H.E.S.S. detection of the PSR J0855-4644 nebula

1. Problem Statement

The region surrounding the supernova remnant (SNR) RX J0852.0−4622 (commonly known as Vela Junior) is a complex source of TeV $\gamma$-ray emission. Located near the pulsar PSR J0855−4644, previous observations by the High Energy Stereoscopic System (H.E.S.S.) identified the SNR shell but could not definitively separate potential emission from the pulsar's wind nebula (PWN) due to:

• Spatial Overlap: The PWN is projected onto the bright shell of the Vela Junior SNR.
• Analysis Limitations: Previous 2D morphological and spectral analyses lacked the sensitivity to distinguish the distinct spectral and morphological signatures of the PWN from the SNR shell.
• Unresolved Hypothesis: While PSR J0855−4644 is a known X-ray and radio PWN, its existence as a TeV source remained unconfirmed, with previous data suggesting only a possible spectral hardening that lacked statistical significance.

2. Methodology

The authors employed advanced analysis techniques on a comprehensive dataset to resolve the complex source structure:

• Data Collection: Utilized H.E.S.S. observations from January 2004 to February 2021 (approx. 130 hours on-source). This included data from the initial 4-telescope array and 55 hours from the upgraded camera system (post-2016).
• 3D Binned Likelihood Analysis: Unlike previous 2D analyses, the team performed a 3D spectro-morphological analysis using the Gammapy framework. This approach simultaneously models the spatial distribution and energy spectrum of the source and background.
  • Region: A $6^\circ \times 6^\circ$ analysis region centered on Vela Junior.
  • Energy Binning: 13 logarithmic bins between 0.7 and 100 TeV.
  • Background Modeling: Used a 3D Field-of-View (FoV) background model adapted to specific observation conditions.
• Iterative Modeling: The team started with a background-only hypothesis and iteratively added parametric components (spatial templates and spectral functions) until residuals were statistically insignificant ($<4\sigma$).
  • Components: They modeled the SNR shell (using a template and a uniform shell model) and added Gaussian components to capture substructures.
• Multi-wavelength Joint Fitting: A one-zone leptonic model was applied to jointly fit X-ray data (from XMM-Newton, 2–8 keV) and $\gamma$-ray data (H.E.S.S.) using the Naima and Sherpa packages to constrain the electron population and magnetic field.

3. Key Contributions

• First TeV Detection of PSR J0855−4644 PWN: This paper reports the first unambiguous detection of the pulsar wind nebula associated with PSR J0855−4644 at TeV energies.
• Source Separation: Successfully disentangled the PWN emission from the bright Vela Junior SNR shell, a feat not achieved in previous H.E.S.S. publications (e.g., H.E.S.S. Collaboration et al. 2018b).
• Advanced 3D Analysis: Demonstrated the efficacy of 3D likelihood analysis in resolving complex, overlapping Galactic sources where standard 2D methods fail.
• Broadband Modeling: Provided a consistent physical model linking X-ray synchrotron emission and TeV inverse-Compton emission for this specific PWN.

4. Key Results

Morphological and Spectral Findings

• New Component (Component E): A distinct, extended $\gamma$-ray source was identified coincident with PSR J0855−4644.
  • Significance: Detected at 12.2$\sigma$.
  • Extension: Gaussian width $\sigma_E = 0.14^\circ \pm 0.02^\circ$ (significantly larger than the X-ray PWN size of $\sim 45''$).
  • Position: $l = 266.91^\circ \pm 0.03^\circ$, $b = -0.94^\circ \pm 0.02^\circ$, matching the pulsar location.
• Spectral Properties:
  • Component E exhibits a hard power-law spectrum with an index $\Gamma_E = 1.81 \pm 0.07$.
  • This spectrum is statistically distinct from the surrounding SNR shell components (which have softer indices and exponential cutoffs).
  • No significant spectral cutoff was found for Component E up to $\sim 10$ TeV.
• Energy-Dependent Morphology: Flux maps show that while the SNR shell dominates at lower energies ($<4.8$ TeV), the emission at energies $>10.3$ TeV is dominated by the hard-spectrum component coincident with the pulsar.

Physical Parameters (Joint Fit)

Using a one-zone leptonic model (Synchrotron + Inverse-Compton):

• Electron Spectral Index: $\alpha = 1.88 \pm 0.01$ (consistent with particle acceleration at a termination shock).
• Magnetic Field: A lower limit of $B \geq 1.6 \, \mu\text{G}$ was derived. This low value is consistent with the expansion of the nebula and radiative losses as particles move away from the termination shock.
• Luminosity & Efficiency:
  • TeV luminosity ($1-10$ TeV): $L_{1-10\text{TeV}} \approx 3.2 \times 10^{32} \, \text{erg s}^{-1}$ (assuming $d=900$ pc).
  • Conversion efficiency: $\sim 0.03\%$ of the pulsar's spin-down power ($\dot{E} = 1.1 \times 10^{36} \, \text{erg s}^{-1}$), which is typical for evolved PWNe.

5. Significance

• Confirmation of PWN Population: This detection adds PSR J0855−4644 to the growing population of known TeV PWNe, confirming that high-spin-down pulsars in the Galactic plane are efficient TeV accelerators.
• Methodological Benchmark: The study validates the use of 3D analysis techniques for resolving complex, overlapping sources in the Galactic plane, offering a template for future analyses of similar regions (e.g., the Galactic Center or other SNR-PWN complexes).
• Physical Insight: The discrepancy between the X-ray and TeV sizes of the nebula provides direct evidence of electron cooling and transport effects, while the joint spectral fit constrains the magnetic field environment of the nebula, supporting standard PWN evolution models.
• Fermi-LAT Context: The hard spectrum explains why this PWN was not detected by Fermi-LAT in the GeV range; the flux drops steeply at lower energies, making it undetectable against the background in that band.

In conclusion, the paper establishes PSR J0855−4644 as a confirmed TeV pulsar wind nebula, resolving a long-standing ambiguity in the Vela Junior region through advanced 3D analysis and multi-wavelength modeling.