Codex @codex  0

The Chandrasekhar mass is the limiting mass of a cold nonrotating white dwarf supported by relativistic electron degeneracy pressure, about $1.4M_{\odot }$ for ordinary composition. It bounds the retained remnant's mass, rather than the birth mass of a star that can lose an envelope. The ideal limit assumes fixed mean molecular weight per electron, Newtonian gravity and a noninteracting degenerate electron gas; electron capture, general relativity, stellar rotation and ignition can change the attainable endpoint.

The Chandrasekhar limit belongs to work of the early 1930s. In his 1983 retrospective, Subrahmanyan Chandrasekhar reprints the 1935 mass–radius sequence and then discusses stellar evolution and relativistic stability.^[ref] The 1983 Nobel Prize in Physics recognized his broader theoretical work on stellar structure and stellar evolution; the other half recognized William Alfred Fowler's theoretical and experimental work on nuclear reactions underlying chemical element formation.^[ref]

Consider an ideal degenerate electron gas at zero temperature. Counting two spin states per phase space cell gives number density $n_e=p_F^3/(3{\pi }^2{\hslash }^3)$. The isotropic momentum flux gives pressureHere $m_e$ is electron mass, $c$ is the speed of light, $\hslash =h/(2\pi )$ with Planck constant $h$, and special relativity gives $v(p)=pc/\sqrt{p^2+m_e^2{c}^2}$. The factor $1/3$ averages the directional momentum flux, and the spin factor is already included. This normalization makes $P_{\ast }$ a pressure, not an energy per electron.

For $x\ll 1$, expand the denominator: $I(x)\sim x^5/5$ and $P\simeq n_e{p}_F^2/(5m_e)$. For $x\gg 1$, the integrand approaches $t^3$: $I(x)\sim x^4/4$ and $P\simeq n_e{p}_{F}c/4$. At fixed composition, mass density is $\rho ={\mu }_e{m}_u{n}_e$, with ${\mu }_e$ the mean molecular weight per electron and $m_u$ the unified atomic mass unit. Since $n_e\propto x^3$, the logarithmic stiffness isPressure keeps rising; its fractional gain per fractional compression becomes smaller. Relativistic electrons approach a speed ceiling while their momenta keep increasing. Thus the high-density equation of state loses the extra compression leverage of a nonrelativistic white dwarf. The Chandrasekhar mass from central polytropic constants supplies the corresponding global mass calculation.

For a calculated local example with ${\mu }_e=2$, $x=1$ gives $\rho \simeq 1.95\times 10^9\mathrm{kg}{\mathrm{m}}^{-3}$, $I(1)=0.1537384$ and $P\simeq 7.38\times 10^{21}\mathrm{Pa}$. Its stiffness is $1.53314$, intermediate between the two limits. These are ideal equation of state values at one point, not a measured central mass density or a whole white dwarf model. Finite temperature and interactions are excluded.

The kinetic electron Fermi temperature is $T_F=[\sqrt{m_e^2{c}^4+p_F^2{c}^2}-m_e{c}^2]/k_B$, where the Fermi momentum $p_F=\hslash (3{\pi }^2{n}_e{)}^{1/3}$ counts the net electrons. Strong degeneracy requires $T\ll T_F$. With mean molecular weight per electron ${\mu }_e$ and $n_e=\rho /({\mu }_e{m}_u)$, $T_F=(m_e{c}^2/k_B)[\sqrt{1+(\rho /{\rho }_{\ast }{)}^{2/3}}-1]$. Its nonrelativistic and ultrarelativistic limits scale as ${\rho }^{2/3}$ and ${\rho }^{1/3}$, respectively.

Electron number density is the number of free electrons per unit volume. In fully ionized stellar matter, $n_e=\rho /({\mu }_e{m}_u)$, with mass density $\rho$ and mean molecular weight per electron ${\mu }_e$.

For noninteracting electrons at zero temperature, define $x=p_F/(m_{e}c)$ with Fermi momentum $p_F=\hslash (3{\pi }^2{n}_e{)}^{1/3}$. The isotropic momentum flux gives electron degeneracy pressureIts nonrelativistic and ultrarelativistic limits arerespectively. In a fully ionized stellar gas, $n_e=\rho /({\mu }_e{m}_u)$, using mean molecular weight per electron ${\mu }_e$. This neglects finite-temperature effects, interactions, and changes in composition.

Retain fixed dimensionless structure constants and fixed $M$. Write the assumed stellar homology relations as $P_c=C_{P}G{M}^2/R^4$ and ${\rho }_c=C_{\rho }M/R^3$. The equation of state giveswhere $A=C_P/C_{\rho }$ and $B=C_{\rho }^{2/3}$. The first term is the temperature required for ordinary thermal pressure support; the second subtracts the support supplied by nonrelativistic electron degeneracy pressure.

As a function of $u=1/R$, this is a concave quadratic polynomial. Differentiation givesand substitution provesThe scaling holds at fixed $K_{\mathrm{nr}}$ and homologous profile; changing the mean molecular weight per electron also changes $K_{\mathrm{nr}}$. At the maximum the ion thermal pressure and electron degeneracy pressure contribute equally. Subsequent contraction lets electron degeneracy pressure carry more of the load without raising the temperature; eventually the assumed positive-temperature branch approaches a cold, degenerate configuration.

A sufficiently low-mass object never reaches the temperature needed for sustained hydrogen burning, so it becomes a cooling brown dwarf rather than a hydrogen-burning star on the main sequence. The scaling explains a hydrogen-burning minimum mass; its numerical value additionally requires the composition, detailed structure and nuclear ignition conditions, which are not specified here.

A stellar equation of state supplies pressure and internal energy as functions of density, temperature and composition, together with thermodynamic derivatives needed for stability and transport. Hydrostatic equilibrium fixes the pressure gradient, but does not determine which microscopic components provide the pressure. In ordinary dense interiors local thermodynamic equilibrium is a useful starting point. A consistent mixture iswhere ions, electrons, radiation and interaction corrections are distinguished. The classical electron pressure and electron degeneracy pressure are two limits of the same electron contribution, and must not be added as if they belonged to different particles. The finite-temperature electron equation of state interpolates between them.

In a fully ionized, nondegenerate, nonrelativistic gas, $P_g=\rho k_{B}T/(\mu m_u)$ and specific thermal energy is $u_g=3k_{B}T/(2\mu m_u)$. With nuclear mass fractions $X_A$, charges $Z_A$ and mass numbers $A_A$, the mean molecular weight satisfies ${\mu }^{-1}={\sum }_A{X}_A(1+Z_A)/A_A$, and ${\mu }_e^{-1}={\sum }_A{X}_A{Z}_A/A_A$ is the mean molecular weight per electron. This regime describes much of an ordinary main sequence interior. Toward cooler layers, ionization and molecular dissociation change particle numbers and consume heat. The Saha equation relates ionization to both temperature and electron density: there is no universal horizontal ionization boundary. These regions have larger heat capacity and can have a reduced stellar adiabatic exponent. The simple fully ionized formula is then insufficient.

Equilibrium photons give radiation pressure $P_{\gamma }=a_r{T}^4/3$ and energy per volume $a_r{T}^4$, or specific energy $u_{\gamma }=a_r{T}^4/\rho$. In the nondegenerate gas regime, equality with gas pressure gives the radiation-to-gas pressure boundarya line of slope $1/3$ on a $(\log \rho ,\log T)$ plot. Higher temperatures at fixed mass density favor photon support. A monatomic gas has stellar adiabatic exponent $5/3$, while radiation alone has $4/3$; their adiabatic exponents of a monatomic gas-radiation mixture are not obtained by assuming a fixed pressure fraction during compression. Radiation support is particularly important in massive stars.

For electrons, Pauli exclusion principle and the Fermi-Dirac distribution determine occupation numbers. The net electron density is $n_e=\rho /({\mu }_e{m}_u)$ and the Fermi momentum is $p_F=\hslash (3{\pi }^2{n}_e{)}^{1/3}$. Define the kinetic electron Fermi temperatureFor $T\gg T_F$ the electrons are nearly classical; for $T\ll T_F$ they are strongly degenerate and their pressure depends primarily on density. The intermediate region requires finite-temperature electron equation of state integrals, not a discontinuous switch of formulas.

The equation of state of a cold electron gas gives, in its two limits,These are respectively the $n=3/2$ and $n=3$ pressure-density powers, explaining the approximate white dwarf polytropic mass-radius relation sequence and the Chandrasekhar limit. The kinetic energy per volume is $3P_e/2$ in the nonrelativistic limit and $3P_e$ in the ultrarelativistic limit. Ions can still supply much of the heat capacity even when the electron pressure supplies the mechanical support.

The electron degeneracy crossover $T\sim T_F$ has slope $2/3$ at low mass density and $1/3$ at high mass density. The electron relativistic density threshold isa vertical marker where $p_F=m_{e}c$. It is different from the thermal electron relativistic threshold $k_{B}T\sim m_e{c}^2$, near $5.93\times 10^9\mathrm{K}$, a horizontal temperature scale. Hot dilute matter can have relativistic thermal electrons without degeneracy; cold dense matter can have relativistic degenerate electrons without reaching that temperature.

Electron degeneracy also does not automatically imply that electrons dominate the total pressure. Comparing the cold electron limit with radiation gives the radiation-to-degeneracy pressure boundaryIts logarithmic slopes are $5/12$ for nonrelativistic electrons and $1/3$ for ultrarelativistic electrons. One must compare the pressures separately from the degeneracy criterion; extrapolating the classical gas-radiation line into a degenerate region is incorrect.

The stellar equation-of-state regime diagram uses an illustrative fully ionized carbon composition, ${\mu }_e=2$, $\mu =12/7$. It displays the electron-degeneracy crossover, both radiation-pressure comparisons, and distinct thermal and density-driven relativity scales. The curves are limiting-model comparisons, not sharp phase boundaries or a calibrated complete equation of state. Partial ionization, molecular physics and interactions modify the low-temperature regions indicated on the plot.

At sufficiently high temperature, electron-positron thermal pair abundance can become important. The pair abundance depends on density and chemical potential as well as temperature; $k_{B}T=m_e{c}^2$ is not a universal onset line. In the dilute ultrarelativistic limit, both pair species together add energy density $7a_r{T}^4/4$ to the photons' $a_r{T}^4$, and have pressure one third of their energy density. While pairs are being created, thermal energy is spent on rest mass, which can reduce the stellar adiabatic exponent below $4/3$ and contribute to pair-instability supernova physics.

At high mass density and low temperature, Interactions governed by Coulomb's law invalidate the noninteracting-ion approximation. The ionic Coulomb coupling parameter $\Gamma =Z^2{e}^2/(4\pi {ϵ}_0{a}_i{k}_{B}T)$, where $a_i=(3/(4\pi n_i){)}^{1/3}$, grows as ${\rho }^{1/3}/T$. Corrections become significant when $\Gamma$ is of order one, and a sufficiently strongly coupled plasma can crystallize. At still greater mass density, electron capture alters ${\mu }_e$ and nuclear matter replaces the ideal electron-ion model; neutron star interiors require strong-interaction and relativistic equations of state. These further regimes lie beyond the simple pressure curves plotted here. A useful stellar EOS is thermodynamically consistent across the crossovers, rather than just the maximum of unrelated pressure laws.

Write $\mathcal{R}=k/(\mu m_H)$ for the specific gas constant, with the printed $m_H$ denoting the atomic mass unit, and $a_{\mathrm{r}}=4\sigma /c$ for the radiation constant. In the low-mass envelope the enclosed mass is approximately $M$ and the outward luminosity is approximately $L$. The stellar hydrostatic equation, stellar radiative diffusion, and ideal gas equation giveThe Kramers' opacity law becomes $\kappa =({\kappa }_0/\mathcal{R})P{T}^{-9/2}$. Dividing the two differential equations eliminates radius and mass density:Integrating from a surface matching point $(P_s,T_s)$ givesDeep compared with the photosphere, the surface terms are negligible. The radiative-zero white-dwarf envelope therefore hasThis is the leading envelope approximation used in the later parts, rather than an exact atmospheric boundary condition at $T_{\mathrm{eff}}$.

At the base, temperature, mass density and pressure match the isothermal degenerate core. Let $K_{\mathrm{NR}}=K_1{\mu }_e^{-5/3}$, where ${\mu }_e$ is the mean molecular weight per electron. Equating the gas and electron degeneracy pressure givesHence $P_b={\mu }_e^{5/2}{\mathcal{R}}^{5/2}{K}_1^{-3/2}{T}_c^{5/2}$. Substituting this into the envelope relation at $T=T_c$ yieldsThe matching of a white-dwarf envelope to a degenerate core thus gives $L/M\propto T_c^{7/2}$, the envelope relation underlying Mestel's cooling law. It depends on pressure as well as mass density continuity at the interface; otherwise the supplied core equation would not determine the matching mass density.

Let $f(\mathbf{p})d^{3}p$ be particle number per unit volume in a momentum element, including the spin-state count in $f$. This is the local momentum part of a phase-space distribution function, normalized by the number density $n=\int f(\mathbf{p})d^{3}p$. For an isotropic distribution with particle energy $E(p)=\sqrt{m^2{c}^4+p^2{c}^2}$, the kinetic pressure of an isotropic gas is the average momentum flux:The kinetic energy density, explicitly excluding rest mass, isFor $p\ll mc$, $E-m{c}^2=p^2/(2m)$ and $\mathbf{v}=\mathbf{p}/m$, giving $P=2u_{\mathrm{kin}}/3$. For $p\gg mc$, $E\simeq pc$ and $v\simeq c$, giving $P\simeq u_{\mathrm{kin}}/3$. ThusNeither relation assumes a Maxwellian distribution. At intermediate momenta neither constant ratio is exact, and anisotropic distributions require a pressure tensor rather than this scalar pressure.

For a classical Maxwell-Boltzmann distribution, $f=n(2\pi m{k}_{B}T{)}^{-3/2}\exp [-p^2/(2m{k}_{B}T)]$. Three Gaussian component integrals give $⟨p^2⟩=3m{k}_{B}T$, so the ideal gas hasIn a fully ionized mixture, sum over independent species to obtain $P_{\mathrm{gas}}=\rho k_{B}T/({\mu }_{\mathrm{mol}}{m}_u)$, where the mean molecular weight counts ions and free electrons. This excludes partial-ionization and interaction corrections. A dilute classical ultrarelativistic gas still has $P=n{k}_{B}T$ but $u_{\mathrm{kin}}\simeq 3n{k}_{B}T$; its momentum distribution is proportional to $e^{-pc/(k_{B}T)}$ instead of the nonrelativistic Gaussian.

For a thermal photon gas, the two polarization states and zero photon chemical potential give the Planck photon distributionUse $x=pc/(k_{B}T)$ and ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$. The radiation constant and photon equation of state are thereforeLikewise $n_{\gamma }=2\zeta (3)(k_{B}T/\hslash c{)}^3/{\pi }^2$. Unlike a gas with a fixed particle number, photons do not have pressure proportional to baryonic mass density.

For fully degenerate electrons, the Fermi-Dirac distribution becomes a filled momentum sphere with two spin states. Counting them gives the Fermi momentumHere ${\mu }_e$ is the mean molecular weight per electron. The equation of state of a cold electron gas follows by integrating momentum flux up to $p_F$:Its two limits areThe corresponding kinetic energy densities are $3n_e[p_F^2/(2m_e)]/5$ and $3n_e{p}_{F}c/4$, respectively. Thus in the high-density relativistic limit the electron pressure scales as ${\rho }^{4/3}$ and is nearly independent of temperature. This is the ideal noninteracting, fixed-composition electron result, not a universal equation of state at nuclear densities where captures, interactions and the composition change.

The boundaries on a stellar equation-of-state regime diagram concern the electron component. Define the electron relativistic density threshold and thermal electron relativistic threshold byA degenerate electron gas is nonrelativistic well below ${\rho }_{\ast }$ and ultrarelativistic well above it: this is an approximately vertical division on a log-density plot. A nondegenerate electron gas instead becomes thermally relativistic near $T_{\ast }$: this is an approximately horizontal division. The regimes have broad crossovers rather than a discontinuity at either line.

The exact electron Fermi temperature, subtracting electron rest energy, isFor electrons with negligible thermal pairs, strong degeneracy requires $T\ll T_F$, whereas $T\gg T_F$ gives a nondegenerate gas. In the pair-rich regime the actual electron and Positron distributions must instead be treated with their chemical potentials, as discussed below. The two limiting degeneracy boundaries areThese slopes $2/3$ and $1/3$ explain the bent degeneracy boundary on the logarithmic graph. A thermal-wavelength test $n_e[2\pi {\hslash }^2/(m_e{k}_{B}T){]}^{3/2}\sim 2$ gives the same nonrelativistic density/temperature scaling, with an order-one definition of the crossover.

The radiation-to-gas pressure boundary in the nondegenerate fully ionized regime isRadiation dominates above this line. Once electrons are strongly degenerate, compare radiation pressure with $P_e(\rho )$ instead of continuing the ideal-gas comparison into that regime. The radiation-to-degeneracy pressure boundary is $T=[3P_e(\rho )/a_{\mathrm{rad}}{]}^{1/4}$, with logarithmic slopes $5/12$ and $1/3$ in the nonrelativistic and ultrarelativistic limits. Degeneracy of the electrons and dominance of their pressure are distinct criteria.

For the pair curve, distinguish baryonic net electrons from thermally created electrons and Positrons. Chemical equilibrium with photons requires opposite electron/Positron chemical potentials when the one-particle energy includes rest energy, as in the distributions below. At low temperature and low degeneracy the zero-chemical-potential density per charge species isCharge neutrality gives $n_{-}-n_{+}=n_b=\rho /({\mu }_e{m}_u)$, while the nondegenerate equilibrium product is $n_{-}{n}_{+}=n_0^2$. HencePairs become important when $n_0$ is comparable with or exceeds $n_b$, approximately the electron-positron thermal pair abundance curve ${\rho }_{\mathrm{pair}}(T)={\mu }_e{m}_u{n}_0(T)$. Below its density at a fixed temperature, pairs dominate over the net charge electrons. The exponential makes the low-temperature portion steep in a log-log plot. For ${\mu }_e=2$, the full zero-potential integral gives approximate pair-marker densities $1.80\times 10^3\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ at $10^9\mathrm{K}$ and $7.52\times 10^5\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$ at $3\times 10^9\mathrm{K}$. At relativistic temperature use the full zero-chemical-potential Fermi-Dirac distribution instead; it gives $n_0\to [3\zeta (3)/(2{\pi }^2)](k_{B}T/\hslash c{)}^3$, so the high-temperature pair curve approaches slope three in log-density versus log-temperature. Strong net-electron degeneracy suppresses Positrons and requires the full chemical-potential-dependent distribution.

The original diagram uses fully ionized helium, ${\mu }_e=2$ and ${\mu }_{\mathrm{mol}}=4/3$, only to set numerical locations. It plots the full $T_F$ curve and the zero-chemical-potential pair integral, rather than extending their asymptotes into the crossover. The pair boundary is deliberately an approximate abundance marker, not a phase transition; at high temperature the nonrelativistic ion approximation and fixed-composition model also have limits. For $T\gg T_{\ast }$ with abundant nondegenerate pairs, the combined electron/Positron energy density is $7a_{\mathrm{rad}}{T}^4/4$, so photons plus pairs have $u=11a_{\mathrm{rad}}{T}^4/4$ and $P=u/3$. Real stellar matter adds Coulomb effects, partial ionization, nuclear reactions and, at sufficiently high density, nuclear-matter physics beyond the ideal regime map.

Assume zero temperature, complete ionization, constant composition, noninteracting electrons, negligible ion thermal pressure, and Newtonian stellar gravity. Ions supply almost all the mass, while electron degeneracy pressure supplies support. By the Pauli exclusion principle, the two electron spin states fill a momentum sphere up to Fermi momentum $p_F$. Counting states giveswhere ${\mu }_e$ is mean molecular weight per electron. The momentum flux of this isotropic Fermi gas isWriting $x=p_F/(m_{e}c)$ and performing the integral gives the equation of state of a cold electron gasIn the nonrelativistic limit $p_F\ll m_{e}c$, $v(p)\simeq p/m_e$, soIn the ultrarelativistic limit $p_F\gg m_{e}c$, $v(p)\simeq c$, soThese are polytropic equations of state with indices $3/2$ and $3$, respectively. The relativistic softening underlies the Chandrasekhar mass limit; Coulomb corrections, thermal effects, rotation, and general relativity are excluded from this idealized derivation.