{"id":"41bf3d8b-ccd5-4492-92bb-2e30756a1595","arxiv_id":"2607.11784","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Equilibrium E×B flow shear can break the Dimits-state zonal flows of ITG turbulence, producing a non-monotonic transport response with a sharp heat-flux increase between weak and strong shear.","lead":"Gyrokinetic simulations show that a sheared plasma flow can destroy the self-organized zonal flows that normally suppress turbulence, causing heat transport to jump sharply instead of falling. The result suggests spherical-tokamak rotation may be set by heat input rather than by momentum injection alone, contradicting the standard picture that flow shear always improves confinement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Geometric incompatibility argument hinges on unproven fixed finite ℓmax; if ℓmax varies with ω⊥ or box size, the predicted instability interval may be an artifact.","rationale":"The paper presents a strong simulation evidence base: two independent models (GENE gyrokinetics and a reduced fluid model), PVG term removal, a partial global check, and multiple MAST-U discharges. The non-monotonic transport response is a robust numerical observation. The geometric argument is an elegant formalization, but its key premise — a finite, fixed, physically determined ℓmax — is explicitly acknowledged as unexplained in Section III. This is precisely the condition that makes λ < 1 and creates the claimed instability interval. If ℓmax is instead set by the simulation box or by the imposed shear itself, the formal argument loses its predictive content, even if the phenomenon persists. The reader's weakest_assumption already identifies this same issue, so the conditional verdict remains appropriate. No change in verdict is needed; the concern reinforces the need for the supplementary estimates to be independently scrutinized and for a dedicated box-size/ω⊥-dependence test to be reported.","tokens_in":15572,"tokens_out":5972,"duration_ms":55486,"concrete_test":"In the idealized Miller geometry at R0/LTi=15 (below Dimits threshold), run local GENE simulations at ω⊥=0 and at ω⊥≈0.5ωc with radial box sizes Lx=50, 100, 200, 400 ρs. From the time-averaged profiles of total shear (e.g., Fig. 2e–h), measure the maximum width of contiguous regions where the total shear has one sign and magnitude near ωc (ℓmax) and the minimum width (ℓmin). Test two things: (1) Does ℓmax at ω⊥=0 saturate with Lx, or does it scale with Lx? (2) Does ℓmax at ω⊥≈0.5ωc differ by more than the run-to-run variability from its ω⊥=0 value? If ℓmax scales with Lx or changes with ω⊥, Eq. (4)'s fixed-λ interval is not supported; if ℓmax saturates and is ω⊥-independent, compute λ from these widths and compare with the observed onset of the transport rise (ω⊥_onset/ωc from Fig. 1) to validate the quantitative prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formal claim is the compatibility condition ω⊥/ωc ≤ λ = (ℓmax−ℓmin)/(ℓmax+ℓmin) (Eq. 4), which defines the predicted instability interval. The derivation of Eq. (4) is algebraically correct given the piecewise-constant ±ωc profile and the inequalities ℓ+ ≤ ℓmax and ℓ− ≥ ℓmin. However, the argument's predictive power rests entirely on ℓmax and ℓmin being fixed, finite, and independent of the imposed shear ω⊥ and of the radial domain size. Section III explicitly states 'the mechanism setting ℓmax remains to be established' and offers only a speculative 'finite radial mean free path of avalanches/ferdinons'. If ℓmax is not a genuine physical scale — e.g., if it grows with the radial box size in local flux-tube simulations (where radial boundary conditions enforce periodic zonal patterns) or if it depends on ω⊥ because the turbulence itself changes with total shear — then λ is not a fixed threshold, and the claimed interval λ < ω⊥/ωc < 1 may be an artifact of the chosen simulation domain rather than a generic instability. The observed finite zonal scale in the ω⊥=0 saturated state is reported as independent of box size, but this is only a statement for ω⊥=0 and is not documented in the main text (Supplemental [50]). This is the weakest link in the otherwise well-supported central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a non-monotonic dependence of ion heat transport on imposed equilibrium E×B flow shear in local gyrokinetic ITG simulations: weak shear leaves the Dimits state roughly intact, intermediate shear destroys the self-generated zonal flows and produces a sharp transport increase, and strong shear quenches turbulence. The authors interpret this as a geometric incompatibility: if the Dimits state requires alternating zonal-shear layers of widths between ℓmin and ℓmax, an imposed shear ω⊥ is compatible only for ω⊥/ωc ≤ λ=(ℓmax−ℓmin)/(ℓmax+ℓmin), leaving a destabilizing interval λ<ω⊥/ωc<1. A reduced 2D fluid model reproduces the phenomenology. Gyrokinetic simulations of MAST-U discharges find the inferred rotation shear at or below the onset of the transport increase, suggesting rotation can be limited by heat injection.","tokens_in":15990,"tokens_out":5885,"duration_ms":60136,"significance":"If the mechanism holds, it overturns the standard expectation that equilibrium flow shear monotonically suppresses ITG transport and identifies a concrete constraint on spherical-tokamak operation, with possible broader implications for shear-driven zonal/mean-flow interactions. The paper's strengths are its systematic numerical evidence: the GENE scans cover multiple gradients and explicitly remove PVG to rule out the obvious alternative mechanism; the fluid model shows the effect is not tied to kinetic or toroidal details; the MAST-U analysis covers six discharges and includes kinetic electrons, electromagnetic fluctuations, and collisions; and a global GENE run is cited as a partial check. The main weakness is that the central geometric criterion, Eq. (4), is not yet a predictive theory because ℓmin and ℓmax are estimated from the same simulations and the scale ℓmax is admittedly not understood.","major_comments":[{"comment":"The compatibility condition is derived from assumptions (i)–(iv), but its predictive content depends on ℓmin and ℓmax being fixed, finite, and independent of ω⊥ and of the simulation domain. The manuscript states that 'the mechanism setting ℓmax remains to be established' and offers a speculative finite-mean-free-path argument; the quoted estimates are made from the ω⊥=0 saturated Dimits state in the Supplemental Material. Thus λ is calibrated on the same simulations Eq. (4) is supposed to explain. If ℓmax changes with ω⊥ — e.g., because the avalanche/ferdinon mean free path depends on total shear — or if it grows with the radial box size in the flux-tube geometry, the predicted interval is not a robust prediction. Please provide direct measurements of ℓ±(ω⊥) and ℓmax(Lx), an independent estimate of ℓmax, and a comparison of Eq. (4) with the measured onset in Fig. 1 without free adjustme","section":"Section III, Eqs. (2)–(4)"},{"comment":"The derivation idealizes the total shear profile as piecewise constant equal to ±ωc and counts N discrete shear regions. The actual profiles in Fig. 2(e–h) are smooth and asymmetric, with no clear square-wave structure. It is not demonstrated that the inequalities ℓ+≤ℓmax and ℓ−≥ℓmin apply to the actual profiles or that N is well-defined. Please show that the square-wave idealization is conservative, for example by checking Eq. (3) against the measured widths in the GK and fluid runs, or derive the criterion for continuous profiles.","section":"Section III, assumption (i) and Fig. 4"},{"comment":"The simulations establishing the effect (Fig. 1 and the fluid model) are local, gradient-driven, and radially periodic; the global gyrokinetic check mentioned in footnote [75] is not documented. Because the geometric argument invokes a finite ℓmax and the local flux-tube domain imposes radial periodicity on zonal flows, the possibility that ℓmax is influenced by the periodic box (or that profile relaxation changes the zonal response) needs to be addressed with quantitative evidence. Please report the global-run setup and result, or provide a dedicated finite-domain convergence study of the onset shear.","section":"Section II and footnote [75]"}],"minor_comments":[{"comment":"The panel labeling is inconsistent: (a,c) are time traces while (b),(d) are flux-versus-shear plots, but the caption reads 'Time traces ... (a,c), and (b) ion heat flux and (d) ion toroidal angular momentum flux ...'. Please clarify, e.g., 'Panels (a) and (c): time traces; panels (b) and (d): fluxes versus flow shear.'","section":"Figure 5 caption"},{"comment":"The symbol λ is introduced but not given a name; later it is referred to as a threshold. Consider defining it explicitly as the 'compatibility threshold' to avoid confusion with the plasma micro-scales.","section":"Section III, Eq. (4)"},{"comment":"The estimates of ℓmin and ℓmax, the box-size-independence check, and the derivation of the mean-flow-shear terms in the fluid model are all relegated to the Supplement. Since the finite-ℓmax assumption underpins the central theoretical claim, the box-size-independence result should appear in the main text or at least be described with a quantitative summary.","section":"Supplemental Material [50]"},{"comment":"The remark that y-axes are '(same for the pairs of simulations with equal radial box size)' is unclear. Please specify the radial box size in each panel or state explicitly where the normalization changes.","section":"Figure 2 caption"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper deserves a serious referee, and the central result — that imposed E×B shear can destroy the Dimits-state zonal flows and sharply increase transport — is credible. The non-monotonic transport response is demonstrated in gyrokinetic simulations, reproduced in a reduced fluid model, and supported by the MAST-U comparisons. The PVG-removal checks and the bistability check are well done and rule out the most obvious alternative explanations.\n\nThe geometric incompatibility argument (Section III) is an appealing way to frame the mechanism, but the predictive threshold λ depends on ℓmin and ℓmax, and ℓmax's mechanism is explicitly left open. Those constants are inferred from the same simulations the argument is meant to explain, so the predicted instability interval is partially calibrated rather than independently derived. The stress-test worry that ℓmax may vary with ω⊥ or box size is a fair concern, but not fatal: the paper reports in the supplemental that the zero-shear zonal scale is box-size independent, which gives the argument some footing. Still, the threshold λ is only as good as the assumption that ℓmax is fixed. That should be a main referee question, along with an error bar on λ.\n\nThe idealized scans use adiabatic electrons; the authors acknowledge this and provide kinetic-electron MAST-U runs, so it is a limitation but not a hidden one. The experimental section relies on inferred rotation shear with modest statistics; the \"at or just below threshold\" claim is suggestive rather than decisive.\n\nThe paper is honest about what it does not know — it states plainly that \"the mechanism setting ℓmax remains to be established\" — and the main physics (zonal-flow breakdown) is well supported by the simulation phenomenology. The fluid model is a nice touch, showing the effect is not an artifact of toroidal geometry or kinetic details.\n\nThis is a paper for plasma transport theorists and tokamak rotation specialists, and it deserves referee time. I would send it out and ask for a dedicated test of ℓmax's dependence on ω⊥ and box size, plus a sensitivity scan of λ to the estimates of ℓmin/ℓmax.","headline":"Imposed E×B shear can break the Dimits state and sharply increase ITG transport; the simulation evidence is solid, though the geometric threshold theory leans on constants fitted from the same runs.","tokens_in":16430,"tokens_out":1367,"would_cite":true,"duration_ms":15078,"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":"Imposed E×B flow shear can destroy the zonal flows that regulate plasma turbulence, setting off a sharp rise in heat transport.","keywords":["ion-temperature-gradient turbulence","zonal flows","E×B flow shear","Dimits regime","gyrokinetic simulations","turbulent transport","spherical tokamaks","flow-shear destabilization"],"falsifier":"Run a flux-tube gyrokinetic or fluid simulation below the Dimits threshold and scan the imposed E×B shear finely across 0 < ω⊥/ωc < 1; if the time-averaged heat flux never rises above its zero-shear value before being quenched, the claimed destabilization window does not exist. Equivalently, a controlled tokamak rotation-shear scan with fixed profiles that shows no local heat-flux peak would contradict the mechanism.","tokens_in":15502,"feed_emoji":"🌀","tokens_out":6234,"duration_ms":59518,"temperature":0.7,"pith_summary":"This paper establishes that in the low-transport state known as the Dimits regime, an externally imposed E×B flow shear does not simply suppress turbulence: once the imposed shear approaches the strength of the self-generated zonal flows, it can destroy them, causing the turbulent heat flux to rise sharply. The mechanism is geometric: the Dimits state needs alternating bands of flow shear with widths in a fixed range, and an imposed shear cannot be accommodated once it exceeds a threshold set by the ratio of the minimum and maximum allowed band widths. A reduced fluid model reproduces the same behavior, showing the effect is not an artifact of kinetic or toroidal details. Simulations of spherical tokamak discharges find the experimentally inferred rotation shear sitting at or just below the threshold, suggesting that toroidal rotation in such devices is limited primarily by heat injection rather than by momentum injection alone. This overturns the textbook expectation that equilibrium flow shear monotonically improves confinement.","feed_headline":"Flow shear can destroy the zonal flows that tame plasma turbulence","feed_subtitle":"In the Dimits regime, a window of imposed shear sharply raises heat transport — and tokamak rotation sits near its edge.","key_machinery":"The load-bearing object is a simple inequality, the compatibility condition ω⊥/ωc ≤ λ ≡ (ℓmax−ℓmin)/(ℓmax+ℓmin), where ω⊥ is the imposed E×B shear, ωc the critical shear of the Dimits state, and ℓmin and ℓmax the minimum and maximum radial widths that a zonal-shear band can have while still suppressing turbulence. The argument counts how many alternating shear regions of total shear ±ωc can fit across a domain of width L; when the imposed shear is too large, no integer number of bands satisfies both width bounds, so the self-organized pattern breaks down. This geometric argument carries the paper's explanation for why transport increases in the interval λ < ω⊥/ωc < 1. A second key ingredient","core_discovery":"The central claim is that imposed equilibrium E×B flow shear can destabilize the Dimits state of ion-temperature-gradient turbulence. In this low-transport state, self-organized zonal flows—large-scale bands of perpendicular flow—regulate the turbulence. The paper shows that the turbulent eddies respond only to the total perpendicular shear, the sum of imposed and zonal shear. Weak imposed shear is absorbed by a reorganization of the zonal-flow pattern, but when the imposed shear becomes comparable to the intrinsic zonal shear, the alternating zonal-shear regions can no longer satisfy the required width bounds, the zonal flows break down, and heat transport rises sharply before being quenche","pith_inferences":["If the compatibility condition is generic, existing databases of flow-shear scans in other tokamaks could be re-examined for a heat-flux peak below the quench threshold; a local maximum in transport versus rotation shear would be a direct experimental signature.","The same geometric incompatibility might appear in planetary atmospheres or oceans, where externally forced mean zonal winds interact with self-organized zonal jets; idealized beta-plane simulations with imposed large-scale shear could test whether the zonal jet pattern breaks down analogously.","The empirical finding that the Dimits shift grows at low safety factor and tight aspect ratio implies the destabilization window widens in such geometries, which may make the effect more prominent in future compact tokamaks or spherical devices."],"forward_implications":["In the Dimits regime, equilibrium flow shear is not a monotonic confinement knob: there is a window of imposed shear in which heat transport rises sharply before larger shear quenches the turbulence.","The mechanism is independent of kinetic effects and toroidal geometry, since a minimal two-dimensional fluid model reproduces it; analogous combinations of self-organized and imposed shear in other turbulent systems may show the same breakdown.","In spherical tokamaks, steady-state rotation shear can be pinned below the destabilization threshold, so toroidal rotation is set by the balance of heat and momentum injection rather than by momentum diffusivity alone.","Machines operating in this regime face a heat-flux hill: to reach the strongly suppressed high-shear state, enough heat must be injected to sustain the profiles through the enhanced transport."],"fun_headline_variants":["Shear kills zonal flows, spiking heat transport","Imposed flow shear can shatter turbulence-regulating zonal flows","When shear overpowers zonal flows, heat transport rises sharply","How equilibrium shear can destroy the Dimits state","Flow shear breaks zonal-flow barrier, ramps up transport"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The geometric argument assumes that the Dimits state always requires alternating shear bands whose widths stay between a fixed minimum and a fixed maximum that do not depend on the imposed shear; the paper notes the mechanism setting the maximum width is not yet established.","fun_headline_variants_meta":{"raw":{"variants":["Shear kills zonal flows, spiking heat transport","Imposed flow shear can shatter turbulence-regulating zonal flows","When shear overpowers zonal flows, heat transport rises sharply","How equilibrium shear can destroy the Dimits state","Flow shear breaks zonal-flow barrier, ramps up transport"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000174,"raw_usage":{"total_tokens":1067,"prompt_tokens":641,"completion_tokens":426,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":385,"completion_tokens_details":{"reasoning_tokens":343}},"tokens_in":385,"tokens_out":426,"duration_ms":5071,"temperature":1.0,"reasoning_tokens":343,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T06:44:34.148629+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a flux-tube gyrokinetic or fluid simulation below the Dimits threshold and scan the imposed E×B shear finely across 0 < ω⊥/ωc < 1; if the time-averaged heat flux never rises above its zero-shear value before being quenched, the claimed destabilization window does not exist. Equivalently, a controlled tokamak rotation-shear scan with fixed profiles that shows no local heat-flux peak would contradict the mechanism.","supporting_citations":[],"review_version":2}