Frequency and wavelength peaks of Planck's law (source code)

= Frequency and wavelength peaks of Planck's law
{title2=$\lambda_{\max}T=b,\quad\nu_{\max}/T=x_\nu k_B/h$}

At a fixed <temperature>, the <frequency> maximizing $B_\nu$ and the <wavelength> maximizing $B_\lambda$ in <Planck's law> are different ways of describing the same <black-body radiation>. They are maxima of different densities. The inverse relation between <frequency> and <wavelength> does not preserve those maxima because its interval-conversion factor varies with <wavelength>.

For a short extremum calculation, substitute $x=h\nu/(k_BT)$ into $B_\nu$, leaving a positive constant times $x^3/(e^x-1)$. Substitute $x=hc/(\lambda k_BT)$ into $B_\lambda$, leaving a positive constant times $x^5/(e^x-1)$. Setting each <derivative> to zero gives
$$
3(1-e^{-x_\nu})=x_\nu,\qquad 5(1-e^{-x_\lambda})=x_\lambda.
$$
Select the positive roots: $x_\nu=2.821439372\ldots$ and $x_\lambda=4.965114232\ldots$. Although $x=0$ solves the rearranged equations, it is an endpoint, not either peak. For $m=3$ or $5$, the function $m(1-e^{-x})-x$ initially increases, its <derivative> $me^{-x}-1$ crosses zero once, and it eventually becomes negative. Hence it has exactly one positive root. Each density tends to zero at both ends, so this root gives its unique maximum. The <bisection method> brackets these roots in $(1,3)$ and $(1,5)$.

Consequently <Wien's displacement law> gives
$$
\lambda_{\max}T=\frac{hc}{k_Bx_\lambda}=b\approx2.897771955\times10^{-3}\,\mathrm{m\,K},
\qquad \frac{\nu_{\max}}T=\frac{k_Bx_\nu}{h}\approx5.878925758\times10^{10}\,\mathrm{Hz\,K^{-1}}.
$$
These constants agree with NIST's tabulation. <Fundamental Physical Constants — Complete Listing by National Institute of Standards and Technology and others>{ref} The <wavelength> peaks at $3000$, $4500$, and $6000\,\mathrm K$ are $965.924$, $643.949$, and $482.962\,\mathrm{nm}$. Doubling <temperature> halves this peak <wavelength>.

At $6000\,\mathrm K$, the <frequency> peak is $\nu_{\max}\approx352.736\,\mathrm{THz}$, whose converted <wavelength> is $c/\nu_{\max}\approx0.849907\,\mu\mathrm m$. The $B_\lambda$ peak is instead $0.482962\,\mu\mathrm m$, corresponding to about $620.737\,\mathrm{THz}$. These numbers are compatible: equal-width <frequency> bins and equal-width <wavelength> bins partition the same radiation differently. Specifying which density was maximized is essential when using a peak to infer <temperature>.

\Image[/planck-frequency-and-wavelength-peaks.png]
{title=One $6000\,\mathrm K$ spectrum, two spectral-coordinate peaks}
{description=The left panel plots the <Planck function> per terahertz and marks its <frequency> peak. The right panel plots the same <black-body radiation> per micrometre and marks its <wavelength> peak. Dashed lines locate the counterpart of the other panel's peak. The converted <frequency> peak at $0.849907\,\mu\mathrm m$ is not the $B_\lambda$ peak at $0.482962\,\mu\mathrm m$: the <change of variables> changes density heights through $c/\lambda^2$.}