E. Aliu, S. Archambault, A. Archer, W. Benbow, R. Bird, Jonathan Biteau, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, J. V. Cardenzana, M. Cerruti, Xuhui Chen, L. Ciupik, M. P. Connolly, W. Cui, H. J. Dickinson, J. D. Eisch, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, A. Flinders, P. Fortin, L. Fortson, A. Furniss, G. H. Gillanders, S. Griffin, J. Grube, G. Gyuk, M. Huetten, Nils Hakansson, J. Holder, T. B. Humensky, C. A. Johnson, P. Kaaret, P. Kar, N. Kelley-Hoskins, M. Kertzman, D. Kieda, M. Krause, M. J. Lang, A. Loo, G. Maier, S. McArthur, A. McCann, K. Meagher, P. Moriarty, R. Mukherjee, T. Nguyen, D. Nieto, R. A. Ong, A. N. Otte, D. Pandel, N. Park, V. Pelassa, A. Petrashyk, M. Pohl, A. Popkow, Elisa Pueschel, J. Quinn, K. Ragan, P. T. Reynolds, G. T. Richards, E. Roache, C. Rulten, M. Santander, G. H. Sembroski, K. Shahinyan, A. W. Smith, D. Staszak, Igor O. Telezhinsky, J. V. Tucci, J. Tyler, A. Varlotta, S. Vincent, S. P. Wakely, O. M. Weiner, A. Weinstein, Alina Wilhelm, D. A. Williams, B. Zitzer, M. Chernyakova, M. S. E. Roberts

• The binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259-63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to beThe binary millisecond radio pulsar PSR J1023+0038 exhibits many characteristics similar to the gamma-ray binary system PSR B1259-63/LS 2883, making it an ideal candidate for the study of high-energy nonthermal emission. It has been the subject of multiwavelength campaigns following the disappearance of the pulsed radio emission in 2013 June, which revealed the appearance of an accretion disk around the neutron star. We present the results of very high energy (VHE) gamma-ray observations carried out by the Very Energetic Radiation Imaging Telescope Array System before and after this change of state. Searches for steady and pulsed emission of both data sets yield no significant gamma-ray signal above 100 GeV, and upper limits are given for both a steady and pulsed gamma-ray flux. These upper limits are used to constrain the magnetic field strength in the shock region of the PSR J1023+0038 system. Assuming that VHE gamma rays are produced via an inverse Compton mechanism in the shock region, we constrain the shock magnetic field to be greater than similar to 2 G before the disappearance of the radio pulsar and greater than similar to 10 G afterward.…