Plasma Beta

Since beta (β) = p / (B²/2μ₀) and pressure p = n k T, we can express density as n = (β B²)/(2 μ₀ k T). Substituting into the triple product n·T·τ_E gives n·T·τ_E = (β B² τ_E)/(2 μ₀ k). Thus, for a fixed magnetic field B and confinement time τ_E, higher beta directly increases the triple product, bringing the plasma closer to ignition. A factor of 2 in beta doubles the triple product, which is why achieving high beta is not just an efficiency metric but a direct path to reducing the required size or field for ignition.

Yes, there is an absolute upper bound known as the 'Bishop's limit' or the 'ideal MHD beta limit' for a given safety factor profile. For a toroidal plasma, ideal MHD stability theory predicts that the maximum beta is limited by the ballooning instability. The theoretical limit scales as β_max ∝ (a/R) · (1/q²) · (1 + κ²) for a given pressure profile, but this is typically higher than the empirical Troyon limit. The absolute maximum for any stable configuration is often cited as β ~ 40-50% for a spherical tokamak, but in practice, the limit is much lower due to neoclassical tearing modes and resistive wall modes, which arise from finite resistivity and non-ideal effects.

The Troyon beta limit is often expressed as β_N = β / (I_p/(a B)) ≈ 2.8 to 3.5 for conventional aspect ratios (A > 3). However, this limit scales with elongation and inversely with aspect ratio. For spherical tokamaks (A < 1.5), the stability improves because the magnetic field line curvature is more favorable (the magnetic well deepens) and the plasma can achieve higher elongation. The empirical scaling shows that the normalized beta limit increases roughly as β_N ∝ 1/A, so a spherical tokamak with A=1.5 can achieve β_N up to 6-7, corresponding to volume-averaged beta values of 20-40% because the denominator (I_p/(aB)) is also smaller. This makes spherical tokamaks attractive for compact, high-beta designs.

Volume-averaged beta is calculated from the global plasma pressure (obtained from diamagnetic loops) and the known magnetic field. To measure the local beta profile, one needs spatially resolved density and temperature profiles: n(r) from interferometry or Thomson scattering, T_e(r) from Thomson scattering or ECE, and T_i(r) from charge-exchange spectroscopy or neutron spectrometry. Then p(r) = n_e(r) k T_e(r) + n_i(r) k T_i(r). The local beta is β(r) = p(r) / (B(r)²/(2μ₀)). This profile is crucial because instabilities like NTMs depend on the local pressure gradient (dp/dr), not just the global β. The ratio of the peak beta to the average beta (the 'peaking factor') is a key parameter in stability analysis.

The bootstrap current is proportional to the pressure gradient (∇p) and inversely proportional to the collisionality. In a high-beta plasma, the pressure gradient is steep, and the bootstrap current density can reach a significant fraction of the total current (up to 80% in advanced scenarios). The fraction f_bs = I_bs / I_p scales roughly with β_p (poloidal beta) and the aspect ratio. For steady-state operation, the external current drive can be minimized, reducing the recirculating power. However, if the bootstrap current profile becomes too peaked, it can drive NTMs. Therefore, optimizing the pressure profile to maximize bootstrap current while remaining MHD-stable is a central design challenge for reactors like ITER's steady-state mission.

Global beta (β) is the volume-averaged ratio, while local beta (β(r)) is the radial profile. The fusion power density S(r) ∝ n(r)² ⟨σv⟩(T(r)) E, and since n(r) ∝ β(r) B² / T(r), the local fusion power is determined by the local beta and temperature. A peaked beta profile (high core beta) yields a high central power density, but may be unstable. The global beta gives a macroscopic figure of merit for the whole device, but the local beta at the pedestal and core determines both the stored energy and the stability limits. For example, the pedestal beta determines the H-mode edge barrier, which in turn sets the global confinement.

Plasma rotation, particularly toroidal rotation, can stabilize resistive wall modes (RWMs) by Doppler-shifting the mode frequency and creating a 'rotational stabilization' effect. The rotation shear (change in rotation speed with radius) also suppresses turbulence, which indirectly allows a higher pressure gradient. However, the rotation itself is driven by neutral beam torque or intrinsic mechanisms, and it dissipates via viscosity. In ITER, the rotation is expected to be lower, making active feedback control of RWMs necessary. The rotation also affects the beta limit through the Mach number, which modifies the pressure balance. Higher rotation generally allows higher beta before instability, but it introduces additional engineering complexity (e.g., momentum confinement).

ITER's baseline scenario aims for a volume-averaged beta of about 2.3% (β_N ≈ 1.8) to achieve Q=10, operating safely below the Troyon limit of β_N ≈ 2.8. This conservative beta ensures robust stability against NTMs and RWMs. The advanced scenario, designed for longer pulses and higher Q, targets a higher β_N ≈ 2.5–3.0, corresponding to β ≈ 2.8%, achieved by optimizing the current profile and using active NTM stabilization. This requires operating closer to the stability boundary and relies on EC current drive for NTM suppression. The choice between scenarios reflects a trade-off between performance and risk.

A perfectly conducting wall ideally stabilizes the kink mode, allowing beta to exceed the free-boundary Troyon limit indefinitely. In reality, walls have finite resistivity (time constant τ_w = μ₀ σ d², where d is the wall thickness and σ is conductivity). For typical fusion devices, τ_w is on the order of 10–100 ms. The RWM grows on the timescale of τ_w; if the plasma rotation is faster than the growth rate, the mode is stabilized. If rotation slows down (e.g., during a disruption or NBI off), the RWM can grow, leading to disruption. Active stabilization with feedback coils is used to artificially maintain stability, allowing operation at beta values up to the ideal wall limit (β_N up to ~4) even when rotation is low.

In stellarators, beta is limited not by current-driven modes (since there is no net current) but by the equilibrium response to pressure. The primary limit is the 'finite-beta equilibrium shift' or the 'Shafranov shift' which distorts the vacuum magnetic surfaces. At high beta, the plasma's pressure gradient causes a large outward shift, which can degrade the magnetic well and lead to ballooning instability. Additionally, the confinement of fast particles can be degraded. For optimized stellarators like W7-X, which are designed to minimize the Shafranov shift, the beta limit is predicted to be about 3-5% (volume-averaged) for the standard configuration. However, advanced optimization may push this to ~6-7%, but still lower than spherical tokamaks due to the lack of a stabilizing toroidal current and the more complex 3D shaping.

The maximum beta is limited by stability, but the resulting fusion power (which scales as β² B⁴) determines the neutron and heat flux to the first wall. For a given fusion power, higher beta allows a smaller device, but the wall loading (MW/m²) increases. The divertor must handle the exhaust heat, which scales with the total fusion power. If beta is too high for a given B and size, the neutron flux may exceed material limits (~5 MW/m² for steady-state). Therefore, the beta limit is not just a physics limit but a coupled engineering constraint: the reactor must be designed with sufficient wall area and cooling to handle the high heat and neutron fluxes that come with high-beta operation. This often forces a compromise—operating at a lower beta to keep wall loading manageable, or increasing the size to reduce wall loading, which increases cost.

The total plasma pressure includes both the thermal component (from bulk ions and electrons) and the fast-ion component (from energetic particles). Fast ions have a distinct pressure profile and can be highly anisotropic. While the total beta is the sum (β_total = β_thermal + β_fast), the instability driving mechanisms differ: thermal pressure drives ballooning modes (which depend on the pressure gradient), while fast ions can drive Alfvén eigenmodes (AE) and fishbones, which are not captured by ideal MHD beta limits. The fast-ion beta (β_f = p_f/(B²/2μ₀)) is typically a few percent and can stabilize or destabilize certain modes. For example, a high β_f from alpha particles can drive toroidal Alfvén eigenmodes (TAEs) that expel fast ions, reducing heating efficiency. Thus, reactor designs must account for the fast-ion beta profile separately, which is why integrated modeling includes both thermal and fast particle kinetics.