A Method for Testing Diffusive Shock Acceleration and Diffusion Propagation of 1-100 TeV Cosmic Electron with Multi-wavelength Observation of Geminga Halo and Pulsar Wind Nebula

This paper develops a method to test diffusive shock acceleration and diffusion propagation theories using multi-wavelength observations of the Geminga pulsar wind nebula and halo, finding current data consistent with the models but limited by wide energy bins, while anticipating that future high-precision observations will enable rigorous validation of the energy-dependent diffusion coefficient up to sub-PeV energies.

The Gist

Imagine the universe as a giant, chaotic highway where tiny particles called cosmic rays (mostly electrons in this story) zoom around at nearly the speed of light. For decades, scientists have had a theory about how these particles get so fast: Diffusive Shock Acceleration (DSA).

Think of DSA like a cosmic ping-pong game. Imagine a shockwave (like a sonic boom from a supersonic jet) moving through space. Particles bounce back and forth across this shockwave, hitting it from both sides. Every time they cross, they get a little "kick" of energy, speeding up faster and faster until they become ultra-high-energy cosmic rays.

This theory works great for low-energy particles (like those near our Sun), but for the super-fast, high-energy ones (trillions of times more energetic), we haven't been able to prove it works directly. It's like having a perfect recipe for a cake but never having tasted the final product.

The Geminga "Cosmic Lighthouse"

Enter Geminga, a dead star (pulsar) about 800 light-years away. It's like a cosmic lighthouse spinning in the dark. As it spins, it shoots out a wind of particles. When this wind hits the slower gas of space, it creates a massive "termination shock"—the perfect ping-pong table for our cosmic particles.

Around Geminga, there is a giant, invisible bubble of high-energy particles called a halo. This halo is so big it covers a huge chunk of the sky. Because the particles are so energetic, they glow in gamma rays (a type of light we can't see with our eyes, but telescopes can).

The Detective Work: Two Clues

The authors of this paper acted like detectives trying to solve a mystery: "Does the DSA theory actually explain how Geminga's halo works?"

They used two different clues to check the theory:

1. The Energy Spectrum (The "Speedometer"): They looked at how much energy the particles have. If the DSA theory is right, the number of particles at different speeds should follow a specific pattern, like a smooth slide.
2. The Shape (The "Footprint"): They looked at how spread out the halo is. If the particles diffuse (spread out) too slowly, the halo would be a tight, small ball. If they diffuse too fast, it would be a huge, fuzzy cloud.

The Big Discovery: A "Traffic Jam" in Space

Here is the tricky part. The scientists found that the particles near Geminga behave strangely.

• The Analogy: Imagine a highway where cars (particles) usually drive at a steady speed. But near Geminga, there's a massive traffic jam caused by a "slow-diffusion zone." The particles are stuck in a local neighborhood, unable to escape quickly.
• The Twist: The data showed that these particles are trapped in this traffic jam up to a certain speed (around 100 TeV). But once they get faster than that, the traffic jam suddenly clears up, and they can zoom away freely again.

The paper shows that the DSA theory (the ping-pong game) combined with this diffusion model (the traffic jam clearing up) perfectly matches what the telescopes (HAWC and Fermi-LAT) are seeing. The "speedometer" reading and the "footprint" shape of the halo tell the same story.

Why This Matters

For a long time, we could only test these theories with low-energy particles. This paper is like finally getting a test drive for the high-speed version of the theory.

• The Result: The theory holds up! The way particles accelerate and spread out in the Geminga halo matches the predictions of the standard model.
• The Catch: The data we have right now is a bit like looking at a blurry photo. The "traffic jam" (the diffusion coefficient) changes very rapidly above 100 TeV, and our current telescopes can't see the details clearly enough to be 100% sure about the exact shape of that change.

The Future

The authors are excited because future telescopes (like LHAASO) will act like high-definition cameras. They will be able to see the "traffic jam" in sharp detail. If they confirm that the diffusion coefficient (how fast particles spread) spikes up exactly as predicted above 100 TeV, it will be the "smoking gun" proof that our understanding of how the universe accelerates particles is correct.

In a nutshell: The scientists used the Geminga pulsar as a giant laboratory to test how cosmic particles get supercharged. They found that the standard "ping-pong" acceleration theory works, but only if you assume the particles get stuck in a local traffic jam that suddenly opens up at very high speeds. It's a major step forward in understanding the high-energy universe.

Technical Summary

Here is a detailed technical summary of the paper "A Method for Testing Diffusive Shock Acceleration and Diffusion Propagation of 1-100 TeV Cosmic Electron with Multi-wavelength Observation of Geminga Halo and Pulsar Wind Nebula."

1. Problem Statement

Diffusive Shock Acceleration (DSA) and diffusion propagation are foundational components of the standard cosmic ray model. While DSA is robustly validated in the MeV energy range (e.g., solar processes), its direct validation at sub-PeV energies (1 TeV to several hundred TeV) remains elusive.

• The Challenge: Existing models often rely on extrapolation. There is a need to test whether DSA and specific diffusion mechanisms (specifically the "SNR-induced" slow-diffusion model) hold true for high-energy electrons in the 1–100 TeV range.
• The Target: The Geminga pulsar and its surrounding wind nebula (PWN) and extended gamma-ray halo. Geminga is an ideal laboratory because it exhibits a PWN with a termination shock (TS) and a large-scale halo where high-energy electrons are confined, offering multi-wavelength data (X-ray from the PWN, GeV/TeV from the halo).

2. Methodology

The authors developed a self-consistent validation method linking the acceleration mechanism (DSA) with particle propagation (diffusion) using multi-wavelength data from Fermi-LAT (GeV) and HAWC (TeV).

A. Theoretical Framework

1. Electron Acceleration (DSA):

   • Modeled electron acceleration at the Termination Shock (TS) of the Geminga PWN.
   • Solved the 1D diffusion-convection equation for a bounded upstream region (finite size $a$), moving beyond the classic infinite upstream assumption.
   • Derived electron spectra under two boundary conditions:
     • Zero-particle condition: Electrons captured by the pulsar at the far upstream.
     • Zero-streaming condition: Electrons reflected/cooled in the far upstream (analogous to solar wind TS).
   • Established the relationship between the accelerated electron spectrum $f_0(p)$ and the upstream diffusion coefficient $D_P$.

2. Electron Propagation:

   • Modeled electron diffusion from the PWN into the surrounding halo.
   • Assumed a broken power-law diffusion coefficient $D(E)$ to account for potential changes in turbulence regimes:
     $$D(E) = \begin{cases} D_c (E/E_c)^{\delta_1} & E < E_c \\ D_c (E/E_c)^{\delta} & E \ge E_c \end{cases}$$
   • Calculated the electron density $N(E, r, t)$ in the halo considering injection, diffusion, and energy losses (Inverse Compton scattering).

3. Radiation Modeling:

   • X-ray: Synchrotron radiation from the PWN (constrained by Chandra data).
   • Gamma-ray: Inverse Compton (IC) scattering off CMB and interstellar dust photons (constrained by Fermi-LAT and HAWC data).

B. Consistency Check & Fitting

• Linking PWN and Halo: The authors derived a relationship between the diffusion coefficients in the PWN ($D_P$) and the halo ($D_H$) based on electron densities required to produce observed X-ray and Gamma-ray fluxes. They found $D_P \approx D_H$ (or $D_P$ slightly lower) is consistent with observations.
• Joint Fitting: Used Markov Chain Monte Carlo (MCMC) to simultaneously fit:
  1. The gamma-ray energy spectrum (Fermi-LAT + HAWC).
  2. The morphological profile of the halo (HAWC angular distribution).
  3. The diffusion coefficient at 100 TeV ($D_{100}$) measured by HAWC.
• Parameters Fitted: Break energy ($E_c$), upstream scale ($u_1a$), diffusion indices ($\delta, \delta_1$), spectral slope ($\sigma$), and conversion efficiency ($\epsilon$).

3. Key Contributions

• Unified Validation Method: Proposed a rigorous framework to test DSA and diffusion theories simultaneously by requiring a single set of parameters to explain both the energy spectrum (acceleration physics) and the spatial morphology (propagation physics) of the Geminga halo.
• Bounded Upstream Model: Introduced a more realistic treatment of the finite upstream region of the pulsar wind termination shock, deriving specific spectral cutoffs dependent on the diffusion coefficient.
• Multi-wavelength Synthesis: Successfully integrated X-ray (PWN structure), GeV (Fermi-LAT), and TeV (HAWC) data to constrain the diffusion coefficient across a wide energy range (10 GeV to >100 TeV).

4. Key Results

• Consistency with DSA: The theoretical models (DSA + Diffusion) are consistent with experimental observations. The best-fit parameters yield a gamma-ray spectrum that matches Fermi-LAT and HAWC data well.
• Diffusion Coefficient Behavior:
  • The diffusion coefficient $D(E)$ exhibits a rapid increase above a break energy of approximately 100 TeV.
  • Below 100 TeV, diffusion is slow (suppressed), consistent with the "slow-diffusion zone" hypothesis.
  • Above 100 TeV, $D(E)$ rises sharply, approaching the typical Galactic interstellar medium (ISM) diffusion value.
  • The diffusion index $\delta$ in the high-energy regime is significantly steeper than the standard Kolmogorov index ($\delta \approx 13-23$ vs. 1.33), supporting the SNR-induced turbulence model where turbulence cannot be excited for wave vectors smaller than the electron Larmor radius.
• Boundary Conditions: The zero-streaming condition provided a better fit (lower $\chi^2$) than the zero-particle condition, suggesting electrons are reflected/cooled in the far upstream rather than captured.
• Morphology: The model successfully reproduces the observed angular profile of the Geminga halo in HAWC energy bins (17.8–31.6 TeV and 31.6–316 TeV).
• Magnetic Field: The analysis implies a magnetic field strength of $\sim 20 \mu G$ in the Geminga PWN. If the field is stronger ($\sim 44 \mu G$), the diffusion coefficients in the PWN and halo would be nearly identical.

5. Significance and Future Outlook

• Validation of Theory: This work provides strong evidence supporting the DSA mechanism and the SNR-induced diffusion model for cosmic rays in the sub-PeV range, elevating direct testing from MeV to sub-PeV energies.
• Model Specificity: The confirmation specifically validates the acceleration and propagation theory; it does not rule out other acceleration mechanisms (like magnetic reconnection), but similar tests would be valuable for those models.
• Future Prospects:
  • Current limitations include wide energy bins in morphological data, preventing high-precision tests of the energy-dependent diffusion coefficient.
  • Future observations by LHAASO-KM2A and improved HAWC data are expected to provide higher precision, particularly to confirm the rapid increase of the diffusion coefficient above 100 TeV.
  • The findings suggest that a "cutoff" or break in the diffusion coefficient may be a universal feature of high-energy particle acceleration sources (PeVatrons), with the break energy scaling with the source's maximum acceleration capability.

In conclusion, the paper establishes a robust methodology for using the Geminga system as a cosmic laboratory, demonstrating that the standard DSA and diffusion models, when applied with a bounded upstream and energy-dependent diffusion, successfully explain the complex multi-wavelength phenomenology of the Geminga halo.