---
title: Gyrokinetic Vlasov–Maxwell System Overview
url: https://www.emergentmind.com/topics/gyrokinetic-vlasov-maxwell-system
type: topic
---

# Gyrokinetic Vlasov–Maxwell System Overview

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 \)) [1411.1790]. 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 \( P_{\rm gy} = -\delta {\mathcal K}/\delta E \) 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 [1411.1790, 1606.06652, 1301.6066, 1702.04747, 2108.09384, 1907.11204].

## 2. Noncanonical Poisson Bracket and Casimir Invariants

The dynamic evolution is governed by a noncanonical Poisson bracket. For functionals \( F \), \( G \) of \( (f, D, B) \),
\[
\{F, G\}_{\rm GVM} = \sum_{s=1}^{N_s} \int_{TQ} B_{\rm gy}\left( d\left(\frac{\delta F}{\delta f_s}\right) - \frac{4\pi e_s}{c} T^*\left(\frac{\delta F}{\delta D}\right), d\left(\frac{\delta G}{\delta f_s}\right) - \frac{4\pi e_s}{c} T^*\left(\frac{\delta G}{\delta D}\right) \right) f_s \\
+ 4\pi c \int_Q \left( \frac{\delta F}{\delta D} \wedge * \frac{\delta G}{\delta B} - \frac{\delta G}{\delta D} \wedge * \frac{\delta F}{\delta B} \right),
\]
where \( B_{\rm gy} \) is defined by the gyrocenter symplectic form, and \( T^* \) is the pullback from Q to \( TQ \) [1411.1790, 1606.06652].

Casimir invariants take the form:
- **Particle-label Casimirs:** \( C_h[f] = \sum_s \int_{TQ} h(F_s) \Omega_s \), with \( h \) arbitrary.
- **Field Casimirs:** Any functional of \( \delta D + 4\pi P_{\rm gy} \) is a Casimir, corresponding to the conservation of Gauss’s law \( \nabla \cdot D - 4\pi\rho_{\rm gy} \) [1411.1790, 1301.6066].

The bracket is verified to be antisymmetric and to satisfy the Jacobi identity, ensuring a bona fide Hamiltonian system [2108.09384, 1606.06652].

## 3. Equations of Motion and Electromagnetic Closure

The evolution equations consist of:
- **Gyrokinetic Vlasov equation:** For each \( s \),
  \[
  \partial_t f_s + \{f_s, H_{\rm GVM}\} = 0,
  \]
  or, in the characteristic form in phase space (gyrocenter or guiding-center coordinates),
  \[
  \frac{\partial F_s}{\partial t} + \nabla \cdot (F_s \dot{\mathbf X}) + \partial_{p_\parallel}(F_s \dot p_\parallel) = 0,
  \]
  where trajectories are generated by the Hamiltonian via the noncanonical Poisson bracket [1411.1790, 1604.03538, 1907.11204].

- **Maxwell's equations (with gyrokinetic constitutive relations):**
  \[
  \nabla \cdot D = 4\pi \sum_s \int dP\, q_s\, F_s,
  \]
  \[
  \nabla \times H - \frac{1}{c} \frac{\partial D}{\partial t} = \frac{4\pi}{c}(J_{\rm gc} + J_{\rm coll}),
  \]
  where the polarization and magnetization terms are systematically derived from the gyrocenter Hamiltonian, and \( J_{\rm coll} \) arises from collisional metriplectic extensions [2203.08464].

The polarization and magnetization enter Maxwell’s equations via
\[
\mathbf D = \mathbf E_1 + 4\pi \mathbf P, \qquad \mathbf H = \mathbf B_1 - 4\pi \mathbf M,
\]
where explicit expressions are given, typically as moments over the gyrocenter distribution [1703.07582, 1807.04367].

## 4. Energy, Momentum, and Entropy Conservation

The system possesses an exact energy conservation law, derivable by Noether’s theorem from the gyrokinetic action:
\[
\frac{d}{dt} \int d^3x\, \mathcal E = 0, \qquad \mathcal{E} = \sum_s \int d^6Z\,F\,H_{\rm gy} + \frac{1}{8\pi}\left(|\nabla\Phi_1|^2 + |\mathbf{B}_0+\epsilon\nabla\times\mathbf{A}_1|^2\right),
\]
valid in both the collisionless [1411.1790, 1702.04747, 1604.03538] and metriplectic collisional [2203.08464] frameworks.

Canonical and toroidal angular momentum are also exactly conserved in an axisymmetric background:
\[
\frac{d\mathcal P_\varphi }{dt} = 0,
\]
where \( \mathcal P_\varphi \) includes both gyrocenter and field angular momentum [1907.11204, 2108.09384, 1604.03538, 1702.04747].

In metriplectic (collisional) frameworks, the entropy functional,
\[
\mathcal S[F] = -\sum_s \int dZ\, F_s \ln(F_s/B_{\parallel s}^*),
\]
is a Casimir of the Poisson bracket and is strictly non-decreasing under the metric (dissipative) bracket, ensuring monotonic entropy production [2203.08464].

## 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 \( \epsilon = \rho/L \) [1604.03538, 1909.05036]. This reduction yields a gyrokinetic Hamiltonian
\[
H_{\rm gy} = H_0 + \epsilon H_1 + \epsilon^2 H_2 + \ldots,
\]
with \( H_0 = p_\parallel^2/2m + \mu B_0 \), \( H_1 \) encoding field perturbations (e.g., \( e\langle \Phi_1 \rangle \)), and \( H_2 \) containing polarization and ponderomotive terms [1702.04747, 1807.04367].

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) [1703.07582]. Ion and electron dynamics are treated consistently, allowing for variations such as maximal ordering for fusion devices [1909.05036].

## 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 [1611.03064, 1301.6066]. 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 [1703.07582]. The geometric framework also enables large-eddy simulations (“α-models”) and energy–Casimir-based nonlinear stability analyses [1301.6066].

## 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 \( \mathbf J_{\rm coll} \) ensuring local charge conservation,
  \[
  \partial_t F_s + \{F_s, H\} + C[F_s] = 0,
  \]
  \[
  \frac{d\mathcal H}{dt} = 0, \qquad \frac{d\mathcal S}{dt} \ge 0
  \]
  [2203.08464].
- **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 [2108.09384].

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: [1411.1790], [1606.06652], [2108.09384], [1301.6066]  
- Field-particle variational methods and models: [1807.04367], [1702.04747], [1703.07582], [1604.03538], [1907.11204]  
- Metric (dissipative) extension: [2203.08464]  
- Gauge-free structure: [2108.09384]  
- Structure-preserving discretization: [1611.03064]  
- Maximal-ordering reductions: [1909.05036]

Source: https://www.emergentmind.com/topics/gyrokinetic-vlasov-maxwell-system