{"id":"86382ea1-51ad-443f-891d-19090454c472","arxiv_id":"2412.10027","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Global electromagnetic gyrokinetic simulations show that an internal transport barrier forms only when the minimum safety factor sits at a low-order rational value, q=2, and that electron-driven zonal currents help lock in this condition.","lead":"This paper simulates plasma turbulence in fusion devices and identifies two ingredients needed to form an internal transport barrier: a safety factor minimum at a rational value (q=2) and electron dynamics that flatten the magnetic field profile. The result is a step toward predicting and controlling improved confinement regimes in tokamaks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Section 6.1's decisive qmin=2 vs 2.03 comparison is not fully simulated: the qmin=2.03 flux-driven run is stopped before quasi-steady state and the no-ITB conclusion is carried by an extrapolated core profile, so the claimed necessity of rational qmin remains undemonstrated.","rationale":"The reader's verdict identifies the same weakness, and the manuscript itself flags it in Section 6.1 ('has not reached quasi steady-state yet'). I do not see a stronger objection. The gradient-driven results (Figs. 3, 5) independently show a sharp transport reduction for qmin=2 versus 2.03, so the mechanism is plausible; the flux-driven qmin=2 run's power balance is converged; the qmin=2.01 run shows the q-profile being pulled toward 2.0, and Appendix A checks the df0/dt approximation. These elements support the overall qualitative picture. What is missing is only the completed qmin=2.03 flux-driven comparison, which is the key counterfactual for 'necessary.' Because the conclusion is conditional on an extrapolated final state, the appropriate verdict is unchanged (conditional acceptance), not rejection: the claimed mechanism is coherent and supported by several auxiliary simulations, but the central necessity claim should be verified with a converged run before being asserted as established.","tokens_in":24895,"tokens_out":4791,"duration_ms":49210,"concrete_test":"Continue the Section 6.1 qmin=2.03 flux-driven EM run to an operational quasi-steady state, defined as the radially integrated turbulent loss in Fig. 19 matching the injected power to within 10% over the core (s<0.4), and then compare the final ion temperature profile with the dashed 'future' curve in Fig. 17. If a computationally cheaper proxy is needed, restart from the end state with the source strength increased by a factor that preserves the target on-axis temperature, so the core relaxes on a shorter time scale. If the converged qmin=2.03 profile develops a local ion-gradient buildup near s≈0.4-0.6 or remains above the qmin=2 temperature profile, the claimed necessity of qmin=2 is weakened; if it falls to or below the dashed extrapolation and stays smooth, the central claim is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central assertion is that qmin close to a low-order rational value (2) is a necessary ingredient for ITB formation. The decisive evidence is the flux-driven comparison in Section 6.1, where qmin=2 forms an ion ITB and qmin=2.03 is said not to. The qmin=2.03 case had not reached quasi-steady state when stopped; the text states this explicitly, and the 'no ITB' conclusion is carried by a dashed black 'future' curve in Fig. 17, built from the power-balance imbalance in Fig. 19. The extrapolation is conservative and internally coherent: core turbulent losses exceed the source, so the core temperature is expected to drop. But the central necessity claim is supported by a predicted, not simulated, final state. A second embedded assumption is that the same heating source (taken from the qmin=2 gradient-driven calibration) is the right common comparator; if the qmin=2.03 run were instead run until its own power balance converged, a barrier or a different corrugation could in principle develop. Thus the sharp contrast is not yet directly demonstrated, only inferred from a non-converged trajectory. This is a correctness risk, not a contradiction: the remaining uncertainty can be settled by continuing the simulation or by a cheaper proxy that accelerates relaxation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports global electromagnetic gyrokinetic simulations with the ORB5 code, examining internal transport barrier (ITB) formation in reversed-shear tokamaks with the minimum safety factor qmin near 2. The paper compares three electron models (adiabatic, hybrid, and fully kinetic) and finds that kinetic electron dynamics is essential: it produces strong zonal flows, temperature profile corrugation, and zonal current sheets that flatten the local safety factor profile. In flux-driven simulations, the authors report that qmin = 2 leads to an ion-channel ITB, while qmin = 2.03 does not, and they attribute this difference to turbulent eddy self-interaction enabled by the rational qmin. Additional simulations show that the width of the q-flattening region scales between rho* and rho_i, that the input power scaling is close to gyroBohm for the two system sizes considered, and that an initial qmin = 2.01 profile evolves toward qmin = 2.0 under the action of turbulence-induced zonal currents.","tokens_in":25199,"tokens_out":3733,"duration_ms":42000,"significance":"If the main claims hold, this work provides a credible causal chain linking a rational minimum safety factor to ITB formation through kinetic-electron-driven zonal flows and zonal currents, extending earlier flux-tube results to global flux-driven simulations. The systematic comparison of adiabatic, hybrid, and fully kinetic electron models is valuable, and the Appendix A check of the df0/dt|_0=0 assumption strengthens confidence in the numerical setup. The q-flattening mechanism is qualitatively consistent with the companion flux-tube paper (Ref. [1]), and the reported system-size dependence is a useful step toward extrapolating ITB physics to reactor scale. However, the central claim that 'qmin close to a lowest order rational value is a necessary ingredient' is currently supported by a single non-converged comparison, so the significance is contingent on closing that gap.","major_comments":[{"comment":"The decisive comparison between qmin = 2 and qmin = 2.03 is not completed: the text explicitly states that the qmin = 2.03 case 'has not reached quasi steady-state yet', and the no-ITB conclusion is carried by the dashed black 'future' curve in Fig. 17, which is an extrapolation based on the power imbalance in Fig. 19. Because this comparison is the primary evidence for the abstract's claim that rational qmin is a 'necessary ingredient', the claim is not yet demonstrated. Please extend the qmin = 2.03 run to a genuine quasi-steady state, or provide a quantitative convergence check (e.g., time traces of the core temperature and power balance showing approach to a steady state), and report how the central conclusion would change if the extrapolation is inaccurate.","section":"Section 6.1, Figs. 17-19"},{"comment":"The heating source for the qmin = 2.03 run is taken from the qmin = 2 gradient-driven calibration. The extrapolated cooling of the core relies on the assumption that the turbulent losses exceeding the source in the s in [0.2, 0.3] region persist long enough to lower the core temperature. However, if the qmin = 2.03 equilibrium were instead simulated toward its own power balance with a different source profile or amplitude, a different outcome cannot be excluded a priori. Please justify the common source choice or test the sensitivity of the conclusion to the source strength and shape.","section":"Section 6.1, Fig. 19"},{"comment":"The term 'necessary ingredient' overstates the logical force of a single pair of simulations. The paper demonstrates that, in this specific setup, qmin = 2 is sufficient to produce an ITB and qmin = 2.03 is not (if the extrapolation holds); it does not establish that no other qmin value or alternative mechanism can produce an ITB. Please soften the claim to state that rational qmin is a necessary ingredient in the parameter regime studied here, or explicitly discuss the limited scope of the necessity claim.","section":"Abstract and Section 7"}],"minor_comments":[{"comment":"The sentence 'The modified q-profile shown in Figure 15 has been computed using equation 5' appears to be a citation error: Eq. (5) defines the effective heat diffusivity chi, while the q-modification is computed via Eq. (7)-(8). Please correct the cross-reference.","section":"Section 5.2, text near Fig. 15"},{"comment":"The conclusion of 'almost GyroBohm scaling' is based on only two values of rho* with different levels of power balance; please add an explicit caveat that the scaling is preliminary and may depend on the chosen comparison temperature and on the non-converged status of the qmin = 2.03 case.","section":"Section 6.2, Fig. 24-25"},{"comment":"The phrase 'solves the full-f Vlasov equation in spite of the delta-f splitting' is awkward and potentially confusing; it would be clearer to say that the code evolves the full distribution function while using a control-variate delta-f splitting for noise reduction.","section":"Section 2, numerical setup"},{"comment":"There are several typographical errors and inconsistent hyphenation, for example 'Ècole Polytechnique Féedérale', 'V olčokas', and the varying use of 'quasi-steady state' versus 'quasi steady-state'. A careful proofreading pass is recommended.","section":"Throughout"},{"comment":"The dashed black curve is described as an extrapolation, but the method used to construct it (beyond 'smoothing of R/LT' and the power imbalance) is not fully specified; please describe the extrapolation procedure or add a note on its uncertainty.","section":"Figure 17"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically substantive and the numerical machinery appears appropriate for the questions addressed. The main concern is that the central necessity claim rests on a non-converged simulation, which the authors themselves acknowledge. The self-citation to Refs. [1,16,17] is justified because those works provide the flux-tube foundation for the eddy self-interaction concept. I would be willing to accept a revised version that either completes the qmin = 2.03 simulation or explicitly re-scopes the claim to the parameter regime simulated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a solid, useful paper that does something new—global electromagnetic flux-driven gyrokinetic simulations showing an ion ITB forming when qmin sits at the low-order rational value 2, with the q-profile flattened by turbulence-generated zonal currents. The gradient-driven sections are convincing: the factor-of-two transport reduction at qmin=2 versus 2.03 is shown clearly for both hybrid and fully kinetic electron models, and the parallel eddy self-interaction picture carries over from the group's flux-tube work. The qmin=2.01 case is a nice additional piece showing partial self-interaction evolving into full self-interaction via the q-flattening feedback. That part is arguably the strongest evidence for the mechanism because it is a simulated transition, not just a static contrast.\n\nWhere the paper goes soft: the decisive flux-driven contrast in Section 6.1. The qmin=2.03 run was stopped before reaching quasi-steady state, and the conclusion that no ITB forms rests on an extrapolated core temperature profile. The paper is honest about this—it says so in the text—but the abstract's \"necessary ingredient\" claim is stronger than what the simulation actually demonstrates. A longer run of the 2.03 case, or a cheaper proxy that accelerates relaxation, would settle it. Without that, the necessity claim is inference from a non-converged trajectory, not a direct demonstration. Also, single runs without uncertainty estimates mean the sharp sensitivity could be partly numerical; and the rho* scaling uses two points only, so \"almost GyroBohm\" is a hint, not a solid scaling law.\n\nNone of this sinks the paper. The mechanism is plausible, the gradient-driven evidence is solid, and the qmin=2 flux-driven barrier is clearly visible and close to power balance. The flaw is localized and fixable. I would send this to a serious referee, with the request that the authors either run the qmin=2.03 case to quasi-steady state or qualify the necessity claim as an inference. A good referee's time would be well spent on this.\n\nBest,\n\n[Your name]","headline":"A serious global gyrokinetic study with a genuinely new flux-driven ITB result, but the headline 'necessary ingredient' claim leans on an unconverged qmin=2.03 run and needs a longer run or a softer claim.","tokens_in":25770,"tokens_out":2083,"would_cite":true,"duration_ms":22268,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["52.35.Ra","52.55.Fa","52.65.Tt"],"model":"deepseek-v4-flash","headline":"Reversed-shear tokamaks form internal transport barriers only when the safety-factor minimum sits close to a low-order rational value (q=2), because that lets turbulent eddies close on themselves and sustain a strong sheared zonal flow.","keywords":["internal transport barrier","reversed magnetic shear","gyrokinetic simulation","zonal flows","zonal currents","eddy self-interaction","safety factor","turbulence suppression"],"falsifier":"Run the $q_{\\min}=2.03$ flux-driven simulation long enough that the radially integrated turbulent losses equal the injected power everywhere in the core, then compare the converged ion temperature and temperature-gradient profiles with the $q_{\\min}=2$ run; if the core gradient equilibrates at a similar level, the necessary-ingredient claim fails.","tokens_in":24726,"feed_emoji":"🌀","tokens_out":10500,"duration_ms":98592,"temperature":0.7,"pith_summary":"This paper tries to establish a causal condition for internal transport barrier (ITB) formation in reversed-shear tokamaks: the minimum of the safety factor profile, $q_{\\min}$, must lie on or very near a low-order rational value (2 in these simulations), so that turbulent eddies can wrap around the torus and 'bite their own tail.' Using global electromagnetic gyrokinetic simulations with kinetic electrons, the authors show that including electron dynamics is essential even when the dominant instability is the ion temperature gradient mode; kinetic electrons produce strong zonal-flow shearing that tears eddies apart and zonal current sheets that flatten and lower the q profile around the rational surface. In flux-driven runs with $q_{\\min}=2$ an ion-channel transport barrier develops at two radii around the q minimum, while with $q_{\\min}=2.03$ no barrier forms. The result matters because it points to a concrete testable rule for triggering ITBs in fusion devices and explains why experiments preferentially see barriers when $q_{\\min}$ hits an integer.","feed_headline":"Rational qmin is a necessary ingredient for tokamak ITBs","feed_subtitle":"Flux-driven simulations form an ion barrier at qmin=2, none at 2.03, because eddies bite their own tail.","key_machinery":"The load-bearing mechanism is parallel eddy self-interaction at a low-order rational surface. In a low-shear region, turbulent eddies become extremely elongated along the magnetic field; when $q_{\\min}$ is exactly 2, they extend far enough to close on themselves, 'biting their own tail' and locking turbulence to that surface. This self-interaction, enabled only by kinetic electron dynamics, produces two coupled structures: a strong zonal $E\\times B$ shearing dipole that de-correlates eddies and quenches transport, and a steady zonal parallel electron current that, through the perturbed vector potential, flattens and locally shifts the safety factor. The flattening creates a positive feedback loop when $q_{\\min}$ is close enough to the rational value (within roughly $\\Delta q \\sim \\rho^* q_0/n$), dragging $q_{\\min}$ toward 2 and strengthening self-interaction.","core_discovery":"On the paper's own terms, the central discovery is that rational $q_{\\min}$ is a necessary ingredient for ITB formation in reversed-shear configurations. With adiabatic electrons, the transport coefficient shows only a smooth monotonic increase with $q_{\\min}$ and no special response at $q_{\\min}=2$; once kinetic (drift-kinetic) electrons are included, $q_{\\min}=2$ produces a global reduction of heat transport, strong radial profile corrugation, and a steady dipole of $E\\times B$ shearing around the zero-shear surface, whereas $q_{\\min}=2.03$ behaves like a non-rational value. The electromagnetic response adds zonal parallel currents that persist in time and, through Ampère's law, modify the safety factor: the q profile flattens around $q_{\\min}$, an effect that is weak when self-interaction is already complete but decisive when it is only partial. In flux-driven simulations the $q_{\\min}=2$ case forms an internal transport barrier in the ion channel at inner and outer radii around $q_{\\min}$; the $q_{\\min}=2.03$ case, run to a state where core losses still exceed input power, is extrapolated to lose core temperature and not form a barrier. A case starting at $q_{\\min}=2.01$ shows partial eddy self-interaction evolving into complete self-interaction as zonal currents drag q down to 2.0, in agreement with flux-tube predictions. The width of the flattened q region scales between $\\rho^*$ and $\\rho_i$, and the power needed for equal on-axis temperatures is close to gyro-Bohm in the two system sizes considered.","pith_inferences":["This suggests a possible control strategy: modest current drive that nudges $q_{\\min}$ onto an integer could trigger an ITB even without a strong shear reversal, and the resulting barrier might then be self-sustaining via the q-flattening feedback.","Because $q_{\\min}=2.01$ evolves toward 2 while $q_{\\min}=2.03$ does not, the outcome may be history-dependent for a narrow range of $q_{\\min}$, with the same nominal target producing a barrier or not depending on the path taken.","If the $\\rho^*$-dependent tolerance is generic, larger devices (smaller $\\rho^*$) would need tighter control of $q_{\\min}$ to exploit rational-surface triggering, which could be tested by scanning $q_{\\min}$ offsets across device sizes.","A testable extension would be to map the threshold curve $\\Delta q_{\\max}(\\rho^*)$ in global simulations and see whether it quantitatively matches the flux-tube estimate $\\Delta q \\sim \\rho^* q_0/n$."],"forward_implications":["If the claim is right, placing $q_{\\min}$ on a low-order rational value is a necessary condition for an ion-channel ITB in reversed-shear tokamaks, not merely a helpful detail.","Kinetic electron response must be retained in simulations of ITB onset: adiabatic-electron models miss the $q_{\\min}=2$ sensitivity entirely in this scenario.","Turbulent self-generated currents can move $q_{\\min}$ downward toward a rational value, so the q profile and the turbulence state are self-consistently coupled: a plasma starting at $q_{\\min}=2.01$ can transition to full self-interaction and form a barrier.","The sensitivity window for self-interaction widens with $\\rho^*$, so smaller devices (larger $\\rho^*$) should tolerate larger deviations of $q_{\\min}$ from the rational value.","Power requirements for equal on-axis temperatures in this barrier regime are close to gyro-Bohm: going from the smaller to the larger device size costs only about 20% more power in the two cases considered."],"supporting_citations":[{"why":"Establishes the flux-tube result that turbulence-generated parallel currents flatten the safety factor at rational surfaces in low-shear tokamaks, the effect that the global runs here reproduce.","marker":"[1]"},{"why":"Shows that ultra-long turbulent eddies and magnetic topology control internal transport barrier triggering, providing the self-interaction picture used throughout.","marker":"[16]"},{"why":"Quantifies how eddy self-interaction depends on how far $q_{\\min}$ sits above the rational value, including the estimate $\\Delta q \\sim \\rho^* q_0/n$ used to interpret $q_{\\min}=2.01$.","marker":"[17]"},{"why":"Reports gyrokinetic simulations with kinetic electrons showing large off-axis minimum-q profile corrugations and global transport reduction at $q_{\\min}=2$, the earlier result this work extends to electromagnetic flux-driven runs.","marker":"[13]"},{"why":"Shows with adiabatic electrons that transport stays smooth across a minimum-q surface, explaining why adiabatic models miss the rational-$q_{\\min}$ effect.","marker":"[12]"},{"why":"Connects non-adiabatic passing-electron dynamics to fine radial structures near low-order rational surfaces, supporting the electron-dynamics mechanism.","marker":"[14]"},{"why":"Argues that parallel eddy length is set by critical balance in ITG turbulence, underpinning why eddies can become long enough to self-interact near rational q.","marker":"[34]"}],"fun_headline_variants":["Rational qmin: the on-switch for tokamak ITBs","Tokamak ITBs require qmin=2; electrons set the stage","Eddy self-interaction at qmin=2 keys ITB formation","Zonal currents at rational qmin build transport barriers","qmin=2 plus electron dynamics flips on ITBs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The $q_{\\min}=2.03$ flux-driven case was stopped before reaching quasi-steady state, so the claim that no ITB forms at $q_{\\min}=2.03$ depends on the assumption that the core temperature continues to drop as its yet-unbalanced power losses suggest.","fun_headline_variants_meta":{"raw":{"variants":["Rational qmin: the on-switch for tokamak ITBs","Tokamak ITBs require qmin=2; electrons set the stage","Eddy self-interaction at qmin=2 keys ITB formation","Zonal currents at rational qmin build transport barriers","qmin=2 plus electron dynamics flips on ITBs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00034,"raw_usage":{"total_tokens":2047,"prompt_tokens":1292,"completion_tokens":755,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":908,"completion_tokens_details":{"reasoning_tokens":664}},"tokens_in":908,"tokens_out":755,"duration_ms":8611,"temperature":1.0,"reasoning_tokens":664,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:25:39.765378+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the $q_{\\min}=2.03$ flux-driven simulation long enough that the radially integrated turbulent losses equal the injected power everywhere in the core, then compare the converged ion temperature and temperature-gradient profiles with the $q_{\\min}=2$ run; if the core gradient equilibrates at a similar level, the necessary-ingredient claim fails.","supporting_citations":[{"cited_title":"Since the first discovery in JET [ 2], a safety factor profile with a reversed shear region has been used to generate Internal Transport Bar riers (ITB) in several devices","cited_arxiv_id":null,"evidence_quote":"Establishes the flux-tube result that turbulence-generated parallel currents flatten the safety factor at rational surfaces in low-shear tokamaks, the effect that the global runs here reproduce."},{"cited_title":"Hugon, B.Ph","cited_arxiv_id":null,"evidence_quote":"Shows that ultra-long turbulent eddies and magnetic topology control internal transport barrier triggering, providing the self-interaction picture used throughout."},{"cited_title":"Joffrin, C.D","cited_arxiv_id":null,"evidence_quote":"Quantifies how eddy self-interaction depends on how far $q_{\\min}$ sits above the rational value, including the estimate $\\Delta q \\sim \\rho^* q_0/n$ used to interpret $q_{\\min}=2.01$."},{"cited_title":"Our simulations setup is inspired by Cyclone Base Case param eters but excludes density gradients","cited_arxiv_id":null,"evidence_quote":"Reports gyrokinetic simulations with kinetic electrons showing large off-axis minimum-q profile corrugations and global transport reduction at $q_{\\min}=2$, the earlier result this work extends to electromagnetic flux-driven runs."},{"cited_title":"Candy, R","cited_arxiv_id":null,"evidence_quote":"Shows with adiabatic electrons that transport stays smooth across a minimum-q surface, explaining why adiabatic models miss the rational-$q_{\\min}$ effect."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Connects non-adiabatic passing-electron dynamics to fine radial structures near low-order rational surfaces, supporting the electron-dynamics mechanism."},{"cited_title":"McMillan, S","cited_arxiv_id":null,"evidence_quote":"Argues that parallel eddy length is set by critical balance in ITG turbulence, underpinning why eddies can become long enough to self-interact near rational q."}],"review_version":1}