Double-Adiabatic Equations of State for Relativistic Plasmas

arXiv:2603.25669 · astro-ph.HE, physics.plasm-ph · Submitted 2026-03-26 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Next we'll be talking about the paper "Double-Adiabatic Equations of State for Relativistic Plasmas".

Jocelyn: The paper was written by A. Wierzchucka, P. J. Bilbao, A. G. R. Thomas, D. A. Uzdensky and A. A. Schekochihin from Rudolf Peierls Centre for Theoretical Physics, University of Oxford and Merton College, Oxford and Lady Margaret Hall, Oxford and Gérard Mourou Center for Ultrafast Optical Science, University of Michigan.

Vera: Stay tuned as we take you through the paper and discuss its implications.

Paper discussion segment 1: Vera: Jocelyn, I just opened this new preprint on arXiv called "Double-Adiabatic Equations of State for Relativistic Plasmas," and the title alone makes my head spin. It’s coming from a heavy-hitting group at Oxford and Michigan, including Agnieszka Wierzchucka and Dmitry Uzdensky.

Jocelyn: That sounds like some serious high-energy physics, Vera! When you say "double-adiabatic," are we talking about something that affects how we interpret the light coming from those distant sources?

Vera: It definitely does, because it’s looking at plasmas moving at nearly the speed of light. Most of our standard models assume these gases behave in a very simple, predictable way.

Jocelyn: So if they're saying the "equations of state" are different for relativistic plasmas, does that mean our current simulations might be missing something huge?

Subrahmanyan: You've hit on the exact problem, Jocelyn. In these extreme environments like pulsar magnetospheres or black hole accretion disks, the particles aren't just bumping into each other randomly like a normal gas. They are actually spiraling around magnetic field lines so fast that their motion becomes highly organized and anisotropic.

Vera: I see where you're going, Subrahmanyan, but how does that change what we actually see through a telescope?

Subrahmanyan: Well, if the pressure isn't the same in all directions—what they call anisotropy—the whole fluid behaves differently when it gets squeezed or stretched. This paper is trying to provide the mathematical rules for that specific, messy behavior in a relativistic setting.

Jocelyn: That sounds like a nightmare for someone trying to model a jet from an Active Galactic Nucleus! Are they claiming the old rules just don't apply at all?

Subrahmanyan: They aren't saying the old rules are wrong, but rather that they are incomplete for these high-energy regimes. If we use the wrong equations, our models of how these jets evolve might be fundamentally skewed.

Vera: It’s fascinating because it bridges the gap between simple fluid theory and the much more complex kinetic reality. We're going to look at what they actually found in their results next.

Jocelyn: I wonder if they actually found something concrete or just more math!

Subrahmanyan: They found something very concrete indeed, Jocelyn.

Paper discussion segment 2: Vera: We’ve established that "Double-Adiabatic Equations of State for Relativistic Plasmas" is tackling the math of how extreme gases behave. Now I'm looking at their summary, and they actually managed to derive exact functional forms for these pressures.

Jocelyn: So they aren't just making an educated guess? They actually found a way to link pressure directly to density and magnetic field strength?

Vera: They did, but it's much more complicated than the simple power laws we use in non-relativistic physics. In the standard non-relativistic case, you have these neat relations like pressure being proportional to density times the magnetic field.

Jocelyn: But for these relativistic ones, they say it depends on the anisotropy itself? That seems like a massive jump in complexity for a modeler.

Subrahmanyan: It is a huge jump, but it's physically necessary because of how the Lorentz factor works. In a relativistic plasma, the momentum integration involves that non-linear gamma factor, which means you can't just treat the parallel and perpendicular pressures as independent simple variables anymore.

Vera: I was reading their results section, and they found that if you have an ultra-relativistic plasma that is almost isotropic, the scaling actually changes to something like density to the power of four-fifths.

Jocelyn: That is a huge difference from the five-thirds or four-thirds we usually see! If I'm looking at a pulsar wind nebula, does this mean my estimates for its internal energy could be way off?

Subrahmanyan: Potentially, yes. They show that depending on whether the perpendicular pressure is much larger than the parallel pressure, or vice versa, you get entirely different scaling laws. For example, in one extreme case, they found a logarithmic correction that shows up when the parallel pressure dominates.

Vera: It's incredible how such a subtle change in particle motion leads to these totally different mathematical landscapes. Let's talk about how they actually proved this wasn't just theoretical math.

Jocelyn: Please tell me they did some actual testing!

Subrahmanyan: They certainly did.

Paper discussion segment 3: Vera: This is where it gets really impressive, because the authors didn't just stay at a chalkboard; they actually ran these massive Particle-in-Cell simulations to verify everything. They used a code called OSIRIS to simulate an electron-positron plasma being squeezed in a box.

Jocelyn: So they actually built a digital version of a relativistic plasma and watched it compress? How much did the simulation results match their new equations?

Vera: Jocelyn, it was nearly perfect! They showed that as the density increased due to compression, the true pressures followed those exact complex curves they derived.

Jocelyn: That must have been a massive computational undertaking. Did they run into any limits where their theory stopped working?

Subrahmanyan: They actually did, and that's one of the most honest parts of the paper. The double-adiabatic theory only works as long as certain symmetries hold, specifically when things like the first adiabatic invariant are conserved.

Vera: Right, they noted that once you hit a certain level of compression, these kinetic instabilities—like the mirror or firehose instabilities—start to kick in.

Jocelyn: And those instabilities must break all that beautiful math they just worked so hard on!

Subrahmanyan: Exactly. Once the plasma becomes unstable, it starts creating its own microscale magnetic fluctuations, which scatters the particles and destroys the very symmetry that makes these equations work. The simulation showed the pressure evolution departing from their theory right at that onset point.

Vera: It’s actually quite a relief to see them define exactly where the model is valid and where it breaks down. They've basically given us a map for when to use this new tool and when we need to switch to something much more expensive and complex.

Jocelyn: It sounds like they've provided a much-needed bridge between simple fluid models and full-scale kinetic simulations.

Subrahmanyan: It really is a vital contribution for anyone modeling high-energy environments.

Conclusion: Vera: We are running out of time, but I am so glad we got to discuss "Double-Adiabatic Equations of State for Relativistic Plasmas." This paper really feels like a major step forward for high-energy astrophysics.

Jocelyn: I agree, Vera. For anyone trying to model the most violent events in the universe, these new scaling laws are going to be essential tools. It’s going to change how we interpret the data from our next generation of telescopes.

Subrahmanyan: It really does provide a way to include relativistic effects in large-scale simulations without needing a supercomputer for every single step. It makes the whole field of relativistic magnetohydrodynamics much more robust.

Vera: We'll definitely be watching how this is used in future studies of magnetic reconnection and synchrotron cooling.

Jocelyn: Thanks for joining us, Subrahmanyan! You always help us see the big picture behind these equations.

Subrahmanyan: My pleasure; it's a thrilling time to be working on these high-energy problems.

Vera: Well, that's all for today! We'll catch you next time with another fascinating paper from the arXiv. Goodbye!

Jocelyn: Bye everyone! See you at the next one!

Rudolf Peierls Centre for Theoretical Physics, University of Oxford · Merton College, Oxford · Lady Margaret Hall, Oxford · Gérard Mourou Center for Ultrafast Optical Science, University of Michigan

astro-ph.HE, physics.plasm-ph

Submitted: 2026-03-26

Updated: 2026-09-16

Comments: 22 pages, 4 figures, accepted to JPP

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 88/100

The gist: The paper presents "a general first-principle formalism to derive adiabatic laws using the symmetries of the system," which recovers "the adiabatic equation of state P ∝ nΓ for isotropic plasmas

Terminology

Summary

The paper presents a general first-principle formalism to derive adiabatic laws using the symmetries of the system, which recovers the adiabatic equation of state P ∝ nΓ for isotropic plasmas and the double-adiabatic equations of state for collisionless, magnetised plasmas. The authors extend the latter to the relativistic regime, finding that their exact functional form depends on the pressure anisotropy and is not a simple power law.

Key aspects of the research include:

(

Derivation Methodology:

The authors abandon the standard moment-based derivation of the evolution equations and instead [formulate] a theory of adiabatic invariance focused on the system’s symmetries. By using phase-space volume conservation properties, they demonstrate that imposing certain symmetries leads to a self-similar evolution of the distribution function.

Findings for Isotropic Plasmas:

For isotropic plasmas, the study recovers the adiabatic equation of state P ∝ nΓ, where the adiabatic index Γ ranges from 5/3 in the non-relativistic regime to 4/3 in the ultra-relativistic limit. The authors note that in general, however, P depends not solely on density but also on the form of the initial distribution function f0 (p).

Findings for Relativistic Gyrotropic Plasmas:

For collisionless, magnetised plasmas—where the first adiabatic invariant, mu, is conserved—the authors derive new double-adiabatic equations of state. They find that for an ultra-relativistically hot plasma with an isotropic reference distribution, the exact form depends on the pressure anisotropy. The asymptotic limits are identified as:

(

1. Large perpendicular pressure (Δ ≫ 1):

P⊥ ∝ nB1/2 and P ∥ ∝ n3 /B5/2.

2. Almost isotropic plasma (Δ ≪ 1):

P⊥ ∝ (nB)4/5 and P ∥ ∝ (n3 /B2)4/5.

3. Large parallel pressure (Δ ≈ −1):

P⊥ ∝ B2 ln(n′2 /B′3) and P ∥ ∝ n3 /B5/2.

The authors note that these evolution equations agree with the initial discussion of the relativistic CGL equations in [35], up to the logarithmic correction in the second P⊥ expression, and they highlight that their results differ drastically from the non-relativistic theory [17], in which the same double-adiabatic relations hold for all anisotropies.

Validation and Applications:

The theoretical predictions were confirmed using two-dimensional particle-in-cell (PIC) simulations incorporating a large-scale compressive flow, where they found excellent agreement with our theoretical predictions. The authors state that these equations of state can be the EoS for relativistic collisionless plasmas and are "crucial for modelling large-scale high-energy astrophysical systems by means of a fluid theory, e.g., to determine regions where kinetic processes like the firehose and mirror instabilities are excited in high-beta plasmas. Furthermore, they suggest applicability to complex high-energy astrophysical processes, such as the evolution of plasmoid structures formed during magnetic reconnection and phenomena like synchrotron cooling."

Limitations:

The authors acknowledge that the distribution function is only gyrotropic if the coordinate r refers to the gyrocentre of the particle, which implies the results are valid only in the limit of vanishing gyroradius. Additionally, double adiabaticity requires neglect of heat fluxes and collisions, and their current derivation is restricted to the case of a plasma with a non-relativistic bulk flow.

Improvements for AI systems

To improve AI systems using the findings from this paper, we must move away from treating physical environments as static or purely statistical entities and instead implement models that respect the fundamental conservation laws of phase-space geometry.

The following are specific improvements for AI architectures, particularly those used in high-fidelity physical simulations (Physics-Informed Neural Networks - PINNs), autonomous control of plasma-based systems, and astrophysical modeling.

  1. Improvement: Transition from Scalar/Tensor Regression to Symmetry-Preserving Lagrangian Manifold Learning

Instead of training AI to predict pressure values directly (which often leads to unphysical energy gains or losses), the architecture should be redesigned to learn the evolution of the underlying distribution function through a symmetry-constrained latent space.

  • What the improved AI can do: In high-energy simulations, the AI will ensure that any predicted state maintains Liouville’s Theorem (phase-space volume conservation). This prevents numerical heating in long-term simulations of relativistic jets or accretion flows, allowing for stable, multi-million-hour virtual experiments without the simulation drifting into unphysical states.
  1. Improvement: Implementation of Anisotropy-Aware Physics-Informed Neural Networks (A-PINNs)

Current PINNs often use a single adiabatic index (e.g., 5/3 or 4/3). This paper proves that in relativistic regimes, the adiabatic index is not a constant but a dynamic function of the pressure anisotropy parameter:

While the AI is modeling a plasma, it will automatically adjust its internal stiffness (the relationship between density and pressure) based on the local magnetic field strength and parallel/perpendicular pressure ratios.

  • What the improved AI can do: It can accurately predict when a plasma system will undergo a phase transition or instability (like the Mirror or Firehose instabilities). For an AI controlling a fusion reactor or a plasma-based propulsion system, this means it can predict instability onset with much higher precision, enabling preemptive corrective actions before the physical system becomes unstable.
  1. Improvement: Non-Power-Law Adaptive Loss Functions for Relativistic Dynamics

Standard AI loss functions used in fluid dynamics often assume power-law relationships (e.g., scaling as a simple exponent). This paper demonstrates that relativistic EoS follow complex forms, such as logarithmic corrections in the parallel pressure limit:

The loss function of the AI should be reformulated to include these specific transcendental and logarithmic terms derived from the double-adiabatic equations.

  • What the improved AI can do: The AI will exhibit high-fidelity scaling. When simulating extreme astrophysical events (like magnetic reconnection in black hole environments), it will not oversimplify the physics. It will capture the subtle, non-linear pressure evolutions that current AI models miss, providing researchers with accurate data for interpreting telescope observations from events like blazars or pulsar wind nebulae.
  1. Improvement: Hybrid Kinetic-Fluid Surrogate Modeling

The paper provides a bridge between expensive kinetic (PIC) simulations and cheap fluid models. We can use this to create Surrogate Models that use the double-adiabatic EoS as a structural prior.

  • What the improved AI can do: It enables Real-Time Kinetic Intelligence. An AI could run at the speed of a fluid simulation but provide the accuracy of a particle-in-cell (PIC) simulation. This allows for real-time, high-fidelity digital twins of relativistic plasma environments, which is essential for designing next-generation high-energy physics experiments or analyzing massive datasets from cosmic ray observatories.

Abstract

The adiabatic equation of state P proportional to n Γ describes the pressure evolution of highly collisional, isotropic plasmas in terms of their density, providing a possible closure of the fluid moment hierarchy in the absence of heat fluxes and dissipation. An analogous closure exists for collisionless, magnetised plasmas, whose pressure tensor is anisotropic with respect to the magnetic field, and the closure is therefore double adiabatic, prescribing the evolution of the parallel and perpendicular pressures in terms of the magnetic-field strength and density. Here, we present a general first-principle formalism to derive adiabatic laws using the symmetries of the system. With this theory we recover the adiabatic equation of state P proportional to n Γ for isotropic plasmas and the double-adiabatic equations of state for collisionless, magnetised plasmas. We extend the latter to the relativistic regime, finding that their exact functional form depends on the pressure anisotropy and is not a simple power law. Our double-adiabatic equations of state describe simple geometries, like magnetic mirrors or compressed homogeneous plasmas, as well as complex high-energy astrophysical processes, such as the evolution of plasmoid structures formed during magnetic reconnection.

Related papers