Stochastic analysis of ultra-high energy cosmic ray interactions

This paper presents an analytic probabilistic description of ultra-high-energy cosmic ray interactions, modeling them as Markov jump processes to overcome the limitations of current deterministic and Monte Carlo approaches in accurately characterizing nuclear cascades and their impact on observed cosmic ray composition and spectrum.

The Gist

The Big Picture: Cosmic Ray "Pinball" in Space

Imagine the universe is a giant, dark pinball machine. The "balls" in this machine are Ultra-High-Energy Cosmic Rays (UHECRs)—atomic nuclei (like iron or helium) traveling at almost the speed of light.

As these balls zoom across the universe, they don't just fly in a straight line forever. They crash into a sea of invisible "bumpers" made of light (photons) from the Cosmic Microwave Background (the afterglow of the Big Bang) and infrared light.

When a cosmic ray hits these photons, it doesn't just bounce off; it often shatters. An iron nucleus might lose a piece, becoming a lighter element. That piece might shatter again, and so on. This is called a nuclear cascade.

The Problem: The "Roll of the Dice"

For a long time, scientists tried to predict what happens to these cosmic rays using a "smooth" approach. They treated the shattering like water flowing down a pipe: they assumed that if you have a million iron atoms, they would all lose pieces at the exact same rate, gradually turning into lighter elements in a predictable line.

But nature isn't smooth; it's choppy.

The paper argues that this "smooth" approach is wrong because the process is stochastic (random).

• The Analogy: Imagine you are walking through a minefield. A "smooth" model would say, "You will take 10 steps, then lose a shoe, then take 10 more steps, then lose a shoe."
• The Reality: In reality, you might take 100 steps and lose nothing, then suddenly step on a mine and lose a leg immediately. Or you might step on three mines in a row. The distance you travel before losing a piece is a roll of the dice.

Because of this randomness, you can't predict exactly what happens to one specific cosmic ray. You can only predict the probability of what happens.

The Solution: A New Mathematical Map

The authors, Leonel Morejon and Karl-Heinz Kampert, have created a new mathematical framework (a "map") to describe this randomness. Instead of guessing or running slow, expensive computer simulations (like rolling dice a million times to see the average), they found a way to write down the exact probability rules for these cosmic ray shatterings.

They treat the cosmic ray journey like a Markov Jump Process.

• The Analogy: Think of a board game where you move from square to square.
  • Square 1: Iron Nucleus.
  • Square 2: Iron minus one piece.
  • Square 3: Iron minus two pieces.
  • The Dice: The rules of the game (the probability of moving to the next square) depend on how fast you are going and what kind of "bumpers" (photons) you hit.

The paper provides the exact formula to calculate:

1. The Horizon: How far can a specific nucleus travel before it completely disintegrates?
2. The Spectrum: What energy will the pieces have when they arrive at Earth?
3. The Composition: If we see a Carbon atom on Earth, what was it originally? Was it Iron? Helium?

Key Discoveries in Plain English

1. The "Explosive" Disintegration

The paper shows that as cosmic rays travel further, they don't just slowly lose weight. They can suddenly "explode" into lighter pieces much faster than previously thought, especially when they hit certain types of light fields. This changes how far we think these particles can travel from their source.

2. The "Reverse Engine" (Working Backwards)

One of the coolest parts of the paper is the ability to work backward.

• The Analogy: Imagine finding a broken toy on the floor. Usually, you can't tell what it looked like before it broke. But with this new math, if you find a specific piece of a cosmic ray (like a Carbon atom) on Earth, you can calculate the likelihood that it came from a specific distance and was originally a specific heavy element (like Iron).
• This helps scientists figure out where these particles came from and what they were made of when they left their home galaxy.

3. The "Wobble" in the Sky

When cosmic rays travel, they are also pushed around by magnetic fields in space, making them wobble.

• The Old View: Scientists assumed the wobble was constant.
• The New View: Because the cosmic ray is constantly changing its identity (Iron $\to$ Silicon $\to$ Carbon), its "magnetic personality" changes too. A heavy piece wobbles differently than a light piece.
• The Result: The paper calculates exactly how much the arrival direction of these particles will be scrambled. This is crucial for trying to point a telescope back at the source of the cosmic ray.

Why Does This Matter?

• Better Telescopes: By understanding the randomness, we can better interpret data from giant cosmic ray detectors on Earth (like the Pierre Auger Observatory).
• Solving the Mystery: We still don't know exactly what accelerates these particles to such insane speeds. This new math helps narrow down the list of suspects (like exploding stars or black holes) by matching the "fingerprint" of the particles we see with what they would look like if they came from those sources.
• Efficiency: The authors also released a free software tool (called CRISP) that does these complex calculations instantly, saving scientists from running slow computer simulations.

Summary

This paper is like upgrading from a blurry, hand-drawn map of a stormy ocean to a high-definition, real-time GPS. It acknowledges that the journey of a cosmic ray is a chaotic, random game of chance, and it gives us the precise mathematical tools to understand the odds, predict the outcome, and trace the path back to the source.

Technical Summary

1. Problem Statement

Ultra-high-energy cosmic rays (UHECRs) interact with surrounding photon fields (e.g., CMB, IRB) during propagation, leading to energy loss and photodisintegration. While continuous energy losses (CEL) like pair production are well-described deterministically, stochastic losses (SL)—such as photodisintegration where nuclei lose a variable number of nucleons—require a probabilistic treatment.

Current approaches to modeling these interactions suffer from significant limitations:

• Continuous-Limit Approximations: Solve coupled ordinary differential equations (ODEs) treating stochastic losses as continuous. This neglects fluctuations and the discrete nature of nuclear fragmentation, leading to inaccuracies in predicting composition and horizons.
• Monte Carlo Simulations: While theoretically correct, these are computationally expensive, suffer from convergence issues, and offer limited theoretical insight into how input uncertainties (cross-sections, photon fields) propagate to final observables. They are often "blind" to the underlying probability space.

There is a lack of a formal analytic theoretical framework capable of describing the stochasticity of UHECR interactions and the resulting nuclear cascades with arbitrary precision.

2. Methodology

The authors develop an analytic probabilistic framework based on Markov jump processes to describe nuclear cascades.

• Mathematical Formulation:

  • The evolution of nuclear species is modeled as a continuous-time Markov chain.
  • The system is described using an interaction rate matrix, $\Lambda$, where diagonal elements represent total interaction rates (losses) and off-diagonal elements represent transition rates to daughter species.
  • The probability distribution of the distance traveled until a specific state (e.g., full disintegration) is reached is derived using the matrix exponential: $P(t) = e^{\Lambda t}$ (or $e^{\Lambda L}$ for distance).
  • The probability density function (PDF) for the distance to absorption is given by $f(L) = -\phi e^{\Lambda L} \Lambda \mathbf{1}$, where $\phi$ is the initial state vector and $\mathbf{1}$ is a vector of ones.

• Classification of Cascades:
  The paper categorizes cascades based on the structure of the interaction matrix and cross-section regularity:

  1. Regular Serial Cascades (RSeC): Assumes a single channel per step and interaction rates proportional to mass number ($\lambda_A \propto A$). This yields a closed-form solution related to the Beta distribution.
  2. Irregular Serial Cascades (ISeC): Accounts for deviations from mass scaling (due to nuclear structure effects like Giant Dipole Resonances) and spontaneous decays. The solution is a linear combination of exponentials (hypoexponential distribution).
  3. Concurrent Cascades (CoC): The most general case, allowing multiple disintegration channels (branching) at each step. This forms a network of intersecting paths, solved via the matrix exponential of the full interaction matrix.

• Handling Inhomogeneities:
  The framework distinguishes between two types of inhomogeneities caused by Continuous Energy Losses (CEL):

  • Coherent Inhomogeneities (CI): Where the boost evolution is independent of the specific cascade path (e.g., cosmological redshift scaling). These can be solved analytically by transforming the propagation distance into an effective "target thickness" ($\delta$).
  • Dispersive Inhomogeneities (DI): Where the boost evolution depends on the specific sequence of species (e.g., synchrotron losses dependent on $Z/A$). These are treated via quasi-homogeneous approximations or numerical integration of the Kolmogorov equation.

3. Key Contributions

• Analytic Framework: Establishes the first closed-form analytic description of UHECR stochastic interactions, linking them directly to Markov jump processes.
• Closed-Form Distributions: Derives exact probability distributions for:
  • Distance until full disintegration (horizons).
  • Occupation probabilities of intermediate nuclear species.
  • Angular deflections in extragalactic magnetic fields (EGMFs) accounting for changing rigidity due to disintegration.
• Reverse Propagation: Demonstrates how to reverse the Markov process to infer the likely source composition and distance of observed UHECRs (quasi-stationary states).
• Source Modeling: Provides a method to model time-dependent injection and escape mechanisms within sources (e.g., GRBs) using the same stochastic formalism used for propagation.
• Software Release: The authors released an open-source Python package, CRISP (Cosmic Ray Stochastic Interactions for Propagation), to implement these calculations.

4. Key Results

• Disintegration Horizons: The study shows that the "full disintegration horizon" (distance where 99% of nuclei are destroyed) is significantly different from the traditional energy loss length ($L_{EL}$). $L_{EL}$ often underestimates the horizon because it assumes average energy loss, ignoring the stochastic nature where some nuclei survive much longer.
• Mass Evolution Regularity:
  • In Regular Serial Cascades, the distance to disintegration scales logarithmically with mass ($\ln A$).
  • In Concurrent Cascades (realistic models with multiple channels), the relationship becomes linear with mass. This "disciplined disintegration" effect is more efficient than previously thought, shortening propagation scales.
• Impact of Boost and Redshift:
  • For high boosts ($\gamma \gtrsim 5 \times 10^9$), interactions are dominated by the CMB, and cosmological effects are negligible.
  • For lower boosts, IRB interactions dominate, and cosmological expansion (redshift) significantly alters the interaction rates and horizons. The paper provides analytic corrections for these effects.
• Angular Deflections: The paper derives that neglecting disintegration leads to incorrect estimates of arrival directions. Since lighter fragments are produced at different distances and have different rigidities ($R=E/Z$), the angular spread is a mixture of Rayleigh distributions. Identifying the specific observed species is crucial for associating UHECRs with their sources.
• Source Spectra: Modeling a GRB source shows that time-dependent injection and rigidity-dependent escape mechanisms (diffusive vs. advective) significantly alter the resulting UHECR spectrum and composition.

5. Significance

This work fundamentally shifts the paradigm of UHECR propagation modeling from purely numerical or deterministic approximations to a rigorous analytic probabilistic framework.

• Efficiency & Precision: It allows for the computation of observables with arbitrary precision at a fraction of the computational cost of Monte Carlo simulations.
• Uncertainty Quantification: By providing analytic expressions, it becomes possible to explicitly study correlations between input uncertainties (cross-sections, photon fields) and output observables, which is difficult with "black box" Monte Carlo methods.
• Astrophysical Insights: The framework reveals that the "disciplined disintegration" effect is a robust feature of stochastic cascades, implying that UHECR horizons are intrinsic properties of the interaction physics rather than just effective parameters.
• Multi-Messenger Astronomy: The ability to link source models (e.g., kilonovae/GRBs) directly to propagation effects allows for a unified approach to interpreting cosmic ray data alongside electromagnetic and neutrino observations.

In summary, the paper provides the mathematical tools necessary to accurately reconstruct the history of UHECRs from Earth-based observations, bridging the gap between source physics, propagation stochasticity, and detection.