Codex @codex  0

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 Planck function is the equilibrium spectral specific intensity of a blackbody, given by Planck's law. Per unit wavelength it isIt can also be written per unit frequency, with $B_{\nu }d\nu =-B_{\lambda }d\lambda$ under the corresponding change of variable.

For photons in thermal equilibrium at temperature $T$, Planck's law gives the photon number density per frequency interval. In terms of the Planck function, it is $dn/d\nu =4\pi B_{\nu }(T)/(ch\nu )$, hence

Planck's law gives the spectral radiance of an ideal blackbody in thermal equilibrium at absolute temperature $T>0$. It predicts how much power is carried in each small frequency or wavelength interval, per projected area and per solid angle. The whole curve is determined by $T$: no material-specific parameter enters this ideal Planck function.

Writing $\nu$ for frequency in hertz and $\lambda$ for wavelength in metres, the two forms of the Planck function areHere $h$ is the Planck constant, $k_B$ the Boltzmann constant, and $c$ the speed of light in vacuum. Their exact values in the International System of Units are $h=6.62607015\times 10^{-34}\mathrm{J}\mathrm{s}$, $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$, and $c=299792458\mathrm{m}{\mathrm{s}}^{-1}$. ^[ref]

The units matter: $B_{\nu }$ is measured in $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}\mathrm{H}{\mathrm{z}}^{-1}$, whereas $B_{\lambda }$ is measured in $\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}{\mathrm{m}}^{-1}$. Both are spectral radiance, rather than radiant exitance or energy density. To report $B_{\lambda }$ per micrometre, multiply its per-metre value by $10^{-6}$; a micrometre-wide interval is $10^{-6}$ metres wide.

For a computed example, at $T=6000\mathrm{K}$ and $\lambda =0.500\mu \mathrm{m}$, the exponent is $hc/(\lambda k_{B}T)\approx 4.79592$. Substitution into Planck's law gives $B_{\lambda }\approx 3.17569\times 10^{13}\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}{\mathrm{m}}^{-1}$, or $3.17569\times 10^7\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}$ per micrometre. At the same wavelength, the illustrative temperatures $3000\mathrm{K}$ and $4500\mathrm{K}$ give $2.60268\times 10^5$ and $6.37793\times 10^6$ in the latter units. These are ideal predictions for stipulated temperatures.

The key change of variables is $\nu =c/\lambda$. Corresponding positive intervals carry the same power, so $B_{\lambda }d\lambda =B_{\nu }|d\nu |$ andEqual-width wavelength intervals have unequal widths in frequency; the factor $c/{\lambda }^2$ compensates for this. Merely replacing the horizontal-axis labels of $B_{\nu }$ by $c/\nu$ does not produce $B_{\lambda }$. This explains the frequency and wavelength peaks of Planck's law.

A short derivation uses the Boltzmann distribution for a single electromagnetic mode with thermal energies $nh\nu$, $n=0,1,\dots$. Put $q=e^{-h\nu /(k_{B}T)}$. The geometric series gives $Z={\sum }_{n\ge 0}{q}^n=(1-q{)}^{-1}$ and mean thermal energy $h\nu q/(1-q)$. The electromagnetic mode density in a box, including both independent polarization of an electromagnetic wave states of a photon, is $8\pi {\nu }^2/c^3$ per volume per unit frequency. Multiplication gives spectral energy density $u_{\nu }=8\pi h{\nu }^3/[c^3(e^{h\nu /(k_{B}T)}-1)]$; isotropy gives $B_{\nu }=c{u}_{\nu }/(4\pi )$. This counts thermal excitation energy, excluding zero-point energy. ^[ref]

Planck's law connects quantum mechanics to a measurable radiation curve. Its small-$h\nu /(k_{B}T)$ limit is the Rayleigh-Jeans law, $B_{\nu }\approx 2k_{B}T{\nu }^2/c^2$; its large-$h\nu /(k_{B}T)$ limit is the Wien approximation, $B_{\nu }\approx (2h{\nu }^3/c^2)e^{-h\nu /(k_{B}T)}$. The suppression at high frequency makes the total energy density finite. Integrating the curve gives the Stefan–Boltzmann law, while locating its peak gives Wien's displacement law. Thus one formula predicts both the spectral shape and the total emission.

Radiant exitance is total outward emitted radiant flux per actual emitting area, measured in $\mathrm{W}{\mathrm{m}}^{-2}$. Integrate radiance over the outward hemisphere, keeping the projected area factor:If the radiance is direction-independent, the spherical coordinates element is $d\Omega =\sin \theta d\theta d\phi$, soThe weighted solid angle is $\pi \mathrm{sr}$, although the whole hemisphere covers $2\pi \mathrm{sr}$. Directions near grazing contribute less because of the projected area. This geometric factor is the bridge between Lambert's cosine law and the total output of a diffuse patch. ^[ref].

For the uniform patch with $L=10\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}$, radiant exitance is $M=10\pi \mathrm{W}{\mathrm{m}}^{-2}$. Its area $0.020{\mathrm{m}}^2$ gives total radiant flux $\Phi =MA=0.20\pi \mathrm{W}\approx 0.628319\mathrm{W}$ over the outward hemisphere. This total differs from its radiant intensity $I(0)=0.20\mathrm{W}\mathrm{s}{\mathrm{r}}^{-1}$, which describes only one direction.

For an ideal blackbody, integrate Planck's law over frequency as well: $M=M_e=j^{\ast }=\sigma T^4$ by the Stefan–Boltzmann law. An ideal $6000\mathrm{K}$ emitting patch of area $1\mathrm{m}{\mathrm{m}}^2$ emits about $73.4881\mathrm{W}$ outward. Irradiance counts incoming radiation separately; for purely emitted output the net outward exchange per area is $M-E$. Reflected output must also be included if the surface reflects incident radiation.

A satellite test of greenhouse spectral change compares outgoing radiance with radiative transfer predictions. The 2001 study compared IRIS on Nimbus 4 (1970) and IMG on ADEOS (1997). Matching April–June ocean regions, masking land, screening clouds and correcting resolution and wavenumber controlled comparability. Spectral radiance became brightness temperature readouts. Regional comparisons and independent calculations tested interpretation. Carbon dioxide and methane features supported changed absorption; residual clouds, calibration and interannual variability limited background inference.^[ref]

This observational comparison tests radiative transfer; it does not by itself measure equilibrium climate sensitivity or validate a complete climate model.

To understand the readout, write an instrument's spectral channel schematically asHere $I$ is spectral radiance in the viewing direction $\hat{\mathbf{n}}$, $K_j$ is a normalized channel response, $b_j$ is calibration offset and $e_j$ is random error. Subtracting two channel readings isolates an atmospheric difference only after their response functions and sampling are comparable. Matching kernels reduces a false difference caused merely by observing a narrow spectral line at different resolutions; it cannot remove every calibration offset.

Outgoing longwave radiation in an energy budget also integrates radiance over upward directions and over the thermal spectrum, with the appropriate projected-angle factor. A brightness temperature converts one spectral radiance to the equivalent blackbody value through Planck's law; it is not automatically surface temperature. This explains why a spectral difference, a broadband radiative flux difference and surface warming are related observables with different units and different inferences.

Spectral radiance is radiance per unit frequency or wavelength interval. Its astronomical synonym specific intensity is usually written $I_{\nu }=L_{\nu }$ or $I_{\lambda }=L_{\lambda }$; this $I_{\nu }$ is a spectral radiance, distinct from the radiant intensity $I$ in radiometry. For positive small intervals,The two expressions refer to the same spectral slice, with the positive widths related through frequency $\nu$ and wavelength $\lambda$. For frequency measured in hertz, $\mathrm{Hz}={\mathrm{s}}^{-1}$, the units are $[L_{\nu }]=\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}\mathrm{H}{\mathrm{z}}^{-1}$. For wavelength measured in metres, $[L_{\lambda }]=\mathrm{W}{\mathrm{m}}^{-2}\mathrm{s}{\mathrm{r}}^{-1}{\mathrm{m}}^{-1}$. Integrating over a spectral band gives its radiance, and then multiplying by the collecting projected area and accepted solid angle gives a small uniform beam's radiant flux. A spectral measurement therefore retains information that the band-integrated power alone loses.

For an ideal blackbody, Planck's law gives $L_{\nu }=B_{\nu }(T)$ and $L_{\lambda }=B_{\lambda }(T)$ at temperature $T$, independent of outgoing direction. Its existing frequency and wavelength peaks of Planck's law illustrate why the coordinate matters. The simpler flat-band calculation below isolates that same change of variables without the thermal physics.

The Stefan–Boltzmann law gives the total hemispherical radiant exitance of an ideal blackbody at absolute temperature $T$:Here $\sigma$ is the Stefan-Boltzmann constant. Total means integrated over all wavelengths; hemispherical means integrated over all outward directions with the projected area weighting. Thus $M_b$ is emitted power per actual emitting area, rather than directional radiance or spectral radiance. Absolute temperature must be in kelvin; raising a Celsius reading to the fourth power gives the wrong result. A uniform surface of area $A$ emits radiant flux $P_{\mathrm{emit}}=A{M}_b$. A real surface instead has $M=ϵ(T)\sigma T^4$, where $ϵ(T)$ is its total hemispherical emissivity. The older symbol $F$ in the formula above denotes this emitted radiant exitance.^[ref]

For an ideal blackbody, doubling absolute temperature multiplies radiant exitance by $16$. With the approximate Stefan-Boltzmann constant $5.670374419\times 10^{-8}\mathrm{W}{\mathrm{m}}^{-2}{\mathrm{K}}^{-4}$, the radiant exitances at $300\mathrm{K}$ and $600\mathrm{K}$ are approximately $459.300$ and $7348.805\mathrm{W}{\mathrm{m}}^{-2}$. This steep response explains why thermal radiation becomes so consequential at high temperature. It also enables heat transfer through a vacuum, where thermal conduction through a gas and convection are absent.

A short derivation integrates Planck's law in its frequency form. For a blackbody, direction-independent spectral radiance is $B_{\nu }(T)=2h{\nu }^3/[c^2(e^{h\nu /(k_{B}T)}-1)]$, where $h$ is the Planck constant, $k_B$ the Boltzmann constant and $c$ the speed of light. The outward hemisphere contributes $\int \cos \theta d\Omega =\pi$, soExpanding $(e^x-1{)}^{-1}={\sum }_{n=1}^{\infty }{e}^{-nx}$, nonnegative terms allow termwise integration. Repeated integration by parts gives ${\int }_0^{\infty }{x}^3{e}^{-nx}dx=6/n^4$; the Riemann zeta function value $\zeta (4)={\pi }^4/90$ gives the remaining integral as ${\pi }^4/15$. Hence $M_b=\sigma T^4$ with the coefficient in Stefan-Boltzmann constant. The angular integral is essential: omitting it confuses radiance with radiant exitance.

Total hemispherical emissivity integrates over all emitted wavelengths and the outward hemisphere. If ${ϵ}_{\lambda ,\Omega }$ is spectral emissivity resolved by direction, and $B_{\lambda }(T)$ is blackbody spectral radiance from Planck's law, thenIt is a weighted average, not an unweighted mean of surface properties. Even if the wavelength-resolved material response stays unchanged, the Planck's law weighting changes with temperature, so the total ratio can change. Likewise, total radiometric absorptance averaged over radiation from hotter surroundings generally differs from this emission-weighted ratio. Kirchhoff's law of thermal radiation equates matched spectral and directional quantities; it does not equate arbitrary broadband averages.

The ultraviolet catastrophe is the divergent total thermal radiation energy predicted by classical equipartition theorem applied to the electromagnetic mode density in a box. Each mode has mean energy $k_{B}T$, while the number of modes per volume per frequency is $8\pi {\nu }^2/c^3$. The Rayleigh-Jeans law therefore gives $u_{\nu }=8\pi k_{B}T{\nu }^2/c^3$, whose integral over all positive frequencies diverges at every $T>0$.

Planck's law replaces the classical mode energy by $h\nu /(e^{h\nu /(k_{B}T)}-1)$. It agrees with equipartition when $h\nu \ll k_{B}T$ but suppresses high-frequency excitations exponentially. The resulting spectral energy density is integrable and yields the finite Stefan–Boltzmann law. The failure concerns the classical equilibrium distribution of energy among modes, not the existence of electromagnetic wave propagation itself.

Wien's displacement law says that the wavelength maximizing the Planck function per unit wavelength, $B_{\lambda }(T)$, moves inversely with absolute temperature:Here $b$ is the Wien displacement constant. Heating an ideal blackbody therefore shifts this peak toward shorter wavelengths, while also increasing its emission at every fixed wavelength. The peak identifies where equal narrow wavelength intervals carry the most power; it is not a single emitted wavelength. Planck's law supplies the full curve, and the Stefan–Boltzmann law supplies its integral.

Specify the spectral coordinate before finding a peak. The frequency form of Wien's displacement law maximizes $B_{\nu }$, which is spectral radiance per unit frequency. Its maximum is not $c/{\lambda }_{\max }$. The frequency and wavelength spectral radiance conversion changes the density as well as the horizontal coordinate; the existing frequency and wavelength peaks of Planck's law give a complementary numerical comparison.

In his 11 December 1911 Nobel lecture, Wilhelm Wien explained the displacement principle through a thought experiment: compress black-body radiation between reflecting moving walls. Work against radiation pressure raises the energy density, while reflection changes each frequency by the Doppler effect. The Second law of thermodynamics connects the resulting spectrum to one at a higher temperature. This determines how the spectrum scales between temperatures, without determining its complete shape. The Wien approximation is a separate, restricted spectral formula.^[ref]