pith. sign in

arxiv: 2606.12390 · v1 · pith:DT7XO2CVnew · submitted 2026-06-10 · 🌌 astro-ph.HE

Time-dependent cosmic-ray escape from wind bubbles: hard spectra formation

Lukas Merten , Sophie Aerdker , Enrico Peretti This is my paper

Pith reviewed 2026-06-27 08:35 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords cosmic rayswind bubblesdiffusive shock accelerationtime-dependent transportescaping spectraturbulence modelsstochastic differential equations
0
0 comments X

The pith

Wind bubbles produce cosmic-ray escape spectra harder than the standard E to the minus two during the wind-driven phase.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper solves the time-dependent spherical transport equation for cosmic rays injected at the termination shock of a wind bubble and escaping at the moving forward shock, which is modeled as a free boundary. It finds that the escaping spectra can harden beyond the E to the minus two slope expected from diffusive shock acceleration because the escape boundary expands outward over time. A reader would care because this changes the expected output from these systems and could affect multi-messenger signals and the grammage particles accumulate inside the bubble. The hardening occurs specifically in the wind-driven phase and is accompanied by possible strong suppression at the lowest energies that depends on the choice of turbulence model. The calculation uses stochastic differential equations to integrate the advection and diffusion processes.

Core claim

During the wind driven phase, the downstream escaping spectra from wind bubbles can be harder than ∼E^{-2}, the conventional expectation from diffusive shock acceleration, because cosmic rays are continuously injected at the termination shock and propagate by advection and diffusion until they reach the time-dependent position of the forward shock treated as a free escape boundary.

What carries the argument

The time-dependent forward shock treated as a free escape boundary in a one-dimensional spherical advection-diffusion transport model solved via stochastic differential equations.

If this is right

  • The initial energy spectrum can be significantly suppressed at the lowest energies, depending on the turbulence model chosen.
  • Low-energy particles experience efficient confinement inside the bubble.
  • This confinement produces observable effects in multi-messenger radiation and the cosmic-ray grammage accumulated within the bubble.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same time-dependent escape boundary could alter spectra from other expanding shocks whose forward edges move outward on comparable timescales.
  • If wind bubbles contribute substantially to the galactic cosmic-ray population, their hardened output would shift the composite spectrum at certain energies.
  • Gamma-ray or neutrino observations targeting individual bubbles could directly test the predicted low-energy suppression feature.

Load-bearing premise

The forward shock functions purely as a free escape boundary whose position evolves with time, with no reflection or re-acceleration of particles occurring there.

What would settle it

A measurement of the escaping cosmic-ray spectrum from a wind bubble in its wind-driven phase that shows a spectral index of -2 or softer would contradict the reported hardening.

Figures

Figures reproduced from arXiv: 2606.12390 by Enrico Peretti, Lukas Merten, Sophie Aerdker.

Figure 1
Figure 1. Figure 1: Set up of the wind bubble model where particles are injected with a pre-accelerated spectrum at the termination shock (TS) and the forward shock (FS) forms is assumed as the free escape boundary. Both shocks propagate with different speeds so that the distance between them is changing over time. The radii of the termination and forward shocks expanding in a medium of density n0 over time can be estimated b… view at source ↗
Figure 2
Figure 2. Figure 2: shows exemplarily the time evolution of the escaping spectrum in a scenario of Bohm diffusion (δ = 1), constant downstream magnetic field (α = 0), and strong magnetization (ϵB = 0.1). It is clearly visible that low energy particles with (E < 1 TeV) are completely depleted and remain trapped in the downstream region between the two shocks. Above ∼ 10 TeV the escaping spectrum is slightly harder than the inj… view at source ↗
Figure 3
Figure 3. Figure 3: Same as fig. 2 but for Kraichnan diffusion (δ = 3/2). While the low energy depletion is stationary the maximal energy is still increasing over time [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Same as fig. 2 but for Kolmogorov diffusion (δ = 5/3). The maximal energy is strongly time dependent, but the low energy depletion has almost vanished. by roughly one order of magnitude when the radial spectral index α is reduced by one. Since the flux at the high energy cut off—which is independent of the downstream magnetic field—is almost constant, this results in harder spectra for weaker magnetic fiel… view at source ↗
Figure 5
Figure 5. Figure 5: Escaping spectra at t = 1 Myr for different radial dependencies of the downstream magnetic field strength α. The other model parameters are fixed to Bohm diffusion (δ = 1) and low magnetization (ϵB = 0.01). 4. DISCUSSION We conclude that there are physically motivated configurations for which the energy spectra of escaping particles are significantly altered from those injected at the termination shock, du… view at source ↗
Figure 6
Figure 6. Figure 6: Escaping spectra at t = 1 Myr for different diffusion models (δ). The other model parameters are fixed to a radial downstream magnetic field strength dependence of α = −1 and low magnetization (ϵB = 0.01). in the Galactic magnetic field. We assume isotropic diffusion, and in a first approach neglect the structured Galactic magnetic field, which can in principle impact the arrival times significantly (see e… view at source ↗
Figure 7
Figure 7. Figure 7: Expected spectrum of CRs observed at Earth for a wind bubble that is located at a distance of r = 3 kpc. For other diffusion scenarios, not shown here, the results look very similar. Where the energy dependent transport from the wind bubbles to Earth can lead to an additional depletion of low energy particles depending on the observation time and distance between source and Earth. It should be noted, that … view at source ↗
Figure 8
Figure 8. Figure 8: The positions of a purely advective particle are shown in comparison to the two shock positions (dotted lines). For comparison the time scales in the ”frozen” scenario are marked with dash-dotted lines. as discussed above (α = −1, δ = 3/2, and ϵB = 0.1). The number density is highest at the current termination shock position, which is where the CRs are continuously injected. The particle density decreases … view at source ↗
Figure 9
Figure 9. Figure 9: Approximation of the position of the fastest 1% percent of purely diffusive particles at different energies. This shows that low energy particle will not escape the wind bubble purely diffusive. Furthermore, the low energy cut off can be estimated from the intersection of the diffusive tracer position with the forward shock position. Aerdker, S., Habegger, R., Merten, L., Zweibel, E., & Becker Tjus, J. 202… view at source ↗
Figure 10
Figure 10. Figure 10: CR number density for a constant downstream magnetic field strength (α = 0), Kraichnan diffusion coefficient (δ = 3/2) and large magnetization (ϵB = 0.1). The time evolution is color coded with bright yellow corresponding to earlier and dark blue colors corresponding to later times. The positions of the termination RTS and forward RFS shock are marked with a solid and dashed line, respectively. 10 2 10 1 … view at source ↗
Figure 11
Figure 11. Figure 11: Escaping spectra at early times for a constant downstream magnetic field strength (α = 0). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Escaping spectra at early times for a constant downstream magnetic field strength (α = 0) including the cut off. From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization. 10 2 10 1 10 0 10 1 10 2 10 3 10 4 E [TeV] 10 0 10 1 10 2 10 3 d N/d E E 2 [a.u.] t=0.010 Myr t=0.015 Myr t=0.030 Myr t=0.100 Myr 10 2 10 1 10 0 10 1 10 2 10 3 1… view at source ↗
Figure 13
Figure 13. Figure 13: Escaping spectra at early times for a decreasing downstream magnetic field strength (α = −1). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization [PITH_FULL_IMAGE:figures/full_fig_p021_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Escaping spectra at late times for a decreasing downstream magnetic field strength (α = −1). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization. 10 2 10 1 10 0 10 1 10 2 10 3 10 4 E [TeV] 10 0 10 1 10 2 10 3 10 4 d N/d E E 2 [a.u.] t=0.1 Myr t=0.5 Myr t=1.0 Myr t=2.5 Myr t=5.0 Myr 10 2 10 1 10 0 10 1 10 2 10 3 10 4 E [TeV] 1… view at source ↗
Figure 15
Figure 15. Figure 15: Escaping spectra at late times for a constant downstream magnetic field strength (α = 0). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization [PITH_FULL_IMAGE:figures/full_fig_p022_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Escaping spectra at late times for a constant downstream magnetic field strength (α = 0) including the cut off. From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization. 10 2 10 1 10 0 10 1 10 2 10 3 10 4 E [TeV] 10 0 10 1 10 2 10 3 10 4 d N/d E E 2 [a.u.] t=0.1 Myr t=0.5 Myr t=1.0 Myr t=2.5 Myr t=5.0 Myr 10 2 10 1 10 0 10 1 10 2 … view at source ↗
Figure 17
Figure 17. Figure 17: Escaping spectra at late times for a decreasing downstream magnetic field strength (α = −1). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.1 magnetization [PITH_FULL_IMAGE:figures/full_fig_p023_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Escaping spectra at late times for a decreasing downstream magnetic field strength (α = −1) cut off included. From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.1 magnetization. 10 2 10 1 10 0 10 1 10 2 10 3 10 4 E [TeV] 10 0 10 1 10 2 10 3 10 4 d N/d E E 2 [a.u.] t=0.1 Myr t=0.5 Myr t=1.0 Myr t=2.5 Myr t=5.0 Myr 10 2 10 1 10 0 10 1 10 2 10 … view at source ↗
Figure 19
Figure 19. Figure 19: Escaping spectra at late times for a strongly decreasing downstream magnetic field strength (α = −2). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization [PITH_FULL_IMAGE:figures/full_fig_p024_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Escaping spectra at late times for a strongly decreasing downstream magnetic field strength (α = −2). From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.01 magnetization. 0 20 40 60 80 100 120 140 160 r [pc] 10 27 10 26 10 25 10 24 10 23 10 22 10 21 10 20 10 19 10 18 N [a.u.] t=0.0 Myr t=0.5 Myr t=1.0 Myr t=2.5 Myr t=4.5 Myr 0 20 40 60 80 10… view at source ↗
Figure 21
Figure 21. Figure 21: Number density at late times for a constant downstream magnetic field strength (α = 0) including the cut off. From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.1 magnetization [PITH_FULL_IMAGE:figures/full_fig_p025_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: Number density at late times for a constant downstream magnetic field strength (α = −1) including the cut off. From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.1 magnetization. 0 20 40 60 80 100 120 140 160 r [pc] 10 27 10 26 10 25 10 24 10 23 10 22 10 21 10 20 10 19 10 18 N [a.u.] t=0.0 Myr t=0.5 Myr t=1.0 Myr t=2.5 Myr t=4.5 Myr 0 20 40 … view at source ↗
Figure 23
Figure 23. Figure 23: Number density at late times for a constant downstream magnetic field strength (α = −2) including the cut off. From left to right: Bohm, Kraichnan, and Kolmogorov diffusion. From top to bottom: Small ϵB = 0.01 and large ϵB = 0.1 magnetization [PITH_FULL_IMAGE:figures/full_fig_p026_23.png] view at source ↗
read the original abstract

Overview: Wind-driven bubbles are dynamic systems that can accelerate cosmic rays, depending on their physical properties, up to very high energies. We investigate how a time-dependent description of the particle transport may impact the escaping cosmic-ray flux. Model: The wind bubble system is modeled as spherically symmetric. Cosmic rays are continuously injected at the position of the termination shock and propagate through advection and diffusion until the escape at the time-dependent position of the forward shock, which is treated as a free escape boundary. Methods: The one-dimensional spherical time-dependent transport equation is solved by transforming it into the corresponding set of stochastic differential equations, and integrated with a modified version of the open source cosmic-ray propagation framework CRPropa. Results: We find that, during the wind driven phase, the downstream escaping spectra from wind bubbles can be harder than $\sim E^{-2}$, the conventional expectation from diffusive shock acceleration. Depending on the turbulence model the initial energy spectrum can be significantly suppressed at lowest energies, which could be an observable feature to distinguish between different turbulence realizations. This effect could lead to an efficient confinement of low energy particles, potentially leading to observable implication in terms of multi-messenger radiation and cosmic-ray accumulated grammage within the bubble.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript models cosmic-ray transport in spherically symmetric wind bubbles by solving the time-dependent 1D transport equation via stochastic differential equations in a modified CRPropa framework. Particles are injected continuously at the termination shock, advected and diffused, and escape only when they reach the outward-moving forward shock, which is imposed as a free-escape boundary. The central claim is that escaping spectra during the wind-driven phase can be harder than ~E^{-2}, with additional low-energy suppression that depends on the turbulence model and possible implications for multi-messenger signals and grammage.

Significance. If the numerical result is robust, the work supplies a concrete time-dependent mechanism capable of producing spectra harder than standard DSA expectations from wind bubbles, with potential observational signatures in low-energy suppression. The SDE implementation for evolving boundaries is a methodological asset that enables the time-dependent treatment.

major comments (1)
  1. [Model section] Model section (description of forward-shock boundary): the forward shock is treated strictly as a time-dependent absorbing free-escape boundary. This choice directly controls the residence-time distribution that produces the reported hardening; the manuscript does not test alternatives such as partial reflection, finite shock thickness, or re-acceleration at the forward shock. Because the skeptic note identifies this as the load-bearing assumption, a sensitivity study with at least one modified boundary condition is required to confirm that the E^{-2} violation survives.
minor comments (2)
  1. The abstract refers to 'downstream escaping spectra' while the geometry is spherically symmetric; a brief clarification of what 'downstream' denotes in this context would improve readability.
  2. The turbulence-model dependence is stated qualitatively; adding a short table or figure panel that quantifies the low-energy suppression for the two turbulence realizations would strengthen the claim of an observable discriminator.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive review and for recognizing the potential significance of the time-dependent mechanism and the SDE implementation. We address the single major comment below.

read point-by-point responses

  1. Referee: [Model section] Model section (description of forward-shock boundary): the forward shock is treated strictly as a time-dependent absorbing free-escape boundary. This choice directly controls the residence-time distribution that produces the reported hardening; the manuscript does not test alternatives such as partial reflection, finite shock thickness, or re-acceleration at the forward shock. Because the skeptic note identifies this as the load-bearing assumption, a sensitivity study with at least one modified boundary condition is required to confirm that the E^{-2} violation survives.

    Authors: We agree that the strictly absorbing free-escape boundary at the forward shock is a central modeling choice that shapes the residence-time distribution and the resulting spectral hardening. This boundary condition is physically motivated by the expectation that particles crossing the forward shock enter the interstellar medium and are no longer confined to the bubble. Nevertheless, we acknowledge that the robustness of the E^{-2} violation under alternative treatments has not been demonstrated. In the revised manuscript we will add a sensitivity study implementing at least one modified boundary condition (e.g., a partially reflecting boundary with a tunable reflection coefficient) and will show whether the reported hardening persists. The new results will be presented in an expanded Model section or dedicated appendix. revision: yes

Circularity Check

0 steps flagged

Numerical integration of transport equation produces emergent hardening; no circular reduction

full rationale

The paper solves the one-dimensional spherical time-dependent transport equation numerically by conversion to stochastic differential equations and integration in CRPropa. Cosmic-ray injection occurs at the termination shock and escape is imposed at the time-dependent forward-shock position treated as a free-escape boundary. The reported spectral indices harder than E^{-2} are outputs of this integration under the stated boundary conditions and advection-diffusion physics. No parameter is fitted to the target spectrum, no output is defined in terms of itself, and no load-bearing step reduces by construction to a self-citation or ansatz imported from prior work by the same authors. The model assumptions are explicit inputs; the hardening is an emergent numerical result rather than a definitional identity.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

Because only the abstract is available, the ledger is populated from the modeling choices explicitly named: spherical symmetry, continuous injection at the termination shock, free-escape outer boundary, and two unspecified turbulence models.

free parameters (1)
  • turbulence model parameters
    The abstract states that results depend on the choice of turbulence model; these parameters control low-energy suppression and are therefore free parameters of the simulation.
axioms (2)
  • domain assumption Spherically symmetric geometry
    The wind bubble system is modeled as spherically symmetric.
  • domain assumption Forward shock as free escape boundary
    Escape occurs at the time-dependent position of the forward shock treated as a free escape boundary.

pith-pipeline@v0.9.1-grok · 5747 in / 1344 out tokens · 12495 ms · 2026-06-27T08:35:16.595622+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

299 extracted references · 222 canonical work pages · 45 internal anchors

  1. [1]

    Science , year = 2005, month = sep, volume = 309, pages =

    Voyager 1 in the Foreshock, Termination Shock, and Heliosheath. Science , year = 2005, month = sep, volume = 309, pages =. doi:10.1126/science.1117569 , adsurl =

    work page doi:10.1126/science.1117569 2005

  2. [2]

    Physics of the Inner Heliosheath , year = 2006, series =

    Low-energy ions near the termination shock. Physics of the Inner Heliosheath , year = 2006, series =. doi:10.1063/1.2359308 , adsurl =

    work page doi:10.1063/1.2359308 2006

  3. [3]

    Science , year = 2005, month = sep, volume = 309, pages =

    Voyager 1 Explores the Termination Shock Region and the Heliosheath Beyond. Science , year = 2005, month = sep, volume = 309, pages =. doi:10.1126/science.1117684 , adsurl =

    work page doi:10.1126/science.1117684 2005

  4. [4]

    , year = 2008, month = jul, volume = 454, pages =

    An asymmetric solar wind termination shock. , year = 2008, month = jul, volume = 454, pages =. doi:10.1038/nature07022 , adsurl =

    work page doi:10.1038/nature07022 2008

  5. [5]

    Physics of the Inner Heliosheath , year = 2006, series =

    Termination shock particle spectral features. Physics of the Inner Heliosheath , year = 2006, series =. doi:10.1063/1.2359310 , adsurl =

    work page doi:10.1063/1.2359310 2006

  6. [6]

    , keywords =

    Modeling of the Heliospheric Interface, Magnetic Field, and Cosmic-Ray Transport. , keywords =. doi:10.1086/512848 , adsurl =

    work page doi:10.1086/512848

  7. [7]

    Advances in Space Research , year = 2009, month = may, volume = 43, pages =

    The acceleration of Anomalous Cosmic Rays by stochastic acceleration in the heliosheath. Advances in Space Research , year = 2009, month = may, volume = 43, pages =. doi:10.1016/j.asr.2009.02.010 , adsurl =

    work page doi:10.1016/j.asr.2009.02.010 2009

  8. [8]

    , archivePrefix = "arXiv", eprint =

    A Magnetic Reconnection Mechanism for the Generation of Anomalous Cosmic Rays. , archivePrefix = "arXiv", eprint =. doi:10.1088/0004-637X/709/2/963 , adsurl =

    work page doi:10.1088/0004-637x/709/2/963

  9. [9]

    , archivePrefix = "arXiv", eprint =

    A Model of Acceleration of Anomalous Cosmic Rays by Reconnection in the Heliosheath. , archivePrefix = "arXiv", eprint =. doi:10.1088/0004-637X/703/1/8 , adsurl =

    work page doi:10.1088/0004-637x/703/1/8

  10. [10]

    , keywords =

    The Theory of Anomalous Cosmic Rays. , keywords =

  11. [11]

    Astrophysics and Space Science Library , year = 2009, series =

    Nonlinear Cosmic Ray Diffusion Theories. Astrophysics and Space Science Library , year = 2009, series =. doi:10.1007/978-3-642-00309-7 , adsurl =

    work page doi:10.1007/978-3-642-00309-7 2009

  12. [12]

    Amsterdam: North-Holland, 1990, edited by Ginzburg, V.L

    Astrophysics of cosmic rays. Amsterdam: North-Holland, 1990, edited by Ginzburg, V.L. , year = 1990, adsurl =

    1990

  13. [13]

    1997 , chapter=

    Complex Variables: Introduction and Applications , author=. 1997 , chapter=

    1997

  14. [14]

    , archivePrefix = "arXiv", eprint =

    Transport of Cosmic-Ray Protons in Intermittent Heliospheric Turbulence: Model and Simulations. , archivePrefix = "arXiv", eprint =. doi:10.1088/0004-637X/781/2/93 , adsurl =

    work page doi:10.1088/0004-637x/781/2/93

  15. [15]

    Anomalous Transport of High-energy Cosmic Rays in Galactic Superbubbles. I. Numerical Simulations. , keywords =. doi:10.1088/0004-637X/749/2/178 , adsurl =

    work page doi:10.1088/0004-637x/749/2/178

  16. [16]

    Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications. Phys. Rep. , year = 1990, volume = 195, pages =

    1990

  17. [17]

    Annals of Mathematical Statistics , year = 1958, volume = 29, pages =

    A Note on the Generation of Random Normal Deviates. Annals of Mathematical Statistics , year = 1958, volume = 29, pages =

    1958

  18. [18]

    , year = 2006, month = jan, volume = 439, pages =

    The scaling laws of human travel. , year = 2006, month = jan, volume = 439, pages =. doi:10.1038/nature04292 , adsurl =

    work page doi:10.1038/nature04292 2006

  19. [19]

    Journal of the American Statistical Association , volume =

    A Method for Simulating Stable Random Variables , author =. Journal of the American Statistical Association , volume =. 1976 , publisher =

    1976

  20. [20]

    Physics Letters A , year = 1998, month = feb, volume = 239, pages =

    A fractional diffusion equation to describe L \'e vy flights. Physics Letters A , year = 1998, month = feb, volume = 239, pages =. doi:10.1016/S0375-9601(97)00947-X , adsurl =

    work page doi:10.1016/s0375-9601(97)00947-x 1998

  21. [21]

    Soviet Journal of Experimental and Theoretical Physics , year = 1995, month = nov, volume = 81, pages =

    Stochastic transport and fractional derivatives. Soviet Journal of Experimental and Theoretical Physics , year = 1995, month = nov, volume = 81, pages =

    1995

  22. [22]

    , keywords =

    Anisotropic Three-Dimensional Focused Transport of Solar Energetic Particles in the Inner Heliosphere. , keywords =. doi:10.1088/0004-637X/709/2/912 , adsurl =

    work page doi:10.1088/0004-637x/709/2/912

  23. [23]

    , eprint =

    Anomalous transport and particle acceleration at shocks. , eprint =

  24. [24]

    , year = 2001, month = feb, volume = 95, pages =

    Anomalous Cosmic Rays: Messengers from the Outer Heliosphere. , year = 2001, month = feb, volume = 95, pages =

    2001

  25. [25]

    8th International Conference of Numerical Modeling of Space Plasma Flows (ASTRONUM 2013) , year = 2014, series =

    Numerical Modeling of Anomalous Transport of Cosmic Rays. 8th International Conference of Numerical Modeling of Space Plasma Flows (ASTRONUM 2013) , year = 2014, series =

    2013

  26. [26]

    Handbook of stochastic methods for physics, chemistry and the natural sciences

  27. [27]

    and De Fabritiis, G

    Gorenflo, R. and De Fabritiis, G. and Mainardi, F. Discrete random walk models for symmetric L\'evy-Feller diffusion processes. Physica A: Statistical Mechanics and its Applications. 1999

    1999

  28. [28]

    , keywords =

    Scaling Theory for Cross-Field Transport of Cosmic Rays in Turbulent Fields. , keywords =. doi:10.1088/0004-637X/711/2/997 , adsurl =

    work page doi:10.1088/0004-637x/711/2/997

  29. [29]

    Numerical Approximation of a Fractional-In-Space Diffusion Equation, I , author =

  30. [30]

    Physica A Statistical Mechanics and its Applications , eprint =

    Stochastic calculus for assets with non-Gaussian price fluctuations. Physica A Statistical Mechanics and its Applications , eprint =. doi:10.1016/S0378-4371(02)00803-8 , adsurl =

    work page doi:10.1016/s0378-4371(02)00803-8

  31. [31]

    Physics of Plasmas , year = 2014, month = feb, volume = 21, number = 2, eid =

    A gedankenexperiment for anomalous diffusion in a charge-fluctuating dusty plasma. Physics of Plasmas , year = 2014, month = feb, volume = 21, number = 2, eid =. doi:10.1063/1.4865408 , adsurl =

    work page doi:10.1063/1.4865408 2014

  32. [32]

    , keywords =

    Velocity Correlation and the Spatial Diffusion Coefficients of Cosmic Rays: Compound Diffusion. , keywords =. doi:10.1086/308492 , adsurl =

    work page doi:10.1086/308492

  33. [33]

    , eprint =

    Fractional Schr \"o dinger equation. , eprint =. doi:10.1103/PhysRevE.66.056108 , adsurl =

    work page doi:10.1103/physreve.66.056108

  34. [34]

    , archivePrefix = "arXiv", eprint =

    Superdiffusion of Cosmic Rays: Implications for Cosmic Ray Acceleration. , archivePrefix = "arXiv", eprint =. doi:10.1088/0004-637X/784/1/38 , adsurl =

    work page doi:10.1088/0004-637x/784/1/38

  35. [35]

    , keywords =

    Sub- and superdiffusive molecular displacement laws in disordered porous media probed by nuclear magnetic resonance. , keywords =. doi:10.1103/PhysRevE.74.066309 , adsurl =

    work page doi:10.1103/physreve.74.066309

  36. [36]

    , keywords =

    Non-Gaussian statistics and superdiffusion in a driven-dissipative dusty plasma. , keywords =. doi:10.1103/PhysRevE.78.046403 , adsurl =

    work page doi:10.1103/physreve.78.046403

  37. [37]

    , keywords =

    Competition between subdiffusion and L \'e vy flights: A Monte Carlo approach. , keywords =. doi:10.1103/PhysRevE.75.056702 , adsurl =

    work page doi:10.1103/physreve.75.056702

  38. [38]

    Journal of Computational and Applied Mathematics , keywords =

    Fox H functions in fractional diffusion. Journal of Computational and Applied Mathematics , keywords =

  39. [39]

    ACM Trans

    Matsumoto, Makoto and Nishimura, Takuji , title =. ACM Trans. Model. Comput. Simul. , issue_date =. 1998 , issn =. doi:10.1145/272991.272995 , acmid =

    work page doi:10.1145/272991.272995 1998

  40. [40]

    Journal of Computational and Applied Mathematics , keywords =

    Finite difference approximations for fractional advection-dispersion flow equations. Journal of Computational and Applied Mathematics , keywords =

  41. [41]

    Meerschaert and Charles Tadjeran

    Mark M. Meerschaert and Charles Tadjeran. Finite difference approximations for two-sided space-fractional partial differential equations. Applied Numerical Mathematics. 2006. doi:10.1016/j.apnum.2005.02.008

    work page doi:10.1016/j.apnum.2005.02.008 2006

  42. [42]

    , year = 2000, month = dec, volume = 339, pages =

    The random walk's guide to anomalous diffusion: a fractional dynamics approach. , year = 2000, month = dec, volume = 339, pages =. doi:10.1016/S0370-1573(00)00070-3 , adsurl =

    work page doi:10.1016/s0370-1573(00)00070-3 2000

  43. [43]

    Journal of Physics A Mathematical General , year = 2004, month = aug, volume = 37, pages =

    TOPICAL REVIEW: The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics. Journal of Physics A Mathematical General , year = 2004, month = aug, volume = 37, pages =. doi:10.1088/0305-4470/37/31/R01 , adsurl =

    work page doi:10.1088/0305-4470/37/31/r01 2004

  44. [44]

    EPL (Europhysics Letters) , keywords =

    Deriving fractional Fokker-Planck equations from a generalised master equation. EPL (Europhysics Letters) , keywords =. doi:10.1209/epl/i1999-00279-7 , adsurl =

    work page doi:10.1209/epl/i1999-00279-7

  45. [45]

    Planetary and Space Science , year = 1965, month = jan, volume = 13, pages =

    The passage of energetic charged particles through interplanetary space. Planetary and Space Science , year = 1965, month = jan, volume = 13, pages =. doi:10.1016/0032-0633(65)90131-5 , adsurl =

    work page doi:10.1016/0032-0633(65)90131-5 1965

  46. [46]

    , keywords =

    Evidence of Superdiffusive Transport of Electrons Accelerated at Interplanetary Shocks. , keywords =. doi:10.1086/525523 , adsurl =

    work page doi:10.1086/525523

  47. [47]

    Astrophysics and Space Sciences Transactions , year = 2008, month = jul, volume = 4, pages =

    Observations of anomalous transport of energetic electrons in the heliosphere. Astrophysics and Space Sciences Transactions , year = 2008, month = jul, volume = 4, pages =. doi:10.5194/astra-4-27-2008 , adsurl =

    work page doi:10.5194/astra-4-27-2008 2008

  48. [49]

    , keywords =

    Ion Superdiffusion at the Solar Wind Termination Shock. , keywords =. doi:10.1088/0004-637X/693/2/L118 , adsurl =

    work page doi:10.1088/0004-637x/693/2/l118

  49. [50]

    , keywords =

    Nonclassical Transport and Particle-Field Coupling: from Laboratory Plasmas to the Solar Wind. , keywords =. doi:10.1007/s11214-013-9966-9 , adsurl =

    work page doi:10.1007/s11214-013-9966-9

  50. [51]

    , eprint =

    Anomalous transport of cosmic ray electrons. , eprint =

  51. [52]

    , keywords =

    A new theory for perpendicular transport of cosmic rays. , keywords =. doi:10.1051/0004-6361:20077260 , adsurl =

    work page doi:10.1051/0004-6361:20077260

  52. [53]

    Journal of Physics A Mathematical General , year = 2007, month = sep, volume = 40, pages =

    Generalized compound transport of charged particles in turbulent magnetized plasmas. Journal of Physics A Mathematical General , year = 2007, month = sep, volume = 40, pages =. doi:10.1088/1751-8113/40/36/017 , adsurl =

    work page doi:10.1088/1751-8113/40/36/017 2007

  53. [54]

    Physics Today , keywords =

    Fractional Kinetics. Physics Today , keywords =. doi:10.1063/1.1535007 , adsurl =

    work page doi:10.1063/1.1535007

  54. [55]

    Fractional Calculus and Applied Analysis , volume=

    The space-fractional diffusion-advection equation: Analytical solutions and critical assessment of numerical solutions , author=. Fractional Calculus and Applied Analysis , volume=. 2014 , publisher=

    2014

  55. [56]

    , keywords =

    Sign for Super-diffusive Transport of Energetic Ions Associated with a Coronal-mass-ejection-driven Interplanetary Shock. , keywords =. doi:10.1088/2041-8205/731/2/L34 , adsurl =

    work page doi:10.1088/2041-8205/731/2/l34 2041

  56. [57]

    Journal of Computational Physics , year = 2007, month = jan, volume = 220, pages =

    A second-order accurate numerical method for the two-dimensional fractional diffusion equation. Journal of Computational Physics , year = 2007, month = jan, volume = 220, pages =. doi:10.1016/j.jcp.2006.05.030 , adsurl =

    work page doi:10.1016/j.jcp.2006.05.030 2007

  57. [58]

    , keywords =

    Quasi-ballistic and superdiffusive transport for impulsive solar particle events. , keywords =. doi:10.1051/0004-6361/201016278 , adsurl =

    work page doi:10.1051/0004-6361/201016278

  58. [59]

    Compound and Perpendicular Diffusion of Cosmic Rays and Random Walk of the Field Lines. I. Parallel Particle Transport Models. , keywords =. doi:10.1086/507415 , adsurl =

    work page doi:10.1086/507415

  59. [60]

    , keywords =

    From L \'e vy Walks to Superdiffusive Shock Acceleration. , keywords =. doi:10.1088/0004-637X/778/1/35 , adsurl =

    work page doi:10.1088/0004-637x/778/1/35

  60. [61]

    , keywords =

    Superdiffusive and Subdiffusive Transport of Energetic Particles in Solar Wind Anisotropic Magnetic Turbulence. , keywords =. doi:10.1086/502676 , adsurl =

    work page doi:10.1086/502676

  61. [62]

    Advances in Space Research , year = 2012, month = jun, volume = 49, pages =

    Anomalous particle transport in the heliosphere. Advances in Space Research , year = 2012, month = jun, volume = 49, pages =. doi:10.1016/j.asr.2011.10.022 , adsurl =

    work page doi:10.1016/j.asr.2011.10.022 2012

  62. [63]

    , keywords =

    Scale-invariant motion in intermittent chaotic systems. , keywords =. doi:10.1103/PhysRevE.47.851 , adsurl =

    work page doi:10.1103/physreve.47.851

  63. [64]

    , keywords =

    A Generalized Nonlinear Guiding Center Theory for the Collisionless Anomalous Perpendicular Diffusion of Cosmic Rays. , keywords =. doi:10.1088/0004-637X/716/1/671 , adsurl =

    work page doi:10.1088/0004-637x/716/1/671

  64. [65]

    , keywords =

    Modeling the Modulation of Galactic and Jovian Electrons by Stochastic Processes. , keywords =. doi:10.1088/0004-637X/735/2/83 , adsurl =

    work page doi:10.1088/0004-637x/735/2/83

  65. [66]

    Computer Physics Communications , year = 2012, month = mar, volume = 183, pages =

    A stochastic differential equation code for multidimensional Fokker-Planck type problems. Computer Physics Communications , year = 2012, month = mar, volume = 183, pages =. doi:10.1016/j.cpc.2011.11.014 , adsurl =

    work page doi:10.1016/j.cpc.2011.11.014 2012

  66. [67]

    , keywords =

    Time-Dependent Acceleration of Interstellar Pickup Ions at the Heliospheric Termination Shock Using a Focused Transport Approach. , keywords =. doi:10.1088/0004-637X/693/1/534 , adsurl =

    work page doi:10.1088/0004-637x/693/1/534

  67. [68]

    A Hitch-hiker's Guide to Stochastic Differential Equations: Solution Methods for Energetic Particle Transport in Space Physics and Astrophysics

    A Hitch-hiker's Guide to Stochastic Differential Equations. Solution Methods for Energetic Particle Transport in Space Physics and Astrophysics. , keywords =. doi:10.1007/s11214-017-0351-y , archivePrefix =. 1703.06192 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1007/s11214-017-0351-y

  68. [69]

    , keywords =

    Computation of cosmic-ray acceleration by Ito's stochastic differential equations. , keywords =

  69. [70]

    , keywords =

    A more accurate numerical scheme for diffusive shock acceleration. , keywords =. doi:10.1111/j.1365-2966.2010.17868.x , archivePrefix =. 1103.3049 , primaryClass =

    work page doi:10.1111/j.1365-2966.2010.17868.x 2010

  70. [71]

    Physics of Plasmas , year = 2022, volume =

    A nonlinear model of diffusive particle acceleration at a planar shock. Physics of Plasmas , year = 2022, volume =. doi:10.1063/5.0094667 , adsurl =

    work page doi:10.1063/5.0094667 2022

  71. [72]

    , keywords =

    Calculation of Diffusive Shock Acceleration of Charged Particles by Skew Brownian Motion. , keywords =. doi:10.1086/309429 , adsurl =

    work page doi:10.1086/309429

  72. [73]

    On the Origin of the Cosmic Radiation. Phys. Rev. , year = 1949, volume =. doi:10.1103/PhysRev.75.1169 , adsurl =

    work page doi:10.1103/physrev.75.1169 1949

  73. [74]

    , keywords =

    The structure of cosmic ray shocks. , keywords =

  74. [75]

    , keywords =

    The acceleration of cosmic rays in shock fronts - II. , keywords =. doi:10.1093/mnras/182.3.443 , adsurl =

    work page doi:10.1093/mnras/182.3.443

  75. [76]

    , keywords =

    The acceleration of cosmic rays in shock fronts - I. , keywords =. doi:10.1093/mnras/182.2.147 , adsurl =

    work page doi:10.1093/mnras/182.2.147

  76. [77]

    doi:10.1086/182658 , Journal =

    Particle acceleration by astrophysical shocks. , keywords =. doi:10.1086/182658 , adsurl =

    work page doi:10.1086/182658

  77. [78]

    , keywords =

    CRPropa 3.2 - an advanced framework for high-energy particle propagation in extragalactic and galactic spaces. , keywords =. doi:10.1088/1475-7516/2022/09/035 , archivePrefix =. 2208.00107 , primaryClass =

    work page doi:10.1088/1475-7516/2022/09/035 2022

  78. [79]

    Energy Spectrum of Fast Second Order Fermi Accelerators as Sources of Ultra-High-Energy Cosmic Rays

    Energy spectrum of fast second order Fermi accelerators as sources of ultra-high-energy cosmic rays. Astroparticle Physics , keywords =. doi:10.1016/j.astropartphys.2018.04.004 , archivePrefix =. 1612.03675 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.astropartphys.2018.04.004 2018

  79. [80]

    1951 , publisher=

    It. 1951 , publisher=

    1951

  80. [81]

    , keywords =

    Cosmic-Ray Spectra at Spherical Termination Shocks. , keywords =. doi:10.1086/375054 , adsurl =

    work page doi:10.1086/375054

Showing first 80 references.