Double-Adiabatic Equations of State for Relativistic Plasmas

This paper presents a general first-principle formalism based on system symmetries to derive adiabatic equations of state for both isotropic and anisotropic plasmas, specifically extending the double-adiabatic closure for collisionless, magnetised plasmas to the relativistic regime where the exact functional form depends on pressure anisotropy rather than following a simple power law.

The Gist

Imagine you are trying to predict how a crowd of people moves through a city. If everyone is bumping into each other constantly (like a high-pressure, crowded market), they tend to move as a single, smooth fluid. You can predict their behavior with simple rules: if you squeeze the crowd into a smaller space, the pressure goes up in a predictable way. In physics, this is called an adiabatic equation of state.

But now, imagine a different scenario: a massive, empty stadium where people are running around on their own, rarely touching, but all following the same invisible magnetic "tracks" laid out on the floor. This is a collisionless plasma, like the stuff found in space near black holes or in the solar wind. Here, the simple rules break down because the people (particles) aren't bumping into each other; they are spiraling around the magnetic tracks.

This paper is about figuring out the new, more complex rules for how this "space crowd" behaves when they are moving at speeds close to the speed of light (relativistic).

The Old Rules vs. The New Reality

The Old Rules (Non-Relativistic):
For decades, physicists have used a set of rules called the CGL equations (named after Chew, Goldberger, and Low). Think of these as a "Double-Adiabatic" law.

• The Analogy: Imagine the particles are like spinning tops.
  • Rule 1: If you squeeze the magnetic field lines closer together, the tops spin faster (perpendicular pressure goes up).
  • Rule 2: If you stretch the space they are in, they slow down their forward motion (parallel pressure changes differently).
• In the slow-speed world, these rules are simple power laws. It's like saying, "If you double the density, the pressure goes up by exactly 8 times." Simple and clean.

The Problem:
When these particles move near the speed of light, things get weird. Their mass effectively increases, and their energy behaves differently. The old "simple power law" rules stop working. The relationship between density, magnetic field, and pressure becomes messy and depends on how the particles are distributed.

The Paper's Big Idea: Symmetry as a Compass

The authors (Wierzchucka, Bilbao, et al.) didn't just try to brute-force the math. Instead, they looked at the symmetries of the system.

• The Metaphor: Imagine a spinning dancer. No matter how fast she spins, her shape looks the same from the side (symmetry). The authors realized that even in the chaotic relativistic world, the plasma particles still have these hidden symmetries. They are "gyrotropic" (spinning symmetrically around magnetic lines) and "parity symmetric" (looking the same if you flip them forward or backward).
• By focusing on these symmetries and the fact that the "volume" of the phase space (a fancy way of saying the total space of all possible positions and speeds) is conserved, they derived a general formula for how the pressure evolves.

The Surprising Discovery: It's Not Just a Simple Formula

When they applied this new theory to ultra-fast (ultra-relativistic) plasmas, they found something shocking: There is no single simple formula.

In the old, slow world, the rule was always the same. In the relativistic world, the rule changes depending on the balance of the pressure:

1. If the particles are mostly spinning sideways (Perpendicular pressure is huge): The pressure grows in one specific way.
2. If the particles are mostly shooting forward (Parallel pressure is huge): The pressure grows in a completely different way, involving logarithms (a mathematical curve that grows slowly).
3. If the particles are balanced (Isotropic): It follows a middle-ground rule.

The Analogy:
Think of the old rule as a recipe that says, "Add 2 cups of flour for every cup of water."
The new relativistic rule says, "Add 2 cups of flour if the water is cold, but add 3 cups if the water is hot, and if the water is boiling, you need to add a pinch of salt and wait 5 minutes." The recipe depends entirely on the current state of the mixture.

Why Does This Matter?

Why should a general audience care about the pressure of invisible space particles?

1. Understanding the Universe's Engines: Many of the most energetic events in the universe—like black hole jets, pulsar winds, and magnetic reconnection (where magnetic field lines snap and reconnect, releasing massive energy)—happen in this relativistic, collisionless regime. To simulate these on a computer, scientists need a "rulebook" (an Equation of State) to tell the computer how the plasma reacts. This paper provides that rulebook.
2. Predicting Instabilities: Just like a stretched rubber band can snap, plasmas can become unstable. If the pressure gets too lopsided, the plasma explodes into turbulence (instabilities like the "firehose" or "mirror" instability). The authors found that in the relativistic world, these plasmas might be more stable than we thought. It takes much more squeezing to trigger these explosions in fast-moving plasmas than in slow ones.
3. New Physics: They confirmed their theory using supercomputer simulations (Particle-in-Cell simulations), which act like a digital wind tunnel. The results matched their new, complex formulas perfectly, proving that the universe really does follow these tricky, state-dependent rules.

The Takeaway

This paper is a bridge between the simple, intuitive physics of slow-moving gases and the chaotic, high-speed physics of the relativistic universe.

The authors showed that when particles move near the speed of light, the "double-adiabatic" laws aren't simple, one-size-fits-all equations. Instead, they are shape-shifters. The way the plasma pressure evolves depends on whether the particles are spinning wildly or shooting straight ahead. By understanding these nuances, we can finally build better models to explain how the most energetic objects in the cosmos work, from the jets of black holes to the radiation belts around Earth.

Technical Summary

1. Problem Statement

The paper addresses the challenge of modeling collisionless, magnetized plasmas in the relativistic regime.

• Context: In non-relativistic, collisionless plasmas, the pressure tensor is anisotropic (different parallel $P_\parallel$ and perpendicular $P_\perp$ pressures relative to the magnetic field). The evolution of these pressures is traditionally described by the Chew-Goldberger-Low (CGL) equations (the "double-adiabatic" approximation), which rely on the conservation of two adiabatic invariants: the magnetic moment ($\mu$) and the longitudinal invariant ($J$).
• The Gap: While the non-relativistic CGL equations provide a simple power-law closure ($P_\perp \propto nB$, $P_\parallel \propto n^3/B^2$), their relativistic generalization is non-trivial. In the relativistic regime, the definition of pressure involves the Lorentz factor $\gamma$, which is a non-linear function of both parallel and perpendicular momenta. Consequently, $P_\perp$ and $P_\parallel$ depend on the full momentum distribution, not just their respective components.
• Limitation of Previous Work: Existing relativistic extensions (e.g., modified pressures) often fail to provide a sufficient fluid closure because they evolve "modified" pressures rather than the physical pressures required for the fluid momentum equation. Furthermore, the exact functional form of the relativistic equation of state (EoS) was previously unknown, particularly regarding its dependence on pressure anisotropy.

2. Methodology

The authors developed a first-principles formalism based on phase-space symmetries rather than standard moment expansions of the Vlasov equation.

• Symmetry-Based Derivation:
  • They treat the distribution function $f$ as the density of an incompressible fluid in 6D phase space.
  • They impose specific symmetries on the distribution function in the "peculiar momentum" frame (the frame moving with the bulk flow):
    1. Gyrotropy: Cylindrical symmetry around the magnetic field vector $\mathbf{B}$.
    2. Parity Symmetry: Symmetry along the direction of $\mathbf{B}$.
  • They utilize Liouville's theorem (conservation of phase-space volume) and the conservation of the first adiabatic invariant ($\mu$, implying conservation of perpendicular phase-space area) to derive how the distribution function evolves self-similarly under compression/expansion.
• Derivation of EoS:
  • By mapping the initial distribution $f_0$ to the evolved distribution $f(t)$ using the derived scaling laws for momentum components ($\tilde{p}_\perp$ and $\tilde{p}_\parallel$), they calculate the pressure moments.
  • They specifically analyze the ultra-relativistic limit ($p \gg mc$) assuming an initially isotropic distribution.
• Numerical Validation:
  • The theoretical predictions were tested using 2D Particle-in-Cell (PIC) simulations using the OSIRIS code.
  • Simulations involved electron-positron plasmas compressed both perpendicular and parallel to a uniform magnetic field to induce varying degrees of pressure anisotropy.

3. Key Contributions

• General Formalism: Established a symmetry-based method to derive adiabatic laws that applies to both isotropic and anisotropic plasmas, recovering known non-relativistic limits and extending them to the relativistic regime.
• Exact Relativistic EoS: Derived the exact functional form of the double-adiabatic equations for ultra-relativistic plasmas. Crucially, they demonstrated that the EoS is not a simple power law in the general case; it depends explicitly on the pressure anisotropy.
• Asymptotic Limits: Provided analytical solutions for three distinct regimes:
  1. Near Isotropy ($P_\perp \approx P_\parallel$): $P_\perp \propto (nB)^{4/5}$ and $P_\parallel \propto (n^3/B^2)^{4/5}$.
  2. Strong Perpendicular Anisotropy ($P_\perp \gg P_\parallel$): $P_\perp \propto n B^{1/2}$ and $P_\parallel \propto n^3 / B^{5/2}$.
  3. Strong Parallel Anisotropy ($P_\parallel \gg P_\perp$): $P_\perp \propto B^2 \ln(n^2/B^3)$ and $P_\parallel \propto n^3 / B^{5/2}$.
• Generalized Adiabatic Indices: Introduced generalized adiabatic indices ($\Gamma_\perp, \Gamma_\parallel$) that vary continuously with the anisotropy parameter $A = B'^3/n'^2$, allowing for a fluid-like closure that adapts to the plasma state.
• Eulerian Formulation: Reformulated the Lagrangian results into an Eulerian form suitable for implementation in numerical fluid simulations, including non-ideal effects like cooling.

4. Results

• Theoretical Agreement: The derived equations (Eqs. 51 and 52) successfully describe the evolution of $P_\perp$ and $P_\parallel$ in the absence of heat fluxes. The results differ drastically from non-relativistic theory, where the power laws are constant regardless of anisotropy.
• Simulation Validation:
  • PIC simulations of compressing boxes showed excellent agreement with the theoretical predictions for $P_\perp$ and $P_\parallel$ as functions of density $n'$ and magnetic field $B'$.
  • The simulations confirmed the transition from the "near-isotropic" scaling ($4/5$ power) to the "highly anisotropic" scalings as the compression proceeded.
• Instability Onset: The simulations demonstrated the breakdown of the double-adiabatic approximation when kinetic instabilities (Mirror instability for perpendicular compression, Firehose instability for parallel compression) are triggered. This occurs when the pressure anisotropy exceeds specific thresholds, leading to magnetic energy growth and scattering that violates the conservation of adiabatic invariants.
• Stability Implications: The results suggest that relativistic plasmas may be more stable against anisotropy-driven instabilities than non-relativistic plasmas. For instance, under perpendicular compression, the relativistic $\beta_\perp$ decreases, requiring significantly larger compressions to trigger the mirror instability.

5. Significance

• Fluid Modeling of High-Energy Astrophysics: The derived EoS provides a crucial "closure" for fluid models of relativistic, collisionless plasmas. This enables the simulation of large-scale astrophysical systems (e.g., pulsar wind nebulae, black hole accretion flows, AGN jets) without resorting to computationally expensive kinetic simulations for the entire system.
• Magnetic Reconnection: The equations offer a model for particle energization within plasmoids formed during relativistic magnetic reconnection, a key mechanism for non-thermal particle acceleration.
• Synchrotron Cooling: The framework can be coupled with synchrotron cooling processes to model complex phenomena like the synchrotron firehose instability, which creates two-phase media in cooling plasmas.
• Theoretical Foundation: By moving away from moment-based derivations to a symmetry-based approach, the paper clarifies the fundamental conditions (gyrotropy, parity, phase-space conservation) required for double-adiabatic behavior, making the theory robust and applicable to various systems with frozen-in vector fields.

In summary, this work bridges the gap between kinetic theory and fluid modeling for relativistic plasmas, providing a rigorous, symmetry-derived equation of state that captures the complex, anisotropy-dependent evolution of pressure in high-energy astrophysical environments.