Post-Newtonian Boltzmann Equation
- The post-Newtonian Boltzmann equation is a weak-field expansion of relativistic kinetic theory that incorporates 1/c² corrections in gravitational potentials and phase-space measures.
- It unifies collisionless and collisional dynamics through formulations like the Vlasov, BGK, and Marle models, ensuring consistency with 1PN conserved quantities.
- Applications include hydrodynamics, Jeans instability analysis, spherical stellar systems, and multicomponent self-gravitating media, highlighting critical reductions in Jeans mass.
The post-Newtonian Boltzmann equation is the weak-field, slow-motion expansion of relativistic kinetic theory for a one-particle distribution function in a gravitational field. In collisionless problems it is the post-Newtonian Vlasov equation for self-gravitating systems; in collisional problems it is coupled either to the full Boltzmann collision operator or to model operators such as BGK/Marle. Its characteristic feature is that the Newtonian streaming operator is supplemented by terms of order and higher involving scalar and vector gravitational potentials, corrected phase-space measures, and post-Newtonian source terms for the gravitational field equations. The formalism has been developed at $1$PN, $2$PN, and $2.5$PN orders and applied to hydrodynamics, dissipative transport, Jeans instability, spherical stellar systems, Brownian motion, and multicomponent self-gravitating media (Kremer, 2020, Kremer, 2021, Kremer et al., 2023, Kremer, 30 Apr 2026).
1. Weak-field kinetic framework
In the standard post-Newtonian hierarchy, Newtonian order gives the classical Boltzmann or Vlasov equation, $1$PN retains corrections of order , $2$PN retains corrections of order , and $2.5$PN introduces the first odd-order, radiative terms, with the $2.5$PN Boltzmann equation and equilibrium Maxwell–Jüttner distribution determined up to order $1$0 (Kremer, 2020, Kremer, 30 Apr 2026). The dynamical variable is the one-particle distribution function, written either as $1$1 in momentum variables or as $1$2 in three-velocity variables.
A common $1$3PN metric choice is Chandrasekhar’s form
$1$4
where $1$5 is the Newtonian gravitational potential, $1$6 is a scalar $1$7PN potential, and $1$8 is a vector $1$9PN potential (Kremer, 2021). Some treatments use the notation
$2$0
with $2$1 playing the analogous roles (Kremer, 2022). At $2$2PN and $2$3PN order, additional tensor and radiative potentials such as $2$4, and $2$5 enter the metric and the kinetic equation (Kremer, 2020, Kremer, 30 Apr 2026).
The post-Newtonian Boltzmann equation is therefore not a single fixed formula but a family of weak-field kinetic equations whose explicit appearance depends on approximation order, variable choice, and whether collisions are retained. A plausible implication is that the term designates a framework rather than one canonical operator.
2. Collisionless $2$6PN equation and corrected conserved quantities
For collisionless self-gravitating systems, the central object is the $2$7PN collisionless Boltzmann equation obtained by expanding the general relativistic Liouville operator. In Chandrasekhar-type notation one widely used explicit form is
$2$8
which is the equation used in post-Newtonian Jeans analyses (Kremer, 2021). Equivalent formulations may appear with a prefactor multiplying the Newtonian streaming term, but all retain only terms through order $2$9.
A major correction to earlier literature was established in “Kinetic Theory of Collisionless Self-Gravitating Gases: Post-Newtonian Polytropes” (Agón et al., 2011). That work argued that the $2.5$0PN collisionless equation previously used by Rezania and Sobouti is not consistent with the actual $2.5$1PN equations of motion, because the quantities identified there as energy and angular momentum are not true integrals of motion. The corrected $2.5$2PN collisionless Boltzmann equation is consistent with the $2.5$3PN test-particle dynamics, the correct conserved energy
$2.5$4
and, for spherical systems, the corrected angular momentum
$2.5$5
The same paper states that Jeans theorem remains valid at $2.5$6PN order: static equilibria are functions of the true integrals of motion (Agón et al., 2011).
This correction resolved a common misconception that the post-Newtonian collisionless Boltzmann equation can be constructed by a naive relativistic modification of the Newtonian operator. In the corrected formulation, consistency with the microscopic equations of motion is a necessary criterion.
3. Gravitational coupling, invariant measures, and equilibrium distributions
The post-Newtonian Boltzmann equation is coupled to post-Newtonian field equations through velocity moments of the distribution function. In kinetic form, the energy-momentum tensor is written as
$2.5$7
or, in equivalent species-dependent notation,
$2.5$8
with the invariant measure expanded consistently at $2.5$9PN order. One representative measure expansion is
$1$0
and the particle four-velocity has the $1$1PN form
$1$2
(Kremer, 2021, Kremer et al., 2023).
The corresponding Poisson-type field equations couple the kinetic moments to the gravitational potentials. In one common notation,
$1$3
while in multicomponent collisionless systems the sources are sums over species-specific moments of $1$4 (Kremer, 2021, Kremer et al., 2023).
Equilibrium is described by post-Newtonian Maxwell–Jüttner distributions. For a stationary background with vanishing hydrodynamic velocity, one $1$5PN form is
$1$6
with $1$7 the Maxwellian part (Kremer et al., 2023). In peculiar-velocity notation $1$8, the $1$9PN Maxwell–Jüttner distribution acquires corrections involving 0, and 1, and at 2PN the equilibrium distribution contains quartic, sextic, and octic peculiar-velocity terms together with 3, and the 4PN potentials 5 (Kremer, 2020). For non-perfect fluids, the Grad 14-moment distribution has also been expanded to 6PN order, yielding dissipative corrections linear in the pressure deviator, heat flux, and dynamic pressure (Kremer, 2021).
These kinetic ingredients support macroscopic balance laws. In the 7PN non-equilibrium theory, the conserved mass density
8
satisfies
9
and analogous $2$0PN balance equations exist for mass-energy density and momentum density, with explicit viscous-stress and heat-flux terms (Kremer, 2021).
4. Jeans instability and multicomponent self-gravitating media
A central application of the post-Newtonian Boltzmann equation is gravitational instability. In the one-component collisionless case, the perturbative construction uses a homogeneous background, the Jeans swindle, and plane-wave perturbations
$2$1
After solving the linearized kinetic equation and evaluating Gaussian velocity integrals, the resulting dispersion relation yields a post-Newtonian threshold
$2$2
so the critical wavelength is reduced and the Jeans mass is smaller than in the Newtonian theory. With the estimate $2$3, the paper gives
$2$4
which makes explicit that the $2$5PN correction shifts the instability threshold toward collapse at smaller masses (Kremer, 2021).
The multicomponent generalization was developed for a self-gravitating baryon–dark-matter mixture in “A self-gravitating system composed of baryonic and dark matter analysed from the post-Newtonian Boltzmann equations” (Kremer et al., 2023). There the kinetic model is built from two one-particle distribution functions,
$2$6
one for baryons and one for dark matter, each satisfying its own collisionless $2$7PN Boltzmann equation and coupled only through the common gravitational potentials. The perturbed field equations are sourced by moments of $2$8 and $2$9, and the dimensionless dispersion relation depends on the density ratio 0, the velocity-dispersion ratio 1, the post-Newtonian parameter 2, and the background potential 3. The paper shows that both dark matter and 4PN corrections reduce the Jeans mass relative to the standard one-component Newtonian criterion (Kremer et al., 2023).
Dark matter enters that framework in three distinct ways: it has its own distribution function obeying its own collisionless 5PN Boltzmann equation; it contributes separately to the source terms in the post-Newtonian field equations through its own energy-momentum tensor moments; and its equilibrium density 6 and velocity dispersion 7 define the reference Jeans scale (Kremer et al., 2023). The same paper applies the formalism to Bok globules, using 8, 9, and inferred values $2.5$0, and states that the model correctly reproduces the observational stability data (Kremer et al., 2023).
5. Spherical systems, Jeans equations, and post-Newtonian polytropes
For stationary, spherically symmetric stellar systems, the collisionless $2.5$1PN Boltzmann equation can be rewritten in spherical coordinates and reduced by taking velocity moments. In “Post-Newtonian Jeans Equation for Stationary and Spherically Symmetrical Self-Gravitating System” (Kremer, 2022), the radial moment yields the final $2.5$2PN Jeans equation
$2.5$3
with anisotropy parameter
$2.5$4
The same analysis derives $2.5$5 from the angular moment equation and couples the Jeans equation to the scalar Poisson equations
$2.5$6
(Kremer, 2022).
The closure used in that stationary problem is a $2.5$7PN Maxwell–Jüttner equilibrium distribution, from which the required higher velocity moments are evaluated (Kremer, 2022). Applied to a Hernquist model with a central black hole, the post-Newtonian scalar potential $2.5$8 is found to be positive and increasing toward small radius, and the paper states that it attenuates the rise of the velocity-dispersion profile relative to the purely Newtonian prediction (Kremer, 2022).
The corrected collisionless $2.5$9PN Boltzmann equation has also been used to construct post-Newtonian polytropes. For the ergodic distribution
$2.5$0
the resulting field equations are
$2.5$1
which give the Newtonian polytrope equation together with its post-Newtonian correction (Agón et al., 2011). The same work notes that sufficiently strong relativistic corrections can produce negative mass density in outer regions for some parameter choices, and interprets such cases as unphysical and indicative of the limits of the $2.5$2PN approximation (Agón et al., 2011).
6. Collisional, dissipative, and stochastic extensions
The post-Newtonian Boltzmann equation is not restricted to collisionless dynamics. In “Relaxation-Time Model for the Post-Newtonian Boltzmann Equation” (Kremer, 2023), the collision term is replaced by the BGK/Marle model
$2.5$3
and the Chapman–Enskog method is used to obtain the non-equilibrium distribution function, the hydrodynamic balance equations, and constitutive laws for viscous stress and heat flux. The paper identifies Newtonian transport coefficients
$2.5$4
and gives the $2.5$5PN-corrected coefficients
$2.5$6
so both shear viscosity and thermal conductivity depend explicitly on the Newtonian gravitational potential $2.5$7 (Kremer, 2023).
A more general non-equilibrium theory retaining the full Boltzmann collision operator was developed in “Post-Newtonian non-equilibrium kinetic theory” (Kremer, 2021). That paper combines the $2.5$8PN Boltzmann equation, the Eckart decomposition,
$2.5$9
and the Grad 14-moment distribution expanded to post-Newtonian order. It derives a $1$00PN Maxwell–Enskog transfer equation and the $1$01PN hydrodynamic equations for conserved mass density, mass-energy density, and momentum density. The dynamic pressure $1$02 is retained in the Grad distribution but drops out of the final $1$03PN hydrodynamics for a rarefied monatomic gas because its contribution is of order $1$04 (Kremer, 2021).
The collisional framework has also been reduced to a stochastic kinetic description. “Fokker-Planck equation for the Brownian motion in the post-Newtonian approximation” (Kremer, 15 Jul 2025) starts from a $1$05PN Boltzmann equation for a dilute binary mixture, assumes $1$06 and $1$07, and derives a post-Newtonian Fokker–Planck equation for the heavy particles by expanding the collision integral to second order in the small velocity increment. The derivation retains post-Newtonian corrections not only in the streaming operator but also in the invariant measure, invariant collision flux, equilibrium Maxwell–Jüttner distributions, and cross-section dependence on the flux. For hard spheres, the $1$08PN friction coefficient depends explicitly on $1$09, and the quoted $1$10PN extension depends on $1$11, $1$12, and $1$13 (Kremer, 15 Jul 2025). The same paper develops Newtonian and post-Newtonian linear stability analyses and states that for perturbation wavelengths smaller than the Jeans wavelength two propagating modes and one nonpropagating mode appear, while for wavelengths larger than the Jeans wavelength the perturbation either grows or decays, with the growing branch corresponding to instability (Kremer, 15 Jul 2025).
7. Higher-order post-Newtonian structure and regime of validity
Beyond $1$14PN, “Post-Newtonian Kinetic Theory” (Kremer, 2020) derives the Boltzmann equation in the second post-Newtonian approximation directly from the proper-time evolution of the one-particle distribution function along the particle worldline. The resulting $1$15PN equation contains higher-order contributions from $1$16, and $1$17, together with the $1$18PN Maxwell–Jüttner equilibrium distribution and the $1$19PN moment equations for mass density, mass-energy density, and momentum density (Kremer, 2020). An important structural point made there is that combining the mass and mass-energy equations yields the internal-energy equation in the first post-Newtonian approximation, so some conservation laws at a given PN order require higher-order kinetic information (Kremer, 2020).
At $1$20PN, “Boltzmann equation in the $1$21-post-Newtonian approximation” (Kremer, 30 Apr 2026) extends the kinetic theory through order $1$22 in harmonic coordinates. The collisionless streaming operator then contains, in addition to the Newtonian and conservative even-parity PN terms, odd-order radiative potentials such as $1$23, $1$24, $1$25, and $1$26, which are tied to time derivatives of the inertia tensor and represent radiation-reaction effects (Kremer, 30 Apr 2026). The same paper shows that these potentials also enter the equilibrium Maxwell–Jüttner distribution and the higher-order moments $1$27 and $1$28. It further notes that although the Boltzmann and metric expansions include $1$29PN terms, the final total-energy conservation law derived from the available moments does not yet display the full radiative energy-loss term; higher-order metric information would be required (Kremer, 30 Apr 2026).
Across the literature, several assumptions recur: weak gravitational fields, slow bulk and particle motions, truncation at a fixed post-Newtonian order, and either collisionless evolution or a near-equilibrium closure such as BGK or Grad. Many instability calculations also assume a stationary homogeneous background, the Jeans swindle, and plane-wave perturbations (Kremer, 2021, Kremer et al., 2023). A common misconception is that the post-Newtonian Boltzmann equation is only a $1$30PN collisionless construct. The published developments instead include collisionless and collisional theories, dissipative closures, stochastic reductions, and explicit extensions to $1$31PN and $1$32PN order (Kremer, 2020, Kremer, 2021, Kremer, 2023, Kremer, 15 Jul 2025, Kremer, 30 Apr 2026).