skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An In-depth Study of Gamma Rays from the Starburst Galaxy M82 with VERITAS
Abstract Assuming Galactic cosmic rays originate in supernovae and the winds of massive stars, starburst galaxies should produce very-high-energy (VHE;E > 100 GeV) gamma-ray emission via the interaction of their copious quantities of cosmic rays with the large reservoirs of dense gas within the galaxies. Such VHE emission was detected by VERITAS from the starburst galaxy M82 in 2008–09. An extensive, multiyear campaign followed these initial observations, yielding a total of 254 hr of good-quality VERITAS data on M82. Leveraging modern analysis techniques and the larger exposure, these VERITAS data show a more statistically significant VHE signal (∼6.5 standard deviations,σ). The corresponding photon spectrum is well fit by a power law (Γ = 2.3 ± 0.3stat ± 0.2sys), and the observed integral flux isF(>450 GeV) = (3.2 ± 0.6stat ± 0.6sys) × 10−13cm−2s−1, or ∼0.4% of the Crab Nebula flux above the same energy threshold. The improved VERITAS measurements, when combined with various multiwavelength data, enable modeling of the underlying emission and transport processes. A purely leptonic scenario is found to be a poor representation of the gamma-ray spectral energy distribution (SED). A lepto-hadronic scenario with cosmic rays following a power-law spectrum in momentum (indexs ≃ 2.25) and with significant bremsstrahlung below 1 GeV provides a good match to the observed SED. The synchrotron emission from the secondary electrons indicates that efficient nonradiative losses of cosmic-ray electrons may be related to advective escape from the starburst core.  more » « less
Award ID(s):
2200857 2411023 2310028 2411760 2310158 2411860 2209437 2310002 2209605 2110497 2111531 2110737
PAR ID:
10600409
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Corporate Creator(s):
Publisher / Repository:
American Astronomical Society
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
981
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
189
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source, which includes very high energy (VHE, $$E\gt $$ 100 GeV) $$\gamma$$-ray data taken by the MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov telescopes) and H.E.S.S. (High Energy Stereoscopic System) imaging Cherenkov telescopes. The discovery of VHE $$\gamma$$-ray emission happened during a high state of $$\gamma$$-ray activity in July 2016, observed by many instruments from radio to VHE $$\gamma$$-rays. We identify four states of activity of the source, one of which includes VHE $$\gamma$$-ray emission. Variability in the VHE domain is found on daily time-scales. The intrinsic VHE spectrum can be described by a power law with index $$3.27\pm 0.44_{\rm stat}\pm 0.15_{\rm sys}$$ (MAGIC) and $$3.39\pm 0.58_{\rm stat}\pm 0.64_{\rm sys}$$ (H.E.S.S.) in the energy range of 55–300 and 120–500 GeV, respectively. The broadband emission cannot be successfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the data set well and a proton-synchrotron-dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the broad-line region to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be a flat spectrum radio quasar (FSRQ), in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL Lac and FSRQ objects. 
    more » « less
  2. Abstract Observations of GeV gamma-ray emission from the well-studied mixed-morphology supernova remnant (SNR) W44 by Fermi-Large Area Telescope and AGILE imply that it is a site of significant cosmic-ray acceleration. The spectral energy distribution (SED) derived from the GeV data suggests that the gamma-ray emission likely originates from the decay of neutral pions generated by cosmic-ray interactions. It is essential to measure the SED of W44 in the X-ray and very-high-energy (VHE) gamma-ray bands to verify the hadronic origin of the emission and to gauge the potential contributions from leptonic emission. We report an upper limit of the nonthermal X-ray flux from W44 of 5  × 10−13erg cm−2s−1in the 0.5–8.0 keV band based on  ∼300 ks of XMM-Newton observations. The X-ray upper limit is consistent with previously estimated hadronic models, but in tension with the leptonic models. We estimate the VHE flux upper limit of  ∼1.2  × 10−12erg s−1cm−2in the 0.5–5.0 TeV range from W44 using data from the Very Energetic Radiation Imaging Telescope Array System. Our nondetection of W44 at VHE wavelengths is in agreement with observations from other imaging atmospheric Cherenkov telescopes and is perhaps consistent with the evolutionary stage of the SNR. 
    more » « less
  3. Aims.Mrk 421 was in its most active state around early 2010, which led to the highest TeV gamma-ray flux ever recorded from any active galactic nuclei (AGN). We aim to characterize the multiwavelength behavior during this exceptional year for Mrk 421, and evaluate whether it is consistent with the picture derived with data from other less exceptional years. Methods.We investigated the period from November 5, 2009, (MJD 55140) until July 3, 2010, (MJD 55380) with extensive coverage from very-high-energy (VHE;E > 100 GeV) gamma rays to radio with MAGIC, VERITAS,Fermi-LAT,RXTE,Swift, GASP-WEBT, VLBA, and a variety of additional optical and radio telescopes. We characterized the variability by deriving fractional variabilities as well as power spectral densities (PSDs). In addition, we investigated images of the jet taken with VLBA and the correlation behavior among different energy bands. Results.Mrk 421 was in widely different states of activity throughout the campaign, ranging from a low-emission state to its highest VHE flux ever recorded. We find the strongest variability in X-rays and VHE gamma rays, and PSDs compatible with power-law functions with indices around 1.5. We observe strong correlations between X-rays and VHE gamma rays at zero time lag with varying characteristics depending on the exact energy band. We also report a marginally significant (∼3σ) positive correlation between high-energy (HE;E > 100 MeV) gamma rays and the ultraviolet band. We detected marginally significant (∼3σ) correlations between the HE and VHE gamma rays, and between HE gamma rays and the X-ray, that disappear when the large flare in February 2010 is excluded from the correlation study, hence indicating the exceptionality of this flaring event in comparison with the rest of the campaign. The 2010 violent activity of Mrk 421 also yielded the first ejection of features in the VLBA images of the jet of Mrk 421. Yet the large uncertainties in the ejection times of these unprecedented radio features prevent us from firmly associating them to the specific flares recorded during the 2010 campaign. We also show that the collected multi-instrument data are consistent with a scenario where the emission is dominated by two regions, a compact and extended zone, which could be considered as a simplified implementation of an energy-stratified jet as suggested by recentIXPEobservations. 
    more » « less
  4. 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 PWNB-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. 
    more » « less
  5. 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. 
    more » « less