Is it possible to prove the Blackbody radiation law using the fluctuation-dissipation theorem?
Has it been done, or is there some reason why it wouldn't work? I would appreciate if you could point me towards some resources.
I know that Boltzmann first derived it on purely thermodynamic grounds (although that didn't give him the constant of proportionality in $u = a T^4$). That can be obtained with the following derivation, based on the quantisation of energy:
- Start from the Boltzmann probability:
$$p(n) = \frac{ \displaystyle e^{-E_n / k_B T}}{ \displaystyle \sum_{n=0}^{\infty}{e^{-E_n / k_B T}}} \, , $$
where $E_n = n \hbar \omega $
- Mean energy per mode:
$$\overline{E} = \sum_{n=0}^{\infty}{E_n \, p(n)} = \frac{\displaystyle \sum_{n=0}^{\infty} E_n \, \displaystyle e^{-E_n / k_B T}}{ \displaystyle \sum_{n=0}^{\infty}{e^{-E_n / k_B T}}} \, .$$
Using summation rules, we find:
$$\overline{E} = \frac{\hbar \omega}{e^{\hbar \omega / k_B T} -1} \, .$$
- Number of modes per frequency interval is:
$$n(\omega) d \omega = \frac{\omega^2 k_B T}{\pi^2 c^3} d \omega \, .$$
- Hence, the energy density of radiation is:
$$u(\omega) d\omega = \frac{\omega^2 k_B T}{\pi^2 c^3} \frac{\hbar \omega}{e^{\hbar \omega / k_B T} -1} d \omega \, , $$
which is the Planck distribution function.
- Integrate to obtain the Stefan-Boltzmann Law:
$$u = a\, T^4$$ .
Best Answer
The calculation is rather technical and it is usually only carried out when you want to involve the response function (dielectric function) of the material or work at short distances (Wikipedia: near field radiative heat transfer).
For a paper you can check D.Polder and M. Van Hove, Phys.Rev.B 1971 (doi:10.1103/PhysRevB.4.3303), based on Rytov calculations. You can also find a version of it in L.Novotny and B.Hetch "Principles of Nano-Optics" Cambridge.
The idea is the following:
$$\langle j_\alpha(x,t) [j'_\beta(x',t')]^*\rangle\propto\operatorname{Im}(\epsilon)\frac{\hbar\omega^2}{\exp(\hbar\omega/k_{\rm B}T)-1}\delta(x-x')\delta(\omega-\omega')\delta_{\alpha\beta}$$ where $\alpha,\beta=x,y,z$, $\vec{j}$ are the currents, $\epsilon$ is the dielectric function, $\omega$ is frequency and the rest is the usual physical constants.
$$ E_\alpha=\int \mathrm{d}^3r\int \mathrm{d}\omega' \; G_{\alpha\beta}j_\beta(\mathbf r,\mathbf r ',\omega')j_\beta(\mathbf{r}',\omega') $$
where $G_{\alpha,\beta}$ is the dyadic Green tensor for the electric field.
Calculate the Poyting vector, $\langle S\rangle =\langle \vec E \times \vec H\rangle $ where $\vec{H}$ is the magnetic field strenght. Relate each to the currents and use the fluctuation-dissipation relation. Take into account the symmetries of the problem
For the case of two semi-infinite spaces separated by a distance $d$ (for the final formulas, see Wikipedia article on NFRHT), you find that if $\epsilon\to 1$ and $d\to\infty$, you recover that the flux is $\sigma (T^4_1-T^4_2)$, where $T_1,T_2$ are the the temperatures of each semi-space.
Interestingly for some materials this formalism can also recover the emissivity of gray bodies and more complicated responses.