{"id":"f27bf4d4-738a-4819-992d-47da3b7b4e0b","arxiv_id":"2507.16990","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Two bidirectional solar energetic particle events show sunward and anti-sunward proton streams with different release times and path lengths, attributed to shock and flare acceleration inside small-scale flux ropes.","lead":"Solar Orbiter observed two rare solar energetic particle events in which protons streamed both toward and away from the Sun at the same time, with two clear velocity dispersion signatures. The authors argue that both events occurred inside small-scale magnetic flux ropes, that one proton population came from a flare and the other from a CME shock, and that the measured particle paths can help map the shape of the flux rope.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Event 2's 4.62 au S2 path length may be inflated by energy-dependent S1 background contamination, weakening the flux-rope topology claim.","rationale":"The central claim rests on two pillars: (i) two SEP populations with opposite anisotropies originate from flare acceleration and CME-shock acceleration, and (ii) the path-length asymmetry within a flux rope constrains the rope's topology. Pillar (ii) is the less secure. For event 1, the SFR is identified by PyGS, the S1 release time matches the type III burst, and EAS shows bidirectional electrons, so the first event carries independent weight. For event 2, the SFR is only inferred from wavelet helicity and the anti-sunward S2 onset is visually estimated on top of the S1 background; the paper itself flags the S2 path length as less robust. The specific failure mode worth testing is energy-dependent onset bias: S1's decaying, velocity-dispersed background contaminates low-energy S2 channels more than high-energy channels, which systematically delays the low-energy onset estimates and steepens the VDA slope. That would produce a spuriously large path length. If so, the 4.62 au value, which is then used to argue for a long flux-rope leg, cannot support the topology claim. A controlled injection-and-recovery test on the actual count rates settles this concern. I agree with the reader's weakest-assumption identification; the proposed test is the concrete way to see whether the concern lands. No independent support, such as machine-checked derivations or released analysis code for the VDA, is provided, so the path length estimates remain the least externally validated part of the analysis.","tokens_in":16778,"tokens_out":6120,"duration_ms":72285,"concrete_test":"Use the event 2 data to quantify the VDA bias. Take the measured S1 velocity-dispersed background and inject a synthetic S2 population with a known path length near the nominal Parker value of about 1.5 au and a release time about 30 minutes after S1; add Poisson noise at the observed count levels; run the same visual/averaged-onset procedure or an automated CUSUM on background-subtracted intensities; and compare the recovered path length. If the recovered S2 path length is biased toward the reported 4.62 au, the paper's topology constraint must be revised. If it recovers the injected value, the concern is resolved. A cheaper cross-check is to recompute the S2 VDA after subtracting an extrapolated S1 background model from each energy channel and see whether the inferred path length changes by more than about 1 au.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the VDA for event 2's S2. The paper itself concedes in Section 3 that the anti-sunward onset was visually estimated and that under strong scattering the path length might be less robust, especially for S2. This matters because the 4.62 au path length is the main quantitative support for the flux-rope-topology interpretation in event 2, where the SFR is only inferred from wavelet helicity and not confirmed by PyGS. The likely bias is not random: the S2 anti-sunward onset is measured on top of the decaying, velocity-dispersed S1 injection. Low-energy S2 channels are the ones most contaminated by S1 background, so their apparent onsets can be systematically delayed relative to high-energy channels. In a VDA fit of onset time versus 1/v, an energy-dependent late bias inflates the slope, and hence the inferred path length, while also shifting the release-time intercept. If this is what produces 4.62 au, the asymmetry between S1 and S2 path lengths, and therefore the inferred long field-line path inside the flux rope, is an artifact of background contamination rather than a topological constraint.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports two solar energetic particle (SEP) events observed by Solar Orbiter on July 13, 2023, and July 23, 2024, in which the proton anisotropy reveals two oppositely directed streaming populations during the onset phase. For each event the authors identify two separate velocity-dispersion tracks, perform velocity dispersion analysis (VDA) to derive release times and path lengths, and combine the path lengths with in situ magnetic-field diagnostics (automated flux-rope identification via PyGS and wavelet magnetic helicity) to argue that a small-scale flux rope channels flare-accelerated anti-sunward particles along a shorter leg and CME-shock-accelerated sunward particles along a longer leg. The paper concludes that the sunward population is delayed, has a harder spectrum, and is more intense, and that the inferred path-length asymmetry constrains the flux-rope topology.","tokens_in":16868,"tokens_out":10306,"duration_ms":107330,"significance":"Event 1 is a convincing and rare observation: two clean velocity-dispersion tracks with opposite anisotropies, a type III radio burst coincident with the prompt release, a PyGS-identified flux rope, and bidirectional suprathermal electrons. If the interpretation is correct, the work provides a direct observational link between small-scale flux ropes and bidirectional SEP transport, and it proposes a new way to use simultaneous injections to constrain flux-rope topology. The paper is generally careful in stating its assumptions and in acknowledging the weaker constraints for event 2; the use of the public PyGS package and the explicit caveat about the event-2 path length are positive features. However, the strength of the event-2 evidence is currently overstated in the abstract, and the quantitative path-length asymmetry in event 2 rests on a visually estimated onset under strong scattering, so the central claim is only partly supported by the second event.","major_comments":[{"comment":"The 4.62 au S2 path length is the main quantitative support for the path-length asymmetry in event 2 and enters the flux-rope topology argument, yet it is derived from a visually estimated anti-sunward onset in a period of strong scattering and high background. The manuscript itself states that \"the path length for this event might be less robust, especially for the S2.\" Because the S2 onset is measured on top of a decaying, velocity-dispersed S1 population, an energy-dependent late bias in the low-energy onset times would steepen the VDA slope and inflate the inferred path length; the present analysis offers no test of this bias. Please provide a quantitative robustness check, for example a background-subtracted onset fit, a fit restricted to the least contaminated high-energy channels, or a Monte Carlo over plausible onset choices; if such a check is not possible, state explicitly in the abstract and conclusions that the event-2 path-length asymmetry is only qualitative and does not by itself constrain flux-rope topology.","section":"§3, Event 2 (Fig. 3 right panels) and Appendix B, Eq. (B.1)"},{"comment":"The abstract says that \"in both cases, small-scale flux ropes were identified in situ,\" but for event 2 the PyGS algorithm did not reliably detect a flux rope and the identification rests on wavelet helicity with thresholds |σm|≥0.6 and σr≤−0.5, described in the text as a \"potentially Alfvénic SFR.\" This overstates the certainty for event 2. The abstract and the conclusions should carry the same qualification as the body of the paper, namely that the event-2 flux rope is inferred rather than identified.","section":"Abstract and §3 (Event 2)"},{"comment":"The spectral arguments are made without uncertainties: the text contrasts a −2.45 spectral index for S1 with −1.24 for S2 and states that the latter is \"close to the limit of theoretical spectral index −1,\" but no fit uncertainties, energy ranges, or goodness-of-fit values are quoted. Since the spectral hardness is one of the three main discriminants between shock and flare origin, please report the fit errors and the fitting range for each spectrum, or downgrade the spectral claim to a qualitative statement.","section":"§3 (Event 2) and Fig. 4"},{"comment":"The 30-minute release delay is interpreted as shock-acceleration time using assumed values of 1.43 R⊙ for the shock formation height, 3 R⊙ for the release height, and 750 km/s for the average shock speed. These parameters are not measured for these events, and the resulting 24-minute estimate is therefore only a plausibility argument; the delay could also reflect energy-dependent release or transport effects. Please frame this paragraph explicitly as a consistency check and state how the conclusion would change if the assumed heights and speeds were varied over the ranges quoted in the cited statistical studies.","section":"§3 (release-time interpretation)"}],"minor_comments":[{"comment":"The text refers to \"solid lines\" and \"purple lines\" in Fig. 3, while the figure caption describes \"dashed black and magenta lines\"; please make the line-style and color terminology consistent between text and figure.","section":"§3 and Fig. 3"},{"comment":"The notation L0,sun and L0,asun is confusing: L0,sun is derived from the 1.75 au path of the anti-sunward streaming population, whereas the phrase \"sunward streaming\" elsewhere refers to the 3.04 au population. Please rename these quantities so that the mapping between telescopes, streaming directions, and path lengths is explicit.","section":"Appendix B"},{"comment":"The comparison of derived path lengths with the nominal Parker spiral should quote the actual solar wind speed used, and its uncertainty, for each event; Fig. 1 assumes 400 km/s while Appendix B quotes 300 km/s, and the conclusion that the paths are \"substantially greater\" depends on this choice.","section":"§3 and Fig. 1"},{"comment":"Please clarify the radio-burst timing convention: the caption says the vertical line is shifted 7.9 and 7.2 minutes earlier, while the text says the burst began at approximately 06:41 UT from the Sun; state explicitly whether the quoted onset is the time at the Sun or the observed time at the spacecraft.","section":"§3 and Fig. 3 caption"},{"comment":"There are two occurrences of \"consisting with\" where \"consistent with\" is meant; please correct these in the paragraph discussing the acceleration-time estimate.","section":"§3"}],"recommendation":"major_revision","confidential_remarks":"This is a two-event case study in which event 1 is much stronger than event 2. The main risk is overclaiming in the abstract and conclusions: the event-2 flux rope is only inferred, and the 4.62 au path length is not robust under strong scattering. If the authors can either quantify the background-bias effect on the event-2 S2 onset or explicitly downgrade the event-2 path-length and topology claims, the paper would be publishable as a case report. I do not see a circularity problem in the VDA or the flux-rope argument."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about arXiv:2507.16990.\n\nFirst, the observation is genuinely new: two SEP events with two distinct velocity dispersion tracks and opposite anisotropies, both associated with small-scale flux ropes. I don't recall a previous report of that combination, and event 1 is fairly clean—clear sunward and anti-sunward telescope signatures, a type III burst coincident with the first release, and a PyGS-identified flux rope with bidirectional electrons. The interpretation—prompt flare-accelerated anti-sunward protons and delayed shock-accelerated sunward protons—is plausible, and the authors argue sensibly against the mirror-effect alternative.\n\nSecond, the weak spot is event 2's sunward-streaming track (S2), specifically the 4.62 au path length. That track comes from a visually estimated onset in the anti-sunward telescope, on top of a decaying S1 background. The stress-test concern that energy-dependent background contamination inflates the VDA slope is plausible, and the paper's one-line caveat that the S2 path length is less robust does not quantify the bias. This is a legitimate referee question. It does not sink the paper, though, because the quantitative flux-rope geometry in Appendix B uses only event 1, which already shows a substantial path-length asymmetry (1.75 vs 3.04 au) and is much cleaner. Event 2 is supporting evidence, not the foundation.\n\nOther soft spots are minor: spectral indices have no uncertainties, and the shock-acceleration timing argument uses representative shock parameters rather than measured ones. These are consistency checks, and the authors are transparent about the assumptions.\n\nThis is a two-event case study with public data and established tools, so its reach is limited to heliospheric SEP physics and solar wind structure. Within that subfield it is a useful, honest contribution. I would send it to a serious referee: the main pressure points are a sensitivity analysis for event 2's onset estimates and error bars on the spectra, but the paper deserves referee time.","headline":"A genuinely new two-event SEP observation with a plausible flux-rope interpretation; event 2's long path length is the weak link but not fatal.","tokens_in":17580,"tokens_out":8821,"would_cite":true,"duration_ms":81731,"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":"Bidirectional particle bursts trace a flux rope's long field lines.","keywords":["solar energetic particles","bidirectional anisotropy","velocity dispersion analysis","small-scale magnetic flux ropes","magnetic helicity","Solar Orbiter","CME-driven shock acceleration","solar flares"],"falsifier":"Re-measure event 2's anti-sunward onset from background-subtracted spectra with an automated change-point detector; if a shift of roughly ten minutes or more removes the sunward/anti-sunward path-length asymmetry or brings the lengths back to the nominal Parker value, the flux-rope-guided long-path interpretation for that event fails. A future bidirectional event inside a flux rope with near-equal release times, intensities, and spectral indices for both streams would instead support the mirror-reflection alternative.","tokens_in":1797,"feed_emoji":"🧲","tokens_out":2163,"duration_ms":103698,"temperature":0.7,"pith_summary":"Solar energetic particle events usually stream outward along open magnetic field lines, but this paper reports two rare events, observed by Solar Orbiter, in which protons streamed simultaneously toward and away from the Sun. In both events the two streams arrived with separate velocity-dispersion signatures, so their release times and path lengths could be derived independently. The paper argues that the promptly released anti-sunward stream came from flare acceleration, while the delayed, harder-spectrum, more intense sunward stream came from CME-driven shock acceleration, and that a small-scale magnetic flux rope encountered in situ guided the two populations along different, unusually long field lines. If this is right, bidirectional events become a direct way to measure the topology and central axis length of small flux ropes in the inner heliosphere.","feed_headline":"Bidirectional particle bursts trace a flux rope's long field lines","feed_subtitle":"Separate sunward and anti-sunward tracks reveal the rope's long, twisted field lines.","key_machinery":"The argument runs on two diagnostics in parallel. Velocity dispersion analysis fits onset time as a function of inverse particle speed, giving release time and path length separately for the sunward and anti-sunward populations, while the first-order anisotropy tells which population dominates at each energy and time. On the magnetic side, an automated Grad-Shafranov-based flux-rope identification routine and a wavelet-derived normalized reduced magnetic helicity locate the small-scale flux rope, with cross-helicity and residual-energy checks distinguishing a twisted rope from Alfvenic fluctuations. The link between the two is the path-length excess: the fitted lengths exceed the nominal Parker spiral length, and that excess is assigned to helical winding around the flux rope's axis, allowing an estimate of the rope's central axis length.","core_discovery":"The paper reports two solar energetic particle events in which protons arriving at Solar Orbiter showed two distinct velocity-dispersion tracks with opposite first-order anisotropies at the same time. In each event, the promptly arriving anti-sunward stream matched the timing of a type III radio burst, so it is assigned to flare acceleration; the delayed, more intense, spectrally harder sunward stream is assigned to a CME-driven shock that takes roughly half an hour to become an efficient accelerator. Small-scale magnetic flux ropes were identified in situ during both events, and the velocity dispersion analysis implies path lengths of 1.75 au and 3.04 au in event 1 and 1.51 au and 4.62 au in event 2, all substantially longer than the nominal Parker spiral length. The paper concludes that the flux rope's closed, twisted field lines guide the two populations along different, longer paths, and that simultaneous velocity dispersion from two distinct sources can directly constrain the flux rope's topology and central axis length.","pith_inferences":["Beyond the paper, if small flux ropes are as frequent in the slow solar wind as the cited surveys indicate, some solar energetic particle events with anomalously long path lengths or unusual anisotropies may be connectivity effects of these ropes rather than scattering effects.","A testable extension is to search archival events for two velocity-dispersion tracks and check whether the path-length ratio of sunward to anti-sunward streams correlates with flux-rope twist, as the helical winding picture would predict.","The spectral asymmetry proposed here implies a statistical prediction: sunward shock-accelerated streams should cluster near the strong-compression spectral index near minus one, while flare-related anti-sunward streams should be softer, which could be checked across a larger sample.","A rival explanation, mirror reflection, would predict similar release times, intensities, and spectral indices for the two streams; future multipoint observations inside a well-resolved flux rope could test this cleanly by comparing fluxes at the two feet."],"forward_implications":["If the interpretation holds, a single spacecraft can estimate the central axis length of a small flux rope by combining velocity-dispersion path lengths with the observed helical twist, as done in Appendix B for event 1.","The near-zero prompt release of the anti-sunward stream and the roughly 30-minute delayed release of the sunward stream become a usable discriminator between flare and CME-shock acceleration inside a closed magnetic structure.","Small-scale flux ropes, not only large magnetic clouds, must be included in solar energetic particle transport models to account for sunward streams and path lengths longer than nominal Parker spirals.","The rarity of clean bidirectional events is partly a detection effect: the first injection must be much weaker than the second, or its anisotropy signature masks the second source."],"supporting_citations":[{"why":"Supplies the Solar Orbiter EPD ion measurements, together with Wimmer-Schweingruber et al. 2021, from which the anisotropy and velocity-dispersion tracks are derived.","marker":"Rodríguez-Pacheco et al. 2020"},{"why":"Documents the EPD instrument response and calibration used to interpret the directional proton fluxes.","marker":"Wimmer-Schweingruber et al. 2021"},{"why":"Gives the Poisson-CUSUM change-point method used to set the onset times that feed the velocity dispersion analysis.","marker":"Lucas 1985"},{"why":"Provides the automated Grad-Shafranov flux-rope identification algorithm used to detect the event-1 small-scale flux rope.","marker":"Hu et al. 2018"},{"why":"Supplies the Grad-Shafranov reconstruction method used to map the flux rope cross-section and infer its twist.","marker":"Hu & Sonnerup 2002"},{"why":"Defines the normalized reduced magnetic helicity wavelet diagnostic used to identify the event-2 helical structure.","marker":"Zhao et al. 2020"},{"why":"Describes the mirror-reflection alternative interpretation that the paper must exclude to attribute the sunward stream to a second source.","marker":"Li et al. 2020"},{"why":"Provides the proton release-height statistics used to argue that a delayed sunward release is compatible with CME-driven shock acceleration.","marker":"Kouloumvakos et al. 2015"}],"fun_headline_variants":["Bidirectional SEPs unveil flux rope's longer particle paths","Flux rope's field lines channel dual solar particle streams","Rare two-way particle bursts expose flux rope topology","Sunward and anti-sunward protons trace flux rope's shape","Two acceleration sources send particles along flux rope paths"],"cache_read_input_tokens":19584,"weakest_assumption_plain":"The load-bearing premise is that the anti-sunward onset time in the second event, which had to be estimated by eye against a bright, scattered background, is accurate enough for velocity dispersion analysis to yield the 4.62 au path length that underlies the flux-rope argument.","fun_headline_variants_meta":{"raw":{"variants":["Bidirectional SEPs unveil flux rope's longer particle paths","Flux rope's field lines channel dual solar particle streams","Rare two-way particle bursts expose flux rope topology","Sunward and anti-sunward protons trace flux rope's shape","Two acceleration sources send particles along flux rope paths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000503,"raw_usage":{"total_tokens":2498,"prompt_tokens":1026,"completion_tokens":1472,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":642,"completion_tokens_details":{"reasoning_tokens":1391}},"tokens_in":642,"tokens_out":1472,"duration_ms":12385,"temperature":1.0,"reasoning_tokens":1391,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:59:34.508492+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure event 2's anti-sunward onset from background-subtracted spectra with an automated change-point detector; if a shift of roughly ten minutes or more removes the sunward/anti-sunward path-length asymmetry or brings the lengths back to the nominal Parker value, the flux-rope-guided long-path interpretation for that event fails. A future bidirectional event inside a flux rope with near-equal release times, intensities, and spectral indices for both streams would instead support the mirror-reflection alternative.","supporting_citations":[{"cited_title":"F., Mason , G","cited_arxiv_id":null,"evidence_quote":"Supplies the Solar Orbiter EPD ion measurements, together with Wimmer-Schweingruber et al. 2021, from which the anisotropy and velocity-dispersion tracks are derived."},{"cited_title":"F., Janitzek , N","cited_arxiv_id":null,"evidence_quote":"Documents the EPD instrument response and calibration used to interpret the directional proton fluxes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Poisson-CUSUM change-point method used to set the onset times that feed the velocity dispersion analysis."},{"cited_title":"L., Zank , G","cited_arxiv_id":null,"evidence_quote":"Defines the normalized reduced magnetic helicity wavelet diagnostic used to identify the event-2 helical structure."},{"cited_title":"2015, , 580, A80","cited_arxiv_id":null,"evidence_quote":"Provides the proton release-height statistics used to argue that a delayed sunward release is compatible with CME-driven shock acceleration."}],"review_version":1}