--- title: Bootstrap Current Fractions in Fusion Plasmas url: https://www.emergentmind.com/topics/bootstrap-current-fractions type: topic --- # Bootstrap Current Fractions in Fusion Plasmas Bootstrap Current Fractions The bootstrap current fraction quantifies the proportion of plasma current in a magnetically confined system (tokamak, stellarator, or helical device) generated by neoclassical transport processes, primarily through particle drifts and collisional dynamics, as opposed to externally driven currents. This fraction, typically denoted \( f_{\rm bs} \), is critical for both device optimization and predictive modeling, affecting the equilibrium, stability, and confinement properties of fusion-grade plasmas. The accurate computation of bootstrap current fractions requires sophisticated kinetic models, careful treatment of drift- and collisional operators, and, in non-axisymmetric systems, enforcement of exact parallel momentum conservation. Bootstrap current fractions are also of interest in the conformal bootstrap of quantum field theory, where they measure the relative OPE weight carried by conserved current sectors. ## 1. Definitions and Formulation Bootstrap current fraction is defined as the ratio of the bootstrap current to the total current on a flux surface: \[ f_{\rm bs} = \frac{j_{\rm bs}}{j_{\rm tot}} \] or, in an integrated form, \[ f_{\rm bs} = \frac{I_{\rm bs}}{I_{\rm tot}} \] where \( j_{\rm bs} \) is the local bootstrap current density and \( j_{\rm tot} \) is the total (typically toroidal) current density. The precise kinetic definition, as encoded in drift-kinetic models for non-axisymmetric systems, is \[ j_{\rm bs} = e\,n\,\bigl(U_{\parallel,i}-U_{\parallel,e}\bigr)\frac{\langle B\rangle}{\langle B^2\rangle} \] with flux-surface-averaged parallel flows \( U_{\parallel,a} \) for each species \( a \) [1611.03210]. In conformal field theoretic bootstrap, the analogue, for mixed correlators involving a scalar and a conserved current \( J_\mu \), is the OPE “current fraction”—i.e., the ratio of contributions from conserved spin-1 currents to the total conformal block sum [1911.05747, 2512.20803]: \[ \mathrm{fraction}_J(u,v) = \frac{ \sum_{\Delta,J=1}\lambda_{\phi\bar\phi J}^2\,G_{\Delta,1}(u,v) }{ \sum_{\rm all}\lambda^2_{\mathit{Q},\ell}\,G_{\Delta,\ell}(u,v) } \] ## 2. Theoretical Models and Calculation Methodologies Reliable calculation of bootstrap fractions requires both kinetic theory and global plasma modeling. - **Local neoclassical theory:** Formulas such as those of Sauter et al. (1999) express \( j_{\rm bs} \) as a linear combination of normalized density and temperature gradients, weighted by dimensionless coefficients \( L_{31}, L_{32}, L_{34} \) that interpolate across the banana, plateau, and Pfirsch-Schlüter regimes: \[ j_{\rm bs} = -\,\frac{c\,I}{B\langle B^2\rangle} [ L_{32}\frac{dT_e}{d\psi} + L_{34}\frac{dT_i}{d\psi} ] \] [1207.1795]. - **Full drift-kinetic solvers:** For non-axisymmetric (stellarators/helical) systems, it is essential to solve the drift-kinetic equation with exact parallel momentum conservation: \[ v_{\parallel}\,\mathbf{b}\cdot\nabla f_{1a} + \mathbf{v}_d\cdot\nabla_\perp f_{0a} = \mathcal{C}_a[f_{1a}] \] and ensure the surface-averaged momentum balance is respected [1611.03210, 2510.27513]. - **Code benchmarks:** Three principal computational approaches are compared: - DKES: Monoenergetic local drift-kinetic solver employing only pitch-angle scattering. - ZOW: δf Monte Carlo code, energy-dependent, includes tangential drift and momentum-conserving like-species operators. - PENTA: DKES transport coefficients with Sugama-Nishimura momentum correction. ZOW and PENTA, which enforce exact momentum balance, return consistent results, while DKES (without correction) can overestimate the bootstrap-fraction by 600–800% [1611.03210]. - **Strong-gradient/pedestal theory:** In regions where the density/temperature scale length is comparable to the ion orbit width (e.g., tokamak pedestal), standard local neoclassical predictions must be amended by global Fokker-Planck or δf codes which include finite-orbit-width (FOW) and strong gradient modifications [1207.1795, 1206.5912, 2504.03016]. ## 3. Parameter Dependence and Scaling Laws The bootstrap fraction depends sensitively on collisionality, magnetic geometry, and plasma profiles: - **Collisionality (\( \nu_* \))**: At low \( \nu_* \) (banana regime), \( f_{\rm bs} \) rises as trapped particle effects dominate, and falls off as \( 1/\nu_* \) in the Pfirsch-Schlüter regime. In the \( 1/\nu \) regime for standard stellarators/tokamaks, the bootstrap coefficient \( D_{31} \sim 1/\nu_* \), but in omnigenous or optimized quasi-symmetric fields, \( D_{31} \) saturates or even decreases as \( \sqrt{\nu_*} \) [2510.27513, 1701.02501, 2407.21599]. - **Ripple amplitude and field symmetry:** The presence of magnetic ripple increases the trapped particle fraction and alters the transport scaling. In omnigenous fields (including quasi-isodynamic or piecewise-omnigenous configurations), the average radial drift of all orbits vanishes, suppressing the bootstrap current entirely (\( f_{\rm bs} = 0 \)) [2505.02546]. For generic non-axisymmetric fields, alignment of field maxima and equivalent ripple conditions are necessary to minimize bootstrap offsets and attain the Shaing-Callen limit [2407.21599]. - **Gradient strength:** The magnitude and alignment of density and temperature gradients, as well as the flow profiles (e.g., strong radial electric fields \( E_r \)), directly drive the bootstrap fraction. In strong-gradient regimes, FOW corrections can reduce or enhance \( f_{\rm bs} \) by up to 30%, particularly in the edge/pedestal region [2504.03016]. - **Aspect ratio and global geometry:** Lower aspect ratios and strong shaping (elongation, triangularity) reduce bootstrap drive compared to the large-aspect-ratio, circular case [2303.00415]. ## 4. Zero Bootstrap Current and Control Strategies Recent advances have enabled the design of stellarator configurations with identically zero bootstrap current for all profiles—a property not present in generic designs. - **Piecewise-omnigenous fields:** In such configurations, the geometric parameter \( \Delta \) defined from the surface partition vanishes, yielding \( j_{\rm bs} \equiv 0 \) independently of collisionality or plasma profiles. This condition is accessible to analytic design and has been numerically validated with high precision. These fields provide reference cases for code benchmarking and serve as a baseline for quantifying symmetry-breaking effects [2505.02546]. - **QA + piecewise-omnigenous ensembles (QA-pwO):** By blending quasi-axisymmetry with piecewise-omnigenous perturbations, it is possible to tune the bootstrap current continuously from well above the axisymmetric value to zero or negative, which is critical for matching divertor operation constraints (such as in island-diverted reactors) or for enhancing bootstrap drive in steady-state tokamaks [2603.20125]. ## 5. Benchmarks, Code Validation, and Physical Limits Robust predictive capability requires cross-verification between analytic models and numerical solvers: | Configuration | Model/Code | f_bs Values (sample) | Agreement/Deviation | |--------------------------|------------------|-------------------------------|-----------------------------| | FFHR-d1 helical reactor | DKES | 0.35 (core) | Overestimates by ~10x | | | ZOW | 0.04 (core) | Consistent with PENTA | | | PENTA | 0.05 (core) | Consistent with ZOW | | Tokamak (M3D-C1) | Sauter/Redl | — | Agreement with NEO, XGCa, SFINCS <3% | | TJ-II stellarator | DKES/NEO-MC | F_bs ≈ 0.9–1.1 | Consistent with experiment | | Stellarator (piecewise-omnigenous) | MONKES | D_{31} crosses zero at Δ=0 | Machine-precision zero | Numerical codes such as MONKES (Legendre-based solver), SFINCS (full 4D global kinetic), and M3D-C1 (extended-MHD with neoclassical closure) enable rapid and accurate estimation of \( D_{31} \), bootstrap fractions, and their optimization in design workflows [2510.27513, 2507.05166, 2407.21599]. ## 6. Experimental Evidence and Practical Implications Experimental validation on devices such as TJ-II demonstrates the physical reality and diagnostic accessibility of bootstrap current fractions: - **TJ-II:** For various density and temperature regimes, calculated and measured total bootstrap currents agree within 10–20%, with the bootstrap contribution often dominating the total plasma current (\( F_{bs} \approx 0.9-1.1 \)), including sign reversals at high density [1108.3721]. - **Tokamak pedestals:** In DIII-D and ASDEX Upgrade-like pedestals, local neoclassical formulae systematically overpredict \( f_{\rm bs} \) by 20–100% in the plateau or strong-gradient regime; global codes including FOW and E_r reduce this discrepancy [1207.1795, 1206.5912, 2504.03016]. - **Optimization and self-consistency:** In device optimization, particularly for quasi-symmetric stellarators, imposing a penalty function for mismatch between equilibrium and analytically predicted bootstrap current profiles yields self-consistent configurations with optimized confinement and controlled \( f_{\rm bs} \) [2205.02914]. ## 7. Generalization to Field Theory: Bootstrap "Current Fractions" in the CFT Context In conformal field theory, current fractions formalize the share of OPE weight in mixed scalar-current correlation functions that is due to conserved current exchange, with direct connections to finite-temperature transport (e.g., conductivity in the 3d O(2) model). - **Numerical bounds (O(2) model):** At the crossing-symmetric point \( u=v=1/4 \), the current fraction is \( \sim 0.16(2) \), while the stress tensor fraction is \( \sim 0.06(1) \) [1911.05747]. - **Large charge bootstrap:** At large global charge, the current fraction in the OPE is determined by the number of Regge trajectories and related to EFT parameters, with closed analytic formulas for one- and two-trajectory solutions [2512.20803]. - **Physics connection:** These fractions directly relate to high-frequency conductivity in 2+1d CFTs via a Kubo formula, with bootstrap-predicted parameters yielding precise agreement with QMC data. --- ## References - [1611.03210] Benchmark of the bootstrap current simulation in helical plasmas - [2505.02546] A new class of optimized stellarators with zero bootstrap current - [2510.27513] Fast and accurate calculation of the bootstrap current and radial neoclassical transport in low collisionality stellarator plasmas - [1911.05747] Mixed Scalar-Current bootstrap in three dimensions - [2507.05166] Bootstrap Current Modeling in M3D-C1 - [2504.03016] Strong gradient effects on neoclassical electron transport and the bootstrap current - [1701.02501] Stellarator bootstrap current and plasma flow velocity at low collisionality - [2303.00415] Gyrokinetic simulations of neoclassical electron transport and bootstrap current generation in tokamak plasmas in the TRIMEG code - [2603.20125] Control of the bootstrap current in approximately quasi-axisymmetric magnetic fields - [1108.3721] Calculation of the bootstrap current profile for the TJ-II stellarator - [2407.21599] On the convergence of bootstrap current to the Shaing-Callen limit in stellarators - [2205.02914] Optimization of quasisymmetric stellarators with self-consistent bootstrap current and energetic particle confinement - [1207.1795] Local and global Fokker-Planck neoclassical calculations showing flow and bootstrap current modification in a pedestal - [1206.5912] Changes to neoclassical flow and bootstrap current in a tokamak pedestal - [2512.20803] Large charge bootstrap with U(1) current probes Source: https://www.emergentmind.com/topics/bootstrap-current-fractions