Stefan-Boltzmann and Wien Displacement Laws
The Planck law hands you the entire spectral distribution of blackbody radiation — a curve giving the intensity at every frequency. Almost nobody carries that whole curve around. Two of its summary numbers do the daily work of astrophysics, and both were found by experiment before Planck explained why they held. The Stefan-Boltzmann law reads off the area under the Planck curve: the total power radiated per unit area, which scales as the fourth power of temperature. Wien's displacement law reads off the location of the summit: the wavelength at which the curve peaks, which slides inversely with temperature. Area and peak — one integral and one derivative of the same function. Their reach is enormous precisely because each collapses a spectrum into a single measurable quantity. The steep means a factor-of-two rise in temperature brightens a surface sixteenfold, which is why a star's effective temperature is defined through it and why the luminosities on the Hertzsprung-Russell diagram fan out over so many decades. The color-temperature relation lets you take a source's temperature from its hue alone: a cool red-giant photosphere, the glare of an O star, and the cosmic microwave background all announce their temperatures by where they peak. And because both and are built purely from , , and times a pure number, measuring them was one of the first routes to the quantum constants themselves. This article derives both laws from the Planck function, verifies every constant, and shows where they earn their keep.
- The Stefan-Boltzmann law is the area under the Planck curve: a blackbody radiates a flux , with built from , , .
- Wien's displacement law is the peak of that curve: , so hotter bodies peak at shorter (bluer) wavelengths.
- Both are exact consequences of one spectrum: integrate the Planck function to get ; differentiate it and set the result to zero to get .
- The comes from four powers of pulled out of the frequency integral; the pure number is from .
- Together they yield effective temperatures, radii, luminosities, and color temperatures — from the Sun (, peak ) to the CMB (, peak ).
- The frequency peak and the wavelength peak sit at different photon energies — a Jacobian trap that catches everyone once.
- Key idea: Integrate the Planck function over frequency and solid angle to get ; differentiate it and set the result to zero to get . Both the scaling and the color-temperature relation are exact properties of one spectrum.
Where this lives today
The two laws are the routine arithmetic of modern observation. Every stellar effective temperature in a survey catalogue — Gaia's billion-plus stars, the SDSS spectroscopic samples — is a Stefan-Boltzmann inversion at heart, and the color indices those surveys report are Wien's law in broadband disguise. Exoplanet science leans on both: a planet's equilibrium temperature balances absorbed stellar flux against its own re-radiation, the first cut on the habitable zone, and the same sets Earth's climate budget. At the cosmological frontier, COBE's FIRAS measured the CMB to be a blackbody at with deviations below one part in — the most perfect blackbody ever observed, and direct confirmation that the early Universe was in thermal equilibrium. Current CMB experiments (the Simons Observatory, CMB-S4) work right at the Wien peak because that is where the signal is strongest against foregrounds. And in stellar interiors, the radiation constant sets the radiation pressure that holds massive stars up.
The physical picture
Picture the Planck spectrum as a single bell-shaped hill whose whole shape is fixed by one dial: the temperature. Two numbers summarize the hill. Its total area — integrated over all frequencies — is the energy the surface pours out each second, the Stefan-Boltzmann quantity. The position of its summit is the Wien quantity. Turn the dial up and two things happen at once: the hill grows dramatically taller (its area climbs as , its peak height even faster, as ) and its summit slides toward shorter wavelengths. The first figure shows this hill for the Sun and marks both numbers on the same curve.
The second face is the raw numerology of the exponents. Why and not ? Substitute in the frequency integral and four factors of fall out — three from the prefactor, one from — leaving a pure number behind. (Photon number density alone scales as ; each photon's mean energy supplies the fourth power.) Wien's law, meanwhile, says the curve is self-similar: written against the Planck shape is universal, so its peak sits at a fixed and the peak wavelength must scale as . The second figure turns that scaling into a straight line and hangs four real blackbodies on it.
Because the color-temperature relation is monotonic, you do not even need to resolve the peak to use it: a single flux ratio through two broadband filters pins the temperature of any thermal source — the basis of photometric temperatures for stars too faint to spectrograph.
The core equations
The headline result is the total emergent flux. Integrate the Planck intensity over all frequencies and over the outward hemisphere and every trace of the spectral shape collapses into one number times .
where: — emergent radiative flux from a blackbody surface (W m⁻²); — Stefan-Boltzmann constant, ; — absolute temperature (K). Double the temperature and the surface brightness rises by . The constant is not fundamental. It is a fixed bundle of the genuinely fundamental constants times the pure number that came out of the integral, packaged with the hemisphere factor.
where: — Boltzmann constant, ; — Planck and reduced Planck constant (); — speed of light. The two forms are identical: substituting turns into . The companion headline is the peak location. In wavelength, the product of the peak wavelength and the temperature is a constant.
where: — wavelength of the peak of ; — Wien displacement constant; is the nonzero root of . Reading: at the peak sits at , blue-green. The peak has a different location if you plot per unit frequency instead of per unit wavelength, because the two are densities in different variables. The frequency peak obeys its own transcendental root.
where: — frequency of the peak of ; is the nonzero root of . Reading: and do not satisfy — see the warning below. Everything above is extracted from one integrand, the Planck function, here in its per-unit-frequency form.
where: — specific intensity per unit frequency (W m⁻² Hz⁻¹ sr⁻¹), isotropic for a blackbody. This is the curve whose area gives (1) and whose summit gives (4). The flux out of a surface is times the frequency-integrated intensity, the being the projected solid angle of the outward hemisphere.
where: — frequency-integrated Planck intensity (W m⁻² sr⁻¹); — the value of over the outward hemisphere. Inside an enclosure the relevant quantity is the isotropic radiation energy density, which carries the same but a different constant.
where: — radiation energy density; — radiation constant. The factor separates a one-sided emergent flux from an isotropic interior energy density. Finally, the astrophysical workhorse: multiply the surface flux by the area of a sphere to get a star's bolometric luminosity, which defines its effective temperature.
where: — bolometric luminosity (W); — stellar radius; — effective temperature, the blackbody temperature reproducing the star's output at radius .
The goal is to start from nothing but the emergent flux of an isotropic blackbody surface and the Planck intensity, carry out the frequency and angle integrals, and arrive at both the law and the closed form — showing exactly where the fourth power comes from (four factors of from the substitution) and where the pure number comes from (a standard Bose-Einstein integral). The constant is built, not quoted.
Step 1. Write the flux emerging from a surface as the intensity weighted by the projection factor and integrated over the outward hemisphere: , where is the frequency-integrated intensity. Why: a detector on the surface sees each ray foreshortened by , and over a hemisphere is exactly , not .
Step 2. Insert the Planck function and non-dimensionalize with , so and : . Why: the contributes three powers of and the contributes a fourth, which is the entire origin of the fourth power — the remaining integral is a pure number.
Step 3. Evaluate the dimensionless integral as a Gamma-times-zeta product: . Why: expanding and integrating term by term gives , and is a standard result.
Step 4. Assemble and identify the constant: . Why: collecting every constant in front of literally defines ; plugging in CODATA returns . *Verified: , and from CODATA constants; matches Carroll & Ostlie Eq. (3.18).*
Go deeper · Derivation B · 1 minWien displacement law by maximizing the Planck curve
The goal is to locate the summit of the Planck spectrum rigorously: switch the intensity to a per-unit-wavelength density using the correct Jacobian, differentiate with respect to wavelength, set the derivative to zero, and reduce the condition to a transcendental equation whose root delivers the displacement constant . Along the way it becomes clear why the wavelength and frequency peaks obey different equations and sit at different photon energies.
Step 1. Convert to with the Jacobian of , namely , giving . Why: intensity is a density, so it must be transformed with the Jacobian, and it is this weighting that will make the wavelength peak differ from the frequency peak.
Step 2. Set . Writing and differentiating the product gives, after clearing common factors, the condition , equivalently . Why: at the maximum the two competing trends — the prefactor pushing the peak to short wavelength and the Bose factor cutting off the short-wavelength end — exactly balance, which is what the derivative-equals-zero condition encodes.
Step 3. Solve the transcendental equation numerically. The trivial root is a minimum; the physical root is , found by iterating to convergence. Why: the equation cannot be solved in closed form, but fixed-point iteration converges quickly to ; the analogous frequency condition has root , and the two differ because of the Jacobian in Step 1.
Step 4. Extract the constant: , so . Why: with fixed and known, the product is a pure constant — the whole content of the displacement law. *Verified: the root of is , and ; the frequency root is giving . Matches Ryden Eq. (2.30)-type values.*
Go deeper · Derivation C · 1 minRadiation energy density and the factor c over 4
The goal is to derive the interior radiation energy density and pin down why its constant differs from the surface constant by exactly . This closes the loop between how much energy sits in a radiation field and how much streams out of a hole in it, supplying the that runs through stellar-interior and cosmology calculations.
Step 1. The isotropic energy density is the intensity integrated over all steradians and divided by : . Why: energy density counts photons flowing in every direction, so it uses the full solid angle and converts the streaming intensity to a stored density by dividing by the propagation speed .
Step 2. Substitute the frequency-integrated Planck intensity from the previous derivation (note so ): . Why: the same spectral integral that built builds ; only the geometric prefactor changes from (one-sided flux) to (all-sky density).
Step 3. Compare the two constants: , and equivalently the emergent flux from a cavity is . Why: the is the standard effusion factor — averaging over the outward hemisphere and normalizing over all directions gives , the same factor that relates gas pressure to flux through a hole. *Verified: and reproduce ; consistent with Rybicki & Lightman Ch. 1.*
Limiting cases
The two laws are the integrated and extremal summaries of one spectrum, so their behavior in each corner of the Planck curve is worth reading off directly.
| Case | Result | Meaning |
|---|---|---|
| Rayleigh-Jeans (low frequency, ) | Classical limit; integrating this alone over all diverges (the ultraviolet catastrophe), so the finite demands the full Planck cutoff. | |
| Wien tail (high frequency, ) | The spectrum falls exponentially, which is what makes the peak sharp and finite; historically Wien's own 1896 formula. | |
| Effective temperature | Defines for any star as the blackbody temperature giving its observed luminosity at radius , even when the true spectrum is only approximately Planckian. | |
| Energy density vs flux | The radiation constant relates interior energy density to ; a factor separates isotropic density from one-sided flux. | |
| Wavelength peak vs frequency peak | vs | Different transcendental roots, so ; the mismatch is the Jacobian , not an error. |
Worked examples
1 · Effective temperature and color of the Sun (fully worked)
Given the solar luminosity and radius , find the effective temperature from the Stefan-Boltzmann law, then predict the peak wavelength from Wien's law.
Step 1. Invert equation (8) for temperature: . Assemble the denominator first: , and multiplying by gives .
Step 2. Divide and take the fourth root: , so . This is the defined effective temperature of the Sun — the standard IAU value is , an exact match.
Step 3. Feed that temperature into Wien's wavelength law (3): . That lands in the blue-green, entirely consistent with the Sun's G2 spectral class. T_eff = 5772 K and λ_max = 502 nm. Note the Sun looks white, not blue-green, because the eye integrates the whole broad curve — the peak is only where the density is largest.
2 · Temperature of the cosmic microwave background (fully worked)
The CMB is the most perfect blackbody known, with . Find its peak in frequency and in wavelength, and confirm they disagree.
Step 1. Frequency peak from equation (4): , squarely in the microwave band that COBE/FIRAS measured.
Step 2. Wavelength peak from equation (3): .
Step 3. Check the Jacobian trap: naively converting the frequency peak gives , not the wavelength peak — the ratio is , exactly as the two roots predict. 160 GHz and 1.06 mm are both correct; they are peaks of different densities.
3 · Radius of a red supergiant from both laws
Betelgeuse has bolometric luminosity and photospheric temperature (peak at , deep red). Estimate its radius.
Step 1. Solve equation (8) for radius: , with and .
Step 2. Denominator , so , i.e. .
Step 3. A star with the same as a cool dwarf but times its luminosity must be enormous — which is exactly why the coolest, most luminous stars occupy the supergiant branch of the Hertzsprung-Russell diagram. R ≈ 800 R_⊙.
Physics of the portrait
Behind the two tidy laws stand four physicists and two decades of argument. The empirical relations came first — a measured slope, a measured shift — and the man who turned the peak into a theorem, and set the stage Planck would later fill, was Wilhelm Wien.
Origins · 1 minWilhelm Wien and how this was discovered
Wilhelm Wien worked at the Physikalisch-Technische Reichsanstalt in Berlin, whose precise blackbody measurements would soon force the quantum revolution. In 1893 he proved the displacement law by a thermodynamic argument: subject cavity radiation to a slow adiabatic compression, track how each wavelength Doppler-shifts off the moving walls and how the temperature rises, and the whole spectral shape must be a universal function of the single combination — so the peak wavelength has to scale as , with no microscopic model. This "displacement" of the peak is the law that bears his name. In 1896 he went further and guessed a full spectral form, , which fit the short-wavelength data but failed in the infrared — the failure that Planck patched in 1900 by interpolating between Wien's exponential tail and the Rayleigh-Jeans slope, forcing him to quantize energy. Wien's law thus sits exactly at the hinge of classical and quantum physics: it is derivable from thermodynamics alone, yet its constant secretly contains . He received the 1911 Nobel Prize "for his discoveries regarding the laws governing the radiation of heat." The total-power side is older: Josef Stefan inferred empirically in 1879, and his student Ludwig Boltzmann derived it thermodynamically in 1884 by treating radiation as a gas with pressure — which is why that law carries both names.
Going further
Texts. Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 3 (blackbody radiation, effective temperature) · Ryden, Introduction to Cosmology, 2nd ed., Ch. 2 and Ch. 8 (CMB blackbody) · Rybicki & Lightman, Radiative Processes in Astrophysics, Ch. 1 (thermal radiation, the factor) · Dodelson, Modern Cosmology, Ch. 1-2 · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 10 (stellar photometry) for the color-temperature machinery in a galaxy-survey context.
Provenance. Equations and constants are grounded in Carroll & Ostlie (Ch. 3), Ryden (Ch. 2), and Rybicki & Lightman (Ch. 1). A verifier should recompute ; confirm the transcendental roots from and from ; reproduce , , and from CODATA ; and check the worked numbers , , , .
Prerequisites. Blackbody Radiation and the Planck Law · Specific Intensity and Radiative Flux · Kirchhoff's Laws and Thermal Emission. Leads to. Local Thermodynamic Equilibrium · The Saha Ionization Equation.
Ask the tutor: Ask me why the Sun's frequency peak lands in the infrared while its wavelength peak is green — the Jacobian that trips up everyone.
Is the Stefan-Boltzmann law only true for a perfect blackbody?
Strictly, is the flux from an ideal blackbody surface. A real ("grey") body radiates , where the emissivity measures how close it comes to ideal. For stars this is folded into the definition of effective temperature: rather than track a wavelength-dependent emissivity, one defines as the temperature of the blackbody emitting the star's actual total flux from its actual radius. So holds exactly, by construction, with all the non-ideal detail absorbed into the single number .
Where does the pure number π⁴/15 actually come from?
It is the value of the dimensionless integral left after the temperature and physical constants are factored out of the Planck integral. Expand the Bose factor as a geometric series and integrate term by term: each term gives , and summing over gives . So the in traces back to the zeta value , not to any geometric circle.
Why is there a factor of π in F = πB but a factor of 4π/c in the energy density?
Both come from integrating over solid angle, but they answer different questions. Flux asks how much energy crosses the surface per second, so each ray is weighted by its projection and integrated over the outward hemisphere: . Energy density asks how much energy is stored per unit volume, weighting every direction equally over the full steradians and dividing by the propagation speed: . Combining the two gives the effusion relation , hence .
The Sun peaks at 502 nm (green) — so why isn't the Sun green?
Because your eye responds to the whole broad spectrum, not to the single wavelength of maximum spectral density. The Planck curve at is wide and only gently peaked; it delivers copious red, green, and blue, and the brain integrates all of it to "white" (sunlight defines white). The peak wavelength is a real feature of the density curve, but it is not "the color" of the source — and it depends on whether you plot per wavelength or per frequency, the Sun's frequency peak sitting near in the infrared.
Can I use Wien's law on a source that isn't a perfect blackbody?
With care. Wien's law locates the peak of the Planck function; a source with strong emission or absorption lines, a non-thermal component, or significant dust reddening will have a peak displaced from . For a stellar photosphere the continuum is close enough to Planckian that the peak (or, more robustly, a broadband color index) gives a good color temperature, though it can differ modestly from the effective temperature. The safe practice is to fit a color index calibrated against models rather than read a single peak off a low-resolution spectrum.
How were σ and b used to measure the fundamental constants?
Turn the logic around. Measure in the laboratory (total radiated power versus temperature) and independently measure (peak position versus temperature), and you have two equations, and , in the unknowns and (with known). Planck did exactly this in 1900: from the fitted radiation data he extracted the first reasonable values of both and , and from with the gas constant he even obtained Avogadro's number. The blackbody spectrum was historically one of the first precision routes to the quantum of action.