{"id":"601f701b-1a23-48e4-b2a7-b40b21c57325","arxiv_id":"2501.03066","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The inverse velocity dispersion observed by Parker Solar Probe in the 5 September 2022 SEP event is reproduced by a simulation where the weak flank of the CME shock gradually accelerates particles to higher energies.","lead":"A new study combines 3D shock modeling with particle transport simulations to explain why Parker Solar Probe saw delayed, slowly arriving energetic particles during a September 2022 solar storm. The delayed arrival and the strange energy pattern are attributed to the spacecraft's magnetic connection to a weak part of the expanding shock that took time to accelerate particles.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"IVD is largely built into the prescribed roll-over momentum (Eq. 3): pc rises with shock radius rs, so delayed high-energy onset at PSP is an input to the model; η=0.35 tuned to PSP makes the match a calibration check, not an independent test of flank acceleration.","rationale":"The paper's central claim is that PSP's inverse velocity dispersion results from gradual acceleration at a weak, strengthening shock flank, and that the PARADISE simulation reproduces the event. The strongest independent evidence would be a model that predicts, not prescribes, the evolution of the accelerated spectrum. Here, Eq. (3) prescribes the roll-over momentum as an increasing function of shock radius, which by construction delays the appearance of high-energy particles at any observer connected to an expanding shock. The IVD is therefore not a prediction of the model; it is a consequence of the chosen source parameterization. The one free parameter, η=0.35, is chosen to optimize PSP agreement, so the successful match reported in Section 3.3 is an exercise in calibration rather than validation. The Solar Orbiter comparison strengthens this concern: intensities are overestimated by about 10x and the spectrum is too hard, which the authors attribute to an overestimated apex shock speed, but which is also consistent with Eq. (3) giving too high a roll-over when the shock is strong. This does not invalidate the physical scenario—the shock is indeed weak at the PSP flank and strengthens over time, and the VDA independently supports a delayed release—but it means the quantitative reproduction cannot be used as decisive support for the mechanism. The proposed sensitivity test would show whether the IVD and onset timing are robust to η and to the rs-dependence; if they are, the conditioning concern is resolved. The reader's conditional verdict already captures this: the paper should be accepted only with the caveat that the source prescription is not independently tested. Hence no change in verdict.","tokens_in":25827,"tokens_out":6478,"duration_ms":63972,"concrete_test":"Run a sensitivity suite with the D⊥=0 PARADISE simulation: (1) vary η over {0.1, 0.2, 0.5, 1.0}; (2) as a control, replace Eq. (3) by a constant roll-over pc = η pinj (no rs dependence) and re-run. Record the PSP onset times and the energy at which the IVD appears. If the IVD transition energy shifts by more than a factor of 2 across η∈[0.2, 0.5], or if the IVD disappears when pc is held constant, the reproduction is calibration-driven and the attribution to slow flank acceleration loses support. Optionally, permit η to depend on local Mfm or ΘBN and check whether the Solar Orbiter discrepancy is reduced while the PSP match is retained.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the time-dependent acceleration at the PSP-connected shock flank, rather than transport or cross-field diffusion, causes the observed IVD. But the particle source in PARADISE is not a time-dependent acceleration solution. The shock distribution is the steady-state DSA spectrum of Eq. (2), with an exponential cutoff whose roll-over momentum pc is prescribed by Eq. (3) as a monotonic function of the shock radial distance rs (through rs/di) and of the local shock compression ratio. Because rs grows as the shock expands, pc increases with time at any given shock point. Consequently, any observer connected to an outward-moving flank will see low-energy particles arrive before high-energy particles; the IVD is essentially an input of the source model, not an emergent outcome. The only free normalization, η, is set to 0.35 to fit PSP observations (Section 3.3). Thus the reported reproduction of the PSP onset, spectrum, and IVD demonstrates that a suitable choice of cutoff evolution can match the data; it does not independently confirm that slow flank acceleration is the physical cause. This is highlighted by the Solar Orbiter comparison (Section 4): the same prescription overproduces intensities by an order of magnitude and gives a harder spectrum, indicating that Eq. (3) is not universally accurate. The constant-η assumption is also unverified, and the conclusion that cross-field diffusion is minor relies on comparing a D⊥≠0 run against a D⊥=0 run with the same η, without re-calibrating η for the diffusive run.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the 5 September 2022 SEP event observed by Parker Solar Probe (PSP) and Solar Orbiter. It combines multi-spacecraft EUV/coronagraph observations, a 3D ellipsoidal shock reconstruction, MAS-derived coronal MHD parameters, and PARADISE focused-transport simulations with a prescribed diffusive-shock-acceleration source spectrum. The simulation reproduces the PSP onset times (Table 1) and the observed inverse velocity dispersion (IVD) above about 1 MeV when a reduction factor eta = 0.35 is applied to the roll-over momentum in Eq. (3). A run with perpendicular diffusion yields worse onset-time agreement, leading the authors to conclude that cross-field diffusion was minor. The same model overestimates Solar Orbiter intensities by roughly an order of magnitude, and the authors attribute this to uncertainties in the roll-over momentum at the shock apex.","tokens_in":26184,"tokens_out":6219,"duration_ms":58664,"significance":"If the central interpretation is correct, the paper would provide a compelling multi-spacecraft case that the delayed SEP release and inverse velocity dispersion at PSP are governed by the time-evolving strength of the connected shock flank rather than by cross-field transport. The study is valuable as a detailed event analysis that integrates 3D shock reconstruction, MHD background modeling, and focused transport, and it is commendably transparent about its calibrated parameter and its disagreement at Solar Orbiter. However, because the delayed high-energy onset is effectively prescribed by Eq. (3) and eta is tuned to the PSP observations, the reported match does not independently confirm the proposed physical mechanism; the central claim therefore requires a careful reframing or additional sensitivity analysis.","major_comments":[{"comment":"The inverse velocity dispersion is not an emergent result of the simulation. In Eq. (3), the roll-over momentum pc is prescribed as a monotonically increasing function of the shock radius rs through rs/di, and the reduction factor eta = 0.35 is explicitly calibrated to match the PSP observations in the simulations presented below. Because pc sets the energy above which the source spectrum is suppressed, the delayed arrival of >1 MeV protons at PSP is built into the time-dependent source. The agreement in Table 1 therefore demonstrates that a suitable choice of cutoff evolution can reproduce the observations, but it does not independently establish that the IVD is caused by slow ongoing acceleration at the shock flank, as stated in the abstract.","section":"§3.3, Eq. (3), Fig. 8, Table 1"},{"comment":"The use of a single calibrated eta is undermined by the Solar Orbiter comparison. With eta = 0.35, the model overproduces Solar Orbiter intensities by roughly an order of magnitude and gives a harder spectrum, and the authors state in §3.3 that 'it is expected that eta varies across the shock surface.' The same paper therefore demonstrates that a constant eta is not accurate and that the PSP match may be specific to the calibrated value. The claim that the model is data-driven is weakened because the key free parameter is tuned to the target dataset, and the authors' own acknowledgment of the Solar Orbiter discrepancy should be given more weight in the interpretation.","section":"§4 and Fig. 11"},{"comment":"The conclusion that cross-field diffusion had a minor effect on this event is not strongly supported. The comparison between the D⊥ = 0 and D⊥ ≠ 0 runs uses the same source calibration, which was optimized for the D⊥ = 0 case; the worsening of the D⊥ ≠ 0 onset-time agreement is therefore not an independent test of transport. Moreover, the authors note that the shift of the IVD transition energy from about 1 MeV to 3 MeV may not be significant given the shock modeling uncertainties. A sensitivity study over the diffusion parameters, rather than a single diffusion on/off comparison, would be needed to justify the 'minor effect' conclusion.","section":"§3.3, Table 1, Fig. 10"}],"minor_comments":[{"comment":"The VDA-derived path length L ~ 6.3 R_sun is used to infer a release altitude of about 9 R_sun, but the VDA assumptions of scatter-free propagation and simultaneous release are not consistent with the extended, time-dependent source model used later in the paper; this interpretation should be framed more cautiously or cross-checked with the simulation.","section":"§2.2"},{"comment":"The cross-reference to 'Figure 10' when discussing the east flank's roll-over energy below the depicted energy channels appears to be a typo; it should likely refer to Figure 9.","section":"§3.3"},{"comment":"The normalization of the source spectrum would benefit from an explicit statement of the units and reference radius for np, so that the source amplitude in Eq. (2) is reproducible from the text.","section":"§3.3, Eq. (2)"},{"comment":"The Acknowledgments contain a sentence thanking 'the anonymous reviewer' for helping improve the paper; this sentence is inappropriate in a submitted manuscript and should be removed or rewritten.","section":"Acknowledgments"},{"comment":"The paper lacks an explicit data and code availability statement; given that the study relies on open-source tools such as PyThea, PARADISE, and MAS/EUHFORIA, a clear statement of the versions used and the availability of input data would improve reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"This is a competent and transparent case study, but the central physical interpretation overreaches the model's evidential support. The IVD is essentially an input through Eq. (3), and eta is calibrated to the very observations the model is claimed to reproduce. The authors should either substantially soften the causal attribution or add sensitivity analyses (e.g., varying eta, varying turbulence parameters) and a demonstration that the flank-acceleration scenario is distinguishable from alternative source evolution models. If the journal emphasizes physical interpretation over data-driven reproduction, these changes are necessary before acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid event study, and the first simulation-based attempt I know to tie the PSP inverse velocity dispersion to delayed acceleration at the expanding shock flank. It is not a clean confirmation of that mechanism. The headline match is partly built into the source prescription, and the authors are transparent enough that a careful reader can see it.\n\nWhat is genuinely new: Cohen et al. reported the IVD observation; this paper adds the full modeling chain—ellipsoid shock reconstruction, MAS corona, EUHFORIA background, PARADISE transport—and shows that a weak flank whose roll-over energy grows with radial distance can produce the delayed high-energy onset and the qualitative IVD. The timing comparison in Table 1 is genuinely good for the D⊥=0 run: most channels within a few minutes of the observed threshold times. That is real evidence the scenario is viable. The shock reconstruction also checks reasonably against in-situ arrival at both PSP and Solar Orbiter, and the paper is honest about the discrepancies in speed and Mach number.\n\nSoft spots, in order of importance. First, the central quantitative result is calibrated rather than predicted. Equation (3) makes pc increase with rs, so any observer connected to an outward-moving flank will see high-energy particles arrive late. The IVD is therefore largely an input of the source model, not an emergent outcome of time-dependent acceleration. The reduction factor η=0.35 is explicitly chosen to match PSP, and the authors admit η likely varies over the shock. They then keep the same η in the D⊥≠0 run, so the cross-field diffusion conclusion rests on a comparison with only one free parameter re-tuned. This does not kill the paper, but it downgrades the claim from \"we explain the IVD\" to \"this mechanism can reproduce the IVD under a tuned cutoff prescription.\" Second, Solar Orbiter is overproduced by an order of magnitude with a harder spectrum. The authors' explanation—overestimated apex shock strength and constant η—is plausible, but it means Eq. (3) is not a reliable quantitative predictor away from the tuned region. Third, the minor role of cross-field diffusion is asserted from timing comparisons without uncertainty propagation; the ±3° connectivity shading is shown in the intensity plots but not folded into the Δt values.\n\nWho this is for: heliospheric SEP modelers and anyone interpreting PSP IS⊙IS onset timing. It deserves serious peer review. A good referee should ask for an explicit statement that the IVD is built into pc(rs), and for a sensitivity run varying η, before the flank-acceleration interpretation is presented as the explanation.","headline":"A credible, well-documented case study whose central IVD reproduction is largely baked into a tuned roll-over momentum prescription, so the flank-acceleration mechanism is plausible but not independently confirmed.","tokens_in":26788,"tokens_out":2126,"would_cite":true,"duration_ms":22710,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The inverse velocity dispersion seen by Parker Solar Probe on 5 September 2022 is attributed to slow acceleration at the expanding shock's weak flank, not to cross-field diffusion.","keywords":["solar energetic particles","coronal mass ejection shocks","inverse velocity dispersion","diffusive shock acceleration","Parker Solar Probe","SEP onset timing","cross-field diffusion","particle transport modeling"],"falsifier":"Measure the proton energy spectrum in the foreshock region upstream of the shock during the roughly 40 minutes before the shock reaches PSP at 17:27 UT on 5 September 2022. The model predicts that the energy at which the spectrum rolls over climbs from well below 1 MeV near onset (16:48 UT) to about 10 MeV at crossing, with higher-energy protons arriving later as the cutoff rises. If the cutoff energy stays roughly constant during that interval, or if a rising cutoff is not accompanied by later arrival of the highest-energy protons, the flank-acceleration explanation is falsified.","tokens_in":25657,"feed_emoji":"☀️","tokens_out":12593,"duration_ms":201350,"temperature":0.7,"pith_summary":"On 5 September 2022, Parker Solar Probe was only 15.4 solar radii from the Sun when a fast coronal mass ejection shock swept past it, and its particle instruments recorded an odd arrival pattern: above about 1 MeV, the higher-energy protons arrived later than the lower-energy ones, the opposite of normal velocity dispersion. The paper sets out to show that this inverse velocity dispersion, together with a roughly 35-minute delay in SEP release relative to Solar Orbiter, came from the shock itself. PSP was magnetically connected to the weak east flank of the expanding shock, where particle acceleration was slow and the maximum accelerated energy rose only gradually as the shock strengthened and moved onto stronger regions. A 3D shock reconstruction combined with a data-driven focused-transport simulation reproduces the observed SEP onset, the evolving energy spectrum, and the inverse velocity dispersion when cross-field diffusion is omitted; adding cross-field diffusion makes the match worse. The authors conclude that the local strength evolution of the connected shock flank governed the delayed release, so release-time interpretations and forecasting need to treat shock acceleration as spatially and temporally inhomogeneous.","feed_headline":"Shock flank, not diffusion, delayed Parker Solar Probe's SEPs","feed_subtitle":"A model reproduces the inverse velocity dispersion above 1 MeV and ties SEP release timing to local shock strength.","key_machinery":"The machinery that carries the argument is the prescription of the proton spectrum injected at each point of the reconstructed shock front: a power law in momentum with an exponential cutoff at a roll-over momentum $p_c$. The cutoff is set by $p_c = \\eta\\, p_{\\mathrm{inj}}\\left(\\frac{\\pi \\epsilon \\sigma}{4\\, r_s/d_i}\\right)^{1/(\\sigma-3)}$, in which $p_{\\mathrm{inj}}$ is the injection momentum, $\\sigma$ the diffusive-shock-acceleration spectral index, $r_s$ the radial distance of the shock point, and $d_i$ the proton inertial length. Because a shock flank that formed only recently has had little time to trap and accelerate particles, its $p_c$ is low; as the shock expands and strengthens, $p_c$ at the PSP-connected point rises from below 1 MeV to about 10 MeV. The constant reduction factor $\\eta = 0.35$ is calibrated against the PSP observations. The time- and space-dependent cutoff, carried to the observer by a focused-transport solver, is what produces the inverse velocity dispersion in the simulated spectrograms. The 3D ellipsoidal shock reconstruction and the coronal and heliospheric MHD background make the prescription specific to this event.","core_discovery":"The central claim is that the SEPs observed at PSP were not released promptly at the Sun and then transported, but were accelerated gradually at the flank of the CME-driven shock to which PSP was magnetically connected. The paper derives release times from velocity dispersion analysis: SEPs reached Solar Orbiter promptly along a Parker-spiral path (release 16:09 ± 4 min), while at PSP the inferred release time is 16:45 ± 1 min with a much shorter path length (about 6.3 $R_\\odot$), indicating release higher in the corona rather than near the solar surface. The shock reconstruction shows PSP's connection point started on a subcritical, oblique flank ($M_{fm} \\approx 1.2$, $\\Theta_{BN} \\approx 32^\\circ$) and only later became supercritical, so the roll-over energy of the accelerated spectrum at that connection point climbed from well below 1 MeV to about 10 MeV by the time the shock crossed the spacecraft. In the simulation this gradual rise of the cutoff energy is what makes protons above ~1 MeV arrive after lower-energy protons. The paper also shows that including cross-field diffusion degrades the timing agreement, and concludes that perpendicular transport had only a minor effect for this near-Sun event.","pith_inferences":["A survey of Parker Solar Probe energetic-particle events should find inverse velocity dispersion preferentially when the spacecraft is connected to a young, weak shock flank; if IVD events occur under other connectivity geometries, the causal link to local flank acceleration would need revision.","The short VDA path length at PSP could be repurposed as a diagnostic of where supercritical conditions first appear along the connected field line, turning a transport-derived quantity into a probe of shock acceleration onset.","The conclusion that cross-field diffusion was minor should not be generalized to all near-Sun events: PSP was inside the steep radial-intensity gradient near the shock, where the direct shock source dominates; farther from the flank, perpendicular transport may matter more.","A position-dependent reduction factor $\\eta$ that decreases at strong shock regions would likely correct the Solar Orbiter overestimate while preserving the PSP match; this is a concrete, testable modification of the model."],"forward_implications":["Release-time estimates from near-Sun spacecraft will be biased late whenever the connected shock flank starts subcritical; the path length derived from velocity dispersion analysis then reflects the height at which the shock first becomes an efficient accelerator rather than a particle path along the Parker spiral.","SEP timing predictions should be made energy-dependent: a weak flank can supply low-energy protons promptly while high-energy protons are delayed by tens of minutes, so forecasts based on a single shock speed or a single release site will miss the observed dispersion pattern.","If cross-field diffusion really was minor for this event, fitting SEP profiles to infer transport coefficients without modeling the spatially varying shock acceleration will overestimate perpendicular diffusion.","The failure of the same injection prescription at Solar Orbiter's strong quasi-parallel connection means the calibration is flank-specific; simulations of other events should adopt position-dependent injection efficiency rather than one global factor."],"supporting_citations":[{"why":"Provides the roll-over momentum formula that fixes the energy cutoff of the injected spectrum, and hence the threshold where the modeled inverse velocity dispersion forms.","marker":"Vainio et al. (2014)"},{"why":"Reports the IS☆IS observations of the event, including the inverse velocity dispersion above ~1 MeV that the simulation is designed to reproduce.","marker":"Cohen et al. (2024)"},{"why":"Supplies the method for estimating 3D shock Mach number and obliquity distributions along the reconstructed wavefront.","marker":"Kouloumvakos et al. (2019)"},{"why":"Presents the focused-transport SEP model used to propagate injected particles from the shock to the observers.","marker":"Wijsen et al. (2019)"},{"why":"Documents the model setup and transport equation used in the simulations.","marker":"Wijsen (2020)"},{"why":"Provides in situ shock parameters and arrival times at PSP and Solar Orbiter used to validate the shock reconstruction and connectivity.","marker":"Trotta et al. (2024)"},{"why":"Earlier argument that evolving shock properties at connected field lines control SEP release timings, which this event is used to test.","marker":"Kouloumvakos et al. (2023)"}],"fun_headline_variants":["Shock flank acceleration, not diffusion, slows SEP arrival","PSP SEP delay points to gradual shock flank acceleration","SEP onset delay explained by slow shock flank acceleration","Inverse velocity dispersion from shock flank, not transport"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on the assumption that a single calibration factor, set to 0.35 and chosen to make the simulated spectra match Parker Solar Probe, correctly describes how the maximum accelerated energy grows at every point of the shock flank that was connected to the spacecraft; if that factor actually varies across the shock surface, or if the steady-state acceleration theory is invalid for this weak oblique flank, the reproduced delay and inverse velocity dispersion could be an artifact of the tuning.","fun_headline_variants_meta":{"raw":{"variants":["Shock flank acceleration, not diffusion, slows SEP arrival","PSP SEP delay points to gradual shock flank acceleration","SEP onset delay explained by slow shock flank acceleration","Inverse velocity dispersion from shock flank, not transport"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1323,"prompt_tokens":1087,"completion_tokens":236,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":703,"completion_tokens_details":{"reasoning_tokens":172}},"tokens_in":703,"tokens_out":236,"duration_ms":3077,"temperature":1.0,"reasoning_tokens":172,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:56:33.296320+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the proton energy spectrum in the foreshock region upstream of the shock during the roughly 40 minutes before the shock reaches PSP at 17:27 UT on 5 September 2022. The model predicts that the energy at which the spectrum rolls over climbs from well below 1 MeV near onset (16:48 UT) to about 10 MeV at crossing, with higher-energy protons arriving later as the cutoff rises. If the cutoff energy stays roughly constant during that interval, or if a rising cutoff is not accompanied by later arrival of the highest-energy protons, the flank-acceleration explanation is falsified.","supporting_citations":[{"cited_title":"2019, , 622, A28, 10.1051/0004-6361/201833958","cited_arxiv_id":null,"evidence_quote":"Presents the focused-transport SEP model used to propagate injected particles from the shock to the observers."}],"review_version":1}