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.

A cryogenic radiometer uses low-temperature operation to improve electrical substitution radiometry through particular material and thermal properties. Lower heat capacity of a copper receiver permits a substantial cavity without proportionately slow relaxation; cold shields reduce unwanted thermal radiation, and heater leads with superconductivity can suppress resistive lead heating. These features do not make radiometric absorptance exactly one or abolish optical-electrical heating non-equivalence.^[ref]

Cold surroundings matter as much as a cold receiver. In a simplified gray-surface model, net outward thermal radiation is $ϵ\sigma A(T^4-T_s^4)$, using emissivity $ϵ$, area $A$, shield temperature $T_s$ and the Stefan-Boltzmann constant $\sigma$. Cold shields reduce the incoming term as well as the receiver's emission; cooling the receiver alone while leaving it exposed to warm surroundings does not remove incoming radiation. An entrance window still provides an optical and thermal path. Vacuum reduces gas-mediated heat transfer, while an intentional radiometer thermal link remains necessary. This model explains specific losses, rather than promising a universal improvement factor or perfect accuracy.

Choose a positive constant $\Sigma$ and the illustrative interface lawHere $T_h$ and $T_c$ are the hot and cold interface-side temperatures, not distant thermometer readings. This is a low-temperature phonon model under restrictive assumptions, not a universal fit to Kapitza's measured exponent or a material-property table. Its coefficient $\Sigma$ is not the Stefan-Boltzmann constant: the carriers in this model are phonons, rather than photons crossing a vacuum.

The fourth power has a simple counting explanation. In the three-dimensional low-temperature Debye model, a linear phonon dispersion relation gives a density of states proportional to ${\omega }^2$. A phonon has energy $\hslash \omega$ and Bose-Einstein distribution occupation $[\exp (\hslash \omega /k_{B}T)-1{]}^{-1}$. If propagation speeds and the angle-averaged transmission introduce no additional frequency scale over the occupied modes, the outgoing energy flux has the formExtending the upper integration limit assumes occupied modes lie well below the cutoff. Balancing the two outgoing fluxes, with coefficients satisfying zero net heat flux at equal temperatures, gives the chosen difference of fourth powers. Frequency-dependent transmission, different excitations or additional transport channels can change this conclusion.

The finite-drop area-specific thermal resistance follows by factoring a difference of powers:Near thermal equilibrium, write $T_h=T+\delta T$ and $T_c=T$. Expansion gives $q=4\Sigma T^3\delta T+O((\delta T{)}^2)$, so the differential area-specific Kapitza resistance isAt a finite bias, varying only the hot side gives $(\partial q/\partial T_h{)}^{-1}=1/(4\Sigma T_h^3)$; this is generally different from $r_{K,\sec }$. Stating what is held fixed makes a differential resistance reproducible.

For an original normalized example, choose reference scales $T_0>0$ and $A_0>0$, and setFor contact area $A=2A_0$, the total heat transfer rate is $\dot{Q}=2.1472A_0\Sigma T_0^4$. The finite-drop resistance satisfies $\Sigma T_0^3{r}_{K,\sec }=0.2/1.0736=125/671\approx 0.18629$. Using the equilibrium tangent at $T_0$ instead would predict $j=4(0.2)=0.8$, illustrating why a substantial temperature rise requires the nonlinear law. Doubling the common equilibrium temperature reduces $R_K$ by eight in this chosen model.

A series check can be done without guessing any material properties. Stipulate a hot-side slab with dimensionless areal bulk resistance $b_s=\Sigma T_0^3{L}_s/k_s=0.10$, and neglect the cold-side bulk drop. At the same $j=1.0736$, the hot reservoir must have $T_H/T_0=1.2+b_{s}j=1.30736$. The imposed reservoir drop splits into $0.10736$ in the bulk and $0.20$ at the interface. Applying the fourth-power law directly to $T_H$ would assign the bulk drop to the interface and overestimate the heat flux.

Finally, transfer heat at rate $Aq$ between maintained interface-side temperatures $T_h>T_c>0$. The hot side loses entropy at rate $Aq/T_h$, and the cold side gains it at rate $Aq/T_c$. With no interface entropy accumulation, the entropy production per unit area isFor the normalized example it is ${\dot{S}}_{\mathrm{prod}}/(A\Sigma T_0^3)=1.0736(1-1/1.2)=0.178933\dots$. The same energy crosses both sides in a steady state, while total entropy increases, consistently with the Second law of thermodynamics. A finite temperature jump during heat transfer therefore expresses dissipation, not a failure of energy conservation.

Ice-ocean albedo feedback is the ocean-surface case of ice-albedo feedback: loss of bright sea ice, often covered by snow, exposes darker water, lowers surface albedo and increases absorbed sunlight. The extra stored ocean heat content can promote further melting and delay refreezing. Reversing the sequence can amplify cooling. This is a positive climate feedback, meaning that it reinforces the initial change; the word positive does not mean beneficial, and it does not guarantee runaway warming.

For unchanged downward solar irradiance $S_{\downarrow }$ on a newly exposed opaque patch, the extra absorbed heat flux isFor illustrative local surface albedos $0.60$ and $0.10$, and $S_{\downarrow }=200\mathrm{W}{\mathrm{m}}^{-2}$, the patch absorbs an additional $100\mathrm{W}{\mathrm{m}}^{-2}$. This is a local sunlight calculation, not a whole-planet radiative forcing: newly exposed area, illumination, clouds and season determine its global contribution. With no incident sunlight this particular heating pathway vanishes. NASA's notes describe the contrast between reflective ice and snow and darker water. ^[ref]

To isolate how a positive climate feedback competes with Planck feedback, retain the assumptions of the zero-dimensional planetary energy balance and prescribe the smooth broadband Bond albedoThe exponential function makes a smooth transition from a bright cold state near $0.6$ to a darker warm state near $0.3$. The midpoint $240\mathrm{K}$ and width $3\mathrm{K}$ are chosen to illustrate the dynamics. They are not a seawater freezing law, an observed global relation, or a calibration of modern Earth. No evolving ice thickness or latent heat is represented.

An equilibrium satisfies $F(T_{\ast })=0$. For a small perturbation $\delta T$, linear stability follows fromSince $C>0$, $F^{\prime }<0$ gives an asymptotically stable equilibrium, while $F^{\prime }>0$ gives an unstable equilibrium. With $x=(T-240\mathrm{K})/(3\mathrm{K})$, the derivative isThus the optical term $-S{\alpha }^{\prime }/4$ is destabilizing, while the thermal radiation term $-4\sigma T^3$ is stabilizing. Their sum, not the optical term alone, decides the linear stability.

For the same $S=1360\mathrm{W}{\mathrm{m}}^{-2}$ and Stefan-Boltzmann constant used above, the bisection method on the intervals $(220,223)$, $(238,241)$ and $(253,255)$, with temperature in kelvins, gives:

Unrounded numerical zeros of a function have $|F(T_{\ast })|<10^{-10}\mathrm{W}{\mathrm{m}}^{-2}$; the table rounds them for readability. In the positive-temperature domain the sign of $F$ alternates between heating, cooling, heating and cooling across these three equilibria. The middle equilibrium separates the basins of attraction of the two stable equilibria. For example, a starting temperature of $235\mathrm{K}$ cools toward the cold state, whereas $245\mathrm{K}$ warms toward the warm state: this is bistability at identical external illumination.

The same model illustrates hysteresis. Its equilibrium curve can be written as $S(T)=4\sigma T^4/[1-\alpha (T)]$. The two turning points, where $4[1-\alpha (T)]+T{\alpha }^{\prime }(T)=0$, are saddle-node bifurcations: approximately $(T,S)=(232.637\mathrm{K},1567.801\mathrm{W}{\mathrm{m}}^{-2})$ and $(245.893\mathrm{K},1250.481\mathrm{W}{\mathrm{m}}^{-2})$. For illumination strictly between these thresholds there are two stable equilibria and one unstable equilibrium. Under sufficiently slow increase of $S$, the cold branch ends at the upper threshold; under sufficiently slow decrease, the warm branch ends at the lower threshold. The ensuing warming and cooling transitions differ, so the state depends on its history. Finite sweep rates can delay the transitions; these thresholds describe this model's static equilibrium curve, not predicted tipping points for modern Earth.

The striking feature is that a smooth albedo law can produce distinct long-time states and abrupt transitions, although no discontinuity was inserted into the optical law. It is useful for understanding why positive climate feedback can amplify a response while stabilizing thermal radiation still limits it. Here $0.3<\alpha (T)<0.6$ bounds the incoming heating, whereas outgoing thermal radiation grows as $T^4$; high enough temperatures therefore cool rather than run away indefinitely.

Historically, M. I. Budyko's 1969 paper coupled changing albedo from expanding ice cover to a latitude-dependent temperature and heat transfer calculation, showing how that coupling amplified the effect of changed sunlight. ^[ref] William D. Sellers's 1969 paper constructed an annual energy-budget model in 10-degree latitude belts, including planetary albedo, atmospheric infrared transmission and atmospheric and oceanic exchange. ^[ref] These primary papers provide historical examples of ice-sensitive energy-budget modeling; the smooth single-temperature example above is a separate illustration, not a reproduction of either paper's equations or predictions.

The vertical hydrostatic equilibrium equation balances the gradient of the total pressure against gravity. In a thin disk with Keplerian rotation, expanding the central gravitational field about the midplane gives vertical acceleration $-{\Omega }^{2}z$, so $p_z=-\rho {\Omega }^{2}z$. Disk self-gravity is neglected in this approximation.

The flux-gradient equation is local thermal balance: the upward radiative flux increases with height because differential rotation dissipates energy. For a Newtonian effective dynamic viscosity $\mu$, the volumetric heating rate is $\mu (r{\Omega }^{\prime }{)}^2$, which becomes $9\mu {\Omega }^2/4$ for Keplerian rotation. The symbol $\mu$ here is dynamic viscosity, not kinematic viscosity; the latter is $\nu =\mu /\rho$.

The radiative diffusion equation applies to an optically thick medium in local thermodynamic equilibrium. Its diffusion coefficient is $c/(3\kappa \rho )$, and the equilibrium radiation energy density is $u_{\mathrm{rad}}=4\sigma T^4/c$. HenceThe minus sign means energy flows down the temperature gradient. Here $\sigma$ is the Stefan-Boltzmann constant, $\kappa$ the opacity per unit mass, and $c$ the speed of light. Constant Thomson scattering opacity is the adopted transport approximation; local thermodynamic equilibrium is an additional assumption behind the radiation equation of state.

Finally the equation of state adds ideal gas gas pressure $p_g=\rho kT/({\mu }_m{m}_H)$ and isotropic radiation pressure $p_r=u_{\mathrm{rad}}/3=4\sigma T^4/(3c)$. The Boltzmann constant is $k$, and ${\mu }_m{m}_H$ is the mean mass per gas particle. Both contributions supply the total pressure entering vertical hydrostatic equilibrium.

The first equation is vertical hydrostatic equilibrium: the upward pressure gradient balances the central mass's vertical gravitational force. For $|z|\ll r$, the point-mass acceleration is $-GMz/(r^2+z^2{)}^{3/2}\simeq -{\Omega }_K^{2}z$, where Keplerian rotation gives ${\Omega }_K^2=GM/r^3$. The disc's self-gravity and vertical inertia are neglected.

The second equation is local thermal balance. The divergence of the vertical radiative flux equals viscous dissipation. In the convention of this alpha disk, $\nu =\alpha p/(\rho \Omega )$, and the dominant Keplerian shear has $rd\Omega /dr=-3\Omega /2$. ThereforeThe radiative flux increases from zero at the symmetric midplane to a positive outgoing value at the upper surface. Radial heat transport, thermal advection and rapid storage of thermal energy are excluded from this local balance.

The third equation is radiative diffusion, equivalentlyIt applies to optically thick material, with $\kappa$ the opacity per unit mass and $\sigma$ the Stefan-Boltzmann constant. Here constant electron-scattering opacity represents Thomson scattering; the formula uses a diffusive, nearly thermal radiation field. Positive upward flux requires temperature to decrease upwards.

The final equation is the ideal gas pressure law, with mean molecular weight ${\mu }_m$ in proton-mass units. Writing $\mathcal{R}=k/({\mu }_m{m}_p)$ gives $p=\mathcal{R}\rho T$. The pressure is gas pressure because radiation pressure is assumed negligible. The four equations express vertical support, local heating balance, diffusive cooling and the gas equation of state.

Use geometrized units with $\hslash =k_B=1$ initially, and write $f(r)=1-2M/r$, $h(r)=1-Q^2/(Mr)$. The given nonextremality assumption makes $h_H=h(2M)=1-Q^2/(2M^2)>0$. Thus the static Killing vector field has a simple zero of its norm atThe exterior $r>2M$ is static. For radial null geodesics, $dt/dr=\pm 1/f$, independently of $h$. With $v=t+\int dr/f$, the metric iswhich is regular at $r_H$ since $h_H>0$. The outgoing null equation is $dr/dv=f/2$. Inside the future horizon these rays move to smaller radius, while outside they can escape to infinity. Hence the zero is the event horizon. The inner zero of $h$ is not another regular Killing horizon: the Killing norm there diverges rather than vanishes. The strict charge bound matters when making the regular extension.

For the Euclidean black-hole regularity condition, continue $t=i{t}_E$, consistent with the sign of the given quantum time-translation operator. Set $x=r-2M$. To leading order $f=x/(2M)$ and $h=h_H$, so the radial–time part of the Euclidean metric isIntroduce the radial proper length coordinate $\rho$ by $x=h_H{\rho }^2/(8M)$. ThenThe first two terms are a polar plane with angle $t_E/(4M)$. Smoothness at its origin requiresA different period would make a conical singularity. The time coordinate is normalized to proper time at infinity, so no further asymptotic redshift factor is needed. The cancellation of $h_H$ proves the charge independence of the Euclidean horizon period under a regular radial conformal factor; this argument would fail at $h_H=0$.

To justify the thermal interpretation, analytic continuation of the time-evolution operator $U(t)=e^{iHt}$ gives $U(i{t}_E)=e^{-H{t}_E}$. Its matrix element is an imaginary-time propagation kernel. Identifying the initial and final configurations and integrating over them takes its trace, giving the Euclidean path integral $Z=\mathrm{Tr}{e}^{-\beta H}$. In an energy eigenbasis each state has weight $e^{-\beta E}$, exactly the Boltzmann factor of a canonical ensemble. Thus the geometrically required period is the inverse Hawking temperature. With dimensions restored, an imaginary time period $\Delta t_E$ gives $\beta =\Delta t_E/\hslash =1/(k_{B}T)$. ConsequentlyHere $\mathcal{M}$ is mass in ordinary units; for a neutral hole $M=G\mathcal{M}/c^2$ and the physical imaginary-time period is $8\pi G\mathcal{M}/c^3$. Bosonic thermal fields are periodic and fermionic fields antiperiodic over that same Euclidean period. The relation to the canonical density operator, not just dimensional reasoning, identifies the temperature.

For an evaporation estimate take a neutral, isolated Schwarzschild black hole and approximate its radiated power with the Stefan–Boltzmann law. Its area is $A=16\pi G^2{\mathcal{M}}^2/c^4$, and the photon Stefan-Boltzmann constant is ${\sigma }_{SB}={\pi }^2{k}_B^4/(60{\hslash }^3{c}^2)$. Substituting the temperature givesThe black-hole evaporation speeds up as the mass decreases. Integrating, rather than dividing the initial rest energy by the initial power, givesThe corresponding initial Hawking temperature is about $6\times 10^{-8}$K. This is the Stefan-Boltzmann estimate of Schwarzschild evaporation time at the level of accuracy requested. Greybody factors and other emitted species alter the numerical coefficient, but not the leading cubic scaling of a large isolated hole. The estimate assumes negligible accretion or ambient radiation absorption; a stellar-mass hole in the present cosmic radiation bath is not currently losing mass at that isolated-hole rate. Finally the semiclassical formula does not describe the Planck-scale endpoint: extrapolating it to zero mass is an estimate of the macroscopic lifetime, not a proof of what the final quantum object is.

Let $h_i,h_s$ denote ice and snow thickness, let $z$ increase downwards, and count external atmospheric fluxes as positive into the surface. The Maykut-Untersteiner model is a time-dependent column model, rather than a single thickness formula: heat conduction, heat storage, internal absorption of sunlight and moving freezing/melting boundaries must be solved together. In an ice-following description its energy equation can be writtenwhere $e_i$ is specific enthalpy, $S$ is the prescribed ice salinity, and $I$ is the downward penetrating shortwave flux. The material derivative includes coordinate motion if a moving numerical grid is used. A common attenuation form is $I=I_0{e}^{-\kappa z}$ within a homogeneous ice layer, with the origin taken at that layer's top. The corresponding source is $\kappa I$. Snow-covered cases have much less penetration; absorbed surface radiation must not also be counted as internal heating.

The liquid brine fraction contributes latent heat to the effective specific heat capacity: for the illustrative liquidus approximation ${\varphi }_b\simeq -\mu S/T$ with Celsius $T<0$ in the mixed-phase range $0<{\varphi }_b<1$, $c_{\mathrm{eff}}={\partial }_T{e}_i\simeq c_0+L_f\mu S/T^2$. The precise ice conductivity and heat capacity are salinity- and temperature-dependent. This is why the full transient temperature profile need not be linear. At the snow-ice interface impose continuity of temperature and conductive heat flux; at the ice bottom impose the seawater freezing temperature $T_f$.

Write $I_0=i_0(1-\alpha )F_{\mathrm{SW}}$ and define the atmospheric input absorbed at the surface byHere $\alpha$ is the surface albedo, $\sigma$ is the Stefan-Boltzmann constant, and $F_H,F_E$ are signed sensible and latent turbulent heat fluxes. With upward conductive flux $q_c=k{\partial }_{z}T$, the cold, nonmelting surface satisfies $Q_A+q_{c,s}=0$. Once the surface reaches its melting temperature, constrain it there and use excess energy for ablation:The mass-loss rate ${\dot{m}}_s$ first removes snow at thickness rate ${\dot{m}}_s/{\rho }_s$; after the snow is gone it removes ice at rate $m_i={\dot{m}}_s/{\rho }_i$. Snow thickness also changes through prescribed snowfall and any separately represented compaction or sublimation.

At the bottom, the Stefan condition equates latent heat released by growth to upward conductive removal minus upward ocean heat flux $F_o$:Here $g_b$ is positive for basal freezing and negative for basal melting, and $L_b$ is the appropriate effective latent heat. These equations, the incident fluxes and snowfall determine the seasonal thickness trajectory. Annual sea-ice thermodynamic equilibrium means a periodic temperature field and $h_i(t+\mathcal{T})=h_i(t)$, soWith constant $L_b$, this is equivalently $\int (q_{c,b}-F_o-{\rho }_i{L}_b{m}_i)dt=0$. An “equilibrium thickness” must refer to the mean or a specified phase of this seasonal cycle.

A useful winter reduction neglects heat storage and solar penetration and treats each conductivity as constant. The two layers act as thermal resistances in series:This displays the dependence on snow and ice thickness, but is not a replacement for the full seasonally forced Maykut-Untersteiner model.

Increasing ocean heat flux subtracts directly from basal growth and increases bottom ablation. On a stable perennial branch it lowers annual equilibrium thickness, and sufficiently large input can eliminate a year-round ice cover. With zero oceanic input, removing this heat sink can allow much thicker thermodynamic ice; it does not remove the seasonal surface energy balance.

Additional annual snowfall has competing effects. Larger $h_s/k_s$ suppresses winter heat conduction and basal growth, favoring thinner ice. But persistent bright snow raises albedo, consumes energy when it melts and delays exposure of dark bare ice, suppressing summer ice ablation. The original column calculations show near compensation for moderate snowfall and substantial thickening when much more snow survives into summer. Thus annual equilibrium thickness need not decrease with snowfall, even though greater snow thickness at fixed winter temperatures always reduces the conductive growth rate. The balance also depends on snowfall timing and the treatment of snow-to-ice conversion. These seasonal responses are illustrated by the original model calculations.

Reducing surface albedo increases absorbed shortwave radiation. It promotes surface ablation, internal heating and, where sunlight reaches the ocean, basal ablation; the ice-albedo feedback reinforces the reduction in equilibrium ice thickness. Increasing albedo has the opposite effect. This sensitivity is strongest when sunlight is available and dark surfaces are exposed, rather than during polar night.

For the observed bottom melt, use ${\rho }_i\simeq 917\mathrm{kg}{\mathrm{m}}^{-3}$, $L_f\simeq 3.34\times 10^5\mathrm{J}\mathrm{k}{\mathrm{g}}^{-1}$ and approximately 122 days. The basal melt heat flux equivalent isThe integrated latent heat is about $6.1\times 10^8\mathrm{J}{\mathrm{m}}^{-2}$. This is energy used for melting, not automatically the complete ocean-to-ice heat input: the basal balance gives $F_o={\rho }_i{L}_b(-g_b)+q_{c,b}$, with sensible heat and salinity corrections if needed.

A likely major source is solar heating of an unusually large open-water fraction. Dark water absorbs much more sunlight than snow-covered ice, and warm water transfers that stored heat to floe bottoms by mixing and lateral transport. Thin or ponded ice also transmits more sunlight. Wind-driven divergence and export expose additional water, while advection of warm Pacific or Atlantic water and mixing of subsurface heat can contribute regionally. The supplied melt amount alone cannot partition these sources. The 2007 Beaufort Sea heat-budget analysis found that locally absorbed solar heat was sufficient to account for the enhanced basal melt. The latent-heat equivalent is about $60\mathrm{W}{\mathrm{m}}^{-2}$; increased solar absorption by exposed ocean is a well-supported principal explanation.

The first relation is vertical hydrostatic equilibrium. In a thin disc, the central object's gravity near the midplane is $g_z=-{\Omega }^{2}z$; the pressure gradient balances it. The second is an equation of state combining ideal gas pressure with isotropic radiation pressure. Write $\mathcal{R}=k/({\mu }_m{m}_p)$, so $p_g=\mathcal{R}\rho T$ and $p_r=4\sigma T^4/(3c)$; $\sigma$ is the Stefan-Boltzmann constant and ${\mu }_m$ the mean molecular weight in proton-mass units.

The third is local energy balance: the increase of upward radiative flux with height equals viscous dissipation. Here $\mu$ is dynamic viscosity, not mean molecular weight. A Keplerian disc has shear $rd\Omega /dr=-3\Omega /2$, so its heating rate is $\mu (r{\Omega }^{\prime }{)}^2=9\mu {\Omega }^2/4$. The fourth is optically thick radiative diffusion, with opacity per unit mass $\kappa$ and diffusion directed down the temperature gradient. Midplane symmetry gives $F(0)=0$ and oppositely directed fluxes from the two faces. The four relations express vertical force balance, gas–radiation thermodynamics, local heating, and diffusive cooling.

The standard thin-disk dissipation flux, summed over both faces, follows by substituting the Keplerian accretion disk profile into the viscous heating rate:The Stefan–Boltzmann law gives $\mathcal{C}=2\sigma T^4$ because there are two emitting faces. Local radiative equilibrium consequently gives the effective-temperature profile of a zero-torque diskHere $\sigma$ is the Stefan-Boltzmann constant. The scale $T_{\mathrm{in}}$ is not the temperature exactly at the inner boundary: the zero-torque inner boundary condition makes that formal temperature zero. Maximizing $T^4$ shows that the maximum effective temperature of a zero-torque disk occurs at $r=49r_{\ast }/36$, with $T_{\max }\simeq 0.488T_{\mathrm{in}}$.

At equal central mass and accretion rate, characteristic effective temperatures scale as $r_{\ast }^{-3/4}$. Thus, comparing corresponding values of $r/r_{\ast }$,The neutron-star disk is about 180 times hotter. The Planck law and Wien displacement law move its characteristic emission to about 180 times higher frequency, or 180 times shorter wavelength. A white dwarf disk commonly emits in optical and ultraviolet bands, while the hotter neutron star disk can emit in X-rays. Absolute bands require an actual accretion rate; the relative shift follows directly from the stated scaling.

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.

Write $\kappa =R_d/c_p$, where $R_d$ is the specific gas constant and $c_p$ the specific heat capacity at constant pressure. Conservation of potential temperature on a dry adiabatic process givesLet $Q=(1-{\alpha }_c)S$ be the absorbed stellar irradiance, and let $T_0=(Q/\sigma {)}^{1/4}$, with $\sigma$ the Stefan-Boltzmann constant. The transparent gas contributes no thermal radiation. For a cloud with emissivity $e$, Kirchhoff's law of thermal radiation gives thermal absorptivity $e$; assume thermal transmissivity $1-e$ and negligible thermal reflection, as in the usual one-layer grey body model. The outgoing radiant flux is the transmitted surface emission plus upward cloud emission:ThusNo factor of four is needed: $S$ is the incident irradiance of the specified column, rather than a stellar-beam irradiance that still needs a spherical average.

For fixed ${\alpha }_c,p_c,p_s$, increasing $e$ monotonically warms the ground when $0<q<1$. At $e=0$, $T_g=T_0$; at $e=1$, $T_c=T_0$ and $T_g=T_0/q$. A cloud at the surface, $q=1$, produces no such warming. Comparing with a genuinely cloud-free column also requires including its different albedo: relative to $(S/\sigma {)}^{1/4}$, the ratio is $[(1-{\alpha }_c)/(1-e+e{q}^4){]}^{1/4}$. Reflection can outweigh the greenhouse effect.

The governing balance here is the whole-column radiation budget. Imposing an additional purely radiative surface balance would incorrectly remove the convective heat transfer that the prescribed dry adiabatic process may require.

Interpret habitability as the ability to sustain life with broadly terrestrial biochemistry. Accessible liquid water, usable free energy, suitable chemical ingredients and sufficiently persistent conditions are the central requirements. Earth-like geography, a breathable atmosphere and actual habitation are not part of that definition.

A viable local physical environment. The relevant temperature and pressure must permit liquid water where the organisms live. For a surface habitat this normally requires an appropriate atmosphere and surface temperature; the roughly 6.1 millibar triple-point pressure of pure water is a minimum phase constraint, not a guarantee of a temperate surface. High salinity, extreme acidity and low water activity can make apparently wet environments unsuitable for known life. A subsurface ocean beneath ice can provide liquid conditions even when the surface is frozen, so the entire body need not lie in the conventional surface-water habitable zone.

A suitable orbit and rotation. For a surface-water planet, orbital distance must give an absorbed radiation budget compatible with liquid conditions after the body's Bond albedo, atmosphere and greenhouse effect are included. The simple planetary equilibrium temperature estimateuses orbital irradiance $S$, Bond albedo $A$ and the Stefan-Boltzmann constant $\sigma$, assumes unit thermal emissivity and averages over the whole sphere. It is not a prediction of surface temperature beneath a atmosphere with a greenhouse effect. Excessive absorbed energy risks a runaway greenhouse; inadequate surface heating can freeze surface water, although heat from the interior or tidal heating can maintain a subsurface habitat.

The orbit must persist without collision or ejection for a biologically useful interval. Large orbital eccentricity can produce extreme variations in heating, but neither exactly circular motion nor Earth's particular obliquity is essential. Rotation period and obliquity govern illumination and seasons; atmospheric and ocean heat flux determine how severe those variations become. Tidal locking does not automatically make a body uninhabitable, because sufficient heat transport can moderate the day–night contrast. A natural satellite also needs a stable orbit around its primary and a tolerable combination of irradiation and tidal heating: enough heating may maintain an ocean, whereas excessive heating can make the environment hostile. These are constraints on the body's location and path, without requiring a particular kind of surrounding planetary system.

An adequate chemical inventory and access to it. Known life uses carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, together with catalytic trace elements. A rocky planet supplies mineral surfaces and reservoirs, but total inventory is not enough: chemical weathering, water-rock interaction or other transport must make nutrients accessible. A water-rich body with an ocean isolated from its rocky interior by a thick high-pressure ice layer could have difficulty recycling some nutrients; this is a possible limitation, not proof of sterility.

A usable energy supply. Photosynthesis can exploit light in an illuminated habitat. A dark subsurface habitat instead needs sustained chemical disequilibrium, for example oxidized and reduced species brought together by water-rock interaction, hydrothermal circulation or radiolysis. Warm water at complete chemical equilibrium is not, by itself, an energy supply for sustained metabolism.

Retention, protection and longevity. The body's mass, radius, surface gravity and thermal history affect retention of water and an atmosphere. Escape velocity must be considered together with gas temperature, molecular mass, atmospheric escape and replenishment: there is no universal Earth-mass threshold. An atmosphere, ice or rock overburden can protect life from damaging ionizing radiation and short-wavelength light. An intrinsic planetary magnetic field can modify charged-particle exposure and escape, but it is neither a universal requirement nor a guarantee of retention.

Finally, conditions must last long enough. Stored formation heat, radiogenic heating and possibly tidal heating can maintain circulation and chemical exchange. Plate tectonics can recycle materials and support a carbonate-silicate cycle, but it is not established as a necessary condition for every habitable body: conditional long-term carbon cycling is possible on a stagnant-lid planet as well.^[ref] The necessary end result is persistence of a suitable environment, rather than a checklist requiring every feature of the modern Earth.

Define the requirement before listing planetary properties. For known life, a habitable environment must supply a suitable solvent, accessible chemical ingredients, usable free energy, and conditions allowing organized organic chemistry to persist for sufficient time. Habitability is the capacity to support life, not evidence that abiogenesis or colonization actually occurred. Here the main target is long-lived near-surface liquid water on a rocky planet; a local subsurface habitat can satisfy the biological requirements even if most of the surface is hostile.

Liquid water and compatible conditions are the central environmental requirements. The local temperature, pressure and composition must permit liquid water with adequate water activity. The pure-water triple point, approximately $273.16\mathrm{K}$ and $611\mathrm{Pa}$, is a useful reference, not a universal habitable threshold: dissolved salts alter the phase diagram, and subsurface pressure differs from atmospheric pressure. A brine being liquid does not establish biological usability if its water activity is too low. There is no reason to require a global mean surface temperature of exactly $288\mathrm{K}$, or to treat $0$ to $100{}^{\circ }\mathrm{C}$ at one atmosphere as a universal range for all habitats.

Composition, nutrients and energy must be accessible together. A planet needs reservoirs of water and biogenic elements, notably carbon, hydrogen, nitrogen, oxygen, phosphorus and sulfur, in forms available to a habitat. Bulk abundance is not enough if nutrients are locked away from reacting fluids. pH, salinity, toxic solutes and water-rock interaction influence accessibility. Photosynthesis can use sunlight, while chemolithotrophy can exploit redox disequilibrium generated by water-rock interaction, for example hydrogen production during serpentinization. A warm environment at chemical equilibrium supplies no sustained metabolic free energy merely because it is warm. Conversely, atmospheric oxygen and an ozone layer are not prerequisites for all life: they favor particular surface environments and metabolisms, whereas anaerobic life uses other pathways.

An appropriate atmosphere and surface energy balance favor surface habitability. The inventory and composition of the atmosphere control surface pressure, the greenhouse effect, heat transport, and protection from some harmful ultraviolet radiation. Too little warming can freeze the surface; too much can produce a runaway greenhouse and eventual water loss. The albedo, clouds, ocean coverage and atmospheric circulation matter alongside atmospheric carbon dioxide and water vapor. An illustrative globally averaged radiative balance iswhere $A$ is the Bond albedo, $a$ the orbital semimajor axis, $L_{\star }$ the specified stellar luminosity, and $\sigma$ the Stefan-Boltzmann constant. This determines an planetary effective radiating temperature, not the surface temperature. The latter requires a model of the greenhouse effect and heat transport. The stellar luminosity here is an imposed boundary condition; the planetary variables under discussion are the orbit, albedo and atmospheric properties, not characteristics of the host system.

Mass and interior properties govern retention and replenishment. Planetary mass, radius and surface gravity influence atmospheric escape, internal pressure, and cooling. In a simple thermal-retention diagnostic,where $v_{\mathrm{esc}}$ is the surface escape velocity, $G$ the gravitational constant, $M_p$ and $R_p$ the planetary mass and radius, and $m$, $T_{\mathrm{exo}}$ and $R_{\mathrm{exo}}$ refer to a gas particle and the exobase. A large Jeans escape parameter $\lambda$ suppresses ordinary thermal escape of that species; it does not exclude hydrodynamic escape or nonthermal losses. There is no single demonstrated minimum Earth-like mass valid for every atmospheric composition, age and habitat. Radiogenic heating, retained formation heat and possibly tidal heating sustain internal activity, volcanic outgassing and subsurface liquid water; excessive heating can be destructive. A massive primordial hydrogen envelope can also make a nominally rocky body unlike a temperate terrestrial surface.

The orbit must supply a tolerable irradiation history. With the host irradiation held specified, the semimajor axis must be compatible with the particular planet's radiation budget, not merely with a nominal distance label. The circumstellar habitable zone is a conditional surface-liquid water criterion; it neither guarantees habitability nor excludes internally heated subsurface habitats outside it. Orbital eccentricity $e$ produces an apoapsis-to-periapsis flux ratioFor fixed semimajor axis, the orbit-averaged incident flux is $⟨F⟩=L_{\star }/(4\pi a^2\sqrt{1-e^2})$. Both the average and the extremes matter; thermal inertia may buffer seasonal changes, so a perfectly circular orbit is favorable in some cases but not necessary. Sustained orbital stability matters more than requiring $e=0$.

Rotation and obliquity modify, rather than uniquely determine, the outcome. Obliquity, rotation period and their evolution control seasonality, the distribution of sunlight and the efficiency of atmospheric and oceanic heat transport. Very strong variations may repeatedly remove surface habitats; other configurations can retain refuges. Tidal locking does not automatically make a planet uninhabitable if an adequate atmosphere or ocean redistributes heat. Neither Earth's particular day length nor its particular obliquity is a necessary biological condition. Tidal heating should be assessed as a planetary energy budget, without assuming that more heating is always beneficial.

Long-lived recycling is favorable; its mechanism need not be unique. Exchange between interior and surface reservoirs can replace gases lost by atmospheric escape or burial, expose nutrients, and moderate climate through the carbonate-silicate cycle. Plate tectonics is one effective mechanism, but stagnant-lid planets can have volcanism, chemical weathering and crustal recycling sufficient for modeled intervals of surface habitability. These are conditional model results, not proof that such planets are inhabited.^[ref]

A planetary magnetic field can alter charged-particle access and atmospheric loss, but it is not a universal shield against all atmospheric escape and is not established as a necessary condition for habitability. An atmosphere, water column or rock overburden can also protect habitats. Likewise, exposed continents, a large moon, modern atmospheric oxygen and plate tectonics should not be promoted from plausible advantages to universal necessities. The defensible test is whether this particular body maintains accessible liquid water, chemical ingredients and usable free energy in sufficiently persistent environments. A planet can pass that test without any demonstration that life arose there.

The effective temperature is defined by the Stefan–Boltzmann law applied to a star's total luminosity:Here $R$ is the adopted photospheric radius and ${\sigma }_{\mathrm{SB}}$ is the Stefan-Boltzmann constant. This defines the temperature of a black body radiating the same total flux per unit surface area; it does not assert that every wavelength of the stellar spectrum follows a perfect Planck spectrum.

The effective temperature is defined by $\boxed{L=4\pi R_{\ast }^2\sigma T_{\mathrm{eff}}^4}$, with bolometric luminosity $L$, photospheric radius $R_{\ast }$ and Stefan-Boltzmann constant $\sigma$. It is the blackbody temperature giving the same total surface flux, not a claim that each spectral line forms at that temperature.

A photon gas is a collection of photons described statistically as radiation. Here the equilibrium model is homogeneous, isotropic electromagnetic radiation in vacuum at temperature $T>0$: each photon has energy $E=c|\mathbf{p}|=h\nu$, with momentum $\mathbf{p}$ and frequency $\nu$, and there are two independent transverse states of electromagnetic polarization. An absorbing cavity with walls at a common temperature provides a way to establish this thermal equilibrium. The bulk formulas assume cavity dimensions much larger than the typical wavelength $hc/(k_{B}T)$, so the electromagnetic mode density in a box can replace individual cavity modes by a continuum. ^[ref]

The crucial constraint is whether photon emission and photon absorption, or other number-changing processes, can equilibrate the photon population. When they can, photon number is not an independently conserved quantity and the equilibrium photon chemical potential is ${\mu }_{\gamma }=0$. For each polarization and frequency, the mean occupation number then follows the Bose-Einstein distribution:If number-changing processes are ineffective but number-conserving interactions redistribute energy, a constrained Bose-Einstein distribution can instead have $f=[\exp ((h\nu -{\mu }_{\gamma })/(k_{B}T))-1{]}^{-1}$ with chemical potential ${\mu }_{\gamma }\le 0$ fixed by the number constraint. Thus thermal equilibrium with respect to energy exchange alone does not establish ${\mu }_{\gamma }=0$. The zero value can also persist in an already prepared Planck photon distribution during suitable subsequent evolution; continual photon emission is not required to preserve it. ^[ref]

Combining $f$ with the two-polarization density of states, $g_{\nu }/V=8\pi {\nu }^2/c^3$, gives the Planck photon distribution and the spectral photon energy density:Here $B_{\nu }$ is the spectral radiance in Planck's law; $u_{\nu }$ is stored energy density per unit frequency, not outgoing power per area. Their different units are connected by the full-sphere solid angle and the speed of light $c$. The Boltzmann constant is $k_B$, and $h$ is the Planck constant. ^[ref]

The equilibrium equation of state is particularly simple:The radiation constant $a_R$ is the same constant called $a_{\mathrm{rad}}$ in black-body radiation; $\sigma$ is the Stefan-Boltzmann constant. An outline of the energy density calculation substitutes $x=h\nu /(k_{B}T)$ and uses ${\int }_0^{\infty }{x}^3/(e^x-1)dx={\pi }^4/15$. The radiation pressure follows because the angular average of the squared directional component of photon momentum is $1/3$ for isotropy; the thermodynamic derivation of blackbody radiation pressure supplies an independent route. These formulas describe the radiation contribution, not the cavity walls. The isotropic massless relation $p=u/3$ survives a different isotropic spectrum, whereas the coefficient in $u=a_R{T}^4$ requires the zero-chemical potential Planck photon distribution.

For a computed example with stipulated $T=2.725\mathrm{K}$, the photon energy density is $u\simeq 4.17174\times 10^{-14}\mathrm{J}{\mathrm{m}}^{-3}$ and the radiation pressure is $p\simeq 1.39058\times 10^{-14}\mathrm{Pa}$. The photon number density in the International System of Units isHere $\zeta$ is the Riemann zeta function. The mean photon energy differs from the energy at the peak of a spectral curve: counting photons and weighting them by energy emphasize different frequencies. The calculations use the exact International System of Units values of the Planck constant, Boltzmann constant and speed of light. ^[ref]

A motivating observation is the cosmic microwave background spectrum measured by the Far Infrared Absolute Spectrophotometer (FIRAS) aboard COBE. Its apparatus used a polarizing Michelson interferometer, with one input viewing a roughly seven-degree sky patch and the other an internal blackbody. Scanning mirrors changed the optical path length difference; heat-sensitive detectors recorded the varying power, and Fourier transform spectroscopy recovered the spectral difference. Operation outside the terrestrial atmosphere removed atmospheric emission and absorption, while shielding and cooling the optics in a helium cryostat suppressed the instrument's own radiation. ^[ref]

Absolute calibration came from moving an external blackbody into the sky aperture, filling the same beam used for the sky. This was an iron-loaded epoxy absorber shaped as a re-entrant cone. Repeated reflections made it a close approximation to a blackbody, while thermometers supplied its temperature. Mather and colleagues explicitly evaluated reflections, temperature gradients, radiation leaking through the material and leakage around its edge. These are limits on how faithfully the calibrator realizes thermal equilibrium radiation, rather than details that can be eliminated merely by calling it black. ^[ref]

During sky observations the internal blackbody temperature was adjusted to nearly null the spectral difference, reducing gain and drift errors. Fixsen and colleagues' 1996 analysis used calibrated spectra over wavenumbers $2$ to $21\mathrm{c}{\mathrm{m}}^{-1}$, about $60$ to $630\mathrm{GHz}$. They separated the sky's spatially uniform component, its dipole and Galactic emission, and varied sky cuts and foreground templates as controls. The resulting cosmic microwave background spectrum agreed with a single Planck function, with reported RMS deviations below 50 parts per million of its peak brightness. This is a residual statistic over the analyzed spectrum, not a pointwise relative error bound at every frequency. ^[ref]

The observation motivates using a photon gas as a physical model rather than merely a formal counting exercise. Its interpretation has a clear boundary: agreement with Planck's law constrains departures from an equilibrium spectral shape within the calibrated band and foreground assumptions; it does not directly measure radiation pressure or entropy density of a photon gas. Nor does it establish that today's cosmic microwave background is continually re-equilibrating with matter. A previously thermal spectrum can retain its Planck function shape as expansion of the universe lowers every photon frequency in the same proportion. The usefulness of the model is that one measured spectral temperature, together with its equilibrium assumptions, predicts energy density, radiation pressure, photon number density and entropy density of a photon gas consistently.

The Stefan-Boltzmann constant converts the fourth power of absolute temperature into blackbody radiant exitance:Here $k_B$ is the Boltzmann constant, $h$ the Planck constant and $c$ the speed of light. In the current International System of Units, these three constants have exact defining values, so $\sigma$ is exact as an expression in them and $\pi$. The displayed finite decimal is an approximation, not an exact decimal or a measured uncertainty. The 2022 recommended-constants table lists an ellipsis and an exact uncertainty for $\sigma$.^[ref]