{"id":"1a37b345-2f51-4371-b09e-92a2e743aab9","arxiv_id":"2508.13835","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In 2.5D MHD simulations of breakout CMEs, increasing background poloidal field strength suppresses eruptions, and the growth rate of absolute net current helicity best discriminates failed, single, and multiple eruptions.","lead":"This paper runs 2.5D magnetohydrodynamic simulations of solar coronal mass ejections and finds that a stronger background magnetic field suppresses eruptions, while the rate at which magnetic helicity builds up separates failed, single, and multiple eruptions. A smart generalist should read it because that growth rate could become a practical forecasting indicator for space weather.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ANCH-slope criterion is collinear with the imposed shear amplitude, so the simulations do not establish ANCH growth rate as the crucial determinant of eruption likelihood.","rationale":"The reader's conditional verdict is appropriate. The most load-bearing problem is not the numerical threshold shift alone but that the central claim attributes causation to ANCH growth rate when the experimental design cannot separate that quantity from the imposed shear amplitude. This is a logical gap rather than a simple precision issue, and it applies even if every run is perfectly converged. I therefore would keep the paper at CONDITIONAL: the qualitative trend that stronger background poloidal fields suppress breakout CMEs is clean and physically plausible, and the ANCH-slope diagnostic is worth testing further, but the headline claim that ANCH growth rate is the crucial factor is not established by three single runs with one varied input. The proposed constant-v0 ensemble would test whether ANCH slope carries predictive power after removing collinearity with v0, and a resolution repeat would establish whether the 0.086-unit slope separation is numerically meaningful. I do not see grounds for rejection: the simulations use a credible code, the setup is described in enough detail to reproduce, and the authors explicitly acknowledge the resolution dependence, which is a limitation statement that should count in their favor. The missing piece is a demonstration that the diagnostic is not merely a proxy for shear amplitude.","tokens_in":15021,"tokens_out":6848,"duration_ms":79256,"concrete_test":"Generate a small constant-v0 ensemble: fix the maximum shear velocity at 36.2 km/s and vary the background poloidal field in 0.1 G steps around 2.2 G, covering both eruptive and failed outcomes. If the ANCH slope separates the eruptive from the failed cases within this fixed-v0 set, the criterion is not merely a proxy for shear amplitude; if the slopes overlap across outcomes, the central claim must be weakened to a parameterization of shear-driven eruptions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that ANCH growth rate be a controlling variable, not just a correlate. In §3.2 the three eruption scenarios are generated by changing only the maximum shear velocity v0 in Eq. (8): 35.8, 36.2, and 38.8 km/s. The reported ANCH slopes (0.412, 0.498, 0.645 simulation units) therefore increase monotonically with the prescribed helicity-injection rate. Any monotone diagnostic of injected helicity, such as magnetic energy slope or TUCH slope, would order the three cases the same way; the curves in Fig. 9 are effectively a relabeled version of the input parameter. The abstract's statement that ANCH growth rate 'determines the likelihood' of eruption is a causal claim, but the design has a single varied parameter and one run per value, so causation and parameter dependence cannot be separated. The numerical-resistivity caveat in §2 compounds this: the regime boundaries sit 0.4-2.6 km/s apart, and the authors state that the required v0 changes with resolution, so even the quantitative slope thresholds are not shown to be converged. In addition, the practical forecast link is indirect: SHARP ANCH is a photospheric surface quantity, whereas the simulated ANCH is a coronal volume integral over B·J (Eqs. 9-11); no mapping between the two is supplied. The result is a plausible diagnostic hypothesis, not an established physical criterion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents 2.5D MHD simulations of breakout CMEs using MPI-AMRVAC. The setup combines a background dipole field with a triple-arcade quadrupolar configuration; after relaxation, a time-dependent shear flow is imposed at the base of the central arcade. Varying the maximum shear velocity (35.8, 36.2, and 38.8 km/s) yields three outcomes: failed eruption, single eruption, and multiple eruptions. The authors also vary the polar background field strength from 2.2 G to 2.5 G and report that stronger poloidal fields suppress eruptions. Tracking global magnetic parameters (excess magnetic energy, total unsigned current helicity, and absolute net current helicity), they conclude that the growth rate of absolute net current helicity (ANCH) is the crucial factor determining eruption likelihood, and they connect this to the excess of weak CMEs in Solar Cycle 24.","tokens_in":15383,"tokens_out":6040,"duration_ms":65204,"significance":"If established, the ANCH-slope criterion would be a practically valuable forecasting precursor, since ANCH-like quantities are available from SHARP photospheric magnetograms. The poloidal-field result is also physically interesting and offers a plausible mechanism for the Solar Cycle 24 weak-CME excess through reduced background field strength. The numerical setup follows established practice (breakout model, AMR, GLM divergence control), and the supplementary movies are a useful resource. However, the central causal claim is not yet supported: the evidence rests on three single runs in which the proposed predictor is controlled by the same input parameter (shear amplitude), and the eruption thresholds are admitted to be resolution dependent. The paper is a worthwhile exploratory study, but the abstract overstates what the simulations demonstrate.","major_comments":[{"comment":"The ANCH-slope criterion is collinear with the input shear amplitude. The three scenarios are generated by changing only v0 in Eq. (8) (35.8, 36.2, and 38.8 km/s), and the reported ANCH slopes (0.412, 0.498, and 0.645 in simulation units) increase monotonically with v0. Any monotone diagnostic of injected helicity, such as the slope of magnetic energy or TUCH, would order the three cases in the same way, so the design does not separate the proposed predictor from the control parameter. To support the causal claim, the authors should vary other physical parameters (e.g., background field, heating, or arcade geometry) at fixed v0, or include runs with overlapping ANCH slopes and different outcomes, and compare the discriminative power of ANCH slope against energy slope and TUCH slope in a quantitative analysis.","section":"§3.2, Fig. 9, Eq. (8)"},{"comment":"The eruption thresholds are resolution dependent. The manuscript states that increasing resolution lowers the maximum shearing velocity needed to initiate an eruption, yet the three regime boundaries are separated by only 0.4-2.6 km/s. No convergence study is provided, so the quoted thresholds of 35.8, 36.2, and 38.8 km/s, and hence the associated ANCH slope thresholds, may shift with grid resolution. Without a resolution study, the quantitative slope values in Fig. 9 cannot be distinguished from numerical artifacts. This is a load-bearing issue because the central claim depends on these threshold separations.","section":"Section 2, final paragraph; §3.2"},{"comment":"The forecast link to observations is not established. The simulated ANCH is a coronal volume integral of B·J (Eq. 11), while SHARP ANCH is a photospheric surface quantity. No mapping, unit conversion, or comparison of time evolution between the two is supplied, and the reported slopes are given only in unspecified simulation units. The paper should either provide a forward model connecting the simulated volume quantity to the observable photospheric quantity, or explicitly limit the claim to a qualitative diagnostic hypothesis rather than a ready-to-use forecasting metric.","section":"Section 4, Eqs. (9)-(11)"},{"comment":"The poloidal-field result is presented mainly through field-line snapshots and qualitative statements. Quantitative eruption metrics for the four Bpol values, such as flux-rope height, velocity, magnetic energy, or helicity budgets, are not reported, and each case is a single run. Since increasing Bpol changes the equilibrium field connectivity (Fig. 1), the claim that a marginal (<5%) change in dipole strength controls eruption likelihood needs quantitative support and a sensitivity check. A quantitative comparison, for example of the maximum central-arcade height or the time of flux-rope formation against Bpol, would make the conclusion much more robust.","section":"§3.1, Fig. 2"}],"minor_comments":[{"comment":"The text reports velocities for three CMEs in the multiple-eruptions case, but the figure caption and legend appear to show only two curves (first and second CMEs). Please clarify the correspondence and label all curves consistently.","section":"Fig. 6"},{"comment":"The definition of t0 is confusing: it is first described as the steady-state time (~200 h) and then reset to zero. Please state explicitly that Eq. (8) is evaluated after the reset and remove the ambiguity.","section":"Section 2, Eq. (8)"},{"comment":"The abstract says the ANCH growth rate 'determines the likelihood' of CME eruptions, while Section 4 says it 'can serve as the most effective indicator.' These are different strengths of claim; please harmonize them and avoid causal wording unless the additional simulations recommended above are performed.","section":"Abstract and Section 4"},{"comment":"The footnote '* Released on March, 1st, 2021' appears to be a leftover from an earlier version and should be removed or updated.","section":"Page 1 footnote"},{"comment":"The first author name appears as 'Nitin V ashishtha' with a missing space or title; please check the author block for formatting errors.","section":"Author block"}],"recommendation":"major_revision","confidential_remarks":"The central claim in the abstract is not supported by the current experimental design because the proposed predictor is controlled by the same input parameter that defines the scenarios, and the thresholds are admitted to be resolution dependent. I would require additional simulations (at minimum a resolution study and runs that break the collinearity between ANCH slope and shear amplitude) or a substantially softened claim. The paper fits the journal's scope and the numerical work appears competent, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The poloidal-field scan (2.2–2.5 G, a <5% change) is a real, clean result: it suppresses the eruption and supports the idea that weaker background fields in Cycle 24 allow more weak CMEs. The ANCH-slope comparison is a legitimate extension of Talpeanu et al. 2020, but the causal claim attached to it is not supported by the evidence.\n\nWhat is actually new and good: the narrow poloidal-field scan is something I have not seen in prior breakout-model work. The method is standard AMRVAC, with careful boundary conditions and a shear profile motivated by observed coronal velocities. The authors also transparently note that the shear threshold depends on grid resolution, which is more honest than most papers in this area.\n\nSoft spots, in order of importance. First, the three eruption scenarios differ only by v0 (35.8, 36.2, 38.8 km/s), one run each. The reported ANCH slopes (0.412, 0.498, 0.645) increase monotonically with v0. Any monotone diagnostic of injected helicity—magnetic energy slope, TUCH slope—would order the runs the same way. So the design cannot separate “ANCH growth rate determines eruption” from “the prescribed shear amplitude determines both.” That is not fatal if framed as a candidate predictor, but the abstract says “establish,” which is too strong. Second, there is no resolution study or error bars for the slopes, and the authors admit the regime boundaries move with resolution; slopes separated by 0.1–0.2 simulation units may not survive. Third, the poloidal-field result is mostly field-line figures; no quantitative eruption heights, speeds, or energy release are reported for the 2.3–2.5 G cases. Fourth, the forecast link to SHARP is unquantified: SHARP ANCH is a photospheric surface quantity, while the simulated ANCH is a coronal volume integral of B·J, and no mapping is supplied.\n\nAlso note an internal inconsistency: the abstract says they “establish that the growth rate of absolute net current helicity is the crucial factor,” whereas the summary says it “can serve as the most effective indicator.” The second is defensible; the first is not.\n\nBottom line: the paper deserves a serious referee. It is a plausible diagnostic hypothesis, not an established physical criterion. I would ask for a resolution study, multiple runs per parameter, quantitative poloidal-scan metrics, and softer causal language. The poloidal-field suppression result alone is worth publishing.","headline":"A solid 2.5D breakout CME parameter study whose ANCH-growth-rate claim overreaches a single-run design; the poloidal-field suppression result is the more defensible new piece.","tokens_in":15891,"tokens_out":2634,"would_cite":false,"duration_ms":28030,"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 growth rate of absolute net current helicity, not its magnitude, determines whether a sheared solar arcade erupts as a CME.","keywords":["coronal mass ejection","breakout model","magnetohydrodynamics","current helicity","absolute net current helicity","solar cycle 24","poloidal magnetic field","flux rope eruption"],"falsifier":"Take one of the three shear speeds, for example 36.2 km/s, and rerun the 2.5D breakout simulation at twice and four times the adaptive-mesh refinement resolution; if the eruption classification flips or the ordering of ANCH slopes (multiple greater than single greater than failed) changes, the central claim that ANCH growth rate controls eruption likelihood is not robust.","tokens_in":14876,"feed_emoji":"☀️","tokens_out":7044,"duration_ms":67361,"temperature":0.7,"pith_summary":"This paper uses 2.5-dimensional magnetohydrodynamic simulations of the solar corona to ask what controls whether a sheared magnetic arcade erupts as a coronal mass ejection. Working within the breakout model, the authors vary two things: the strength of the background poloidal (global dipole) magnetic field and the amount of helicity injected by shearing the base of a central arcade. They find that increasing the background dipole field by less than five percent is enough to suppress flux-rope formation and eruption, and that the same shearing mechanism produces failed, single, or multiple eruptions depending on the shear speed. Their central claim is that the growth rate of the absolute net current helicity in the arcade, not just its magnitude, is the quantity that separates eruptive from non-eruptive cases. If this carries over to observations, tracking how fast a region's net current helicity is rising could be a practical eruption precursor, and weaker global solar magnetic fields would naturally produce more weak CMEs, as was seen in Solar Cycle 24.","feed_headline":"Fast net helicity growth signals a CME, simulations show","feed_subtitle":"A 2.5D breakout model ties eruption chance to the rise rate of a region's absolute net current helicity.","key_machinery":"The load-bearing mechanism is breakout reconnection in a multipolar quadrupole-like arcade: shear flow at the base of the middle arcade builds the azimuthal magnetic field and magnetic pressure; the expanding arcade flattens the overlying X-point, reconnection removes the restraining field, and a flux rope can escape. The diagnostic that carries the argument is the absolute net current helicity, the magnitude of the volume-integrated $\\mathbf{B}\\cdot\\mathbf{J}$ over the arcade region, and in particular its time derivative before the first breakout. The background dipole field is the control parameter that sets the height and connectivity of the arcade, turning a less-than-five-percent field increase into a qualitative change in eruption outcome.","core_discovery":"On its own terms, the paper establishes that the breakout CME outcome is set by the balance between helicity injection at the base and the confinement supplied by the background poloidal magnetic field. In simulations with a polar field of 2.2 G, the sheared central arcade forms a flux rope and erupts; at 2.3 G the arcade rises but no flux rope forms; at 2.4–2.5 G even the rise is suppressed, despite a configuration that is otherwise unchanged. With the polar field fixed at 2.2 G, maximum shear speeds of 35.8, 36.2, and 38.8 km s$^{-1}$ yield a failed eruption, a single eruption, and multiple eruptions. Computing the absolute net current helicity $\\mathrm{ANCH} = |\\sum_i \\int \\mathbf{B}\\cdot\\mathbf{J}\\,dV|$ across the arcade, the authors find that the slope of ANCH up to the first breakout reconnection orders the three cases: 0.645 (multiple), 0.498 (single), 0.412 (failed) in simulation units. In the multiple-eruption case, the second eruption is preceded by a steeper second rise (0.801), so the rate of increase, not the peak value, is the discriminator. The authors conclude that the ANCH growth rate is the crucial factor determining the likelihood of eruption.","pith_inferences":["Inference: the paper leaves implicit that observational predictors built on SHARP magnetograms should use the time derivative of absolute net current helicity, for example a 6 to 24 hour slope, instead of the snapshot value; this is a direct testable extension.","Inference: because the regime boundaries sit on shear-speed differences of 0.4 to 2.6 km/s and depend on grid resolution, higher-resolution or three-dimensional simulations will likely shift the exact threshold speeds even if the slope ordering survives.","Inference: if ANCH growth rate is the controlling factor, then helicity injection history, not just total injected helicity, should enter flare and CME forecasting; an active region that gains helicity quickly should be weighted more heavily than one that gains it slowly.","Inference: connecting this to Solar Cycle 24, the paper's mechanism suggests the weak-CME excess is physical rather than purely a detection-cadence artifact; a direct observational check would compare ANCH slopes from magnetograms in Cycle 24 and Cycle 23 active regions."],"forward_implications":["In an active region, a time series of absolute net current helicity should separate eruptive from non-eruptive cases by slope rather than by instantaneous value, giving a possible space-weather precursor.","Reducing the background poloidal field strength by a few percent makes eruptions easier at the same shear input, consistent with the observed excess of weak CMEs in Solar Cycle 24 under a weakened global magnetic field.","All breakout CMEs in the investigated domain are slow CMEs, below about 220 km/s, suggesting the mechanism is naturally a source of weak, slow CMEs rather than fast ones.","The same physical mechanism can produce failed, single, or multiple eruptions under shear-speed differences of only a few km/s, so predicting an eruption requires more than knowing the magnetic configuration at one instant."],"supporting_citations":[{"why":"Supplies the breakout model that the simulation implements.","marker":"Antiochos et al. 1999"},{"why":"Provides the multipolar arcade vector potential and setup used to build the coronal configuration.","marker":"van der Holst et al. 2007"},{"why":"Earlier breakout simulations with a similar arcade and background-field configuration that this paper extends.","marker":"Zuccarello et al. 2012"},{"why":"Similar simulation framework and shear-induced eruptions whose results this paper's slight-variation findings echo.","marker":"Talpeanu, D.-C. et al. 2020"},{"why":"Observational link between reduced solar and heliospheric pressure in Cycle 24 and more weak CMEs, the context the poloidal-field result supports.","marker":"Gopalswamy et al. 2015"},{"why":"Observational evidence that weaker polar fields allowed more eruptions in Cycle 24.","marker":"Petrie 2015"},{"why":"The rival observational-bias explanation that the authors' physical mechanism argues against.","marker":"Wang & Colaninno 2014"},{"why":"Defines the SHARP data products used to compute helicity-based predictors observationally.","marker":"Bobra et al. 2014"}],"fun_headline_variants":["Net helicity growth rate sets CME eruption odds","Helicity rise speed, not peak, predicts CME eruptions","Poloidal field strength and helicity injection steer CMEs","Helicity growth rate sorts failed, single, and multiple CMEs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The magnetic reconnection in the simulation is numerical in origin, tied to the grid spacing rather than to a physical resistivity model, so the tiny shear-speed thresholds that define the eruption regimes could shift with resolution.","fun_headline_variants_meta":{"raw":{"variants":["Net helicity growth rate sets CME eruption odds","Helicity rise speed, not peak, predicts CME eruptions","Poloidal field strength and helicity injection steer CMEs","Helicity growth rate sorts failed, single, and multiple CMEs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001091,"raw_usage":{"total_tokens":4657,"prompt_tokens":1142,"completion_tokens":3515,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":758,"completion_tokens_details":{"reasoning_tokens":3444}},"tokens_in":758,"tokens_out":3515,"duration_ms":27270,"temperature":1.0,"reasoning_tokens":3444,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:10:58.939054+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take one of the three shear speeds, for example 36.2 km/s, and rerun the 2.5D breakout simulation at twice and four times the adaptive-mesh refinement resolution; if the eruption classification flips or the ordering of ANCH slopes (multiple greater than single greater than failed) changes, the central claim that ANCH growth rate controls eruption likelihood is not robust.","supporting_citations":[{"cited_title":"P., Meliani, Z., & Poedts, S","cited_arxiv_id":null,"evidence_quote":"Earlier breakout simulations with a similar arcade and background-field configuration that this paper extends."},{"cited_title":"2015, Astrophysical Journal Letters, 804, doi: 10.1088/2041-8205/804/1/L23","cited_arxiv_id":null,"evidence_quote":"Observational link between reduced solar and heliospheric pressure in Cycle 24 and more weak CMEs, the context the poloidal-field result supports."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Observational evidence that weaker polar fields allowed more eruptions in Cycle 24."},{"cited_title":"M., & Colaninno, R","cited_arxiv_id":null,"evidence_quote":"The rival observational-bias explanation that the authors' physical mechanism argues against."}],"review_version":2}