TeV Gamma-ray Observations Research Articles

Upper limit on the axion-photon coupling from Markarian 421

Markarian421 is a well-known nearby BL Lac blazar at the redshift z=0.031. Many previous works were investigated to constrain the axion-photon coupling from its TeV gamma-ray observations, showing the upper limit on the coupling constant gaγ≲2.0×10−11GeV−1 for the axion mass [5.0×10−10eV≲ma≲5.0×10−7eV]. While in this work, we obtain a more stringent upper limit on the axion-photon coupling from the 1038 days gamma-ray observations of the blazar Markarian421. The long-term VHE gamma-ray spectra are measured by the collaborations Large Area Telescope on board NASA's Fermi Gamma-ray Space Telescope (Fermi-LAT) and High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory from 2015 June to 2018 July. We show the best-fit spectral energy distributions (SEDs) of Markarian421 under the null and axion hypotheses. Then we set the axion-photon limit in the {ma,gaγ} plane. The 99% C.L. upper limit set by Markarian421 is gaγ≲4.0×10−12GeV−1 for the axion mass [1.0×10−9eV≲ma≲1.0×10−8eV]. It is the most stringent upper limit in this axion mass region.

HAWC Study of the Ultra-high-energy Spectrum of MGRO J1908+06

Abstract We report TeV gamma-ray observations of the ultra-high-energy source MGRO J1908+06 using data from the High Altitude Water Cherenkov Observatory. This source is one of the highest-energy known gamma-ray sources, with emission extending past 200 TeV. Modeling suggests that the bulk of the TeV gamma-ray emission is leptonic in nature, driven by the energetic radio-faint pulsar PSR J1907+0602. Depending on what assumptions are included in the model, a hadronic component may also be allowed. Using the results of the modeling, we discuss implications for detection prospects by multi-messenger campaigns.

TeV Instrumentation: Current and Future

During the last 20 years, TeV astronomy has turned from a fledgling field, with only a handful of sources, into a fully-developed astronomy discipline, broadening our knowledge on a variety of types of TeV gamma-ray sources. This progress has been mainly achieved due to the currently operating instruments: imaging atmospheric Cherenkov telescopes, surface arrays and water Cherenkov detectors. Moreover, we are at the brink of a next generation of instruments, with a considerable leap in performance parameters. This review summarizes the current status of the TeV astronomy instrumentation, mainly focusing on the comparison of the different types of instruments and analysis challenges, as well as providing an outlook into the future installations. The capabilities and limitations of different techniques of observations of TeV gamma rays are discussed, as well as synergies to other bands and messengers.

Investigating γ -ray halos around three HAWC bright sources in Fermi -LAT data

Numerous extended sources around Galactic pulsars have shown significant $\gamma$-ray emission from GeV to TeV energies, revealing hundreds of TeV energy electrons scattering off of the underlying photon fields through inverse Compton scattering (ICS). HAWC TeV gamma-ray observations of few-degree extended emission around the pulsars Geminga and Monogem, and LAT GeV emission around Geminga, suggest that systems older than 100 kyr have multi-TeV $e^\pm$ propagating beyond the SNR-PWN system into the interstellar medium. Following the discovery of few $\gamma$-ray sources by HAWC at energies E$>100$~TeV, we investigate the presence of an extended $\gamma$-ray emission in Fermi-LAT data around the three brightest sources detected by HAWC up to 100 TeV. We find an extended emission of $\theta_{68} = 1.00^{+0.05}_{-0.07}$ deg around eHWC J1825-134 and $\theta_{68} = 0.71\pm0.10$ deg eHWC J1907+063. The analysis with ICS templates on Fermi-LAT data point to diffusion coefficient values which are significantly lower than the average Galactic one. When studied along with HAWC data, the $\gamma$-ray Fermi-LAT data provide invaluable insight into the very high-energy electron and positron parent populations.

Probing the Properties of the Pulsar Wind in the Gamma-Ray Binary HESS J0632+057 with NuSTAR and VERITAS Observations

Abstract HESS J0632+057 is a gamma-ray binary composed of a compact object orbiting a Be star with a period of about 315 days. Extensive X-ray and TeV gamma-ray observations have revealed a peculiar light curve containing two peaks, separated by a dip. We present the results of simultaneous observations in hard X-rays with NuSTAR and in TeV gamma-rays with VERITAS, performed in 2017 November and December. These observations correspond to the orbital phases ϕ ≈ 0.22 and 0.3, where the fluxes are rising toward the first light-curve peak. A significant variation of the spectral index from 1.77 ± 0.05 to 1.56 ± 0.05 is observed in the X-ray data. The multiwavelength spectral energy distributions (SED) derived from the observations are interpreted in terms of a leptonic model, in which the compact object is assumed to be a pulsar and nonthermal radiation is emitted by high-energy electrons accelerated at the shock formed by the collision between the stellar and pulsar wind. The results of the SED fitting show that our data can be consistently described within this scenario, and allow us to estimate the magnetization of the pulsar wind at the location of the shock formation. The constraints on the pulsar wind magnetization provided by our results are shown to be consistent with those obtained from other systems.

Measurement of the Crab Nebula Spectrum Past 100 TeV with HAWC

Abstract We present TeV gamma-ray observations of the Crab Nebula, the standard reference source in ground-based gamma-ray astronomy, using data from the High Altitude Water Cherenkov (HAWC) Gamma-Ray Observatory. In this analysis we use two independent energy estimation methods that utilize extensive air shower variables such as the core position, shower angle, and shower lateral energy distribution. In contrast, the previously published HAWC energy spectrum roughly estimated the shower energy with only the number of photomultipliers triggered. This new methodology yields a much-improved energy resolution over the previous analysis and extends HAWC’s ability to accurately measure gamma-ray energies well beyond 100 TeV. The energy spectrum of the Crab Nebula is well fit to a log-parabola shape with emission up to at least 100 TeV. For the first estimator, a ground parameter that utilizes fits to the lateral distribution function to measure the charge density 40 m from the shower axis, the best-fit values are (TeV cm2 s)−1, , and . For the second estimator, a neural network that uses the charge distribution in annuli around the core and other variables, these values are (TeV cm2 s)−1, , and β = 0.06 ± 0.01 ± 0.02. The first set of uncertainties is statistical; the second set is systematic. Both methods yield compatible results. These measurements are the highest-energy observation of a gamma-ray source to date.

Constraints on axion-like particle properties with TeV gamma-ray observations of Galactic sources

Axion-like particles (ALPs) can oscillate into photons and vice versa in the magnetic fields, which provides a feasible explanation of the observed transparency of the Universe for TeV gamma rays. The coupling gaγ for the ALP mass ma⩽ 10−7 eV has been tightly constrained by the Fermi-LAT and H.E.S.S. observations of some extragalactic sources. In this work, we search for the ALP oscillation effect with the H.E.S.S. observation data of some TeV sources in the Galactic plane for the first time. No significant ALP signal is found. The constraints on the ALP parameters, though sensitively depend on the Galactic magnetic field (GMF) model, are more stringent than the limits set by the CAST experiment at ma ∼ 0.1μ eV. The upcoming CTA observations of these Galactic TeV sources are able to improve the constraints substantially.

Searching for secluded dark matter with H.E.S.S., Fermi-LAT, and Planck

Short-lived mediators are often used to describe dark matter interactions with Standard Model particles. When the dark matter mass is heavier than the mass of the mediator, it may self-annihilate into short-lived mediators, and in some cases this might be the dominant annihilation channel. This scenario is known as secluded dark matter. We use Fermi-LAT observations of dwarf spheroidal galaxies, H.E.S.S. data from the Galactic center, and Planck measurements of the Cosmic Microwave Background to constrain secluded dark matter. We explore the interplay between these experiments and we assess the impact of the mediator mass on our bounds, an often overlooked yet very important point. In particular, we exclude pair -annihilation cross-sections greater or on the order of σ v ∼ 4 × 10−27 cm3/s for dark matter masses around 10 GeV and greater or on the order of σ v ∼ × 10−25 cm3/s for dark matter masses around a TeV . Our findings supersede previous constraints which use Fermi-LAT data, and constitute the first limits on secluded dark sectors using the H.E.S.S. telescope. We also show that one can fit TeV gamma-ray observations from H.E.S.S. with secluded dark matter annihilations, with the mediator mass impacting the best-fit dark matter particle mass. Our findings indicate that any assessment of secluded dark sectors in the context of indirect detection significantly depends on the choice of the mediator mass.

Improved bound on isotropic Lorentz violation in the photon sector from extensive air showers

Cosmic rays have extremely high particle energies (up to $10^{20} \; \text{eV}$) and can be used to search for violations of Lorentz invariance. We consider isotropic nonbirefringent Lorentz violation in the photon sector for the case of a photon velocity larger than the maximum attainable velocity of the standard fermions. Up to now, Earth-based bounds on this type of Lorentz violation have been determined from observations of TeV gamma rays. Here, we elaborate on a novel approach to test Lorentz invariance with greatly improved sensitivity. This approach is based on investigating extensive air showers which are induced by cosmic-ray particles in the Earth's atmosphere. We study the impact of two Lorentz-violating decay processes on the longitudinal development of air showers, notably the atmospheric depth of the shower maximum $X_\text{max}$. Specifically, the two Lorentz-violating decay processes considered are photon decay into an electron-positron pair and modified neutral-pion decay into two photons. We use Monte Carlo simulations performed with the CONEX code which was extended to include these two Lorentz-violating decay processes at a magnitude allowed by the best previous Earth-based bound. Compared to standard physics, these Lorentz-violating decay processes reduce the average $X_\text{max}$ for showers with primary energies above $10^{18}\;\text{eV}$ by an amount that is significantly larger than the average resolution of current air shower experiments. Comparing the simulations of the average $X_\text{max}$ to observations, new Earth-based bounds on this type of Lorentz violation are obtained, which are better than the previous bounds by more than three orders of magnitude. Prospects of further studies are also discussed.

Axion-like particles and recent observations of the cosmic infrared background radiation

The CIBER collaboration released their first observational data of the Cosmic IR background (CIB) radiation, which has significant excesses at around the wavelength $\sim$ 1 $\mu$m compared to theoretically-inferred values. The amount of the CIB radiation has a significant influence on the opaqueness of the Universe for TeV gamma-rays emitted from distant sources such as AGNs. With the value of CIB radiation reported by the CIBER experiment, through the reaction of such TeV gamma-rays with the CIB photons, the TeV gamma-rays should be significantly attenuated during propagation, which would lead to energy spectra in disagreement with current observations of TeV gamma ray sources. In this article, we discuss a possible resolution of this tension between the TeV gamma-ray observations and the CIB data in terms of axion [or Axion-Like Particles (ALPs)] that may increase the transparency of the Universe by the anomaly-induced photon-axion mixing. We find a region in the parameter space of the axion mass, $m_a \sim 7 \times 10^{-10} - 5 \times 10^{-8}$eV, and the axion-photon coupling constant, $1.2 \times 10^{-11} {\rm GeV}^{-1} \lesssim g_{a\gamma} \lesssim 8.8 \times 10^{-10} {\rm GeV}^{-1}$ that solves this problem.

The ARGO-YBJ legacy to next generation wide field-of-view experiments

The ARGO-YBJ experiment has been in stable data taking for more than 5 years at the YangBaJing Cosmic Ray Observatory (Tibet, P.R. China, 4300 m a.s.l., 606 g/cm 2 ). With a duty-cycle greater than 86% the detector collected about 5×10 11 events in a wide energy range, from few hundreds GeV up to about 10 PeV. High altitude location and detector features make ARGO-YBJ capable of investigating a wide range of important issues in Cosmic Ray and Astroparticle Physics by imaging the front of atmospheric showers with unprecedented resolution and detail. In this contribution some of the latest physics results obtained by ARGO-YBJ in gamma-ray astronomy and in cosmic ray physics are summarized. The prospects of TeV gamma-ray observations with new ground-based wide field-of-view detectors are presented.

NONTHERMAL RADIATION OF YOUNG SUPERNOVA REMNANTS: THE CASE OF CAS A

The processes responsible for the broad-band radiation of the young supernova remnant Cas A are explored using a new code which is designed for a detailed treatment of the diffusive shock acceleration of particles in nonlinear regime. The model is based on spherically symmetric hydrodynamic equations complemented with transport equations for relativistic particles. Electrons, protons and the oxygen ions accelerated by forward and reverse shocks are included in the numerical calculations. We show that the available multi-wavelength observations in the radio, X-ray and gamma-ray bands can be best explained by invoking particle acceleration by both forward and reversed shocks. Although the TeV gamma-ray observations can be interpreted by interactions of both accelerated electrons and protons/ions, the measurements by Fermi LAT at energies below 1 GeV give a tentative preference to the hadronic origin of gamma-rays. Then, the acceleration efficiency in this source, despite the previous claims, should be very high; 25% of the explosion energy (or approximately $3\cdot 10^{50}$ erg) should already be converted to cosmic rays, mainly by the forward shock. At the same time, the model calculations do not provide extension of the maximum energy of accelerated protons beyond 100 TeV. In this model, the acceleration of electrons is dominated by the reverse shock; the required $10^{48}$ erg can be achieved under the assumption that the injection of electrons (positrons) is supported by the radioactive decay of $^{44}$Ti.

VERITAS OBSERVATIONS OF THE MICROQUASAR CYGNUS X-3

We report results from TeV gamma-ray observations of the microquasar Cygnus X-3. The observations were made with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) over a time period from 2007 June 11 to 2011 November 28. VERITAS is most sensitive to gamma rays at energies between 85 GeV to 30 TeV. The effective exposure time amounts to a total of about 44 hours, with the observations covering six distinct radio/X-ray states of the object. No significant TeV gamma-ray emission was detected in any of the states, nor with all observations combined. The lack of a positive signal, especially in the states where GeV gamma rays were detected, places constraints on TeV gamma-ray production in Cygnus X-3. We discuss the implications of the results.

Nonthermal emission properties of the northwestern rim of supernova remnant RX J0852.0-4622

The supernova remnant (SNR) RX J0852-4622 (Vela Jr., G266.6-1.2) is one of the most important SNRs for investigating the acceleration of multi-TeV particles and the origin of Galactic cosmic rays because of its strong synchrotron X-ray and TeV gamma-ray emission, which show a shell-like morphology similar to each other. Using the XMM-Newton archival data consisting of multiple pointing observations of the northwestern rim of the remnant, we investigate the spatial properties of the nonthermal X-ray emission as a function of distance from an outer shock wave. All X-ray spectra are well reproduced by an absorbed power-law model above 2 keV. It is found that the spectra show gradual softening from a photon index 2.56 in the rim region to 2.96 in the interior region. We show that this radial profile can be interpreted as a gradual decrease of the cutoff energy of the electron spectrum due to synchrotron cooling. By using a simple spectral evolution model that includes continuous synchrotron losses, the spectral softening can be reproduced with the magnetic field strength in the post-shock flow to less than several tens of uG. If this is a typical magnetic field in the SNR shell, gamma-ray emission would be accounted for by inverse Compton scattering of high-energy electrons that also produce the synchrotron X-ray emission. Future hard X-ray imaging observations with Nustar and ASTRO-H and TeV gamma-ray observations with the Cherenkov Telescope Array (CTA) will allow to us to explore other possible explanations of the systematic softening of the X-ray spectra.

Very-high energy gamma-ray astronomy

Very-high energy (VHE) gamma quanta contribute only a minuscule fraction - below one per million - to the flux of cosmic rays. Nevertheless, being neutral particles they are currently the best "messengers" of processes from the relativistic/ultra-relativistic Universe because they can be extrapolated back to their origin. The window of VHE gamma rays was opened only in 1989 by the Whipple collaboration, reporting the observation of TeV gamma rays from the Crab nebula. After a slow start, this new field of research is now rapidly expanding with the discovery of more than 150 VHE gamma-ray emitting sources. Progress is intimately related with the steady improvement of detectors and rapidly increasing computing power. We give an overview of the early attempts before and around 1989 and the progress after the pioneering work of the Whipple collaboration. The main focus of this article is on the development of experimental techniques for Earth-bound gamma-ray detectors; consequently, more emphasis is given to those experiments that made an initial breakthrough rather than to the successors which often had and have a similar (sometimes even higher) scientific output as the pioneering experiments. The considered energy threshold is about 30 GeV. At lower energies, observations can presently only be performed with balloon or satellite-borne detectors. Irrespective of the stormy experimental progress, the success story could not have been called a success story without a broad scientific output. Therefore we conclude this article with a summary of the scientific rationales and main results achieved over the last two decades.

TeV gamma-ray observations of Markarian 421 using TACTIC during 2009–10

We have observed the blazar Markarian 421 with the TACTIC γ-ray telescope at Mt. Abu, India, from 22 November 2009 to 16 May 2010 for 265 h. Detailed analysis of the data so recorded revealed the presence of a TeV γ-ray signal with a statistical significance of 12.12σ at Eγ ⩾ 1 TeV. We have estimated the time-averaged differential energy spectrum of the source in the energy range 1.0–16.44 TeV. The spectrum fits well with the power-law function of the form (dF/dE = f0E−Γ) with f0 = (1.39 ± 0.239) × 10−11 cm−2 s−1 TeV−1 and Γ = 2.31 ± 0.14.

OBSERVATION OF TeV GAMMA RAYS FROM THE CYGNUS REGION WITH THE ARGO-YBJ EXPERIMENT

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

TeV AND MULTI-WAVELENGTH OBSERVATIONS OF Mrk 421 IN 2006-2008

We report on TeV gamma-ray observations of the blazar Mrk 421 (redshift of 0.031) with the VERITAS observatory and the Whipple 10m Cherenkov telescope. The excellent sensitivity of VERITAS allowed us to sample the TeV gamma-ray fluxes and energy spectra with unprecedented accuracy where Mrk 421 was detected in each of the pointings. A total of 47.3 hrs of VERITAS and 96 hrs of Whipple 10m data were acquired between January 2006 and June 2008. We present the results of a study of the TeV gamma-ray energy spectra as a function of time, and for different flux levels. On May 2nd and 3rd, 2008, bright TeV gamma-ray flares were detected with fluxes reaching the level of 10 Crab. The TeV gamma-ray data were complemented with radio, optical, and X-ray observations, with flux variability found in all bands except for the radio waveband. The combination of the RXTE and Swift X-ray data reveal spectral hardening with increasing flux levels, often correlated with an increase of the source activity in TeV gamma-rays. Contemporaneous spectral energy distributions were generated for 18 nights, each of which are reasonably described by a one-zone SSC model.

12 mm line survey of the dense molecular gas towards the W28 field TeV gamma-ray sources

We present 12mm Mopra observations of dense molecular gas towards the W28 supernova remnant (SNR) field. The focus is on the dense molecular gas towards the TeV gamma-ray sources detected by the H.E.S.S. telescopes, which likely trace the cosmic-rays from W28 and possibly other sources in the region. Using the NH3 inversion transitions we reveal several dense cores inside the molecular clouds, the majority of which coincide with high-mass star formation and HII regions, including the energetic ultra-compact HII region G5.89-0.39. A key exception to this is the cloud north east of W28, which is well-known to be disrupted as evidenced by clusters of 1720MHz OH masers and broad CO line emission. Here we detect broad NH3, up to the (9,9) transition, with linewidths up to 16 km/s. This broad NH3 emission spatially matches well with the TeV source HESS J1801-233 and CO emission, and its velocity dispersion distribution suggests external disruption from the W28 SNR direction. Other lines are detected, such as HC3N and HC5N, H2O masers, and many radio recombination lines, all of which are primarily found towards the southern high-mass star formation regions. These observations provide a new view onto the internal structures and dynamics of the dense molecular gas towards the W28 SNR field, and in tandem with future higher resolution TeV gamma-ray observations will offer the chance to probe the transport of cosmic-rays into molecular clouds.

MAGIC TeV gamma-ray observations of Markarian 421 during multiwavelength campaigns in 2006

The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of 2006 April 28, 29, and 2006 June 14. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from 2006 April 22 to 30 and on 2006 June 14. We inferred light curves with sub-day resolution and night-by-night energy spectra. A strong gamma-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis between (0.92+-0.11)10^-10 cm^-2 s^-1 (0.57 Crab units) and (3.21+-0.15)10^-10 cm^-2 s^-1 (2.0 Crab units) in 2006 April. There is a clear indication for intra-night variability with a doubling time of 36+-10(stat) minutes on the night of 2006 April 29, establishing once more rapid flux variability for this object. For all individual nights gamma-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing flux, as expected in leptonic acceleration models.