• Very High Energy outburst of Markarian 501 in May 2009(1608.01569)

Aug. 4, 2016 astro-ph.HE
The very high energy (VHE; E $>$ 100 GeV) blazar Markarian 501 was observed between April 17 and May 5 (MJD 54938--54956), 2009, as part of an extensive multi-wavelength campaign from radio to VHE. Strong VHE $\gamma$-ray activity was detected on May 1st with Whipple and VERITAS, when the flux (E $>$ 400 GeV) increased to 10 times the pre-flare baseline flux ($3.9{\times 10^{-11}}~{\rm ph~cm^{-2}~s^{-1}}$), reaching five times the flux of the Crab Nebula. This coincided with a decrease in the optical polarization and a rotation of the polarization angle by 15$^{\circ}$. This VHE flare showed a fast flux variation with an increase of a factor $\sim$4 in 25 minutes, and a falling time of $\sim$50 minutes. We present the observations of the quiescent state previous to the flare and of the high state after the flare, focusing on the flux and spectral variability from Whipple, VERITAS, Fermi-LAT, RXTE, and Swift combined with optical and radio data.
• VERITAS Deep Observations of the Dwarf Spheroidal Galaxy Segue 1(1202.2144)

July 7, 2015 hep-ex, astro-ph.HE
The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant $\gamma$-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are $\mathrm{<\sigma v >^{95% CL} \lesssim 10^{-23} cm^{3} s^{-1}}$, improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of two for dark matter particle masses $\mathrm{m_{\chi}\gtrsim 300 GeV}$. The lower limits on the decay lifetime are at the level of $\mathrm{\tau^{95% CL} \gtrsim 10^{24} s}$. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.
• Observations of the unidentified gamma-ray source TeV J2032+4130 by VERITAS(1401.2828)

Jan. 16, 2014 astro-ph.HE
TeV J2032+4130 was the first unidentified source discovered at very high energies (VHE; E $>$ 100 GeV), with no obvious counterpart in any other wavelength. It is also the first extended source to be observed in VHE gamma rays. Following its discovery, intensive observational campaigns have been carried out in all wavelengths in order to understand the nature of the object, which have met with limited success. We report here on a deep observation of TeV J2032+4130, based on 48.2 hours of data taken from 2009 to 2012 by the VERITAS (Very Energetic Radiation Imaging Telescope Array System) experiment. The source is detected at 8.7 standard deviations ($\sigma$) and is found to be extended and asymmetric with a width of 9.5$^{\prime}$$\pm1.2^{\prime} along the major axis and 4.0^{\prime}$$\pm$0.5$^{\prime}$ along the minor axis. The spectrum is well described by a differential power law with an index of 2.10 $\pm$ 0.14$_{stat}$ $\pm$ 0.21$_{sys}$ and a normalization of (9.5 $\pm$ 1.6$_{stat}$ $\pm$ 2.2$_{sys}$) $\times$ 10$^{-13}$TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. We interpret these results in the context of multiwavelength scenarios which particularly favor the pulsar wind nebula (PWN) interpretation.
• HESS J0632+057 is the only gamma-ray binary known so far whose position in the sky allows observations with ground-based observatories both in the northern and southern hemispheres. Here we report on long-term observations of HESS J0632+057 conducted with the VERITAS and H.E.S.S. Cherenkov Telescopes and the X-ray Satellite Swift, spanning a time range from 2004 to 2012 and covering most of the system's orbit. The VHE emission is found to be variable, and is correlated with that at X-ray energies. An orbital period of $315 ^{+6}_{-4}$ days is derived from the X-ray data set, which is compatible with previous results, $P = (321 \pm 5$) days. The VHE light curve shows a distinct maximum at orbital phases close to 0.3, or about 100 days after periastron passage, which coincides with the periodic enhancement of the X-ray emission. Furthermore, the analysis of the TeV data shows for the first time a statistically significant ($> 6.5 \sigma$) detection at orbital phases 0.6--0.9. The obtained gamma-ray and X-ray light curves and the correlation of the source emission at these two energy bands are discussed in the context of the recent ephemeris obtained for the system. Our results are compared to those reported for other gamma-ray binaries.
• Multiwavelength Observations of The TeV Binary LS I +61 303 with VERITAS, Fermi-LAT and Swift-XRT During a TeV Outburst(1310.7913)

Oct. 29, 2013 astro-ph.HE
We present the results of a multiwavelength observational campaign on the TeV binary system LS I +61 303 with the VERITAS telescope array (>200 GeV), Fermi-LAT (0.3-300 GeV), and Swift-XRT (2-10 keV). The data were taken from December 2011 through January 2012 and show a strong detection in all three wavebands. During this period VERITAS obtained 24.9 hours of quality selected livetime data in which LS I +61 303 was detected at a statistical sig- nificance of 11.9 sigma. These TeV observations show evidence for nightly variability in the TeV regime at a post-trial significance of 3.6 sigma. The combination of the simultaneously obtained TeV and X-ray fluxes do not demonstrate any evidence for a correlation between emission in the two bands. For the first time since the launch of the Fermi satellite in 2008, this TeV detection allows the construction of a detailed MeV-TeV spectral energy distribution from LS I +61 303. This spectrum shows a distinct cutoff in emission near 4 GeV, with emission seen by the VERITAS observations following a simple power-law above 200 GeV. This feature in the spectrum of LS I +61 303, obtained from overlapping observations with Fermi-LAT and VERITAS, may indicate that there are two distinct populations of accelerated particles producing the GeV and TeV emission.
• Long term observations of B2 1215+30 with VERITAS(1310.6498)

Oct. 24, 2013 astro-ph.HE
We report on VERITAS observations of the BL Lac object B2 1215+30 between 2008 and 2012. During this period, the source was detected at very high energies (VHE; E > 100 GeV) by VERITAS with a significance of $8.9\sigma$ and showed clear variability on time scales larger than months. In 2011, the source was found to be in a relatively bright state and a power-law fit to the differential photon spectrum yields a spectral index of $3.6 \pm 0.4_{\mathrm{stat}} \pm 0.3_{\mathrm{syst}}$ with an integral flux above 200 GeV of $(8.0 \pm 0.9_{\mathrm{stat}} \pm 3.2_{\mathrm{syst}}) \times 10^{-12}\, \mathrm{cm}^{-2} \mathrm{s}^{-1}$. No short term variability could be detected during the bright state in 2011. Multi-wavelength data were obtained contemporaneous with the VERITAS observations in 2011 and cover optical (Super-LOTIS, MDM, Swift-UVOT), X-ray (Swift-XRT), and gamma-ray (Fermi-LAT) frequencies. These were used to construct the spectral energy distribution (SED) of B2 1215+30. A one-zone leptonic model is used to model the blazar emission and the results are compared to those of MAGIC from early 2011 and other VERITAS-detected blazars. The SED can be well reproduced with model parameters typical for VHE-detected BL Lacs.
• Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil
• Multiwavelength observations and modeling of 1ES 1959+650 in a low flux state(1307.6772)

July 25, 2013 astro-ph.HE
We report on the VERITAS observations of the high-frequency peaked BL Lac object 1ES 1959+650 in the period 2007-2011. This source is detected at TeV energies by VERITAS at 16.4 standard deviation (sigma) significance in 7.6 hours of observation in a low flux state. A multiwavelength spectral energy distribution (SED) is constructed from contemporaneous data from VERITAS, Fermi-LAT, RXTE PCA, and Swift UVOT. Swift XRT data is not included in the SED due to a lack of simultaneous observations with VERITAS. In contrast to the orphan gamma-ray flare exhibited by this source in 2002, the X-ray flux of the source is found to vary by an order of magnitude, while other energy regimes exhibit less variable emission. A quasi-equilibrium synchrotron self-Compton model with an additional external radiation field is used to describe three SEDs corresponding to the lowest, highest, and average X-ray states. The variation in the X-ray spectrum is modeled by changing the electron injection spectral index, with minor adjustments of the kinetic luminosity in electrons. This scenario produces small-scale flux variability of order >~2 in the HE (E>1 MeV) and VHE (E>100 GeV) gamma-ray regimes, which is corroborated by the Fermi-LAT, VERITAS, and Whipple 10m telescope light curves.
• Discovery of TeV Gamma-ray Emission Toward Supernova Remnant SNR G78.2+2.1(1305.6508)

May 28, 2013 astro-ph.HE
We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0.23^{\circ} \pm 0.03^{\circ} (stat)+0.04^{\circ}_{-0.02}^{\circ}(sys) and its spectrum is well-characterized by a differential power law (dN/dE = N_0 \times (E/TeV)^{-\Gamma}) with a photon index of {\Gamma} = 2.37 \pm 0.14 (stat) \pm 0.20 (sys) and a flux normalization of N0 = 1.5 \pm 0.2 (stat) \pm 0.4(sys) \times 10^-12 ph TeV^{-1} cm^{-2} s^{-1}. This yields an integral flux of 5.2 \pm 0.8 (stat) \pm 1.4 (sys) \times 10^-12 ph cm^{-2} s^{-1} above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the supernova remnant shock.
• Constraints on Cosmic Rays, Magnetic Fields, and Dark Matter from Gamma-Ray Observations of the Coma Cluster of Galaxies with VERITAS and Fermi(1208.0676)

Aug. 3, 2012 astro-ph.CO, astro-ph.HE
Observations of radio halos and relics in galaxy clusters indicate efficient electron acceleration. Protons should likewise be accelerated, suggesting that clusters may also be sources of very high-energy (VHE; E>100 GeV) gamma-ray emission. We report here on VHE gamma-ray observations of the Coma galaxy cluster with the VERITAS array of imaging Cherenkov telescopes, with complementing Fermi-LAT observations at GeV energies. No significant gamma-ray emission from the Coma cluster was detected. Integral flux upper limits at the 99% confidence level were measured to be on the order of (2-5)*10^-8\ ph. m^-2 s^-1 (VERITAS, >220 GeV} and ~2*10^-6 ph. m^-2 s^-1 (Fermi, 1-3 GeV), respectively. We use the gamma-ray upper limits to constrain CRs and magnetic fields in Coma. Using an analytical approach, the CR-to-thermal pressure ratio is constrained to be < 16% from VERITAS data and < 1.7% from Fermi data (averaged within the virial radius). These upper limits are starting to constrain the CR physics in self-consistent cosmological cluster simulations and cap the maximum CR acceleration efficiency at structure formation shocks to be <50%. Assuming that the radio-emitting electrons of the Coma halo result from hadronic CR interactions, the observations imply a lower limit on the central magnetic field in Coma of (2 - 5.5) muG, depending on the radial magnetic-field profile and on the gamma-ray spectral index. Since these values are below those inferred by Faraday rotation measurements in Coma (for most of the parameter space), this {renders} the hadronic model a very plausible explanation of the Coma radio halo. Finally, since galaxy clusters are dark-matter (DM) dominated, the VERITAS upper limits have been used to place constraints on the thermally-averaged product of the total self-annihilation cross section and the relative velocity of the DM particles, <\sigma v>. (abr.)
• VERITAS Observations of the Nova in V407 Cygni(1205.5287)

May 23, 2012 astro-ph.GA, astro-ph.HE
We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.
• VERITAS Observations of day-scale flaring of M87 in 2010 April(1112.4518)

Dec. 19, 2011 astro-ph.CO
VERITAS has been monitoring the very-high-energy (VHE; >100GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +- 0.16) X 10^{-11} cm^{-2} s^{-1} at energies above 350GeV. In 2010 April, VERITAS detected a flare from M87 with peak flux of (2.71 +- 0.68) X 10^{-11} cm^{-2} s^{-1} for E>350GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90^{+0.22}_{-0.15} days. The shortest detected exponential rise time is three times as long, at 2.87^{+1.65}_{-0.99} days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M87 and the constraints they provide on the physical size of the emitting region.
• Multiwavelength Observations of the Previously Unidentified Blazar RX J0648.7+1516(1110.5949)

Nov. 9, 2011 astro-ph.HE
We report on the VERITAS discovery of very-high-energy (VHE) gamma- ray emission above 200 GeV from the high-frequency-peaked BL Lac object RXJ0648.7+1516 (GBJ0648+1516), associated with 1FGLJ0648.8+1516. The photon spectrum above 200 GeV is fit by a power law dN/dE = F0(E/E0)-{\Gamma} with a photon index {\Gamma} of 4.4 {\pm} 0.8stat {\pm}0.3syst and a flux normalization F0 of (2.3 {\pm}0.5stat {\pm} 1.2sys) {\times} 10-11 TeV-1cm-2s-1 with E0 = 300 GeV. No VHE vari- ability is detected during VERITAS observations of RXJ0648.7+1516 between 2010 March 4 and April 15. Following the VHE discovery, the optical identifica- tion and spectroscopic redshift were obtained using the Shane 3-m Telescope at the Lick Observatory, showing the unidentified object to be a BL Lac type with a redshift of z = 0.179. Broadband multiwavelength observations contemporaneous with the VERITAS exposure period can be used to sub-classify the blazar as a high-frequency-peaked BL Lac (HBL) object, including data from the MDM ob- servatory, Swift-UVOT and XRT, and continuous monitoring at photon energies above 1 GeV from the Fermi Large Area Telescope (LAT). We find that in the absence of undetected, high-energy rapid variability, the one-zone synchrotron self-Compton model (SSC) overproduces the high-energy gamma-ray emission measured by the Fermi-LAT over 2.3 years. The SED can be parameterized sat- isfactorily with an external-Compton or lepto-hadronic model, which have two and six additional free parameters, respectively, compared to the one-zone SSC model.
• VERITAS Observations of the Unusual Extragalactic Transient Swift J164449.3+573451(1107.1738)

July 8, 2011 astro-ph.HE
We report on very-high-energy ($>$100 GeV) gamma-ray observations of Swift J164449.3+573451, an unusual transient object first detected by the {\it Swift} Observatory and later detected by multiple radio, optical and X-ray observatories. A total exposure of 28 hours was obtained on Swift J164449.3+573451 with VERITAS during 2011 March 28 -- April 15. We do not detect the source and place a differential upper limit on the emission at 500 GeV during these observations of $1.4 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ (99% confidence level). We also present time-resolved upper limits and use a flux limit averaged over the X-ray flaring period to constrain various emission scenarios that can accommodate both the radio-through-X-ray emission detected from the source and the lack of detection by VERITAS.