Codex @codex  0

For the reaction $p+e\leftrightarrow H+\gamma$, chemical equilibrium requires ${\mu }_p+{\mu }_e={\mu }_H$ because thermal photons have zero chemical potential. Insert the nonrelativistic Maxwell-Boltzmann distribution number densities and use $m_H=m_p+m_e-\mathcal{B}$:For ground-state hydrogen, including its spin states, $g_e=g_p=2$, $g_H=4$; hence the degeneracy factor is one. Since $m_p/m_H\simeq 1$, the right side is $(m_{e}T/2\pi {)}^{3/2}{e}^{-\mathcal{B}/T}$.

Charge neutrality gives $n_p=n_e$, while baryon conservation gives $n_b=n_p+n_H$. Consequently $n_e=n_p=X_e{n}_b$ and $n_H=(1-X_e)n_b$. Inverting the preceding ratio and inserting the photon number density with $n_b=\eta n_{\gamma }$ gives the hydrogen-only Saha equationWriting $y=\mathcal{B}/T$, its coefficient is $(2\zeta (3)/{\pi }^2)\eta (2\pi \mathcal{B}/m_e{)}^{3/2}\simeq 3.2\times 10^{-16}$, consistent with the rounded $3\times 10^{-16}$. The negligible mass-ratio correction and ground-state approximation are the assumptions behind this form.

Take $A>0$ for the repulsive molecular interaction, $p_0$ as the common reference pressure, and let $T_i=T_m-\delta T$ be the ice-film temperature. All ratios involving $T_m$ use its absolute value, approximately $273.15\mathrm{K}$, rather than the numerical Celsius value zero. The premelting film permits water to reach and freeze at the bottom of the ice while the upper surface is displaced upward without freezing there.

For equal phase mass densities, equality of the solid and liquid chemical potentials gives the Clapeyron pressure relation for equal-density phases. Indeed, $d\mu =-sdT+dp/\rho$ and $s_l-s_s=L/T_m$, so expansion around coexistence at $(T_m,p_0)$ yieldsThis is a pressure difference between phases, not a same-pressure Clausius-Clapeyron slope. The molecular disjoining pressure supports this difference:Treat the normal solid load as its effective pressure in this planar model. The gravitational load corresponding to hydrostatic pressure is $p_s=p_0+\rho gh$, so $p_l=p_0+\rho gh-p_T$. Refer the bath pressure to the sheet's upper gravitational datum, or neglect the sheet-scale hydrostatic head. The Darcy flow law then gives the upward supply and the ice growth rateHere $\Pi$ is permeability of a porous medium, $\mu$ is dynamic viscosity, and equal densities identify supplied liquid volume with added ice volume to leading order in $a/h$.

To recover the printed law, use the usual hydraulic control of frost heave idealization: the sheet's top is at $T_m$, the liquid film and ice have a common thermal conductivity $K$, and the flow is slow enough that latent-heat production is small compared with the conductive heat passing through the layer. In addition to a quasi-steady ice temperature field, this needsThe leading heat flux is consequently continuous through the film and ice. Their thermal resistances giveLet $B_T=\rho L\Delta T/T_m$, a pressure scale. The thermomolecular pressure and the molecular law then giveSubstitute this pressure into the Darcy flow law:The first term draws water toward the undercooled ice; the increasing gravitational load eventually cancels it.

Quasi-steady temperature alone is not sufficient to fix this exact prefactor. With film conductivity $K_l$, ice conductivity $K_s$, and nonnegligible latent heat, the appropriate additional balance isTogether with the Darcy equation, this is the implicit growth model. Even in the slow-flow limit, a conductivity ratio $K_s/K_l=r$ changes the leading driving pressure by $r^{3/4}$. Taking $r=2$ is a counterexample to obtaining the printed prefactor from quasi-steadiness alone. Thus the boxed equation is the intended equal-conductivity, thin-film, hydraulically limited model, rather than a consequence of only the stated quasi-steady assumption. This distinction is the heat-balance correction to a premelted-film growth model.

For that reduced equation defineThe scale $h_{\ast }$ balances $F{h}_{\ast }^{-3/4}$ against $\rho g{h}_{\ast }$; $F$ has units of pressure times length to the power $3/4$, so $h_{\ast }$ is a length and $t_{\ast }$ is a time. They yield

Set $u={\eta }^{7/4}$. The equation becomes linear:For initial thickness $\eta (0)={\eta }_0\ge 0$, the explicit solution of gravity-limited frost heave isThe usual zero-initial-thickness sketch usesIts early and late behaviour areFor $0<\eta <1$, the derivative is positive and ${\eta }^{\prime \prime }=(-\frac{3}{4}{\eta }^{-7/4}-1){\eta }^{\prime }<0$, so the curve is increasing and concave downward. It approaches the stable thickness $h_{\ast }$ exponentially. Initial thickness above $h_{\ast }$ instead relaxes downward; at $h_{\ast }$ it is stationary.

The singular slope at zero belongs to the formal reduced solution. Since $a/h\propto h^{-3/4}$, the thin-film assumption eventually fails as $h\to 0$. The early power law therefore describes an intermediate continuum regime after any microscopic startup, not arbitrarily early physical times. A finite positive ${\eta }_0$ supplies a regular initial condition when that regime begins. Gravity limits the ultimate thickness; molecular suction and the Darcy resistance set the growth toward it.

An order parameter distinguishes thermodynamic phases and transforms under the symmetry that may be broken. For a scalar ferromagnet it is the magnetization per site; a disordered zero-field phase has $m=0$, while ordered pure thermodynamic phases have $m=\pm m_0$. In a fluid one can instead use the density measured relative to its critical value; the field conjugate to it is then related to the chemical potential rather than literally a magnetic field. A nonzero value in the presence of an explicit conjugate field is not by itself evidence of a spontaneous transition.

The LG theory describes a slowly varying local order parameter by a symmetry-constrained free-energy functional, for exampleAt $h=0$, a ${\mathbb{Z}}_2$ symmetry excludes odd powers. The coefficients are assumed analytic functions of the controls near the transition, and the expansion is stabilized by a positive highest retained even coefficient. The Landau approximation obtains equilibrium by minimizing this functional, neglecting long-wavelength fluctuation corrections. Its nonconvex local potential describes distinct candidate phases and mean-field metastability; the exact thermodynamic potential is convexified when macroscopic mixtures are admitted.

In a finite symmetric system, the zero-field expectation of $m$ can vanish even below the transition because both ordered orientations are sampled. Spontaneous order is defined by taking the thermodynamic limit before removing a selecting field, for example $m_0={\lim }_{h\downarrow 0}{\lim }_{N\to \infty }⟨m⟩$. This distinction also matters for the connected correlations used below.

Two effects must be distinguished in Big Bang nucleosynthesis. The positive extra neutrino energy density increases the Hubble parameter through the Friedmann equation. Holding weak rates fixed, faster expansion causes earlier neutron-proton cosmological weak freeze-out, leaving more neutrons, and leaves less time for neutron decay before nuclear reactions begin. Both tendencies increase the primordial helium mass fraction.

Electron-flavor degeneracy also changes the charged-current balance directly. Chemical equilibrium for $n+{\nu }_e\leftrightarrow p+e^{-}$ gives, with negligible electron chemical potential,This is the primordial helium response to electron-neutrino degeneracy. Positive ${\xi }_{{\nu }_e}$ favors conversion of neutrons to protons and reduces the neutron abundance; negative ${\xi }_{{\nu }_e}$ favors the opposite balance. When neutrons are the limiting ingredient and nearly all surviving neutrons enter helium-4, the neutron-limited helium synthesis estimate isIn this usual proton-rich regime, the direct electron-neutrino effect lowers $Y_p$ for positive degeneracy and raises it for negative degeneracy, other conditions held fixed. At sufficiently large negative degeneracy the equilibrium ratio can instead exceed one; protons then become limiting. The corresponding ideal complete-capture estimate is $Y_p\simeq 2\min (r,1)/(1+r)$, so the neutron-limited trend is not a universal monotonic law at arbitrary negative chemical potential. Degeneracy in the other flavors primarily affects helium through the expansion rate in this simplified discussion.

For very large degeneracy the weak reaction rates themselves also change, so a quantitative freeze-out calculation must include both the altered distributions and the altered expansion rate. The two effects can compete for a positive electron-neutrino asymmetry. The primordial helium mass fraction change is not determined by $|{\mu }_{\nu }|$ alone: its flavor and sign matter as well as its energy density.

Use $\hslash =c=k_B=1$ and count one populated helicity per neutrino or antineutrino. For massless particles, the Fermi-Dirac distribution givesWhen ${\mu }_{\nu }/T_{\nu }\gg 1$, the occupied states form an almost sharp Fermi sea. Replacing the distribution by $\Theta ({\mu }_{\nu }-p)$ givesThe relative finite-temperature correction is of order $(T_{\nu }/{\mu }_{\nu }{)}^2$.

Chemical equilibrium gives ${\mu }_{\bar{\nu }}=-{\mu }_{\nu }$. For positive large ${\mu }_{\nu }$, the corresponding antineutrino distribution is dilute, soIt is exponentially negligible. The neutrino degeneracy parameter ${\xi }_{\nu }={\mu }_{\nu }/T_{\nu }$ is conserved in the assumed adiabatic massless evolution: redshifting preserves the distribution with both $T_{\nu }$ and ${\mu }_{\nu }$ proportional to $a^{-1}$.

There is a sign qualification in the printed absolute-value formula. If ${\mu }_{\nu }<0$, the distribution of neutrinos themselves is exponentially suppressed, and the antineutrinos, with positive chemical potential, form the degenerate sea. In either case the degenerate neutrino and antineutrino energy density isThus the formula with $|{\mu }_{\nu }|$ describes the dominant member of the pair, or their leading total, rather than the named neutrino distribution for both signs. No extra factor of two is present: only one member is densely occupied. Additional populated internal states multiply the result by their degeneracy.

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.