Cosmic-ray electrons and positrons (CREs) at GeV-TeV energies are a unique probe of our local Galactic neighborhood. CREs lose energy rapidly via synchrotron radiation and inverse-Compton scattering ...processes while propagating within the Galaxy, and these losses limit their propagation distance. For electrons with TeV energies, the limit is on the order of a kiloparsec. Within that distance, there are only a few known astrophysical objects capable of accelerating electrons to such high energies. It is also possible that the CREs are the products of the annihilation or decay of heavy dark matter (DM) particles. VERITAS, an array of imaging air Cherenkov telescopes in southern Arizona, is primarily utilized for gamma-ray astronomy but also simultaneously collects CREs during all observations. We describe our methods of identifying CREs in VERITAS data and present an energy spectrum, extending from 300 GeV to 5 TeV, obtained from approximately 300 hours of observations. A single power-law fit is ruled out in VERITAS data. We find that the spectrum of CREs is consistent with a broken power law, with a break energy at 710±40stat±140syst GeV.
We present constraints on the annihilation cross section of weakly interacting massive particles dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results ...are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of ∼230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is 1.35×10−23 cm3 s−1 at 1 TeV for the bottom quark (bb¯) final state, 2.85×10−24 cm3 s−1 at 1 TeV for the tau lepton (τ+τ−) final state and 1.32×10−25 cm3 s−1 at 1 TeV for the gauge boson (γγ) final state.
Abstract
Flat-spectrum radio quasars (FSRQs) are the most luminous blazars at GeV energies but only rarely emit detectable fluxes of TeV gamma rays, typically during bright GeV flares. We explore the ...gamma-ray variability and spectral characteristics of three FSRQs that have been observed at GeV and TeV energies by Fermi-LAT and VERITAS, making use of almost 100 hr of VERITAS observations spread over 10 yr: 3C 279, PKS 1222+216, and Ton 599. We explain the GeV flux distributions of the sources in terms of a model derived from a stochastic differential equation describing fluctuations in the magnetic field in the accretion disk and estimate the timescales of magnetic flux accumulation and stochastic instabilities in their accretion disks. We identify distinct flares using a procedure based on Bayesian blocks and analyze their daily and subdaily variability and gamma-ray energy spectra. Using observations from VERITAS, as well as Fermi, Swift, and the Steward Observatory, we model the broadband spectral energy distributions of PKS 1222+216 and Ton 599 during very high energy (VHE)–detected flares in 2014 and 2017, respectively, strongly constraining the jet Doppler factors and gamma-ray emission region locations during these events. Finally, we place theoretical constraints on the potential production of PeV-scale neutrinos during these VHE flares.
We present the low-temperature Geiger-mode characteristics of GaN (gallium nitride) p-i-n avalanche photodiodes (APDs). The devices have a breakdown voltage of −95V and a temperature dependence of ...<inline-formula> <tex-math notation="LaTeX">\sim </tex-math></inline-formula>0.0159 ± 0.0034 V/K near 300K. The room-temperature (300 K) dark-count rate (DCR) is 23.8 MHz for a <inline-formula> <tex-math notation="LaTeX">75\times 75\,\,\mu \text{m}^{2} </tex-math></inline-formula> device biased at 1 % overvoltage. The DCR halves when lowering the temperature by 50°C. Based on the temperature-dependent characteristics of the DCR, we identify band-to-band tunneling as the dominant DCR generation mechanism. At 4.65 % overvoltage and 375 nm, the photon detection efficiency (PDE) is 0.82 %- limited by a low breakdown probability of 1.7 %. We discuss the measurement setup and the method to extract count rates, which is based on the Poisson distribution of the time intervals between Geiger-mode breakdowns of the APD. The setup includes a custom circuit to bias the diode and amplify its signals, a steady ultraviolet (UV) light source, and a system to control the temperature of the APD with a thermoelectric element in the range from −40 to 20°C.
Abstract
The Galactic Center (GC) region hosts a variety of powerful astronomical sources and rare astrophysical processes that emit a large flux of nonthermal radiation. The inner 375 pc × 600 pc ...region, called the Central Molecular Zone, is home to the supermassive black hole Sagittarius A*, massive cloud complexes, and particle accelerators such as supernova remnants (SNRs). We present the results of our improved analysis of the very-high-energy gamma-ray emission above 2 TeV from the GC using 125 hr of data taken with the Very Energetic Radiation Imaging Telescope Array System imaging-atmospheric Cerenkov telescope between 2010 and 2018. The central source VER J1745–290, consistent with the position of Sagittarius A*, is detected at a significance of 38 standard deviations above the background level (38
σ
), and we report its spectrum and light curve. Its differential spectrum is consistent with a power law with exponential cutoff, with a spectral index of
, a flux normalization at 5.3 TeV of
TeV
−1
cm
−2
s
−1
, and cutoff energy of
TeV. We also present results on the diffuse emission near the GC, obtained by combining data from multiple regions along the GC ridge, which yield a cumulative significance of 9.5
σ
. The diffuse GC ridge spectrum is best fit by a power law with a hard index of 2.19 ± 0.20, showing no evidence of a cutoff up to 40 TeV. This strengthens the evidence for a potential accelerator of PeV cosmic rays being present in the GC. We also provide spectra of the other sources in our field of view with significant detections, composite SNR G0.9+0.1, and HESS J1746–285.
High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it ...is a young, type Ia SNR that has been well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hr of VERITAS and 84 months of Fermi-LAT observations, which represent about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is 2.92 0.42stat 0.20sys. It is also softer than the spectral index in the GeV energy range, 2.14 0.09stat 0.02sys, measured in this study using Fermi-LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi-LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously.
We present results from deep observations toward the Cygnus region using 300 hr of very high energy (VHE) γ-ray data taken with the VERITAS Cerenkov telescope array and over 7 yr of high-energy γ-ray ...data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse γ-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 J2021.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 γ-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula flux. We interpret these observations in a multiwavelength context and present the most detailed γ-ray view of the region to date.
We present a new measurement of the energy spectrum of iron nuclei in cosmic rays from 20 TeV 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 dFdE=f0·(EE0)−γ over the full energy range, with γ=2.82±0.30(stat)−0.27+0.24(syst) and f0=(4.82±0.98(stat)−2.70+2.12(syst))×10−7 m−2 s−1 TeV−1 sr−1 at E0=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.
On 2017 September 22, the IceCube Neutrino Observatory reported the detection of the high-energy neutrino event IC 170922A, of potential astrophysical origin. It was soon determined that the neutrino ...direction was consistent with the location of the gamma-ray blazar TXS 0506+056 (3FGL J0509.4+0541), which was in an elevated gamma-ray emission state as measured by the Fermi satellite. Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations of the neutrino/blazar region started on 2017 September 23 in response to the neutrino alert and continued through 2018 February 6. While no significant very-high-energy (VHE; E > 100 GeV) emission was observed from the blazar by VERITAS in the two-week period immediately following the IceCube alert, TXS 0506+056 was detected by VERITAS with a significance of 5.8 standard deviations ( ) in the full 35 hr data set. The average photon flux of the source during this period was (8.9 1.6) × 10−12 cm−2 s−1, or 1.6% of the Crab Nebula flux, above an energy threshold of 110 GeV, with a soft spectral index of 4.8 1.3.
We present a search for magnetically broadened gamma-ray emission around active galactic nuclei (AGNs), using VERITAS observations of seven hard-spectrum blazars. A cascade process occurs when ...multi-TeV gamma-rays from an AGN interact with extragalactic background light (EBL) photons to produce electron-positron pairs, which then interact with cosmic microwave background photons via inverse-Compton scattering to produce gamma-rays. Due to the deflection of the electron-positron pairs, a non-zero intergalactic magnetic field (IGMF) would potentially produce detectable effects on the angular distribution of the cascade emission. In particular, an angular broadening compared to the unscattered emission could occur. Through non-detection of angularly broadened emission from 1ES 1218+304, the source with the largest predicted cascade fraction, we exclude a range of IGMF strengths around 10−14 G at the 95% confidence level. The extent of the exclusion range varies with the assumptions made about the intrinsic spectrum of 1ES 1218+304 and the EBL model used in the simulation of the cascade process. All of the sources are used to set limits on the flux due to extended emission.