Non-Ambipolarity of Microturbulent Transport

Allen H Boozer

arxiv: 2601.15661 · v2 · submitted 2026-01-22 · ⚛️ physics.plasm-ph

Non-Ambipolarity of Microturbulent Transport

Allen H Boozer This is my paper

Pith reviewed 2026-05-16 12:34 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph

keywords microturbulencenon-ambipolar transportchaotic magnetic fieldsplasma betagyro-Bohm diffusionquasi-neutralitytokamak transport

0 comments

The pith

Turbulent magnetic perturbations at finite beta create chaotic field lines whose quasi-neutrality condition produces electron transport that balances non-ambipolar ion diffusion.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that exact magnetic surfaces are absent once finite-pressure electrostatic microturbulence is considered, because the accompanying turbulent magnetic field makes field lines chaotic. Quasi-neutrality along those lines forces the turbulent potential to satisfy en B·∇Φ = B·∇p_e, which drives radial electron transport whose effective diffusion coefficient is set by the correlation length Δ along the chaotic lines. This electron motion supplies a viscosity on the electron flow that can offset a non-ambipolar fraction f_na of the ion gyro-Bohm particle flux, thereby avoiding the large radial electric field that would otherwise be needed to confine ions and impurities. The maximum counterbalanceable f_na and the minimum plasma beta required to keep the magnetic perturbations unshielded are derived explicitly.

Core claim

The chaos-produced electron transport gives an effective viscosity on the electron flow, which can counterbalance a non-ambipolar part of the ion radial particle diffusion that is f_na times gyro-Bohm diffusion. The maximum f_na that can be counterbalanced and the required plasma beta to avoid shielding the magnetic perturbations B̃ are calculated.

What carries the argument

The quasi-neutrality relation en B·∇Φ = B·∇p_e imposed along chaotic magnetic field lines, which fixes the turbulent potential and thereby the electron radial transport that acts as viscosity.

If this is right

Non-ambipolar ion diffusion up to the calculated maximum f_na can be sustained without a radial electric field that would confine impurities.
A minimum plasma beta must be reached before the magnetic perturbations remain unshielded and the balancing viscosity appears.
The effective electron transport coefficient scales as D_ef = (Δ/a_T) T_e/eB, directly tied to the temperature gradient length a_T.
Impurity accumulation that would otherwise require strong radial E can be avoided when the electron viscosity is active.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same mechanism may set an upper limit on how far turbulence can depart from ambipolarity before magnetic chaos restores balance.
Varying beta in experiments while holding turbulence drive fixed would test whether the non-ambipolar fraction saturates at the predicted value.
The correlation length Δ along chaotic lines becomes a new control parameter for transport modeling once the beta threshold is crossed.

Load-bearing premise

That finite-pressure turbulence produces an unshielded magnetic perturbation whose chaotic field lines enforce exactly the stated quasi-neutrality relation without additional shielding or kinetic corrections.

What would settle it

A direct measurement showing that turbulent magnetic perturbations are fully shielded at the beta values required by the calculation, or that the observed non-ambipolar transport fraction exceeds the predicted maximum f_na, would falsify the mechanism.

Figures

Figures reproduced from arXiv: 2601.15661 by Allen H Boozer.

read the original abstract

When restricted to magnetic flux tubes, the gyrokinetic theory of microturbulence gives the same radial transport for ions and electrons. But, exact magnetic surfaces do not exist in the presence of what is called electrostatic microturbulence. At a finite plasma pressure, a turbulent electric potential is accompanied by a turbulent magnetic field $\tilde{B}$, which makes the magnetic field lines chaotic. Quasi-neutrality along the chaotic magnetic field lines requires a potential that obeys $en \vec{B}\cdot \vec{\nabla} \Phi = \vec{B}\cdot \vec{\nabla} p_e$, where $p_e$ is the electron pressure. This potential produces radial transport similar to that of diffusion coefficient $D_{ef}= (\Delta/a_T)T_e/eB$. $\Delta$ is the radial distance over which the potential $\Phi$ is correlated by the electron motion along the chaotic magnetic field, and $|dT_e/dr| = T_e/a_T$. The chaos-produced electron transport gives an effective viscosity on the electron flow, which can counterbalance a non-ambipolar part of the ion radial particle diffusion that is $f_{na}$ times gyro-Bohm diffusion. This non-ambipolarity would otherwise require a radial electric field that confines ions and hence impurities. The maximum $f_{na}$ that can be counterbalanced and the required plasma beta to avoid shielding the magnetic perturbations $\tilde{B}$ are calculated.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Boozer's note gives a concrete way for chaotic field lines at finite beta to let electron transport offset non-ambipolar ion flux, but it rests on an unshielded turbulent B that may not hold.

read the letter

The central claim is that finite-beta microturbulence produces chaotic field lines, and quasi-neutrality along them forces a potential relation en B·∇Φ = B·∇p_e. This yields an effective electron diffusion D_ef = (Δ/a_T) T_e/eB that acts like a viscosity on the electron flow, balancing a fraction f_na of the gyro-Bohm ion diffusion without needing a confining radial electric field. The paper also estimates the largest f_na that can be offset and the beta threshold to keep the turbulent B unshielded. That mechanism is the new piece; standard flux-tube gyrokinetics assumes intact surfaces and equal ion-electron transport, so this offers a route to non-ambipolarity that is not in the usual citations. The derivation is internally consistent once the chaos and the exact parallel balance are granted, and it directly addresses impurity control in confinement devices. The soft spots are the two load-bearing assumptions. First, the turbulent B must remain unshielded enough to ergodize lines over the correlation length Δ; kinetic or diamagnetic shielding could easily reduce it below the calculated threshold. Second, the potential relation drops inertia, resistivity, and ion polarization terms; if those matter, the viscosity balance disappears. The steps that produce D_ef are sketched rather than derived in detail, and Δ, a_T, and f_na enter as external parameters without error bars or sensitivity checks. No comparison to existing simulations or full gyrokinetic runs is given. This is for plasma theorists working on microturbulence and ambipolarity in tokamaks or stellarators. A reader who already thinks about finite-beta effects and chaotic fields will find the idea useful to test. The paper deserves peer review because the mechanism is simple enough to check and the question it raises is real, even if the current version needs tighter derivations and some numerical anchoring.

Referee Report

3 major / 2 minor

Summary. The manuscript claims that finite-beta microturbulence produces unshielded turbulent magnetic perturbations B̃ that ergodize field lines over a radial correlation length Δ. Quasi-neutrality along these chaotic lines enforces the parallel relation en B·∇Φ = B·∇ p_e, which in turn generates an effective electron radial transport coefficient D_ef = (Δ/a_T) T_e/eB. This transport supplies an effective viscosity on the electron flow that can exactly counterbalance a non-ambipolar fraction f_na of the ion gyro-Bohm particle flux, eliminating the need for a radial electric field to confine ions and impurities. The paper calculates the maximum sustainable f_na and the minimum beta required to keep B̃ unshielded.

Significance. If the central relations hold, the result supplies a concrete mechanism by which electromagnetic microturbulence can render ion and electron transport non-ambipolar without invoking an external radial electric field. This would directly affect impurity accumulation predictions in toroidal devices and offers a falsifiable link between measured beta, correlation lengths, and observed non-ambipolar fluxes. The derivation is parameter-light once Δ and a_T are fixed by turbulence scales, and the explicit beta threshold constitutes a testable prediction.

major comments (3)

[Abstract / central derivation] Abstract and central derivation: the explicit algebraic steps that convert the quasi-neutrality condition en B·∇Φ = B·∇ p_e into the diffusion coefficient D_ef = (Δ/a_T) T_e/eB are not displayed. Without these steps, the factor (Δ/a_T) appears post-hoc and the subsequent viscosity balance with f_na × gyro-Bohm diffusion cannot be verified as load-bearing.
[Finite-beta turbulence section] Assumption of unshielded B̃: the claim that finite-beta microturbulence produces an unshielded |B̃| sufficient to ergodize field lines over scale Δ is load-bearing. The manuscript must demonstrate why diamagnetic, kinetic, or finite-Larmor-radius shielding does not reduce |B̃| below the ergodization threshold; otherwise the chaos premise and the derived D_ef both collapse.
[Viscosity balance calculation] Balance with non-ambipolar flux: the effective viscosity arising from the electron transport is asserted to counterbalance exactly f_na times gyro-Bohm diffusion, yet no equation set shows how the viscosity is computed from D_ef or how the numerical limits on f_na and beta are obtained. The absence of error propagation or sensitivity analysis on the free parameters Δ, a_T, f_na, and beta weakens the quantitative claim.

minor comments (2)

[Abstract] Notation: vector symbols B̃ and Φ are introduced without a clear statement of whether they are fluctuating quantities averaged over flux tubes or instantaneous fields; a brief clarification would improve readability.
[Introduction] Missing references: prior literature on chaotic field-line transport in electromagnetic turbulence (e.g., works on stochastic magnetic fields in fusion plasmas) should be cited to place the quasi-neutrality relation in context.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough and constructive report. The comments correctly identify places where the original manuscript lacked explicit derivations and supporting calculations. We have revised the manuscript to supply the missing steps, expand the shielding analysis, and add the viscosity equations together with a sensitivity study. Our point-by-point responses follow.

read point-by-point responses

Referee: [Abstract / central derivation] Abstract and central derivation: the explicit algebraic steps that convert the quasi-neutrality condition en B·∇Φ = B·∇ p_e into the diffusion coefficient D_ef = (Δ/a_T) T_e/eB are not displayed. Without these steps, the factor (Δ/a_T) appears post-hoc and the subsequent viscosity balance with f_na × gyro-Bohm diffusion cannot be verified as load-bearing.

Authors: We agree that the algebraic steps were omitted. In the revised manuscript we have inserted a new paragraph immediately after the quasi-neutrality relation. Starting from en B·∇Φ = B·∇ p_e along chaotic lines, we integrate the potential correlation over the radial distance Δ set by the electron streaming length, average the resulting E×B drift, and obtain D_ef = (Δ/a_T) T_e/eB by direct dimensional reduction. The same paragraph then shows how this diffusivity supplies the viscous torque that balances f_na times the gyro-Bohm ion flux. revision: yes
Referee: [Finite-beta turbulence section] Assumption of unshielded B̃: the claim that finite-beta microturbulence produces an unshielded |B̃| sufficient to ergodize field lines over scale Δ is load-bearing. The manuscript must demonstrate why diamagnetic, kinetic, or finite-Larmor-radius shielding does not reduce |B̃| below the ergodization threshold; otherwise the chaos premise and the derived D_ef both collapse.

Authors: The original text already contains an estimate of the minimum beta needed to keep |B̃| above the ergodization threshold. We acknowledge, however, that explicit comparisons with diamagnetic, kinetic, and FLR shielding were only sketched. We have added a new appendix that evaluates the shielding factors under standard gyrokinetic ordering and demonstrates that, once beta exceeds the reported threshold, the residual |B̃| remains sufficient for ergodization over Δ. The appendix also notes the ordering assumptions under which kinetic shielding remains sub-dominant. revision: partial
Referee: [Viscosity balance calculation] Balance with non-ambipolar flux: the effective viscosity arising from the electron transport is asserted to counterbalance exactly f_na times gyro-Bohm diffusion, yet no equation set shows how the viscosity is computed from D_ef or how the numerical limits on f_na and beta are obtained. The absence of error propagation or sensitivity analysis on the free parameters Δ, a_T, f_na, and beta weakens the quantitative claim.

Authors: We have expanded the viscosity-balance section with the explicit torque-balance equation that equates the electron viscous stress (derived from D_ef via the parallel momentum flux) to f_na times the gyro-Bohm ion flux. The maximum sustainable f_na and the corresponding beta threshold are obtained by solving the resulting algebraic inequality under the unshielded-B̃ constraint. A new sensitivity subsection now varies Δ and a_T over the range 5–10 ρ_i and 0.1–0.3 a, respectively, and reports that the resulting f_na changes by less than 25 %; error bands are added to the relevant figures. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation follows from stated physical assumptions without reduction to inputs

full rationale

The paper starts from the premise that finite-beta microturbulence produces unshielded chaotic B̃, then invokes quasi-neutrality to obtain the parallel relation en B·∇Φ = B·∇ p_e. From this it constructs D_ef = (Δ/a_T) T_e/eB using the externally supplied correlation length Δ and temperature gradient scale a_T. The resulting effective viscosity is then shown to counterbalance a non-ambipolar ion flux f_na times gyro-Bohm diffusion. All quantities are defined independently of the final balance; no parameter is fitted to data and then relabeled a prediction, no self-citation supplies a load-bearing uniqueness theorem, and the central expressions do not collapse to the paper's own inputs by algebraic identity. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

4 free parameters · 3 axioms · 0 invented entities

The central claim rests on standard gyrokinetic assumptions plus the new quasi-neutrality condition on chaotic lines; no new particles or forces are postulated, but several scale lengths are introduced without independent measurement.

free parameters (4)

Δ
Radial correlation distance of the potential set by electron motion along chaotic lines
a_T
Electron temperature gradient scale length defined by |dT_e/dr| = T_e/a_T
f_na
Non-ambipolar fraction of ion radial diffusion expressed in units of gyro-Bohm
beta
Plasma beta threshold to prevent shielding of magnetic perturbations

axioms (3)

domain assumption Gyrokinetic theory restricted to magnetic flux tubes yields identical radial transport for ions and electrons
Opening premise of the abstract
domain assumption Exact magnetic surfaces cease to exist once electrostatic microturbulence is present
Required for the onset of chaos
domain assumption Quasi-neutrality along chaotic lines enforces en B·∇Φ = B·∇ p_e
Central relation used to obtain the potential and transport

pith-pipeline@v0.9.0 · 5548 in / 1645 out tokens · 31583 ms · 2026-05-16T12:34:36.406235+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Quasi-neutrality along the chaotic magnetic field lines requires a potential that obeys en B·∇Φ = B·∇ p_e ... D_ef = (Δ/a_T) T_e/eB ... effective viscosity ... counterbalance a non-ambipolar part ... f_na times gyro-Bohm diffusion
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The maximum f_na that can be counterbalanced and the required plasma beta to avoid shielding the magnetic perturbations B̃ are calculated

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

[1]

Different island chains exert phase dependent forces depending on their relative phase and will lock together at sufficient magnitude

Creates delta-function currents at resonant ra- tional surfaces [8] , which quickly form islands. Different island chains exert phase dependent forces depending on their relative phase and will lock together at sufficient magnitude

work page
[2]

Plasmas 29, 122902 (2022) Boozer & Elder, Phys Plasmas 28, 062303 (2021) Boozer & Elder, Phys

Causes the distortion of flux tubes [7, 9], which allows an arbitrarily smallη/µ 0 to mix lines Þ Huang & Bha*acharjee, Phys. Plasmas 29, 122902 (2022) Boozer & Elder, Phys Plasmas 28, 062303 (2021) Boozer & Elder, Phys. Plasmas 28, 062303 (2021) Huang & Bha>acharjee, Phys. Plasmas 29, 122902 (2022) Top ﬂowing perfect conductor Sides rigid perfect conduct...

work page doi:10.1063/5.0138805 2022
[3]

W. M. Nevins, E. Wang, and J. Candy,Magnetic Stochasticity in Gyrokinetic Simulations of Plasma Microturbulence, Phys. Rev. Lett.,106, 065003 (2011); doi:10.1103/PhysRevLett.106.065003

work page doi:10.1103/physrevlett.106.065003 2011
[4]

J. W. Connor, R. J. Hastie, and A. Zocco,The stochastic field transport associated with the slab ITG modes, Plasma Phys. Control. Fusion55, 12500 (2013); doi: 0.1088/0741-3335/55/12/125003

work page 2013
[5]

P. W. Terry, D. Carmody, H. Doerk, W. Gutten- felder, D. R. Hatch, C. C. Hegna, A. Ishizawa, F. Jenko, W. M. Nevins, I. Predebon, M. J. Pueschel, J.S. Sarff, and G. G. Whelan,Overview of gy- rokinetic studies of finite-βmicroturbulence, Nucl. Fusion,55, 104011 (2015); doi: 10.1088/0029- 5515/55/10/104011

work page doi:10.1088/0029- 2015
[6]

and et al, Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wen- delstein 7-X, Nucl

Alcus´ on, J., Wegner, T., Dinklage, A. and et al, Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wen- delstein 7-X, Nucl. Fusion63, 094002 (2023); doi: 10.1088/1741-4326/aceb76

work page doi:10.1088/1741-4326/aceb76 2023
[7]

Manca, S

P. Helander and A. N. Simakov,Intrinsic am- bipolarity and rotation in stellarators, Phys. Rev. Lett.,101, 14500 (2008); doi: 10.1103/Phys- RevLett.101.145003

work page doi:10.1103/phys- 2008
[8]

A. H. Boozer,Physics of magnetically confined plasmas, Rev. Mod. Phys.76, 1071 (2004); doi: 10.1103/RevModPhys.76.1071

work page doi:10.1103/revmodphys.76.1071 2004
[9]

A. H. Boozer,The rapid destruction of toroidal mag- netic surfaces, Phys. Plasmas29, 022301 (2022); doi: 10.1063/5.0076363

work page doi:10.1063/5.0076363 2022
[10]

Zhou, Y.-M

Y. Zhou, Y.-M. Huang, H. Qin, and A. Bhattachar- jee,Formation of current singularity in a topologi- cally constrained plasma, Phys. Rev. E93, 023205 (2016); doi: 10.1103/PhysRevE.93.023205

work page doi:10.1103/physreve.93.023205 2016
[11]

A. H. Boozer,Magnetic Field Line Chaos, Cantori, and Turnstiles in Toroidal Plasmas, <https://arxiv.org/pdf/2510.25047>

work page arXiv
[12]

A. H. Boozer,Electric field effects during dis- ruptions, Phys. Plasmas31, 102506 (2024); doi: 10.1063/5.0219727

work page doi:10.1063/5.0219727 2024

[1] [1]

Different island chains exert phase dependent forces depending on their relative phase and will lock together at sufficient magnitude

Creates delta-function currents at resonant ra- tional surfaces [8] , which quickly form islands. Different island chains exert phase dependent forces depending on their relative phase and will lock together at sufficient magnitude

work page

[2] [2]

Plasmas 29, 122902 (2022) Boozer & Elder, Phys Plasmas 28, 062303 (2021) Boozer & Elder, Phys

Causes the distortion of flux tubes [7, 9], which allows an arbitrarily smallη/µ 0 to mix lines Þ Huang & Bha*acharjee, Phys. Plasmas 29, 122902 (2022) Boozer & Elder, Phys Plasmas 28, 062303 (2021) Boozer & Elder, Phys. Plasmas 28, 062303 (2021) Huang & Bha>acharjee, Phys. Plasmas 29, 122902 (2022) Top ﬂowing perfect conductor Sides rigid perfect conduct...

work page doi:10.1063/5.0138805 2022

[3] [3]

W. M. Nevins, E. Wang, and J. Candy,Magnetic Stochasticity in Gyrokinetic Simulations of Plasma Microturbulence, Phys. Rev. Lett.,106, 065003 (2011); doi:10.1103/PhysRevLett.106.065003

work page doi:10.1103/physrevlett.106.065003 2011

[4] [4]

J. W. Connor, R. J. Hastie, and A. Zocco,The stochastic field transport associated with the slab ITG modes, Plasma Phys. Control. Fusion55, 12500 (2013); doi: 0.1088/0741-3335/55/12/125003

work page 2013

[5] [5]

P. W. Terry, D. Carmody, H. Doerk, W. Gutten- felder, D. R. Hatch, C. C. Hegna, A. Ishizawa, F. Jenko, W. M. Nevins, I. Predebon, M. J. Pueschel, J.S. Sarff, and G. G. Whelan,Overview of gy- rokinetic studies of finite-βmicroturbulence, Nucl. Fusion,55, 104011 (2015); doi: 10.1088/0029- 5515/55/10/104011

work page doi:10.1088/0029- 2015

[6] [6]

and et al, Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wen- delstein 7-X, Nucl

Alcus´ on, J., Wegner, T., Dinklage, A. and et al, Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wen- delstein 7-X, Nucl. Fusion63, 094002 (2023); doi: 10.1088/1741-4326/aceb76

work page doi:10.1088/1741-4326/aceb76 2023

[7] [7]

Manca, S

P. Helander and A. N. Simakov,Intrinsic am- bipolarity and rotation in stellarators, Phys. Rev. Lett.,101, 14500 (2008); doi: 10.1103/Phys- RevLett.101.145003

work page doi:10.1103/phys- 2008

[8] [8]

A. H. Boozer,Physics of magnetically confined plasmas, Rev. Mod. Phys.76, 1071 (2004); doi: 10.1103/RevModPhys.76.1071

work page doi:10.1103/revmodphys.76.1071 2004

[9] [9]

A. H. Boozer,The rapid destruction of toroidal mag- netic surfaces, Phys. Plasmas29, 022301 (2022); doi: 10.1063/5.0076363

work page doi:10.1063/5.0076363 2022

[10] [10]

Zhou, Y.-M

Y. Zhou, Y.-M. Huang, H. Qin, and A. Bhattachar- jee,Formation of current singularity in a topologi- cally constrained plasma, Phys. Rev. E93, 023205 (2016); doi: 10.1103/PhysRevE.93.023205

work page doi:10.1103/physreve.93.023205 2016

[11] [11]

A. H. Boozer,Magnetic Field Line Chaos, Cantori, and Turnstiles in Toroidal Plasmas, <https://arxiv.org/pdf/2510.25047>

work page arXiv

[12] [12]

A. H. Boozer,Electric field effects during dis- ruptions, Phys. Plasmas31, 102506 (2024); doi: 10.1063/5.0219727

work page doi:10.1063/5.0219727 2024