Geodesic Acoustic Modes (GAM)

• GAMs are oscillatory, axisymmetric zonal flows in toroidal plasmas driven by field-line curvature and plasma compressibility, regulating drift-wave turbulence.
• Their linear behavior is characterized by a frequency scaling as ω_GAM ~ c_s/R, with modifications from collisionless damping, finite radial structure, and plasma shaping.
• Nonlinear interactions, energetic particle effects, and electromagnetic influences lead to complex GAM evolution, impacting turbulence regulation and transport barrier formation.

Geodesic Acoustic Modes (GAMs) are oscillatory, axisymmetric ($m = n = 0$) zonal flows unique to toroidal plasmas, arising from the interplay of field-line curvature, plasma compressibility, and the up–down antisymmetric ($m=1$) pressure (density/temperature) perturbations. They regulate drift-wave turbulence, impact anomalous transport, and are critical to the formation of transport barriers in tokamak and stellarator experiments. GAMs are characterized by a frequency $\omega_{\rm GAM} \sim c_s/R$, where $c_s = \sqrt{T_e/m_i}$ is the sound speed and $R$ is the major radius; they exhibit a broad range of physical phenomena including collisionless (Landau) damping, energetic-particle-driven instabilities (EGAMs), nonlinear interactions with turbulence, intermittency, and a range of kinetic, collisional, and geometric effects.

1. Linear Theory: Dispersion Relations and Key Dependencies

The fundamental frequency and damping of GAMs are captured by both fluid and kinetic models. In a large-aspect-ratio, circular tokamak, the canonical linear GAM frequency for electrostatic, collisionless, adiabatic-electron plasmas is

$$\omega_{\rm GAM}^2 = \frac{T_i + T_e}{m_i R^2}\left( \frac{7}{4} + \frac{T_e}{T_i} \right),$$

with damping arising from ion and electron Landau resonance. Including finite radial structure (finite $k_r$) and shaping (elongation) leads to modified expressions involving transit harmonics and finite orbit width (FOW) (Biancalani et al., 2017, Novikau et al., 2017).

Key Physics Parameters

Parameter	Effect on $\omega_{\rm GAM}$	Effect on damping $\gamma$
Safety factor $q$	Increases $m=1$0 (quadratic scaling); at large $m=1$1, corrections are $m=1$2	Damping increases rapidly with $m=1$3 up to $m=1$4; saturates for larger $m=1$5 (all-transit-harmonic regime)
$m=1$6	Small $m=1$7: negligible, but important at larger $m=1$8; higher harmonics required	Damping $m=1$9 for FOW; at large $\omega_{\rm GAM} \sim c_s/R$0: algebraic decay for damping (Chen et al., 2017)
Elongation $\omega_{\rm GAM} \sim c_s/R$1	Reduces $\omega_{\rm GAM} \sim c_s/R$2 for moderate $\omega_{\rm GAM} \sim c_s/R$3	Damping increases with $\omega_{\rm GAM} \sim c_s/R$4
$\omega_{\rm GAM} \sim c_s/R$5	Increases both frequency and damping
Kinetic electrons	Saturate $\omega_{\rm GAM} \sim c_s/R$6 at high $\omega_{\rm GAM} \sim c_s/R$7; increase damping (factor 2–10)	Adds electron Landau damping

Fully kinetic, drift-kinetic, and gyrokinetic treatments reveal that electrons, FOW, and shaping significantly affect the linear spectrum and damping characteristics (Novikau et al., 2017, Biancalani et al., 2016, Biancalani et al., 2017).

2. GAM Excitation, Nonlinear Evolution, and Intermittency

GAMs are spontaneously excited by nonlinear interactions—primarily Reynolds and Maxwell stress transfer from drift-wave turbulence. Parametric instability and modulational-interaction formalisms lead to spontaneous GAM drive equations of the form

$\omega_{\rm GAM} \sim c_s/R$8

In global and local turbulence models, GAMs typically manifest as limit-cycle oscillations, intermittent bursts, or chirped wave packets (Hager et al., 2011). The nonlinear growth rate due to turbulence can be explicitly calculated, and the dynamical structure can enter regimes of soliton-like propagation, convective turbulence spreading, and spatiotemporal intermittency (Chen et al., 2021, Liu et al., 2020).

A unified two-field theory captures coupled drift-wave and GAM evolution:

• In the linear phase, parametric-modulational equations reproduce standard three-wave instability conditions and enhanced group velocities (joint DW-GAM packets can propagate orders-of-magnitude faster than linear group velocity).
• In the nonlinear regime, soliton-like GAM–DW structures self-organize, producing convective turbulence penetration and core–edge coupling (Chen et al., 2021).

Numerical evidence confirms that “caviton-instanton” transitions and local reversals in drift-wave energy flow trigger intermittent, radially propagating GAM bursts, consistent with experimental observations in ASDEX-U, DIII-D, and T-10 (Liu et al., 2020).

3. Energetic Particle Effects and EGAM Physics

GAMs can be driven unstable by energetic particle (EP) populations (EGAMs), notably via:

• Inverse Landau resonance with positive velocity-space gradient of the EP distribution.
• Absence of pitch-angle instability threshold for non–fully-slowed-down EP beams, introducing strong “simple-pole” drives (in contrast to weaker logarithmic terms for fully-slowed-down distributions) (Cao et al., 2015).
• EGAM frequencies can exceed the local thermal-GAM frequency and exhibit rapid growth (instability timescales far shorter than slowing-down time).

Dispersion relations incorporating EPs (Cao et al., 2015, Novikau et al., 2019, Chen et al., 2017, Qiu et al., 2018):

• Reproduce the observed features of fast EGAM bursts and frequency trends in LHD and ASDEX-Upgrade.
• Show that both the “logarithmic” (classical) and “pole” (threshold-free) terms in the EP response are crucial—pole terms dominating when the EP distribution is incomplete (not fully slowed).
• Demonstrate that at short wavelengths, the EGAM drive and damping both scale algebraically with $\omega_{\rm GAM} \sim c_s/R$9 and EP drift frequency, rather than exponentially.

MPR (Mode-Particle-Resonance) diagnostics in global gyrokinetic simulations quantitatively resolve energy transfer between EGAMs, bulk ions, and EP/thermal electrons. Typical heating rates indicate dominant ion heating, with kinetic electron effects reducing EGAM amplitude and associated heating (Novikau et al., 2019).

4. Electromagnetic, Collisional, and Geometric Effects

Electromagnetic (EM) GAMs

• Including EM fluctuations modifies both the restoring force and the electron pressure/parallel current response.
• The electromagnetic GAM frequency becomes sensitive to electron temperature gradients (not present in the electrostatic limit), with the EM branch corresponding to larger radial wavelengths (up to $c_s = \sqrt{T_e/m_i}$0), consistent with high-frequency oscillations in devices like TCABR (Sgalla, 2014).
• m=1 and m=2 poloidal magnetic sidebands coexist in EM-GAMs, with the m=1 amplitude strongly enhanced by drift-coupled harmonics and potentially comparable to or dominating m=2 (Xie et al., 2022). Both scale with $c_s = \sqrt{T_e/m_i}$1, and must be incorporated when $c_s = \sqrt{T_e/m_i}$2.

Collisionality and Plasma Beta

• Collisions up-shift the GAM frequency at low $c_s = \sqrt{T_e/m_i}$3 and down-shift at high $c_s = \sqrt{T_e/m_i}$4.
• Finite $c_s = \sqrt{T_e/m_i}$5 reduces $c_s = \sqrt{T_e/m_i}$6 in all regimes due to magnetic compressibility.
• Coupled system including surface-averaged density and temperature evolution leads to emergent low-frequency branches and can exacerbate experimental–theoretical discrepancies in observed GAM frequency, indicating that higher $c_s = \sqrt{T_e/m_i}$7 coupling and nonlinear effects play non-negligible roles (Singh et al., 2015).

Plasma Geometry

• Shaping (elongation, triangularity, X-point symmetry) influences fundamental frequency, growth/damping rates, and even the directionality and pulse-structure of GAM activity, including selection by up–down asymmetry and preference for specific $c_s = \sqrt{T_e/m_i}$8 signs in single-null configurations, leading to pulsed GAM activity (Hager et al., 2011).

Second Harmonic and Infinite-m Coupling

• Retaining $c_s = \sqrt{T_e/m_i}$9 sidebands in fluid/gyrofluid models systematically increases the GAM frequency and modifies the nonlinear excitation, particularly in shaped plasmas or in the presence of strong gradients (Anderson et al., 2014).
• Infinite $R$0-chains (matrix continued-fraction solutions) demonstrate that neglect of high-$R$1 couplings overestimates experimental GAM frequencies by up to a factor of two in low-collisionality discharges; convergence with experimental measurements is only obtained when these couplings are included (Singh et al., 2015).

5. Decay, Damping, and Phase Mixing

The lifetime and amplitude decay of GAMs in realistic plasmas are dominated by the interplay between:

• Landau (collisionless) damping, which exponentially suppresses GAMs via wave–particle resonance (explicit formulas given in, e.g., (Biancalani et al., 2016)).
• Phase mixing (continuum damping) arising from equilibrium profile gradients (notably the ion temperature profile $R$2), which increases $R$3 in time and thus strongly accelerates the decay rate by pumping the mode into regimes of higher Landau damping (Biancalani et al., 2016, Novikau et al., 2017).
• Combined, these effects can reduce the GAM half-decay time by over an order of magnitude compared to collisionless damping alone, imposing stringent limits on the persistence and transport-regulatory role of GAMs in steep-gradient or strongly inhomogeneous pedestals.

Nonlinear predator–prey dynamics between GAMs and background microturbulence (e.g., ITG, ETG modes) yield elevated turbulence saturation and feedback, with high-frequency GAM branches possible through electron-scale ETG turbulence coupling (Anderson et al., 2012).

6. Anisotropy, Rotation, and Kinetic/Geometric Extensions

• Analytical and gyrokinetic treatments with bi-Maxwellian ions ($R$4) reveal that increasing ion temperature anisotropy ($R$5) raises both $R$6 and $R$7, especially in high $R$8 regimes (Sama et al., 2022, Ren, 2015).
• Toroidal rotation adds further up-shift of frequency and reduction of damping, with explicit dependence on the Mach number and anisotropy threshold (Ren, 2015).
• Kinetic extensions allow for realistic experimental geometry, mass ratio, elongation, and nonuniformity, enabling robust comparison and cross-validation against multi-code gyrokinetic benchmarks (Biancalani et al., 2017, Novikau et al., 2017).

Key foundational and review works:

• "Cross-code gyrokinetic verification and benchmark on the linear collisionless dynamics of the geodesic acoustic mode" (Biancalani et al., 2017)
• "Linear gyrokinetic investigation of the geodesic acoustic modes in realistic tokamak configurations" (Novikau et al., 2017)
• "Kinetic theory of geodesic acoustic modes in toroidal plasmas: a brief review" (Qiu et al., 2018)
• "Geodesic Acoustic Modes with poloidal mode couplings ad infinitum" (Singh et al., 2015)
• "Nonlinear dynamics of energetic-particle driven geodesic acoustic modes in ASDEX Upgrade" (Novikau et al., 2019)

Conclusion

GAMs embody a broad suite of physics effects central to turbulence regulation and transport in toroidal plasmas. Their frequency, damping, and stability properties are conditioned by kinetic effects, plasma geometry, energetic particle content, electromagnetic dynamics, collisionality, and strong nonlinear coupling to turbulence. The theoretical apparatus for GAMs is now sufficiently mature to interpret and predict most experimental features within the constraints of kinetic and global toroidal geometry, but continuing discrepancies (notably frequency down-shifts and robustness in low-collisionality regimes) motivate further work in nonlocal/nongyrotropic closure, high-$R$9 spectral resolution, and inclusion of full geometric and collisional complexity.