Gyrokinetic Vlasov–Maxwell System Overview

• Gyrokinetic Vlasov–Maxwell system is a set of equations derived by averaging fast cyclotron motions to describe plasma dynamics in reduced phase space.
• Its Hamiltonian formulation using a noncanonical Poisson bracket guarantees conservation of energy, momentum, and a hierarchy of Casimir invariants.
• The framework supports structure-preserving discretization and incorporates extensions such as collisional effects, variational principles, and gauge-free formulations.

The gyrokinetic Vlasov–Maxwell system is the self-consistent set of equations governing the evolution of weakly collisional, strongly magnetized plasmas with turbulence and electromagnetic fluctuations on scales comparable to the Larmor radius. It arises from a systematic dynamical reduction, in which the fast cyclotron (gyro) motion is averaged out to yield a description of plasma dynamics in five-dimensional (guiding center) or six-dimensional (gyrocenter) phase space. The system is formulated as a noncanonical Hamiltonian field theory featuring a nontrivial Poisson bracket, a rich hierarchy of Casimir invariants, and exact energy and momentum conservation properties. Extensions include collisional effects, metric dissipation, variational formulations (Euler–Poincaré and metriplectic), and amenability to structure-preserving discretization.

1. Hamiltonian Structure and Variational Principles

The Hamiltonian formulation of the gyrokinetic Vlasov–Maxwell (GVM) system is constructed on the phase-space manifold

$$Q \times \text{(fields)} = \{ (f_s, D, B) \mid s = 1, \ldots, N_s \},$$

where $f_s$ is the species-s gyrocenter distribution on $TQ$, $D$ the electric displacement 2-form, and $B$ the magnetic 2-form (with $dB = 0$) (Burby et al., 2014). The total Hamiltonian functional is

$$H_{\rm GVM}[f, D, B] = {\mathcal K}[f, \hat E, B] + \frac{1}{8\pi} \int_Q (\hat E \wedge *\hat E + B \wedge *B),$$

where $\hat E = \hat E(f, D, B)$ is determined by the constitutive relation

$$D = \hat E - 4\pi \frac{\delta {\mathcal K}}{\delta E}.$$

Here, $${\mathcal K}[f, E, B] = \sum_s \int_{TQ} K_s(E, B) f_s$$ is the total gyrocenter kinetic energy, and $f_s$0 is the gyrokinetic polarization.

This Hamiltonian structure is rigorously derived either through direct variational principles (Eulerian or constrained Eulerian), noncanonical brackets, or via lifted Poisson brackets through Lie-transform perturbation theory (Burby et al., 2014, Brizard et al., 2016, Squire et al., 2013, Brizard, 2017, Brizard, 2021, Brizard, 2019).

2. Noncanonical Poisson Bracket and Casimir Invariants

The dynamic evolution is governed by a noncanonical Poisson bracket. For functionals $f_s$1, $f_s$2 of $f_s$3,

$f_s$4

where $f_s$5 is defined by the gyrocenter symplectic form, and $f_s$6 is the pullback from Q to $f_s$7 (Burby et al., 2014, Brizard et al., 2016).

Casimir invariants take the form:

• Particle-label Casimirs: $f_s$8, with $f_s$9 arbitrary.
• Field Casimirs: Any functional of $TQ$0 is a Casimir, corresponding to the conservation of Gauss’s law $TQ$1 (Burby et al., 2014, Squire et al., 2013).

The bracket is verified to be antisymmetric and to satisfy the Jacobi identity, ensuring a bona fide Hamiltonian system (Brizard, 2021, Brizard et al., 2016).

3. Equations of Motion and Electromagnetic Closure

The evolution equations consist of:

• Gyrokinetic Vlasov equation: For each $TQ$2,

$TQ$3

or, in the characteristic form in phase space (gyrocenter or guiding-center coordinates),

$TQ$4

where trajectories are generated by the Hamiltonian via the noncanonical Poisson bracket (Burby et al., 2014, Tronko et al., 2016, Brizard, 2019).

• Maxwell's equations (with gyrokinetic constitutive relations):

$TQ$5

$TQ$6

where the polarization and magnetization terms are systematically derived from the gyrocenter Hamiltonian, and $TQ$7 arises from collisional metriplectic extensions (Hirvijoki et al., 2022).

The polarization and magnetization enter Maxwell’s equations via

$TQ$8

where explicit expressions are given, typically as moments over the gyrocenter distribution (Tronko et al., 2017, Brizard, 2018).

4. Energy, Momentum, and Entropy Conservation

The system possesses an exact energy conservation law, derivable by Noether’s theorem from the gyrokinetic action: $TQ$9 valid in both the collisionless (Burby et al., 2014, Brizard, 2017, Tronko et al., 2016) and metriplectic collisional (Hirvijoki et al., 2022) frameworks.

Canonical and toroidal angular momentum are also exactly conserved in an axisymmetric background: $D$0 where $D$1 includes both gyrocenter and field angular momentum (Brizard, 2019, Brizard, 2021, Tronko et al., 2016, Brizard, 2017).

In metriplectic (collisional) frameworks, the entropy functional,

$D$2

is a Casimir of the Poisson bracket and is strictly non-decreasing under the metric (dissipative) bracket, ensuring monotonic entropy production (Hirvijoki et al., 2022).

5. Gyrocenter Coordinate Reductions and Model Hierarchy

The systematic two-step coordinate transformation (guiding-center, then gyrocenter) employs near-identity (Lie or polynomial) transforms to remove gyrophase dependence order-by-order in the small parameter $D$3 (Tronko et al., 2016, Zoni et al., 2019). This reduction yields a gyrokinetic Hamiltonian

$D$4

with $D$5, $D$6 encoding field perturbations (e.g., $D$7), and $D$8 containing polarization and ponderomotive terms (Brizard, 2017, Brizard, 2018).

Different ordering limits—electrostatic, electromagnetic/low-β, full-β, long-wavelength (second-order FLR), and full-FLR—produce a hierarchy of models implemented in major codes (ORB5, GENE) (Tronko et al., 2017). Ion and electron dynamics are treated consistently, allowing for variations such as maximal ordering for fusion devices (Zoni et al., 2019).

6. Structure-Preserving Discretization and Numerical Applications

Structure-preserving numerical schemes are developed by discretizing fields using compatible finite-element exterior calculus (FEEC) bases and the distribution using macroparticles. The resulting finite-dimensional Hamiltonian system inherits the noncanonical bracket and all continuous invariants (energy, Casimirs, Gauss constraints), enabling symplectic time integration and long-term invariance preservation (Burby, 2016, Squire et al., 2013). Discrete Poisson brackets, Dirac reductions, and geometric integrators are employed to guarantee fidelity to the analytic theory.

Applications extend to code verification, as in the benchmark efforts for ORB5 and GENE, enabling parametric β-scans and turbulence studies (Tronko et al., 2017). The geometric framework also enables large-eddy simulations (“α-models”) and energy–Casimir-based nonlinear stability analyses (Squire et al., 2013).

7. Metriplectic and Gauge-Free Extensions

Recent developments include:

• Metriplectic structure: A fully self-consistent metriplectic (Hamiltonian+metric) structure for collisional gyrokinetic Vlasov–Maxwell–Landau theory, with a new collisional current $D$9 ensuring local charge conservation,

$B$0

$B$1

(Hirvijoki et al., 2022).

• Gauge-free formulations: Using only field strengths (E, B) in the Hamiltonian and Poisson bracket achieves manifest gauge invariance and facilitates the construction of structure-preserving integrators (Brizard, 2021).

These extensions enhance the fidelity and physical transparency of gyrokinetic field theories and enable the rigorous simulation of kinetic plasma turbulence and transport.

---

Key References:

• Hamiltonian foundation and Poisson bracket: (Burby et al., 2014, Brizard et al., 2016, Brizard, 2021, Squire et al., 2013)
• Field-particle variational methods and models: (Brizard, 2018, Brizard, 2017, Tronko et al., 2017, Tronko et al., 2016, Brizard, 2019)
• Metric (dissipative) extension: (Hirvijoki et al., 2022)
• Gauge-free structure: (Brizard, 2021)
• Structure-preserving discretization: (Burby, 2016)
• Maximal-ordering reductions: (Zoni et al., 2019)