{"id":"30e066ba-2c14-42ac-bd7f-fb7e1dd0ba2a","arxiv_id":"2412.13895","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Full-orbit simulations in a turbulent Parker spiral show drift reduction factors of 0.2-0.9, much larger (less suppression) than predicted by standard drift-suppression theories at low proton energies.","lead":"This paper uses computer simulations of protons flying through a model of the Sun's magnetic field and turbulence to measure how much turbulence reduces the particles' large-scale drift. The result matters for predicting where solar energetic particles travel in the inner solar system, because the simulations find much weaker drift suppression than earlier theories predicted, meaning drift may shape particle arrival patterns more than assumed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Scatter-model baseline likely underestimates the no-turbulence drift by enforcing isotropy, inflating fs; a scatter-free control is needed to test whether the less-suppression claim survives.","rationale":"I read the paper as a serious full-orbit simulation study with a novel heliospheric turbulence model and a clever use of distribution asymmetry to isolate drift from stochastic spreading. The central quantitative claim rests on the definition of the no-turbulence reference. The reader's weakest-assumption analysis identified the scatter baseline; I agree and sharpen it. The paper's defence of the baseline (τΩ >> 1) is a uniform-field TGK argument and does not account for the effect of scattering on the pitch-angle distribution in the inhomogeneous spiral. Adiabatic focusing in a scattering-free spiral would produce a field-aligned beam with a larger drift velocity than the isotropic distribution enforced by the scatter model; hence vd_scat is a low-side estimate of the no-turbulence drift, and fs is a high-side estimate of the true reduction. A simple scatter-free control run can settle this: if the baseline drift rises, the reported fs values shrink and the discrepancy with theoretical predictions may largely disappear. This is a correctable experimental issue rather than a fatal flaw, so the conditional verdict stands. The paper should be accepted only after such a control simulation is reported; the method itself is a useful contribution regardless of the outcome. I am not claiming the authors are wrong, only that the current control simulation does not implement the intended counterfactual, so the headline quantitative claim is not yet established. The proposed test is inexpensive and directly targets the load-bearing assumption.","tokens_in":16656,"tokens_out":9803,"duration_ms":91548,"concrete_test":"Run the same 100,000-particle, 48-hour, first/last 1-au crossing analysis in the pure Parker spiral of Eq. (13) with no turbulence and no scattering (full-orbit Lorentz force only) for simulation sets 1–8, and compare the median vd_no_turb(Δt) with vd_scat(Δt). If vd_no_turb is systematically larger than vd_scat, refit Eq. (18) using vd_no_turb as the baseline; if the resulting fs values approach the Bieber & Matthaeus and Engelbrecht et al. curves, the paper's central claim fails. For set 4, a 1.5-fold increase in baseline would reduce fs from 0.66 to about 0.34, which would substantially narrow the gap with theory.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that turbulence reduces SEP drifts by only a factor fs = 0.2–0.9, far less than the Bieber & Matthaeus (1997) and Engelbrecht et al. (2017) predictions. The reduction factor is defined relative to the 'scatter' simulation (Sec. 3.2, Eq. 18), in which the Parker spiral is kept turbulence-free but particles undergo ad-hoc isotropic pitch-angle scattering with λ∥ from quasi-linear theory. The paper justifies this baseline by noting that τ = λ∥/v >> 1/Ω, so that gyrophase decorrelation is negligible and the scatter model should reproduce the no-turbulence drift. This argument is incomplete. It addresses only the Taylor-Green-Kubo decorrelation mechanism derived for a uniform background field. In the Parker spiral, scattering also competes with adiabatic focusing. In a genuinely turbulence-free spiral, protons injected isotropically near the Sun become focused along the field as B decreases, developing v∥ >> v⊥; the drift velocity in Eq. (B1) scales as (1/2 v⊥² + v∥²), which for a focused beam approaches v², i.e. up to 1.5 times the isotropic value 2v²/3 used in Eq. (B3). The scatter model prevents this focusing by randomizing the velocity direction, keeping the distribution isotropic and thus suppressing the drift relative to the true no-turbulence limit. The measured fs = vd_turb / vd_scat is therefore inflated: the turbulence is actually reducing the drift more than reported. Quantitatively, for set 4 (100 MeV, δB²/B² = 0.6), vd_scat = 120 km/s versus the isotropic 1-au theory of 160 km/s; if the focused no-turbulence drift is about 1.5 times higher (≈240 km/s), fs drops from 0.66 to ≈0.34, close to the Engelbrecht et al. prediction of 0.29. The Appendix B theoretical drift is itself pitch-angle averaged under an isotropy assumption and cannot serve as an independent baseline. Thus the headline 'much less efficient suppression' may be an artifact of the chosen control.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper investigates how magnetic turbulence reduces the large-scale gradient and curvature drifts of solar energetic protons in the inner heliosphere. The authors use full-orbit test particle simulations in a Parker spiral superposed with a new analytic 2D-slab composite turbulence model, and compare the latitudinal drift measured between first and last crossings of the 1 au sphere in these 'turbulence' runs with that in 'scatter' runs where the same Parker spiral is used but turbulence is replaced by ad-hoc isotropic pitch-angle scattering with a quasi-linear mean free path. They introduce a drift reduction factor fs defined as the fitted ratio of the turbulent drift to the scatter drift. The reported fs values range from 0.2 to 0.9 depending on proton energy and turbulence amplitude, and are substantially larger than the predictions of Bieber & Matthaeus (1997) and Engelbrecht et al. (2017). The authors conclude that drifts are much less suppressed by turbulence than previously thought.","tokens_in":17071,"tokens_out":11420,"duration_ms":97745,"significance":"If the quantitative result holds, the paper challenges the strong drift suppression adopted in many cosmic-ray modulation models and provides a new tool for SEP propagation studies. The study is the first to assess drift reduction in a realistic Parker spiral geometry rather than in a uniform or gradient-only background field, and it makes a useful comparison with earlier test-particle simulations. The method of estimating drift from 1 au crossing statistics is novel, and the data are publicly released. The main caveat is that the reference 'no turbulence' drift is itself model-dependent, so the numerical values of fs carry a systematic uncertainty that is not yet quantified.","major_comments":[{"comment":"The reference drift used to define fs is not validated against a true no-turbulence limit. The scatter model randomises the velocity vector at a rate set by λ∥, and for the parameters in Table 1 the scattering time at 1 au is comparable to the adiabatic focusing timescale in the Parker spiral. In a genuinely turbulence-free spiral, adiabatic focusing would make the pitch-angle distribution anisotropic, with (1/2 v⊥² + v∥²) approaching v² rather than the isotropic 2v²/3 used in Eq. (B3). The paper's justification that τ = λ∥/v ≫ 1/Ω only addresses gyrophase decorrelation, not the pitch-angle distribution. Because vd,scat may therefore underestimate the true no-turbulence drift, fs could be systematically overestimated; for set 4, vd,scat = 120 km/s versus the theoretical isotropic value 160 km/s at 1 au, and the focusing contribution to this discrepancy is not quantified. Please run a scatter-free Parker spiral control using the same first-last crossing analysis, or provide an analytic estimate of the focusing enhancement, to support the reported fs values.","section":"Sec. 2.3, Eq. (16)"},{"comment":"The first-last crossing estimator excludes all particles that cross the 1 au sphere fewer than two times, but the paper does not report the fraction of excluded particles or how this selection differs between the turbulence and scatter simulations. At 10 MeV, where the B+ and B- distributions in Fig. 3(a) overlap substantially, this censorship could bias the median drift. Please report the crossing statistics for each simulation set and test the sensitivity of fs to the minimum-crossing criterion.","section":"Sec. 2.3, Eq. (16)"},{"comment":"The fit forces a single constant fs across all ∆t. Figure 2 suggests that the ratio of the turbulence to scatter drift is not obviously constant in ∆t: the turbulence curve appears flatter than the scatter curve, and both decline with ∆t for different reasons (turbulent decoupling vs. spatial sampling of weaker drift regions). The weighting by particle number mitigates the influence of large-∆t bins, but if the true ratio is time-dependent, the fitted fs is a weighted average whose physical meaning is unclear. Please present fs(∆t) or test the constancy of the ratio, for example by fitting in separate ∆t intervals.","section":"Sec. 3.2, Eq. (18)"}],"minor_comments":[{"comment":"The TGK validity argument contains a reversed inequality: from Eq. (10), fs < 0.90 implies τΩ < 3, not τΩ > 3. The conclusion that ∆t > 100 s is much larger than τ is unaffected, but the logic as written is wrong.","section":"Sec. 3.2"},{"comment":"Typo: 'condider' should be 'consider'.","section":"Sec. 2.3"},{"comment":"The caption refers to panel '(d) 1000 MeV' but the figure has only panels (a)-(c); the 1000 MeV panel is (c).","section":"Fig. 3 caption"},{"comment":"Column 4 header has a double bracket 'km s−1]]'.","section":"Table 1"},{"comment":"Averaging the SQLT mean free path over 2 r⊙ to 1 au is a crude approximation; please state the resulting radial variation or justify that the average is representative.","section":"Sec. 2.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the journal scope and the data release is a strength. The main technical risk is the systematic uncertainty in the reference drift; a scatter-free control would substantially strengthen the claims. The self-citation to Laitinen et al. (2023a) is a normal input-model citation, not circularity. The paper is likely publishable after revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read it. First, this is the first full-orbit test of drift reduction in a Parker spiral with a realistic 2D-slab turbulence model, and it gives concrete numbers: for 10–1000 MeV protons, drift survives at 20–90% of the no-turbulence value, much higher than Bieber–Matthaeus and Engelbrecht et al. predict at low energy. Second, the quantity that drives the headline, fs, is measured relative to a “scatter” simulation that injects isotropic pitch-angle scattering into a Parker spiral without turbulence. That baseline is not the same as a genuinely turbulence-free spiral, and the difference may be exactly where the action is.\n\nWhat the paper does well: the simulations are large (100k particles per set, 100 turbulence realisations), the turbulence model is physically motivated (the 2D component is perpendicular to the Parker spiral), the data are on Zenodo, and the authors are explicit about what is and isn't in the model (no heliospheric current sheet, no motional E-field). The new crossing-based drift estimator is a sensible way to extract a systematic latitudinal drift in the presence of strong stochastic spreading, and the comparison with Minnie et al. and Tautz & Shalchi is honest about parameter mismatch. The paper is clearly written and the limitations are discussed.\n\nThe soft spot is the baseline. The paper justifies the scatter model as a stand-in for “no turbulence” by arguing that the parallel scattering time is much longer than the gyroperiod, so gyrophase decorrelation is negligible. But that argument addresses only the TGK decorrelation mechanism. In a Parker spiral, adiabatic focusing also acts: in a scatter-free field, a particle's pitch-angle distribution narrows as B decreases, and the drift velocity depends on (1/2 v_perp^2 + v_par^2), which grows from 2/3 v^2 toward v^2 as focusing develops. The scatter model randomises the pitch angle and keeps the distribution closer to isotropic, so vd_scat is likely an underestimate of the true no-turbulence drift. That would inflate all the fs values. For the 100 MeV moderate-turbulence case, a rough estimate puts the focused drift about 1.5 times higher than the isotropic value, which would pull fs from 0.66 down to roughly 0.3–0.4, near the Engelbrecht prediction. The paper does not provide a scatter-free control run, and the Appendix B theory is itself pitch-angle averaged under isotropy, so it cannot serve as an independent anchor.\n\nThere are smaller issues: excluding particles with zero or one 1-au crossing could select against the most focused particles, and the 10 MeV result carries a sigma of 0.13, so the low-energy statement rests on a broad distribution. These are secondary.\n\nNet: this is a serious paper worth refereeing, but the headline claim of much weaker suppression than theory is not yet nailed down. The baseline choice needs scrutiny, ideally a scatter-free control or at least a quantitative check of the pitch-angle anisotropy in the scatter runs. If the baseline issue is resolved, the fs values could move substantially. I'd send it to review and ask for that analysis before publication. I'd cite the method and turbulence model, but I'd be cautious about quoting the fs values without caveats.","headline":"First full-orbit test of drift reduction in a Parker spiral with 2D-slab turbulence, but the headline numbers rest on a scatter-model baseline that may understate the true no-turbulence drift and inflate fs.","tokens_in":17727,"tokens_out":6584,"would_cite":true,"duration_ms":63517,"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":"Turbulence in the heliosphere cuts solar-energetic-particle drifts by only 20–90%, far less than theory predicts, so drifts remain a major force in SEP and cosmic-ray transport.","keywords":["solar energetic particles","guiding centre drift","heliospheric turbulence","Parker spiral","test particle simulations","drift reduction factor","cosmic ray transport","pitch-angle scattering"],"falsifier":"Run the same particle energies and turbulence parameters in a full-orbit simulation with a constant background field plus a transverse gradient, following the earlier gradient-field approach: if those simulations give drift reduction factors matching the theoretical curves rather than the paper's $f_s$ values, the Parker-spiral geometry would be the origin of the discrepancy. Observationally, multi-spacecraft measurements of the heliolatitude dependence of ~100 MeV SEP events would test whether latitudinal drifts survive at the level the simulations imply.","tokens_in":16479,"feed_emoji":"☀️","tokens_out":8606,"duration_ms":69612,"temperature":0.7,"pith_summary":"The paper asks how strongly the plasma turbulence of the solar wind suppresses the large-scale guiding-centre drift that the Parker-spiral magnetic field imposes on energetic particles. The authors run full-orbit test-particle simulations of 10, 100 and 1000 MeV protons in a newly developed analytic heliospheric turbulence model, and measure the drift reduction factor $f_s$ against a no-turbulence scattering baseline. They find $f_s$ between 0.2 and 0.9, depending on energy and turbulence amplitude, which is a much weaker suppression than the decorrelation-based theoretical models of 1997 and 2017 predict, especially at low proton energies. If this result holds, guiding-centre drifts are a significant factor in the evolution of solar energetic particle intensities and in cosmic-ray transport in the inner heliosphere.","feed_headline":"Turbulence suppresses solar particle drifts less than theory says","feed_subtitle":"Full-orbit runs put drift reduction at 0.2–0.9, keeping drifts significant for SEP and cosmic-ray transport.","key_machinery":"The central quantity is the drift reduction factor $f_s = \\Omega^2\\tau^2/(1+\\Omega^2\\tau^2)$ from the Taylor-Green-Kubo formalism, where $\\Omega$ is the particle gyrofrequency and $\\tau$ the gyromotion decorrelation timescale; theoretical models estimate $\\tau$ from field-line or particle cross-field diffusion, while the paper extracts $f_s$ directly by comparing two sets of simulations. The turbulence model adds a 2D-slab composite fluctuation field to the Parker spiral, with the dominant 2D component's wave vector and magnetic field vector both normal to the spiral, and the drift is measured by a new method: the median change in colatitude between a particle's first and last crossings of the 1 au sphere, averaged over 100 turbulence realisations.","core_discovery":"The paper reports that in a turbulent Parker-spiral heliosphere, the mean latitudinal drift velocity of energetic protons is reduced to a factor $f_s = 0.2$–$0.9$ of the drift expected in the same field without turbulence, with the least reduction at high energy and low turbulence amplitude: $f_s \\approx 0.9$ for 1000 MeV protons and $\\delta B^2/B^2 = 0.2$, down to $f_s \\approx 0.2$ for 10 MeV protons and $\\delta B^2/B^2 = 0.6$. The suppression is therefore real, but considerably weaker than the values below 0.1 that decorrelation-based theoretical models predict at low energies. The authors conclude that drifts should be retained in models of solar energetic particle propagation and cosmic-ray modulation, at least for protons above about 100 MeV and for heavier ions with larger Larmor radii.","pith_inferences":["If the modest suppression holds for galactic cosmic rays as well, modulation models that currently suppress drift coefficients to match observations may be attributing too much of the suppression to turbulence; the heliospheric current sheet and other large-scale structure may be doing more of the work.","The crossing-based drift-measurement method could be applied to spacecraft data by tracking the centroid of SEP intensity in latitude over successive solar rotations, offering a direct observational check of the simulated reduction factors.","A natural extension is to repeat the analysis with a pre-computed 3D turbulence grid rather than the analytic Fourier-mode model; a large change in $f_s$ would show the analytic model's strict transverse-2D geometry influences the result.","The present model excludes the heliospheric current sheet and the motional electric field, so the quoted $f_s$ values apply to unipolar-field regions; particles crossing the sector boundary may experience different effective drift reduction."],"forward_implications":["Models of solar energetic particle propagation should include guiding-centre drifts; at 100 MeV and above, drift reduction is only about 10–40%, so omitting drifts misplaces particles in heliolatitude and longitude.","The drift reduction factor for 100 MeV protons spans 0.91 at weak turbulence to 0.23 at strong turbulence, meaning solar-cycle variations in turbulence amplitude change how much drift matters.","Since drift speed scales with Larmor radius, heavier ions at the same energy per nucleon experience even less relative suppression than protons, strengthening drift effects in heavy-ion SEP events.","The theoretical drift-reduction models overpredict suppression at low energies, which indicates their decorrelation timescale is too short or is set by a different physical process than the one the models assume."],"supporting_citations":[{"why":"Supplies the theoretical drift reduction model based on a decorrelation timescale, the main comparison target of the paper.","marker":"Bieber & Matthaeus (1997)"},{"why":"Provides the alternative theoretical drift reduction estimate using perpendicular mean free path, also compared against.","marker":"Engelbrecht et al. (2017)"},{"why":"Supplies the analytic heliospheric turbulence model used in the turbulent simulations, including the 2D component transverse to the Parker spiral.","marker":"Laitinen et al. (2023a)"},{"why":"Supplies the scatter simulation approach with ad-hoc pitch-angle scattering used as the no-turbulence baseline.","marker":"Marsh et al. (2013)"},{"why":"Gives the theoretical drift velocity in the Parker spiral used for comparison and for defining the drift baseline.","marker":"Dalla et al. (2013)"},{"why":"Provides the Fourier-mode method used to realise the turbulent magnetic field in the simulations.","marker":"Giacalone & Jokipii (1999)"},{"why":"Provides the prior full-orbit simulation in a gradient field whose drift reduction results are compared with the paper's.","marker":"Minnie et al. (2007)"},{"why":"Provides prior simulation-based drift reduction values in constant background field for comparison.","marker":"Tautz & Shalchi (2012)"}],"fun_headline_variants":["Solar particle drifts resist turbulence more than expected","Drift suppression weaker than theory predicts for solar particles","SEP drifts survive turbulence better than models predict","Drift reduction only 0.2–0.9 in turbulent heliosphere"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The drift reduction factor is measured relative to the scatter simulation, which replaces turbulence with ad-hoc isotropic pitch-angle scattering; if that baseline does not faithfully represent the turbulent run's parallel transport or the no-turbulence drift, the quoted values of $f_s$ would be biased.","fun_headline_variants_meta":{"raw":{"variants":["Solar particle drifts resist turbulence more than expected","Drift suppression weaker than theory predicts for solar particles","SEP drifts survive turbulence better than models predict","Drift reduction only 0.2–0.9 in turbulent heliosphere"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00072,"raw_usage":{"total_tokens":3265,"prompt_tokens":1008,"completion_tokens":2257,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":624,"completion_tokens_details":{"reasoning_tokens":2197}},"tokens_in":624,"tokens_out":2257,"duration_ms":15149,"temperature":1.0,"reasoning_tokens":2197,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:40:44.591866+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same particle energies and turbulence parameters in a full-orbit simulation with a constant background field plus a transverse gradient, following the earlier gradient-field approach: if those simulations give drift reduction factors matching the theoretical curves rather than the paper's $f_s$ values, the Parker-spiral geometry would be the origin of the discrepancy. Observationally, multi-spacecraft measurements of the heliolatitude dependence of ~100 MeV SEP events would test whether latitudinal drifts survive at the level the simulations imply.","supporting_citations":[{"cited_title":"E., Strauss , R","cited_arxiv_id":null,"evidence_quote":"Provides the alternative theoretical drift reduction estimate using perpendicular mean free path, also compared against."},{"cited_title":"S., Dalla , S., Kelly , J., & Laitinen , T","cited_arxiv_id":null,"evidence_quote":"Supplies the scatter simulation approach with ad-hoc pitch-angle scattering used as the no-turbulence baseline."},{"cited_title":"S., Kelly , J., & Laitinen , T","cited_arxiv_id":null,"evidence_quote":"Gives the theoretical drift velocity in the Parker spiral used for comparison and for defining the drift baseline."},{"cited_title":"W., Matthaeus , W","cited_arxiv_id":null,"evidence_quote":"Provides the prior full-orbit simulation in a gradient field whose drift reduction results are compared with the paper's."},{"cited_title":"C., & Shalchi , A","cited_arxiv_id":null,"evidence_quote":"Provides prior simulation-based drift reduction values in constant background field for comparison."}],"review_version":1}