We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels \textit{Swift} Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 to 500 TeV. The measurement makes use of a template-based analysis method, which, for the first time, is applied to the energy reconstruction of iron-induced air showers recorded by the VERITAS array of imaging atmospheric Cherenkov telescopes. The event selection makes use of the direct Cherenkov light which is emitted by charged particles before the first interaction, as well as other parameters related to the shape of the recorded air shower images. The measured spectrum is well described by a power law $\frac{\mathrm{d} F}{\mathrm{d} E}=f_0\cdot \left(\frac{E}{E_0}\right)^{-\gamma}$ over the full energy range, with $\gamma = 2.82 \pm 0.30 \mathrm{(stat.)} ^{+0.24}_{-0.27} \mathrm{(syst.)}$ and $f_0 = \left( 4.82 \pm 0.98 \mathrm{(stat.)}^{+2.12}_{-2.70} \mathrm{(syst.)} \right)\cdot 10^{-7}$m$^{-2}$s$^{-1}$sr$^{-1}$TeV$^{-1}$ at $E_0=50$TeV, with no indication of a cutoff or spectral break. The measured differential flux is compatible with previous results, with improved statistical uncertainty at the highest energies.

HESS J1943+213 is a very-high-energy (VHE; $>$100 GeV) $\gamma$-ray source in the direction of the Galactic Plane. Studies exploring the classification of the source are converging towards its identification as an extreme synchrotron BL Lac object. Here we present 38 hours of VERITAS observations of HESS J1943+213 taken over two years. The source is detected with $\sim$20 standard deviations significance, showing a remarkably stable flux and spectrum in VHE $\gamma$-rays. Multi-frequency very-long-baseline array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with European VLBI Network and detect this component in the 4.6 GHz and the 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, $Fermi$-LAT, $Swift$-XRT, FLWO 48$''$ telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron-self-Compton model. The well-measured $\gamma$-ray peak of the SED with VERITAS and $Fermi$-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary $\gamma$-rays from ultra-high-energy cosmic ray-initiated electromagnetic cascades to the $\gamma$-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and $\gamma$-ray bands. No statistically significant flux or spectral variability is detected.

We present results from deep observations towards the Cygnus region using 300 hours of very-high-energy (VHE) $\gamma$-ray data taken with the VERITAS Cherenkov telescope array and over seven years of high-energy $\gamma$-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse $\gamma$-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130), and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL 2021.0+4031e (the Gamma-Cygni supernova remnant) was examined and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE $\gamma$-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula ux. We interpret these observations in a multiwavelength context and present the most detailed $\gamma$-ray view of the region to date.

Compilation of papers presented by the VERITAS Collaboration at the 34th International Cosmic Ray Conference (ICRC), held July 30 through August 6, 2015 in The Hague, The Netherlands.

The Galactic center is an interesting region for high-energy (0.1-100 GeV) and very-high-energy (E > 100 GeV) gamma-ray observations. Potential sources of GeV/TeV gamma-ray emission have been suggested, e.g., the accretion of matter onto the supermassive black hole, cosmic rays from a nearby supernova remnant (e.g. SgrA East), particle acceleration in a plerion, or the annihilation of dark matter particles. The Galactic center has been detected by EGRET and by Fermi/LAT in the MeV/GeV energy band. At TeV energies, the Galactic center was detected with moderate significance by the CANGAROO and Whipple 10 m telescopes and with high significance by H.E.S.S., MAGIC, and VERITAS. We present the results from three years of VERITAS observations conducted at large zenith angles resulting in a detection of the Galactic center on the level of 18 standard deviations at energies above ~2.5TeV. The energy spectrum is derived and is found to be compatible with hadronic, leptonic and hybrid emission models discussed in the literature. Future, more detailed measurements of the high-energy cutoff and better constraints on the high-energy flux variability will help to refine and/or disentangle the individual models.

In March 2013, a flaring episode from the Crab Nebula lasting ~2 weeks was detected by the Fermi-LAT (Large Area Telescope on board the Fermi Gamma-ray Space Telescope). VERITAS provides simultaneous observations throughout this period. During the flare, the Fermi-LAT detected a 20-fold increase in flux above the average synchrotron flux >100 MeV seen from the Crab Nebula. Simultaneous measurements with VERITAS are consistent with the non-variable long-term average Crab Nebula flux at TeV energies. Assuming a linear correlation between the very-high-energy flux change >1 TeV and the flux change seen in the Fermi-LAT band >100 MeV during the period of simultaneous observations, the linear correlation factor can be constrained to be at most 8.6 * 10^-3 with 95% confidence.

This white paper describes the contributions of the VERITAS gamma-ray telescope array to the CF6-A subgroup of the 'Cosmic Frontier' working group, part of the Snowmass Community Summer Study 2013.

Compilation of papers contributed by the VERITAS Collaboration to the 32nd International Cosmic Ray Conference, held 11-18 August 2011 in Beijing, China.

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including {\gamma}-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

Indirect dark matter searches with ground-based gamma-ray observatories provide an alternative for identifying the particle nature of dark matter that is complementary to that of direct search or accelerator production experiments. We present the results of observations of the dwarf spheroidal galaxies Draco, Ursa Minor, Bootes 1, and Willman 1 conducted by VERITAS. These galaxies are nearby dark matter dominated objects located at a typical distance of several tens of kiloparsecs for which there are good measurements of the dark matter density profile from stellar velocity measurements. Since the conventional astrophysical background of very high energy gamma rays from these objects appears to be negligible, they are good targets to search for the secondary gamma-ray photons produced by interacting or decaying dark matter particles. No significant gamma-ray flux above 200 GeV was detected from these four dwarf galaxies for a typical exposure of ~20 hours. The 95% confidence upper limits on the integral gamma-ray flux are in the range 0.4-2.2x10^-12 photons cm^-2s^-1. We interpret this limiting flux in the context of pair annihilation of weakly interacting massive particles and derive constraints on the thermally averaged product of the total self-annihilation cross section and the relative velocity of the WIMPs. The limits are obtained under conservative assumptions regarding the dark matter distribution in dwarf galaxies and are approximately three orders of magnitude above the generic theoretical prediction for WIMPs in the minimal supersymmetric standard model framework. However significant uncertainty exists in the dark matter distribution as well as the neutralino cross sections which under favorable assumptions could further lower the limits.

M87 is a nearby radio galaxy that is detected at energies ranging from radio to VHE gamma-rays. Its proximity and its jet, misaligned from our line-of-sight, enable detailed morphological studies and extensive modeling at radio, optical, and X-ray energies. Flaring activity was observed at all energies, and multi-wavelength correlations would help clarify the origin of the VHE emission. In this paper, we describe a detailed temporal and spectral analysis of the VERITAS VHE gamma-ray observations of M87 in 2008 and 2009. In the 2008 observing season, VERITAS detected an excess with a statistical significance of 7.2 sigma from M87 during a joint multi-wavelength monitoring campaign conducted by three major VHE experiments along with the Chandra X-ray Observatory. In February 2008, VERITAS observed a VHE flare from M87 occurring over a 4-day timespan. The peak nightly flux above 250GeV was 7.7% of the Crab Nebula flux. M87 was marginally detected before this 4-day flare period, and was not detected afterwards. Spectral analysis of the VERITAS observations showed no significant change in the photon index between the flare and pre-flare states. Shortly after the VHE flare seen by VERITAS, the Chandra X-ray Observatory detected the flux from the core of M87 at a historical maximum, while the flux from the nearby knot HST-1 remained quiescent. Acciari et al. (2009) presented the 2008 contemporaneous VHE gamma-ray, Chandra X-ray, and VLBA radio observations which suggest the core as the most likely source of VHE emission, in contrast to the 2005 VHE flare that was simultaneous with an X-ray flare in the HST-1 knot. In 2009, VERITAS continued its monitoring of M87 and marginally detected a 4.2 sigma excess corresponding to a flux of ~1% of the Crab Nebula. No VHE flaring activity was observed in 2009.

We report on observations of very high-energy gamma rays from the shell-type supernova remnant Cassiopeia A with the VERITAS stereoscopic array of four imaging atmospheric Cherenkov telescopes in Arizona. The total exposure time for these observations is 22 hours, accumulated between September and November of 2007. The gamma-ray source associated with the SNR Cassiopeia A was detected above 200 GeV with a statistical significance of 8.3 s.d. The estimated integral flux for this gamma-ray source is about 3% of the Crab-Nebula flux. The photon spectrum is compatible with a power law dN/dE ~ E^(-Gamma) with an index Gamma = 2.61 +/- 0.24(stat) +/- 0.2(sys). The data are consistent with a point-like source. We provide a detailed description of the analysis results, and discuss physical mechanisms that may be responsible for the observed gamma-ray emission.

We present results from an intensive VERITAS monitoring campaign of the high-frequency peaked BL Lac object 1ES 1218+304 in 2008/2009. Although 1ES 1218+304 was detected previously by MAGIC and VERITAS at a persistent level of ~6% of the Crab Nebula flux, the new VERITAS data reveal a prominent flare reaching ~20% of the Crab. While very high energy (VHE) flares are quite common in many nearby blazars, the case of 1ES 1218+304 (redshift z = 0.182) is particularly interesting since it belongs to a group of blazars that exhibit unusually hard VHE spectra considering their redshifts. When correcting the measured spectra for absorption by the extragalactic background light, 1ES 1218+304 and a number of other blazars are found to have differential photon indices less than 1.5. The difficulty in modeling these hard spectral energy distributions in blazar jets has led to a range of theoretical gamma-ray emission scenarios, one of which is strongly constrained by these new VERITAS observations. We consider the implications of the observed light curve of 1ES 1218+304, which shows day scale flux variations, for shock acceleration scenarios in relativistic jets, and in particular for the viability of kiloparsec-scale jet emission scenarios.

We report the first detection of very-high-energy (VHE) gamma-ray emission above 140 GeV from PKS 1424+240, a BL Lac object with an unknown redshift. The photon spectrum above 140 GeV measured by VERITAS is well described by a power law with a photon index of 3.8 +- 0.5_stat +- 0.3_syst and a flux normalization at 200 GeV of (5.1 +- 0.9_stat +- 0.5_syst) x 10^{-11} TeV^-1 cm^-2 s^-1, where stat and syst denote the statistical and systematical uncertainty, respectively. The VHE flux is steady over the observation period between MJD 54881 and 55003 (2009 February 19 to June 21). Flux variability is also not observed in contemporaneous high energy observations with the Fermi Large Area Telescope (LAT). Contemporaneous X-ray and optical data were also obtained from the Swift XRT and MDM observatory, respectively. The broadband spectral energy distribution (SED) is well described by a one-zone synchrotron self-Compton (SSC) model favoring a redshift of less than 0.1. Using the photon index measured with Fermi in combination with recent extragalactic background light (EBL) absorption models it can be concluded from the VERITAS data that the redshift of PKS 1424+240 is less than 0.66.

Index of the VERITAS contributions to ICRC 2009

The accretion of matter onto a massive black hole is believed to feed the relativistic plasma jets found in many active galactic nuclei (AGN). Although some AGN accelerate particles to energies exceeding 10^12 electron Volts (eV) and are bright sources of very-high-energy (VHE) gamma-ray emission, it is not yet known where the VHE emission originates. Here we report on radio and VHE observations of the radio galaxy M87, revealing a period of extremely strong VHE gamma-ray flares accompanied by a strong increase of the radio flux from its nucleus. These results imply that charged particles are accelerated to very high energies in the immediate vicinity of the black hole.

We report on the results of two coordinated multiwavelength campaigns that focused on the blazar Markarian 421 during its 2006 and 2008 outbursts. These campaigns obtained UV and X-ray data using the XMM-Newton satellite, while the gamma-ray data were obtained utilizing three imaging atmospheric Cerenkov telescopes, the Whipple 10m telescope and VERITAS, both based in Arizona, as well as the MAGIC telescope, based on La Palma in the Canary Islands. The coordinated effort between the gamma-ray groups allowed for truly simultaneous data in UV/X-ray/gamma-ray wavelengths during a significant portion of the XMM-Newton observations. This simultaneous coverage allowed for a reliable search for correlations between UV, X-ray and gamma-ray variability over the course of the observations. Investigations of spectral hysteresis and modeling of the spectral energy distributions are also presented.