🔭 astrophysics

Investigating the emission signatures of pulsar halo candidate HESS J1813-126

This paper investigates HESS J1813-126 as a candidate pulsar halo powered by PSR J1813-1246 by modeling the pulsar's emission with a synchro-curvature framework and testing various particle transport scenarios against multi-wavelength data to derive predictive surface brightness profiles and aperture-dependent signatures for future observational discrimination.

Original authors: Agnibha De Sarkar

Published 2026-02-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Agnibha De Sarkar

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, bustling city. In this city, pulsars are like incredibly fast-spinning lighthouses or power plants. They spin hundreds of times a second, shooting out beams of energy and particles (like tiny, super-fast electrons and positrons) into the surrounding space.

Usually, when these particles get shot out, they get trapped in a magnetic "bubble" right next to the pulsar, creating a glowing cloud called a Pulsar Wind Nebula (PWN). This is like a bright, tight fog right around the lighthouse.

However, sometimes, these particles escape the bubble and drift far away into the interstellar "streets" of the galaxy. As they drift, they bump into light from stars and the cosmic background, getting supercharged and glowing in high-energy gamma rays. This creates a massive, faint, diffuse glow that can stretch for dozens of light-years. Astronomers call this a Pulsar Halo. It's like the faint, wide halo of light you might see around a streetlamp on a foggy night, stretching much further than the bright bulb itself.

The Mystery of HESS J1813-126

The paper investigates a specific cosmic object called HESS J1813-126. For a long time, astronomers have been puzzled by this object. It's a huge, bright spot of gamma rays in the sky. Right in the middle of it sits a pulsar named PSR J1813-1246.

The big question was: Is this giant glow a "Pulsar Halo" created by the pulsar, or is it something else entirely, like a leftover from a supernova explosion?

How the Authors Solved the Puzzle

The researchers acted like cosmic detectives, using three main tools to figure out what was going on:

1. Checking the "Fingerprint" (The Pulsar's Light)
First, they looked closely at the pulsar itself using the Fermi-LAT telescope (which sees gamma rays) and X-ray telescopes. They used a special mathematical model called Synchro-Curvature radiation to explain how the pulsar shines.

The Analogy: Think of the pulsar as a rollercoaster. The electrons are the cars. As they spiral down the track (magnetic field lines), they emit light. The authors showed that the "fingerprint" of the light coming from the pulsar matches perfectly with what we expect from a powerful, spinning engine. This confirmed the pulsar is the "engine" behind the scene.

2. Testing the "Drift" (The Transport Models)
Next, they tried to simulate how the particles escape the pulsar and spread out to create the giant halo. They tested three different theories on how these particles move through the galaxy:

The "Two-Zone" Model (2ZISD): Imagine the particles are in a crowded room (near the pulsar) where they can't move fast, but once they get to the hallway (further out), they can run freely.
The "Ballistic" Model (B2D): Imagine the particles are shot out like bullets, flying straight for a while before slowing down.
The "Anisotropic" Model (AD): Imagine the particles are like fish swimming in a river; they move easily with the current (magnetic field) but struggle to swim across it.

The Verdict:
They compared their simulations with real data from three giant telescopes: H.E.S.S., HAWC, and LHAASO.

The "Two-Zone" and "Anisotropic" models worked beautifully. They predicted exactly how bright the halo should be and how its shape changes depending on how far out you look.
The "Ballistic" model failed. It predicted the glow would be too big and too bright in the outer edges, which didn't match the telescopes' observations.

3. The "Budget" Check (Energy Accounting)
Finally, they did a math check to see if the pulsar had enough "fuel" to power the giant halo.

The Analogy: It's like checking if a car engine has enough horsepower to pull a massive trailer. They calculated the number of particles the pulsar creates versus the number needed to light up the halo.
The Result: The numbers matched! The pulsar is powerful enough to create the halo without needing any "magic" extra energy sources.

The Big Picture Conclusion

The paper concludes that HESS J1813-126 is almost certainly a Pulsar Halo.

Here is why:

No "Fog": There is no tight, bright cloud (PWN) right next to the pulsar, which suggests the particles have escaped.
The Shape: The glow is huge and spread out, exactly like a halo, not a shell from an explosion.
The Match: The math models for how particles drift through space fit the real telescope data perfectly.

Why does this matter?
This discovery helps us understand how pulsars interact with their surroundings. It also helps solve a mystery about "positrons" (anti-electrons) hitting Earth. These halos might be the source of the extra positrons we detect, acting as local factories for antimatter.

In short, the authors have confirmed that this distant cosmic object is a beautiful, massive "halo" of light, powered by a spinning neutron star that has been shooting particles into the galaxy for thousands of years. Future telescopes will be able to take even sharper pictures of this halo to see exactly how the particles are moving, like watching the wind blow through a field of grass.

1. Problem Statement

Extended gamma-ray sources surrounding middle-aged pulsars, known as pulsar halos, are a distinct phenomenon from compact Pulsar Wind Nebulae (PWNe). These halos are hypothesized to be formed by relativistic electron-positron pairs ( $e^\pm$ ) diffusing into the interstellar medium (ISM) and producing gamma-rays via Inverse-Compton (IC) scattering. However, the nature of specific candidates remains uncertain due to the lack of clear boundaries and the difficulty in distinguishing them from shell-type supernova remnants (SNRs) or confined PWNe.

HESS J1813-126, associated with the energetic, radio-quiet pulsar PSR J1813-1246, is a prime candidate. While it exhibits extended TeV morphology and a hard spectrum, its origin requires rigorous testing against alternative explanations (e.g., residual magnetospheric emission or a GeV PWN). The paper aims to:

Model the high-energy emission of the central pulsar (PSR J1813-1246).
Determine if the extended gamma-ray SED of HESS J1813-126 can be consistently explained as a pulsar halo powered by PSR J1813-1246.
Test different particle transport models to predict observable signatures (Surface Brightness Profiles and aperture-dependent flux) for future telescopes.

2. Methodology

The study employs a multi-faceted approach combining observational data analysis with theoretical modeling:

Observational Data Analysis:
- Fermi-LAT: Performed standard and phased analyses of GeV data. The authors searched for off-pulse emission to rule out a GeV PWN. They utilized the Bayesian blocks algorithm to define off-pulse intervals and tested spatial models (RadialDisk vs. RadialGaussian) to determine source extension.
- Multi-wavelength Data: Incorporated X-ray data (XMM-Newton, INTEGRAL) and TeV data from H.E.S.S., HAWC, and LHAASO (including >100 TeV detections).
Pulsar Magnetosphere Modeling (Synchro-Curvature):
- Used the Synchro-Curvature (SC) radiation model to explain the combined X-ray and GeV spectra of PSR J1813-1246. This model accounts for particles spiraling and sliding along curved magnetic field lines in the outer magnetosphere, interpolating between synchrotron and curvature radiation regimes.
Pulsar Halo Transport Modeling:
- Simulated the transport and energy losses of $e^\pm$ injected by the pulsar using time-dependent diffusion-loss equations.
- Tested three distinct transport scenarios using the PHECT code:
  1. Two-Zone Isotropic Suppressed Diffusion (2ZISD): Assumes a suppressed diffusion coefficient in an inner zone and standard ISM diffusion in an outer zone.
  2. Ballistic-to-Diffusion (B2D): Assumes quasi-ballistic propagation initially, transitioning to diffusion.
  3. Anisotropic Diffusion (AD): Accounts for diffusion along and across magnetic field lines, dependent on the viewing angle relative to the field.
- Calculated the resulting IC emission (scattering off CMB, dust, and starlight) and compared it with observational SEDs.
- Generated predictions for Surface Brightness Profiles (SBPs) and aperture-dependent emission by convolving model predictions with Gaussian Point Spread Functions (PSFs) representative of current and future instruments (CTA, SWGO).
Energetic Consistency Check:
- Calculated the pair multiplicities ( $\kappa$ ) for both the magnetospheric SC emission and the halo injection to verify if the pulsar's spin-down power can energetically support both components without exotic pair-creation mechanisms.

3. Key Results

A. Pulsar Properties and GeV Emission

Fermi-LAT Analysis: The source 4FGL J1813.4-1246 is confirmed as the phase-averaged GeV counterpart of PSR J1813-1246. It is point-like (consistent with a pulsar, not an extended nebula) and has a significantly curved spectrum best fit by a PLSuperExpCutoff4 function.
Off-Pulse Search: A faint source (PS J1813.3-1246) was detected during the off-pulse phase. However, its spectrum is a soft power law, and it is point-like. The authors conclude this is likely residual magnetospheric contamination (imperfect gating) rather than a GeV PWN, as a PWN would be spatially extended and have a harder spectrum.
SC Model Fit: The SC framework successfully reproduces the X-ray and GeV SEDs. The model indicates X-rays are produced via synchrotron radiation at low altitudes, while GeV gamma-rays are produced via curvature radiation further out in the magnetosphere.

B. Pulsar Halo Modeling (HESS J1813-126)

The study evaluated the three transport models against multi-instrument data (H.E.S.S., HAWC, LHAASO):

2ZISD Model (Two-Zone Isotropic Suppressed Diffusion):
- Performance: Provides the most consistent description of the data. It reproduces the H.E.S.S. SED within the standard aperture (0.34°) and, crucially, does not overshoot the HAWC and LHAASO fluxes when larger apertures are considered.
- Prediction: Predicts that increasing the aperture primarily enhances the low-energy flux (8–40 TeV) because lower-energy electrons diffuse further before cooling. High-energy flux saturates quickly. This results in an energy-dependent morphology where low-energy emission is extended, and high-energy emission is compact.
B2D Model (Ballistic-to-Diffusion):
- Performance: Can match the H.E.S.S. data within a small aperture but fails at larger apertures. The quasi-ballistic propagation leads to flatter SBPs, causing the integrated flux to overshoot HAWC and LHAASO measurements significantly.
- Issue: Requires a high conversion efficiency ( $\eta \gtrsim 10\%$ ), which is higher than typical expectations for pulsar halos.
AD Model (Anisotropic Diffusion):
- Performance: The results depend heavily on the viewing angle ( $\zeta$ $ζ$ ) relative to the magnetic field.
  - Aligned ( $\zeta = 0^\circ$ ): Produces a nearly isotropic, circular morphology that fits the data well, similar to the 2ZISD case.
  - Misaligned ( $\zeta = 10^\circ$ ): Produces elongated morphologies and flatter SBPs that overshoot the observed flux at wider apertures.
- Implication: The data favors a field-aligned configuration or a spherical morphology, ruling out strong asymmetry.

C. Energetic Consistency

The calculated high-energy pair multiplicities are $\kappa_{mag} \approx 5.9 \times 10^4$ (magnetosphere) and $\kappa_{halo} \approx 2.2 \times 10^2$ (halo).
These values are consistent with standard pulsar physics and indicate that the pulsar's spin-down power is sufficient to supply the high-energy pairs for both the magnetospheric emission and the halo, without requiring additional pair-creation channels.

4. Key Contributions

Confirmation of Halo Candidate: Strongly supports the interpretation of HESS J1813-126 as a pulsar halo powered by PSR J1813-1246, based on the absence of a GeV/X-ray PWN, the extended TeV morphology, and the successful reproduction of the SED by halo models.
Diagnostic Signatures: Provides specific, testable predictions for Surface Brightness Profiles (SBPs) and aperture-dependent flux evolution.
- The 2ZISD and aligned AD models predict that future high-resolution Imaging Atmospheric Cherenkov Telescopes (IACTs) like CTA will observe a significant increase in flux at lower energies as the aperture grows, while high-energy flux remains saturated.
- This contrasts with the B2D model, which predicts a uniform increase across energies that would be ruled out by current data.
Unified Modeling: Successfully links the central engine (via SC radiation) and the extended halo (via diffusion models) within a single energetic framework, demonstrating self-consistency in the pair budget.

5. Significance

This work advances the understanding of pulsar halos by:

Refining Transport Physics: It demonstrates that the 2ZISD model (suppressed diffusion in the inner zone) is currently the most viable explanation for HESS J1813-126, offering a robust alternative to the B2D model which is disfavored by multi-instrument data.
Guiding Future Observations: The derived SBP predictions serve as a roadmap for upcoming instruments (CTA, SWGO). Observing the specific energy-dependent morphological changes predicted by the 2ZISD model will be the "smoking gun" to confirm the diffusion physics governing these halos.
Resolving the PeVatron Mystery: By confirming HESS J1813-126 as a pulsar halo with emission extending to >100 TeV, the paper reinforces the idea that pulsar halos are significant contributors to the Galactic "PeVatron" population, potentially explaining the cosmic-ray positron excess and the origin of ultra-high-energy gamma rays.

In conclusion, the paper establishes HESS J1813-126 as a textbook example of a pulsar halo, providing a coherent physical picture that connects the pulsar's magnetospheric emission to its extended interstellar halo, while offering clear observational tests to distinguish between competing diffusion theories.