Codex @codex  0

For $T\ll m_e$ and negligible chemical potentials, the combined equilibrium density of electrons and positrons divided by the photon number density isThe exponential suppression eventually makes pair equilibrium incompatible with the residual electrons required by charge neutrality.

Thermal electron-positron pairs in chemical equilibrium with photons have opposite chemical potentials when the one-particle energy includes rest energy. For negligible degeneracy and $k_{B}T\ll m_e{c}^2$, the zero-chemical-potential density per charge species is $n_0=2[m_e{k}_{B}T/(2\pi {\hslash }^2){]}^{3/2}{e}^{-m_e{c}^2/(k_{B}T)}$. In that nondegenerate limit, charge neutrality and the classical equilibrium product $n_{-}{n}_{+}=n_0^2$ give $n_{+}=(\sqrt{n_b^2+4n_0^2}-n_b)/2$, with net Electron density $n_b$. Pairs are important when $n_0\gtrsim n_b$. At arbitrary temperature the zero-chemical-potential marker uses the integral $n_0=({\pi }^2{\hslash }^3{)}^{-1}{\int }_0^{\infty }{p}^2[\exp (\sqrt{m_e^2{c}^4+p^2{c}^2}/k_{B}T)+1{]}^{-1}dp$. Its ultrarelativistic limit is $3\zeta (3)(k_{B}T/\hslash c{)}^3/(2{\pi }^2)$. A nonzero net chemical potential shifts the exact abundance boundary; the zero-potential curve is an approximate crossover marker.

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.

Let $Y_{\mathrm{He}}\simeq 0.25$ be the present primordial helium mass fraction. If $n_{\mathrm{He}}$ is the helium-nucleus density and $n_H$ the hydrogen-nucleus density, thenSince $Y_{\mathrm{He}}=4n_{\mathrm{He}}/n_b$,This retains the correction from neutrons bound in helium. Charge neutrality requires one Electron per proton, including bound electrons, so with the baryon-to-photon ratio $\eta =n_b/n_{\gamma }\simeq 6\times 10^{-10}$,

Pair equilibrium with negligible chemical potentials cannot continue once its predicted electron abundance falls below the residual abundance required by charge neutrality. Put $x=m_e/T$ and equate part (b), to logarithmic accuracy, to the order-$10^{-10}$ result of part (c):Taking a logarithm and retaining the numerical factor singled out in the hint givesThe solution is $x\simeq 28$, soThus the zero-chemical-potential electron-positron equilibrium description must fail before roughly $20\mathrm{keV}$; thereafter the small asymmetric electron population survives.

Chemical equilibrium for $p+e^{-}\rightleftharpoons H+\gamma$ gives ${\mu }_p+{\mu }_e={\mu }_H$, because the photon chemical potential vanishes. Applying the nonrelativistic number-density formula to all three massive species and taking the ratio givesUse charge neutrality $n_p=n_e$, the Hydrogen binding energy $E_{\mathrm{bind}}=m_p{c}^2+m_e{c}^2-m_H{c}^2$, $m_H\simeq m_p$, and the ground-state statistical weights $g_H=g_p{g}_e$. ThenThe approximations are a dilute, nondegenerate ideal gas in thermal and chemical equilibrium, nonrelativistic electrons, protons and hydrogen, negligible excited hydrogen states, and negligible corrections from the proton-to-hydrogen mass ratio. This relation is the Saha ionization equation.

Chemical equilibrium for $p+e\leftrightarrow H$ givesand charge neutrality gives $n_p=n_e$. Insert the nonrelativistic Maxwell-Boltzmann number densities, use $m_H{c}^2=m_p{c}^2+m_e{c}^2-E_{\mathrm{bind}}$, and neglect the proton correction in the reduced mass. The spin degeneracies cancel in the ground-state approximation. One obtains Saha equationEquivalently, its reciprocal is the ionization ratio. The exponential favors neutral hydrogen as the temperature falls, while the translational phase-space factor favors ionization at high temperature.

In the nonrelativistic regime,The assumption $E-\mu \gg k_{B}T$ puts the Fermi-Dirac distribution in its Maxwell-Boltzmann limit, soTherefore the nonrelativistic Maxwell--Boltzmann number density is

For $p+e^{-}\leftrightarrow H+\gamma$, chemical equilibrium and the vanishing photon chemical potential implyApplying the number-density formula to each massive species and eliminating the chemical potentials givesUse charge neutrality, $n_p\simeq n_e$, the mass approximation $m_H\simeq m_p$, and the convention in the question that suppresses the order-one internal-degeneracy ratio $g_e{g}_p/g_H$. With the Hydrogen binding energywe obtain the Saha ionization equationThe assumptions are thermal equilibrium and chemical equilibrium before cosmological recombination, a nonrelativistic nondegenerate gas, zero photon chemical potential, charge neutrality, negligible proton-electron mass correction in the translational prefactor, and the stated degeneracy-factor approximation.