We present a study of γ-ray emission from the core-collapse supernova remnant Cas A in the energy range from 0.1 GeV to 10 TeV. We used 65 hr of the Very Energetic Radiation Imaging Telescope Array ...System (VERITAS) data to cover 200 GeV-10 TeV, and 10.8 yr of Fermi-Large Area Telescope (LAT) data to cover 0.1-500 GeV. The spectral analysis of Fermi-LAT data shows a significant spectral curvature around 1.3 0.4stat GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from Fermi-LAT and VERITAS deviates significantly from a simple power law, and it is best described by a power law with a spectral index of 2.17 0.02stat and a cutoff energy of 2.3 0.5stat TeV. These results, along with radio, X-ray, and γ-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed γ-ray spectrum. From modeling of the entire multiwavelength spectrum, a minimum magnetic field inside the remnant of Bmin 150 G is deduced.
Abstract
Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their ...conventional parameter space (100 GeV ≲
M
χ
≲ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM;
M
χ
≳ 100 TeV) has been suggested as an underexplored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS
γ
-ray observatory. With 216 hr of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a
γ
-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle.
Context.
The response of imaging atmospheric Cherenkov telescopes to incident
γ
-ray-initiated showers in the atmosphere changes as the telescopes age due to exposure to light and weather. These ...aging processes affect the reconstructed energies of the events and
γ
-ray fluxes.
Aims.
This work discusses the implementation of signal calibration methods for the Very Energetic Radiation Imaging Telescope Array System (VERITAS) to account for changes in the optical throughput and detector performance over time.
Methods.
The total throughput of a Cherenkov telescope is the product of camera-dependent factors, such as the photomultiplier tube gains and their quantum efficiencies, and the mirror reflectivity and Winston cone response to incoming radiation. This document summarizes different methods to determine how the camera gains and mirror reflectivity have evolved over time and how we can calibrate this changing throughput in reconstruction pipelines for imaging atmospheric Cherenkov telescopes. The implementation is validated against seven years of observations with the VERITAS telescopes of the Crab Nebula, which is a reference object in very-high-energy astronomy.
Results.
Regular optical throughput monitoring and the corresponding signal calibrations are found to be critical for the reconstruction of extensive air shower images. The proposed implementation is applied as a correction to the signals of the photomultiplier tubes in the telescope simulation to produce fine-tuned instrument response functions. This method is shown to be effective for calibrating the acquired
γ
-ray data and for recovering the correct energy of the events and photon fluxes. At the same time, it keeps the computational effort of generating Monte Carlo simulations for instrument response functions affordably low.
Abstract
G106.3+2.7, commonly considered to be a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct (“head” and “tail”) regions in the ...radio band. A discovery of very-high-energy gamma-ray emission (
E
γ
> 100 GeV) followed by the recent detection of ultrahigh-energy gamma-ray emission (
E
γ
> 100 TeV) from the tail region suggests that G106.3+2.7 is a PeVatron candidate. We present a comprehensive multiwavelength study of the Boomerang PWN (100″ around PSR J2229+6114) using archival radio and Chandra data obtained two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi-LAT and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with Γ = 1.52 ± 0.06 up to 20 keV. Contrary to the previous radio study by Kothes et al., we prefer a much lower PWN
B
-field (
B
∼ 3
μ
G) and larger distance (
d
∼ 8 kpc) based on (1) the nonvarying X-ray flux over the last two decades, (2) the energy-dependent X-ray size of the PWN resulting from synchrotron burn-off, and (3) the multiwavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock ∼1000 yr ago. In this case, the head region should be formed by GeV–TeV electrons injected earlier by the pulsar propagating into the low-density environment.
Abstract
Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ∼ 10–100 times greater than ordinary core-collapse supernovae. One popular model to explain the ...enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma-rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar’s spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV–100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA.
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.
The extragalactic background light (EBL), a diffuse photon field in the optical and infrared range, is a record of radiative processes over the universe's history. Spectral measurements of blazars at ...very high energies (>100 GeV) enable the reconstruction of the spectral energy distribution (SED) of the EBL, as the blazar spectra are modified by redshift- and energy-dependent interactions of the gamma-ray photons with the EBL. The spectra of 14 VERITAS-detected blazars are included in a new measurement of the EBL SED that is independent of EBL SED models. The resulting SED covers an EBL wavelength range of 0.56-56 m, and is in good agreement with lower limits obtained by assuming that the EBL is entirely due to radiation from cataloged galaxies.
Blazars are known for their variability on a wide range of timescales at all wavelengths. Most studies of TeV gamma-ray blazars focus on short timescales, especially during flares. With a decade of ...observations from the Fermi-LAT and VERITAS, we present an extensive study of the long-term multiwavelength radio-to-gamma-ray flux-density variability, with the addition of a couple of short-time radio-structure and optical polarization observations of the blazar 1ES 1215+303 (z = 0.130), with a focus on its gamma-ray emission from 100 MeV to 30 TeV. Multiple strong GeV gamma-ray flares, a long-term increase in the gamma-ray and optical flux baseline, and a linear correlation between these two bands are observed over the ten-year period. Typical HBL behaviors are identified in the radio morphology and broadband spectrum of the source. Three stationary features in the innermost jet are resolved by Very Long Baseline Array at 43.1, 22.2, and 15.3 GHz. We employ a two-component synchrotron self-Compton model to describe different flux states of the source, including the epoch during which an extreme shift in energy of the synchrotron peak frequency from infrared to soft X-rays is observed.
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.
Abstract
HESS J0632+057 belongs to a rare subclass of binary systems that emit gamma rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp “primary” ...peak and broader “secondary” peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between 2019 December and 2020 February, when the orbital phase (
ϕ
) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multiwavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from 2019 October to 2020 March, the MDM observatory did not detect significant variation in H
α
and H
β
line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light-curve data with an intrabinary shock model constrained the orbital parameters, suggesting two orbital phases (at
ϕ
D
= 0.13 and 0.37), where the pulsar crosses the Be-disk, as well as phases for the periastron (
ϕ
0
= 0.30) and inferior conjunction (
ϕ
IFC
= 0.75). The broadband X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases.