Lane-Emden Equation

Lane-Emden Equation is a dimensionless second-order nonlinear ordinary differential equation in stellar structure and theoretical astronomy that governs the interior radial density, pressure, and gravitational potential profiles of a spherically symmetric, self-gravitating fluid sphere obeying a polytropic equation of state, P = Kρ^{1 + 1/n}, in Newtonian hydrostatic equilibrium. First formulated in 1870 by the American physicist Jonathan Homer Lane in an investigation of the convective internal temperature and density of the Sun and systematically generalized in 1907 by the Swiss astrophysicist Robert Emden in his monograph Gaskugeln ("Gas Spheres"), the equation reduces the coupled partial differential equations of stellar hydrostatics and Poisson gravity to a single universal curve θ(ξ) parameterized by the polytropic index n. In the 1920s and 1930s, Sir Arthur Stanley Eddington used the n = 3 polytrope to build his "Standard Model" of radiative main-sequence stars (and the Schönberg–Chandrasekhar limit), while Subrahmanyan Chandrasekhar used the transition from n = 1.5 (non-relativistic electron degeneracy) to n = 3 (ultra-relativistic electron degeneracy) to derive the maximum stable mass of a white dwarf star at the Chandrasekhar limit.

Polytropic Equation of State and Newtonian Hydrostatics

Consider a non-rotating, unmagnetized, spherically symmetric fluid configuration in which the isotropic pressure P(r) at radial distance r from the center depends only on the local mass density ρ(r) through the polytropic relation:

P(r) = K ρ(r)^γ = K ρ(r)^{1 + 1/n}

where K is the polytropic constant (determined by specific entropy in an adiabatic convective gas or by universal quantum constants ℏ, m_e, and c in a degenerate Fermi gas), γ ≡ 1 + 1/n is the polytropic exponent, and n = 1/(γ − 1) is the polytropic index.

A spherical shell at radius r in mechanical equilibrium satisfies two coupled differential equations—the equation of hydrostatic equilibrium and mass continuity (or equivalently Poisson's equation for the Newtonian gravitational potential Φ):

1. Hydrostatic Equilibrium: (1 / ρ) dP / dr = −G M(r) / r² = −dΦ / dr
2. Mass Continuity: dM(r) / dr = 4π r² ρ(r)
3. Spherical Poisson Equation: (1 / r²) d/dr [ (r² / ρ) dP/dr ] = −4π G ρ(r)

Dimensionless Lane-Emden Transformations (θ and ξ)

To eliminate the physical constants (G, K, and the central density ρ_c ≡ ρ(0)) from the differential equation, introduce the dimensionless Lane-Emden variable θ(ξ) (sometimes called the polytropic temperature), which normalizes the radial density profile to unity at the stellar center:

ρ(r) = ρ_c [θ(ξ)]ⁿ ⟹ P(r) = P_c [θ(ξ)]ⁿ⁺¹

where P_c = K ρ_c^{1 + 1/n} is the central pressure. Notice that for an ideal gas obeying P = (ρ / μm_u)k_BT with uniform mean molecular weight μ, the temperature scales linearly with θ(ξ): T(r) = T_c θ(ξ). Substituting ρ(r) = ρ_c θⁿ into the left-hand side of the spherical Poisson equation yields:

(1 / ρ) dP / dr = K (n + 1) ρ_c^1/n (dθ / dr)

Next, define the dimensionless radial coordinate ξ scaled by the characteristic length parameter α:

r = α ξ, where α ≡ [ ((n + 1) K ρ_c^{(1/n) − 1}) / (4π G) ]^1/2

Substituting r = αξ and ρ = ρ_c θⁿ into the spherical Poisson equation cancels all dimensional factors and produces the canonical Lane-Emden equation:

(1 / ξ²) d/dξ ( ξ² dθ / dξ ) = −θⁿ or equivalently d²θ / dξ² + (2 / ξ) (dθ / dξ) + θⁿ = 0

Central Boundary Conditions and the Stellar Surface (ξ₁)

Integrating outward from the center of the star (ξ = 0) requires two boundary conditions:

1. Normalized Central Density: By definition, ρ(0) = ρ_c, so θ(0) = 1.
2. Vanishing Central Gravity: Because the enclosed mass M(r) → (4π/3)ρ_cr³ vanishes as r → 0, the gravitational acceleration g(r) = GM(r)/r² ∝ r must vanish at the origin, requiring a flat pressure and density peak at the center: dθ / dξ|_ξ=0 = 0.

Expanding θ(ξ) in a Taylor series around ξ = 0 with these boundary conditions gives the universal near-center behavior for any index n:

θ(ξ) = 1 − (1 / 6)ξ² + (n / 120)ξ⁴ − [ n(8n − 5) / 15,120 ]ξ⁶ + …

As ξ increases outward from zero, θ(ξ) decreases monotonically. For any polytropic index n < 5, θ(ξ) crosses zero at a finite first root ξ = ξ₁ (where θ(ξ₁) = 0), which defines the physical surface of the star where both density ρ(R) and pressure P(R) vanish at physical radius R = αξ₁.

Although the Lane-Emden equation must be integrated numerically (e.g., via fourth-order Runge–Kutta methods) for general values of n, closed-form analytical solutions exist for three specific indices:

1. Incompressible Homogeneous Sphere (n = 0, γ → ∞)

When n = 0, θ⁰ = 1, so the density ρ(r) = ρ_c is constant throughout the interior (representing an incompressible liquid sphere or first-order rocky planet). The differential equation (1/ξ²)d/dξ(ξ²dθ/dξ) = −1 integrates twice directly to give:

θ₀(ξ) = 1 − (ξ² / 6), with surface zero at ξ₁ = √6 ≈ 2.44949

2. Linear Sinc Polytrope (n = 1, γ = 2)

When n = 1, the Lane-Emden equation becomes a linear Helmholtz/spherical Bessel differential equation, d²θ/dξ² + (2/ξ)dθ/dξ + θ = 0. Substituting u(ξ) = ξ·θ(ξ) reduces the equation to the simple harmonic oscillator d²u/dξ² + u = 0. Applying θ(0) = 1 and θ'(0) = 0 yields the spherical sinc function:

θ₁(ξ) = sin(ξ) / ξ, with surface zero at ξ₁ = π ≈ 3.14159

Remarkably, for n = 1 the length scale α = [K / (2πG)]^1/2 is independent of the central density ρ_c. Consequently, all n = 1 polytropes with a given equation-of-state constant K possess the same physical radius R = π[K/(2πG)]^1/2 regardless of their total mass M. This property explains why gas giant planets and brown dwarfs spanning 1 to 80 Jupiter masses—whose interior equation of state is governed by a combination of metallic hydrogen Coulomb interactions (n ≈ 1) and partially degenerate electrons (n = 1.5)—all maintain nearly the same radius (R ≈ 1 R_J ≈ 70,000 km).

3. Plummer Sphere (n = 5, γ = 6/5)

Discovered by Arthur Schuster (1883) and applied by Henry C. Plummer (1911) to model the spatial distribution of stars in globular clusters, the n = 5 case admits the algebraic solution:

θ₅(ξ) = (1 + ξ² / 3)^−1/2

Although θ₅(ξ) > 0 for all finite ξ—meaning the physical radius is infinite (ξ₁ = ∞)—its total integrated mass converges to a finite value because −ξ²dθ₅/dξ = (ξ³/3)(1 + ξ²/3)^−3/2 → √3 as ξ → ∞. For all n > 5, both the radius and the total enclosed mass diverge to infinity.

Integrating the mass continuity equation from ξ = 0 to the stellar surface ξ = ξ₁ using the Lane-Emden equation itself to substitute ξ²θⁿ = −d/dξ(ξ²dθ/dξ) yields the total stellar mass M and radius R:

• Stellar Radius: R = α ξ₁ = [ ((n + 1) K) / (4π G) ]^1/2 ρ_c^{(1 − n) / (2n)} ξ₁
• Total Enclosed Mass: M = 4π α³ ρ_c ( −ξ₁² dθ/dξ|_ξ1 ) = 4π [ ((n + 1) K) / (4π G) ]^3/2 ρ_c^{(3 − n) / (2n)} ( −ξ₁² dθ/dξ|_ξ1 )
• Central-to-Mean Density Ratio: ρ_c / ρ̄ = −ξ₁ / (3 dθ/dξ|_ξ1)

Eliminating the central density ρ_c between the radius R and mass M expressions when K is fixed yields the fundamental polytropic mass–radius relation:

R ∝ K^{n / (3 − n)} M^{(1 − n) / (3 − n)}

This single algebraic exponent governs three key regimes of stellar astrophysics:

1. Non-Relativistic Degenerate Stars (n = 1.5, γ = 5/3): In low-mass white dwarfs (M ≲ 0.6 M_☉) and neutron stars supported by non-relativistic Fermi–Dirac degeneracy pressure (P = K₁ρ^5/3), substituting n = 3/2 gives (1 − n)/(3 − n) = (−1/2)/(3/2) = −1/3, so R ∝ M^−1/3. More massive white dwarfs are physically smaller in volume (M R³ = constant).
2. Eddington Standard Model & Ultra-Relativistic Degeneracy (n = 3, γ = 4/3): When n = 3, the exponent (3 − n) / (2n) on ρ_c in the mass equation vanishes (3 − n = 0)! Consequently, the central density ρ_c cancels out of the mass formula altogether:

M(n = 3) = 4π ( K / πG )^3/2 ( −ξ₁² dθ/dξ|_ξ1 )_n=3, where (−ξ₁² dθ/dξ|_ξ1)_n=3 = 2.018236 In Eddington's radiative standard model (where the ratio β = P_gas/P_tot of gas pressure to total gas-plus-radiation pressure P_rad = aT⁴/3 is constant throughout the star), this fixes K as a function of β, allowing Eddington to derive the mass–luminosity relation of main-sequence stars. In an ultra-relativistic white dwarf where K = K₂ = (ℏc/4)(3π²)^1/3(μ_em_H)^−4/3 is fixed by fundamental constants, the n = 3 mass equation yields a unique maximum mass independent of radius—the Chandrasekhar limit (M_Ch ≈ 1.44 M_☉).

Why does a fully convective star obey an n = 1.5 polytrope?

In a star or planetary interior where energy is transported outward by vigorous subsonic adiabatic convection (such as M-dwarf red dwarfs with M < 0.35 M_☉, pre-main-sequence Hayashi-track protostars, or the outer convective envelope of the Sun), ascending and descending fluid parcels expand and compress along an adiabat with virtually constant specific entropy s. For an ideal monatomic ionized gas with three translational degrees of freedom, the ratio of specific heats is γ_ad = c_P/c_V = 5/3. Equating γ_ad = 5/3 to the polytropic exponent 1 + 1/n gives 1/n = 2/3, or n = 1.5.

What is the physical meaning of the Emden–Chandrasekhar isothermal sphere (n → ∞)?

When n → ∞ (γ = 1), the equation of state becomes isothermal, P = c_s² ρ, with constant sound speed c_s = (k_BT / μm_u)^1/2. Because the standard substitution ρ = ρ_cθⁿ becomes singular as n → ∞, one instead sets ρ(r) = ρ_c exp(−ψ(ξ)), yielding the isothermal Lane-Emden equation: (1/ξ²)d/dξ(ξ²dψ/dξ) = exp(−ψ). At large radii (ξ ≫ 1), its density decays as an inverse-square law ρ(r) ∝ r⁻², producing a flat rotation curve and forms the theoretical basis for Bonnor–Ebert spheres in star-forming molecular cloud cores and simplified dark-matter halo models.

Can the Lane-Emden equation be used for neutron stars?

For low-mass neutron stars where Newtonian gravity and non-relativistic neutron degeneracy (n = 1.5) provide a rough first approximation, the Lane-Emden equation correctly predicts the inverted R ∝ M^−1/3 scaling. However, because a neutron star's compactness parameter GM/(Rc²) ≈ 0.15–0.30 is large, accurate neutron-star structure requires replacing Newtonian hydrostatic equilibrium with the general-relativistic Tolman–Oppenheimer–Volkoff (TOV) equation and a nuclear-matter equation of state.

1. Lane, J. H. (1870). On the theoretical temperature of the Sun, under the hypothesis of a gaseous mass maintaining its volume by its internal heat, and depending on the laws of gases as known to terrestrial experiment. American Journal of Science, s2-50(148), 57–74. https://doi.org/10.2475/ajs.s2-50.148.57
2. Plummer, H. C. (1911). On the problem of distribution in globular star clusters. Monthly Notices of the Royal Astronomical Society, 71(5), 460–470. https://doi.org/10.1093/mnras/71.5.460
3. Chandrasekhar, S. (1931). The maximum mass of ideal white dwarfs. The Astrophysical Journal, 74, 81–82. https://doi.org/10.1086/143324
4. Horedt, G. P. (1986). Seven-digit tables of Lane-Emden functions. Astrophysics and Space Science, 126(2), 357–408. https://doi.org/10.1007/BF00639386