---
title: Geodesic Acoustic Modes (GAM)
url: https://www.emergentmind.com/topics/geodesic-acoustic-modes-gam
type: topic
---

# Geodesic Acoustic Modes (GAM)

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) [1705.06554, 1709.01818].

### Key Physics Parameters

| Parameter      | Effect on $\omega_{\rm GAM}$                       | Effect on damping $\gamma$           |
|----------------|----------------------------------------------------|--------------------------------------|
| Safety factor $q$  | Increases $\omega_{\rm GAM}$ (quadratic scaling); at large $q$, corrections are $O(q^{-2})$ | Damping increases rapidly with $q$ up to $q \sim 3.5$; saturates for larger $q$ (all-transit-harmonic regime) |
| $k_r \rho_i$   | Small $k_r$: negligible, but important at larger $k_r$; higher harmonics required | Damping $\propto (k_r \rho_i)^2$ for FOW; at large $k_r$: algebraic decay for damping ([1709.02085]) |
| Elongation $e$ | Reduces $\omega_{\rm GAM}$ for moderate $e$        | Damping increases with $e$           |
| $T_e/T_i$      | Increases both frequency and damping               |                                     |
| Kinetic electrons | Saturate $\omega_{\rm GAM}$ at high $k_r$; 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 [1709.01818, 1608.03447, 1705.06554].

## 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
\[
\frac{\partial^2 \delta \phi}{\partial t^2} + \omega_{\rm GAM}^2\, \delta\phi = \text{nonlinear forcing} .
\]
In global and local turbulence models, GAMs typically manifest as limit-cycle oscillations, intermittent bursts, or chirped wave packets [1109.3324]. 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 [2106.04131, 2007.00250].

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 [2106.04131].

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 [2007.00250].

## 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) [1510.06105].
- 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 [1510.06105, 1912.07950, 1709.02085, 1801.01622]:
- 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 $k_r$ 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 [1912.07950].

## 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 $\lambda_r \sim 25\,\mathrm{cm}$), consistent with high-frequency oscillations in devices like TCABR [1410.6827].
- 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 [2203.09024]. Both scale with $\beta q^2$, and must be incorporated when $\beta q^2 \gtrsim 1\%$.

### Collisionality and Plasma Beta
- Collisions up-shift the GAM frequency at low $k_r$ and down-shift at high $k_r$.
- Finite $\beta$ reduces $\omega_{\rm GAM}$ 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 $m$ coupling and nonlinear effects play non-negligible roles [1507.03232].

### 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 $k_r$ signs in single-null configurations, leading to pulsed GAM activity [1109.3324].

### Second Harmonic and Infinite-m Coupling
- Retaining $m=2$ 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 [1407.8037].
- Infinite $m$-chains (matrix continued-fraction solutions) demonstrate that neglect of high-$m$ 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 [1512.06066].

## 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., [1608.03447]).
- Phase mixing (continuum damping) arising from equilibrium profile gradients (notably the ion temperature profile $T_i(r)$), which increases $k_r$ in time and thus strongly accelerates the decay rate by pumping the mode into regimes of higher Landau damping [1608.03447, 1709.01818].
- 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 [1203.1174].

## 6. Anisotropy, Rotation, and Kinetic/Geometric Extensions

- Analytical and gyrokinetic treatments with bi-Maxwellian ions ($T_\perp \neq T_\parallel$) reveal that increasing ion temperature anisotropy ($\eta = T_\perp/T_\parallel$) raises both $\omega_{\rm GAM}$ and $|\gamma|$, especially in high $T_e/T_i$ regimes [2209.14874, 1503.01315].
- Toroidal rotation adds further up-shift of frequency and reduction of damping, with explicit dependence on the Mach number and anisotropy threshold [1503.01315].
- Kinetic extensions allow for realistic experimental geometry, mass ratio, elongation, and nonuniformity, enabling robust comparison and cross-validation against multi-code gyrokinetic benchmarks [1705.06554, 1709.01818].

## References

Key foundational and review works:
- "Cross-code gyrokinetic verification and benchmark on the linear collisionless dynamics of the geodesic acoustic mode" [1705.06554]
- "Linear gyrokinetic investigation of the geodesic acoustic modes in realistic tokamak configurations" [1709.01818]
- "Kinetic theory of geodesic acoustic modes in toroidal plasmas: a brief review" [1801.01622]
- "Geodesic Acoustic Modes with poloidal mode couplings ad infinitum" [1512.06066]
- "Nonlinear dynamics of energetic-particle driven geodesic acoustic modes in ASDEX Upgrade" [1912.07950]

## 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-$m$ spectral resolution, and inclusion of full geometric and collisional complexity.

Source: https://www.emergentmind.com/topics/geodesic-acoustic-modes-gam