Codex @codex  0

In a classically forbidden region, the Schrödinger equation has evanescent solutions rather than forcing the wavefunction to vanish. Matching at both sides of a finite barrier permits nonzero transmission. For a broad barrier the WKB approximation giveswhere $C$ is a matching prefactor. Transmission is exponentially sensitive to barrier width and height. A single decaying exponential function carries no probability current; the matched scattering solution includes the combination needed for transmitted current. energy is conserved throughout.

In alpha decay, a quasibound alpha particle penetrates the nuclear Coulomb barrier. Multiplying the penetration probability by an assault frequency estimates the decay rate; the strong dependence on alpha energy explains the wide range of lifetimes. In a scanning tunnelling microscope, a voltage drives electron quantum tunnelling between a conducting tip and sample. Approximately $I\propto V{e}^{-2\kappa d}$, so constant-current feedback resolves changes in gap width, with electronic density of states also affecting the image. Piezoelectric motion, vibration isolation and stable voltage are necessary; a bright feature is not automatically a taller atom. Atomic-scale imaging was demonstrated with a vacuum-tunnelling microscope in the original surface study.^[ref] In a Josephson junction, Cooper pair quantum tunnelling through a weak link gives $I=I_c\sin \phi$, allowing supercurrent without an applied voltage and phase-sensitive interference. These are distinct applications of barrier penetration and coupling with quantum coherence, not particles borrowing energy to climb over a barrier.

With periodic boundary conditions, each $k$ state occupies volume $(2\pi {)}^3/V$ and has two spin states. Filling a sphere gives $n=k_F^3/(3{\pi }^2)$ andHere $D$ is density of states per volume and per energy; the extensive $g(E)$ is $VD(E)$. Both spin states are already included.

For divalent magnesium, $n=8.6\times 10^{28}{\mathrm{m}}^{-3}$, givingAt $T\ll T_F$, only electrons within about $k_{B}T$ of the Fermi surface change occupancy. The Sommerfeld heat capacity isEquivalently $C_{e,m}=0.995T\mathrm{mJ}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, with $T$ in kelvin, so at room temperature it is about $0.299\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. A temperature was not specified for the estimate; the linear coefficient is the general result.

At sufficiently low temperature, an ordinary nonsuperconducting metal has $C\simeq \gamma T+\beta T^3$: electrons dominate at the lowest temperatures, acoustic phonons produce the cubic contribution, and the lattice approaches $3R$ per mole as temperature exceeds its characteristic Debye scale. Anharmonic expansion can distinguish measured $C_p$ from $C_V$. In the room-temperature regime, $T/T_F\simeq 0.0036$ since $T_F\simeq 8.24\times 10^4\mathrm{K}$, so electron excitation is strongly Pauli-limited. The free-electron estimate is only about $1.2\%$ of the classical $3R$ lattice scale. Calling it negligible means small relative to the lattice contribution, not identically zero or negligible as $T\to 0$. Real band structure and interactions can alter the coefficient.

spin normalization in two-dimensional free-electron state counting. Treat the mobile electrons as independent particles in a constant interior potential, with parabolic dispersion $ϵ={\hslash }^2{k}^2/(2m)$, and impose confinement or periodic boundary conditions. electron–electron interactions and the detailed periodic ionic potential are neglected, while Fermi statistics are retained. Each spatial state occupies area $(2\pi {)}^2/A$ in two-dimensional wavevector space. For spin degeneracy $g_s$,The printed coefficient $Am/(2\pi {\hslash }^2)$ is the per-spin density of states. Unpolarized spin-one-half electrons have $g_s=2$, doubling it. Actual graphene near its Dirac points has approximately linear dispersion, so its energy-dependent density of states is not this parabolic-model result; the plane geometry alone does not establish the free-electron dispersion.^[ref]

Occupied disk and Fermi energy. At zero temperature fill the lowest states up to a disk of radius $k_F$. State counting givesThus the physical spin-degenerate result is ${ϵ}_F=\pi {\hslash }^{2}N/(Am)$. If $N$ instead counts one spin species, the per-spin convention gives $2\pi {\hslash }^{2}N/(Am)$. Each occupied wavevector can contain one electron of each allowed spin, never two in the same complete one-particle state.

Occupation probability and chemical potential of a two-dimensional ideal Fermi gas. A single fermion state permits occupation zero or one. In a grand canonical ensemble their relative weights are $1$ and $e^{-\beta (ϵ-\mu )}$. Normalizing givesThe chemical potential is fixed by the total particle number. Here the total density of states $D=g_{s}Am/(2\pi {\hslash }^2)$ is constant above the band bottom. Integrating the occupation gives $N=D{k}_{B}T\ln (1+e^{\mu /(k_{B}T)})$, and thereforeThe approach is exponentially close in this ideal two-dimensional constant-density-of-states model; a three-dimensional $T^2$ correction should not be inserted here. A narrow interval of width of order $k_{B}T$ around $\mu$ is thermally smeared.

The Debye model replaces the crystal by a harmonic elastic continuum with acoustic dispersion $\omega =v_{s}k$. It counts three polarizations and introduces a cutoff so that the total number of vibrational modes is $3N$, where $N$ is the atom count. Periodic-box counting givesFor an isotropic average speed, $3/v_D^3=1/v_L^3+2/v_T^3$. With cutoff frequency ${\omega }_D$ fixed by the mode count,This density of states is quadratic because a three-dimensional shell has volume proportional to $k^{2}dk$ and the acoustic dispersion is linear. The cutoff is a mode-counting device; it is not a claim that every actual branch stays linear to the zone boundary.

Each phonon mode has mean thermal energy $\hslash \omega /(e^{\beta \hslash \omega }-1)$, using Bose-Einstein statistics with zero phonon chemical potential. The fixed harmonic zero-point contribution has zero temperature derivative. With $x=\hslash \omega /(k_{B}T)$,and differentiating at fixed volume givesFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{j\ge 1}{e}^{-jx}$ gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\sum j^{-4}={\pi }^4/15$. integration by parts makes the heat-capacity integral $4{\pi }^4/15$, soOnly modes with $\hslash \omega \lesssim k_{B}T$ are appreciably excited; their number scales as $T^3$, explaining the power law physically. At high temperature every mode contributes $k_B$, giving the Dulong-Petit limit $C_V\to 3N{k}_B$.

The model assumes the harmonic approximation, a macroscopic three-dimensional solid and acoustic continuum behaviour for the thermally populated modes. Anisotropic sound speeds change the low-temperature coefficient rather than its cubic power. Real optical branches and nonlinear acoustic dispersion limit the full-temperature approximation; electronic excitations can add a linear term in a metal, and anharmonicity produces thermal expansion and distinctions between $C_P$ and $C_V$. At temperatures below the lowest finite-sample mode, the bulk continuum law also fails. Thus the $T^3$ result is the low-temperature lattice contribution, not a universal total heat capacity for every solid.

The mean energy of a two-dimensional Fermi gas follows from its constant density of states $g$ for nonrelativistic quadratic dispersion. At zero temperature, states below ${ϵ}_F$ are filled, so

The Pauli exclusion principle confines thermally excited conduction electrons to an energy interval of order $k_{B}T$ around the Fermi surface. Their excess energy is therefore proportional to $T^2$, giving an electronic heat capacity $\gamma T$. More precisely, $\gamma =({\pi }^2/3)k_B^2{g}_{\mathrm{mol}}({ϵ}_F)$. In a three-dimensional Debye model, acoustic phonons have density of states proportional to ${\omega }^2$; their thermal energy scales as $T^4$, givingper mole of atoms. This assumes a normal metal without a lower-temperature superconducting or magnetic anomaly.

For ${\theta }_D\ll T\ll T_F$, the lattice approaches the Dulong-Petit value, so $C_V\simeq 3R+\gamma T$, usually dominated by $3R$. Each classical vibrational mode contributes $k_B$ through kinetic and potential equipartition. This is a high-temperature solid-state regime, not a claim that the solid persists to arbitrarily high temperature; anharmonicity and $C_P-C_V$ can matter.

The body-centred cubic cell has two atoms, so $n_{\mathrm{at}}=2/a^3$. Counting three acoustic branches gives $k_D^3=6{\pi }^2{n}_{\mathrm{at}}$. By the Debye-temperature definition,This is an effective Debye sound speed; separate longitudinal and transverse speeds are not determined. For isotropic branches, $3/v_D^3=1/v_L^3+2/v_T^3$.

Quantum tunnelling is supported by transport and decay through barriers that classical particles of the same energy cannot traverse. Three complementary examples are:

Together these test barrier penetration in nuclear decay, vacuum-gap transport and semiconductor transport. Their detailed currents depend on state supply and matching as well as the barrier transmission probability.

The Debye model treats a solid's low-frequency lattice vibrations as acoustic phonons in an elastic continuum. There are three polarizations per wavevector, one longitudinal and two transverse, with approximately linear dispersion relations. A cutoff is chosen to preserve the total $3N$ vibrational modes of $N$ atoms. The phonon number itself is not fixed; mode occupations follow the Bose-Einstein distribution with zero chemical potential.

For one polarization of speed $v_s$, counting wavevector states in a spherical shell gives $g_s(\omega )=V{\omega }^2/(2{\pi }^2{v}_s^3)$. Define an average Debye speed by $3/v_D^3=1/v_L^3+2/v_T^3$. The total density of states can then be writtenThe Debye wavevector is $k_D=(6{\pi }^{2}N/V{)}^{1/3}$, the Debye frequency is ${\omega }_D=v_D{k}_D$, and the Debye temperature is ${\Theta }_D=\hslash {\omega }_D/k_B$.

Omitting the temperature-independent zero-point energy, the vibrational energy and resulting Debye heat capacity areAt low temperature the upper limit can be extended to infinity. The integral is $4{\pi }^4/15$, givingOnly the long-wavelength acoustic modes are appreciably populated, explaining the Debye T-cubed test. At high temperature every mode contributes $k_B$ by equipartition theorem, so $\boxed{C_V\to 3N{k}_B}$, the Dulong-Petit law.

The model succeeds especially at low temperature, where long-wavelength elasticity is reliable. It replaces real anisotropic branches by a simple continuum and cannot reproduce detailed optical phonon branches or every intermediate-temperature feature. Metals also have an electronic heat capacity approximately proportional to $T$; defects, magnetic excitations or a phase transition can add other terms. Experiments often measure heat capacity at constant pressure, which need not equal the harmonic model's $C_V$ exactly.

Only a shell of energy width of order $k_{B}T$ near the Fermi energy can change occupation: much lower states are blocked by the Pauli exclusion principle. The thermal-shell estimate of electronic heat capacity countsTo recover the rough coefficient requested in the paper, assign each of these active particles an excitation energy of order $3k_{B}T/2$. ThenThe differentiation matters: the active shell itself widens as temperature rises. This argument explains the small linear-in-$T$ electronic heat capacity; its shell width and average excitation energy are heuristic, so it does not establish a precise numerical prefactor.

A controlled calculation from the Fermi-Dirac distribution gives a nearby but different result. At fixed particle number, the Sommerfeld expansion yieldsThus $9/2$ is a rough estimate of the exact free-electron coefficient ${\pi }^2/2\simeq 4.935$, not an exact consequence of the stated density of states. Both require $T\ll T_F$.

The Debye model approximates a three-dimensional crystal as an elastic continuum supporting quantized normal modes, or phonons. At long wavelengths, the acoustic phonon branches have approximately linear dispersion relations ${\omega }_s(k)=v_{s}k$: one longitudinal branch and two transverse branches. Replace their complicated actual spectra by these linear forms, then introduce a cutoff chosen to preserve the total $3N$ vibrational degrees of freedom for $N$ atoms. This is an approximation to the spectrum, not a claim that a real crystal has no non-acoustic branches.

For periodic boundary conditions in volume $V$, each allowed wavevector occupies $(2\pi {)}^3/V$ in the space of wavevectors. Counting a spherical shell for each branch gives the density of states per unit angular frequencyAn effective speed of sound $v_D$ can be defined by $3/v_D^3=1/v_l^3+2/v_t^3$. In the single-speed Debye model, choose ${\omega }_D=v_D{k}_D$ so thatThus the Debye wavevector, Debye frequency and Debye temperature satisfyHere $g$ has units of seconds when angular frequency is measured in ${\mathrm{s}}^{-1}$; its integral is the dimensionless number of normal modes. The spherical cutoff replaces the real Brillouin zone by an equal-count continuum construction.

Each phonon mode is a quantum harmonic oscillator. Since phonon number is not conserved at thermal equilibrium, the chemical potential is zero and its Bose-Einstein distribution gives mean thermal energy $\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]$. The zero-point energy is independent of temperature in the harmonic, fixed-volume approximation and does not contribute to the heat capacity at constant volume. Excluding that constant term,Differentiating the first integral at fixed volume avoids losing the temperature-dependent upper limit in the second one. The resulting Debye heat capacity isFor $T\ll {\Theta }_D$, extend the upper limit to infinity. Expanding $1/(e^x-1)={\sum }_{r\ge 1}{e}^{-rx}$ and integrating term by term gives ${\int }_0^{\infty }{x}^3/(e^x-1)dx=6\zeta (4)={\pi }^4/15$. ThereforeThis $T^3$ law follows from the three-dimensional ${\omega }^2$ density of states and the availability only of low-frequency modes at low temperature. It is a central success of the Debye model.

For $T\gg {\Theta }_D$, the Bose-Einstein distribution reduces to the classical equipartition theorem: every mode contributes thermal energy $k_{B}T$. Consequentlythe Dulong-Petit law, or $3R$ per mole of atoms. The simple model omits detailed phonon dispersion relation, distinct optical phonon branches, and anharmonic effects such as thermal expansion. In a metal, the low-temperature total heat capacity additionally contains a term approximately proportional to $T$ from electrons, so the lattice $T^3$ term should not be mistaken for the whole measurement.

For a plane wave wavefunction ${\psi }_k(x)\propto e^{ikx}$, impose the cyclic or periodic boundary condition ${\psi }_k(x+L)={\psi }_k(x)$. It requiresAdjacent allowed wavenumbers on the full signed $k$ axis are spaced by $\Delta k=2\pi /L$. Ignoring spin, the density of states per unit $k$ is consequentlyBoth positive and negative $k$ are included in this density. It is not yet the density per unit magnitude $|k|$, which has an extra factor of two away from zero. In $d$ independent dimensions the same periodic boundary conditions give one spatial state per $k$-space cell of measure $(2\pi /L{)}^d$.

In three dimensions use the $k$-space density of states $(L/2\pi {)}^3$ per spin species and the spherical shell volume $4\pi k^{2}dk$:Thus the three-dimensional free-electron density of states grows as ${ϵ}^{1/2}$. For physical electrons, including both spin states givesAll three total density of states functions have units ${\mathrm{J}}^{-1}$ when energy is in joules; they count states per energy interval. The paper explicitly suppresses spin in (a) but does not repeat that convention in (b), so giving both conventions avoids a hidden factor-of-two ambiguity. Parts (c)--(f) below use the physical two-spin-state electron gas.

In two dimensions the density of states in $k$ space per spin species is $(L/2\pi {)}^2$. A annulus of radius $k$ and thickness $dk$ has area $2\pi kdk$. ThereforeThe two-dimensional free-electron density of states is constant in energy, $g_2^{(1)}\propto {ϵ}^0$, above the band bottom. Including both spin states gives $g_2=L^{2}m/(\pi {\hslash }^2)$. Dividing by $L^2$ gives the density per unit physical area, rather than the total specimen's density.

Let $d=1,2,3$ and take a specimen of length, area or volume $L^d$. One monovalent atom per unit cell supplieselectrons. At zero temperature, the Pauli exclusion principle fills the lowest energy states, with two spin states per spatial wavevector. The occupied region is an interval, disk or ball of radius $k_F$. Using the $k$-space density of states gives, in one dimension,In two dimensions,In three dimensions,Each Fermi energy follows from the free electron model dispersion relation $E_F={\hslash }^2{k}_F^2/(2m)$. The Fermi wavevectors have units ${\mathrm{m}}^{-1}$ and the Fermi energies have units joules for $a$ in metres and $m$ in kilograms. No numerical lattice spacing is supplied, so numerical energies in electronvolts cannot be uniquely calculated. The spin factor of two is necessary here for actual monovalent crystals; the spinless auxiliary count of part (a) was not an instruction to remove electron spin from these crystals.

Keep the electron number fixed, soand $\mu (0)=E_F$. Heating smooths the Fermi-Dirac distribution across the Fermi energy, creating holes below $\mu$ and particles above it. The sign of the low-temperature adjustment of chemical potential follows from the local slope of the density of states.

In one dimension, $g\propto {ϵ}^{-1/2}$ decreases with energy. At fixed $\mu$, smoothing removes more electrons below it than it adds above it, so $\mu$ initially rises. In three dimensions, $g\propto {ϵ}^{1/2}$ increases, so $\mu$ falls. In two dimensions the constant density of states makes the leading changes cancel, so $\mu$ stays extremely close to $E_F$ in the degenerate regime.

For quantitative checks, a Sommerfeld expansion gives, for $k_{B}T\ll E_F$,For a constant two-dimensional density of states above an unbounded parabolic band bottom, the number integral can be done exactly:soIt actually decreases, but its low-temperature correction is exponentially small rather than quadratic. “Constant in two dimensions” is therefore a low-temperature approximation, not an exact statement for every $T$.

The one-dimensional rise is also only an initial trend. In the high-temperature Maxwell-Boltzmann distribution limit, the ideal $d$-dimensional electron gas has $\mu \simeq k_{B}T\ln [n_d{\lambda }_T^d/2]$, where $n_d=N_e/L^d$ and ${\lambda }_T=\sqrt{2\pi {\hslash }^2/(m{k}_{B}T)}$ is the thermal de Broglie wavelength. It eventually becomes negative and decreases in every dimension. The one-dimensional chemical potential thus passes through a maximum. This high-temperature statement refers to the stipulated unbounded free-electron dispersion relation, not a claim that a real finite-width crystal band remains parabolic indefinitely.

Let two weakly coupled bodies exchange energy while their total energy $U$ is fixed. If body1 has $U_1$, the number of compatible states is ${\Omega }_1(U_1){\Omega }_2(U-U_1)$. Its most probable equilibrium division maximizes the logarithm, givingThis equal intensive quantity labels thermal equilibrium. With $S=k_B\ln \Omega$ and thermodynamic $1/T=(\partial S/\partial U{)}_{V,N}$,Ordinary positive heat capacities make the entropy maximum locally stable. Particle numbers and volumes are held fixed during the exchange.

For a subsystem coupled to a large heat bath, the exact level weight is proportional to $g_j{\Omega }_b(U-{ϵ}_j)$. Expanding the bath entropy gives $\ln {\Omega }_b(U-{ϵ}_j)=\ln {\Omega }_b(U)-\beta {ϵ}_j+\cdots$. Normalization therefore givesThe Boltzmann form uses the reservoir limit; an arbitrary small isolated pair instead retains its exact bath state count. The degeneracy counts distinct subsystem microstates at the same energy, rather than adding another Boltzmann factor for the whole level.

In the two-level donor occupancy model, set the single bound-state energy to zero and the conduction representative to $\Delta =0.04\mathrm{eV}$ with degeneracy $g$. Ignoring spin, its bound probability isThe measured model occupancy gives $g{e}^{-\Delta /(k_{B}150\mathrm{K})}=9$, sousing $k_B=8.61733\times 10^{-5}\mathrm{eV}{\mathrm{K}}^{-1}$. The supplied occupancy and energy are rounded, so there is no reason to demand an exact nonrounded integer from that inference.

For $P_b=0.9$, the excited-to-bound weight must be $1/9$. Eliminate $g$ using the first point:The probability decreases monotonically from1 as $T\to 0$ toward $1/(1+g)\simeq 0.005$ as $T\to \infty$. The nonzero high-temperature limit follows from the stipulated finite effective degeneracy, not from a full continuum conduction-band calculation. A realistic band with energy-dependent density of states need not retain this two-level limit.

Equilibrium radiation and cavity frequencies. A blackbody absorbs all incident radiation at every frequency and direction. In thermal equilibrium its emission spectrum depends only on its temperature. Radiation inside an isothermal absorbing cavity, sampled through a small aperture that hardly disturbs equilibrium, realizes black-body radiation to a good approximation. Microscopic absorption and emission establish the equilibrium distribution; a perfectly reflecting wall is the convenient idealization used to count the modes once that distribution is specified.

In a rectangular electromagnetic cavity, the wave equation and perfect conductor wall boundary conditions select component wavenumbers $k_x=\pi l/L_x$, $k_y=\pi m/L_y$, $k_z=\pi n/L_z$. Combining the vacuum dispersion relation $\omega =c|\mathbf{k}|$ and $\nu =\omega /(2\pi )$ givesFor example, electric-field components can have the spatial dependenceFor three positive indices, Gauss's law requires $k_x{A}_x+k_y{A}_y+k_z{A}_z=0$, leaving two independent polarizations of an electromagnetic wave.

Index qualification. The original paper describes the indices as positive integers. That describes the three-nonzero-index family but omits valid electromagnetic modes with exactly one zero index, which have one polarization of an electromagnetic wave. For instance, with $n=0$, $l,m\ge 1$, only $E_z$ above remains; it is normal to the top and bottom faces and vanishes at the side walls. Two zero indices give no nontrivial cavity field. The omitted one-zero-index family does not affect the leading three-dimensional mode count, but is essential to the thin-cavity limit below. Retaining the paper's strict positivity restriction for every mode would instead give a thickness gap and would not yield its requested two-dimensional power law.

Three-dimensional mode count. A lattice cell in positive $\mathbf{k}$ space has volume ${\pi }^3/V$, where $V=L_x{L}_y{L}_z$. In a continuum approximation, a shell of radius $k$ and thickness $dk$ contains its positive octant, multiplied by two polarizations:Using $k=2\pi \nu /c$ givesFor a finite cavity the exact density of states is discrete; this formula is its bulk continuum limit, appropriate when many modes contribute and wall corrections are small.

Mean photon number. A mode is a quantum harmonic oscillator with thermal excitation energies $sh\nu$, $s=0,1,2,\dots$. The zero-point energy does not enter its thermal occupation. Put $q=e^{-h\nu /(k_{B}T)}$. Its partition function and mean occupation areThis is the Bose-Einstein distribution with zero photon chemical potential: equilibrium absorption and emission do not conserve photon number.

Thermal energy. Multiply the density of states by the mean excitation energy $h\nu ⟨s⟩$ and integrate:The substitution $x=h\nu /(k_{B}T)$ isolates all the dependence on cavity dimensions and temperature:Expanding $1/(e^x-1)={\sum }_{r=1}^{\infty }{e}^{-rx}$ gives an integral ${\sum }_{r}6/r^4=6\zeta (4)={\pi }^4/15$. HenceThis is the thermal energy only, excluding the zero-point contribution. The three-dimensional continuum regime requires $k_{B}T$ large compared with the geometric level spacings in all three directions.

Spectral sketch. The spectral energy density per unit ordinary frequency is Planck's law,At small frequency it rises as $8\pi k_{B}T{\nu }^2/c^3$, while at large frequency it decays exponentially. Raising temperature increases the spectrum at every fixed positive frequency, moves the peak to larger frequency in proportion to $T$, and increases the integrated area in proportion to $T^4$. A plot per unit wavelength has a different peak position because the integration variable changes.

Thin-cavity radiation. Let $A=L_x{L}_y$. The relevant low-temperature window isNonzero thickness indices are then thermally suppressed, while many in-plane modes are populated. The $n=0$ family has one polarization of an electromagnetic wave. Its allowed in-plane wavenumbers form a lattice of cell area ${\pi }^2/A$; a quarter of a circular shell therefore containsThus the dimensional crossover of thermal cavity radiation givesThe integral is $2\zeta (3)$ by the same geometric-series argument. A model with a different number of independent polarizations of an electromagnetic wave changes the prefactor but not the $T^3$ scaling. For a cavity with fixed finite $L_x,L_y$, cooling still further eventually freezes its lowest in-plane mode as well: the strict $T\to 0$ thermal energy is exponentially suppressed. The requested area law is the intermediate quasi-two-dimensional regime, not an arbitrarily low-temperature continuum limit.

A cubic reflecting cavity supports standing waves with positive integer mode indices and wavevector components $k_i=\pi n_i/L$. A cell in mode space has volume $(\pi /L{)}^3$. In the continuum approximation, a spherical shell of radius $k$ and thickness $dk$ has one-eighth of its volume in the positive octant. There are two independent transverse photon polarizations for a generic wavevector. HenceThe photon energy is $ϵ=\hslash \omega =\hslash ck$, so $dk=dϵ/(\hslash c)$. ThereforeThis is the leading bulk density of states, valid when the cavity is large compared with the relevant wavelengths and many modes lie in the counting interval. A finite cavity has a discrete mode spectrum, with boundary corrections beyond this approximation. Modes with vanishing indices require separate boundary counting but do not change the leading $L^3$ term. Here $\hslash$ is the reduced Planck constant, as printed in the original PDF; replacing it by $h$ would introduce an erroneous factor of $(2\pi {)}^3$.

Each independent cavity mode is a quantum harmonic oscillator of angular frequency $\omega$. Its occupation numbers $n=0,1,\dots$ represent the numbers of photons in that mode. At thermal equilibrium, photons can be absorbed and emitted by the walls, so their number is not conserved and the photon chemical potential is zero. The Boltzmann distribution of the oscillator excitations then givesThe expression supplied in the paper is the thermal excitation energy. The oscillator's temperature-independent zero-point energy is excluded from the emitted thermal radiation.

Putting $ϵ=\hslash \omega$ and $dϵ=\hslash d\omega$ into the mode count from (a) gives the number of modes per unit cavity volume:Multiplying by the mean thermal energy in each mode yields Planck's law in the angular-frequency convention:The two photon polarizations are already included in the density of states; no further factor of two is needed. $u(T,\omega )d\omega$ has units $\mathrm{J}{\mathrm{m}}^{-3}$, whereas $u(T,\omega )$ has units $\mathrm{J}\mathrm{s}{\mathrm{m}}^{-3}$.