arxiv: 2604.18703 · v1 · submitted 2026-04-20 · 🌌 astro-ph.HE · astro-ph.IM

Recognition: unknown

SAETASS: Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry

Jos\'e Mar\'ia Garc\'ia-Morillo , Stefano Menchiari , Rub\'en L\'opez-Coto

Authors on Pith no claims yet

Pith reviewed 2026-05-10 03:32 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.IM

keywords astroparticle transportspherical symmetrynumerical solvercosmic ray propagationoperator splittingfinite volume methodsupernova remnantsdiffusion regimes

0 comments

The pith

SAETASS solves the time-dependent transport equation for astroparticles in spherical symmetry with a conservative finite-volume method.

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

The paper presents SAETASS as a lightweight numerical solver for non-thermal particle transport in radially stratified astrophysical settings such as supernova remnants and superbubbles. It uses a finite-volume discretization that enforces exact particle conservation together with an operator-splitting scheme separating advection and losses from diffusion. The code is validated on separate advection, diffusion, and loss problems and then applied to cosmic-ray proton propagation, where it reproduces known steady-state profiles while exposing pre-equilibrium time evolution under Kolmogorov, Kraichnan, and Bohm diffusion. This approach offers an efficient alternative to full three-dimensional galactic propagation codes when spherical symmetry is a reasonable approximation.

Core claim

SAETASS is a time-dependent solver built on a conservative finite-volume framework and a modular operator-splitting architecture; radial advection and continuous momentum losses are advanced with a second-order MUSCL-Hancock scheme while diffusion is integrated implicitly with a batched Crank-Nicolson method, allowing stable treatment of steep gradients, spatial discontinuities, and the regularity condition at the origin.

What carries the argument

The modular operator-splitting architecture that isolates advection and losses (treated by MUSCL-Hancock) from diffusion (treated by implicit Crank-Nicolson) to manage their combined effects while preserving particle number.

If this is right

The solver recovers established steady-state limits for cosmic-ray protons in spherical symmetry.
It reveals pre-equilibrium temporal dynamics that steady-state models miss across three standard diffusion regimes.
The same framework can be applied to other non-thermal particles or to expanding superbubbles with radially varying properties.
Validation tests confirm stability and conservation for pure advection, pure diffusion, and pure loss cases.

Where Pith is reading between the lines

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

The same splitting strategy could be combined with time-dependent magnetic-field evolution to study self-consistent particle transport.
Direct comparison of SAETASS output with selected three-dimensional codes would quantify the error introduced by the spherical-symmetry assumption.
Adding stochastic re-acceleration or spatially varying injection would extend the tool to more realistic acceleration sites.

Load-bearing premise

The operator-splitting procedure accurately combines advection, diffusion, and losses without introducing significant numerical artifacts when all processes act together.

What would settle it

A simulation that combines strong advection, diffusion, and losses in which the total particle number deviates from exact conservation or the solution fails to approach the known analytic steady-state profile.

read the original abstract

In order to model astrophysical environments characterized by radial stratification, such as supernova remnants or expanding superbubbles; correctly understanding the transport of non-thermal particles in astrophysical plasmas is essential. While large-scale Galactic propagation codes exist, they are often optimized for Cartesian or cylindrical geometries and lack the efficiency of one-dimensional spherically symmetric problems. In this work, we present SAETASS (Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry), a novel, open-source numerical tool designed to solve the time-dependent transport equation for astroparticles. The solver is built upon a conservative finite-volume framework that ensures exact particle conservation and numerical stability. To manage the interplay between diverse physical processes, SAETASS employs a modular operator-splitting architecture. Radial advection and continuous momentum losses are treated using a second-order, shock-capturing MUSCL-Hancock scheme, while the diffusive operator is integrated via an implicit, batched Crank-Nicolson algorithm. This approach allows for the robust handling of steep gradients, spatial discontinuities and regularity conditions at the origin. We rigorously validate the code through a suite of tests for pure advection, diffusion and losses. Finally, we demonstrate the solver's capabilities by modelling cosmic-ray proton transport in a real astrophysical scenario. Our results successfully recover established steady-state limits while revealing relevant pre-equilibrium temporal dynamics across Kolmogorov, Kraichnan and Bohm diffusion regimes. SAETASS provides the community with a lightweight, flexible tool for investigating particle acceleration and propagation in complex, radially dependent astrophysical environments.

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.

Desk Editor's Note private letter to a colleague

SAETASS is a practical new open-source 1D spherical solver for astroparticle transport that handles radial geometry cleanly, but its validation for the combined advection-diffusion-losses case is too thin to fully support the pre-equilibrium claims.

read the letter

SAETASS gives the community a lightweight, conservative finite-volume code built specifically for spherical symmetry. The MUSCL-Hancock treatment of advection and losses paired with implicit Crank-Nicolson diffusion is a sensible split for keeping conservation and stability while managing the origin regularity condition. That architecture is genuinely useful for supernova remnants and superbubbles, where full 3D or even 2D Cartesian runs are overkill. The separate tests for pure advection, pure diffusion, and pure losses are standard and reassuring, and recovering the known steady-state limits in the cosmic-ray proton example shows the code at least does not break the basics.

Referee Report

2 major / 2 minor

Summary. The paper presents SAETASS, an open-source finite-volume solver for the time-dependent astroparticle transport equation in spherical symmetry. It uses a conservative framework with MUSCL-Hancock treatment of radial advection and momentum losses combined with implicit Crank-Nicolson diffusion via operator splitting. The code is validated on separate pure-advection, pure-diffusion, and pure-losses tests, then applied to cosmic-ray proton transport in a radially stratified astrophysical scenario, where it recovers known steady-state limits and exhibits pre-equilibrium temporal evolution across Kolmogorov, Kraichnan, and Bohm diffusion regimes.

Significance. If the operator-splitting implementation proves accurate for simultaneous processes, SAETASS would supply a lightweight, particle-conserving 1-D spherical tool that fills a gap between full 3-D Galactic propagation codes and steady-state approximations. The modular architecture, exact conservation property, and demonstration of time-dependent dynamics in realistic environments constitute clear strengths for modeling supernova remnants and superbubbles.

major comments (2)

[Validation tests] Validation section: tests are reported only for pure advection, pure diffusion, and pure losses separately. No combined test with all three processes active simultaneously is shown against an analytic or high-fidelity reference solution in spherical geometry, which is required to substantiate the claim that pre-equilibrium dynamics are faithfully recovered without splitting artifacts when advection, diffusion, and losses compete on comparable timescales.
[Astrophysical demonstration] Astrophysical application section: the cosmic-ray proton demonstration recovers steady-state limits but presents the new pre-equilibrium temporal dynamics without quantitative error norms, convergence studies, or comparison to a reference solution that includes the full operator set in spherical coordinates; this weakens the assertion that the revealed dynamics are numerically reliable.

minor comments (2)

[Abstract] Abstract: the phrase 'rigorously validate the code through a suite of tests for pure advection, diffusion and losses' should be qualified to note that combined-process validation remains to be added.
[Numerical method] Notation: ensure the diffusion coefficient D(p,r) and its regime-specific forms (Kolmogorov, Kraichnan, Bohm) are defined with explicit radial dependence before first use in the results.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful and constructive review of our manuscript on SAETASS. Their comments on validation are well taken, and we address each major point below with the strongest honest response possible.

read point-by-point responses

Referee: Validation section: tests are reported only for pure advection, pure diffusion, and pure losses separately. No combined test with all three processes active simultaneously is shown against an analytic or high-fidelity reference solution in spherical geometry, which is required to substantiate the claim that pre-equilibrium dynamics are faithfully recovered without splitting artifacts when advection, diffusion, and losses compete on comparable timescales.

Authors: We acknowledge the value of a direct combined test. However, the modular operator-splitting design permits verification of each component in isolation, which is standard practice for such codes. The astrophysical demonstration activates all operators simultaneously and recovers the known analytic steady-state solution to high accuracy; any significant splitting artifact would prevent this recovery. We have revised the Validation section to include an explicit discussion of this indirect validation through steady-state recovery and added a grid-convergence study for a case with competing processes. revision: partial
Referee: Astrophysical application section: the cosmic-ray proton demonstration recovers steady-state limits but presents the new pre-equilibrium temporal dynamics without quantitative error norms, convergence studies, or comparison to a reference solution that includes the full operator set in spherical coordinates; this weakens the assertion that the revealed dynamics are numerically reliable.

Authors: The steady-state recovery is already compared quantitatively to the expected analytic spectrum. For the time-dependent pre-equilibrium phase, a full independent reference solution does not exist in the literature for this spherical setup. We have added convergence studies with respect to spatial resolution and time-step size in the revised manuscript, along with L2 error norms for the steady-state phase, confirming that the reported dynamics converge and remain consistent across the three diffusion regimes. revision: partial

Circularity Check

0 steps flagged

Numerical solver with separate-component validation; no derivation chain present

full rationale

The paper describes a finite-volume code using operator splitting for the spherical transport equation. Validation consists of independent tests on pure advection, pure diffusion, and pure losses, followed by a demonstration run on combined physics. No analytic derivations, fitted parameters renamed as predictions, or self-citation load-bearing steps appear in the provided text. The recovery of steady-state limits is a standard code-verification step against known analytic solutions for the separate operators, not a circular reduction. Pre-equilibrium dynamics are outputs of the time-dependent integration, not inputs. This is a standard software-methods paper whose central claims rest on numerical implementation and external benchmarks rather than self-referential logic.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are introduced; the work relies on standard numerical methods for transport equations.

pith-pipeline@v0.9.0 · 5597 in / 1112 out tokens · 52435 ms · 2026-05-10T03:32:08.299101+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

45 extracted references · 13 canonical work pages · 1 internal anchor

[1]

Strong, I.V

A.W. Strong, I.V. Moskalenko and V.S. Ptuskin,Cosmic-ray propagation and interactions in the galaxy,Annual Review of Nuclear and Particle Science57(2007) 285–327

2007
[2]

Grenier, J.H

I.A. Grenier, J.H. Black and A.W. Strong,The nine lives of cosmic rays in galaxies,Annual Review of Astronomy and Astrophysics53(2015) 199–246

2015
[3]

Parker,The passage of energetic charged particles through interplanetary space,Planetary and Space Science13(1965) 9

E. Parker,The passage of energetic charged particles through interplanetary space,Planetary and Space Science13(1965) 9

1965
[4]

Weaver, R

R. Weaver, R. McCray, J. Castor, P. Shapiro and R. Moore,Interstellar bubbles. II. Structure and evolution.,The Astrophysical Journal218(1977) 377

1977
[5]

Bykov, D.C

A.M. Bykov, D.C. Ellison, P.E. Gladilin and S.M. Osipov,Supernovae in compact star clusters as sources of high-energy cosmic rays and neutrinos,Advances in Space Research62(2018) 2764 [1706.01135]

work page arXiv 2018
[6]

Aharonian, R

F. Aharonian, R. Yang and E. de Oña Wilhelmi,Massive stars as major factories of galactic cosmic rays,Nature Astronomy3(2019) 561–567

2019
[7]

Morlino, P

G. Morlino, P. Blasi, E. Peretti and P. Cristofari,Particle acceleration in winds of star clusters, Monthly Notices of the Royal Astronomical Society504(2021) 6096

2021
[8]

Vieu and B

T. Vieu and B. Reville,Massive star cluster origin for the galactic cosmic ray population at very-high energies,MNRAS519(2023) 136 [2211.11625]

work page arXiv 2023
[9]

Strong and I.V

A.W. Strong and I.V. Moskalenko,Propagation of Cosmic-Ray Nucleons in the Galaxy,ApJ 509(1998) 212

1998
[10]

Evoli, D

C. Evoli, D. Gaggero, A. Vittino, G. Di Bernardo, M. Di Mauro, A. Ligorini et al.,Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients,JCAP2017 (2017) 015

2017
[11]

Kissmann,PICARD: A novel code for the Galactic Cosmic Ray propagation problem, Astroparticle Physics55(2014) 37 [1401.4035]

R. Kissmann,PICARD: A novel code for the Galactic Cosmic Ray propagation problem, Astroparticle Physics55(2014) 37 [1401.4035]. – 35 –

work page arXiv 2014
[12]

Maurin,USINE: Semi-analytical models for Galactic cosmic-ray propagation,Computer Physics Communications247(2020) 106942

D. Maurin,USINE: Semi-analytical models for Galactic cosmic-ray propagation,Computer Physics Communications247(2020) 106942

2020
[13]

Telezhinsky, V.V

I. Telezhinsky, V.V. Dwarkadas and M. Pohl,Time-dependent escape of cosmic rays from supernova remnants, and their interaction with dense media,A&A541(2012) A153

2012
[14]

Telezhinsky, V.V

I. Telezhinsky, V.V. Dwarkadas and M. Pohl,Acceleration of cosmic rays by young core-collapse supernova remnants,A&A552(2013) A102 [1211.3627]

work page arXiv 2013
[15]

Brose, I

R. Brose, I. Telezhinsky and M. Pohl,Transport of magnetic turbulence in supernova remnants, A&A593(2016) A20 [1606.04477]

work page arXiv 2016
[16]

Sushch, R

I. Sushch, R. Brose and M. Pohl,Modeling of the spatially resolved nonthermal emission from the Vela Jr. supernova remnant,A&A618(2018) A155 [1807.10551]

work page arXiv 2018
[17]

Brose, I

R. Brose, I. Sushch, M. Pohl, K.J. Luken, M.D. Filipović and R. Lin,Nonthermal emission from the reverse shock of the youngest Galactic supernova remnant G1.9+0.3,A&A627(2019) A166 [1906.02725]

work page arXiv 2019
[18]

Stone, K

J.M. Stone, K. Tomida, C.J. White and K.G. Felker,The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers,The Astrophysical Journal Supplement Series249(2020) 4

2020
[19]

García-Morillo, S

J.M. García-Morillo, S. Menchiari and R. López-Coto,SAETASS: Solver for Astroparticle Equation of Transport Analysis in Spherical Symmetry, Apr., 2026. 10.5281/zenodo.19468484

work page doi:10.5281/zenodo.19468484 2026
[20]

LeVeque,Finite Volume Methods for Hyperbolic Problems, Cambridge Texts in Applied Mathematics, Cambridge University Press (2002)

R.J. LeVeque,Finite Volume Methods for Hyperbolic Problems, Cambridge Texts in Applied Mathematics, Cambridge University Press (2002)

2002
[21]

Brenner and L.R

S.C. Brenner and L.R. Scott,The Mathematical Theory of Finite Element Methods, Springer New York (2008), 10.1007/978-0-387-75934-0

work page doi:10.1007/978-0-387-75934-0 2008
[22]

R.C. Chin,Dispersion and gibbs phenomenon associated with difference approximations to initial boundary-value problems for hyperbolic equations,Journal of Computational Physics18 (1975) 233

1975
[23]

McLachlan and G.R.W

R.I. McLachlan and G.R.W. Quispel,Splitting methods,Acta Numerica11(2002) 341–434

2002
[24]

Blanes, F

S. Blanes, F. Casas and A. Murua,Splitting methods for differential equations,Acta Numerica 33(2024) 1–161

2024
[25]

Toro,Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer, 3rd ed

E.F. Toro,Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer, 3rd ed. (2009), 10.1007/978-3-540-49834-6

work page doi:10.1007/978-3-540-49834-6 2009
[26]

LeVeque,Numerical Methods for Conservation Laws, Birkhäuser Basel, Basel (1990), 10.1007/978-3-0348-5116-9

R.J. LeVeque,Numerical Methods for Conservation Laws, Birkhäuser Basel, Basel (1990), 10.1007/978-3-0348-5116-9

work page doi:10.1007/978-3-0348-5116-9 1990
[27]

Lie and F

S. Lie and F. Engel,Theorie der Transformationsgruppen., American Mathematical Society (1970)

1970
[28]

Trotter,On the product of semi-groups of operators,Proceedings of the American Mathematical Society10(1959) 545

H.F. Trotter,On the product of semi-groups of operators,Proceedings of the American Mathematical Society10(1959) 545

1959
[29]

Strang,On the construction and comparison of difference schemes,SIAM Journal on Numerical Analysis5(1968) 506

G. Strang,On the construction and comparison of difference schemes,SIAM Journal on Numerical Analysis5(1968) 506

1968
[30]

Godunov,A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics,Matematicheskii Sbornik89(1959) 271

S.K. Godunov,A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics,Matematicheskii Sbornik89(1959) 271

1959
[31]

van Leer,Towards the ultimate conservative difference scheme

B. van Leer,Towards the ultimate conservative difference scheme. iv. a new approach to numerical convection,Journal of Computational Physics23(1977) 276

1977
[32]

van Leer,Towards the ultimate conservative difference scheme

B. van Leer,Towards the ultimate conservative difference scheme. v. a second-order sequel to godunov’s method,Journal of Computational Physics32(1979) 101. – 36 –

1979
[33]

van Leer,On the relation between the upwind-differencing schemes of godunov, engquist–osher and roe,SIAM Journal on Scientific and Statistical Computing5(1984) 1

B. van Leer,On the relation between the upwind-differencing schemes of godunov, engquist–osher and roe,SIAM Journal on Scientific and Statistical Computing5(1984) 1

1984
[34]

Harten,High resolution schemes for hyperbolic conservation laws,Journal of Computational Physics49(1983) 357

A. Harten,High resolution schemes for hyperbolic conservation laws,Journal of Computational Physics49(1983) 357

1983
[35]

Zhang and C.-W

X. Zhang and C.-W. Shu,On maximum-principle-satisfying high order schemes for scalar conservation laws,Journal of Computational Physics229(2010) 3091

2010
[36]

Crank and P

J. Crank and P. Nicolson,A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type,Mathematical Proceedings of the Cambridge Philosophical Society43(1947) 50–67

1947
[37]

LeVeque,Finite Difference Methods for Ordinary and Partial Differential Equations, Society for Industrial and Applied Mathematics (2007), 10.1137/1.9780898717839

R.J. LeVeque,Finite Difference Methods for Ordinary and Partial Differential Equations, Society for Industrial and Applied Mathematics (2007), 10.1137/1.9780898717839

work page doi:10.1137/1.9780898717839 2007
[38]

Thomas,Elliptic problems in linear difference equations over a network,Watson Sci

L.H. Thomas,Elliptic problems in linear difference equations over a network,Watson Sci. Comput. Lab. Rept., Columbia University, New York1(1949) 71

1949
[39]

Britz, O

D. Britz, O. Østerby and J. Strutwolf,Damping of crank–nicolson error oscillations, Computational Biology and Chemistry27(2003) 253

2003
[40]

Menchiari, G

S. Menchiari, G. Morlino, E. Amato, N. Bucciantini and M.T. Beltrán,Cygnus ob2 as a test case for particle acceleration in young massive star clusters,Astronomy & Astrophysics686 (2024) A242

2024
[41]

Blumenthal and R.J

G.R. Blumenthal and R.J. Gould,Bremsstrahlung, synchrotron radiation, and compton scattering of high-energy electrons traversing dilute gases,Rev. Mod. Phys.42(1970) 237

1970
[42]

Mannheim and R

K. Mannheim and R. Schlickeiser,Interactions of cosmic ray nuclei,A&A286(1994) 983 [astro-ph/9402042]

work page internal anchor Pith review arXiv 1994
[43]

Krakau and R

S. Krakau and R. Schlickeiser,Pion production momentum loss of cosmic ray hadrons,The Astrophysical Journal802(2015) 114

2015
[44]

Ginzburg,Theoretical Physics and Astrophysics, vol

V.L. Ginzburg,Theoretical Physics and Astrophysics, vol. 99 ofInternational Series in Natural Philosophy, Pergamon Press, Oxford (1979)

1979
[45]

Courant and D

R. Courant and D. Hilbert,Methods of mathematical physics. Volume 2, Wiley-VCH, Weinheim, 2. repr ed. (1989). – 37 –

1989