Anisotropic Diffusion in Pulsar Halos: Interpreting the asymmetric morphology of Geminga and Monogem halos measured by HAWC

This paper utilizes an anisotropic diffusion model to interpret the asymmetric morphologies of the Geminga and Monogem pulsar halos observed by HAWC, revealing distinct magnetic field orientations and a local coherence length of approximately 100 pc, thereby establishing halo morphology as a powerful diagnostic tool for interstellar magnetic turbulence.

The Gist

Imagine the universe is a vast, foggy ocean. In this ocean, there are lighthouses called pulsars. These lighthouses don't just shine light; they shoot out a massive stream of tiny, invisible particles (electrons and positrons) like a high-pressure hose.

Usually, when you spray water into a calm lake, it spreads out in a perfect circle. But when astronomers looked at two specific pulsars, Geminga and Monogem, they didn't see perfect circles. Instead, they saw strange, stretched-out, "comet-like" shapes. The glow around these pulsars was lopsided.

This paper is like a detective story trying to figure out why the glow is lopsided.

The Mystery: Why isn't it a circle?

For a long time, scientists thought the particles were just getting stuck in a thick, sticky soup around the pulsars, slowing them down equally in all directions. But that theory couldn't explain the weird shapes.

The authors of this paper propose a new idea: The "Magnetic Wind" Theory.

Imagine the space around these pulsars isn't empty; it's filled with invisible magnetic "wind" or "currents" (magnetic fields).

• The Old Idea: Particles move randomly in all directions, like a drunk person stumbling in a crowd.
• The New Idea: Particles are like skiers. They can zoom very fast down a slope (along the magnetic field lines), but they are terrible at moving sideways across the slope. They get stuck trying to cross the "tracks."

The "Viewing Angle" Trick

Here is the clever part of the explanation.

Imagine you are looking at a long, straight hallway.

• If you look straight down the hallway, it looks like a tiny dot.
• If you look from the side, it looks like a long, stretched-out line.

The authors suggest that the "magnetic wind" around Geminga and Monogem is pointing almost directly at us (like looking down the hallway). Because the particles are zooming fast toward us and away from us, but moving very slowly sideways, the glow looks squashed and asymmetric from our perspective. It's an optical illusion caused by the angle we are looking at the magnetic field.

The Detective Work

The team used a supercomputer to simulate millions of different scenarios. They asked: "If the magnetic field is tilted at this angle, and the wind is this strong, does the resulting glow look like what the HAWC telescope actually saw?"

They found the "Goldilocks" settings:

1. The Angle: The magnetic fields aren't pointing straight at us, but they are tilted significantly (about 20–30 degrees). This explains the lopsided shape.
2. The "Wind" Strength: They measured something called the "Alfvénic Mach number" (a fancy way of saying how turbulent the magnetic wind is). They found it's relatively calm (about 0.2), meaning the particles are definitely struggling to move sideways.
3. The Neighborhood: They discovered that Geminga and Monogem are neighbors (about 100 light-years apart), but they live in different magnetic neighborhoods. The magnetic fields around them point in slightly different directions, like two houses on the same street having their front doors facing different ways.

Why Does This Matter?

Think of the interstellar magnetic field as the invisible skeleton of our galaxy. We can't see it directly, but we can see how it affects the "glow" of these pulsars.

By studying these weird shapes, the authors are essentially doing galactic MRI scans. They are mapping out the invisible magnetic currents of our galaxy. They found that these magnetic "neighborhoods" stay consistent for about 100 light-years before the direction changes.

The Big Takeaway

This paper tells us that the universe isn't just a random mess of particles. It's structured. The strange, lopsided shapes of the Geminga and Monogem halos aren't mistakes; they are clues.

They prove that:

1. Magnetic fields act like invisible rails, guiding particles in specific directions.
2. We are looking at these cosmic structures from a specific angle that makes them look stretched.
3. By understanding these shapes, we can finally start to map the invisible magnetic highways that crisscross our galaxy.

In short: The universe is playing a game of perspective, and these pulsars are the props that help us see the invisible magnetic skeleton holding it all together.

Technical Summary

Here is a detailed technical summary of the paper "Anisotropic Diffusion in Pulsar Halos: Interpreting the asymmetric morphology of Geminga and Monogem halos measured by HAWC" by Si-Zhe Wu, Chao-Ming Li, and Ruo-Yu Liu.

1. Problem Statement

Recent observations by the High Altitude Water Cherenkov (HAWC) Observatory and the Large High Altitude Air Shower Observatory (LHAASO) have revealed that the TeV gamma-ray halos surrounding the nearby pulsars Geminga and Monogem exhibit asymmetric morphologies.

• The Discrepancy: Standard isotropic diffusion models struggle to explain these asymmetries without invoking ad-hoc assumptions. Previous isotropic fits suggested a diffusion coefficient ($D$) significantly lower ($\sim 10^{28} \text{ cm}^2/\text{s}$) than the Galactic average, implying a highly turbulent local environment.
• The Hypothesis: The anisotropic diffusion model proposes that particle transport is governed by sub-Alfvénic turbulence where particles diffuse preferentially along magnetic field lines, while cross-field propagation is suppressed. If the mean magnetic field is not aligned with the observer's line of sight (LOS), this anisotropy projects onto the sky as an asymmetric halo shape.
• The Goal: This study aims to quantitatively constrain the properties of local interstellar magnetic turbulence (specifically the mean field orientation and Alfvénic Mach number) by modeling the asymmetric morphologies of Geminga and Monogem within the anisotropic diffusion framework.

2. Methodology

The authors employed a comprehensive modeling approach combining particle transport physics with statistical inference:

• Transport Equation: They solved the electron transport equation (Eq. 1) incorporating:
  • Anisotropic Diffusion: Defined by parallel ($D_\parallel$) and perpendicular ($D_\perp$) diffusion coefficients, related by the Alfvénic Mach number ($M_A$) via $D_\perp = M_A^4 D_\parallel$.
  • Energy Losses: Synchrotron and Inverse Compton (IC) scattering off the Cosmic Microwave Background (CMB) and Far-Infrared (FIR) fields, including Klein-Nishina corrections.
  • Injection: A time-dependent injection term following the pulsar's spin-down evolution, treated as a point source.
• Coordinate System: A specific coordinate system was defined where the $z$-axis aligns with the mean magnetic field. The geometry is parameterized by:
  • $\phi$: The inclination angle between the mean magnetic field and the LOS.
  • $\chi$: The azimuthal angle of the field projection relative to Right Ascension.
  • $M_A$: The Alfvénic Mach number.
• Comparison with Observations:
  • Instead of fitting raw surface brightness profiles directly, the authors simulated the 2D intensity map, divided it into four quadrants (matching HAWC's analysis), and calculated the effective diffusion coefficient ($D_{\text{eff}}$) or characteristic diffusion angle ($\theta_d$) for each quadrant.
  • They matched these model-derived $\theta_d$ values to the quadrant-dependent values reported by HAWC (A. Albert et al. 2024).
• Statistical Analysis: A Monte-Carlo Markov Chain (MCMC) method (using emcee) was used to scan the parameter space ($\phi, \chi, M_A$, and pulsar distances) to find the best-fit values and confidence intervals. Two integration radii ($r_{\text{max}} = 100$ pc and $200$ pc) were tested to assess sensitivity to the integration volume.

3. Key Contributions

• Quantitative Constraints on Turbulence: The study provides the first rigorous constraints on the Alfvénic Mach number ($M_A$) and mean magnetic field orientation for both Geminga and Monogem simultaneously within a unified anisotropic framework.
• Resolution of Statistical Tension: Previous studies suggested the magnetic field must be nearly parallel to the LOS ($\phi < 5^\circ$) to explain the lack of observed elongation, a statistically unlikely scenario. This work demonstrates that the observed asymmetry is consistent with larger inclination angles ($\phi \sim 19^\circ$ for Geminga, $\sim 29^\circ$ for Monogem), resolving the tension between non-detection of extreme elongation and the probability of viewing angles.
• Magnetic Coherence Length Estimation: By comparing the derived magnetic field orientations for the two pulsars (which are spatially separated by $\sim 100$ pc), the authors estimate the coherence length of the local interstellar magnetic field.

4. Key Results

• Alfvénic Mach Number ($M_A$): The fitted values for both pulsars are consistently around $M_A \approx 0.2$ (specifically $0.16 - 0.23$). This indicates sub-Alfvénic turbulence, confirming that cross-field diffusion is strongly suppressed, which explains the low effective diffusion coefficients observed.
• Magnetic Field Orientation:
  • Geminga: $\phi \approx 19^\circ$ (for $r_{\text{max}}=100$ pc).
  • Monogem: $\phi \approx 29^\circ$ (for $r_{\text{max}}=100$ pc).
  • The distinct orientations imply that Geminga and Monogem reside in different magnetic coherence regions.
• Coherence Length ($L_c$):
  • For $r_{\text{max}} = 100$ pc, the angle between the best-fit field directions is $\sim 29^\circ$ with minimal confidence contour overlap.
  • For $r_{\text{max}} = 200$ pc, the angle decreases to $\sim 10^\circ$ with partial overlap.
  • This variation suggests a local magnetic field coherence length of $L_c \sim 100$ pc. A shorter length ($\sim 10$ pc) would wash out asymmetry, while a much longer one ($>500$ pc) would result in identical field directions for both pulsars.
• Injection Efficiency ($\eta_e$): The study estimated the fraction of spin-down energy going into relativistic electrons, finding values consistent with previous literature, though sensitive to distance uncertainties.
• Robustness: The perpendicular diffusion coefficient ($D_\perp$) remains stable across different assumed parallel diffusion coefficients ($D_0$), reinforcing the reliability of the turbulence constraints.

5. Significance

• Diagnostic Tool: The paper establishes pulsar halo morphology as a powerful diagnostic tool for probing the microphysics of interstellar magnetic turbulence, specifically the transition from isotropic to anisotropic transport regimes.
• Validation of Anisotropic Model: The successful reproduction of the asymmetric morphologies without invoking extreme local turbulence amplification validates the anisotropic diffusion model as the leading explanation for TeV halos.
• Future Directions: The authors highlight that current single-coherence models are simplifications. Future high-precision morphological measurements are required to test multi-coherence scenarios (where magnetic domains vary stochastically) and to resolve the detailed structure of the interstellar magnetic field. This work sets the stage for more sophisticated modeling of electron distributions in complex magnetic environments.