Codex @codex  0

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.

A diatomic molecule has one internal bond-stretching vibrational coordinate. Near its equilibrium separation, approximate the bond by a spring of stiffness $K$ and use the reduced mass $\mu =m_1{m}_2/(m_1+m_2)$. Its frequency is $\omega =\sqrt{K/\mu }$ and its quantum harmonic oscillator levels are ${ϵ}_r=(r+1/2)\hslash \omega$.

For independent molecules in thermal equilibrium, set $x=\hslash \omega /(k_{B}T)={\theta }_v/T$. The single-molecule canonical partition function is a geometric series:Using $U=-\partial \log (z_v^N)/\partial \beta$, where $\beta =1/(k_{B}T)$, givesPer mole replace $N{k}_B$ by $R$. The zero-point energy is present in $U_v$ but is temperature-independent and contributes nothing to the heat capacity.

For $T\ll {\theta }_v$, occupation of the first excited vibrational state is exponentially suppressed: $C_{V,v}\simeq N{k}_B{x}^2{e}^{-x}\to 0$. The vibrational degree of freedom is frozen out. For $T\gg {\theta }_v$, $C_{V,v}\to N{k}_B$. This recovers classical equipartition: the quadratic kinetic term and quadratic elastic term each contribute $k_B/2$ per molecule.

The vibrational heat capacity of a diatomic gas must be distinguished from the total heat capacity. In a dilute gas with rotations already active and vibration still frozen, translation contributes $3R/2$ and rotation $R$, so molar $C_V\simeq 5R/2$. Once vibration is active it adds $R$, giving $C_V\simeq 7R/2$. For an ideal gas of a fixed number of intact molecules, $C_P=C_V+R$, so the corresponding values are $7R/2$ and $9R/2$. Very low temperatures also freeze rotations, while sufficiently high temperatures bring anharmonicity, dissociation or electronic excitations; the harmonic plateau is restricted to the range where those effects are negligible.

The photoelectric effect, black-body radiation and the Davisson–Germer experiment challenge different parts of a fully classical picture. Classical electromagnetic waves successfully describe propagation, interference and diffraction, while classical particles follow trajectories determined by forces. The difficulty is combining that picture with the observed exchange of energy between radiation and matter, the equilibrium distribution of radiation, and electron scattering. The three experiments expose the roles of discrete energy exchange and coherent matter waves.

In a photoelectric effect experiment, monochromatic light illuminates a metal cathode in vacuum and a second electrode collects emitted electrons. Varying the collecting voltage measures a photocurrent; reversing it until collection ceases determines the stopping potential. The observations distinguish two controls: at fixed frequency above threshold, increasing light intensity increases the emission rate, whereas the largest emitted kinetic energy is determined mainly by frequency. There is a material-dependent photoelectric threshold frequency, and emission above threshold does not require the long charging-up delay predicted by a simple picture of gradual classical energy accumulation at low intensity.

A classical wave carries continuously distributed energy, with intensity proportional to squared field amplitude. On a fully classical absorption picture, increasing intensity would be expected to increase the energy delivered to an electron, and sufficiently long exposure could accumulate the escape energy even below the observed threshold. It does not naturally produce the measured universal straight-line relation between stopping potential and light frequency. In 1905, Albert Einstein instead treated light-energy transfer as localized quanta. His argument distinguished the fastest escaping electrons from electrons losing energy before escape.^[ref]

In the single-photon regime, conservation of energy gives the Einstein photoelectric equationwhere $\Phi$ is the surface work function. The photon supplies energy $h\nu$; the escape cost is paid first, and the remainder sets the upper electron-energy edge. Increasing intensity at fixed frequency increases the number of available photons, so it raises current without raising this edge. Below ${\nu }_0$, one photon cannot pay the escape cost. In an actual measurement, stable surface preparation, suppression of unwanted higher-frequency light, and the photoelectric contact-potential correction are needed to test the slope $h/e$ reliably. Robert Andrews Millikan's experiment used selected mercury spectral lines, freshly exposed metal surfaces and an electrometer to make such a quantitative test.^[ref] The interpretation is for weak single-photon illumination; multiphoton photoemission at sufficiently intense fields is a different regime.

Black-body radiation tests equilibrium physics rather than the escape of individual electrons. A small aperture in an isothermal absorbing cavity approximates an ideal blackbody: radiation entering the hole undergoes repeated interactions with the walls, and the equilibrium spectrum depends on temperature rather than wall material. A spectrally resolving detector outside the aperture measures how emitted power varies with frequency or wavelength. Experimentally the spectrum has a finite peak, is strongly suppressed at short wavelengths, shifts toward shorter wavelengths when the cavity is heated, and has integrated emission proportional to $T^4$.

For a large cavity, the electromagnetic mode density in a box, including two independent transverse polarizations, is $g(\nu )=8\pi {\nu }^2/c^3$ per volume per unit frequency. The classical equipartition theorem assigns mean energy $k_{B}T$ to each mode, leading to the Rayleigh-Jeans lawThis agrees at low frequencies but diverges as frequency increases; integrating it over all frequencies gives infinite energy density at any nonzero temperature. This is the ultraviolet catastrophe, a failure of classical equipartition applied to the electromagnetic modes.

Max Planck's essential modification was that energy at a frequency is exchanged in units $h\nu$. In a modern oscillator description, the thermal excitation levels are $nh\nu$, $n=0,1,\dots$. Using the Boltzmann distribution and summing a geometric series givesMultiplying by the mode density yields Planck's law for the spectral energy density,When $h\nu \ll k_{B}T$, expansion of the denominator recovers the classical law. When $h\nu \gg k_{B}T$, the large excitation cost suppresses occupancy exponentially, removing the divergence. Its integral gives the Stefan–Boltzmann law, while its peak gives Wien's displacement law. Thus quantization repairs the spectral shape, total emission and temperature dependence together. The displayed thermal energy excludes the oscillator's temperature-independent zero-point energy.

The Davisson–Germer experiment supplies the complementary evidence that matter cannot always be treated as classical particles with definite scattering trajectories. In 1927, Clinton J. Davisson and Lester Halbert Germer sent a narrow, accelerated electron beam normally onto a single nickel crystal in vacuum. A movable collector measured the current of scattered electrons as direction and accelerating voltage were varied; a retarding potential selected electrons retaining approximately the incident energy. Sharp angular beams depended on crystal orientation and beam energy, with a prominent feature near $54\mathrm{V}$ and detector colatitude $50^{\circ }$. These are scattering directions and currents, not photographs of electron paths.^[ref]

The de Broglie wavelength predictsAt $U=54\mathrm{V}$, this gives $\lambda \simeq 0.167\mathrm{nm}$, comparable to atomic spacings. A crystalline arrangement supplies regularly spaced scattering sites, so scattered probability amplitudes can add constructively in selected directions. For a suitable family of crystal planes, Bragg's law expresses the condition as $2d\sin {\theta }_B=n\lambda$, with ${\theta }_B$ measured relative to those planes. The detector colatitude is not automatically this angle. The measured angular structure and its variation with electron momentum are described by electron diffraction with the de Broglie wavelength. A classical point-particle deflection model has no phase wavelength $h/p$ and does not supply that interference condition.

The common quantum account preserves the successes of wave propagation and particle detection while changing what is quantized. Photons exchange energy in frequency-dependent units; a thermal distribution of those excitations gives the blackbody spectrum. Electrons have a wavefunction whose amplitudes interfere on scattering from a crystal, while detection transfers localized charge. The Born rule turns squared amplitudes into event probabilities. These results invalidate a fully classical description of matter and radiation; the photoelectric energy law alone does not establish every property of a quantized electromagnetic field, but the three phenomena together provide complementary evidence for the quantum description.

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.

Use the canonical commutation relations $[\hat{x},{\hat{p}}_x]=[\hat{y},{\hat{p}}_y]=i\hslash$, with cross-coordinate commutators zero. For either coordinate,Thus $2\alpha \beta \hslash =1$. To obtain the required Hamiltonian, the ratio of its position operator and momentum operator coefficients must be preserved:where the positive root follows from the requested signs. ConsequentlyThe annihilation operators are dimensionless; $\alpha$ has units of inverse length and $\beta$ has units of inverse momentum. Operators belonging to different coordinates commute automatically. Writing the number operators as ${\hat{N}}_j={\hat{a}}_j^{†}{\hat{a}}_j$, substitution givesThe constant is the sum of the two zero-point energies of this two-dimensional isotropic harmonic oscillator.

For excitation energy $\Delta =\hslash \omega$, the geometric series gives the thermal partition function of a quantum harmonic oscillator with its zero-point energy subtracted:Thus $U\to 0$ at low temperature and approaches $k_{B}T$ at high temperature. Restoring the zero-point term multiplies $Z$ by $e^{-\beta \Delta /2}$ and adds $\Delta /2$ to $U$, without changing the thermal occupation probabilities.

To count electromagnetic modes, put a large cavity in a periodic box of volume $V$. One transverse polarization has wavevector-state density $V/(2\pi {)}^3$; there are two independent polarizations. In the shell $k\to k+dk$,With $\omega =ck$, the electromagnetic mode density in a box becomesThe continuum expression is the large-cavity leading term; using standing-wave boundary conditions gives the same bulk count. Both photon polarizations must be included. This is density per angular frequency, not per ordinary frequency.

Each mode is an oscillator with zero photon chemical potential in thermal equilibrium. Multiplying its mean thermal energy by the mode density per volume gives the angular-frequency Planck energy densityAt small $\hslash \omega /(k_{B}T)$ it reduces to the Rayleigh–Jeans result $k_{B}T{\omega }^2/({\pi }^2{c}^3)$; exponential suppression at high frequency prevents a divergent thermal energy. With $x=\hslash \omega /(k_{B}T)$,The integral is the stipulated ${\pi }^4/15$. This is a volume energy density, $a{T}^4$; for isotropic cavity radiation the outward black-surface flux is $cu/4={\sigma }_{\mathrm{SB}}{T}^4$, with ${\sigma }_{\mathrm{SB}}={\pi }^2{k}_B^4/(60{\hslash }^3{c}^2)$. Distinguishing flux, spectral density and total volume energy avoids extra factors of $c$, $4$ or $2\pi$.

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.

Entropy at absolute zero. The third law of thermodynamics concerns the low-temperature limit of equilibrium states. In the zero-entropy formulation, the entropy of a perfect crystal with a nondegenerate ground state tends to zero as $T\to 0$. The statistical interpretation is $S=k_B\ln \Omega$: a unique accessible ground state has $\Omega =1$. The entropy-difference formulation of the third law of thermodynamics states that entropy differences associated with regular changes of mechanical parameters tend to zero as absolute zero is approached. These statements require equilibrium and appropriate ground-state assumptions; a frozen disordered material or a degenerate ground state can carry residual entropy.

For a phase of fixed composition and pressure,A finite limiting entropy constrains this integral. Ordinary equilibrium heat capacities therefore vanish as temperature approaches zero, with familiar examples $C\propto T^3$ for three-dimensional acoustic phonons and $C\propto T$ for a normal degenerate Fermi liquid. A constant positive heat capacity extrapolated to zero would give a logarithmically divergent entropy difference and is impossible in this regime. Discontinuous transitions and anomalous limits need separate treatment rather than extending a high-temperature approximation indefinitely.

Gibbs free energy. For fixed particle numbers,Thus, when $S\to 0$, the Gibbs free energy approaches its zero-temperature value with a horizontal tangent:The equality $G(0,p)=H(0,p)$ uses the enthalpy $H=U+pV$; neither limiting value needs to be zero. The low-temperature Gibbs free-energy expansion makes the approach quantitative. For $C_p=a{T}^r$ with $r>0$ and $S(0)=0$,The second derivative is ${\partial }_T^{2}G=-C_p/T\le 0$, although it need not tend to zero. For a reaction between such equilibrium phases, ${\partial }_T\Delta G=-\Delta S\to 0$, and $\Delta G$ and $\Delta H$ approach one another. If a genuine residual entropy $S_0$ is retained instead, $G=G(0,p)-S_{0}T+\cdots$ has limiting slope $-S_0$. A horizontal tangent cannot be asserted without the zero-entropy condition.

Why helium remains a liquid. Helium atoms attract each other weakly and are exceptionally light. Their zero-point energy is large compared with the binding that would localize them in a solid. Both helium-4 and helium-3 can therefore remain liquid down to the lowest temperatures at low pressure; substantial applied pressure is needed to solidify them. This is not a violation of the third law. A quantum liquid can have a unique low-entropy ground state, and zero-point motion does not itself imply thermal entropy. The third law does not require every substance to crystallize at atmospheric pressure.

A helium-4 atom has integer spin zero and is a boson. Below the lambda point near $2.17\mathrm{K}$ at saturated vapor pressure, its liquid is superfluid. Its strong interactions mean that it is not quantitatively an ideal Bose gas, and the condensate fraction need not equal the superfluid fraction. The two-fluid model of superfluidity separates an entropy-carrying normal component of a superfluid from the superfluid component, whose entropy is zero. At sufficiently low temperature the dominant thermal excitations are acoustic phonons: their linear dispersion in three dimensions gives $U-U_0\propto V{T}^4$, $C_V\propto V{T}^3$ and $S\propto V{T}^3$. Higher-energy rotons are exponentially suppressed. Thus the normal component of a superfluid and entropy disappear as $T\to 0$, while superfluid flow and quantum motion can persist.

A helium-3 atom has spin $1/2$ and is a fermion. Well below its Fermi temperature, but above its superfluid transition, liquid helium-3 is a degenerate interacting Fermi liquid. Only quasiparticles close to the Fermi surface are thermally excited, giving $C_V\propto T$ and $S\propto T$. At temperatures of order a millikelvin, helium-3 also becomes superfluid through formation of Cooper pairs of fermions. This is not Bose-Einstein condensation of individual helium-3 atoms. Pairing reorganizes the low-energy excitation spectrum; in the fully gapped bulk B phase, the fermionic quasiparticle contribution to heat capacity is exponentially suppressed at temperatures well below the pairing gap. Gapless collective modes can still give thermal terms following a power law. The entropy nevertheless tends to zero. The normal Fermi-liquid law cannot be extrapolated through its pairing transition unchanged.^[ref]

Helium-3 also illustrates the importance of nuclear spin entropy. Over an intermediate low-temperature range, localized spins in the solid can retain approximately $k_B\ln 2$ of entropy per atom, while the degenerate liquid has much less. Reversible conversion of liquid to the higher-entropy solid then absorbs heat; adiabatic compression along the melting curve can cool the sample. This is Pomeranchuk cooling. At still lower temperature, interactions between the nuclear spins remove this nearly free-spin entropy, so a constant $k_B\ln 2$ cannot be extrapolated to absolute zero.^[ref]

Liquefaction and cooling methods. A gas must be below its liquid-gas critical point to form a distinct liquid phase by compression. Compression above the critical temperature alone cannot achieve ordinary liquid-vapour coexistence. Practical liquefaction removes the compression heat, cools the compressed gas in a heat exchanger and uses expansion to lower its temperature further.

A throttling valve operates approximately at constant specific enthalpy when kinetic-energy changes and heat leakage are negligible. The Joule-Thomson coefficient isusing consistent molar quantities, or consistent total quantities for a fixed amount. A pressure drop cools the gas if ${\mu }_{\mathrm{JT}}>0$. An ideal gas has zero Joule-Thomson coefficient; real intermolecular interactions are essential. In particular, room-temperature helium warms on throttling, so it must first be brought into a region of positive Joule-Thomson coefficient. Expansion through an engine or turbine can instead remove work and produce approximately isentropic cooling.

A regenerative liquefier sends the cold returning gas through a counterflow heat exchanger to precool the incoming compressed gas before its next expansion. Expansion engines and external precooling can supply the initial cooling; repeated operation then produces liquid. Modern helium liquefaction commonly combines work-producing expansion with a final throttle. ^[ref]

Once liquid exists, pumping away vapour lowers its saturation pressure and boiling temperature. The escaping vapour carries latent heat away, cooling the remaining liquid. At still lower temperatures, a dilution refrigerator exploits the different entropy carried by helium-3 in concentrated and dilute helium mixtures; adiabatic demagnetization and Pomeranchuk cooling provide other routes.

Approaching zero. The unattainability formulation of the third law states that regular cooling operations with finite resources cannot reach absolute zero in finitely many stages. As the entropy available to remove becomes small and thermal couplings weaken, each method loses its ability to take the final step. An isentropic cooling trajectory starting at nonzero entropy cannot reach a common zero-entropy endpoint through regular finite changes of parameters. A fully general finite-time claim also requires assumptions about available controls and heat-transfer rates. Quantum liquids and refrigeration methods can approach $0\mathrm{K}$ while remaining consistent with vanishing equilibrium entropy; they do not provide a finite-operation route to absolute zero.

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}$.

Zero-point energy is the ground-state energy above the classical minimum, remaining even at zero temperature. For a harmonic oscillator it is $\hslash \omega /2$; simultaneous vanishing position and momentum fluctuations is forbidden by the uncertainty principle.

A concrete test is the Casimir interaction: changing conducting boundaries changes electromagnetic ground-state modes and produces a measurable force. Ideal parallel plates have pressureA torsion-pendulum experiment measured the attraction of closely spaced conducting surfaces, testing this boundary-dependent interaction.^[ref] Such measurements test changes in ground-state energy, not an arbitrary absolute additive vacuum energy; real-material and thermal corrections must be included when comparing precision data.

The classical minimum is $E_{\mathrm{cl}}=0$, attained at $x=p=0$. The quantum oscillator has $E_0=\hslash \omega /2$, its zero-point energy; simultaneous zero position and momentum uncertainties are forbidden.

Let $\ell =\sqrt{\hslash /(m\omega )}$ and $\xi =x/\ell$. The first three normalized harmonic oscillator eigenfunctions areThey have respectively zero, one and two nodes, with even, odd and even parity. Their phases may be reversed without changing the state.^[ref]

On a smooth wavefunction $f$, the position operator and momentum operator giveThus the canonical commutation relation is $\boxed{[X,P]=i\hslash I}$. Since $X$ and $P$ are self-adjoint operators and $m,\omega$ are real, the adjoint operator of $a_{-}$ isExpanding the products in the stated order is essential:Subtracting gives $\boxed{[a_{-},a_{+}]=2m\hslash \omega I}$. These identities initially hold on smooth rapidly decaying functions, which contain the quantum harmonic oscillator energy eigenstates used below.

The factorization also proves a lower bound: for a normalized energy eigenstate $H\psi =E\psi$,Therefore $E\ge E_0:=\hslash \omega /2$. The commutator product rule, applied to $H=a_{+}{a}_{-}/(2m)+E_0$, givesSo a nonzero application of the creation operator $a_{+}$ raises the energy eigenvalue by $\hslash \omega$, while the annihilation operator $a_{-}$ lowers it by the same amount.

Begin with any energy eigenstate, whose existence is assumed. Repeated lowering must terminate: otherwise it would eventually produce a nonzero energy eigenstate below $E_0$, contradicting the lower bound. More explicitly, if $(E-E_0)/(\hslash \omega )=N+r$ with integer $N\ge 0$ and $0\le r<1$, the first $N$ lowering steps are nonzero. If $r>0$, one further step is also nonzero by its positive squared norm and has energy $E_0+(r-1)\hslash \omega <E_0$, an impossibility. Thus $r=0$ and the last state satisfies $a_{-}{\psi }_0=0$. This provesEvery level exists as well: $\Vert a_{+}{\psi }_n{\Vert }^2=2m(E_n+\hslash \omega /2)=2m\hslash \omega (n+1)>0$, so repeated raising from ${\psi }_0$ never vanishes. The same factorization bounds the energy expectation of any normalized state in the one-dimensional harmonic oscillator form domain below by $E_0$, proving that the terminated state really is a ground state.

To find its probability density, solve the first-order equationIts solutions are ${\psi }_0=C_0\exp [-m\omega x^2/(2\hslash )]$. With $\alpha =m\omega /\hslash >0$, the supplied Gaussian integral givesTaking $C_0=(\alpha /\pi {)}^{1/4}$ fixes an irrelevant global phase. ConsequentlyThe ground-state energy is positive even though the potential minimum is zero: it is the oscillator's zero-point energy.

The density operator is diagonal in the oscillator number basis, so $⟨a^2⟩=⟨a^{†2}⟩=0$. With $\bar{n}=⟨a^{†}a⟩$ and $a{a}^{†}=a^{†}a+1$,HenceEquality holds even at zero temperature, when each contribution is $\hslash \omega /4$. Classical equipartition alone would miss this zero-point energy.

On the unbounded line the momentum density of states in a reference length $L$ is $Ldk/(2\pi )$. Count both signs of $k$, each with zero-point energy $|k|/2$ and the same frequency regulator:The extensive divergent term agrees exactly with c. Subtract it before taking the cutoff away:The force tending to change the interval length isIts negative sign means attraction: separating the boundaries increases this negative vacuum-energy difference toward zero. The bulk subtraction must be made before interpreting the force, otherwise differentiating an unphysical cutoff-dependent bulk energy would contaminate it.

This is the one-dimensional scalar Dirichlet Casimir force. The source's name “Casimir-Polder” normally refers to a retarded atom–atom or atom–surface dispersion force, whose geometry and power law differ. No atom polarizability or electromagnetic mode multiplicity enters the scalar result derived here.

For one oscillator, set $x=\hslash \Omega /(k_{B}T)$. Its quantum-oscillator energy and heat capacity areAt low temperature, $U\to \hslash \Omega /2$ and $C_V\sim k_B{x}^2{e}^{-x}\to 0$: excitations freeze out but zero-point energy remains. At high temperature, $U\sim k_{B}T$ and $C_V\to k_B$, recovering equipartition. The zero-point term does not contribute to heat capacity.

Ensemble and applications. A grand canonical ensemble describes a system exchanging both energy and particles with a large reservoir. Its fixed control parameters are $T,V,\mu$, while microstates with different $N$ have probabilitiesWeak coupling allows the system's energy levels to remain meaningful; equilibrium fixes the reservoir temperature and chemical potential. The grand potential is $\Phi =-k_{B}T\ln \Xi =U-TS-\mu N$, with $N$ now the mean population. It generates entropy, pressure and mean particle number by differentiation. This ensemble is convenient for gases in contact with particle reservoirs, independently occupied quantum levels and surface adsorption; its particle-number fluctuations are part of the description, not an error in counting.

Classical single-level limit. For the single-level grand partition function in the dilute limit, put $x_k=e^{-\beta ({ϵ}_k-\mu )}$. Empty and singly occupied levels contribute $1$ and $x_k$. In the dilute classical limit $x_k\ll 1$, higher occupancies contribute only at order $x_k^2$ or beyond:This truncation does not impose the Pauli principle on a classical gas. A fermion state has exactly $1+x_k$, a boson state has $(1-x_k{)}^{-1}$, and Maxwell-Boltzmann occupancy weights give $e^{x_k}$; they all agree to this order. In particular, $\ln {\Xi }_k=x_k+O(x_k^2)$, with the classical Maxwell-Boltzmann result exactly $x_k$.

grand potential and chemical potential. With $z=e^{\beta \mu }$, summing the single-level logarithms and using the stated density of states givesSubstitute $ϵ=k_{B}T{x}^2$, for whichCombining this with the density-of-states prefactor gives $\ln \Xi =zV/{\lambda }^3$. ThusDifferentiating at fixed $T,V$,Similarly $p=-({\partial }_V\Phi {)}_{T,\mu }=k_{B}Tz/{\lambda }^3$, so $pV=N{k}_{B}T$. This is the grand-canonical chemical potential of a classical ideal gas. The condition $N{\lambda }^3/V\ll 1$ states the classical regime. The supplied density of states counts one internal state; extra internal degeneracies would have to be included consistently rather than inserted at only one stage.