{"id":"0374dae7-85f6-4020-a686-51d3d36885fc","arxiv_id":"2511.04051","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Fishbone-driven zonal flows suppress electromagnetic ITG turbulence and reduce ion heat transport to near-neoclassical levels in global gyrokinetic simulations of EAST-like plasmas.","lead":"This paper uses the first global gyrokinetic simulations to couple fishbone instability and ion-temperature-gradient turbulence in a tokamak, finding that fishbone-driven zonal flows create radial electric field shear that suppresses the turbulence and cuts ion heat transport to near-neoclassical levels. If correct, it provides a mechanism for fishbone-induced internal transport barriers and could affect predictions for ITER-like burning plasmas.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zonal-flow shearing rate may be dominated by unresolved fine-scale radial structure; convergence of the second-derivative shear is unshown.","rationale":"The reader's weakest assumption focused on numerical convergence and marker noise damping ITG turbulence. My concern is more specific: the fishbone-driven zonal-flow shearing rate, which is the mechanism responsible for suppression, is computed from a second radial derivative of a zonal potential that is explicitly described as having fine-scale radial structure. Second derivatives amplify short-wavelength components, and with only 100 radial grid points, the fine structure may be at the grid scale, making the shear overestimated. This is load-bearing because the suppression claim depends on the shearing rate exceeding the ITG growth rate; if that threshold is an artifact of under-resolution, the central mechanism fails. The imposed-Er control does not resolve this because it uses the same uncertain Er profile. The paper's assertion of meticulous convergence is not backed by a scan, and the ECE validation does not test the zonal-flow shear. I agree partially with the reader: we both identify numerical fidelity as the weak point, but I pinpoint the specific derived quantity that is most vulnerable. The verdict remains CONDITIONAL, as the paper otherwise has strong internal consistency (the imposed-Er run reproduces coupled suppression) and experimental correlation, but the quantitative threshold needs verification.","tokens_in":8823,"tokens_out":6174,"duration_ms":60330,"concrete_test":"Run the fishbone-only GTC case at Nψ=200 (and/or 1000 markers/cell) and recompute the zonal potential δφ00 and ω_E×B profile from Fig. 2(c). If the peak shearing rate drops below γ_ITG or the fine-scale radial structure loses amplitude when the grid is refined, the suppression mechanism is a numerical artifact. Also plot the radial power spectrum of δφ00 and confirm it decays before the grid Nyquist wavenumber; if it is cut off at the grid scale, the second derivative is unreliable. If the shearing rate remains >γ_ITG by the same margin at higher resolution, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central suppression mechanism rests on the fishbone-driven zonal-flow shearing rate ω_E×B ∝ ∂²δφ00/∂ψ² exceeding the ITG growth rate (Fig. 2c). The authors report a 'fine-scale radial structure' in the zonal potential, attributed to the zonal electron density response, but provide no evidence that this structure is resolved. With Nψ=100 grid points over r/a=0–0.63 (Δr/a≈0.0063), the second radial derivative amplifies grid-scale oscillations; if the fine-scale structure is not converged (or is marker noise), the computed shearing rate could be artificially large, making the suppression of ITG turbulence (δφ from 1e-3 to 1e-4) appear stronger than physical. The imposed-Er control (Fig. 4d) uses the same Er profile, so it inherits the same risk. The one-line claim 'Numerical convergence studies are meticulously achieved' is not substantiated with a scan, and the ECE validation (Fig. 2d) checks the fishbone temperature perturbation, not the zonal-flow shear. Thus the key quantitative threshold — ω_E×B > γ_ITG — is not robustly established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports global gyrokinetic GTC simulations of electromagnetic ITG turbulence and the n=1 fishbone instability using EAST discharge #93890 equilibrium. The central claim is that, in coupled nonlinear simulations, fishbone-driven zonal flows/radial electric fields at saturation produce E×B shearing rates that exceed the ITG linear growth rate, suppress ITG fluctuations from δφ~1e-3 to ~1e-4, and reduce the ion heat conductivity close to the neoclassical level. The mechanism is isolated by three related simulations: ITG+EPs, ITG+EPs with an imposed fishbone-driven Er, and fully coupled fishbone+ITG. Comparisons with ECE, ONETWO power balance, and Doppler reflectometry on EAST discharges are presented as experimental support.","tokens_in":9130,"tokens_out":6111,"duration_ms":57243,"significance":"If the central suppression mechanism holds, the paper establishes a new cross-scale channel by which a macroscopic EP-driven mode improves core ion confinement, with implications for ITER hybrid scenarios and for the interpretation of fishbone-associated ITBs in EAST and other devices. The study is a genuine computational first (global gyrokinetic coupling of fishbone and ITG) and has a clean control structure: the imposed-Er run and coupled run agree, and the fishbone-only ECE validation anchors the mesoscale part. The main weakness is that the correctness of the load-bearing shear-rate threshold rests on unresolved numerical convergence claims.","major_comments":[{"comment":"The central threshold ω_E×B > γ_ITG is computed from ∂²δφ00/∂ψ² on an Nψ=100 grid over r/a=0–0.63 (Δr/a≈0.0063). The text asserts 'Numerical convergence studies are meticulously achieved' but no scan is shown. Since the zonal potential has a fine-scale radial structure (attributed to the zonal electron density response), the second radial derivative can be dominated by grid-scale oscillations, and marker noise or insufficient radial resolution could artificially inflate the shearing rate. The imposed-Er control run in Fig. 4(d) uses the same Er profile and therefore inherits this risk. ECE validation in Fig. 2(d) checks the fishbone temperature perturbation, not the zonal-flow shear. Please provide a convergence study in Nψ (and, if possible, in marker number and time step) demonstrating that the peak shearing rate and the resulting transport reduction are robust. This is load-bearing fo","section":"Simulation setups / Fig. 2(c)"},{"comment":"Figure 6 is presented as direct evidence that fishbone-driven zonal flows suppress turbulence, but the Doppler reflectometry measurement shows a Doppler shift during fishbone bursts and a correlated reduction of turbulence intensity; a Doppler shift alone does not uniquely identify the zonal-flow shearing rate or its radial structure. Without either a quantitative shearing-rate estimate from the DR data or a measurement of the zonal-flow radial wavelength, Fig. 6 remains correlative. This does not invalidate the simulation-based conclusion, but the abstract's phrasing 'agrees well with experimental observations' should be softened or supported by a more quantitative comparison.","section":"Fig. 6 / Experimental evidence"},{"comment":"The reduction of χi to 'close to the neoclassical level' is presented without error bars or sensitivity analysis; the three curves in Fig. 4(b) come from single realizations. A convergence/robustness statement (e.g., variation with grid resolution, marker number, or time averaging window) would strengthen the quantitative claim.","section":"Fig. 4(b) / transport comparison"}],"minor_comments":[{"comment":"Typos and garbled words should be corrected: 'ITER-likc', 'cross-cale', 'marcoscopic', 'scries', 'bule line', 'conductvity', and 'balck' in the Fig. 6 caption.","section":"Throughout"},{"comment":"Define the notation 'n=m=0' and clearly state the variables in the shearing-rate expression; an equation number for ω_E×B would help reproducibility.","section":"Fig. 2(a) and Fig. 2(c)"},{"comment":"The numerical values for χi appear inconsistently ('~0.2 m2/s' and '~0.3 m2/s') near the comparison with ONETWO; please unify the notation and label each curve explicitly in the text.","section":"Fig. 4(b) and text"},{"comment":"The claim of 'first cross-scale simulations' should be qualified relative to previous AE–ITG coupling studies (e.g., Refs. [14,16]); clarify that the novelty is specifically the fishbone–ITG coupling.","section":"Introduction"},{"comment":"Axes and units in Fig. 6 are not fully labeled; the Doppler-shift trace needs its vertical scale and relation to zonal-flow shear stated in the caption.","section":"Fig. 6"}],"recommendation":"major_revision","confidential_remarks":"The main gate for publication is the convergence evidence. If the authors can supply a resolution scan (radial grids, marker number, time step) showing that the fine-scale zonal-flow shearing rate and the χi reduction are stable, the paper would be a strong candidate for acceptance. The correlative experimental evidence is secondary and should be framed as supporting, not conclusive."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Zhuo — quick read of Ma et al. The thing to know: the central claim — fishbone-driven zonal flows shear apart ITG eddies and cut ion heat transport to neoclassical levels — is supported by a cleanly designed set of GTC simulations, including a control where the fishbone Er is imposed on ITG and reproduces the suppression. That control is the strongest piece of evidence. The claim about being first is credible: previous work hypothesized this channel (refs 28, 29) but didn't couple the two in a global nonlinear run.\n\nWhat is new and good: the coupled n=1 fishbone + n≈25 ITG simulation with EPs; the demonstration that the fishbone saturates first and drives zonal flows larger than ITG's self-generated ones; the volume-averaged turbulence intensity traces; and the transport coefficients compared across three cases. The experimental comparison is correlational, but the ECE envelope matching the fishbone-only simulation and the Doppler reflectometry rotation during fishbone bursts are reasonable supporting evidence.\n\nWhere it's soft. The main issue is convergence. The paper states \"Numerical convergence studies are meticulously achieved\" but shows no scan. With Nψ=100 over r/a 0–0.63 (Δr/a ≈ 0.0063), the shearing rate ω_E×B requires a second radial derivative of δφ00, which amplifies grid-scale oscillations. The fishbone-driven Er explicitly has \"fine-scale radial structure.\" If that structure is partly marker noise or under-resolved, the computed shearing rate could be inflated, and the ω_E×B > γ_ITG threshold would not be robust. The imposed-Er control inherits this because it uses the same profile. This is fixable — a resolution scan, showing δφ00 and ω_E×B converge — but without it the quantitative claim isn't closed. The transport reduction itself (δφ from 1e-3 to 1e-4) is a large qualitative effect, so even if the shear rate is off by a factor of two, the mechanism likely survives; it's the threshold comparison that needs care.\n\nAlso, no error bars on χi or fluctuation amplitudes, and marker noise isn't discussed. The experimental validation is correlational (turbulence intensity drops during fishbone bursts) and doesn't directly measure shearing. Minor: a fair number of self-citations in the mechanism lineage, but they are the relevant prior work.\n\nConclusion: this is a serious paper and deserves a real referee. I'd send it to peer review; it's a genuine advance, and the main soft spot is an absent convergence scan, not a misfitted parameter. I'd probably cite it.","headline":"First coupled global gyrokinetic fishbone+ITG simulation with a credible zonal-flow suppression mechanism; the quantitative shearing-rate claim needs a convergence scan before fully buying it.","tokens_in":9582,"tokens_out":1804,"would_cite":true,"duration_ms":16681,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Fishbone-driven zonal radial electric fields suppress electromagnetic ITG turbulence and cut ion heat transport near neoclassical levels.","keywords":["fishbone instability","ITG turbulence","zonal flows","gyrokinetic simulation","multiscale interaction","ion thermal transport","internal transport barrier","tokamak plasmas"],"falsifier":"Run the coupled fishbone-plus-ITG simulation with doubled radial resolution and quadrupled marker count; if the roughly 100-fold drop in ITG amplitude and the χ_i reduction to about 0.2 m²/s weaken or disappear, the suppression is a numerical artifact rather than the claimed physical mechanism.","tokens_in":8727,"feed_emoji":"⚛️","tokens_out":5626,"duration_ms":51080,"temperature":0.7,"pith_summary":"The paper tries to establish that a macroscopic instability, the fishbone, can actively improve ion confinement in a tokamak by suppressing microscopic turbulence. Using global gyrokinetic simulations that for the first time treat the fishbone and ITG turbulence together, it finds the fishbone nonlinearly drives zonal radial electric fields whose shear rate exceeds the ITG growth rate, breaking up turbulent eddies. As a result, ion heat conductivity drops from roughly 1.5 m²/s to about 0.2 m²/s, near the neoclassical collisional floor. The paper argues this zonal-flow channel is the dominant cross-scale coupling mechanism and matches experimental observations of turbulence suppression during fishbone bursts. If correct, it identifies a favorable multiscale effect on confinement in fusion plasmas.","feed_headline":"Fishbone bursts suppress turbulence and slash ion heat loss","feed_subtitle":"A macroscopic instability creates zonal shears that quench microscopic eddies, potentially explaining improved core confinement.","key_machinery":"The central object is the zonal radial electric field, the axisymmetric (n=0) component of the perturbed electrostatic potential that the fishbone excites through a beat-driven process. This field produces an E×B shearing rate well above the ITG linear growth rate, and its fine radial structure decorrelates the ITG eddies. The paper proves the zonal flow is the active mediator by imposing the extracted fishbone-driven electric field on an ITG-only simulation and recovering nearly the same turbulence suppression as in the fully coupled run.","core_discovery":"The paper reports the first global gyrokinetic simulations that simultaneously resolve the macroscopic fishbone (toroidal mode n=1) and microscopic electromagnetic ITG turbulence (n≈25) in a tokamak. The fishbone saturates first because of its higher growth rate and drives zonal radial electric fields through a beat process; the resulting E×B shearing rate exceeds the ITG linear growth rate and suppresses the turbulence. The ion heat conductivity falls to about 0.2 m²/s, close to neoclassical values, while the fishbone itself is nearly unaffected because ITG fluctuations are much weaker. The simulation's electron-temperature perturbation envelope matches electron cyclotron emission measureme","pith_inferences":["If the shearing mechanism is general, other macroscopic modes that efficiently drive zonal flows could be used to actively control turbulent transport; the paper's logic does not require the mode to be a fishbone specifically.","The fine radial structure of the fishbone-driven zonal field suggests the suppression is localized; tailoring the safety-factor profile or energetic-particle drive could position that shear layer exactly where a transport barrier is desired.","A direct experimental test would measure the shearing rate of the fishbone-induced Doppler shift and the turbulence amplitude across a single burst; the suppression should track the zonal-flow amplitude, not the fishbone mode amplitude itself.","The multi-burst experimental trace hints at cumulative suppression; a dedicated experiment with repeated fishbone bursts could determine whether each burst deepens the barrier or the effect saturates."],"forward_implications":["Turbulence and fishbone cannot be treated independently: transport models that ignore the zonal-flow coupling will miss a dominant suppression channel.","Fishbone bursts in the plasma core can lower ion heat transport to near-neoclassical levels, offering a candidate mechanism for internal transport barriers observed after fishbone activity.","The interaction is largely one-directional: the fishbone regulates ITG turbulence, while the ITG leaves the fishbone's growth rate and saturation amplitude essentially unchanged.","In a burning plasma, alpha-particle-driven fishbone could similarly regulate core turbulence, but net confinement depends on competing energetic-particle transport caused by the fishbone.","The experimental Doppler-reflectometry signature—rotation and reduced fluctuation amplitude during bursts—can serve as a diagnostic marker for this cross-scale mechanism."],"fun_headline_variants":["Fishbone-driven zonal shears quench ITG turbulence","Macro fishbone suppresses micro turbulence, cuts heat loss","First global simulation: fishbone tames turbulence via zonal fields","Tokamak fishbone generates shears that slash ion heat transport","Cross-scale win: fishbone's zonal fields silence eddies"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The suppression claim rests on the coupled gyrokinetic simulation resolving both the fishbone and ITG turbulence with enough fidelity that the large fishbone-driven zonal flow does not quench the ITG eddies through numerical noise rather than physical shearing.","fun_headline_variants_meta":{"raw":{"variants":["Fishbone-driven zonal shears quench ITG turbulence","Macro fishbone suppresses micro turbulence, cuts heat loss","First global simulation: fishbone tames turbulence via zonal fields","Tokamak fishbone generates shears that slash ion heat transport","Cross-scale win: fishbone's zonal fields silence eddies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000222,"raw_usage":{"total_tokens":1233,"prompt_tokens":628,"completion_tokens":605,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":372,"completion_tokens_details":{"reasoning_tokens":517}},"tokens_in":372,"tokens_out":605,"duration_ms":6025,"temperature":1.0,"reasoning_tokens":517,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T23:46:29.137824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the coupled fishbone-plus-ITG simulation with doubled radial resolution and quadrupled marker count; if the roughly 100-fold drop in ITG amplitude and the χ_i reduction to about 0.2 m²/s weaken or disappear, the suppression is a numerical artifact rather than the claimed physical mechanism.","supporting_citations":[],"review_version":1}