Equilibria for the Vlasov-Maxwell system related to plasma confinement

Technical Summary: Equilibria for the Vlasov-Maxwell System Related to Plasma Confinement

Problem Statement
The paper investigates the existence of semi-explicit, physically relevant equilibria for the two-species Vlasov-Maxwell (VM) system. The study focuses on a cylindrically and vertically symmetric setting ($x, y \in \mathbb{R}$, $z \in \mathbb{T}$), where the particle densities $f^\pm$, electric field $E$, and magnetic field $B$ exhibit specific invariance properties. The primary goal is to construct electrically neutral ($E=0, \rho=0$), magnetically driven equilibria that correspond to three canonical plasma confinement configurations: $z$-pinch, $\theta$-pinch, and screw-pinch. Unlike previous works that often rely on fluid closure assumptions, this study utilizes a kinetic description where equilibria are determined self-consistently from particle distributions and electromagnetic fields, making the VM setting suitable for rigorous stability analysis.

Methodology
The authors reduce the existence problem of these equilibria to a set of second-order ordinary differential equations (ODEs) in the radial variable $r$. The methodology proceeds as follows:

1. Symmetry Reduction: By assuming invariance in the azimuthal ($\theta$) and vertical ($z$) directions, the VM system simplifies. The densities $f^\pm$ are expressed as functions of the particle trajectory invariants: energy $e^\pm$, the $z$-component of linear momentum $p_z^\pm$, and the $z$-component of angular momentum $\ell_z^\pm$.
2. Ansatz Construction: The authors propose specific functional forms for the distribution functions $\mu^\pm$ (typically exponential functions of the invariants) that ensure electrical neutrality ($\rho=0$).
3. Field Reduction: Substituting these densities into the Maxwell equations (specifically Ampère's law) reduces the field equations to nonlinear second-order ODEs for the magnetic vector potential components ($A_z$, $A_\theta$, or a hybrid $A_\eta$).
4. Dimensionless Scaling: Physical constants (masses $m^\pm$, charges $q^\pm$, thermal speeds $\sigma^\pm$, etc.) are grouped into dimensionless parameters (e.g., $d$, $\lambda$, $a$) to streamline the ODEs.
5. Analysis of ODEs: The resulting ODEs are analyzed for existence, boundary conditions (regularity at $r=0$ and decay at $r \to \infty$), and asymptotic behavior. In some cases, explicit solutions are found; in others, numerical simulations and qualitative mechanical energy analysis are employed.

Key Contributions and Results

The paper presents four main theorems corresponding to different equilibrium types:

• Theorem 1.1 ($z$-pinch, constant drift):

  • Configuration: Magnetic field points in the $\theta$-direction; current flows in the $z$-direction.
  • Distribution: $\mu^\pm$ depends on $e^\pm$ and $p_z^\pm$ with a linear term in $p_z^\pm$, resulting in a constant drift velocity $\pm \sigma^\pm$ in the $z$-direction.
  • Result: The field equation reduces to $(-\partial_{\bar{r}}^2 - \bar{r}^{-1}\partial_{\bar{r}})\bar{A}_z = d e^{\bar{A}_z}$. For a specific choice of the dimensionless constant $d=8$, an explicit solution (the Bennett solution) exists: $\bar{A}_z = -2\log(1+e^D\bar{r}^2) + D$.
  • Decay: The density decays as $r^{-4}$, and the magnetic field decays as $r^{-1}$. The authors note that while the current vanishes at infinity, the magnetic field fails to be in $L^2(\mathbb{R}^3)$, consistent with the non-existence of nontrivial finite-energy equilibria in $\mathbb{R}^3$ (citing [GS84]).

• Theorem 1.2 ($z$-pinch, variable drift):

  • Configuration: Similar to 1.1 but with a quadratic dependence on $p_z^\pm$ in the exponent, leading to a variable drift velocity dependent on $r$.
  • Result: The field equation becomes $(-\partial_{\bar{r}}^2 - \bar{r}^{-1}\partial_{\bar{r}})\bar{A}_z = -d \bar{A}_z e^{-\bar{A}_z^2/2}$.
  • Analysis: Mechanical energy analysis shows that solutions satisfying regularity at $r=0$ cannot decay exponentially to zero at infinity. Instead, $\bar{A}_z$ diverges logarithmically ($\sim \log \bar{r}$) as $\bar{r} \to \infty$.
  • Decay: This logarithmic divergence leads to a "super-power" decay rate for the current and number densities (faster than any polynomial power), specifically $n^\pm \sim \bar{r}^{-C^2 \log \bar{r}}$ (with $C \approx 1.67$ in simulations). The drift velocity, given by $\mp \bar{A}_z(\bar{r})\sigma_\pm/2$, diverges as $r \to \infty$. This contrasts with Theorem 1.1, where the drift velocity is constant.

• Theorem 1.3 ($\theta$-pinch):

  • Configuration: Magnetic field points in the $z$-direction; current flows in the $\theta$-direction.
  • Distribution: $\mu^\pm$ depends on $e^\pm$ and $\ell_z^\pm$.
  • Result: The field equation for $\bar{A}_\theta$ is $(-\partial_{\bar{r}}^2 - \bar{r}^{-1}\partial_{\bar{r}} + \bar{r}^{-2})\bar{A}_\theta = -d \frac{\bar{r}^2 \bar{A}_\theta}{(1+\bar{r}^2)^{3/2}} \exp\left[-\frac{\bar{r}^2 \bar{A}_\theta^2}{2(1+\bar{r}^2)}\right]$.
  • Decay: Solutions satisfying boundary conditions behave as $\bar{A}_\theta \sim \bar{r}$ as $\bar{r} \to \infty$. Consequently, the magnetic field $B_z$ approaches a non-zero constant, but the particle densities and current decay at a near-Gaussian rate ($\sim e^{-C^2 \bar{r}^2/2}$), indicating strong confinement despite the non-decaying magnetic field.

• Theorem 1.4 (Screw-pinch):

  • Configuration: A generalization interpolating between $z$- and $\theta$-pinches using a pitch parameter $k$. The current flows helically.
  • Method: The authors introduce a rotated basis $(e_\eta, e_\sigma)$ aligned with the helical symmetry. The distribution depends on the invariant $\ell_z^\pm + k p_z^\pm$.
  • Result: The field equation reduces to a single ODE for $\bar{A}_\eta$.
  • Decay: Similar to the $\theta$-pinch, the magnetic field does not decay to zero, but the densities and current exhibit Gaussian-like decay, ensuring strong localization of particles.

Extensions Discussed

• External Fields (Section 8): The authors discuss how including an external field allows for solutions where the magnetic field converges exponentially to zero, addressing the finite-energy limitation of the self-contained cases.
• Relativistic Case (Section 9): The paper outlines how the non-relativistic results can be adapted to the relativistic case by modifying the energy term in the distribution function, though the resulting ODEs would need to be re-derived.

Significance and Claims
The paper claims to provide the first construction of semi-explicit, electrically neutral Vlasov-Maxwell equilibria for $z$-pinch, $\theta$-pinch, and screw-pinch configurations within a cylindrically and vertically symmetric geometry. The authors emphasize that these solutions are "semi-explicit": while the distribution functions are given in closed form, the magnetic potentials are solutions to nonlinear ODEs that generally require numerical calculation or qualitative analysis, though explicit solutions exist for specific parameter choices (e.g., the Bennett solution).

The primary significance of this work, as stated by the authors, is to provide a foundation for future rigorous stability analysis. They intend to use these specific equilibria to test linear stability and instability frameworks developed by Lin & Strauss [LS07, LS08] and Lin [Lin22], which reduce stability questions to spectral conditions of associated scalar problems. The authors explicitly note that the stability of these specific semi-explicit equilibria has not yet been thoroughly investigated and that their current work is a necessary precursor to such stability studies. They also acknowledge that obtaining more complex equilibria depending on multiple variables would require more abstract methods (e.g., bifurcation theory), positioning their radial solutions as a starting point for such analyses.