{"id":"9cc25169-bb1c-40b9-8ad3-b2e5de5c7a65","arxiv_id":"2604.18703","paper_version":2,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"SAETASS solves the time-dependent astroparticle transport equation in spherical symmetry using a conservative finite-volume, operator-splitting framework, validated against analytical solutions and applied to cosmic-ray proton transport in a stellar-wind bubble.","lead":"SAETASS is an open-source Python solver for time-dependent cosmic-ray transport in 1D spherical symmetry, using operator splitting with MUSCL-Hancock advection and Crank-Nicolson diffusion. It gives astrophysicists a lightweight, validated alternative to large Galactic propagation codes for radially stratified systems like supernova remnants and stellar-wind bubbles.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Validation tests each operator in isolation; no evidence of a manufactured solution exercising all operators simultaneously with stiff coupling, leaving pre-equilibrium dynamics unverified.","rationale":"The reader correctly identified operator-splitting error in stiff coupled regimes as the load-bearing concern. My stress-test confirms this and sharpens it: the specific gap is the absence of a coupled manufactured solution that exercises all operators simultaneously. The isolated-operator validations (Sections 3.1–3.4, Appendix D) and the steady-state recovery in Section 4 do not constrain transient accuracy in the regime where the novel scientific claims (pre-equilibrium dynamics) live. However, this is a validation-completeness gap rather than evidence of incorrectness — the methods themselves (MUSCL-Hancock, Crank-Nicolson, Strang splitting) are standard and sound. The verdict remains CONDITIONAL because the tool appears well-constructed but the evidence for the transient-dynamics claim is incomplete. The reader's LOW confidence due to limited text availability is also justified — the critical Sections 2.3–2.6, 3, and 4 are not in the provided text, so this assessment is necessarily partial. If the full text reveals a coupled manufactured solution or a splitting-error convergence study, the concern would be substantially mitigated.","tokens_in":4267,"tokens_out":1888,"duration_ms":412926,"concrete_test":"Construct a manufactured solution for the full coupled transport equation (advection + diffusion + losses + source simultaneously active) in spherical symmetry with stiff parameter choices mimicking the Section 4 shock regime (e.g., u ~ 10^8 cm/s, D ~ 10^24 cm²/s/Bohm, loss timescale ~ 10^3 yr, domain ~ 10 pc). Drive all operators with a known analytic source term and verify that the numerical solution converges at the expected splitting order (2nd order for Strang) as the splitting timestep is refined. If the convergence rate degrades below 2nd order or the error exceeds the tolerance used for the isolated-operator tests, the pre-equilibrium dynamics in Section 4 are not quantitatively reliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim has two parts: (1) the solver is correct, and (2) it reveals physically meaningful pre-equilibrium dynamics. Part (1) is supported by Section 3 tests, but from the TOC and available text, these validate advection (3.1), diffusion (3.2), losses (3.3), and time dependence (3.4) separately. Appendix D (manufactured solutions) covers advection with space-time dependence (D.1) and diffusion with non-stationary coefficient (D.2) — again individually. There is no indication of a manufactured solution that simultaneously exercises advection + diffusion + losses + sources with cross-coupled stiff terms. The operator-splitting error is O(Δt²) for Strang splitting when each sub-problem is solved exactly, but when sub-stepping is used (Section 2.6) and the sub-problems are themselves discretized, the effective order and error magnitude depend on the ratio of timescales. In the astrophysical application (Section 4), near the termination shock, advection can be transonic while diffusion is simultaneously large (Bohm regime) and momentum losses are rapid — a regime where splitting errors are largest. The steady-state limit is recovered (which only tests that the fixed point of the split operator matches the continuous fixed point), but steady-state recovery does not constrain transient accuracy. The pre-equilibrium temporal dynamics — the novel scientific output — could therefore carry unquantified splitting errors. The reader's concern is essentially correct; the sharpening is that the gap is specifically the absence of a coupled manufactured solution, not just generic splitting risk.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"The manuscript presents SAETASS, an open-source numerical solver for the time-dependent astroparticle transport equation in 1D spherical symmetry. The solver uses a conservative finite-volume framework with operator splitting: radial advection and momentum losses are handled by a second-order MUSCL-Hancock scheme, diffusion by an implicit batched Crank-Nicolson algorithm, and sources are treated separately. The code is validated against analytical solutions for individual operators (advection, diffusion, losses) and against manufactured solutions for time-dependent cases (Appendix D). An astrophysical application to cosmic-ray proton transport in a stellar-wind bubble recovers known steady-state limits and presents pre-equilibrium temporal dynamics across Kolmogorov, Kraichnan, and Bohm diffusion regimes. The code is released as open-source.","tokens_in":4841,"tokens_out":780,"duration_ms":91805,"significance":"The paper addresses a genuine gap: general-purpose Galactic propagation codes are not optimized for spherically symmetric, radially stratified environments, and existing analytical treatments are limited to steady-state. A lightweight, open-source, modular solver with exact conservation properties and validated convergence rates is a useful community resource. The open-source release and the modular architecture (Appendix A) are commendable. The application to a real astrophysical scenario (Section 4) demonstrates practical utility beyond toy problems. The main scientific novelty — pre-equilibrium temporal dynamics — depends on the accuracy of the operator-splitting approach in coupled, stiff regimes, which is the central point requiring scrutiny.","major_comments":[{"comment":"The validation suite (Section 3 and Appendix D) tests each operator in isolation: advection (3.1, D.1), diffusion (3.2, D.2), losses (3.3), and time dependence (3.4). There is no manufactured solution or benchmark that simultaneously exercises advection + diffusion + losses + sources with cross-coupled terms. The astrophysical application (Section 4) involves a regime near the termination shock where advection is transonic, diffusion can be large (Bohm), and momentum losses are rapid — precisely where operator-splitting errors are largest. The recovery of the steady-state limit (Section 4) only confirms that the fixed point of the split operator matches the continuous fixed point; it does not constrain transient accuracy. Since the pre-equilibrium temporal dynamics are the novel scientific output, the absence of a coupled validation case is a load-bearing gap. The authors should add at a","section":null}],"minor_comments":[{"comment":"The abstract contains a semicolon where a comma is expected: 'expanding superbubbles; correctly understanding' should read 'expanding superbubbles, correctly understanding'.","section":null},{"comment":"The TOC entries for Sections 4 and 5 are garbled with replacement characters in the provided text; these should be verified before final submission.","section":null},{"comment":"The conclusion (Section 5) lists planned features including AMR, parallelization, and gamma-ray emission utilities. This is informative but could be shortened to keep the conclusion focused on what SAETASS currently supports.","section":null},{"comment":"References to specific equations and figures in Sections 3 and 4 could not be verified from the available text; the authors should ensure all cross-references are complete and point to the correct objects.","section":null}],"recommendation":"major_revision","confidential_remarks":"The reader's report and stress-test note both flag the same core issue — the absence of a coupled manufactured solution exercising all operators simultaneously. On inspection of the available text, this concern is well-founded and load-bearing for the paper's central novelty claim. The paper is otherwise a solid methods paper with a useful open-source contribution, and the issue should be addressable by adding one or two coupled test cases with quantified splitting error. I would encourage the authors rather than discourage them; the foundation is sound."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive report. The central concern — that the validation suite tests operators in isolation but lacks a coupled benchmark exercising advection + diffusion + losses + sources simultaneously — is well taken. We agree this is a load-bearing gap given that the pre-equilibrium temporal dynamics are the novel scientific output. We will add a coupled manufactured-solution test and discuss its implications for splitting-error control in the revised manuscript.","responses":[{"response":"The referee is correct on all counts. Our current validation suite (Section 3 and Appendix D) tests each operator independently, and the method-of-manufactured-solutions tests in Appendix D address time dependence but only for individual operators (advection in D.1, diffusion in D.2). We do not present a test that simultaneously exercises advection, diffusion, momentum losses, and sources with cross-coupled terms. The referee is also correct that the steady-state recovery in Section 4 constrains only the fixed point of the split scheme, not transient accuracy — and transient accuracy is precisely what the pre-equilibrium dynamics claim to demonstrate. This is a genuine gap in the manuscript as it stands. We will add a coupled manufactured-solution test in which we construct an analytical source term S(r,p,t) chosen so that the full transport equation with non-trivial advection velocity V(r,t), diffusion coefficient D(r,p,t), and momentum loss rate b(p) admits a known time-dependent solution f(r,p,t). This exercises all operators simultaneously and in cross-coupled fashion. We will report convergence rates and confirm that second-order accuracy is maintained in the coupled regime. We will also discuss the splitting timestep strategy (Section 2.6) in light of this test, explicitly addressing how the substepping control handles regimes where operator-splitting errors would otherwise be largest. The revised manuscript will include this as a new validation subsection (Section 3.5 or a new Appendix D entry) and we will reference it explicitly when discussing the reliability of the pre-equilibrium results in Section 4.","revision_made":"yes","referee_comment":"The validation suite tests each operator in isolation; there is no manufactured solution or benchmark that simultaneously exercises advection + diffusion + losses + sources with cross-coupled terms. The astrophysical application involves a regime near the termination shock where advection is transonic, diffusion can be large (Bohm), and momentum losses are rapid — precisely where operator-splitting errors are largest. The recovery of the steady-state limit only confirms that the fixed point of the split operator matches the continuous fixed point; it does not constrain transient accuracy. Since the pre-equilibrium temporal dynamics are the novel scientific output, the absence of a coupled validation case is a load-bearing gap. The authors should add a coupled validation case."}],"tokens_in":3882,"tokens_out":576,"duration_ms":90584,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"SAETASS is an open-source 1D spherical transport solver for cosmic rays — a tool paper, not a new physical result. The numerical methods (MUSCL-Hancock for advection/losses, implicit Crank-Nicolson for diffusion, Strang-type operator splitting) are all standard. What is new is the specific combination, the spherical geometry handling, and the open-source release. That is a legitimate contribution if the validation holds up.","headline":"Tool paper with solid individual-operator validation but a real gap in coupled testing","tokens_in":4987,"tokens_out":538,"would_cite":false,"duration_ms":19575,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"New solver cracks time-dependent cosmic-ray transport in spherical symmetry","keywords":["cosmic-ray transport","spherical symmetry","finite-volume method","operator splitting","MUSCL-Hancock","Crank-Nicolson","diffusion","astroparticle physics"],"falsifier":"If the operator-splitting errors in coupled regimes (e.g., simultaneous strong advection, diffusion, and rapid losses near a shock) prove larger than the substepping strategy can control, the transient pre-equilibrium dynamics would be quantitatively unreliable even when the steady-state limit is correct.","tokens_in":4366,"feed_emoji":"🌀","tokens_out":1045,"duration_ms":20902,"temperature":0.7,"pith_summary":"The paper introduces SAETASS, an open-source numerical solver for the fully time-dependent astroparticle transport equation in one-dimensional spherical symmetry. The core claim is that by combining a conservative finite-volume framework with operator splitting — a shock-capturing MUSCL-Hancock scheme for advection and momentum losses, and an implicit Crank-Nicolson method for diffusion — the solver achieves exact particle conservation and numerical stability even across steep gradients and discontinuities. The authors validate each operator against analytical solutions and manufactured solutions, then apply the tool to cosmic-ray proton transport in a wind-blown bubble, recovering known steady-state profiles while also exposing the transient pre-equilibrium dynamics that steady-state models cannot capture. The central object is the transport equation itself, solved in spherical coordinates with radial and momentum dependence, and the central mechanism is the modular splitting of physical operators with substepping control.","feed_headline":"New solver captures cosmic-ray transient dynamics in spherical systems","feed_subtitle":"SAETASS solves the full time-dependent transport equation in 1D spherical symmetry, revealing pre-equilibrium phases that steady-state tools","key_machinery":"Conservative finite-volume discretization in 1D spherical symmetry; operator splitting with substepping; MUSCL-Hancock scheme for advection and continuous momentum losses; implicit batched Crank-Nicolson for diffusion; method of manufactured solutions for time-dependent validation.","core_discovery":"SAETASS demonstrates that the time-dependent astroparticle transport equation in spherical symmetry can be solved accurately and conservatively by decomposing it into separately treated operators — hyperbolic advection and losses via MUSCL-Hancock, parabolic diffusion via implicit Crank-Nicolson — stitched together by a substepping strategy. Applied to a real astrophysical scenario, the solver reproduces established steady-state limits across Kolmogorov, Kraichnan, and Bohm diffusion regimes, while revealing energy-dependent transient phases of proton transport that static analytical models miss entirely.","pith_inferences":["If pre-equilibrium dynamics are physically significant, then gamma-ray or neutrino observations of young or rapidly evolving systems (e.g., young supernova remnants, freshly formed wind bubbles) may carry temporal signatures that steady-state models misattribute to spatial or spectral structure.","The operator-splitting approach could in principle be extended to include stochastic reacceleration or catastrophic losses (spallation, escape) as additional operators, though the splitting-error budget in regimes with stiff cross-coupling would need separate validation.","Comparing SAETASS time-dependent results against 3D codes in regimes where both apply would serve as a useful cross-check: agreement would validate the spherical-symmetry reduction, while discrepancies would quantify the cost of the symmetry assumption."],"forward_implications":["Time-dependent modelling of cosmic-ray transport in supernova remnants, stellar-wind bubbles, and pulsar wind nebulae is now accessible without the overhead of full 3D Galactic propagation codes.","Pre-equilibrium temporal dynamics — the transient behaviour before a system reaches steady state — can be quantified for radially stratified astrophysical environments, potentially altering interpretations of gamma-ray observations that assume steady state.","The open-source, modular architecture allows the community to extend the solver with additional physics such as hadronic and leptonic interactions, adaptive mesh refinement, or pre-computed evolving background profiles.","Energy-dependent transient phases across different diffusion regimes (Kolmogorov, Kraichnan, Bohm) can be compared directly within a single framework, enabling systematic study of how diffusion-law assumptions affect time-dependent signatures."],"fun_headline_variants":["Spherical transport solver reveals cosmic-ray transient phases","Operator-split solver captures cosmic-ray transients in spherical systems","New spherical solver exposes pre-equilibrium cosmic-ray transport dynamics","Time-dependent solver captures cosmic-ray dynamics in spherical symmetry","Modeling cosmic-ray transport in radial astrophysical systems with SAETASS"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The paper assumes that sequentially applying advection, diffusion, losses, and source operators — rather than solving them as a fully coupled system — introduces errors small enough to be controlled by substepping, even in regimes where multiple processes act stiffly and simultaneously, such as near a shock front with rapid momentum losses and strong diffusion.","fun_headline_variants_meta":{"raw":{"variants":["Spherical transport solver reveals cosmic-ray transient phases","Operator-split solver captures cosmic-ray transients in spherical systems","New spherical solver exposes pre-equilibrium cosmic-ray transport dynamics","Time-dependent solver captures cosmic-ray dynamics in spherical symmetry","Modeling cosmic-ray transport in radial astrophysical systems with SAETASS"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":1358,"prompt_tokens":609,"completion_tokens":749,"prompt_tokens_details":null},"tokens_in":609,"tokens_out":749,"duration_ms":34580,"temperature":1.0,"reasoning_tokens":809,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-05T10:33:40.902061+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the operator-splitting errors in coupled regimes (e.g., simultaneous strong advection, diffusion, and rapid losses near a shock) prove larger than the substepping strategy can control, the transient pre-equilibrium dynamics would be quantitatively unreliable even when the steady-state limit is correct.","supporting_citations":[],"review_version":2}