Codex @codex  0

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.

The useful high-density core bound from TOV compactness follows without choosing the unknown high-density equation of state. Let $a$ be the radius at which the decreasing energy density first reaches ${\rho }_0$, and put $\mu =m(a)$, $p_0=p({\rho }_0)$. The Tolman–Oppenheimer–Volkoff equation gives $p^{\prime }<0$ away from a nonempty centre; since $dp/d\rho >0$, energy density decreases outward. HenceFor the local stellar compactness, if $z\ge 0$, $1-z+\sqrt{1+z}\le 2$, since its derivative is $-1+(2\sqrt{1+z}{)}^{-1}<0$. Applying the supplied compactness inequality at the core boundary therefore gives $\mu /a<4/9$. Combining the two bounds yieldsThe pressure-dependent form further restricts the allowed core-boundary data:In particular, positivity of its right side gives $6\pi a^2{p}_0<3$, so $a^2<1/(2\pi p_0)$ when $p_0>0$. These restrictions depend only on known quantities, not on the equation of state above ${\rho }_0$.

The intended stellar-mass argument now attaches the known low-density envelope to each admissible pair $(a,\mu )$ with pressure $p_0$. Its Tolman–Oppenheimer–Volkoff equation involves only the known equation of state. For a physical low-density envelope with a uniformly bounded mass on this bounded set of core data, taking the supremum of the envelope-plus-core masses gives an upper mass limit independent of any high-density continuation. Stars with ${\rho }_c\le {\rho }_0$ are also determined entirely by the known equation of state. Unknown dense-core physics cannot evade the displayed core bounds.

A qualification is necessary for a claim about the entire star: the printed monotonicity assumptions alone do not guarantee a finite maximum total mass for an arbitrary known low-density equation of state. An explicit low-density polytropic mass divergence occurs for $p=K{\rho }^{5/4}$ with fixed $K>0$. In the low-density Newtonian limit, this is a stellar polytrope of index $n=4$, whose Lane-Emden equation has a finite first zero. With $\rho ={\rho }_c{\vartheta }^4$ and $r=\ell \xi$, hydrostatic balance fixeswhere ${\mu }_4>0$ is the dimensionless Lane-Emden surface mass constant. Thus $M\to \infty$ as ${\rho }_c\to 0$, although $M/R\propto {\rho }_c^{1/4}\to 0$. The first Lane-Emden zero is simple. The continuous dependence of an ODE solution on parameters therefore preserves this finite surface for sufficiently small relativistic parameter. The relativistic correction is controlled by $p_c/{\rho }_c=K{\rho }_c^{1/4}\to 0$, so the corresponding regular finite-radius Tolman–Oppenheimer–Volkoff equation solutions have the same limiting scalings. They never sample $\rho >{\rho }_0$, and satisfy the local positivity and monotonicity assumptions in their nonvacuum interiors. This shows exactly why a suitable bounded-envelope hypothesis, or an appropriate restriction to physically stable stellar models, must accompany the intended total-mass conclusion. The compactness inequality supplies a bound on the unknown core, not a bound on the mass of every conceivable dilute envelope.

A regular central expansion of the TOV equations requires $m(0)=0$. Once ${\rho }_c$ is chosen, the barotropic stellar equation of state fixes $p_c=p({\rho }_c)$. Smoothness at the centre forces the radial first derivatives of $p$ and $\rho$ to vanish. Expanding the Tolman–Oppenheimer–Volkoff equation givesSince $dp/d\rho >0$, the equation of state is locally invertible, and these coefficients determine the corresponding energy density expansion as well. Integrating the regular solution outward gives the surface at the first zero of the pressure, when such a finite surface exists. The ordinary differential equation uniqueness theorem then fixes the interior and its exterior mass from ${\rho }_c$.

The additive constant ${\Phi }_c$ is a time-coordinate normalization, not a second physical parameter. Matching to the Schwarzschild metric with time normalized at infinity fixes it. Therefore the regular stellar solutions, modulo time rescaling, have one physical parameter: ${\rho }_c$. This statement uses a suitably smooth fixed equation of state and refers to central densities for which the solution actually forms a regular star; not every central energy density need give a finite-radius configuration.

The spherical stellar structure equations apply to a planet as well, with its appropriate equation of state:The last equation is the local first law, including gravitational contraction through the specific entropy change. For an isolated non-burning planet, ${ϵ}_{\mathrm{nuc}}\simeq 0$, soIf the convective interior has nearly uniform specific entropy $s(t)$, positive outgoing luminosity gives $\dot{s}=-L_{\mathrm{int}}/\int Tdm<0$. This entropy loss of an isolated convective planet describes its secular cooling.

During assembly, the gravitational energy scale $G{M}^2/R$ provides a source of heat. Integrating hydrostatic equilibrium gives the virial theorem, $\Omega +3\int PdV=0$, neglecting the surface-pressure term. For a monatomic ideal gas this becomes $2U+\Omega =0$, hence $E=U+\Omega =\Omega /2$. Contraction makes $\Omega$ more negative, releases radiation, and heats part of the gas. A characteristic formation temperature has scale $k_{B}T\sim \mu m_{H}GM/R$ up to structure factors. The retained fraction depends on the accretion shock, producing hot and cold starts of a giant planet.

Formation releases heat, and an isolated planet subsequently loses specific entropy and intrinsic luminosity as it approaches a cooler, more degenerate state. The structure equations alone do not fix one initial temperature, and the early contraction of an ideal gas can increase its central temperature while its total energy decreases. “Cools” therefore need not mean that every depth becomes colder at every instant. Once degeneracy limits contraction, the thermal reservoir declines more directly; the emergent effective temperature usually falls along a fixed-mass cooling track.

This makes young giant planets attractive for exoplanet direct imaging: their thermal near/mid-infrared emission can be far brighter than that of old planets of the same mass. Imaging and spectroscopy measure projected separation, astrometry, luminosity, an exoplanet emission spectrum, temperatures and atmospheric molecular and exoplanet cloud signatures. Radius follows from luminosity and fitted temperature; mass and initial specific entropy generally require age/evolutionary information, or an independent dynamical mass measurement. A detected luminosity alone is not a model-independent planet mass.

The four equations describe, respectively, vertical hydrostatic equilibrium in the central object's gravity, the vertical divergence of radiative flux balancing viscous dissipation, optically thick radiative diffusion, and the equation of state combining ideal gas gas pressure with radiation pressure. Here $\mu$ is dynamic viscosity, not the mean molecular weight ${\mu }_m$; $\kappa$ is the opacity, $\sigma$ the Stefan-Boltzmann constant, and $\Omega$ the local orbital frequency.

In the radiation-pressure limit, $p=4\sigma T^4/(3c)$, so radiative diffusion givesCompare this with vertical hydrostatic equilibrium. Wherever $\rho \ne 0$,At fixed radius, constant opacity makes $dF/dz=c{\Omega }^2/\kappa$. The viscous dissipation equation therefore yieldsindependent of height within the optically thick radiation-supported layer. The result assumes the supplied gravity, heating and diffusion equations; optically thin surface corrections are outside that approximation.

Introduce the thermal expansion coefficient ${\alpha }_T=-{\rho }_0^{-1}(\partial \rho /\partial T)$, so the linear equation of state is $\rho ={\rho }_0[1-{\alpha }_T(T-T_0)]$. With gravity $-g\hat{\mathbf{z}}$, the resting conductive state of Rayleigh-Bénard convection satisfies the steady heat equation and hydrostatic pressure balance. Consequentlyand, up to an arbitrary constant,Indeed $P_b^{\prime }=-g{\rho }_b$, where ${\rho }_b={\rho }_0[1-{\alpha }_T\Delta T(1-z/d)]$.

The Boussinesq approximation retains temperature-dependent mass density in buoyancy while replacing it by ${\rho }_0$ in inertial coefficients; it requires $|{\alpha }_T\Delta T|\ll 1$. Subtracting the conductive state of Rayleigh-Bénard convection and dropping products of perturbations gives the dimensional Linearized Boussinesq equationsThe minus sign in the temperature equation comes from ${\mathbf{u}}_{\ast }\cdot {\nabla }_{\ast }{T}_b=-w_{\ast }\Delta T/d$.

Use the thermal-diffusion scaling of a convection layerThe Rayleigh number and Prandtl number areThus the nondimensional perturbation equations areHere $\nu$ is kinematic viscosity, $\kappa$ is thermal diffusivity, and positive $\mathrm{Ra}$ denotes destabilizing heating from below for ${\alpha }_T>0$.

For a quasistatic spherical star with fixed total mass, the stellar structure equations in radius areHere $s$ is the specific entropy of a moving mass shell, and the entropy term accounts for release or storage of thermal and gravitational energy. In thermal equilibrium it vanishes. An equation of state, opacity, energy-generation law and composition evolution close the equations. The absence of mass loss fixes the outer mass coordinate; it does not require the star to be in thermal equilibrium.

The radiative temperature gradient needed to carry the full luminosity isHere the actual temperature gradient is $\nabla =d\log T/d\log P$. For uniform composition the Schwarzschild criterion says that a radiative layer is stable if ${\nabla }_{\mathrm{rad}}<{\nabla }_{\mathrm{ad}}$, with marginal stability at equality, where the adiabatic temperature gradient is $(\partial \log T/\partial \log P{)}_s$. Use $\nabla ={\nabla }_{\mathrm{rad}}$ there. If ${\nabla }_{\mathrm{rad}}>{\nabla }_{\mathrm{ad}}$, convection carries some of the flux. Efficient convection gives $\nabla \simeq {\nabla }_{\mathrm{ad}}$; inefficient surface convection needs a transport prescription and can be superadiabatic. For composition gradients, use the Ledoux criterion, with threshold ${\nabla }_{\mathrm{ad}}+\varphi {\nabla }_{{\mu }_{\mathrm{mol}}}/{\delta }_{\mathrm{th}}$ instead. Here ${\nabla }_{{\mu }_{\mathrm{mol}}}=d\log {\mu }_{\mathrm{mol}}/d\log P$, $\varphi =(\partial \log \rho /\partial \log {\mu }_{\mathrm{mol}}{)}_{P,T}$ and ${\delta }_{\mathrm{th}}=-(\partial \log \rho /\partial \log T{)}_{P,{\mu }_{\mathrm{mol}}}$, using the mean molecular weight ${\mu }_{\mathrm{mol}}$.

Above the thin burning shell, take constant luminosity $L$, constant mean molecular weight ${\mu }_{\mathrm{mol}}$, ideal-gas pressure $P=\mathcal{R}T\rho$ with $\mathcal{R}=k_B/({\mu }_{\mathrm{mol}}{m}_u)$, and the Kramers opacity law $\kappa ={\kappa }_0\rho T^{-7/2}$. Then the temperature equation isTo implement the printed assumption that $P,\rho$ and $M_r$ are all power laws while retaining mass conservation, write $M_r\propto r^m$, $\rho \propto r^b$, $P\propto r^p$ and $T\propto r^d$. The mass equation gives $m=b+3$; the ideal gas law gives $d=p-b$; hydrostatic equilibrium gives $d=m-1$. The diffusion equation gives $d-1=2b-2-13d/2$. Solving this linear system gives the self-gravitating Kramers power-law envelope:The coefficients can also be matched consistently. If $r_c$ is the core boundary and $M_c$ its enclosed mass,Radiative transfer fixes the constant luminosity throughwhose right-hand side is independent of radius for these exponents. Thus this is a solution of all four envelope equations, not just a dimensional estimate.

Matching the boundary temperature to the isothermal core and extrapolating the idealized power law to the specified photosphere givesFor reference $M_R/M_c=(R/r_c{)}^{1/11}=2500^{1/10}\simeq 2.19$, so the envelope mass is not negligible in this particular solution. The given $M_c\simeq M_{\odot }$ fixes the core-radius normalization $r_c=(11/42)G{M}_c/(\mathcal{R}10^7\mathrm{K})$; a numerical absolute radius needs the unspecified molecular weight. The exponent ratio also gives $\nabla =d/p=5/21<2/5$, so the fully ionized monatomic version of this idealization is convectively stable.

A different common approximation neglects the envelope's self-gravity and sets $M_r\simeq M_c$ in the force equation. It cannot obey the exact mass equation with a nonzero density and exactly constant $M_r$. Under that additional approximation, the same calculation instead givesThis is the constant-core-mass power-law opacity radiative envelope limit, not the full power-law solution with changing enclosed mass. Stating which approximation is used resolves the otherwise different numerical radius ratios.

A compact nearly isothermal core, a thin luminosity-producing shell and a much larger cool envelope describe a shell-burning red giant. The full power-law profile is a deliberately simple model. The Kramers opacity law and an ionized ideal gas are not reliable at a $4000\mathrm{K}$ photosphere; partial ionization, other opacity sources and a convective envelope usually modify real cool giants. The quoted radius ratios are extrapolations within the stipulated model.

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.

For $N$ particles of a classical monatomic ideal gas, the equation of state and internal energy areOnly translational degrees of freedom contribute; the particle number is fixed throughout the thermodynamic cycle.

Substitute $v=3b\bar{v}$, $T=(8a/(27bk))\bar{T}$, and $p=(a/(27b^2))\bar{p}$ into the equation of state. After cancellation of the material parameters,This is the Van der Waals law of corresponding states: fluids described by this model share the same reduced equation, although real fluids need not obey it exactly.

The liquid-gas critical point is the endpoint of the gas-liquid coexistence curve where the distinct phases merge. For this equation of state it is the horizontal inflection of an isotherm, so $p_v=p_{vv}=0$:Dividing these conditions gives $v_c=3b$, then substitution givesAll lie in the physical domain $v>b$.

For a thin spherical shell, the radial pressure force per unit volume is $-dP/dr$. By the shell theorem, the inward gravitational acceleration is $GM(r)/r^2$, because exterior spherical shells contribute no force. Balancing force densities in hydrostatic equilibrium givesThe second equation is mass conservation; it closes the pressure-support equation once an equation of state is supplied.

The exact mixture equation of state is $\rho =(1-\varphi ){\rho }_f+\varphi {\rho }_p$. Keeping first-order terms in the small particle volume fraction and thermal expansion givesThe omitted term is ${\rho }_0\varphi \theta$. The downward solid-volume particle deposition flux is $J_v=W_s\varphi$; the corresponding particle mass flux is $J_m={\rho }_p{W}_s\varphi$. Uniform vertical mixing and constant layer depth implyFor $W_s>0$, set $r=W_s/H$. Thenand henceThe heated particle-laden layer loses its stable density stratification when this contrast vanishes. For $\beta >0$, solving for the neutral-buoyancy time givesThe layer is neutral at $T_s$ and statically unstable for $t>T_s$. The subsequent uniformly mixed lower-layer solution cannot represent the resulting overturning.

If $W_s=0$, then $\varphi ={\varphi }_0$, $\theta =\beta {\varphi }_{0}t$, and the continuous limit is $T_s=\gamma /\beta$. If $\beta =0$, the layer stays denser than its surroundings at every finite time and becomes neutral only asymptotically when $W_s>0$; thus $T_s=\infty$. These limiting cases must replace the printed formula when its denominator vanishes.

Use the shallow-water approximation, a hydrostatic vertically mixed layer, a deep motionless ambient, and no ambient fluid entrainment. Particle loss changes reduced gravity but changes total layer volume only at the discarded dilute-particle order. Likewise thermal expansion changes the equation of state while the leading Boussinesq approximation retains volume conservation. DefineDepth integration of mass conservation, horizontal momentum balance, particle transport, and the heating law givesHere the integrated excess hydrostatic pressure is $b{h}^2/2$. There is no particle source from the bed. Bed drag and the particle contribution to inertia are neglected at this order. The equivalent variable-buoyancy shallow water equations areThe term $h{b}_x/2$ is necessary when heating or particle loss makes the reduced gravity vary horizontally.

For $v=(h,u,\varphi ,\theta {)}^T$, the quasilinear system $v_t+A(v)v_x=S(v)$ hasIts characteristic polynomial factors asThus the characteristic curves have slopesOn the two repeated contact characteristic curves, the particle and heating laws are $D\varphi =-W_s\varphi /h$ and $D\theta =\beta \varphi$. Writing $c_g=\sqrt{bh}$ and $D_{\pm }={\partial }_t+(u\pm c_g){\partial }_x$, the gravity-wave compatibility relations areFor $h>0$ and $b>0$, gravity-wave right eigenvectors can be taken as $(1,\pm c_g/h,0,0{)}^T$. The contact eigenspace has dimension two: $\delta u=0$ and $b\delta h+(h/2)\delta b=0$, with independent $\delta \varphi ,\delta \theta$. Therefore $A$ has a full set of real eigenvectors.

The system is hyperbolic throughout strict static stability, but not strictly hyperbolic because the contact speed is repeated. At $b=0$ the eigenvalues coalesce and this diagonalizability is lost; when $b<0$ the gravity-wave speeds become complex, reflecting the loss of a statically stable lower layer.