GeV gamma-ray emission in the field of the shell-type supernova remnant Vela Jr revisited

Using 15 years of Fermi-LAT data and new eROSITA constraints, this study revises the GeV gamma-ray emission of the Vela Jr supernova remnant as primarily shell-like with minimal pulsar wind nebula contribution, and finds that a hybrid lepton-hadron model best explains the multi-wavelength spectrum, suggesting a hadronic origin for the GeV emission while the TeV emission remains predominantly leptonic.

The Gist

Imagine the universe as a giant, cosmic construction site. In this neighborhood, there's a massive, expanding bubble of debris left behind by a star that exploded thousands of years ago. This bubble is called a Supernova Remnant, and the specific one we are talking about is named Vela Jr (also known as RX J0852.0-4622).

Think of Vela Jr as a cosmic "shrapnel shell." When the star exploded, it sent a shockwave rippling outward, sweeping up gas and dust like a snowplow. This shockwave is a powerful accelerator, smashing particles together and speeding them up to nearly the speed of light.

This paper is like a team of cosmic detectives (astronomers) using a giant space camera called Fermi-LAT to take a new, ultra-sharp photo of this explosion site. They've been watching it for 15 years, which is like watching a movie in slow motion to catch every detail.

Here is what they found, explained simply:

1. The "Ghost" in the Machine

Inside this expanding shell, there is a leftover "engine" from the explosion: a rapidly spinning neutron star (a pulsar) called PSR J0855-4644. Usually, these pulsars are like powerful lighthouses, blasting out energy and creating a glowing cloud (a Pulsar Wind Nebula) that can be very bright.

The astronomers wanted to know: Is the bright light we see coming from the expanding shell, or is it just the pulsar's lighthouse shining through?

• The Analogy: Imagine a campfire (the shell) and a bright flashlight (the pulsar) sitting right next to it. From far away, it's hard to tell if the light is coming from the fire or the flashlight.
• The Discovery: By using a new, high-resolution map of the shell (based on X-ray data from a telescope called eROSITA), they realized the light matches the shape of the campfire perfectly. The flashlight (the pulsar) is there, but it's not the main source of the light in this specific view. The "shell" is doing the heavy lifting.

2. The Energy Puzzle: Two Types of Particles

The big mystery in astronomy is: What is making this light?
There are two main suspects:

• Suspect A (The Electron): Tiny, lightweight particles (electrons) that get zapped by the shockwave and glow when they hit magnetic fields. This is the "Leptonic" theory.
• Suspect B (The Proton): Heavy, nuclear particles (protons) that get smashed into gas clouds, creating a burst of energy. This is the "Hadronic" theory.

For a long time, scientists argued over which suspect was guilty. Some thought it was all electrons; others thought it was all protons.

• The Analogy: Think of it like a crime scene. The "Electron" suspect leaves a trail of footprints (synchrotron radiation) in the mud. The "Proton" suspect leaves behind a pile of smashed bricks (gamma rays from collisions).
• The Discovery: The astronomers built a computer model to test both theories. They found that both suspects are guilty, but they committed the crime at different times of the day.
  • The "Low Energy" Light (GeV): In the lower energy range, the light is a mix. About two-thirds comes from the electrons, but a significant chunk (one-third) comes from the protons smashing into gas. It's a "hybrid" crime.
  • The "High Energy" Light (TeV): When they looked at the most energetic light (the high-energy gamma rays), it was almost entirely the work of the electrons. The protons were mostly quiet here.

3. Why This Matters

This paper is important because it solves a long-standing debate.

• Before: Scientists were arguing, "Is it electrons or protons?"
• Now: They can say, "It's both, but the mix changes depending on how energetic the light is."

It's like realizing that a smoothie isn't just made of strawberries or just blueberries; it's a specific blend where the flavor changes depending on how much you stir it.

The Takeaway

The astronomers used 15 years of data to create a clearer picture of Vela Jr. They confirmed that:

1. The explosion shell is the main source of the light, not the pulsar next door.
2. The light is a "hybrid" mix of two different types of particles.
3. The heavy particles (protons) are doing a lot of the work in the lower-energy range, proving that supernova remnants are indeed the factories that create the high-energy cosmic rays that constantly rain down on Earth.

In short, they've taken a blurry, confusing photo of a cosmic explosion and turned it into a sharp, high-definition 3D movie that tells us exactly how the universe accelerates particles to incredible speeds.

Technical Summary

1. Problem Statement

The origin of $\gamma$-ray emission from the shell-type supernova remnant (SNR) RX J0852.0-4622 (Vela Jr) remains a subject of debate. While the remnant exhibits a clear shell morphology in TeV $\gamma$-rays (detected by H.E.S.S.) and X-rays, the mechanism driving the GeV emission is unclear. Two primary scenarios exist:

• Leptonic: Inverse Compton (IC) scattering of relativistic electrons off ambient photon fields.
• Hadronic: Proton-proton ($pp$) collisions producing neutral pions ($\pi^0$) that decay into $\gamma$-rays.

Previous studies have yielded conflicting results regarding the hadronic contribution, with some suggesting a mixed origin and others favoring a purely leptonic scenario. Furthermore, the potential contribution of the embedded pulsar PSR J0855-4644 and its Pulsar Wind Nebula (PWN) to the extended GeV emission requires clarification. This study aims to revisit the GeV emission using an expanded dataset (15 years of Fermi-LAT data) to refine spatial modeling, spectral constraints, and multi-wavelength (MWL) physical modeling.

2. Methodology

Data Acquisition and Processing:

• Fermi-LAT: The authors analyzed 15 years of Pass 8 data (August 2008 – May 2023) using the P8R3_SOURCE event class.
• Region of Interest (ROI): A $14^\circ \times 14^\circ$ region centered on the remnant.
• Background Modeling: Included the 4FGL-DR4 source catalog, Galactic diffuse emission (gll_iem_v07.fits), and isotropic extragalactic emission.
• X-ray Data: Utilized eROSITA eRASS1 data (1–8 keV) to construct an independent X-ray template for the entire remnant shell.
• Distance: Adopted a reference distance of $1.41 \pm 0.14$ kpc based on recent Gaia-linked measurements (Suherli et al. 2025).

Spatial Analysis:

• Generated Test Statistic (TS) maps in the 5–500 GeV range.
• Tested multiple spatial templates via binned likelihood fits:
  1. Uniform disk (standard Fermi model).
  2. Gaussian template.
  3. eROSITA X-ray template.
  4. H.E.S.S. TeV shell template.
  5. Masked H.E.S.S. template (H.E.S.S. shell with the PWN region masked).
• Model selection was performed using the Test Statistic (TS) and Akaike Information Criterion (AIC).

Spectral and Physical Modeling:

• Spectral Fitting: Tested Power Law (PL), LogParabola, PLSuperExpCutoff (PLEC), and Broken Power Law models.
• MWL Modeling: Used the Naima package to model the broadband SED (Radio to TeV) using two scenarios:
  1. Pure Leptonic: Synchrotron (Radio/X-ray) + IC (GeV/TeV).
  2. Hybrid Lepton-Hadron: Synchrotron + IC + $\pi^0$-decay (Hadronic).
• Constraints included radio data (Duncan & Green 2000), eROSITA X-ray SED, Fermi-LAT GeV points, and H.E.S.S. TeV points.

3. Key Contributions

1. Refined Spatial Modeling: The study demonstrates that the GeV morphology is best described by a masked H.E.S.S. shell template (Model 6), rather than simple geometric shapes or the unmasked H.E.S.S. template. This confirms that the embedded PWN contributes negligibly to the extended GeV flux.
2. Independent X-ray Template: Construction of a new, uniform eROSITA shell template provides a robust, independent constraint on the synchrotron emission, improving upon previous pointed X-ray observations.
3. High-Precision Spectral Fit: The analysis provides a more precise photon index for the GeV spectrum ($1.77 \pm 0.03$) and confirms a smooth connection to the TeV spectrum, reducing systematic uncertainties in joint fits.
4. Quantitative Hadronic Fraction: By fitting the broadband SED with a hybrid model, the authors quantify the relative contributions of leptonic and hadronic processes across different energy bands, moving beyond qualitative arguments.

4. Key Results

Spatial Distribution:

• The GeV emission aligns closely with the TeV shell and the eROSITA X-ray shell.
• The masked H.E.S.S. template yielded the highest $\Delta$AIC (124 units better than the standard disk model), indicating a superior fit.
• Azimuthal profile analysis shows that while the GeV and TeV emissions are inhomogeneous along the shell, the GeV data does not show the significant enhancement at the PWN location seen in TeV data, supporting the conclusion that the PWN does not dominate the GeV emission.

Spectral Properties:

• The spectrum is best fitted by a hard Power Law with a photon index of $\Gamma = 1.77 \pm 0.03$.
• The total $\gamma$-ray flux (0.1–500 GeV) is $(6.49 \pm 0.01_{stat}) \times 10^{-8}$ ph cm$^{-2}$ s$^{-1}$.
• The luminosity is estimated at $(2.47 \pm 0.49) \times 10^{33}$ erg s$^{-1}$.
• Joint Fermi-LAT/H.E.S.S. fits confirm a smooth spectral connection with a cutoff energy of $E_{cut} \approx 6.7$ TeV.

Physical Modeling (Leptonic vs. Hybrid):

• Pure Leptonic Model: Can reproduce the overall broadband shape but is statistically inferior.
• Hybrid Model: Provides a significantly better statistical description ($\Delta$AIC = 204 in favor of the hybrid model).
  • GeV Band (0.1–300 GeV): The hadronic ($\pi^0$-decay) component contributes 34%, while IC contributes 66%.
  • TeV Band (0.3–30 TeV): The emission remains predominantly leptonic, with hadronic contribution dropping to 8%.
• Magnetic Field: The hybrid model implies an average magnetic field of $B \approx 6.9 \, \mu$G, consistent with a wind-blown cavity scenario.
• Energy Budget: The total energy in cosmic ray protons ($W_p \approx 1.66 \times 10^{49}$ erg) exceeds that of electrons ($W_e \approx 3.2 \times 10^{48}$ erg) by a factor of ~5, consistent with standard SNR acceleration efficiency expectations.

5. Significance

This study resolves ambiguities regarding the origin of $\gamma$-rays from Vela Jr by leveraging a decade and a half of Fermi-LAT data and new eROSITA X-ray constraints.

• Confirmation of Mixed Origin: The work provides strong statistical evidence for a mixed-origin scenario, where the GeV emission is a combination of leptonic and hadronic processes, while the TeV emission is dominated by leptonic IC scattering.
• PWN Isolation: It clarifies that the energetic pulsar PSR J0855-4644 and its PWN do not significantly contaminate the extended GeV shell emission, allowing for a cleaner study of the SNR shock acceleration.
• Cosmic Ray Acceleration: The derived energy budget and diffusion constraints suggest that Vela Jr is an efficient accelerator of protons, with suppressed diffusion within the shell, supporting the paradigm that SNRs are the primary sources of Galactic cosmic rays.
• Methodological Advance: The successful application of the masked H.E.S.S. template and the independent eROSITA template sets a precedent for analyzing complex, overlapping high-energy sources in the Galactic plane.