Gamma-ray Observations Research Articles

Testing Lorentz invariance of electrons with LHAASO observations of PeV gamma-rays from the Crab Nebula

The Large High Altitude Air Shower Observatory (LHAASO) recently reported the detection of gamma-ray emissions with energies up to 1.1PeV from the Crab Nebula. Using the absence of vacuum Cherenkov effect by inverse-Compton electrons, we improve previous bounds to linear-order Lorentz invariance violation (LV) in the dispersion relations of electrons by 104 times. We show that the LV effect on electrons is severely constrained, compatible with certain type of LV as expected by some models of quantum gravity (QG), such as the string/D-brane inspired space-time foam. We argue that such models are supported by the Crab Nebula constraints from the LHAASO observations, as well as various LV phenomenologies for photons to date.

Measurement of Galactic 26Al with the Compton Spectrometer and Imager

Abstract The Compton Spectrometer and Imager (COSI) is a balloon-borne compact Compton telescope designed to survey the 0.2–5 MeV sky. COSI’s energy resolution of ∼0.2% at 1.8 MeV, single-photon reconstruction, and wide field of view make it capable of studying astrophysical nuclear lines, particularly the 1809 keV γ-ray line from decaying Galactic 26Al. Most 26Al originates in massive stars and core-collapse supernova nucleosynthesis, but the path from stellar evolution models to Galaxy-wide emission remains unconstrained. In 2016, COSI had a successful 46 day flight on a NASA superpressure balloon. Here, we detail the first search for the 1809 keV 26Al line in the COSI 2016 balloon flight using a maximum-likelihood analysis. We find a Galactic 26Al flux of (8.6 ± 2.5) × 10−4 ph cm−2 s−1 within the Inner Galaxy (∣ℓ∣ ≤ 30°, ∣b∣ ≤ 10°) with 3.7σ significance above background. Within uncertainties, this flux is consistent with expectations from previous measurements by SPectrometer on INTEGRAL (SPI) and the Compton Telescope on the Compton Gamma-Ray Observatory (COMPTEL). This analysis demonstrates COSI’s powerful capabilities for studies of γ-ray lines and underscores the scientific potential of future compact Compton telescopes. In particular, the next iteration of COSI as a NASA Small Explorer satellite has recently been approved for launch in 2025.

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.

Non-destructive 3D imaging method using muonic X-rays and a CdTe double-sided strip detector

Elemental analysis based on muonic X-rays resulting from muon irradiation provides information about bulk material composition without causing damage, which is essential in the case of precious or otherwise unreachable samples, such as in archeology and planetary science. We developed a three-dimensional (3D) elemental analysis technique by combining the elemental analysis method based on negative muons with an imaging cadmium telluride double-sided strip detector (CdTe-DSD) designed for the hard X-ray and soft gamma-ray observation. A muon irradiation experiment using spherical plastic samples was conducted at the Japan Proton Accelerator Research Complex (J-PARC); a set of projection images was taken by the CdTe-DSD, equipped with a pinhole collimator, for different sample rotation angles. The projection images measured by the CdTe-DSD were utilized to obtain a 3D volumetric phantom by using the maximum likelihood expectation maximization algorithm. The reconstructed phantom successfully revealed the 3D distribution of carbon in the bulk samples and the stopping depth of the muons. This result demonstrated the feasibility of the proposed non-destructive 3D elemental analysis method for bulk material analysis based on muonic X-rays.

Physical properties of the molecular cloud, N4, in SS 433: Evidence for an interaction of molecular cloud with the jet from SS 433

Abstract We conducted observations and analyses of the molecular cloud, N4, which is located at ∼40 pc from SS 433 and the same line of sight as that of the radio shell, in 12CO(J = 1–0), 12CO(J = 3–2), 13CO(J = 3–2), and grand-state OH emissions. N4 has a strong gradient of the integrated intensity of 12CO(J = 1–0, 3–2) emission at the northern, eastern, and western edges. The main body of N4 also has a velocity gradient of ∼0.16 km s−1 (20″)−1. A velocity shift by up to 3 km s−1 from the systemic velocity at ∼49 km s−1 is detected at only the northwestern part of N4. The volume density of the molecular hydrogen gas and the kinematic temperature are estimated at eight local peaks of 12CO(J = 1–0) and 13CO(J = 3–2) emissions by the RADEX code. The calculated $n_{\rm (H_2)}$ is an order of 103 cm−3, and Tk ranges from ∼20 to ∼56 K. The mass of N4 is estimated to be ∼7300 M⊙. The thermal and turbulent pressures in N4 are estimated to be ∼105 K cm−3 and ∼107 K cm−3, respectively. The relation of the thermal and turbulent pressures in N4 tends to be similar to that of the molecular clouds in the Galactic plane. However, these values are higher than those in the typical molecular clouds in the Galactic plane. Several pieces of circumstantial evidence representing the physical properties of N4 and comparison with the data of infrared and X-ray radiation suggest that N4 is interacting with a jet from SS 433. However, no gamma-ray radiation is detected toward N4. Compared to the previous study, it is hard to detect the gamma-ray radiation by cosmic-ray proton origin due to the low sensitivity of the current gamma-ray observatories. No OH emission was detected toward N4 due to the low sensitivity of the observation and antenna beam dilution.

The RTApipe framework for the gamma-ray real-time analysis software development

In the multi-messenger era, coordinating observations between astronomical facilities is mandatory to study transient phenomena (e.g. Gamma-ray bursts) and is achieved by sharing information with the scientific community through networks such as the Gamma-ray Coordinates Network. The facilities usually develop real-time scientific analysis pipelines to detect transient events, alert the astrophysical community, and speed up the reaction time of science alerts received from other observatories. We present in this work the RTApipe framework, designed to facilitate the development of real-time scientific analysis pipelines for present and future gamma-ray observatories. This framework provides pipeline architecture and automatisms, allowing the researchers to focus on the scientific aspects and integrate existing science tools developed with different technologies. The pipelines automatically execute all the configured analyses during the data acquisition. This framework can be interfaced with science alerts networks to perform follow-up analysis of transient events shared by other facilities. The analyses are performed in parallel and can be prioritised. The workload is highly scalable on a cluster of machines. The framework provides the required services using containerisation technology for easy deployment. We present the RTA pipelines developed for the AGILE space mission and the prototype of the SAG system for the ground-based future Cherenkov Telescope Array observatory confirming that the RTApipe framework can be used to successfully develop pipelines for the gamma-ray observatories, both space and ground-based.

First observations of solar disk gamma rays over a full solar cycle

The solar disk is among the brightest $\ensuremath{\gamma}$-ray sources in the sky. It is also among the most mysterious. No existing model fully explains the luminosity, spectrum, time variability, and morphology of its emission. We perform the first analysis of solar disk $\ensuremath{\gamma}$-rays over a full 11 year solar cycle, utilizing a powerful new method to differentiate solar signals from astrophysical backgrounds. We produce: (i) a robustly measured spectrum from 100 MeV to 100 GeV, reaching a precision of several percent in the 1--10 GeV range, (ii) new results on the anticorrelation between solar activity and $\ensuremath{\gamma}$-ray emission, (iii) strong constraints on short-timescale variability, ranging from hours to years, and (iv) new detections of the equatorial and polar morphologies of high-energy $\ensuremath{\gamma}$-rays. We find no significant energy dependence in the time variability of solar disk emission, indicating that strong magnetic-field effects close to the solar surface, rather than modulation throughout the heliosphere, must primarily control the flux and morphology of solar disk emission.

A minimal model for understanding secondary cosmic rays

We take a phenomenological approach in a minimal model to understand the spectral intensity of secondary cosmic-ray particles like positrons, antiprotons, Lithium, Beryllium and Boron. Our analysis shows that cosmic rays at ∼ GeV energies pass through a significant amount of matter in regions surrounding the sources. This grammage decreases with increasing cosmic-ray energy and becomes negligible beyond ∼100 GeV. During the subsequent propagation in the interstellar medium cosmic rays of all energies up to ∼105 GeV/n pass through about 1–2 g cm−2 of matter before leaking into the intergalactic medium. It is in the interstellar medium that the bulk of the positrons and antiprotons are generated. Also cosmic-ray nuclei like C, N, and O at all energies generate additional amounts of Li, Be and B nuclei with a spectrum similar to those of C, O etc. The implications of these findings of the minimal model to the observations of gamma rays and also the importance of spatial and temporal discreteness of cosmic-ray sources for modeling cosmic-ray propagation are briefly pointed out.

A two-zone model as origin of hard TeV spectrum in extreme BL lacs

ABSTRACT The emission of the so-called extreme BL Lacs poses challenges to the particle acceleration models. The hardness of their spectrum, ≲ 2, in the high-energy band demands unusual parameters using the standard one-zone synchrotron self-Compton (SSC) model with a deficient magnetized plasma. Some authors use either two-zone or hadronic/lepto-hadronic models to relax these atypical values. In this work, we present a lepto-hadronic two-zone model to explain the multi-wavelength observations of the six best-known extreme BL Lacs. The very-high-energy gamma-ray observations are described by the photo-hadronic processes in a blob close to the AGN core and by SSC and external inverse Compton-processes in an outer blob. The photo-hadronic interactions occur when accelerated protons in the inner blob interact with annihilation line photons from a sub-relativistic pair plasma. The X-ray observations are described by synchrotron radiation from the outer blob. The parameter values found from the description of the spectral energy distribution for each object with our phenomenological model are similar to each other, and lie in the typical range reported in BL Lacs.

Constraining blazar heating with the 2 ≲z≲ 3 Lyman-α forest

ABSTRACTThe intergalactic medium (IGM) acts like a calorimeter recording energy injection by cosmic structure formation, shocks and photoheating from stars and active galactic nuclei. It was recently proposed that spatially inhomogeneous TeV-blazars could significantly heat up the underdense IGM, resulting in patches of both cold and warm IGM around z ≃ 2−3. The goal of this study is to compare predictions of different blazar heating models with recent observations of the IGM. We perform a set of cosmological simulations and carefully compute mock observables of the Lyman-α (Lyα) forest. We perform a detailed assessment of different systematic uncertainties which typically impact this type of observables and find that they are smaller than the differences between our models. We find that our inhomogeneous blazar heating model is in good agreement with the Lyα line properties and the rescaled flux probability distribution function at high redshift (2.5 < z < 3) but that our blazar heating models are challenged by lower redshift data (2 < z < 2.5). Our results could be explained by HeII reionization although state-of-the-art models fall short on providing enough heating to the low-density IGM, thus motivating further radiative transfer studies of inhomogeneous HeII reionization. If blazars are indeed hosted by group-mass haloes of $2\times 10^{13}\, \rm {M}_\odot$, a later onset of blazar heating in comparison with previous models would be favoured, which could bring our findings here in agreement with the evidence of blazar heating from local gamma-ray observations.

Constraints on Cosmic Rays Population in the Radio Halo of the M87 Galaxy from Gamma-Ray Observations

We present results of our search for extended high-energy emission around the M87 galaxy, which is the central galaxy in the Virgo cluster. The M87 galaxy harbors extremely heavy supermassive black hole, {{M}_{{BH}}} > 4 times {{10}^{9}}{kern 1pt} {{M}_{ odot }} and for the last 100 Myr experiences a period of increased activity that resulted in formation of extended lobes, which are prominent radio- and X-ray sources. We perform a search in 13 yr of the Fermi-LAT data in three energy bins: 0.1–1, 1–10, and 10–100 GeV. No significant detection of extended emission was achieved in any bin, with weak indication of a signal ( sim {kern 1pt} 2.5sigma ) in the last bin (10–100 GeV). We used γ-ray observations to put limits on cosmic rays population in the lobes, constraining their fraction X in comparison to the thermal energy of the hot gas. Depending on the value of spectral index, for electrons {{X}_{e}} < (0.1{-} 0.2) and protons were constrained at {{X}_{p}} < (0.2{-} 0.5) level. We confirm that cosmic rays is a subdominant component of the lobes, which are dominated by thermal gas.

The Composition and Power of the Jet of the Broad-line Radio Galaxy 3C 120

Abstract We calculated the electron–positron pair-production rate at the base of the jet of 3C 120 due to collisions of photons from the hot accretion flow using the measurement of its average soft gamma-ray spectrum by the Compton Gamma Ray Observatory. We found that this rate approximately equals the flow rate of leptons emitting the observed synchrotron radio-to-IR spectrum of the jet core, calculated using the extended jet model following Blandford &amp; Königl. This coincidence shows the jet composition is likely to be pair dominated. We then calculated the jet power in the bulk motion of ions and found it greatly exceeds that achievable by the magnetically arrested disk scenario for the maximum black hole spin unless the jet contains mostly pairs. Next, we found that the magnetic flux through the synchrotron-emitting jet equals the maximum poloidal flux that can thread the black hole. Finally, we compared two estimates of the magnetization parameter at the onset of the synchrotron emission and found they are in agreement only if pairs dominate the jet content.

Implications of gamma-ray photon measurements at the LHAASO on Lorentz symmetry * *Supported in part by the Natural Science Foundation of China (11875053, 12035016). Liu Yuxuan acknowledges the support from the National Postdoctoral Program for Innovative Talents (BX2021303), funded by China Postdoctoral Science Foundation

The Large High Altitude Air Shower Observatory (LHAASO) has reported the measurement of photons with high energies of up to 1.42 PeV from twelve gamma-ray sources. We are concerned with the implications of the LHAASO data on the fate of Lorenz symmetry at such high energy levels; thus, we consider the interaction between gamma rays and photons in the cosmic microwave background (CMB) and compute the optical depth, mean free path, and survival probability of photons from these gamma-ray sources. Employing the threshold value predicted by standard special relativity, the lowest survival probability for observed gamma ray photons is found to be approximately 0.60, which is fairly high and implies that abundant photons with energies above the threshold may reach the Earth without Lorentz symmetry violation. We conclude that it is unreasonable to argue that Lorentz symmetry would be violated using current observations at the LHAASO.

Calculated Event Rates for Axion Detection via Atomic and Nuclear Processes

The possibility of detection of 5.5 MeV and 14.4 keV solar axions by observing axion-induced nuclear and atomic transitions is investigated. The presence of nuclear transitions between spin orbit partners can be manifested by the subsequent deexcitation via gamma ray emissions. The transition rates can also be studied in the context of radiative axion absorption by a nucleus. The elementary interaction is obtained in the context of the axion-quark couplings predicted by existing axion models. Then, these couplings will be transformed to the nucleon level utilizing reasonable existing models, which lead to effective transition operators. Using these operators, we calculate the needed nuclear matrix elements employing wave functions obtained in the context of the nuclear shell model. With these ingredients, we discuss possibilities of experimental observation of the axion-induced nuclear gamma rays. In the second part, we will examine the axion-induced production of X-rays (axion-photon conversion) or ionization from deeply bound electron orbits. In this case, the axion electron coupling is predicted by existing axion models; no renormalization is needed. The experimental signal is the observation of directly produced electrons and/or the emission of hard X-rays and Auger electrons, following the deexcitation of the final atom. Critical discussion is made on the experimental feasibility of detecting the solar axions by using multiton scale NaI detectors.

Gamma-Ray Cosmology and Tests of Fundamental Physics

The propagation of gamma-rays over cosmological distances is the subject of extensive theoretical and observational research at GeV and TeV energies. The mean free path of gamma-rays in the cosmic web is limited above 100 GeV due to the production of electrons and positrons on the cosmic optical and infrared backgrounds. Electrons and positrons cool in the intergalactic medium while gyrating in its magnetic fields, which could cause either its global heating or the production of lower-energy secondary gamma-rays. The energy distribution of gamma-rays surviving the cosmological journey carries observed absorption features that gauge the emissivity of baryonic matter over cosmic time, constrain the distance scale of ΛCDM cosmology, and limit the alterations of the interaction cross section. Competitive constraints are, in particular, placed on the cosmic star-formation history as well as on phenomena expected from quantum gravity and string theory, such as the coupling to hypothetical axion-like particles or the violation of Lorentz invariance. Recent theoretical and observational advances offer a glimpse of the multi-wavelength and multi-messenger path that the new generation of gamma-ray observatories is about to open.

Formula omitted]: A new variable for γ/h discrimination in large gamma-ray ground arrays

In this letter, a new strategy to enhance the discrimination of high energy gamma rays from the huge charged cosmic rays background in large cosmic rays ground arrays is presented. This strategy is based on the introduction of a new simple variable, Pγhα, which combines the probability of tagging muons and/or very energetic particles in each single array station. The discrimination power of this new variable, particularly important for and above multi-TeV energies, is illustrated for a few specific examples in the case of a hypothetical water Cherenkov detector cosmic ray array, both in the case of low and high particle stations occupancy. The results are very encouraging and hopefully will be demonstrated in the present and future gamma-ray Observatories.

Constraints on ultracompact minihalos from the extragalactic gamma-ray background observation

Ultracompact minihalo (UCMH) is a special type of dark matter halo with a very steep density profile that may form in the early Universe seeded by an overdense region or a primordial black hole. Constraints on its abundance give valuable information on the power spectrum of primordial perturbation. In this work, we update the constraints on the UCMH abundance in the Universe using the extragalactic gamma-ray background observation. Comparing to previous works, we adopt the updated Fermi-LAT extragalactic gamma-ray background measurement and derive constraints based on a full consideration of the astrophysical contributions. With these improvements, we place constraints on UCMH abundance 1--2 orders of magnitude better than previous results. With the background components considered, we can also attempt to search for possible additional components beyond the known astrophysical contributions.

Deep optical study of the mixed-morphology supernova remnant G 132.7+1.3 (HB3)

ABSTRACT We present optical CCD images of the large supernova remnant (SNR) G 132.7+1.3 (HB3) covering its full extent for the first time, in the emission lines of Hα +[N ii], [S ii], and [O iii], where new and known filamentary and diffuse structures are detected. These observations are supplemented by new low-resolution long-slit spectra and higher resolution images in the same emission lines. Both the flux-calibrated images and spectra confirm that the optical emission originates from shock-heated gas since the [S ii]/Hα &amp;gt; 0.4. Our findings are also consistent with the recently developed emission-line ratio diagnostics for distinguishing SNRs from H ii regions. A multiwavelength comparison among our optical data and relevant observations in radio, X-rays, gamma-rays and CO bands, provided additional evidence on the interaction of HB3 with the surrounding clouds and clarified the borders of the SNR and the adjacent cloud. We discuss the supernova (SN) properties and evolution that led to the current observables of HB3 and we show that the remnant has most likely passed at the pressure driven snowplow phase. The estimated SN energy was found to be (3.7 ± 1.5) × 1051 erg and the current SNR age (5.1 ± 2.1) × 104 yr. We present an alternative scenario according to which the SNR evolved in the wind bubble cavity excavated by the progenitor star and currently is interacting with its density walls. We show that the overall mixed morphology properties of HB3 can be explained if the SN resulted by a Wolf−Rayet progenitor star with mass $\sim \! 34 \rm ~M_{\odot }$.

Atmospheric Electron Spatial Range Extended by Thundercloud Electric Field Below the Relativistic Runaway Electron Avalanche Threshold

AbstractCosmic‐ray shower, hitting the atmosphere all the time, includes high‐energy electrons. Electric fields in thunderstorms accelerate these electrons so that they can reach farther distances with producing secondary electrons. The strong field above 0.284 MV m−1 causes the Relativistic Runaway Electron Avalanches (RREA) process with an exponential increase of electrons. Even with lower electric fields below this RREA threshold which are usually measured in thunderstorms, an increase in electrons' kinetic energy and production of secondary electrons can occur without avalanche processes, known as Modification Of Spectra (MOS) of the cosmic‐ray shower. Both phenomena are related to observed gamma ray glow, which is bremsstrahlung emission from accelerated electrons in thunderstorms. We calculate an extended range and kinetic energy of electrons below the RREA threshold using an analytical evaluation and Geant4 Monte Carlo simulations. For example, at 0.280 MV m−1 slightly below the RREA threshold, 3 MeV electrons gain 100% of their original kinetic energy through the passage of the continuous slowing down approximation (CSDA) range defined at the null electric field, which is consistent with results in previous studies as the transition to the RREA process. Even at a lower field of 0.260 MV m−1, the energy recovery by the field allows an injected 20 MeV electron beam to extend its average electron range from 73 m, in the null field configuration, to 562 m. We have performed analytical calculation and Geant4 simulations of atmospheric electron behavior for initial electron energy of 0.1–100 MeV and an electric field below 0.290 MV m−1.

Consistent explanation for the cosmic-ray positron excess in p-wave Breit\u2013Wigner enhanced dark matter annihilation

Dark matter (DM) annihilation in the galactic halo can be enhanced relative to that in the early Universe due to the Breit–Wigner enhancement, if the DM particles annihilate through a narrow resonance. Although the s-wave Breit–Wigner enhancement can provide a consistent explanation for both the observed cosmic-ray (CR) positron excess and the DM thermal relic density, it is severely constrained by the observations of gamma rays from dwarf spheroidal satellite galaxies (dSphs) and the cosmic microwave background (CMB), which have relatively lower allowed DM annihilation cross section and typical DM velocities than that in the galactic halo. Furthermore, in the s-wave Breit–Wigner enhancement, the case where the resonance mass is below a threshold (twice the DM mass) is ruled out due to the monotonically increasing annihilation cross section with decreasing DM velocity. In this work, we consider Breit–Wigner enhanced p-wave DM annihilation. We explore the parameter regions which can simultaneously account for the CR positron excess and DM thermal relic density without violating the constraints from dSphs gamma rays and CMB. We show that the velocity-dependent cross section in this scenario can peak around the typical DM velocity in the galactic halo for the resonance mass both above and below the threshold. Moreover, the highly suppressed annihilation cross section at extremely low DM velocity can evade the constraints from dSphs gamma rays and CMB easily, which results in larger allowed parameter regions than that in the s-wave case.