🔭 astrophysics

Numerical investigation of particle acceleration at interplanetary shocks: diffusive and superdiffusive scenarios

This paper presents a numerical investigation demonstrating that incorporating superdiffusive transport into a first-order Fermi acceleration model for interplanetary shocks not only reproduces theoretical energy spectra but also successfully matches observed particle fluxes and accelerates particles to energies consistent with ACE spacecraft data more efficiently than standard diffusive models.

Original authors: Giuseppe Prete, Gaetano Zimbardo, Silvia Perri

Published 2026-04-17

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Giuseppe Prete, Gaetano Zimbardo, Silvia Perri

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 space between the planets not as an empty void, but as a bustling, chaotic highway filled with invisible particles (like protons and electrons) zipping around at incredible speeds. Sometimes, a massive "traffic jam" or a sudden wall of compressed gas, called a shock wave, ripples through this highway. These shocks are like cosmic accelerators, capable of boosting these tiny particles to energies that would make a particle collider on Earth jealous.

For decades, scientists have tried to figure out exactly how these particles get so fast. The traditional theory, called Diffusive Shock Acceleration (DSA), is like a game of "pinball" in a very predictable way. Imagine a ball bouncing between two moving walls. Every time it hits a wall, it gets a little push. In this old model, the ball bounces randomly but locally, taking small steps. It's a slow, steady climb up the energy ladder.

However, this new paper by Prete, Zimbardo, and Perri suggests that the universe might be playing a slightly different, more chaotic game. They propose that sometimes, these particles don't just take small steps; they take giant leaps.

The Two Ways Particles Move

To understand the difference, let's use two metaphors:

The Drunkard's Walk (Normal Diffusion): Imagine a person walking home from a bar. They are stumbling randomly, taking small, short steps in random directions. Sometimes they move toward home, sometimes away. It takes a long time to get anywhere because they keep getting stuck in local loops. This is the "old" theory.
The Super-Express (Superdiffusion): Now, imagine that same person, but occasionally, they get a magical boost that lets them teleport 10 miles forward in a single jump. They still stumble a bit, but those giant jumps allow them to cover ground much faster and return to the "starting line" (the shock wave) much more quickly. This is the Superdiffusion or Lévy Walk proposed in this paper.

What the Scientists Did

The researchers built a virtual simulation (a digital playground) to test these two ideas.

They created a digital shock wave.
They injected a crowd of "seed" particles (starting with a modest energy of 70 keV, which is like a low-energy runner).
They let these particles bounce back and forth across the shock. Every time a particle crossed the shock, it gained a little bit of speed (energy), just like a surfer catching a wave.
They ran two versions of the simulation: one where particles only took small steps (Normal Diffusion) and one where they were allowed to take giant, random leaps (Superdiffusion).

The Big Discovery

When they looked at the results, the "Super-Express" model (Superdiffusion) looked much more like reality than the "Drunkard's Walk."

Speed: The particles in the superdiffusive model got accelerated to high energies much faster. Because they could take those giant leaps, they didn't get swept away by the "wind" of the shock as easily. They could jump back upstream, hit the shock again, and get another boost.
The Shape of the Data: When the scientists compared their simulation results to real data from a spacecraft called ACE (which orbits Earth and watches the solar wind), the Superdiffusion model was a perfect match.
- The "Normal Diffusion" model predicted that particle numbers would drop off quickly and smoothly as you moved away from the shock (like a gentle hill).
- The Real Data (and the Superdiffusion model) showed a "power-law" drop-off. This is like a steep cliff that slowly flattens out. It means there are more high-energy particles hanging around far away from the shock than the old theory predicted.

Why Does This Matter?

Think of it like this: If you are trying to predict when a storm of dangerous radiation will hit a satellite or an astronaut, you need to know how fast those particles travel.

If you use the old model, you might think the storm is slow and predictable.
If you use the new model, you realize the storm can arrive faster and with more intense energy than expected because the particles are taking "express lanes" through space.

The Bottom Line

This paper confirms that in the chaotic environment of space, particles don't just shuffle along; they sometimes make giant, chaotic leaps. This "superdiffusive" behavior acts like a turbocharger for cosmic accelerators, allowing particles to reach dangerous energy levels much quicker than we previously thought.

By understanding this, scientists can build better models to predict space weather, protecting our satellites and future astronauts from the unexpected "boosts" of the universe. It turns out, in the cosmic highway, the fastest way to get somewhere isn't always a steady walk—it's sometimes a giant, random leap.

1. Problem Statement

Interplanetary (IP) shocks are primary accelerators of energetic particles in the heliosphere. While Diffusive Shock Acceleration (DSA) is the standard theoretical framework, it assumes normal (Brownian) diffusion where particle transport is characterized by a finite mean free path and a Gaussian distribution of step sizes. However, in-situ observations often reveal particle flux profiles that deviate from DSA predictions, such as:

Upstream fluxes that follow a power-law rise rather than an exponential rise.
Energy spectra that are "harder" (flatter) than the standard $E^{-2}$ or $E^{-3}$ predictions.
Acceleration timescales that appear faster than DSA allows.

These discrepancies suggest that anomalous transport, specifically superdiffusion, may play a critical role. Superdiffusion implies a mean square displacement $\langle \Delta r^2 \rangle \propto t^\alpha$ with $\alpha > 1$ , often modeled via Lévy walks where long free paths are possible. The paper aims to numerically investigate how incorporating superdiffusive transport into a first-order Fermi acceleration mechanism affects particle energy spectra and flux profiles, and to validate these findings against spacecraft observations.

2. Methodology

The authors developed an upgraded 1D numerical test-particle model based on previous work (Prete et al. 2019, 2021). The key features include:

Shock Geometry: An infinite planar shock moving through a plasma. The simulation box spans $L = \pm 10^{10}$ m.
Transport Mechanism: Particle motion is governed by a Langevin-type equation:
$d\mathbf{x}_i = (V_{bulk} + \mathbf{v}_{random})dt_i$
where $V_{bulk}$ is the plasma flow velocity (upstream $V_1$ or downstream $V_2$ ) and $\mathbf{v}_{random}$ represents scattering due to magnetic irregularities.
Scattering Models:
- Normal Diffusion (DSA): Scattering times ( $\tau_s$ ) are fixed, leading to Gaussian random walks.
- Superdiffusion (SSA): Scattering times follow a power-law distribution derived from Lévy walks. The probability density function for free path length $\ell$ and time $\tau_s$ couples space and time ( $\ell = v\tau_s$ ), ensuring a realistic physical description where long jumps require long travel times. The distribution is controlled by a power-law exponent $\mu$ (where $2 < \mu < 3$ for superdiffusion).
Acceleration Mechanism: Unlike previous models that only simulated transport, this version includes momentum gains at every shock crossing (First-order Fermi acceleration).
- Particles gain momentum when crossing from upstream to downstream and vice versa, calculated via Galilean transformations.
- Average momentum gain: $\langle \Delta p \rangle = \frac{2}{3} p \frac{V_1 - V_2}{v}$ .
Simulation Setup:
- Seed Population: Monoenergetic particles injected at $E_0 = 70$ keV.
- Energy Binning: Particles are binned into channels (100–200 keV up to >500 keV) to mimic ACE/LEMS-120 instrument data.
- Parameters: Simulations varied shock compression ratios ( $r=2, 3$ ), scattering times ( $\tau_s$ ), and superdiffusion parameters ( $\mu, \tau_0$ ).

3. Key Contributions

Integration of Acceleration and Anomalous Transport: The study successfully couples a first-order Fermi acceleration mechanism with Lévy walk transport, moving beyond pure transport studies to model the full acceleration process.
Theoretical Validation: The numerical results confirm theoretical predictions that superdiffusion modifies the spectral index $\gamma$ of the energy spectrum. Unlike DSA, where $\gamma$ depends only on the compression ratio $r$ , in SSA, $\gamma$ also depends on the superdiffusion exponent $\mu$ :
$\gamma_{SSA} = \frac{3}{r-1} \frac{\mu-2}{\mu-1} + 1$
Data-Driven Modeling: The authors derived specific model parameters directly from a real spacecraft event (ACE, Dec 14, 2006) to drive the simulations, rather than relying solely on theoretical guesses.

4. Results

A. Parametric Study

Trajectories: In superdiffusive simulations, particles exhibit long "jumps" allowing them to return to the shock from far upstream more frequently than in normal diffusion. This increases the probability of multiple shock crossings and accelerates the energization process.
Energy Spectra:
- Both DSA and SSA simulations converge to theoretical spectral slopes over time.
- Harder Spectra: Superdiffusion produces harder (flatter) energy spectra compared to DSA for the same compression ratio.
- Faster Acceleration: The convergence to the steady-state spectrum occurs more rapidly in the superdiffusive case.
- Compression Ratio: Higher compression ratios ( $r \sim 3$ ) lead to faster convergence and more efficient acceleration.
Spectral Index: The numerical spectral indices $\gamma$ matched the theoretical predictions derived from Eq. (14) across a range of $\mu$ values (2.2 to 2.9).

B. Case Study: ACE Shock (Dec 14, 2006)

The model was tested against observations of an interplanetary shock crossing by the ACE spacecraft at 1 AU.

Upstream Flux Profile:
- DSA Fit: Exponential decay ( $J \propto e^{-b|t|}$ ) yielded a poor fit ( $\tilde{\chi}^2 \approx 0.15-0.20$ ).
- SSA Fit: Power-law decay ( $J \propto |t|^{-\beta}$ ) yielded an excellent fit ( $\tilde{\chi}^2 \approx 0.025-0.029$ ).
- Derived Parameters: The best fit indicated superdiffusive transport with $\beta < 1$ , corresponding to $\mu \approx 2.66$ and $\tau_0 \approx 4$ s.
Flux Comparison:
- Diffusion Model: Failed to reproduce the observed upstream power-law rise and downstream profiles.
- Superdiffusion Model: Successfully reproduced the observed particle fluxes in both upstream and downstream regions across all energy channels.
Downstream Spectrum: The simulated energy spectrum using the derived parameters matched the ACE-measured differential flux ( $dJ/dE \propto E^{-\delta}$ ) remarkably well, confirming the "harder" spectrum predicted by superdiffusion.

5. Significance and Conclusions

Validation of Superdiffusion: The study provides strong numerical evidence that anomalous, superdiffusive transport is a viable and likely dominant mechanism for particle acceleration at interplanetary shocks, explaining observations that standard DSA cannot.
Acceleration Efficiency: Superdiffusion not only produces harder energy spectra but also speeds up the acceleration process, allowing particles to reach high energies more quickly than in normal diffusion scenarios.
Observational Consistency: The model successfully reproduces the specific flux profiles and spectral indices of the 2006 ACE shock event, suggesting that the "scale-free" nature of turbulence in the heliosphere (modeled by Lévy walks) is crucial for understanding energetic particle dynamics.
Future Work: The authors note that while the current model is simplified (1D, test-particle, fixed scattering times), it captures the essential physics. Future work will involve 3D models, energy-dependent scattering times, and self-generated turbulence, utilizing databases from missions like Solar Orbiter.

In summary, this paper demonstrates that incorporating Lévy walk statistics into shock acceleration models resolves discrepancies between standard theory and spacecraft observations, offering a more accurate description of how energetic particles are accelerated in the heliosphere.