Table of Contents
Fetching ...

Kinetic Theory and the Mechanics of Isothermal Gas Spheres: Derivation of the Classical Emden--Chandrasekhar Equation via the Vlasov--Poisson Formalism

Steven D Miller

TL;DR

The paper derives the isothermal Emden--Chandrasekhar framework directly from the stationary Vlasov--Poisson equations by extremizing the Boltzmann entropy, yielding a Maxwell--Boltzmann distribution $\mathcal{F}_E \propto \exp[-\beta(\tfrac{|p|^2}{2m}+m\Phi)]$ and the Poisson--Boltzmann equation for self-gravitating Newtonian gases. Specializing to spherical symmetry recovers the Emden--Chandrasekhar equation, linking kinetic theory, thermodynamics, and stellar structure; the second variation of entropy reproduces the Antonov instability criteria with critical energy $E_c$, while the formalism naturally exhibits negative specific heat and shows uniform density is impossible for isothermal spheres. The work connects kinetic theory to classic stellar models, clarifying the self-consistency between microscopic phase-space dynamics and macroscopic hydrostatic equilibria, and extends the discussion to finite-size truncations (Bonnor--Ebert, King models) and relativistic generalizations via the Einstein--Vlasov system. Overall, the paper provides a unified kinetic-theoretic derivation of foundational isothermal-sphere results and situates them within the broader context of self-gravitating thermodynamics and stellar evolution.

Abstract

We present a derivation of the mechanics of isothermal gas spheres directly from the Vlasov--Poisson equation. By extremising the Boltzmann entropy, we obtain the Maxwell--Boltzmann distribution for a self-gravitating isothermal Newtonian gas, which is a stationary solution of the Vlasov--Poisson system. From this distribution, the corresponding Poisson--Boltzmann equation for the gravitational potential is deduced. The second variation of entropy reproduces the classical Antonov instability criterion: the critical energy is $E_c \simeq -0.335\,\frac{G M^2}{R}$, below which no local entropy maximum exists and the configuration becomes unstable (the so-called "gravothermal catastrophe"). In this work, we assume $E>E_c$, so all equilibria lie on the stable branch, and the Antonov instability does not affect the analysis. Specializing to spherical symmetry, we recover the classical equation for the isothermal gas sphere, as originally studied by Chandrasekhar, which has applications to the formation of red giant stars. We also derive the fundamental equation of hydrostatic equilibrium, the energy integral and the virial theorem directly from the stationary Vlasov--Poisson solution, demonstrating also that an isothermal gas exhibits negative specific heat. Furthermore, we show that an isothermal gas sphere of strictly constant density is an impossibility. This exposition emphasizes some of the deep connections and self-consistency between kinetic theory, statistical mechanics, and stellar structure, while highlighting some formal aspects of classical astrophysical models.

Kinetic Theory and the Mechanics of Isothermal Gas Spheres: Derivation of the Classical Emden--Chandrasekhar Equation via the Vlasov--Poisson Formalism

TL;DR

The paper derives the isothermal Emden--Chandrasekhar framework directly from the stationary Vlasov--Poisson equations by extremizing the Boltzmann entropy, yielding a Maxwell--Boltzmann distribution and the Poisson--Boltzmann equation for self-gravitating Newtonian gases. Specializing to spherical symmetry recovers the Emden--Chandrasekhar equation, linking kinetic theory, thermodynamics, and stellar structure; the second variation of entropy reproduces the Antonov instability criteria with critical energy , while the formalism naturally exhibits negative specific heat and shows uniform density is impossible for isothermal spheres. The work connects kinetic theory to classic stellar models, clarifying the self-consistency between microscopic phase-space dynamics and macroscopic hydrostatic equilibria, and extends the discussion to finite-size truncations (Bonnor--Ebert, King models) and relativistic generalizations via the Einstein--Vlasov system. Overall, the paper provides a unified kinetic-theoretic derivation of foundational isothermal-sphere results and situates them within the broader context of self-gravitating thermodynamics and stellar evolution.

Abstract

We present a derivation of the mechanics of isothermal gas spheres directly from the Vlasov--Poisson equation. By extremising the Boltzmann entropy, we obtain the Maxwell--Boltzmann distribution for a self-gravitating isothermal Newtonian gas, which is a stationary solution of the Vlasov--Poisson system. From this distribution, the corresponding Poisson--Boltzmann equation for the gravitational potential is deduced. The second variation of entropy reproduces the classical Antonov instability criterion: the critical energy is , below which no local entropy maximum exists and the configuration becomes unstable (the so-called "gravothermal catastrophe"). In this work, we assume , so all equilibria lie on the stable branch, and the Antonov instability does not affect the analysis. Specializing to spherical symmetry, we recover the classical equation for the isothermal gas sphere, as originally studied by Chandrasekhar, which has applications to the formation of red giant stars. We also derive the fundamental equation of hydrostatic equilibrium, the energy integral and the virial theorem directly from the stationary Vlasov--Poisson solution, demonstrating also that an isothermal gas exhibits negative specific heat. Furthermore, we show that an isothermal gas sphere of strictly constant density is an impossibility. This exposition emphasizes some of the deep connections and self-consistency between kinetic theory, statistical mechanics, and stellar structure, while highlighting some formal aspects of classical astrophysical models.

Paper Structure

This paper contains 22 sections, 12 theorems, 262 equations.

Key Result

Theorem 3.4

Any steady-state (time-independent) solution of the collisionless Boltzmann equation can be expressed as a function of the integrals of motion of the system. That is, if $I_1, I_2, \dots$ are integrals of motion, then

Theorems & Definitions (29)

  • Definition 3.1: Phase space and distribution function
  • Definition 3.2
  • Remark 3.3
  • Theorem 3.4: Jeans Theorem
  • Lemma 3.5: Stationary Maxwell--Boltzmann Distribution
  • proof
  • Theorem 3.6: Maxwell--Boltzmann solution of the stationary Vlasov--Poisson equation
  • proof
  • Corollary 3.7
  • proof
  • ...and 19 more