Line-of-sight magnetic-field propagation effects on axion-like particle constraints from GRB 221009A

This paper demonstrates that while host-galaxy and Milky Way magnetic fields have only a mild impact on axion-like particle constraints derived from GRB 221009A, the properties of the intergalactic magnetic field constitute the dominant astrophysical uncertainty, significantly altering exclusion limits and introducing oscillatory features in the ALP mass-coupling plane.

The Gist

The Big Picture: A Cosmic Light Show and a Ghostly Particle

Imagine the universe is a giant, dark room filled with invisible fog. This fog is made of ancient light (the Extragalactic Background Light). When a super-bright flash of light (a Gamma-Ray Burst, or GRB) shoots through this room, the fog usually eats the light up. The brighter the light and the farther it travels, the more of it gets swallowed.

However, in October 2022, astronomers saw a flash called GRB 221009A that was so incredibly bright and energetic that it shouldn't have survived the trip. It was like seeing a candle flame from the other side of the galaxy without it getting extinguished.

This paper asks: How did the light survive?

One popular theory suggests the light didn't travel as light the whole way. Instead, it might have turned into a ghostly, invisible particle called an Axion-Like Particle (ALP), slipped through the fog untouched, and then turned back into light just before hitting Earth.

The authors of this paper wanted to test this theory. But to do that, they had to map out the "wind" (magnetic fields) the light traveled through, because the wind controls how easily the light can turn into a ghost and back again.

The Journey: Three Stops on the Road

The light from the explosion traveled through three distinct neighborhoods, each with its own "wind" (magnetic field):

1. The Host Galaxy (The Neighborhood): The galaxy where the explosion happened.
2. Intergalactic Space (The Open Ocean): The vast, empty space between galaxies.
3. The Milky Way (Our Backyard): The magnetic field of our own galaxy, right before the light hits Earth.

The Main Discovery: The "Open Ocean" is the Wild Card

The researchers ran simulations to see how different models of these magnetic fields would change the rules of the game.

• The Neighborhood and Backyard (Host Galaxy & Milky Way): They tried many different maps for the magnetic fields in these two areas. The result? It didn't matter much. Whether they used a rough map or a detailed map for these two spots, the answer about the ghost particles stayed roughly the same.

  • Analogy: Imagine driving through two towns. Whether the street signs are slightly crooked or perfectly straight, you still get to the highway in about the same way.

• The Open Ocean (Intergalactic Magnetic Field): This is where things got wild. The magnetic field in the empty space between galaxies is very poorly understood. It's like trying to drive across an ocean where the currents change direction randomly and nobody knows how strong they are.

  • Analogy: Imagine trying to sail across the ocean. If the currents are weak and steady, you can predict your path. But if the currents are chaotic, strong, and unpredictable, your destination changes completely depending on how you model the water.

The paper's main conclusion: The uncertainty in the "Open Ocean" (Intergalactic Magnetic Field) is the biggest problem. If we don't understand this magnetic field better, we can't be sure if the "ghost particle" theory is actually true or if we are just guessing.

The Three Zones of the Ghost Particle

The paper found that the behavior of these ghost particles depends heavily on their "weight" (mass). They identified three distinct zones:

1. The Heavy Zone (High Mass):

   • If the ghost particles are too heavy, they are too stubborn to change shape. They stay as light the whole time.
   • Result: The "fog" eats them up, just like normal physics predicts. The magnetic fields don't help them survive.

2. The Middle Zone (Medium Mass):

   • This is the "Goldilocks" zone. The particles are light enough to turn into ghosts, but heavy enough that the process is sensitive to the "wind."
   • Result: The survival rate starts to wiggle up and down like a heartbeat. The graph of the results looks like a jagged, oscillating line. It's chaotic and depends entirely on the specific details of the magnetic fields in the "Open Ocean."

3. The Light Zone (Low Mass):

   • If the particles are extremely light, they are very easy to turn into ghosts. They slip through the fog effortlessly.
   • Result: The "wiggles" disappear, and the graph becomes smooth. However, because they survive so well, the rules become very strict: if the ghost particles existed, they would have to be very weakly connected to light, or we would have seen even more light than we did.

Why This Matters

The authors aren't saying ghost particles definitely exist or don't exist. They are saying: "We cannot draw a final conclusion yet because we don't know enough about the magnetic fields in deep space."

They used a specific mathematical tool (a "minimum-χ2 fit") to compare their predictions with the actual data from the LHAASO telescope. They found that while the magnetic fields in our own galaxy and the host galaxy are manageable, the Intergalactic Magnetic Field is the dominant source of error.

The Takeaway

To solve the mystery of the super-bright GRB 221009A and prove or disprove the existence of these ghostly particles, scientists need to stop guessing about the magnetic currents in the deep, empty space between galaxies. Until we map that "Open Ocean" better, the map to the truth remains incomplete.

Technical Summary

Technical Summary: Line-of-sight magnetic-field propagation effects on axion-like particle constraints from GRB 221009A

Problem Statement
The detection of over 5,000 very-high-energy (VHE) photons from GRB 221009A by LHAASO, with energies extending up to 18 TeV, presents a tension with standard astrophysical expectations. Conventional models predict strong attenuation of such photons via pair production ($\gamma\gamma \to e^+e^-$) on the extragalactic background light (EBL). While this tension could stem from systematic uncertainties or limitations in EBL modeling, it also motivates searches for new physics, specifically axion-like particles (ALPs). ALPs can enhance the transparency of the universe by converting into photons in host galaxies, traversing intergalactic space as ALPs (evading EBL attenuation), and reconverting in the Milky Way. However, the efficiency of this photon-ALP oscillation mechanism is highly sensitive to the magnetic field environments along the line of sight. Previous studies utilizing GRB 221009A to constrain ALP parameters ($m_a$, $g_{a\gamma\gamma}$) have not systematically examined how uncertainties in the host-galaxy magnetic field (HGMF), intergalactic magnetic field (IGMF), and Milky Way magnetic field (MWMF) affect the derived exclusion limits.

Methodology
The authors revisit ALP constraints using the LHAASO spectrum of GRB 221009A (specifically the time interval $T_0 + 230$ s to $T_0 + 300$ s, covering 0.23–11.61 TeV). The analysis employs a rigorous propagation framework:

1. Formalism: The beam evolution is described by a Schrödinger-like equation for the photon-ALP state vector, incorporating the mixing matrix which accounts for ALP mass, plasma frequency, QED vacuum polarization, CMB dispersion, and photon absorption. The evolution through multiple magnetized domains is calculated via the ordered product of transfer matrices.
2. Magnetic Field Models:
   • Host Galaxy (HGMF): Three models are tested: a simplified homogeneous field, a phenomenological proxy mimicking the Milky Way, and a realistic model (Model 3) based on HST observations of the host galaxy and NGC 891-like simulations.
   • Milky Way (MWMF): The analysis compares the Jansson–Farrar (JF12, JF12b, JF12c) models against the newer Unger–Farrar (UF) models, selecting the UF "base" model as a benchmark.
   • Intergalactic Space (IGMF): The IGMF is modeled as a sequence of randomly oriented magnetized domains interspersed with voids. Two realizations are considered: the Domain-Like Sharp-Edges (DLSHE) model and the Domain-Like Smooth-Edges (DLSME) model. The study varies the domain length ($L_{dom}$), filling fraction ($f$), and field strength ($B_{IG}$).
3. Statistical Analysis: Constraints are derived using a minimum-$\chi^2$ fit to the observed spectrum at the 95% confidence level. The intrinsic GRB spectrum is modeled as a log-parabolic function with nuisance parameters ($J_0, a, b$). The photon survival probability ($P_{ALP}$) is computed for various ALP parameters and magnetic field configurations to determine the upper limits on $g_{a\gamma\gamma}$ as a function of $m_a$.

Key Contributions and Results
The paper systematically quantifies the impact of magnetic field uncertainties on ALP exclusion contours in the $m_a$–$g_{a\gamma\gamma}$ plane:

• Host and Galactic Field Dependence: The choice of HGMF and MWMF models induces only moderate variations in the exclusion limits. In the low-mass region ($m_a \lesssim 3 \times 10^{-8}$ eV), limits are approximately mass-independent, varying by a factor of roughly 1.5 across different models. Model 3 (HGMF) and the UF base model (MWMF) are identified as robust benchmarks.
• Dominance of IGMF Uncertainty: The IGMF constitutes the dominant astrophysical uncertainty. Varying IGMF parameters leads to significant changes in the exclusion contours, including the generation of pronounced oscillatory features.
• Three Mass Regimes: The inclusion of the IGMF reveals three distinct regimes in the exclusion limits:
  1. Large Mass ($m_a \gtrsim 10^{-8}$ eV): Photon-ALP conversion in the IGMF is inefficient due to a short oscillation length ($l_{osc} \ll L_{dom}$). The constraints closely match those obtained without the IGMF ($B_{IG}=0$).
  2. Intermediate Mass ($10^{-10} \text{ eV} \lesssim m_a \lesssim 10^{-8}$ eV): The oscillation length becomes comparable to the domain size ($l_{osc} \sim L_{dom}$). Multi-domain propagation generates terms involving products of sine and cosine functions of the phase $\omega = \pi L_{dom}/l_{osc}$. This results in exclusion contours with characteristic oscillatory structures dependent on $m_a$ and $L_{dom}$.
  3. Small Mass ($m_a \lesssim 10^{-10}$ eV): The oscillation length becomes large and nearly independent of mass ($l_{osc} \gg L_{dom}$). The oscillatory features disappear, but the IGMF still significantly enhances the high-energy photon survival probability by allowing ALPs to bypass EBL attenuation. Consequently, the exclusion limits become stronger (lower $g_{a\gamma\gamma}$) than in the $B_{IG}=0$ case, with a smooth dependence on mass.
• Domain Length Sensitivity: The onset of the oscillatory region is directly tied to the IGMF coherence scale ($L_{dom}$). Increasing $L_{dom}$ shifts the transition to higher masses.

Significance
The paper concludes that realistic modeling of line-of-sight magnetic fields, particularly the poorly constrained IGMF, is essential for deriving robust ALP constraints from GRB 221009A. While host-galaxy and Galactic field uncertainties are subdominant, the stochastic properties of the IGMF (field strength, coherence scale, and filling fraction) can drastically alter the derived limits. The presence of oscillatory features in the exclusion contours serves as a visible imprint of the IGMF's stochastic nature. The authors emphasize that future gamma-ray searches for ALPs must account for these propagation effects to avoid misinterpreting astrophysical uncertainties as new physics or vice versa.