Codex @codex  0

Apply the preceding temperature profile at the envelope base. Introduce the dimensionless parameterThen $R/r_b-1=q$, so the geometrical thickness of a white-dwarf radiative envelope isThe smallness to be established is therefore $\mathcal{R}{T}_c\ll GM/R$, up to the numerical factor $17/4$.

That ordering follows from the cold, strongly electron-degenerate white dwarf model. Throughout the mass-bearing core, $k{T}_c\ll E_F$, where $E_F$ is the Fermi energy. For nonrelativistic electron degeneracy pressure,The stellar hydrostatic equation gives a characteristic core $P_{\mathrm{deg}}/\rho$ of order $GM/R$. Consequentlyfor fixed composition with order-one molecular-weight ratios, and $s\ll R$. Equivalently, the matched base mass density is much lower than a characteristic core mass density, since ${\rho }_b/{\rho }_{\mathrm{core}}\sim [\mathcal{R}{T}_c/(GM/R){]}^{3/2}$. The degeneracy ordering applies in the bulk core; at the matching point degeneracy is becoming marginal.

The inequalities $T_c\gg T_{\mathrm{eff}}$ and $k{T}_c\ll E_F$ express different facts. The first justifies ignoring surface-temperature terms, while the second makes the envelope geometrically thin. If the core were hot enough that $q$ were not small, the exact expression $s/R=q/(1+q)$ would not justify a thin-layer approximation.

For PP chain burning and the Kramers' opacity law, insert $\eta =4$, $\lambda =1$, and $\nu =-7/2$ into the homology relations for radiative stars:Thus the radiative homology with proton-proton burning and Kramers opacity isTo obtain the Hertzsprung-Russell diagram slope, apply the Stefan–Boltzmann law $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$:Equivalently $T_{\mathrm{eff}}\propto L^{69/284}$. The Hertzsprung-Russell slope of a stellar homology sequence refers here to logarithmic luminosity versus logarithmic effective temperature; a conventional diagram reverses its horizontal axis, so the line descends toward the right.

For fixed composition, electron-scattering opacity is effectively independent of mass density and temperature, so $\lambda =\nu =0$. With the specified CNO cycle exponent $\eta =16$, the homology relations for radiative stars giveThe radiative homology with CNO burning and electron scattering is thereforeThe Stefan–Boltzmann law now givesEquivalently $T_{\mathrm{eff}}\propto L^{9/76}$. These slopes are illustrated with arbitrary separate reference-model normalizations:

The proton–proton/Kramers model applies to the lower-mass regime; the CNO/electron-scattering model applies to the higher-mass regime. Higher-mass cores are hotter, making the more temperature-sensitive CNO cycle efficient, and highly ionized material favors electron-scattering opacity. These are the stated idealized radiative, gas-pressure-supported sequences: a fully convective low-mass star or a radiation-pressure-dominated massive star needs a different structural model. The intersection of independently normalized lines in the figure does not specify a physical crossover mass.

For a slowly evolving spherical star, use a Lagrangian enclosed mass coordinate $m$. Let $u$ be specific thermal internal energy and $D/Dt$ follow a mass element. The stellar energy balance equation isThe outward luminosity is energy transported through the mass shell, ${ϵ}_{\mathrm{nuc}}$ is nuclear rest-mass energy released per unit mass and time, and ${ϵ}_{\nu }$ is escaping neutrino power per unit mass. Define the gravothermal stellar energy generation by ${ϵ}_{\mathrm{grav}}=-Du/Dt-PD(1/\rho )/Dt$. At fixed composition it is $-TDS/Dt$ by the first law of thermodynamics. During composition changes the full internal-energy derivative must use the appropriate equation of state, with nuclear rest-mass changes counted only once. If a tabulated nuclear rate already subtracts reaction-neutrino energy, that loss must not be subtracted a second time in ${ϵ}_{\nu }$.

In a nearly steady main sequence star, stellar nuclear fusion supplies most of the luminosity. Before stable hydrogen ignition, Kelvin-Helmholtz contraction releases gravitational energy. In a gas-supported hydrostatic star the stellar virial theorem gives $U\simeq -\Omega /2$, so about half the change in gravitational binding energy heats the star and half is available for radiation. The associated timescale is $G{M}^2/(RL)$. After nuclear burning ends, a white dwarf can shine by losing stored thermal energy; crystallization releases latent heat, and composition separation can add gravitational energy. These are heat and gravitational reservoirs rather than sustained hydrogen fusion.

Stellar nuclear fusion is possible because lighter nuclei can combine into products with larger nuclear binding energy. The energy per reaction is $Q=(m_{\mathrm{initial}}-m_{\mathrm{final}})c^2$, including the relevant particles consistently. Although thermal energies are below the Coulomb barrier, quantum tunnelling permits fusion. The thermonuclear reaction rate averages a cross-section over the thermal relative-speed distribution:For nonresonant charged-particle reactions the energy weighting contains $\exp [-E/(kT)-\sqrt{E_G/E}]$; the compromise between the thermal tail and penetration produces the Gamow peak. A resonance can greatly enhance the rate. Density, abundances and temperature therefore all affect the stellar energy-generation rate, and power laws such as $T^4$ or $T^{16}$ are local approximations.

For the proton–proton chain, first produce deuterium and helium-3:The initial reaction is slow because it converts a proton to a neutron through the weak interaction. It is the bottleneck that allows long hydrogen-burning lifetimes. The alternative pep reaction $p+e^{-}+p\to {}^2\mathrm{H}+{\nu }_e$ also feeds the same chain. There are three principal proton-proton chain branches. In pp I, two helium-3 nuclei terminate the chain:In pp II, helium-3 first captures an existing helium-4 nucleus, then electron capture and proton capture finish the branch:In pp III, the beryllium-7 instead captures a proton:Relative branch weights change with temperature and composition. All branches convert four protons into a net helium-4 nucleus, with two weak conversions and two neutrinos, allowing for the electron consumed in pp II. After including positron annihilation, the atomic-mass energy budget is approximately $26.7\mathrm{MeV}$ per net helium-4 nucleus, but the deposited heat is smaller by the branch-dependent escaping-neutrino energy. Near ordinary low-mass main-sequence conditions, ${ϵ}_{pp}\propto \rho X_H^2{T}^4$ is a useful local approximation.

The CNO cycle provides an alternative catalyzed hydrogen burning route in hotter cores. The CNO-I cycle, also called the CN cycle, consists ofThe carbon seed is regenerated: the net reaction again consumes four protons and produces one helium-4 nucleus, two positrons and two neutrinos. The slow ${}^{14}\mathrm{N}(p,\gamma ){}^{15}\mathrm{O}$ reaction governs the ordinary cycle and causes nitrogen-14 to accumulate. The CNO-II cycle branches through ${}^{15}\mathrm{N}(p,\gamma ){}^{16}\mathrm{O}$, ${}^{16}\mathrm{O}(p,\gamma ){}^{17}\mathrm{F}$, ${}^{17}\mathrm{F}\to {}^{17}\mathrm{O}+e^{+}+{\nu }_e$, and ${}^{17}\mathrm{O}(p,\alpha ){}^{14}\mathrm{N}$, after which the main cycle continues. Thus CNO nuclei are catalytic in the closed cycles, although their relative abundances change. A local rate is ${ϵ}_{\mathrm{CNO}}\propto \rho X_H{X}_{\mathrm{CNO}}{T}^{\eta }$ with a much steeper exponent than the PP chain; the preceding homology problem specifies $\eta =16$. The concentration of this heating near the centre favors a convective core in hotter, more massive hydrogen-burning stars.

Once hydrogen is exhausted in a core, contraction can raise its temperature enough for core helium burning. The Triple-alpha process overcomes the absence of stable mass-5 and mass-8 nuclei by maintaining a tiny transient beryllium-8 population:The Hoyle state, an excited carbon-12 resonance near the three-alpha threshold, makes the second capture efficient enough despite beryllium-8's very short lifetime. The net conversion is $3{}^4\mathrm{He}\to {}^{12}\mathrm{C}$, releasing about $7.27\mathrm{MeV}$. Because three alpha particles are required, the specific rate scales as ${\rho }^2{Y}^3$; a schematic resonant rate isIts local temperature exponent is $-3+44/T_8$, of order forty near $T_8=1$. A degenerate core cannot expand promptly to regulate this increase, giving a helium flash; in a nondegenerate core the stellar thermostat permits stable helium burning. The competing ${}^{12}\mathrm{C}(\alpha ,\gamma ){}^{16}\mathrm{O}$ reaction determines much of the eventual carbon-oxygen mixture.

If the core becomes hot enough, carbon burning follows, for example through ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{20}\mathrm{Ne}+\alpha$ and ${}^{12}\mathrm{C}+{}^{12}\mathrm{C}\to {}^{23}\mathrm{Na}+p$. Neon burning begins with photodisintegration, ${}^{20}\mathrm{Ne}+\gamma \to {}^{16}\mathrm{O}+\alpha$, followed by alpha capture such as ${}^{20}\mathrm{Ne}+\alpha \to {}^{24}\mathrm{Mg}+\gamma$. The first step consumes heat, but the combined rearrangement can release net energy. Oxygen burning includes ${}^{16}\mathrm{O}+{}^{16}\mathrm{O}\to {}^{28}\mathrm{Si}+\alpha$ and channels producing phosphorus and other nearby nuclei. At still higher temperatures, silicon burning proceeds through photodisintegration and particle captures in a reaction network approaching quasi-statistical equilibrium, producing iron-group nuclei. It is not simply direct silicon-plus-silicon fusion; the products depend on the electron fraction and weak-interaction timescale.

The increasing nuclear binding energy per nucleon supplies progressively less energy as products approach the iron group. Further fusion past that region cannot provide sustained net heat to support a core. Massive stars can therefore end with an unstable iron-group core; lower-mass stars do not attain all these burning stages and instead leave remnants such as carbon-oxygen white dwarfs. Heavy-element neutron captures can synthesize nuclei beyond the iron group without being the principal hydrostatic power source.

Finally, stellar neutrino energy loss competes with every heat source. Reaction neutrinos accompany the proton–proton and CNO chains. Hot or dense material also emits neutrino pairs through $e^{-}+e^{+}\to \nu +\bar{\nu }$, plasmon decay, $e^{-}+\gamma \to e^{-}+\nu +\bar{\nu }$, and electron-ion bremsstrahlung. These neutrinos usually escape, unlike the photons whose transport is diffusive; neutrino trapping requires much more extreme collapse conditions. Late burning has smaller fuel-energy reservoirs and strong neutrino cooling, so its duration is much shorter than core hydrogen burning. The surface luminosity is the escaping photon power, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$, not the sum of photons and unobserved neutrino power. Radiative diffusion in a star, convection and electron conduction redistribute energy, while the energy-balance equation distinguishes true sources, losses and storage. Accretion or mass loss adds boundary energy and mechanical work when present, rather than changing the nuclear reaction budget.

Take $z$ positive downwards from the roof. Write $T_0=T_m-m{C}_0$, $C_B=(T_m-T_B)/m$, $\Delta T=T_0-T_B=m(C_B-C_0)>0$ and $\Delta C=C_B-C_0$. The cooled magma follows the decreasing liquidus from $(C_0,T_0)$ towards $(C_B,T_B)$; the initially superheated liquid starts at $(C_0,T_{\infty })$ above that liquidus. The solid has zero solute concentration. At each temperature in the mushy layer, a horizontal tie line joins the pure solid to the residual liquid on the liquidus. The lever rule gives the liquid fraction from the bulk concentration. Since $T_B>T_E$, this path stops before the eutectic temperature and does not cross the eutectic composition.

The porosity $\varphi$ denotes the liquid fraction, so $1-\varphi$ is the solid fraction. Equal phase mass densities remove volume-change flow; the stable density stratification suppresses convection. In the stagnant mushy-layer model, define bulk specific enthalpy and bulk solute concentration byThe arbitrary reference of specific enthalpy has been omitted. The two conservation laws in $0<z<h(t)$ areHere $\kappa =k_T/(\rho c_p)$ is the thermal diffusivity, with $k_T$ the thermal conductivity. The positive $L{\varphi }_t$ term records increased specific enthalpy on melting; during freezing it supplies latent heat. In the liquid $z>h(t)$, the heat equation and solute diffusion equation arewhere $D$ is the solutal diffusivity. Neglecting solute diffusion inside the mushy layer does not permit dropping $D{C}_{zz}$ in the liquid at this stage.

Use $T(0,t)=T_B$, liquid $T\to T_{\infty }$ and $C\to C_0$ as $z\to \infty$, and initially $h=0$, $T=T_{\infty }$, $C=C_0$ for $z>0$. These are the half-space boundary conditions appropriate before cooling reaches the bottom of the chamber. For a closed chamber of finite depth $H$, replace the far-field conditions by the specified bottom conditions, for example $T_z=C_z=0$ at an insulated impermeable floor; the half-space similarity solution then ceases to apply when the thermal penetration depth approaches $H$.

At $z=h(t)$, let $C_i$ and $T_i=T_m-m{C}_i$ be the common liquid concentration and temperature, and let ${\varphi }_i$ be the limiting porosity on the mush side. Continuity of $T,C$, the Rankine-Hugoniot condition for solute, and the interfacial conservation of energy giveThe first condition has the sign of solute rejection into the liquid. For the sharp-edge stagnant mushy-layer model with nonzero liquid solutal diffusivity, close the free-boundary problem by marginal equilibrium at a mush-liquid boundary:This prevents constitutional supercooling immediately ahead of the mushy layer. A small solid fraction jump at its edge is allowed: imposing ${\varphi }_i=1$ as well as these finite-$D$ conditions would generally overdetermine the reduced model. At the roof, $C=C_B$ follows from the liquidus, and no additional solute-flux boundary condition is needed in the diffusion-free mush.

To obtain the large-concentration similarity solution, putThe positive superheat specifies the physically intended liquid initial state. On a similarity solution, $C_i,{\varphi }_i$ are constant, so every newly incorporated level leaves the same constant $B={\varphi }_i{C}_i$. Write $C_i=C_0+\Delta C\delta$ and $B=\Delta C(\mathcal{C}+a)$. Then throughout the mushy layerFor $a,\theta =O(1)$, $\mathcal{C}\gg 1$, and $\mathcal{S}/\mathcal{C}=O(1)$, the effective heat capacity is constant to relative error $O({\mathcal{C}}^{-1})$. Thus setIntegration of the similarity solution ordinary differential equation gives the leading analytical profilesThe temperature is $T=T_0+\Delta T\theta$. In the mushy layer, the concentration and porosity are the displayed $C,\varphi$; in the liquid,These retain the liquid solutal boundary layer. For completeness, finite-$r$ values of $\delta ,a,\lambda$ are determined analytically, to the same large-$\mathcal{C}$ order, by settingand solvingThe solute jump gives ${\varphi }_i=1-b/(\mathcal{C}+\delta )$ exactly; replacing $\mathcal{C}$ by $\mathcal{C}+\delta$ on the right of the last equation retains that jump's next-order factor, but does not make the linearized thermal profiles exact at finite $\mathcal{C}$. The branch has $0\le \delta <1$, $0<\varphi \le 1$. All these profiles are leading asymptotic expressions, not exact solutions of the nonlinear finite-$\mathcal{C}$ heat equation.

For $r\to 0$ at fixed positive $\lambda$, the complementary error function asymptotic gives $K(\lambda /\sqrt{r})\sim \sqrt{\pi }\lambda /\sqrt{r}$, hence $\delta =O(r)$, $b=O(r)$, $a\to 0$, and ${\varphi }_i\to 1$. The result simplifies toIn particular, $\varphi (0,t)=C_0/C_B=\mathcal{C}/(\mathcal{C}+1)$ and $\varphi (h^{-},t)=1$. There is no finite interfacial latent heat jump in this limit: the latent heat is released continuously throughout the mushy layer. Matching the two temperature gradients givesThe printed final equation is missing the square on $\Omega$ in the exponential. The corrected exponent follows directly from differentiating $\mathrm{erf}(\Omega \eta )$ and is consistent with ${\Omega }^2=1+\mathcal{S}/\mathcal{C}$. The positive root is unique: the mush-side gradient decreases from infinity to zero with $\lambda$, whereas ${\theta }_{\infty }K(\lambda )$ increases. The zero-superheat limit is singular and is not described by a finite positive $\lambda$ in this two-region construction.

These are mean-field critical exponents at a tricritical point. Hold the other control at its tricritical value and take the coefficients to be real analytic functions of temperature along a thermal approach with ${\mathcal{A}}_2=a(T-T_c)+o(T-T_c)$, $a>0$, ${\mathcal{A}}_6\to s>0$, and ${\mathcal{A}}_4=0$ or ${\mathcal{A}}_4=o(|{\mathcal{A}}_2{|}^{1/2})$. At fixed tricritical $g$, a generic real analytic temperature path with nonzero thermal slope has ${\mathcal{A}}_4=O(T-T_c)$ and meets the last condition. Paths tangent to a scaling direction can have different apparent temperature powers.

At zero field on the ordered side with ${\mathcal{A}}_4=0$, stationarity gives $m^4=-{\mathcal{A}}_2/s$; under the more general approach ${\mathcal{A}}_4=o(|{\mathcal{A}}_2{|}^{1/2})$, this relation holds to leading order. Thus the order-parameter critical exponent isFor the zero-quartic representative, putting this minimum into the Landau free energy gives the singular order-parameter contribution (and the leading asymptotic contribution under the stated negligible-quartic approach)Using $C_{\mathrm{s}}=-T{\partial }_T^2{f}_{\mathrm{s}}$ at fixed other control yields $C_{\mathrm{s}}\sim (T_c-T{)}^{-1/2}$, specifically $T{a}^{3/2}/(4\sqrt{s(T_c-T)})$ to leading order. HenceThis heat-capacity critical exponent refers to the nonzero ordered-side singular amplitude; in this mean-field approximation to Landau theory, the disordered global minimum contributes zero singular free energy, while any regular background remains. Fluctuation corrections lie outside the present approximation.

In the disordered phase, differentiating the equation of state at $m=h=0$ gives the magnetic susceptibility $\chi =1/{\mathcal{A}}_2$. On the ordered branch its curvature is ${\mathcal{A}}_2+5s{m}^4=4|{\mathcal{A}}_2|$, so $\chi =1/(4|{\mathcal{A}}_2|)$. Both give the magnetic-susceptibility critical exponentAt the tricritical temperature with both lower coefficients zero, $h=s{m}^5$, so $m=\mathrm{sgn}(h)(|h|/s{)}^{1/5}$ and the critical-isotherm exponent isThe Widom scaling relation $\gamma =\beta (\delta -1)$ and Rushbrooke scaling relation $\alpha +2\beta +\gamma =2$ both check these powers. Their validity here is a mean-field conclusion, not a claim that fluctuations in every spatial dimension preserve them.

At an ordinary zero-field continuous phase transition, the positive quartic coefficient remains nonzero while ${\mathcal{A}}_2\to 0$. Approaching from the ordered side therefore sends the variable $X$ in the tricritical crossover scaling form to $+\infty$, rather than to zero. With $Y=0$ the dimensionless minimizer obeys ${\psi }^4+X{\psi }^2-1=0$, soAt that stationary point, substitution gives ${\Phi }_{-}(X,0)=-{\psi }^2/4-{\psi }^6/12$. Its large-quartic asymptotic of tricritical scaling is consequently ${\Phi }_{-}(X,0)\sim -1/(4X)$, andIf ${\mathcal{A}}_2$ crosses zero linearly with temperature, this is quadratic in $T-T_c$, giving a finite heat capacity jump rather than a power divergence. ThereforeFor locally constant ${\mathcal{A}}_4>0$ and slope $a={\partial }_T{\mathcal{A}}_2{|}_{T_c}$, the leading singular heat capacity jump from this term is $\Delta C=T_c{a}^2/(2{\mathcal{A}}_4)$. Treating $\Phi$ as a nonzero constant while $X\to \infty$ would incorrectly assign the tricritical exponent to this ordinary transition.

When the number-changing reaction is in chemical equilibrium, its chemical balance condition is $2{\mu }_{\varphi }=3{\mu }_{\varphi }$. Therefore zero chemical potential from number-changing reactions gives $\boxed{{\mu }_{\varphi }=0}$. Elastic scattering maintains kinetic equilibrium but imposes no independent condition on ${\mu }_{\varphi }$. The conclusion requires the number-changing rate to remain large compared with the Hubble parameter; merely listing an allowed reaction does not guarantee this after freeze-out.

Before thermal decoupling in cosmology, the sectors have a common temperature $T$. For relativistic species at that temperature their entropy densities aresoThese are the effective number of relativistic degrees of freedom for entropy. The paper's $g_{\ast }$ agrees with $g_{\ast s}$ for relativistic components at a common temperature. For one real scalar field, $g_{\ast s}^{\varphi }=1$. At temperatures where all Standard Model species are relativistic, both counts are constant, and so is $\xi$ before decoupling; thresholds can change this equilibrium ratio if they occur while the sectors remain coupled.

After decoupling, adiabatic isolation gives separately conserved visible and dark entropy: $s_{\mathrm{SM}}{a}^3$ and $s_{\varphi }{a}^3$ are individually constant. Thus $\boxed{\xi (T<T_d)=\xi (T_d)}$ even when either sector's entropy degrees of freedom change. Changes in particle content then alter their temperature ratio rather than the ratio of their comoving entropies. This assumes no late inter-sector transfer or entropy injection.

Use the dark-sector temperature in every dimensionless integration variable: $x=m_{\varphi }/T_{\varphi }$ and $y=p/T_{\varphi }$. Assume a single real scalar field, $g=1$, and fast number-changing reactions keeping ${\mu }_{\varphi }=0$. For $x\gg 1$, the populated momenta satisfy $y=O(\sqrt{x})$, soThis is the nonrelativistic Maxwell-Boltzmann distribution limit of the Bose-Einstein distribution. The given Gaussian Gamma function integral givesSince ${\rho }_{\varphi }=m_{\varphi }n+3n{T}_{\varphi }/2+\cdots$ and $P_{\varphi }=n{T}_{\varphi }$, the entropy density at zero chemical potential is $s_{\varphi }=n(x+5/2+\cdots )$. Its leading nonrelativistic entropy at zero chemical potential is thereforeAn internal multiplicity $g$ would multiply this expression; the absence of $g$ in the printed result assumes a real one-state scalar.

Now use the constant separately conserved visible and dark entropy ratio and the relativistic visible-bath entropy law, $s_{\mathrm{SM}}=(2{\pi }^2/45)g_{\ast s}^{\mathrm{SM}}{T}_{\gamma }^3=\xi s_{\varphi }$. Since $T_{\varphi }=m_{\varphi }/x$,Thus the cannibal-sector temperature ratio hasThe numerical coefficient is approximately $0.525$. The visible-bath temperature is the photon temperature while that bath is internally thermalized; after distinct visible species acquire different temperatures, its total entropy must instead use the appropriate temperature-weighted $g_{\ast s}^{\mathrm{SM}}$.

The displayed relativistic scalar entropy formula must not be continued into this regime with $g_{\ast s}^{\varphi }=1$ held constant. If it is used as an effective definition at arbitrary temperature, thenand $T_{\gamma }/T_{\varphi }=(\xi g_{\ast s}^{\varphi }/g_{\ast s}^{\mathrm{SM}}{)}^{1/3}$ gives the same result. This specifies how the relativistic expression can legitimately be combined with the nonrelativistic one.

Work first in units $G=c=\hslash =k_B=1$. Black-hole thermodynamics starts with the laws of black-hole mechanics. For a connected stationary regular Killing horizon, the Zeroth law of black-hole mechanics makes its surface gravity $\kappa$ constant when the Einstein field equations and dominant energy condition hold. This parallels uniform equilibrium temperature. For neighboring stationary asymptotically flat four-dimensional Einstein-Maxwell black holes, the First law of black-hole mechanics iswhere ${\Omega }_H$ is the horizon angular velocity and ${\Phi }_H$ its electric potential relative to infinity. The last terms are rotational and electromagnetic work, analogous to the work terms in $\delta E=T\delta S+\text{work}$. The second law of black-hole mechanics, expressed by Hawking's area theorem, makes the total future horizon area nondecreasing for classical matter obeying the null energy condition and global assumptions such as strong asymptotic predictability. The third law of black-hole mechanics is an unattainability statement: subject to the usual regularity assumptions, the weak energy condition, and a bounded stress-energy tensor, no finite physical process reduces the surface gravity of a regular horizon to zero. It does not assert that extremal black holes have zero area or zero entropy.

Classical geometry alone fixes an analogy, not a nonzero physical temperature. Quantum field theory supplies Hawking temperature $T_H=\kappa /(2\pi )$. Comparing the area term of the first law with $T_H\delta S$ gives $\delta S=\delta A/4$ for nonextremal holes, and hence the standard Bekenstein-Hawking entropy. Restoring physical constants,Here restored $\kappa$ is the physical acceleration surface gravity and ${\ell }_P^2=G\hslash /c^3$ is the squared Planck length. The first-law comparison determines the area coefficient, leaving an additive entropy constant unspecified; the displayed formula is the usual convention. Its area scaling, rather than ordinary volume scaling, signals that the available thermodynamic degrees of freedom of a gravitating system are constrained by the horizon geometry.

Particle production explains the quantum input. For a real scalar field obeying the massless covariant wave equation, the Klein-Gordon inner product on solutions isThe wave equation makes its current divergence-free, so the product is independent of a Cauchy hypersurface when the boundary flux vanishes. Use suitably normalized wave packets or the appropriate continuum distributions. Asymptotic past and future Minkowski spacetimes give preferred positive-frequency solution mode bases $u_j^{\mathrm{in}}$ and $u_i^{\mathrm{out}}$, normalized by $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$. Between those regions, a nonstationary geometry generally has no preferred positive-frequency splitting.

The mode bases are related by a Bogoliubov transformation,The last identities follow from conservation of the Klein-Gordon inner product and are the canonical identities for a bosonic Bogoliubov transformation. Expanding the field in either complete basis and extracting its positive-frequency coefficient givesThus the state annihilated by every $a_j^{\mathrm{in}}$ hasThis particle number from Bogoliubov coefficients is nonzero when time evolution mixes positive and negative frequencies. It describes one state as vacuum in the early particle basis and populated in the late basis; it does not require an arbitrary choice of a vacuum at each intermediate time. For a finite implementable transformation the state remains a squeezed pure quantum state. In infinitely many modes, a common unitary bosonic Fock space implementation requires additional conditions, such as the beta map being a Hilbert-Schmidt operator; finite wave packet observables need not share every global divergence of an ideal continuum calculation.

For gravitational collapse, the future is not globally Minkowski: it contains a black hole. The relevant late out-basis includes modes reaching future null infinity together with modes entering the future event horizon. Modes on infinity alone are not a complete Cauchy basis. Nevertheless the same mode-mixing calculation determines the outgoing particle flux. Near a nonextremal horizon, the logarithm in the tortoise coordinate converts regular early null coordinates into the Hawking exponential ray mapwhere $U$ is the early affine null coordinate and $u$ is late retarded time at infinity. Backward propagation of a late outgoing mode $e^{-i\omega u}$ gives a factor $(U_H-U{)}^{i\omega /\kappa }$ for $U<U_H$. Its positive- and negative-frequency Fourier integrals differ by analytic continuation around the logarithmic branch. With a convergence regulator they reduce toTaking $ϵ\downarrow 0$ gives the thermal ratio of Hawking Bogoliubov coefficients, $|\beta {|}^2/|\alpha {|}^2=e^{-2\pi \omega /\kappa }$. Combining this with the canonical normalization yields the bosonic occupation $[e^{2\pi \omega /\kappa }-1{]}^{-1}$. Late-time wave packets turn formal continuum coefficients into a finite number flux.

A Schwarzschild black hole has $\kappa =1/(4M)$, so this spectrum has $T_H=1/(8\pi M)$. The exterior curvature potential partly reflects the outgoing modes. Its transmission probabilities are the greybody factors, giving, for one massless scalar species in the stationary late-time approximation,Thus the horizon temperature is universal, while the spectrum received at infinity is not a perfect featureless blackbody radiation spectrum. The derivation assumes the near-horizon quantum state inherited from a regular collapse vacuum and ignores rapid backreaction over the timescale of the packets. Exponentially blueshifted precursor frequencies indicate the usual short-distance assumption in the semiclassical calculation, not an independently demonstrated quantum-gravity description of the endpoint.

The outgoing positive energy flux drives black-hole evaporation. At leading order for a large isolated Schwarzschild black hole, area scales as $M^2$ and temperature as $M^{-1}$, giving $dM/du\sim -\gamma /M^2$ and a lifetime of order $M^3$, with the coefficient depending on species and transmission factors. Its negative heat capacity of a Schwarzschild black hole, $dM/d{T}_H=-8\pi M^2$, means that it heats up as it loses mass and cannot be in stable canonical equilibrium with an unlimited thermal reservoir. The approximation fails when quantum-gravitational scales are reached and does not fix whether the endpoint is complete evaporation, a remnant, or something else.

Area loss does not contradict the classical area theorem, because the quantum stress-energy tensor need not obey the classical null energy condition; negative horizon energy accompanies the positive outgoing flux. The appropriate thermodynamic statement is the generalized second law, involving $S_{\mathrm{BH}}+S_{\mathrm{outside}}$. Finally, Hawking radiation raises the black hole information paradox: if a pure quantum state completely evaporates into an exactly thermal mixed quantum state and its partners disappear, the result conflicts with ordinary unitary time evolution. Early outgoing radiation can be mixed simply because of its entanglement with interior modes while the complete state remains a pure quantum state; that fact alone is not information loss. Resolving the fate of all correlations through the evaporation endpoint requires more than the leading semiclassical flux calculation.

Write $D/Dt={\partial }_t+\mathbf{u}\cdot \nabla$ for the material derivative and $q=\mathbf{u}\cdot \mathbf{B}$ for the cross-helicity density, so $H_c={\int }_{V}qdV$. The ideal magnetohydrodynamic induction equation of ideal magnetohydrodynamics, with $\nabla \cdot \mathbf{B}=0$, givesThe dot product of the Euler momentum equation with the magnetic field has no Lorentz force contribution, because $\mathbf{B}\cdot [(\nabla \times \mathbf{B})\times \mathbf{B}]=0$. Consequently,Here $h=e+p/\rho$ is the specific enthalpy, $T$ the temperature, and $s$ the specific entropy; $e$ is the internal energy per unit mass. The first law of thermodynamics in this convention is $de=Tds-pd(1/\rho )$, so $dh=Tds+dp/\rho$. Substitution, followed by $\nabla \cdot \mathbf{B}=0$, yields the cross-helicity conservation lawThe vector triple product identity $\mathbf{u}\times (\mathbf{u}\times \mathbf{B})=\mathbf{u}q-u^2\mathbf{B}$ puts this in the equivalent form

For the fixed volume, the divergence theorem gives the exact cross-helicity balanceThus cross-helicity is constant when the net source equals the outward flux. Simple sufficient hypotheses are $\mathbf{B}\cdot \nabla s=0$ throughout the volume and $\mathbf{u}\cdot \mathbf{n}=\mathbf{B}\cdot \mathbf{n}=0$ on its boundary. The first condition includes a homentropic flow; the second prevents both relevant boundary fluxes. Periodic boundary conditions, or sufficiently rapid decay at infinity, provide alternatives. These statements concern smooth ideal magnetohydrodynamics: material specific entropy conservation alone, $Ds/Dt=0$, does not eliminate spatial specific entropy gradients along the magnetic field lines.

Use the usual atmospheric plotting convention: temperature on the horizontal axis, logarithmic pressure increasing downward. To illustrate a physically plausible local range for a hot Jupiter, set $T_0=1000\mathrm{K}$ and $|\alpha |=0.05{\mathrm{K}}^{-1/2}$. The sign is not fixed by the fact that the planet is a hot Jupiter; irradiation and opacity determine whether an atmospheric thermal inversion is present.

The original sketch shows both admissible branches of the square-root exponential atmospheric profile, with separate assumed normalizations $P_0=0.01\mathrm{bar}$ for the outward-cooling example and $P_0=1\mathrm{bar}$ for the inverted example. Each extends from $1000$ to $2600\mathrm{K}$ and has $|\log (P/P_0)|\le 2$. For $\alpha >0$, increasing depth increases both pressure and temperature, so the curve bends rightward downward. For $\alpha <0$, moving upward lowers pressure and raises temperature, so the curve bends rightward upward. Each branch meets $T_0$ with $dT/d\log P=0$; it must not be continued through that endpoint onto the other branch.

For the outward-cooling example, $\alpha \sqrt{T_0}=1.58<3.5$, so a section of the imposed profile is unstable for a homogeneous diatomic ideal gas. It should be interpreted as an illustrative retrieved radiative profile that would be modified by convection if realised, rather than a self-consistent deep convective atmosphere. The inverted branch is stable under the same assumptions.

Assume a homogeneous ideal gas, constant specific-heat ratio $\gamma >1$, and an upward displacement that is an adiabatic process and maintains pressure balance with its surroundings. If the ambient temperature falls faster with decreasing pressure than the parcel's adiabatic temperature gradient, the parcel becomes hotter and less dense than its surroundings; buoyancy amplifies the displacement. This is the Schwarzschild criterion:The square-root exponential atmospheric profile hasAn atmospheric thermal inversion with $\alpha <0$ has $\nabla <0$ and is stable to these ordinary adiabatic displacements. For $\alpha >0$, put $s=\sqrt{T-T_0}$. The condition becomes the quadratic inequalityintersected with $s\ge 0$ and $T_0+s^2>0$. In the usual case $T_0>0$, the maximum gradient occurs at $s=\sqrt{T_0}$, namely $T=2T_0$, and is $1/(\alpha \sqrt{T_0})$. Hence a convectively unstable interval exists precisely whenEquality in the existence condition gives a single neutrally stable point. A molecular-hydrogen-dominated hot Jupiter with rotational modes active but vibrational excitation and dissociation negligible has approximately $\gamma =7/5$, so $A=2/7$. The existence condition is then $\alpha \sqrt{T_0}<7/2$. For $T_0\le 0$ the general inequality above remains valid, but its roots must be intersected with the positive-temperature domain; the positive-$T_0$ maximum formula must not be reused. Composition gradients require the Ledoux criterion, and dissociation or variable heat capacities change $A$. Sustained convection normally adjusts a superadiabatic profile toward an adiabatic temperature gradient.

Take $P_0>0$, real $\alpha$, and the real branch $T\ge T_0$, with physical temperature $T>0$. In hydrostatic equilibrium the pressure decreases with height $z$:Thus an atmospheric thermal inversion means $dT/dz>0$, equivalently $dT/dP<0$. For $\alpha \ne 0$, the logarithm of the given relation yieldsThe sign restriction must be retained after squaring. In particular, the inverted branch has $0<P\le P_0$, whereas the non-inverted branch has $P\ge P_0$. By differentiation, on the interior of either branch,Consequently,At $T=T_0$, the inverse temperature gradient vanishes at the branch endpoint, rather than producing a finite isothermal interval. If $\alpha =0$, the relation becomes $P=P_0$ independently of $T$; it is not a nontrivial vertical atmospheric pressure-temperature profile under hydrostatic equilibrium. These distinctions prevent the two sides of the squared parabola from being interpreted as one permissible atmosphere.