{"id":"93ee2127-bc9e-4e2e-b4a8-e3602559e53a","arxiv_id":"2608.02121","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Simulations show that finite-orbit-width effects of runaway electrons, especially curvature drift, stabilize tearing modes and reduce stochastic transport at higher RE energies, counteracting the known destabilizing role of REs.","lead":"This paper simulates how finite drift-orbit width of runaway electrons changes the stability of magnetic islands in a tokamak disruption. Higher-energy runaway electrons are found to stabilize the instability and reduce particle losses, which matters for predicting how runaway beams terminate in ITER.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: standard ∇B cancellation is minor at pitch=0.999; core FOW mechanism supported by current-layer width and transport diagnostics.","rationale":"The paper's quantitative evidence for the FOW mechanism is strong. The reader's weakest_assumption about ∇B cancellation is not the most load-bearing because the ∇B drift is small at pitch=0.999 and the cancellation is a known consistency condition. The central claim is scenario-specific, but that is a limitation, not an error. Verdict remains CONDITIONAL pending artifacts, not due to a scientific flaw.","tokens_in":12020,"tokens_out":14299,"duration_ms":122264,"concrete_test":"Run the 20 MeV case with a fluid RE model that reproduces the same anisotropic pressure tensor but zero orbit width (REs tied to flux surfaces). If the (2,1) growth rate and current-layer width match the kinetic 20 MeV result, the stabilization is due to pressure/equilibrium changes rather than FOW; if they remain close to the low-energy fluid case, the FOW mechanism is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I reviewed the central claim that increasing RE energy stabilizes tearing modes via drift-orbit broadening of the current sheet. The paper provides direct evidence: the (2,1) current layer width increases non-monotonically with energy (Fig. 3), the linear growth rate decreases (Fig. 4), and the Chirikov parameter and transport coefficients decrease (Figs. 8-11). The ∇B cancellation mentioned in Section 2 is a standard consistency requirement between the particle current and the pressure-tensor coupling; at pitch=0.999 the ∇B drift is ~0.2% of the curvature drift, so even an incomplete cancellation would not alter the qualitative result. The main caveat is the absence of a fluid-RE control at matched pressure, but the non-monotonic layer width is a FOW signature that cannot be explained by pressure alone. Thus no load-bearing concern is identified.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents hybrid kinetic-MHD simulations using the JOREK code to study finite-orbit-width (FOW) effects of runaway electrons (REs) on resistive tearing-mode stability and nonlinear dynamics. REs are treated kinetically via full-f relativistic guiding-center Monte Carlo markers coupled to extended MHD. The authors find that increasing RE energy first narrows then broadens the (2,1) current-sheet width (Fig. 3), reduces the linear growth rate (Fig. 4), lowers nonlinear island saturation amplitudes (Fig. 7), and decreases the Chirikov parameter, connection length, and radial diffusion coefficients (Figs. 8, 10, 11). The central claim is that FOW effects stabilize tearing modes by preventing narrow current sheets on rational surfaces, and that this stabilization dominates over enhanced poloidal mode coupling in the investigated scenario.","tokens_in":12293,"tokens_out":8492,"duration_ms":71432,"significance":"If correct, the result is significant because it shows that kinetic FOW effects can qualitatively reverse the destabilizing role of REs predicted by zero-orbit-width (fluid) models, with direct implications for predicting disruption current-quench evolution and benign termination. The study uses multiple complementary diagnostics (layer width, growth rate, Poincare plots, connection length, diffusion, Chirikov parameter) that are internally consistent. The non-monotonic layer width in Fig. 3 is a falsifiable signature that cannot be explained by pressure alone. However, the absence of numerical convergence information and of a zero-orbit-width control at matched pressure leaves the central attribution only partially demonstrated.","major_comments":[{"comment":"The manuscript does not report the number of markers N, grid resolution, time step, or any convergence tests for the hybrid simulations. Since the central claim is a simulation outcome, the lack of these numerical details prevents the reader from assessing whether the observed trends (growth rate decrease, saturation reduction, diffusion decrease) are robust. Please provide the numerical parameters used and a convergence check (at least in marker number and resolution) for the key quantitative results.","section":"§2, Eq. (3); §4, Figs. 4, 7, 10"},{"comment":"The linear growth rate is shown to decrease monotonically with RE energy (Fig. 4), and this is attributed to FOW broadening of the current sheet. However, the RE pressure and the associated Grad-Shafranov shift also increase with energy (Fig. 1(a)), so the stabilizing trend could in part be an equilibrium (pressure/shift) effect rather than a direct FOW current-broadening effect. The 600 keV case in Fig. 3 provides a partial non-monotonic check for the layer width, but the growth rate for this case is not shown. Please provide a fluid-RE (zero orbit width) control at matched pressure, or at least the 600 keV growth rate, to demonstrate that the stabilization is specifically due to FOW broadening.","section":"§4.1, Figs. 3 and 4"},{"comment":"The radial diffusion coefficients are obtained by Gaussian fitting after excluding particles 'strongly affected by the magnetic island region or boundary transport', but the exclusion criterion is not defined quantitatively. The computed D values and the claimed decreasing trend with energy (Fig. 10f) may depend strongly on which particles are removed. Please specify the selection criterion, report the fraction of excluded particles for each energy, and show that the trend in D is robust to reasonable changes of the criterion.","section":"§4.3, Eq. (6) and Fig. 9"}],"minor_comments":[{"comment":"The phrase 'enhanced toroidal mode coupling' describing the (2,1)->(3,1) sideband is a misnomer: the (1,0) equilibrium perturbation couples poloidal harmonics m to m±1 for the same toroidal number n. It should read 'poloidal mode coupling'.","section":"§4.1, after Eq. (5)"},{"comment":"Please clarify what is meant by 'ψ structures' - presumably the perturbed poloidal flux structures of the (2,1) and (3,1) components.","section":"Fig. 5 caption"},{"comment":"The quantity v_RE is used but not defined. Please define it in the text following the equation (e.g., the parallel RE velocity).","section":"Eq. (5)"},{"comment":"The axis label and units of the diffusion coefficient D are not described in the text. Please specify the normalization (e.g., m^2/s) and the radial coordinate used.","section":"§4.3, Fig. 10(f)"},{"comment":"The manuscript does not include a data availability statement or mention whether input files/scripts are available. Given the simulation-based nature of the work, a statement on data/code availability would improve reproducibility.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper presents an interesting and potentially important simulation study. The main concerns are reproducibility (missing numerical parameters/convergence tests) and attribution of the stabilizing effect to FOW rather than equilibrium pressure effects. These are addressable with additional simulations or targeted diagnostics. The paper should be suitable for publication after major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real new result here is that finite-orbit-width (FOW) effects reverse the previously reported destabilizing role of runaway electrons (REs) on tearing modes: as RE energy increases, the linear growth rate drops, saturation amplitudes shrink, and the stochastic region narrows. That is a genuinely new and physically interesting picture, and the paper does a good job of supporting it with multiple diagnostics. The (1,0) current perturbation from orbit displacement, the energy-dependent current-sheet broadening, the Chirikov parameter, connection length, and transport coefficients all point in the same direction. The non-monotonic layer width (narrowing as velocity approaches c, then broadening again with FOW) is a nice signature that is hard to explain by pressure effects alone.\n\nThe hybrid kinetic-MHD model itself is state of the art: full-f Monte Carlo REs coupled self-consistently to JOREK MHD, with the GC equations stated clearly. The ∇B cancellation assumption in Section 2 is a standard consistency requirement, and at pitch=0.999 the ∇B drift is tiny relative to curvature drift, so even an incomplete cancellation would not change the qualitative conclusion. The stress-test note holds up: no load-bearing flaw found.\n\nSoft spots, in proportion: first, no code or data are provided. For a simulation paper claiming a new physical mechanism, convergence checks, numerical settings, and ideally input files are needed for a referee to verify the result. Second, the study is scenario-specific: one circular equilibrium, one q profile, mono-energetic REs. That is fine for a first demonstration, but the authors are appropriately careful about not over-generalizing. Third, there is no fluid-RE control run at matched pressure to fully separate pure FOW effects from pressure-gradient effects; the non-monotonic layer width argues strongly that FOW is doing the work, but a matched-pressure comparison would silence the remaining doubt. Fourth, the Gaussian-fitting transport coefficients exclude island/boundary particles, which makes them effective indicators rather than complete transport measures; the authors do note this.\n\nThe citation pattern is fine; self-citations to prior JOREK hybrid papers are appropriate since the model is described there. This is a credible extension of Refs. [13,14], not a rehash.\n\nWho is this for? Fusion researchers working on RE mitigation and disruption modeling, especially those using MHD codes. It deserves a serious referee. My recommendation: send it to peer review, but ask for the artifacts (JOREK input decks, post-processing scripts, at least convergence plots) before final acceptance.","headline":"A solid scenario-specific simulation study showing that finite-orbit-width effects stabilize runaway-electron-driven tearing modes with increasing energy; the mechanism is plausible and the diagnostics line up, but reproducibility is limited by no artifacts.","tokens_in":12753,"tokens_out":1502,"would_cite":true,"duration_ms":18512,"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":"Higher-energy runaway electrons stabilize the tearing modes of a post-disruption tokamak, because drift-orbit displacement smears out the current sheets that drive the instability.","keywords":["runaway electrons","tearing modes","finite-orbit-width effects","hybrid kinetic-MHD model","tokamak disruption","mode coupling","stochastic magnetic fields","particle transport"],"falsifier":"Run a version of the simulation that keeps the ∇B drift and the RE magnetization current separately (breaking the cancellation) and compare the width of the (2,1) current layer and the (1,0) current component as functions of RE energy. If the layer does not broaden or the (1,0) component does not grow with energy, the proposed stabilisation mechanism is not supported. Alternatively, measure the (1,0) current perturbation in an experiment where the RE energy is scanned, and check that the predicted ordering with energy holds.","tokens_in":12012,"feed_emoji":"⚡","tokens_out":5974,"duration_ms":45049,"temperature":0.7,"pith_summary":"Runaway electrons (REs) created in tokamak disruptions normally make tearing modes more unstable than a purely Ohmic plasma would be, according to earlier fluid-like models. This paper argues that this conclusion changes once the finite width of the particles' drift orbits is included: high-energy REs wander off the magnetic flux surfaces and cannot sustain the narrow current sheets that drive the tearing instability. The orbit deviation also introduces a (1,0) perturbation of the equilibrium current, coupling the dominant (2,1) mode to (3,1) and (1,1) sidebands. For the case studied, the stabilizing effect wins: as RE energy rises from 3 MeV to 30 MeV, islands saturate smaller, stochastic regions shrink, connection lengths grow, and radial particle transport falls. The work thus shows that ignoring drift-orbit width overestimates the threat that RE beams pose to disruption mitigation.","feed_headline":"High-energy runaways stabilize tokamak disruption modes","feed_subtitle":"Drift-orbit smearing of current sheets reverses the known destabilizing effect and cuts particle loss.","key_machinery":"The mechanism is the relativistic guiding-center drift of the REs, dominated by curvature drift, which displaces their orbits from the magnetic flux surfaces by several centimetres at 20 MeV. This displacement produces a strong (1,0) component of the RE current, coupling (m,n) modes to (m±1,n) sidebands, and it broadens the resistive current sheet at the rational surface, which is identified as the stabilizing effect. The model self-consistently couples a full-f Monte-Carlo RE population to an extended-MHD fluid via the RE pressure and by excluding the kinetically evolved RE current from the Ohm's law.","core_discovery":"The central claim is that finite-orbit-width (FOW) effects, caused mainly by the curvature drift of relativistic REs, stabilise the resistive tearing modes in a post-disruption tokamak. In the linear regime, the current layer at the (2,1) rational surface broadens with RE energy because orbits deviate from flux surfaces, reducing the growth rate toward the Ohmic value. In the nonlinear regime, the (2,1) island saturation amplitude, the Chirikov overlap parameter, and the radial diffusion coefficient all decrease as the RE energy is increased, meaning the magnetic field becomes less stochastic and RE transport is suppressed. These results are obtained with a hybrid kinetic-fluid model in whic","pith_inferences":["If the assumed cancellation between the ∇B drift and the RE magnetization current is not exact, the (1,0) perturbation and the current-sheet broadening could change in magnitude or sign, so the stabilisation should be tested against a model that keeps both drift terms.","The clean energy dependence found here suggests an experimental probe: comparing MHD activity in RE beams with different mean energies, e.g., by varying the electric field during the current quench, should show weaker modes and less transport at higher energies.","In a realistic broad energy distribution, low-energy REs would still destabilise while high-energy ones stabilise; the net effect would depend on the spectral weighting, so transport codes should retain the full energy resolution.","The same orbit-width-smoothing principle should apply to other current-driven instabilities and rational surfaces, suggesting that kinetic orbit width is a generic spatial regularisation scale in post-disruption plasmas."],"forward_implications":["The linear growth rate and nonlinear saturation of the (2,1) tearing mode both decrease with RE energy, approaching the pure-Ohmic reference at high energy.","A (1,0) equilibrium-current perturbation appears at high RE energies, generating (m±1,n) sidebands such as the (3,1) mode and enriching the mode spectrum.","The Chirikov parameter, connection length, and radial diffusion coefficients all show reduced stochasticity and transport at higher RE energy.","Fluid RE models that assume zero drift-orbit width overestimate the destabilising effect of REs and would mispredict the nonlinear evolution of a RE beam.","RE beam termination strategies that rely on MHD-triggered stochastic losses may be delayed or weakened for high-energy beams, changing the expected heat loads."],"fun_headline_variants":["Drift-orbit smearing stabilizes runaway tearing modes","Kinetic drift orbits reverse runaway destabilization","Finite-orbit effects tame tokamak disruption modes","Curvature drift widens current sheets, stabilizes runaways"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The model assumes the ∇B drift contribution to the RE current is exactly canceled by the RE magnetization current, so finite-orbit-width effects on the fluid come only from curvature drift; if the cancellation is incomplete, the magnitude or sign of the energy-dependent stabilisation could change.","fun_headline_variants_meta":{"raw":{"variants":["Drift-orbit smearing stabilizes runaway tearing modes","Kinetic drift orbits reverse runaway destabilization","Finite-orbit effects tame tokamak disruption modes","Curvature drift widens current sheets, stabilizes runaways"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000888,"raw_usage":{"total_tokens":3741,"prompt_tokens":887,"completion_tokens":2854,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":2788}},"tokens_in":631,"tokens_out":2854,"duration_ms":15996,"temperature":1.0,"reasoning_tokens":2788,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T14:46:41.947452+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a version of the simulation that keeps the ∇B drift and the RE magnetization current separately (breaking the cancellation) and compare the width of the (2,1) current layer and the (1,0) current component as functions of RE energy. If the layer does not broaden or the (1,0) component does not grow with energy, the proposed stabilisation mechanism is not supported. Alternatively, measure the (1,0) current perturbation in an experiment where the RE energy is scanned, and check that the predicted ordering with energy holds.","supporting_citations":[],"review_version":1}