On the double-adiabatic equations in the relativistic regime
summary
The gist
The gist: This work analytically solves the drift kinetic equation for a homogeneous, magnetized, collisionless plasma and obtains general solutions for time-dependent, gyrotropic distribution
In short
This work analytically solves relativistic drift kinetic equations to extend double-adiabatic evolution equations beyond nonrelativistic limits. It derives new, time-dependent pressure evolution expressions for perpendicular and parallel pressures in relativistic and ultrarelativistic plasmas, providing a general solution applicable to astrophysical phenomena like pulsar winds.
Key concepts
- Drift Kinetic Equation
- This is the fundamental equation describing how the distribution of particles changes over time due to their motion (drift) within a magnetized plasma. It tracks the evolution of the particle distribution function f(t, p⊥, p∥), which depends on momentum components perpendicular and parallel to the magnetic field.
- Double-Adiabatic Equations
- These are equations that describe how two key thermodynamic quantities—perpendicular pressure (P⊥) and parallel pressure (P∥)—evolve independently over time. The paper extends these classic equations to account for relativistic effects, providing a more accurate description of plasma behavior.
- Maxwell-Jüttner Distribution
- This is a specific type of distribution function used to describe the momentum states of particles in a relativistic plasma. It is an extension of the nonrelativistic Maxwellian distribution and is crucial for modeling ultrarelativistic conditions where particle speeds approach the speed of light.
Terminology used across episodes
This episode discusses
- On the double-adiabatic equations in the relativistic regime · Paper Radio
- Double-Adiabatic Equations of State for Relativistic Plasmas · Paper Radio
The paper
On the double-adiabatic equations in the relativistic regime · Read on arXiv
Departamento de Física, Facultad de Ciencias Físicas y Matemıcas, Universidad de Chile · Department of Astronomy, University of Wisconsin--Madison · Department of Physics, University of Wisconsin--Madison
We revisit the double adiabatic evolution equations and extend them to the relativistic and ultrarelativistic regimes. We analytically solve the relativistic, time-dependent drift kinetic equation for a homogeneous, magnetized, collisionless plasma and obtain a solution explicitly dependent on the magnetic field and density variations. In the case of an initial relativistic Maxwellian distribution, a natural extension to an anisotropic Maxwell-Jüttner is obtained. We calculate the moments of this time-dependent solution and obtain analytical expressions for the evolution of the perpendicular and parallel pressures in the ultrarelativistic case. We numerically solve the moment equations in the relativistic case and obtain general expressions for the double-adiabatic equations in this regime. We confirm our results using fully kinetic particle-in-cell simulations of shearing and compressing boxes. Our findings can be readily applied to relativistic species including cosmic-rays and electron-positron pairs, present in astrophysical plasmas like pulsar wind nebulae, astrophysical jets, black hole accretion flows, and Van Allen radiation belts.
DOI: 10.1017/S0022377826102220
Transcript
Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.
Vera: Today's paper: "On the double-adiabatic equations in the relativistic regime".
Jocelyn: The gist: This work analytically solves the drift kinetic equation for a homogeneous, magnetized, collisionless plasma and obtains general solutions for time-dependent, gyrotropic distribution functions in nonrelativistic, relativistic,
Vera: First, who's behind it and why it matters.
Title and authors: Vera: To summarize what they are doing here in "On the double-adiabatic equations in the relativistic regime," they analytically solve the drift kinetic equation to get a time-dependent distribution function f alpha t p⊥ p∥, and then they calculate the moments of this solution to get double-adiabatic evolution equations for P⊥ and P∥.
Jocelyn: They are focusing on solving that drift kinetic equation in three different regimes: nonrelativistic, relativistic, and ultrarelativistic initial distributions. This shows how the physics handles changes from normal speeds up to speeds much closer to the speed of light.
Subrahmanyan: The summary highlights that for a nonrelativistic Maxwellian initial distribution, they recover the classic CGL double-adiabatic equations for T⊥ t and T∥ t. This is just checking if their general method works when things are slow.
Vera: Then they move to the relativistic extension, where they derive new evolution equations for P⊥ and P∥ that depend only on density n(t) and field strength B(t), which is a big step because it decouples the pressure evolution.
Jocelyn: And when they get to the ultrarelativistic case, they obtain analytical expressions for the perpendicular pressure P⊥ (two point one eight) and parallel pressure P∥ (two point two one) that are tractable when kBT over mc squared is greater than one times the speed of light squared, which is what makes those high-speed calculations possible <ref:2603.25594#pg1,obtain analytical expressions for the>.
Subrahmanyan: The paper shows that these analytical expressions for the ultrarelativistic case are also consistent with the general state equations for a relativistic anisotropic plasma from Gedalin in one thousand nine hundred ninety-one showing they hold up across different theoretical frameworks <ref:2603.25594#pg2>.
Vera: So, essentially, they are taking a distribution function evolution and turning it into simple equations for pressure components based on how fast the density and field are changing.
Jocelyn: And this is important because it gives us a way to track how anisotropy develops in these extreme environments without having to solve the full kinetic equation every time.
The paper's summary: Vera: One of the main improvements they propose is that they can generate new, relativistic expressions for the evolution of both perpendicular and parallel pressures by numerically integrating those analytical solutions derived from the drift kinetic equation, like equations two point three three and two point three five.
Jocelyn: They’re essentially using those analytical results to create a new way to get double-adiabatic evolution for the relativistic and ultrarelativistic cases that builds on the classical CGL model by giving it relativistic versions.
Subrahmanyan: This means they've managed to extend the CGL equations, which are useful for constructing fluid models with anisotropic closures to MHD equations, into a regime where particle speeds are significant fractions of light speed.
Vera: They also show that their derived distribution function f t p p is the analytical form of the distribution function that Particle-in-Cell simulations evolve when initialized with a Maxwell-Jüttner distribution, until some instability kicks in.
Jocelyn: So, they’ve connected the abstract kinetic solutions directly to what we see in simulations, which gives us a concrete way to predict how the plasma actually behaves before it breaks down.
Subrahmanyan: This connection is powerful because it suggests that their analytical form for P⊥ (two point three three) and P∥ (two point three five) correctly describes the evolution of P⊥ i and P∥ i in shearing simulations for mildly relativistic temperatures, which is a very specific test case.
Vera: And they also mentioned that this general solution works even if you start with other initial distributions, like power-law or kappa distributions, which are relevant for things like cosmic rays and solar wind studies.
The paper's improvements: Jocelyn: So to wrap up the paper "On the double-adiabatic equations in the relativistic regime," they show that their derived distribution function is a general solution for the time-dependent drift kinetic equation, which allows them to predict pressure evolution based on density and field variations.
Vera: They provide new, relativistic expressions for P⊥ and P∥ when starting with an initial relativistic Maxwellian distribution, and these results are validated by Particle-in-Cell simulations of both shearing and compressing boxes in both mildly relativistic and ultrarelativistic initial temperatures.
Subrahmanyan: The implication is that the relativistic double-adiabatic equations now provide a new closure for extended magnetohydrodynamic MHD theories in the context of collisionless plasmas where heat fluxes are negligible and spatial gradients are large, which is relevant for understanding pulsar wind nebulae or black hole accretion flows.
Jocelyn: They also confirm that their analytical expressions for P⊥ in the ultrarelativistic case are consistent with those from the relativistic moment equations, showing a robust connection across different theoretical approaches.
Vera: So, if you’re looking at how plasma pressure anisotropy evolves in these extreme settings, this paper provides a new tool—the double-adiabatic equations in the relativistic regime—that can be used to model things like jets and radiation belts.
Subrahmanyan: I just think the work confirms that the specific form of these double-adiabatic equations depends on the initial distribution function used, which is a key point for anyone trying to apply this theory to real cosmic ray physics.
Jocelyn: It’s a solid piece of work connecting kinetic theory, simulations, and general fluid models for anisotropic plasmas.
Vera: We're all ready for whatever comes next on the arXiv.
Conclusion: Vera: So we’ve been looking at "On the double-adiabatic equations in the relativistic regime," and essentially, this paper solves the drift kinetic equation to get a general solution for a time-dependent distribution function, which then lets them derive new evolution equations for perpendicular and parallel pressures.
Jocelyn: Right. They’re taking that kinetic equation—the drift kinetic equation—and showing how you can get those double-adiabatic equations for P⊥ and P∥ in the relativistic and ultrarelativistic regimes based only on the density and magnetic field changes.
Subrahmanyan: The big picture here is that this gives us a new way to close extended magnetohydrodynamic theories when you have collisionless plasmas where heat fluxes are negligible but spatial gradients are large, which is super relevant for modeling astrophysical phenomena.
Vera: And they do connect this to simulations, you know? They use Particle-in-Cell simulations of shearing and compressing boxes to confirm that their analytical expressions for P⊥ and P∥ actually describe the pressure evolution seen in those runs across different initial temperatures.
Jocelyn: That’s important because it means we can take those empirical histograms from the simulations, like the faMJ(t, p⊥) distribution, and compare them directly to the analytical marginal distributions they get from their solution.
Subrahmanyan: If this holds up as robustly as it does in shearing simulations for mildly relativistic temperatures, then this method might be a way to predict pressure anisotropies in things like pulsar wind nebulae or jets that we can’t quite reach with standard MHD codes.
Vera: They showed the solution works for several initial conditions, from nonrelativistic Maxwellians all the way up to ultrarelativistic Maxwell-Jüttner distributions, which is pretty comprehensive.
Jocelyn: It really shows how versatile this kinetic approach is; it’s not just a niche thing for one temperature but covers a wide range of particle speeds and initial states.
Subrahmanyan: The paper also suggests using these equations as a new closure relation for those MHD systems, which means you keep the pressure anisotropies in your stress-energy tensor even when you're using fluid approximations.
Vera: So the main implication is that we now have a more rigorous framework for how pressure anisotropy develops in collisionless magnetized plasmas across different relativistic scales.
Jocelyn: It’s a solid piece of work connecting kinetic theory, simulations, and general fluid models for anisotropic plasmas.
Subrahmanyan: We should keep an eye on how this method extends to other initial distributions like power-law or kappa distributions when we look at cosmic ray or solar wind studies.
More episodes
- 2501.00607-Topological Leptogenesis
- 2512.12258-New Massless Spectra from Cosmic String Cusps
- 2609.28637-Breaking de Sitter Symmetry with Cosmic Topology: Cosmological Role of Infrared Vacuum Freedom
- 2603.20409-The impact of gamma-ray propagation effects on indirect dark matter searches
- 2606.00649-Linear causality and stability constraints on relativistic second-order magnetohydrodynamics
- 2605.23802-Astrophysical Parameters of 5056 Open Star Clusters from Bayesian Nested Sampling with PARSEC Isochrones
- 2603.26882-Identifying heating processes in simulations with an entropy-based scheme: A single jet episode in a galaxy cluster
- 2605.31323-Time-reversed stochastic inflation in the quantum well
- 2512.10951-Evidence of galaxy cluster rotation and intrahalo gas motions in the cosmic microwave background
- 2501.12026-Synthetic Modelling of Polarized Dust Emission in Intermediate-Mass YSOs: II: Effects of Radiative Torque Disruption on Dust Grains in Protostellar Jets/Outflows