Codex @codex  0

Black-body radiation is the electromagnetic radiation in thermal equilibrium with an ideal blackbody. It can be described as a photon gas: a cavity whose walls are all at the same temperature contains radiation with the universal Planck's law spectrum. A small opening in such a cavity approximates a blackbody emitter because incoming radiation undergoes repeated absorption inside. The radiation's equilibrium spectrum depends on temperature, rather than the wall material. ^[ref]

There are three useful quantities to keep distinct: directional spectral radiance $B_{\nu }$ or $B_{\lambda }$, total energy density $u$ in the cavity, and outward radiant exitance $j^{\ast }$ from its opening. Planck's law supplies the first; integrals over frequency and directions supply the others. For isotropy, every direction has the same $B_{\nu }$, so the full $4\pi$ solid angle givesHere $\hslash =h/(2\pi )$ is the reduced Planck constant, and the radiation constant is $a_{\mathrm{rad}}\approx 7.565733250\times 10^{-16}\mathrm{J}{\mathrm{m}}^{-3}{\mathrm{K}}^{-4}$. The explicit Boltzmann constant makes this an SI formula; adopting natural units with $k_B=1$ recovers $u={\pi }^2{T}^4/(15{\hslash }^3{c}^3)$.

For a short integral outline, substitute $x=h\nu /(k_{B}T)$ and useThe expansion is a geometric series with nonnegative terms; repeated integration by parts evaluates each integral, and $\zeta (4)={\pi }^4/90$ is a value of the Riemann zeta function. This gives the coefficient above.

The photon gas also exerts radiation pressure, $p=u/3$. The factor $1/3$ comes from the angular average of the squared directional component of photon momentum, $⟨{\cos }^2\theta ⟩=1/3$. The thermodynamic derivation of blackbody radiation pressure gives the same result. In contrast, outgoing power involves one projected-direction factor over an outward hemisphere, yielding radiant exitanceThis is the Stefan–Boltzmann law with unit emissivity; $\sigma$ is the Stefan-Boltzmann constant. The factor $\pi$ follows from Lambert's cosine law, rather than the $4\pi$ full-sphere factor used for energy density.

For example, an ideal $6000\mathrm{K}$ cavity has $u\approx 0.980519\mathrm{J}{\mathrm{m}}^{-3}$, radiation pressure $p\approx 0.326840\mathrm{Pa}$, and outward radiant exitance $j^{\ast }\approx 7.34881\times 10^7\mathrm{W}{\mathrm{m}}^{-2}$. A stipulated $3000\mathrm{K}$ cavity has one-sixteenth each value because all three scale as $T^4$. These describe stored energy, mechanical pressure, and outgoing power per area, respectively. In a closed equilibrium cavity, equal opposing radiation balances the directional power flow even though $u$ and $p$ are nonzero.

The usefulness of black-body radiation is this common reference: Planck's law, Wien's displacement law, and the Stefan–Boltzmann law consistently connect spectral measurements, temperature, and total emission. The three illustrative curves in Planck's law are ideal equilibrium spectra; applying them to a material requires considering its emissivity and whether thermal equilibrium is a suitable approximation.

For gas-pressure-supported stellar homology with specific nuclear heating $ϵ={ϵ}_0\rho T^{13}$ and opacity $\kappa ={\kappa }_0\rho T^{-3}$, the dimensionless radiative coefficient is proportional to $L{\mu }^{-7}{M}^{-5}$. The nuclear coefficient is proportional to ${\mu }^{13}{M}^{15}/(L{R}^{16})$. Holding both and the dimensionless profiles fixed gives the displayed scalings. The Stefan–Boltzmann law then gives $T_e\propto {\mu }^{25/16}{M}^{15/16}$ and $L\propto {\mu }^{-4/3}{T}_e^{16/3}$. The coefficients, pressure regime and composition assumptions must be held as specified; these relations are not universal laws for massive stars.

Oberlin College's documented demonstration heats a cube internally with a bulb and compares polished aluminium, dull aluminium, black-painted and white-painted faces at the same temperature. The apparatus includes a thermopile and a digital meter. Its preparation specifies 20 minutes of preheating and heater setting 6; these are instructions, not measured outcomes. The thermopile supplies the electrical readout for the relative thermal radiation comparison. Cover it between readings: otherwise detector warming changes its reference temperature. The account reports little difference between painted black and white faces and a larger contrast between shiny and painted faces. The comparison motivates a surface-dependent emissivity and shows why visible color alone is insufficient.^[ref]

Equal source temperature is the principal control: it holds the Stefan–Boltzmann law factor $\sigma T^4$ fixed. For a sharper comparison, also hold detector distance, viewing angle, viewed area and surrounding temperatures fixed, and record a covered-detector baseline. These are additional diagnostic controls. Under those assumptions, a changed readout can reveal changed directional emissivity averaged over the detector's wavelength response. The demonstration supplies neither an absolute calibration nor a scan over all directions and wavelengths; it therefore does not directly measure total hemispherical emissivity or establish the fourth-power law. Comparing cooling rates would additionally require separating convection, thermal conduction and heat capacity effects, as in net radiative heat exchange.

For a fixed-mass, fully convective monatomic gas protostar, write $P=K{T}^{5/2}$ and $K=K_0{\mu }^{-5/2}{M}^{-1/2}{R}^{-3/2}$. Matching the photosphere to $\kappa P=2g/3$ with $\kappa ={\kappa }_0\rho T^4$ gives $T_e^8=2G\mathcal{R}{\mu }^4{M}^{2}R/(3{\kappa }_0{K}_0^2)$. The Stefan–Boltzmann law gives $L=\pi ac[2G\mathcal{R}/(3{\kappa }_0){]}^{1/2}{K}_0^{-1}{\mu }^{2}M{R}^{5/2}$. The stellar virial theorem gives total energy $-3G{M}^2/(7R)$, so $L=-3G{M}^2\dot{R}/(7R^2)$. Integrating proves $C=(49\pi ac/6)[2\mathcal{R}/(3{\kappa }_{0}G){]}^{1/2}{K}_0^{-1}$. This holds before stellar nuclear fusion, changing opacity or electron degeneracy pressure invalidates the assumed contraction law.

For a monatomic fully convective ideal gas star, uniform specific entropy gives $P=K{T}^{5/2}$ with $K\propto {\mu }^{-5/2}{M}^{-1/2}{R}^{-3/2}$. Combining this with the stellar surface boundary condition $\kappa P=2g/3$ and $\kappa ={\kappa }_0{\rho }^a{T}^b$ givesAll composition comparisons here hold ${\kappa }_0$ fixed. For $a=1,b=4$, $T_e^8\propto {\mu }^4{M}^{2}R$, and the Stefan–Boltzmann law gives $L\propto {\mu }^{2}M{R}^{5/2}$. The stellar virial theorem and Kelvin-Helmholtz contraction then give $R^{-7/2}-R_0^{-7/2}\propto {\mu }^2{M}^{-1}(t-t_0)$, with a positive coefficient. Once the initial-radius term is negligible, $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$.

A Keplerian accretion disk with a zero-torque inner boundary condition and local blackbody cooling has $T_{\mathrm{eff}}^4=3GM\dot{M}(1-\sqrt{r_{\mathrm{in}}/r})/(8\pi \sigma r^3)$. The Stefan–Boltzmann law uses the flux through one face; the standard thin-disk dissipation flux includes both faces. The scale $T_{\mathrm{in}}=[3GM\dot{M}/(8\pi \sigma r_{\mathrm{in}}^3){]}^{1/4}$ is a normalization, not the actual temperature at the stress-free inner edge.

Assume a passive disk has local isothermal sound speed squared $c_s^2=aT$, disk scale height $H=c_s/{\Omega }_K$, and irradiated height $z_s=fH$ with constant $f$. Its irradiation grazing angle of a disk and Stefan–Boltzmann law give $q^{\prime }=C{q}^8/R^3$ for $q=z_s/R$, where $C=4\pi \sigma (GM/a{)}^4/(L_{\ast }{f}^8)$. Integrating gives $q^{-7}=R_0^2/R^2+B$, with $R_0^2=7C/2$. Where the first term dominates and the thin disk assumptions hold, $q\simeq (R/R_0{)}^{2/7}$, $H\propto R^{9/7}$ and $T\propto R^{-3/7}$. The law cannot be extrapolated into a geometrically thick region using the same grazing-angle approximation.

With negligible heavy-element contribution to the mean molecular weight, $\mu =4/(3+5X)$ and $\kappa \propto 1+X$. At fixed mass, composition-dependent CNO homology gives $L\propto (1+X{)}^{-1}(3+5X{)}^{-4}$ and $R\propto [X(1+X){]}^{1/16}(3+5X{)}^{-9/16}$. The Stefan–Boltzmann law then gives $T_e\propto X^{-1/32}(1+X{)}^{-9/32}(3+5X{)}^{-23/32}$. From the formal $X=1$ reference, declining $X$ initially increases luminosity and effective temperature, unlike the usual redward main-sequence turn-off. A trace CNO catalyst abundance is implicit in the burning coefficient; an exactly metal-free pure-hydrogen star cannot burn through that cycle.

Match a fully convective stellar polytrope of index $3/2$ to the grey-atmosphere convection onset with eighth-power temperature opacity. At the fixed matching optical depth, $T_b\propto T_e$ and $P_b\propto g^{2/3}{T}_e^{-16/3}$. The ideal-gas adiabat $P=K{\rho }^{5/3}$ gives $K\propto T_b^{5/3}{P}_b^{-2/3}\propto T_e^{47/9}{g}^{-4/9}$. The polytropic mass-radius relation gives $K\propto R{M}^{1/3}$; with $g\propto M/R^2$, this yields $T_e^{47}\propto R{M}^7$. The Stefan–Boltzmann law then gives the displayed luminosity. Composition and opacity normalization are held fixed between models.

Match an efficient, monatomic gas, fully convective star to the grey-atmosphere convection onset with thirteenth-power opacity. Its convective relation is $P=K_T{T}^{5/2}$. At the matching temperature $T_b\propto T_e$, the stellar atmosphere gives $K_T^2\propto M{L}^{-1}{T}_e^{-13}\propto M{R}^{13/2}{L}^{-17/4}$ using the Stefan–Boltzmann law. The ideal-gas relation converts this to $P=K_{\rho }{\rho }^{5/3}$ with $K_{\rho }=(\mathcal{R}/\mu {)}^{5/3}{K}_T^{-2/3}$. The polytropic mass-radius relation at index $3/2$ gives $K_{\rho }\propto G{M}^{1/3}R$, hence $K_T^2\propto M^{-1}{R}^{-3}$. Equating the two expressions proves $L^{17/4}\propto M^2{R}^{19/2}$, or $L\propto M^{8/17}{R}^{38/17}$. The entropy coefficient is constant throughout the convective envelope, but can change between models.

For a stellar homology sequence with $R\propto M^X$ and $L\propto M^Y$, the Stefan–Boltzmann law gives $T_{\mathrm{eff}}\propto M^{(Y-2X)/4}$. Eliminating mass gives the displayed Hertzsprung-Russell diagram slope when $Y\ne 2X$. The usual horizontal axis decreases in temperature toward the right, reversing the visual left-to-right sign. The relation is between logarithmic variables and depends on the specified structural assumptions.

Let $\mathcal{R}=k_B/({\mu }_m{m}_p)$ for constant mean molecular weight ${\mu }_m$ in proton-mass units, and let $\sigma$ be the Stefan–Boltzmann constant in the Stefan–Boltzmann law. The gas sound speed used here satisfies $c_s^2=\mathcal{R}T$, and $\Omega$ is the local orbital frequency. For the local opacity approximation $\kappa ={\kappa }_0{\rho }^{1/3}{T}^{10}$, a gas-supported optically thick alpha disk has $\rho \sim \Sigma /H$ and $\mathcal{R}T\sim {\Omega }^2{H}^2$. Integrated viscous heating is $F_{+}\sim \alpha \Sigma {\Omega }^3{H}^2$. Radiative diffusion gives $F_{-}\sim (\sigma /{\kappa }_0)T^{-6}{\Sigma }^{-4/3}{H}^{1/3}\sim (\sigma /{\kappa }_0){\mathcal{R}}^6{\Omega }^{-12}{\Sigma }^{-4/3}{H}^{-35/3}$. Equating the fluxes proves the thickness scaling. At fixed composition and $\alpha$, it gives $H/r\propto r^{53/82}{\Sigma }^{-7/41}$ and $c_s\propto r^{6/41}{\Sigma }^{-7/41}$ in Keplerian rotation. Radial powers therefore require a surface-density profile; they are not determined by opacity alone. The formula is local to the stated opacity regime and small disk aspect ratio.

Take homologous profiles and hold composition, hence mean molecular weight, fixed. The stellar mass conservation equation and stellar hydrostatic equation give ${\rho }_c\propto M/R^3$ and $P_c\propto G{M}^2/R^4$. With ideal gas pressure support, $T_c\propto M/R$. Integrating the specific nuclear rate through homologous profiles then givesThe radiative diffusion in a star equation gives $L_{\mathrm{rad}}\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Here ${\kappa }_c\propto {\rho }_c^{1/2}{T}_c^{-5/2}\propto M^{-2}R$, soEquality of generated and transported luminosity implies $M^7{R}^{-8}\propto M^5{R}^{-1}$, whenceThis is the radiative homology with density-half inverse-five-halves opacity scaling.

The Stefan–Boltzmann law defines the effective temperature by $L=4\pi R^2{\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$. Therefore $T_{\mathrm{eff}}\propto M^{29/28}$. In theoretical Hertzsprung-Russell diagram coordinates, the slope isThe usual graphical convention puts increasing temperature to the left, so the line descends as one moves right. If the horizontal coordinate is $-\log T_{\mathrm{eff}}$, its plotted slope is $-132/29$.

A fully convective, monatomic ideal gas is approximately isentropic through its deep interior. Its adiabatic index is $5/3$, so $P=K{\rho }^{5/3}$ and $T\propto K{\rho }^{2/3}$. The $n=3/2$ stellar polytrope scaling, or direct substitution of the central homology scalings into $K=P_c/{\rho }_c^{5/3}$, gives $K\propto R{M}^{1/3}$.

Near the photosphere, hydrostatic balance expressed in optical depth gives $P_{\mathrm{ph}}\sim g/{\kappa }_{\mathrm{ph}}$, with $g\propto M/R^2$. Treat $T_{\mathrm{ph}}$ as proportional to $T_{\mathrm{eff}}$ and match to the convective adiabat, neglecting a varying surface entropy jump in this homology approximation. Since $\rho \propto P/T$ and ${\kappa }_{\mathrm{ph}}\propto \rho T^{10}\propto P{T}^9$,The adiabat gives $P\propto K^{-3/2}{T}^{5/2}$, so the surface condition yieldsThe Stefan–Boltzmann law consequently gives $L_{\mathrm{surf}}\propto M^{4/7}{R}^{16/7}$. The same specific nuclear law as before still gives $L_{\mathrm{nuc}}\propto M^7{R}^{-8}$. Equating these scalings produces $R^{72/7}\propto M^{45/7}$. Thus the fully convective homology with tenth-power surface opacity isThen $T_{\mathrm{eff}}\propto M^{3/16}$, and the theoretical HR slope isIts steepness expresses the relatively weak change of effective temperature with mass. The graph using $-\log T_{\mathrm{eff}}$ has the negative of this slope.

In stellar homology, write $r=Rx$, $m=Mf(x)$ and $\rho ={\rho }_{c}h(x)$, with common dimensionless profiles. Integrating hydrostatic equilibrium from the surface, where the pressure is negligible, givesThe profile integral is common to the sequence. Similarly, mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^1{x}^{2}h(x)dx$. These prove the homology scaling of central stellar pressure:The ideal gas law consequently gives $T_c\propto \mu M/R$, with the gravitational and gas constants absorbed into the coefficient.

The radiative diffusion in a star equation has a characteristic gradient $T_c/R$. Hence $L\propto R{T}_c^4/({\kappa }_c{\rho }_c)$. Using Kramers' opacity law, and comparing fixed-composition homologous models, this becomesThe volume integral of the stellar energy-generation rate $ϵ={ϵ}_0\rho T^{\nu }$ instead givesThe dimensionless integrals behind these relations do not depend on $M$ along a homologous sequence. In stellar thermal equilibrium the two luminosities agree, proving Kramers homology with a variable nuclear exponent:For the specified proton–proton chain exponent $\nu =7/2$, the mass factor cancels. ThusThe model captures two important features of the lower main sequence: hydrogen burning primarily uses the proton–proton chain, and gas rather than radiation pressure supplies the principal support. A radiative interior and a Kramers' opacity law approximation are useful near the Sun. It is a simplified local model: the Sun has a convective envelope, the coolest main sequence stars become fully convective, and real opacities and burning exponents vary. In particular a literally constant radius is not a realistic law over the entire lower main sequence.

On this proton–proton chain branch $T_c\propto M$, because $R$ and mean molecular weight are fixed. Normalizing to the Sun gives the upper mass estimateFor the specified CNO cycle exponent $\nu =23/2$, the homology relation instead reads $R^{14}\propto M^8$. Substituting into the radiative luminosity relation givesThis branch retains the deliberately idealized fully radiative assumptions, even though real CNO cycle cores are often convective.

For a Hertzsprung-Russell diagram, use the Stefan–Boltzmann law, $L=4\pi R^2\sigma T_{\mathrm{eff}}^4$. On the proton–proton chain branch $T_{\mathrm{eff}}\propto M^{11/8}$, so $L\propto T_{\mathrm{eff}}^4$. On the CNO cycle branch $T_{\mathrm{eff}}\propto M^{57/56}$, giving $L\propto T_{\mathrm{eff}}^{292/57}$. Thus the upper branch rises slightly more steeply toward the hot, luminous upper left in a logarithmic Hertzsprung-Russell diagram. The sketch below matches the approximate branches continuously at their respective burning crossovers; the true crossover has both heating terms and is smooth.

The proton–proton chain begins with reactions between hydrogen nuclei, so its specific heating scales as $X_H^2\rho T^{7/2}$ in this approximation. The catalytic CNO cycle involves hydrogen and CNO nuclei, and scales as $X_H{Z}_{\mathrm{CNO}}\rho T^{23/2}$. Here $X_H$ is the hydrogen mass fraction and $Z_{\mathrm{CNO}}$ the CNO mass fraction, not the total abundance of every metal. ConsequentlyKeeping $X_H$ fixed while dividing $Z_{\mathrm{CNO}}$ by 256 gives the CNO crossover temperatureThe mass estimate uses the unchanged proton–proton chain branch below crossover, together with the stipulated unchanged opacity and mean molecular weight.

The low-CNO proton–proton chain sequence therefore lies on the same locus, but extends to twice the transition mass. Beyond its later crossover the low-CNO CNO cycle sequence has the same logarithmic slope $292/57$ but a different normalization. At a fixed mass where both sequences are CNO-dominated, lowering the nuclear coefficient by 256 givesThus those low-CNO models are smaller, hotter and modestly more luminous at the same mass. At the same effective temperature, their CNO branch is lower in luminosity than the solar-composition CNO branch; these comparisons hold different quantities fixed.

Efficient convection keeps the interior at nearly uniform specific entropy. For the fully ionized monatomic gas, the heat capacity ratio is $5/3$, so the adiabatic process gives $P=K_{\rho }{\rho }^{5/3}$, with $K_{\rho }$ spatially constant. The ideal gas relation $P=\mathcal{R}\rho T/\mu$ then impliesThe constant is spatial, not temporal: loss of specific entropy during contraction changes it. This is a stellar polytrope with polytropic index $n=3/2$. We neglect radiation pressure and take the thin photosphere as an atmospheric matching layer rather than a significant fraction of the stellar mass.

For a fixed dimensionless profile of the Lane-Emden equation, the stellar homology relations have the formwhere $C_{\rho },C_P$ are positive dimensionless structure constants. More explicitly, write $\rho (r)={\rho }_{c}f(x)$ and $m(r)=Mh(x)$ with $x=r/R$ and fixed dimensionless profiles. Mass conservation gives $M=4\pi {\rho }_c{R}^3{\int }_0^{1}f(x)x^{2}dx$. Integrating the hydrostatic pressure support equation, with negligible surface pressure on the interior scale, gives $P_c=(GM{\rho }_c/R){\int }_0^{1}h(x)f(x)x^{-2}dx$. Both integrals are fixed positive numbers, proving the two central scalings. Substituting these relations into $K=P_c/T_c^{5/2}$ givesThus the same polytropic index fixes $K_0$ throughout this homologous sequence, while the dimensional stellar polytrope constant changes.

At the photosphere, whose optical depth is $2/3$, identify the temperature with the effective temperature $T_e$. The ideal gas relation and the interior adiabat give ${\rho }_{\mathrm{ph}}=\mu K{T}_e^{3/2}/\mathcal{R}$. With the prescribed opacity and stellar surface boundary condition,Eliminating $K$ therefore provesUsing the Stefan–Boltzmann law, with $\sigma =ac/4$, gives the luminosity

The gravitational energy of a stellar polytrope is $\Omega =-6G{M}^2/(7R)$ for $n=3/2$. The stellar virial theorem for a monatomic gas gives $2U+\Omega =0$, and hence the total energy is $E=U+\Omega =-3G{M}^2/(7R)$. With no appreciable stellar nuclear fusion or accretion, energy conservation givesCombining the two expressions for $L$ yieldsMultiplication by $-7R^{-9/2}/2$ and integration proves the contraction law with density-linear fourth-power stellar opacity:When the accumulated contraction term dominates the initial-radius term, $R\propto {\mu }^{-4/7}{M}^{2/7}(t-t_0{)}^{-2/7}$. Substitution into the central ideal gas relation giveswhich has the stated $t^{2/7}$ behavior when the time origin is negligible. This late-time approximation still requires the assumed fully convective star, fixed composition and opacity law to remain applicable.

Use Planck units for the following formulas. The black-hole area theorem is a classical statement with hypotheses. Assume the Einstein field equations, the null energy condition and appropriate global future regularity/predictability. For a clean version of the focusing sketch, take the future horizon generators to remain on a regular horizon and to be future complete in affine parameter. The usual strong asymptotic predictability formulation supplies the global causal control needed to exclude the same focusing pathology; an arbitrary spacetime with naked future breakdown is not covered by the theorem.

Let ${\ell }^a$ be an affinely parametrized generator of a smooth portion of the event horizon, and let $\theta =d\log dA/d\lambda$ be its null expansion. The horizon is a null hypersurface, so its generators have zero null twist. In four dimensions the Null Raychaudhuri equation isThe Einstein field equations give $R_{ab}{\ell }^a{\ell }^b=8\pi T_{ab}{\ell }^a{\ell }^b\ge 0$, since the metric terms vanish when contracted with a null vector. The screen-space null shear term is also nonnegative, so ${\theta }^{\prime }\le -{\theta }^2/2$. If ${\theta }_0<0$ at ${\lambda }_0$, integrating this inequality while the congruence remains regular givesIt focuses within affine distance at most $2/|{\theta }_0|$. The vanishing area element is a focal point to an earlier spacelike horizon cut. Past such a point the null geodesic can no longer generate an achronal boundary, but the event horizon is precisely such a boundary. Future completeness and the global hypotheses exclude its escaping this contradiction by ending prematurely at a pathology. Therefore $\theta \ge 0$ everywhere the smooth horizon description applies.

Integrate $d(dA)/d\lambda =\theta dA$ between ordered horizon cuts. Every continuing generator contributes a nondecreasing area element. New generators can enter at past endpoints or merger crease sets, adding area; generators cannot leave a regular future horizon through future endpoints under these hypotheses. Thus, including the total area of disconnected components,This is the future-complete horizon focusing proof of the area theorem. A smooth elementary congruence proof requires the stated regularity; the general theorem also treats the horizon's nonsmooth joining set rather than assuming mergers have a globally smooth horizon.

Hawking radiation changes the situation because the quantum field state in a collapsing geometry is not a classical positive-energy fluid. For a large isolated Schwarzschild black hole,in natural units. Outgoing quantum modes at future null infinity have a nearly thermal occupation, modified by greybody factors; an accompanying negative Killing energy flux near the horizon reduces the hole's mass. On timescales short compared with the evaporation time but long compared with $M$, its geometry is approximately Schwarzschild spacetime with slowly decreasing $M(t)$. It is not an exactly stationary vacuum solution with a parameter changed by hand: the radiation's stress-energy tensor sources the time dependence through the semiclassical Einstein equation. The distant luminosity is carried by the outgoing radiation, while the horizon-area decrease accompanies the inward negative-energy contribution.

As the black hole loses mass, its horizon scale and Bekenstein-Hawking entropy shrink, its temperature rises and its luminosity increases. This is the negative heat capacity of a Schwarzschild black hole. Large astrophysical holes evaporate extremely slowly in isolation; small holes evaporate more rapidly, and increasingly massive particle species become accessible as the temperature rises. An incoming radiation bath or accretion can instead offset the loss, so the pure evaporation law assumes negligible incoming energy. A formal complete evaporation also raises the black hole information paradox; the low-curvature calculation alone does not determine how information or the final state is resolved.

For an elementary luminosity estimate, treat the horizon as emitting photon blackbody radiation. In natural units the Stefan–Boltzmann law has $\sigma ={\pi }^2/60$, givingThe horizon area used here is an estimate of emitting area, not the exact frequency-dependent absorption cross-section. More generally write $P=\alpha /M^2$, where $\alpha >0$ includes the greybody factors and the active particle species. For constant effective $\alpha$, energy balance and direct integration giveHere $t$ is time measured at infinity. The formal evaporation timescale is $M_0^3/(3\alpha )$; with the photon blackbody estimate it is $5120\pi M_0^3$. Restoring constants givesfor that specified estimate. If particle thresholds make $\alpha$ mass dependent, the precise leading law is instead $t={\int }_{M(t)}^{M_0}{m}^{2}dm/\alpha (m)$. The robust result is the inverse-square mass-loss rate for a fixed species regime and the cubic mass scaling of its lifetime. This is the semiclassical cubic mass law for Schwarzschild evaporation. Extrapolating its zero to a definite end state is unjustified once $M$ approaches the Planck mass and curvature requires quantum gravity.

There is no contradiction with the classical area theorem: the renormalized quantum stress-energy tensor can violate the null energy condition used in the focusing inequality. In the evaporation regime,which is permitted when that hypothesis fails. The thermodynamic replacement is the generalized second law, concerning $S_{\mathrm{outside}}+A/(4G\hslash )$ rather than horizon area alone. Outside radiation can carry entropy while the hole's Bekenstein-Hawking entropy decreases.

Buoyancy and the convective criterion. Consider a small fluid parcel displaced outward by $\delta r$. Assume that sound waves keep its pressure equal to the ambient pressure, while heat exchange is slow enough that its displacement is an adiabatic process. Composition is uniform, so the parcel and its surroundings have the same mean molecular weight. For an ideal gas, $P\propto \rho T$, and an adiabatic parcel obeys $P\propto {\rho }^{\gamma }$. If $\delta \ln P<0$ is the pressure change on moving outward, thenThe density excess of the parcel at its new location is thereforeIf $\nabla >(\gamma -1)/\gamma$, this excess is negative: the outward-displaced parcel is lighter and accelerates farther outward. A downward displacement similarly continues downward. Conversely, a smaller gradient gives a restoring buoyancy force. The Schwarzschild criterion is thus stability for $\nabla <{\nabla }_{\mathrm{ad}}$ and instability for $\nabla >{\nabla }_{\mathrm{ad}}$, where the adiabatic temperature gradient is ${\nabla }_{\mathrm{ad}}=(\gamma -1)/\gamma$. For a monatomic gas,A composition gradient would require the additional density effect in the Ledoux criterion; it is absent in this calculation.

Opacity and pressure. The supplied opacity is characteristic of a cool hydrogen-rich stellar atmosphere dominated by negative hydrogen ion opacity. Mostly neutral hydrogen can attach electrons supplied by partially ionized metals; bound-free and free-free absorption involving the resulting $H^{-}$ ions produces a steep positive temperature dependence. The weak density dependence and strong temperature sensitivity can be represented by local power-law fits. At fixed composition the given pressure form is equivalent to $\kappa \propto {\rho }^{1/2}{T}^{17/2}$, because $P\propto \rho T$. This is an approximate fit over the relevant cool-atmosphere regime, not an exact opacity law at all temperatures; molecular opacity or hydrogen ionization eventually changes the regime.

For a geometrically thin atmosphere take the surface gravity of a star to be constant, $g_s=GM/R^2$. Combining $dP/dr=-\rho g_s$ with $d\tau /dr=-\kappa \rho$ gives hydrostatic equilibrium in optical depth,The grey atmosphere relation gives $T^8=T_e^8(3\tau +2{)}^2/16$. With negligible external pressure, $P(0)=0$, integration yieldsAn imposed nonzero external pressure would add $P(0{)}^{3/2}$; the boundary assumption is what fixes the integration constant here.

The factor in parentheses is $3\tau /[2(3\tau +2)]$. HenceThe grey-atmosphere convection onset with eighth-power temperature opacity occurs when this reaches $2/5$, soThe radiative continuation below this boundary would be unstable. Efficient convection instead keeps the deeper star approximately adiabatic.

Matching to the convective interior. Write $P=K{\rho }^{5/3}$. Eliminating density using the ideal gas law givesat fixed composition. Since ${\tau }_c$ is a fixed number, the matching temperature $T_b$ is proportional to $T_e$, while the matching pressure satisfies $P_b^{3/2}\propto g_s{T}_e^{-8}$. ConsequentlyThe polytropic mass-radius relation for index $3/2$ gives $K\propto R{M}^{1/3}$. Raising the matching relation to the ninth power and using $g_s\propto M/R^2$ therefore givesFinally the Stefan–Boltzmann law $L=4\pi R^2{\sigma }_{\mathrm{SB}}{T}_e^4$ yields the Hayashi relation with eighth-power temperature opacity,Composition and the opacity normalization are fixed in this comparison. The entropy coefficient $K$ is spatially constant in each convective model, but need not be the same for stars with different masses or radii.

Mechanical balance and equipotentials. Let $\Phi$ be the Newtonian gravitational potential and $\varpi$ the cylindrical distance from the rotation axis. For rigid angular velocity $\Omega$, a fluid element at rest in the rotating frame has no Coriolis acceleration. The centrifugal acceleration is ${\Omega }^2\varpi \hat{\varpi }$, so hydrostatic equilibrium becomesUsing the Poisson equation for gravity and ${\nabla }^2(x^2+y^2)=4$ gives

On any regular connected equipotential surface, every tangent displacement satisfies $d\varphi =0$ and therefore $dP=0$. Locally $P=P(\varphi )$, with $dP/d\varphi =-\rho$. This makes $\rho$ a function of $\varphi$ as well. Equivalently, taking the curl of the hydrostatic equation gives $\nabla \rho \times \nabla \varphi =0$. Thus $P,\rho$ and ${\nabla }^2\varphi$ are constant on each surface. The last conclusion follows from the modified Poisson equation; it does not constrain the magnitude of the first derivative $|\nabla \varphi |$ to be constant there.

For uniform composition, an equation of state $P=P(\rho ,T)$ determines $T$ uniquely at fixed $P,\rho$ in the usual stellar regime. For the gas-plus-radiation equation, the derivative with respect to $T$ is positive. Hence $T=T(\varphi )$, and the opacity and radiative conductivity are also functions of $\varphi$.

Surface flux and the requested sketch. Define $f(\varphi )=-\chi dT/d\varphi$, positive for the outward radiative flux. ThenOn an outer equipotential surface, $f$ is constant and the outward normal is parallel to $\nabla \varphi$. Its effective gravity is $g=|\nabla \varphi |$, so $F_n=fg$. The Stefan–Boltzmann law $F_n={\sigma }_{\mathrm{SB}}{T}_e^4$ gives the radiative gravity-darkening law,A rapidly rotating star is oblate. Centrifugal support reduces the effective gravity near the equator, which is therefore cooler and dimmer; the poles are hotter. The original sketch below uses a subcritical Roche model of a uniformly rotating star to show this geometry and the corresponding temperature variation. The material temperature on an interior equipotential is not the same quantity as the emergent effective temperature defined by surface flux.

Why purely local radiative balance fails. Differentiating the flux givesOn one equipotential, $f,f^{\prime },{\nabla }^2\varphi$ and the proposed source $\rho ϵ(\rho ,T)$ are all constant, but $|\nabla \varphi {|}^2$ generally varies with latitude. Unless $f^{\prime }=0$ or there is another special cancellation, no single value of $\rho ϵ$ can match this flux divergence at all points on the surface. This is the radiative-equilibrium obstruction in a rotating barotropic star. Slow stellar meridional circulation can transport heat to compensate; it is distinct from abandoning hydrostatic balance at leading order.

Integrated energy balance with circulation. The steady continuity equation is $\nabla \cdot (\rho \mathbf{v})=0$. Combining the entropy–enthalpy relation with hydrostatic balance givesIntegrate the energy equation over $V$ and apply the divergence theorem:Uniform composition makes $h=h(\rho ,T)$ constant on the enclosing equipotential, and $\varphi$ is constant there by definition. Pull $(h+\varphi )$ out of the last integral; the remaining net mass flux is zero by steady continuity. Thus the equipotential luminosity conservation with stellar circulation isThis does not require $\mathbf{v}\cdot \hat{n}=0$ pointwise on $S$: inward and outward circulation may cancel in the net mass flux.

Finally, $f$ is constant on $S$ and the divergence theorem applied to the effective potential givesTherefore $-\chi dT/d\varphi =L/(4\pi Gm-2{\Omega }^{2}V)$. Since $dP/d\varphi =-\rho$,The mass $m$ is the mass enclosed by this surface, not necessarily the total stellar mass. The denominator is positive for the regular hydrostatic surfaces used here; near loss of centrifugal confinement the assumed stellar model must be reconsidered. The derivation uses slow circulation, retaining hydrostatic/equipotential thermodynamics to leading order while allowing a nonzero advective heat flux.