Resolving the origin of the unidentified TeV source HESS J1626-490 as a relic of the ancient cosmic-ray factory SNR G335.2+0.1

arXiv:2507.04407 · astro-ph.HE · Submitted 2025-07-06 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Hi, I'm Vera.

Jocelyn: And I'm Jocelyn, here with Subrahmanyan.

Vera: We're looking at the paper: "Resolving the origin of the unidentified TeV source HESS J1626-490 as a relic of the ancient cosmic-ray factory SNR G335.2+0.1".

Jocelyn: Here's the gist of the paper. The origin of HESS J1626−490 as a relic of SNR G335.2+0.1 is investigated by analyzing multi-wavelength data to test the CR-escape scenario, which suggests that protons accelerated at the SNR escape and illuminate nearby clouds to produce gamma rays Observational Findings The study scrutinized properties of SNR G335.2+0.1 and searched for non-thermal counterparts in radio, X-ray, and GeV gamma-ray bands The distance and age of the SNR are estimated to be 3.3 ± 0.6 kpc and ∼5 kyr X-ray and GeV gamma-ray data analyses revealed an extended GeV gamma-ray emission overlapping with the SNR and H.E.S.S.

Vera: source regions but no X-ray emission The location of the peak of the extended GeV emission changes from near the SNR at ≲ 1 GeV to the H.E.S.S.

Jocelyn: source at > 10 GeV 10 GeV> GeV Gamma-Ray Morphology Fermi-LAT data analysis revealed an apparent excess in a large area containing both the SNR and H.E.S.S.

Vera: source regions The GeV morphology was evaluated by fitting it with a symmetric two-dimensional (2-d) Gaussian function The location of the GeV emission gradually changes with energy, being in good agreement with the SNR in the lower energy band below 2 GeV and consistent with the H.E.S.S.

Title and authors: Vera: Now we're really digging into the specifics of "Resolving the origin of the unidentified TeV source HESS J1626-four hundred ninety as a relic of the ancient cosmic-ray factory SNR G335 point 2+zero point one" <ref:2507.04407#pg1,Resolving the origin of the unidentified TeV source HESS J1626-490 as>. We want to unpack what this title actually means for us beyond just saying it's a leftover source.

Jocelyn: It means they are trying to pinpoint the exact mechanism by which these TeV gamma rays are produced, and whether that mechanism is tied directly to the specific environment of SNR G335 point 2+zero point one <ref:2507.04407#pg1,SNR G335.2+0.1>.

Subrahmanyan: Basically, it’s about testing the CR-escape scenario—the idea that protons get accelerated at an SNR, escape out into space, and then illuminate nearby clouds to make those gamma rays via pi zero decay.

Vera: So, the core question is whether this proton escape and illumination model is physically consistent with all the multi-wavelength data we collected.

Jocelyn: They are testing if the spatial distribution of that cosmic ray interaction matches what we see in the GeV emission maps across different energy bands.

Subrahmanyan: If they can show that the GeV morphology changes with energy exactly as described, it strongly supports a scenario where you have distinct sources contributing at different energies.

Vera: The paper also focuses on modeling how those emissions combine to form the final spectrum, trying to separate what comes from the SNR and what comes from nearby gas clouds.

Jocelyn: And they are using that two-Gaussian model to see if it can actually decompose the observed spectrum into these two distinct parts, which is a way of quantifying that spectral hardening we talked about earlier <ref:2507.04407#pg1>.

Subrahmanyan: Separating the SNR contribution from the cloud contribution requires very precise modeling because those environments are spatially close but physically different, which is a challenge for any theory.

Vera: The improvement they suggest here involves incorporating more detailed CR transport models into their analysis, specifically using solutions to equation B4 with diffusion coefficients that depend on magnetic field strength and turbulence, as described by equations B5 and B13.

Jocelyn: That's a big step because it allows them to simulate a much more realistic spatial distribution for both electrons and protons in the cloud region, which is what we need to get accurate results.

Subrahmanyan: Modeling diffusion coefficients that are functions of magnetic field strength and turbulence means we’re moving away from assuming constant diffusion, which is a major refinement for understanding how cosmic rays move.

The paper's summary: Vera: So, let's look at the actual summary of "Resolving the origin of the unidentified TeV source HESS J1626-four hundred ninety as a relic of the ancient cosmic-ray factory SNR G335 point 2+zero point one" <ref:2507.04407#pg1,Resolving the origin of the unidentified TeV source HESS J1626-490 as>.

Jocelyn: The summary basically says they used multi-wavelength data, including radio continuum maps, X-ray observations from XMM-Newton, and GeV gamma-ray data from Fermi-LAT to constrain this source.

Subrahmanyan: They also found that they estimated the distance and age of SNR G335 point 2+zero point one to be three point three plus or minus zero point six kiloparsecs and about five thousand years old, which gives us a good context for the potential source's origin <ref:2507.04407#pg1,SNR G335.2+0.1>.

Vera: The summary highlights the observations showing an extended GeV gamma-ray emission overlapping with both the SNR and HESS J1626−four hundred ninety regions, but no X-ray emission near the SNR <ref:2507.04407#pg1>.

Jocelyn: And they mention that they looked at how the peak of that extended GeV emission changes with energy, shifting from near the SNR at less than one GeV to being at the HESS J1626−four hundred ninety source when we look above ten GeV <ref:2507.04407#pg1>.

Subrahmanyan: This observation about the energy dependence is a key finding because it links the spatial distribution of gamma rays to different physical scales, which is what we need for structure formation studies.

Vera: The Fermi-LAT data analysis revealed an apparent excess in a large area containing both the SNR and HESS J1626−four hundred ninety regions, and when they fit that GeV morphology with a symmetric two-dimensional Gaussian function across different energy bands <ref:2507.04407#pg1>.

Jocelyn: They found that the location of this GeV emission gradually changes with energy, being in good agreement with the SNR in the lower energy band below two GeV and consistent with the HESS J1626−four hundred ninety source when above eleven GeV <ref:2507.04407#pg1>.

Subrahmanyan: That spectral tracking across energy bands is what helps confirm that we have different physical components contributing to this signal at different energies, which is a strong structural result.

Vera: The observed spectrum of the Gaussian-W shows some spectral hardening at fifty GeV and it smoothly connects to the HESS spectrum, suggesting two components there, like those from the SNR and a molecular cloud.

Jocelyn: To model this broadband spectrum, they found that a combination of emission from escaped cosmic rays and the SNR itself was consistent with observations, assuming protons were accelerated to more than one hundred TeV in the past <ref:2507.04407#pg3>.

Subrahmanyan: So this is where we link the acceleration to the observed gamma-ray shape, suggesting high energy particle acceleration in these old remnants is necessary for this picture.

Vera: They explored both hadronic-dominated and leptonic-dominated models to explain the TeV gamma-ray spectrum, but they chose the hadronic model because it avoided an "unrealistic assumption about the time when electrons escape from the SNR" that some of those leptonic models required.

Jocelyn: The preference for the hadronic model comes with a caveat about how slow cosmic rays are propagating around this remnant, which means we need to be careful about our assumptions on magnetic fields.

Subrahmanyan: That low diffusion coefficient they found is what pushes us toward the hadronic model because it suggests that standard isotropic diffusion models might not describe this region accurately at all.

The paper's improvements: Vera: Now let's talk about the specific improvements suggested by this paper for future work. They are suggesting ways to push the science further based on what they found.

Jocelyn: One big one is incorporating more detailed CR transport modeling, specifically using solutions to equation B4 with diffusion coefficients that are explicitly made functions of magnetic field strength and turbulence, as described by equations B5 and B13.

Subrahmanyan: That’s crucial because it allows them to simulate a much more realistic spatial distribution for both electrons and protons in the cloud region, which is what we need to get accurate results when we talk about cosmic ray propagation.

Vera: And they are doing parameter sensitivity analysis to systematically scan parameters to determine physical properties of CR accelerators, and they found that the SNR has accelerated protons up to more than one hundred TeV in the past <ref:2507.04407#pg3>.

Jocelyn: They also constrained diffusion coefficients like D0, which is shown to be two orders of magnitude lower than the Galactic mean in the hadronic-dominated model.

Subrahmanyan: Constraining that diffusion coefficient helps us understand if these low values are a real physical feature or just an artifact of their modeling choices, and it ties directly into structure formation models.

Vera: They also found that the hadronic model is preferred because the leptonic model is unfavored but not ruled out since there are no observational constraints on the special environment of the magnetic field.

Jocelyn: That leaves us with a bit of uncertainty about whether we should focus our next observations on testing those magnetic field environment constraints.

Subrahmanyan: And for future work, they suggest using sub-PeV gamma rays and neutrinos to test if this SNR was actually a PeVatron, which is a way to probe the very highest energy particle physics.

Conclusion: Vera: So we’ve covered the main points of "Resolving the origin of the unidentified TeV source HESS J1626-four hundred ninety as a relic of the ancient cosmic-ray factory SNR G335 point 2+zero point one" <ref:2507.04407#pg1,Resolving the origin of the unidentified TeV source HESS J1626-490 as>. We've seen how they used multi-wavelength data to link this TeV emission to an older SNR.

Jocelyn: The main implication is that the paper points toward a scenario where protons accelerated at ≳ one hundred TeV in the past, and future sub-PeV gamma ray and neutrino observations could finally test if this remnant was a PeVatron <ref:2507.04407#pg3>.

Subrahmanyan: For cosmology, it’s important because it tests how cosmic rays behave over galactic timescales and whether this remnant acts as a PeVatron.

Vera: It’s interesting stuff because the work shows that both hadronic and leptonic emissions from the cloud can explain the TeV gamma-ray spectrum.

Jocelyn: We're ready to move on to what's next in our research, and we’ve got some big questions about those magnetic fields influencing propagation.

Subrahmanyan: I just think this paper helps bridge the gap between old SNR physics and high-energy particle physics by suggesting that proton acceleration up to ≳ one hundred TeV is plausible <ref:2507.04407#pg3>.

Tomohiko Oka, Wataru Ishizaki, Masaki Mori, Hidetoshi Sano, Hiromasa Suzuki, Takaaki Tanaka

Department of Physical Sciences, Ritsumeikan University Research Organization of Science and Technology, Center for Gravitational Physics, Yukawa Institute for Theoretical Physics, Astronomical Institute Graduate School of Science Tohoku University Faculty of Engineering Gifu University Faculty of Engineering University of Miyazaki Department of Physics Konan University

astro-ph.HE

Submitted: 2025-07-06

Updated: 2025-07-06

Comments: 22 pages, 17 figures, 2 tables; Accepted for publication in ApJ

Journal ref: ApJ, 989, 137 (2025)

DOI: 10.3847/1538-4357/adec80

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 72/100

The gist: The origin of HESS J1626−490 as a relic of SNR G335.2+0.1 is investigated by analyzing multi-wavelength data to test the CR-escape scenario, which suggests that protons accelerated at the SNR

Key concepts

CR-escape scenario: This hypothesis suggests that high-energy protons accelerated within the supernova remnant (SNR) escape its magnetic confinement. These escaping protons then travel through space and interact with nearby gas clouds, producing gamma rays via pion decay ($\pi^0$-decay). This mechanism links the SNR's acceleration to the observed gamma-ray emission.
Hadronic vs. Leptonic Models: The research compared two theoretical frameworks: hadronic (proton-driven) and leptonic (electron/positron-driven). The hadronic model was favored because it avoided unrealistic assumptions about electron escape times. It implies that the SNR has accelerated protons up to 100 TeV in the past.
Diffusion Coefficient: This parameter describes how quickly cosmic rays spread or diffuse through space, influenced by magnetic field strength and turbulence. The study found that a small diffusion coefficient around SNR G335.2+0.1 is two orders of magnitude lower than the average Galactic value, suggesting limitations in current CR propagation models.
Spectral Hardening: The observed gamma-ray spectrum showed a spectral hardening at 50 GeV, connecting smoothly to the HESS spectrum. This feature was modeled as a transition between emissions from the SNR and nearby molecular clouds, rather than a simple power law.

Summary

The origin of HESS J1626−490 as a relic of SNR G335.2+0.1 is investigated by analyzing multi-wavelength data to test the CR-escape scenario, which suggests that protons accelerated at the SNR escape and illuminate nearby clouds to produce gamma rays

Observational Findings

The study scrutinized properties of SNR G335.2+0.1 and searched for non-thermal counterparts in radio, X-ray, and GeV gamma-ray bands The distance and age of the SNR are estimated to be 3.3 ± 0.6 kpc and ∼5 kyr X-ray and GeV gamma-ray data analyses revealed an extended GeV gamma-ray emission overlapping with the SNR and H.E.S.S. source regions but no X-ray emission The location of the peak of the extended GeV emission changes from near the SNR at ≲ 1 GeV to the H.E.S.S. source at > 10 GeV 10 GeV>

GeV Gamma-Ray Morphology

Fermi-LAT data analysis revealed an apparent excess in a large area containing both the SNR and H.E.S.S. source regions The GeV morphology was evaluated by fitting it with a symmetric two-dimensional (2-d) Gaussian function The location of the GeV emission gradually changes with energy, being in good agreement with the SNR in the lower energy band below 2 GeV and consistent with the H.E.S.S. source in the higher energy band above 11 GeV 11 GeV), the Gaussian parameters are consistent with the H.E.S.S. source>

Spectral Analysis and Modeling

The observed spectrum of Gaussian-W shows a spectral hardening at 50 GeV and smoothly connects to the HESS spectrum, suggesting two components, such as those from the SNR and molecular cloud To reproduce the observed broadband spectrum, a combination of emission from escaped CRs and the SNR itself was found to be consistent with observations, assuming protons were accelerated to ≳100 TeV in the past

Hadronic vs. Leptonic Models

The analysis explored both hadronic-dominated and leptonic-dominated models to explain the TeV gamma-ray spectrum The hadronic model was preferred because it avoided an unrealistic assumption about the time when electrons escape from the SNR required by certain leptonic models The hadronic model suggests that the diffusion coefficient around the SNR is two orders of magnitude lower than the Galactic mean

Conclusion

The study concluded that both leptonic and hadronic emissions from the cloud could explain the TeV gamma-ray spectrum, but due to constraints on magnetic field evolution, the hadronic model is preferred The work suggests that this SNR has accelerated protons up to ≳100 TeV in the past and future observations of sub-PeV gamma rays and neutrinos can test whether it was a PeVatron

How it works

The authors analyzed H i and 12CO (J=1−0) line data to identify the cloud associated with the SNR
They tested a CR-escape scenario where protons accelerated at SNR G335.2+0.1 escape from an acceleration site and illuminate nearby clouds, producing gamma rays via π0-decay
They modeled the observed spectrum by combining emissions from the SNR and cloud regions, assuming that CRs were accelerated in the SNR and illuminated the cloud emitting gamma rays
They found that both leptonic and hadronic emissions from the cloud can explain the TeV gamma-ray spectrum
They determined that a small diffusion coefficient around SNR G335.2+0.1 indicates a limitation in modern CR propagation models with isotropic diffusion

How it works

The study utilized multi-wavelength data, including radio continuum maps, X-ray observations from XMM-Newton, and GeV gamma-ray data from Fermi-LAT

The analysis involved fitting GeV gamma-ray maps with symmetric 2-d Gaussian functions across different energy bands to track the emission's location change with energy<ref:2507.

Improvements for AI systems

textbfFundamental CR Propagation Modeling Integration (Eqs. B4-B15): The AI system can now simulate CR transport using a full solution to Eq. B4, incorporating diffusion coefficients that are explicitly made functions of magnetic field strength and turbulence, as described by Eq. B5 and B13, allowing for a more physically realistic spatial distribution of electrons and protons in the cloud region (Eq. B14-B16).
textbfSpectral Modeling with Component Decomposition: The system can now perform a rigorous spectral decomposition using the 2-Gaussian model to explicitly separate emission from the SNR and HESS source, as shown in Fig. 16, and quantify the spectral hardening at 50 GeV by modeling it as a transition between components (SNR and cloud) rather than a simple power law.
textbfPhysically Constrained Hadronic vs. Leptonic Discrimination: The AI can evaluate the viability of emission mechanisms by comparing the resulting models against observational constraints, specifically testing the time when electrons start escaping parameter ξe,inj and finding that the hadronic model is preferred because the leptonic model is unfavored but not ruled out since there are no observational constraints on the special environment of the magnetic field.
textbfParameter Sensitivity Analysis for Accelerator Physics: The system can systematically scan parameters to determine physical properties of CR accelerators, such as finding that the SNR has accelerated protons to ≳100 TeV in the past and constraining diffusion coefficients like D0, which are shown to be two orders of magnitude lower than the Galactic mean in the hadronic-dominated model.
textbfMulti-Wavelength Kinematic Mapping: The system can correlate different observational data sets (H i and 12CO) to map physical structures, such as identifying two cavity-like structures with a similar extent to the SNR shell at specific velocity ranges (VLSR = −22 km s−1 and −45 km s−1), linking these features directly to the SNR's structure.
textbfPlanetary/Neutrino Flux Prediction: The system can predict potential future observational signatures by calculating neutrino spectra using the AAfrag package for both hadronic and leptonic dominated cases, showing that the hadronic-dominated model suggests that the neutrino flux at 1–100 TeV is compatible with the 90% confidence level (CL) sensitivity of KM3NeT.