{"id":"68ae8fb8-4eaf-4fad-bf0c-e137ffed9428","arxiv_id":"2607.29523","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Runaway electrons can form in stellarators during fast temperature collapses if the plasma current is large, but avalanche growth is weaker than in tokamaks due to elongation and lower currents.","lead":"This paper extends the DREAM simulation code to stellarators and scans scenarios in which a stellarator plasma loses its temperature quickly. It finds that runaway electrons can be generated in reactor-scale stellarators with large bootstrap currents, but usually much less severely than in tokamaks.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Plasma current is scanned against a fixed stellarator equilibrium, leaving the paper's quantitative runaway thresholds and avalanche comparison unvalidated.","rationale":"The paper is a well-structured pilot study with a clear central claim: significant runaway generation can occur in stellarators under sufficiently fast, low-temperature, high-current thermal collapse scenarios, but avalanche multiplication is weaker than in tokamaks due to lower currents and stronger elongation. The analysis is careful in many respects: the stellarator extension of DREAM is based on general forms of Faraday's and Ampère's laws, validation against the tokamak version is reported, the profile and boundary-condition sensitivities are explored, and the authors openly list limitations including the single configuration and the prescribed temperature evolution. The reader's weakest-assumption choice—inconsistent current density and fixed magnetic equilibrium—is indeed the most load-bearing issue for the quantitative conclusions. The paper acknowledges this limitation, but it is not merely a caveat: the scanned plasma current is not an independent parameter in a stellarator, and the avalanche dynamics are sensitive to the poloidal flux and geometry that a self-consistent equilibrium would alter. The qualitative possibility claim is robust, and the comparison with an elongated tokamak provides useful physical insight, so a rejection or a substantial downgrade is not warranted. However, because the quantitative thresholds and the 'weaker avalanche' assertion depend on the unphysical fixed-equilibrium scan, the conditional verdict is appropriate. A targeted recomputation with self-consistent equilibria would settle whether the concern materially changes the numbers.","tokens_in":15019,"tokens_out":9518,"duration_ms":88360,"concrete_test":"Recompute the highest-risk simulation points (e.g., Ip = 5 and 10 MA, Tfin = 1 and 10 eV, τ = 0.1 and 2 ms) using stellarator equilibria that are self-consistently solved with DESC for each prescribed current-density profile and a consistent pressure profile, updating the metric coefficients, rotational transform, and poloidal flux boundary condition. Compare the resulting runaway conversion fractions and the Eq. (17) avalanche metric with Figs. 4-6 and 11. If the runaway current changes by more than a factor of ~2 in any case, or if the threshold for >100 kA shifts by more than ~2 MA, the paper should explicitly state that its quantitative thresholds are conditional on the equilibrium model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims—the 'two out of three' conditions, the tTQ < 9 ms threshold, the >100 kA conversion boundaries, and the conclusion that avalanche multiplication is weaker in stellarators than tokamaks—rest on simulations in which the stellarator magnetic equilibrium is kept fixed while the plasma current is varied up to 10 MA (Sec. III, Figs. 1, 3). The paper explicitly states that 'the magnetic configuration is not consistent with the plasma current densities used' and that 'consistency between the plasma pressure and current density is not enforced.' This is not merely a technical simplification: in a reactor-relevant quasi-axisymmetric stellarator, the plasma current is predominantly a bootstrap current determined self-consistently by the pressure gradient and the magnetic geometry. A 10 MA current in a device with a = 1.7 m and R0 = 8.8 m produces a poloidal field on the order of 1 T, comparable to the confining field, so the flux surfaces, rotational transform, and elongation would change substantially. Since the avalanche gain and the proposed metric (Eq. 17) depend explicitly on gθθ/g and V', the fixed-equilibrium assumption can shift the conversion fractions in Figs. 4-6 and the comparison in Figs. 7 and 11. The qualitative finding that runaway generation is possible under fast, low-temperature, high-current conditions is likely robust, because the induced electric field follows from current decay regardless of the exact equilibrium. But the specific thresholds and the quantitative 'weaker avalanche' statements are not yet supported until the equilibrium is updated consistently with the current and pressure.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports a pilot numerical study of runaway-electron generation during temperature collapse in a reactor-scale quasi-axisymmetric stellarator. The authors extend the DREAM disruption code to accept general stellarator equilibria from DESC, using flux-surface-averaged forms of Faraday's and Ampère's laws. They scan initial plasma current I_p = 0.1–10 MA, temperature decay time τ, post-collapse temperature T_fin, current-profile peak, and wall distance b/a, with prescribed exponential temperature decay and self-consistent current/flux evolution. They report conversion fractions up to ~100% for fast collapses; a 'two out of three' condition for significant runaway current; a threshold t_TQ < 9 ms for >1 MA runaway current; and a comparative study showing weaker avalanche multiplication in strongly shaped stellarator/elongated plasmas than in circular tokamak geometry. A metric for avalanche gain (Eq. 17) is proposed.","tokens_in":15363,"tokens_out":6531,"duration_ms":60299,"significance":"If the quantitative results are supported, this is a valuable and timely contribution: it opens a new line of inquiry for stellarator safety, provides a usable code extension, and gives a design-oriented metric. The authors are careful to validate the stellarator model against DREAM tokamak results and to state limitations. The qualitative conclusion that a fast, low-temperature, high-current collapse can produce significant runaway current is plausible and not circular. However, the central quantitative claims are currently not fully validated because the magnetic equilibrium is kept fixed while the plasma current is varied by two orders of magnitude, and because the text gives inconsistent t_TQ thresholds. The paper's value would be substantially increased by recomputing or reframing those thresholds with self-consistent equilibria, or by explicitly presenting the results as proof-of-principle trends rather than as device-specific predictions.","major_comments":[{"comment":"The paper varies I_p from 0.1 to 10 MA while retaining a single fixed equilibrium; the text itself states 'the magnetic configuration is not consistent with the plasma current densities used' and that 'consistency between the plasma pressure and current density is not enforced.' This is a load-bearing issue for the quantitative thresholds. For a = 1.7 m and I_p = 10 MA, the poloidal field μ0 I/(2π a) ≈ 1.2 T, comparable to the confining field, so flux surfaces, rotational transform, and elongation would change substantially. Since the avalanche gain and the proposed metric in Eq. (17) depend on V' and ⟨g_θθ/g⟩, the conversion fractions in Figs. 4–6 and the comparisons in Figs. 7 and 11 could shift. Please recompute with self-consistent equilibria (e.g., DESC equilibria with consistent current and pressure profiles for each I_p) or restrict the quantitative claims to a fixed-background co","section":"Section III, Figs. 4–6 and Eq. (17)"},{"comment":"The t_TQ threshold for significant runaway generation is stated inconsistently: 't_TQ ≲ 2.8 ms' appears in the discussion of Fig. 5, 'If t_TQ < 9 ms ... (I_re ≳ 1 MA)' appears after Fig. 6, and 't_TQ ≲ 4.7 ms' appears in the Discussion. Because this threshold is one of the paper's main quantitative conclusions, the discrepancy must be resolved and tied to a precise definition of 'significant' (an I_re threshold) and to the parameter range (I_p, T_fin). Please correct the typographical/numerical inconsistency and state which figure supports each threshold.","section":"Section IV (Fig. 5 discussion and Fig. 6) vs. Section V"},{"comment":"The conclusion that avalanche multiplication is weaker in stellarators than in tokamaks is based on a single 3D configuration compared with a circular and an elongated tokamak. The elongated tokamak is matched via Eq. (16), but the stellarator equilibrium is not consistent with the 10 MA current, and the matching criterion is a flux-surface average that may not capture the 3D variation. The simulated difference between the stellarator (780 A) and the elongated tokamak (540 A) is modest; the strong reduction is relative to the circular case. The claim should be stated more carefully as 'strongly shaped plasmas, including elongated tokamaks, have much weaker avalanche multiplication than circular tokamaks,' and ideally supported by a self-consistent comparison.","section":"Section IV (Fig. 7) and Section V"}],"minor_comments":[{"comment":"The caption lists 'I_p = 10 MA (black dotted), I_p = 5 MA (purple dashed), and I_p = 10 MA (red solid).' Given the text states I_p ∈ {1, 5, 10} MA, the red solid curve is presumably I_p = 1 MA. Please correct the duplicate 10 MA entry.","section":"Fig. 8 caption"},{"comment":"The caption says 'different curves, of varying colours' without a legend or clear color mapping. A legend or explicit line-style/color table would help the reader identify I_p = 1, 5, 10 MA.","section":"Fig. 5"},{"comment":"Please define the units of the metric explicitly. The numerator includes I(ρ) normalized to 1 MA, while the denominator has dimensions from ⟨g_θθ/g⟩ V'/a; stating the resulting units and the normalization of the integral would improve reproducibility.","section":"Eq. (17)"},{"comment":"Reference [33] points to a personal/public URL for the equilibrium. If possible, provide a permanent repository or DOI, or describe the equilibrium data in sufficient detail for reproduction.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is not fatally flawed; the qualitative conclusion that fast, low-temperature, high-current collapses can generate significant runaways is likely robust. The main issue is that the fixed-equilibrium inconsistency undermines the specific quantitative thresholds as presented, and the t_TQ threshold is stated inconsistently across sections. I would be willing to re-review a revised version that either recomputes the thresholds with self-consistent equilibria or clearly reframes all numerical thresholds as illustrative trends for a fixed background configuration."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThis paper is a worthwhile pilot study: it asks whether runaway electrons could be a problem in reactor-scale stellarators, and gives a cautious 'maybe under some conditions.' The authors extend the DREAM runaway simulation code to stellarator geometry, validate it against a tokamak case, and scan over initial current, post-collapse temperature, and thermal quench time. The main result — that fast, low-temperature collapses with multi-mega-ampere bootstrap currents can produce significant runaway current — is new and plausible. The avalanche metric (Eq. 17) is a nice practical output for configuration optimization.\n\nThe strongest part is the direct comparison of avalanche multiplication between a stellarator and circular/elongated tokamaks with the same current and volume. That cleanly demonstrates that elongation suppresses the avalanche, and gives a physical reason to expect stellarators to be less prone to the avalanche problem than tokamaks. The authors are also honest about the main limitations: one configuration, prescribed temperature decay, no radial transport.\n\nThe soft spot is the one the stress-test note flags: the magnetic equilibrium is held fixed while the plasma current is varied up to 10 MA. That matters, because in a quasi-axisymmetric stellarator the current is predominantly bootstrap, tied to pressure and geometry. A 10 MA current in this device would produce a poloidal field of order 1 T, enough to alter the flux surfaces and rotational transform. Since the avalanche gain depends on the metric components, the quantitative thresholds (tTQ < 9 ms, 'two out of three' conditions, the >100 kA boundaries) should be read as indicative, not precise. The qualitative conclusion that runaway generation is possible under fast, low-T high-Ip conditions is almost certainly robust, because the induced electric field follows from current decay. But 'weaker avalanche in stellarators' is also likely true, though the specific factor could shift with a self-consistent equilibrium.\n\nThere is also a minor internal inconsistency: the two-out-of-three condition gives tTQ <= 2.8 ms in the results section and tTQ <= 4.7 ms in the discussion. Easy to fix, but worth flagging to the authors.\n\nThe paper doesn't ship code or data, which limits reproducibility, though the equilibrium is linked.\n\nOverall: this is a solid pilot study that deserves a serious look. The right referee will ask for a self-consistent equilibrium case (or a sensitivity scan over equilibria) before the thresholds are quoted in design documents. I'd send it to peer review.","headline":"A credible pilot study showing stellarators could produce runaways under fast, cold, high-current collapses, but the quantitative thresholds rely on an equilibrium inconsistent with the imposed currents.","tokens_in":15834,"tokens_out":5821,"would_cite":true,"duration_ms":47921,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.55.Hc"],"model":"deepseek-v4-flash","headline":"A simulation study of stellarator temperature collapses finds that runaway electron generation is possible — especially with high initial current, a fast quench, and low final temperature — but avalanche multiplication is far weaker than in","keywords":["runaway electrons","stellarators","tokamaks","avalanche generation","thermal quench","bootstrap current","plasma disruption","magnetic confinement"],"falsifier":"A direct falsifier would be a measurement of post-collapse runaway current in a stellarator with a known 5–10 MA plasma current and a thermal quench faster than about 5 ms: the paper predicts conversion fractions of tens of percent, so observing less than a few kiloamperes would contradict it.","tokens_in":14914,"feed_emoji":"⚡","tokens_out":3793,"duration_ms":34747,"temperature":0.7,"pith_summary":"The paper asks whether stellarators, which normally lack a large plasma current, can generate relativistic runaway electrons during a sudden temperature collapse. Using an extended simulation code that treats stellarator geometry through generalized Faraday and Ampère laws, the authors scan over initial plasma current, thermal-quench time, and post-collapse temperature. They find significant runaway generation occurs when two of three conditions hold: current above about 5 MA, final temperature below about 20 eV, or a quench faster than a few milliseconds. However, even at high current, the avalanche gain is strongly suppressed by the elongation of stellarator cross-sections compared with circular tokamaks, so reactor-scale stellarators are predicted to be much less prone to damaging runaway beams.","feed_headline":"Stellarators can breed runaways, but far weaker than tokamaks","feed_subtitle":"High current, fast quenches, and low final temperature matter; elongation suppresses the avalanche.","key_machinery":"The central object is the generalized Faraday and Ampère laws written in flux-surface-averaged form valid for any toroidal configuration, implemented in the DREAM simulation code, together with a proposed avalanche-gain metric: the integral over minor radius of the enclosed current divided by the flux-surface-averaged metric factor ⟨gθθ/g⟩V′/a. The metric correlates linearly with the logarithm of the avalanche multiplication factor, letting one estimate and compare runaway avalanche strength across configurations without full simulation.","core_discovery":"The central claim is that runaway electrons, long assumed negligible in stellarators because there is no externally driven current, can be generated in a radiative temperature collapse when the bootstrap or other plasma current is large. The authors implement a fluid model for stellarators in the DREAM code and demonstrate that the avalanche mechanism is exponentially sensitive to current but reduced by plasma shaping: an elongated stellarator produces more than a thousand times less avalanche runaway current than a circular tokamak with the same current and volume. They also propose a simple metric based on Ampère's law that predicts the logarithmic avalanche gain for any configuration, and","pith_inferences":["Our inference: If the avalanche-gain metric is robust, it could be used as a cheap screening tool for proposed stellarator reactor concepts long before full disruption simulations are run.","Our inference: The strong dependence on thermal-quench time suggests that active control of impurity ingress (e.g., by ECRH) may be a more effective runaway mitigation strategy in stellarators than in tokamaks.","Our inference: Since the paper fixes the magnetic equilibria while varying current, a self-consistent equilibrium scan might show the avalanche suppression to be even stronger, because high bootstrap current would also alter the rotational transform and shaping."],"forward_implications":["If the claim is right, stellarator reactor designs should include runaway mitigation only in high-current, fast-quench scenarios, not as a general requirement.","The avalanche-gain metric could become a design constraint in stellarator optimization, favouring configurations with low enclosed current and strong shaping.","Thermal-quench time is the dominant lever: if quenches in large stellarators are slower than about 10 ms (as LHD data suggest), runaway currents stay below roughly 100 kA even at 10 MA initial current.","The comparison with elongated tokamaks implies that tokamak disruption studies should separate the effect of elongation from geometry-specific effects."],"fun_headline_variants":["Stellarators can breed runaways — but far weaker than tokamaks","Runaway electrons: stellarators not immune, but safer","Stellarators: runaways possible, yet thousand-fold weaker","Stellarator runaways: unlikely? Not quite, but less dangerous","Bootstrap current can spark runaways in stellarators"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The simulations hold the stellarator magnetic configuration fixed while scanning the plasma current, and do not enforce consistency between the current density and the plasma pressure; if self-consistent equilibria have different current profiles or poloidal flux, the quantitative runaway fractions could change.","fun_headline_variants_meta":{"raw":{"variants":["Stellarators can breed runaways — but far weaker than tokamaks","Runaway electrons: stellarators not immune, but safer","Stellarators: runaways possible, yet thousand-fold weaker","Stellarator runaways: unlikely? Not quite, but less dangerous","Bootstrap current can spark runaways in stellarators"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000642,"raw_usage":{"total_tokens":2775,"prompt_tokens":710,"completion_tokens":2065,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":454,"completion_tokens_details":{"reasoning_tokens":1977}},"tokens_in":454,"tokens_out":2065,"duration_ms":14211,"temperature":1.0,"reasoning_tokens":1977,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T05:17:29.746882+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct falsifier would be a measurement of post-collapse runaway current in a stellarator with a known 5–10 MA plasma current and a thermal quench faster than about 5 ms: the paper predicts conversion fractions of tens of percent, so observing less than a few kiloamperes would contradict it.","supporting_citations":[],"review_version":1}