🔭 astrophysics

Constraining the slow-diffusion zone size and electron injection spectral index for the Geminga pulsar halo

Using updated HAWC observations of the Geminga pulsar halo within a two-zone diffusion model, this study constrains the slow-diffusion zone size to 30–70 pc and the electron injection spectral index to a maximum of 2.17, finding that a specific combination of these parameters successfully reproduces the AMS-02 positron spectrum at Earth.

Original authors: Kun Fang

Published 2026-04-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Kun Fang

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, there are "pulsars"—dead stars that spin incredibly fast and shoot out beams of energy like cosmic lighthouses. Around these lighthouses, there are clouds of high-speed particles (electrons and positrons) that have escaped the star's immediate neighborhood. These clouds are called Pulsar Halos.

For a long time, scientists knew these halos existed, but they were like trying to understand a foggy morning by looking at a blurry photo. They knew the fog was there, but they couldn't tell exactly how thick it was or how fast the particles were moving through it.

Recently, a giant telescope called HAWC (High-Altitude Water Cherenkov) took a much sharper, higher-resolution picture of the halo around a famous pulsar named Geminga. This new "photo" allowed the author of this paper, Kun Fang, to solve a cosmic mystery: How big is the "slow zone" around Geminga, and how fast are the particles being shot out?

Here is the story of the findings, explained simply:

1. The "Traffic Jam" Analogy (The Slow-Diffusion Zone)

Imagine the space around a pulsar is a highway. Usually, particles (cars) zoom along at high speeds. But right next to Geminga, something strange happens. The space is like a heavy traffic jam. The particles get stuck and move very slowly.

Scientists call this the "Slow-Diffusion Zone." Once the particles escape this traffic jam, they hit the open highway (the rest of the galaxy) and zoom away quickly.

The big question was: How big is this traffic jam?

Too small? If the jam is tiny, the particles escape too fast, and the "fog" (the halo) looks too small and dim.
Too big? If the jam is huge, the particles stay trapped too long, and the energy spectrum (the color of the light) doesn't match what we see.

2. The Detective Work: Measuring the "Fog"

Kun Fang used the new HAWC data to act like a detective, testing different sizes for this traffic jam.

The Lower Limit (The Minimum Size):
The author looked at the shape of the halo (how the light fades as you move away from the center). He found that if the traffic jam were smaller than 30 light-years (about 30 parsecs), the model didn't fit the data at all. The "fog" would disappear too quickly.
- Analogy: It's like trying to explain a giant puddle of water by saying the leak was only a tiny drop. It just doesn't add up. The data says the slow zone must be at least 30 light-years wide.
The Upper Limit (The Maximum Size):
Next, the author looked at the energy of the light (the gamma rays). He knew from physics simulations that the particles being shot out of the pulsar can't be too "hard" (too energetic) or the math breaks down.
He found that if the traffic jam were larger than 70 light-years, the particles would have to be shot out with impossible energy levels to match the observations.
- Analogy: If the traffic jam were a whole city block wide, the cars would have to be driving at the speed of light just to get out, which is physically impossible.

The Verdict: The slow-diffusion zone around Geminga is likely between 30 and 70 light-years wide.

3. The "Firehose" (Electron Injection)

The paper also looked at how the pulsar shoots out these particles. Think of the pulsar as a firehose spraying water.

The "spectral index" is a fancy way of describing the pressure of the hose. Is it a gentle mist (soft spectrum) or a high-pressure jet (hard spectrum)?
The study found that the pulsar is shooting out particles with a specific pressure. If the slow zone is at its smallest (30 light-years), the "pressure" (spectral index) can be at most 2.17. If the zone is bigger, the pressure must be lower.

4. The Big Mystery: The "Positron Excess"

Why do we care? Because there is a mystery in our own solar system.

The Mystery: Detectors on Earth (like AMS-02 on the International Space Station) are finding way more positrons (anti-electrons) than they should. It's like finding a pile of gold coins in your backyard when you only expected a few copper pennies.
The Suspect: Scientists think Geminga might be the culprit, shooting these positrons all the way to Earth.
The Test: Kun Fang ran the numbers. He asked: "If Geminga is the one shooting these particles, and the slow zone is 30 light-years wide, does the spray reach Earth with the right amount of gold coins?"
The Result: Yes! When using the "sweet spot" parameters (Slow zone = 30 light-years, Pressure = 2.17), the predicted positrons hitting Earth match the AMS-02 data perfectly in the 50–500 GeV range.

Summary

This paper is like a cosmic detective story. By using a sharper telescope (HAWC), the author narrowed down the size of the "traffic jam" around the Geminga pulsar to be between 30 and 70 light-years.

This discovery is crucial because it confirms that Geminga is a very strong candidate for the source of the mysterious positron excess we see on Earth. It tells us that the universe has a specific "slow zone" around these stars, and once the particles escape that zone, they travel across the galaxy to reach us, solving a puzzle that has baffled scientists for years.

1. Problem Statement

The study of gamma-ray pulsar halos, particularly the Geminga halo, is crucial for understanding cosmic-ray propagation and the origin of the cosmic-ray positron excess. While previous studies successfully constrained the diffusion coefficient ( $D$ ) around pulsars, they often relied on simplified models (e.g., one-zone diffusion) or fixed assumptions regarding the size of the "slow-diffusion zone" ( $r^*$ ) and the electron injection spectral index ( $p$ ).

The core problems addressed in this paper are:

Parameter Degeneracy: The size of the slow-diffusion zone ( $r^*$ ) and the spectral index of the injected electrons ( $p$ ) are degenerate parameters when fitting observational data. A change in one can mimic the effect of the other.
Lack of Constraints: Previous constraints on $r^*$ were often based on general estimates rather than direct fits to high-precision morphology data.
Positron Excess Interpretation: It remains unclear whether Geminga can account for the positron excess measured by AMS-02 without violating constraints from the gamma-ray spectrum observed by HAWC.

2. Methodology

The author employs a two-zone diffusion model to simulate electron propagation from the Geminga pulsar into the Interstellar Medium (ISM).

Physical Model:
- Diffusion: The model divides space into two regions: a slow-diffusion zone of radius $r^*$ around the pulsar and the standard ISM outside. The diffusion coefficient $D(E, r)$ is defined as:
  $D(E, r) = \begin{cases} D_{100}(E/100 \text{ TeV})^\delta & |r - r_p| \le r^* \\ D_{100, \text{ism}}(E/100 \text{ TeV})^{\delta_{\text{ism}}} & |r - r_p| > r^* \end{cases}$
  where $D_{100, \text{ism}} = 3.4 \times 10^{30} \text{ cm}^2 \text{ s}^{-1}$ and $\delta_{\text{ism}} = 0.38$ .
- Source Term: The electron injection spectrum $q(E_e)$ is modeled as a power law with a super-exponential cutoff: $q(E_e) \propto E_e^{-p} \exp[-(E_e/E_c)^2]$ .
- Propagation: The electron transport is solved using the diffusion-loss equation, accounting for synchrotron and Inverse Compton (IC) energy losses.
Observational Data:
- Morphology: The 1D surface brightness profile (SBP) of the Geminga halo measured by HAWC (High-Altitude Water Cherenkov experiment).
- Spectrum: The gamma-ray energy spectrum of the Geminga halo measured by HAWC (1–100 TeV).
Fitting Strategy:
- Step 1 (Morphology): Fit the SBP to constrain the lower limit of $r^*$ . The Monogem halo contribution is subtracted using a one-zone model.
- Step 2 (Spectrum): Fit the integrated gamma-ray spectrum (within $10^\circ$ ) to constrain the injection spectral index $p$ and the upper limit of $r^*$ .
- Statistical Tools: Goodness-of-fit tests ( $\chi^2$ ) and Likelihood Ratio tests (Wilks' theorem) are used to determine confidence intervals ( $2\sigma$ and $3\sigma$ ).

3. Key Contributions

Refined Two-Zone Modeling: The paper moves beyond simple diffusion coefficient measurements to simultaneously constrain the spatial extent of the slow-diffusion zone ( $r^*$ ) and the injection physics ( $p$ ).
Direct Data-Driven Constraints: Unlike previous works that estimated $r^*$ based on theoretical source sizes, this study derives $r^*$ directly from the goodness-of-fit to the HAWC morphology and spectrum data.
Resolution of Positron Excess Tension: The study demonstrates that specific parameter combinations allowed by HAWC data can simultaneously explain the hard gamma-ray spectrum of Geminga and the positron spectrum measured by AMS-02.

4. Key Results

A. Constraints on Slow-Diffusion Zone Size ( $r^*$ )

Lower Limit ( $r^* \gtrsim 30$ pc): Models with $r^* < 30$ pc are excluded at the 95% confidence level. Small $r^*$ values cause the electron density to drop too sharply at the boundary, resulting in a surface brightness profile that declines too quickly compared to HAWC data.
Upper Limit ( $r^* \lesssim 70$ pc): The upper limit is derived by fitting the gamma-ray spectrum. As $r^*$ increases, the required injection spectral index $p$ to fit the hard HAWC spectrum decreases. Since theoretical simulations and X-ray observations suggest $p \gtrsim 1$ , models with $r^* \ge 70$ pc (which require $p < 1$ ) are disfavored.
Final Constraint: The slow-diffusion zone size around Geminga is constrained to $30 \text{ pc} \le r^* \le 70 \text{ pc}$ .

B. Constraints on Injection Spectral Index ( $p$ )

Inverse Correlation: There is a strong inverse correlation between $r^*$ and the maximum allowed $p$ . A smaller $r^*$ allows for a harder injection spectrum.
Maximum $p$ : When $r^*$ is set to its lower limit (30 pc) and the diffusion index $\delta = 0$ (energy-independent diffusion), the maximum allowed spectral index is $p \le 2.17$ at the $3\sigma$ confidence level.
Dependence on $\delta$ : If the diffusion is energy-dependent (e.g., Kolmogorov $\delta = 1/3$ ), the upper limit on $p$ becomes stricter (e.g., $p \le 1.82$ ).

C. Implications for Positron Excess

AMS-02 Consistency: Using the most optimistic parameters ( $r^* = 30$ pc, $\delta = 0$ , $p = 2.17$ ), the predicted positron flux from Geminga at Earth matches the AMS-02 measurements in the 50–500 GeV range in both spectral shape and absolute flux.
Cutoff Discrepancy: The model predicts a spectral cutoff at a higher energy ( $\sim 800$ GeV) than observed by AMS-02. The author suggests this could be due to higher magnetic fields in the positron's path or that an older pulsar (e.g., PSR B1055-52) contributes significantly to the highest energies.

5. Significance

Validation of Pulsar Origin: The results support the hypothesis that pulsars are the dominant source of the cosmic-ray positron excess, provided the diffusion environment is a slow-diffusion zone of 30–70 pc.
Physical Insight: The derived size of the slow-diffusion zone (30–70 pc) is consistent with the scale of the pulsar's surrounding supernova remnant or the region where escaped electrons amplify magnetic turbulence.
Future Prospects: The paper highlights that the Large High Altitude Air Shower Observatory (LHAASO) will be critical in breaking the degeneracy between $r^*$ , $p$ , and $\delta$ by providing energy-dependent morphology measurements.
Tension with H.E.S.S.: The author notes a potential conflict with recent H.E.S.S. data, which shows a softer spectrum in the 1–10 TeV range within 1 degree, suggesting the need for further verification via LHAASO.

In summary, this work provides the most stringent, data-driven constraints to date on the physical parameters of the Geminga pulsar halo, bridging the gap between gamma-ray morphology, spectral hardness, and cosmic-ray positron observations.