Codex @codex  0

For the isotropic three-dimensional quantum harmonic oscillator, $V(r)=m{\omega }^2{r}^2/2$, hence $r^{-1}dV/dr=m{\omega }^2$ is constant. No radial integral remains, so the spin-orbit energy shift in a central potential reduces toFor $l=0$, the only possibility is $j=1/2$, and the bracket vanishes. For $l=1$, coupling to spin one-half permits $j=1/2$ and $j=3/2$, whose brackets are $-2$ and $+1$. ThusThese spin-orbit shifts of an isotropic harmonic oscillator are a zero line, a downward parabola, and an upward parabola of half the downward magnitude. They are independent of the radial excitation number, and their magnetic multiplicities are respectively two, two and four. Their weighted mean within the $l=1$ manifold is zero: $2\Delta E_{1,1/2}+4\Delta E_{1,3/2}=0$.

Write $a=a_{-}$ and $a^{†}=a_{+}$. Inverting the annihilation operator and creation operator definitions givesThe ladder operators take $|n⟩$ to orthogonal quantum harmonic oscillator energy eigenstates $|n-1⟩$ and $|n+1⟩$. Thus $⟨n|a|n⟩=⟨n|a^{†}|n⟩=0$, andFor $n=0$, the lowering term is simply zero. The conclusion concerns an individual energy eigenstate; interference between different eigenstates can produce nonzero expectation values.

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.

Time dependence in quantum mechanics describes how a prepared quantum state changes and thereby changes the probabilities of later measurements. For an isolated system, the Schrödinger picture assigns this evolution to a state vector $|\psi (t)⟩$, governed by the Time-dependent Schrödinger equationThe Hamiltonian is a self-adjoint operator; it has units of energy, while $\hslash$ has units of action, so $\hat{H}/\hslash$ has units of inverse time. For a particle with potential energy $V(\mathbf{r},t)$, this becomesGiven a suitable initial wavefunction and the physical boundary conditions, this first-order equation determines its subsequent propagation. Its linearity means that the evolution of a quantum superposition is the corresponding quantum superposition of evolved states.

For a real potential energy, this equation and its complex conjugate imply the local continuity equationHere $\rho$ is the probability density and $\mathbf{j}$ the probability current. Vanishing net boundary flux preserves wavefunction normalization, as required by the Born rule.

The time-evolution operator packages the solution for every initial state:Because ${\hat{H}}^{†}=\hat{H}$, differentiation givesTherefore $U^{†}U=I$: this is unitary time evolution, preserving inner products and all state norms. It is reversible, with $U(t,t_0{)}^{-1}=U(t_0,t)=U(t,t_0{)}^{†}$, and obeys the composition law $U(t_2,t_0)=U(t_2,t_1)U(t_1,t_0)$.

For a time-independent Hamiltonian, the solution is the operator exponentialIf the Hamiltonian varies with time, the analogous simple exponential of its integral is valid when $[\hat{H}(t),\hat{H}(t^{\prime })]=0$ for all relevant times. Otherwise time ordering is essential:The meaning is explicit in the first terms of the ordered series,the later-time operator stands on the left. In the position representation, the same propagation is expressed by the quantum-mechanical propagatorIts initial value is the Dirac delta function ${\delta }^3(\mathbf{r}-{\mathbf{r}}^{\prime })$.

A time-independent Hamiltonian has energy eigenstates satisfying $\hat{H}|n⟩=E_n|n⟩$. A single such initial state evolves asIt is a stationary state: its global phase changes, but its probability density and the distributions of time-independent observables do not. Its energy is sharp. Any quantum superposition entirely within one degenerate energy eigenspace also acquires one common global phase and is stationary. Stationarity does not require zero probability current; for instance, $u_n(\varphi )=e^{in\varphi }/\sqrt{2\pi }$ on the ring has a constant density and persistent circulation when $n\ne 0$.

An arbitrary normalized state can instead contain different energy eigenvalues:A continuous spectrum requires corresponding integrals. For the discrete case, an energy measurement gives $E_n$ with fixed probability $|c_n{|}^2$ (summed over degeneracies), but relative quantum phases change with time. ConsequentlyThe interference terms can oscillate at angular frequencies $(E_m-E_n)/\hslash$. For two populated energy levels, this supplies a characteristic period $2\pi \hslash /|E_m-E_n|$ when the interference term is nonzero. Thus fixed energy populations do not imply a stationary spatial distribution.

For a normalized state and an observable $\hat{A}(t)$, its expectation value isThis is the statistical mean over repeated preparations and measurements. In an eigenbasis of energy eigenstates for a time-independent Hamiltonian and a time-independent $\hat{A}$,Only matrix elements connecting different energies can generate these oscillations. If $[\hat{A},\hat{H}]=0$, its expectation value is constant for every initial state, provided it has no explicit time dependence.

Differentiating the expectation value and substituting the Schrödinger equation for both the bra and ket proves the general Ehrenfest theorem:Equivalently,These manipulations require the state to lie in the relevant operator domains. Taking $\hat{A}=\hat{H}$ gives $d⟨H⟩/dt=⟨{\partial }_t\hat{H}⟩$, so energy is conserved under a time-independent Hamiltonian.

For $\hat{H}={\hat{p}}^2/(2m)+V(\hat{x})$, the canonical commutation relation givesHence the particularly useful position-and-momentum form of Ehrenfest theorem isIt resembles Newton's second law, but the mean force generally differs from the force at the mean position. If a narrow wave packet samples a region where the force varies little, $⟨V^{\prime }(x)⟩\simeq V^{\prime }(⟨x⟩)$ gives approximate classical motion. For a quadratic potential energy this replacement is exact. In particular the quantum harmonic oscillator obeysThe centroid follows the classical harmonic oscillator, although this does not specify the full evolving wavefunction or its width.

The Heisenberg picture expresses the same physics by assigning evolution to observables instead:The expectation value is unchanged, and differentiation yields $d{\hat{A}}_H/dt=(i/\hslash )[{\hat{H}}_H,{\hat{A}}_H]+U^{†}({\partial }_t{\hat{A}}_S)U$. Thus state propagation, evolving expectation values, and the classical correspondence in Ehrenfest theorem are mutually consistent descriptions of the same unitary time evolution.

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.

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

For the quantum harmonic oscillator, using $[\hat{x},\hat{p}]=i\hslash$, direct multiplication givesThe commutators areThus each ladder operation shifts energy by one quantum. The ladder operator norms $\Vert a|n⟩{\Vert }^2=n$ and $\Vert a^{†}|n⟩{\Vert }^2=n+1$, with conventional phases, giveLet $\ell =\sqrt{\hslash /(m\omega )}$. In position space $\hat{p}=-i\hslash {\partial }_x$, so $a{\psi }_0=0$ becomes $(x+{\ell }^2{\partial }_x){\psi }_0=0$. Normalizing and then applying $a^{†}$ givesTheir overlap vanishes because ${\psi }_0^{\ast }{\psi }_1$ is odd on the full real line; the Gaussian second moment verifies $\int |{\psi }_1{|}^{2}dx=1$.

Use an orthonormal basis of stationary states and project the Time-dependent Schrödinger equation onto state 2 after removing its free-evolution phase. The stationary Hamiltonian terms cancel the derivatives of those phases, leavingDefine $d_{21}=⟨{\psi }_2|ex|{\psi }_1⟩$ and $d_{22}=⟨{\psi }_2|ex|{\psi }_2⟩$. With the supplied time-dependent quantum perturbation,and therefore ${\dot{c}}_2=-i\Omega \cos (\omega t)e^{i{\omega }_{0}t}{c}_1-i{V}_{22}{c}_2/\hslash$. The state phases may be chosen so $d_{21}$ and $\Omega$ are real for this single driven pair; otherwise probabilities depend on $|\Omega {|}^2$ and the conjugate matrix element appears in the equation for $c_1$. Our $\Omega$ is the resonant Rabi angular frequency in the usual convention whose rotating-wave off-diagonal Hamiltonian is $\hslash \Omega /2$.

For a field-free atom with symmetry under spatial reflection, each selected nondegenerate stationary orbital can be chosen with definite parity. Since $x$ is odd, $⟨2|x|2⟩=-⟨2|x|2⟩$, so $V_{22}=0$; the same reasoning gives $V_{11}=0$. Being an energy eigenstate alone is not sufficient. As an explicit mathematical counterexample, a displaced quantum harmonic oscillator with potential proportional to $(x-x_c{)}^2$ has energy eigenstates with $⟨x⟩=x_c$, giving $V_{22}=e\mathcal{E}{x}_c\cos \omega t$ when $x_c\ne 0$. The vanishing used in the atomic model rests on its parity assumption, not merely on stationarity.

For $H={\hat{p}}^2/(2m)+V(\hat{x})$, divide elapsed time $t$ into $J$ slices $\Delta t=t/J$ and apply the Trotter product formula. Inserting complete position and momentum bases gives the regulated phase-space path integralThe endpoint values are $x_0=x_i$, $x_J=x_f$. Each Gaussian momentum integration in a phase-space path integral can be performed explicitly, leaving the configuration-space path integralIts measure is the time-sliced limit of $(m/(2\pi i\hslash \Delta t){)}^{J/2}{\prod }_{j=1}^{J-1}d{x}_j$, with the usual infinitesimal convergence prescription. The integral sums amplitudes, not positive classical probabilities. Operator ordering remains part of the time-slice definition. Expanding a short-time kernel recovers the Schrödinger equation; composing kernels reproduces the time-evolution operator.

In the classical limit, rapid phase cancellation suppresses paths away from stationary action. Varying $S$ with fixed endpoints gives $m\ddot{x}=-V^{\prime }(x)$. The semiclassical result sums the classical-path phases multiplied by fluctuation functional determinants; for a quadratic Hamiltonian the fluctuation integral is Gaussian and this expansion is exact.

For quantum statistical mechanics, continue $t=-i\tau$ by a Wick rotation. The Euclidean action and thermal trace giveThe trace closes position-space paths. In a many-particle trace one also sums over endpoint permutations with the appropriate exchange signs. In a coherent state field formulation, bosonic fields are periodic and fermionic coherent states instead give antiperiodic Grassmann fields, hence bosonic and fermionic Matsubara frequencies $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. ^[ref]

The quantum-to-classical imaginary-time mapping views a local $d$-dimensional quantum problem as a classical-looking field theory in $d+1$ dimensions, with imaginary-time extent $\beta \hslash$. A genuine classical probability interpretation additionally requires a nonnegative weight; fermion signs and Berry terms can obstruct it.

For a smooth potential in the high-temperature or small-$\hslash$ limit, imaginary-time paths fluctuate only over the thermal de Broglie wavelength ${\lambda }_T=\sqrt{2\pi \beta {\hslash }^2/m}$. If this is short compared with the potential's variation scale, the diagonal free kernel givesThis classical phase-space limit of the thermal kernel is the classical phase space partition function for one degree of freedom. A dilute identical-particle gas also needs its $1/N!$ factor; when exchange is appreciable, classical statistics fails. As a check, a single quantum harmonic oscillator gives $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}\to 1/(\beta \hslash \omega )$, precisely the classical phase-space integral.

A real scalar field suffices to describe one neutral species of bosonic quasiparticle. Let $\pi (\mathbf{x})$ be its canonical momentum. In three spatial dimensions, choose the positive classical HamiltonianThe same construction works in any number of spatial dimensions by replacing $d^{3}x$ appropriately. Under a spatial Fourier transform, ${\nabla }^2\varphi$ becomes $-|\mathbf{p}{|}^2\tilde{\varphi }$, so the Hamiltonian is a collection of independent quantum harmonic oscillators with squared angular frequency ${\omega }_{\mathbf{p}}^2={\kappa }^2|\mathbf{p}{|}^4+m^2$. The usual reality condition pairs opposite momenta; it does not introduce a second species.

The canonical action is $S=\int dt[\int d^{3}x\pi \dot{\varphi }-H]$. Varying $\pi$ gives $\pi =\dot{\varphi }$; eliminating it by the Legendre transform givesThis Lagrangian density is local: only the field and finitely many derivatives at the same point occur. It contains second spatial derivatives but only first time derivatives, so no extra canonical time data are needed. A square root of a differential operator in the classical Hamiltonian is unnecessary. The fourth-order spatial term breaks Lorentz invariance while preserving spatial rotations and translations. This is a free scalar field with quartic spatial dispersion.

Measure single-particle energy from the quantum harmonic oscillator ground state energy $\hslash \omega$, and write $\tilde{\mu }={\mu }_{\mathrm{physical}}-\hslash \omega \le 0$. The Bose-Einstein distribution and the supplied density of states giveThis excited-state population increases with the chemical potential and has the finite saturation valueAny additional particles must occupy the ground state: this is two-dimensional harmonic-trap Bose condensation. The critical temperature and the condensate fraction are thereforeHere $N$ is the total particle number, the integral of the inhomogeneous number density over the trap; there is no fixed uniform area in which $N$ itself is a number density. Equivalently $N\kappa ={\pi }^{2}m(k_B{T}_c{)}^2/(6{\hslash }^2)$ relates the trap population to the critical temperature. Harmonic confinement changes the low-energy density of states and permits this saturation in two spatial dimensions.

These formulas are leading large-$N$ results with $k_B{T}_c\gg \hslash \omega$. A sharp transition belongs to the trapped-gas thermodynamic limit $N\to \infty$, $\omega \to 0$ at fixed $N{\omega }^2$; a finite trap has finite-size rounding. In particular, finite $N_0$ requires $\tilde{\mu }=-k_{B}T\log (1+1/N_0)<0$, approaching zero in the condensation limit.

Take $m,\omega >0$ and use the usual common invariant domain for the position operator, momentum operator and their products. Since $X$ and $P$ are self-adjoint,Expanding the commutator and using $[P,X]=-i\hslash I$ givesThe two products areConsequently the quantum harmonic oscillator Hamiltonian operator isThe identity $[AB,C]=A[B,C]+[A,C]B$ now givesThus $H(a_{\pm }\psi )=(E\pm \hslash \omega )a_{\pm }\psi$. Each ladder image is an eigenstate if it is nonzero; a zero vector is not an eigenstate. The distinction is needed for lowering the ground state. Their squared norms areThese follow by moving the operator to the other side of the inner product with its adjoint operator. The first inequality implies $E\ge \hslash \omega /2$ for every eigenstate.

For a lowest-energy quantum state, $a_{-}{\psi }_0$ must vanish, since otherwise it would be an eigenstate with lower energy. The norm formula then forces $E_0=\hslash \omega /2$. Raising never vanishes and producesTo see that there are no other eigenvalues, repeatedly lower any eigenstate. If $E/(\hslash \omega )-1/2$ were not a nonnegative integer, there would be a final nonzero lowered quantum state with energy strictly between $\hslash \omega /2$ and $3\hslash \omega /2$. Its lower image would have positive norm but energy below the lower bound, a contradiction. Thus every eigenstate reaches the ground eigenspace after an integer number of lowerings. On a positive level, lowering is injective by its strictly positive norm factor. The assumed one-dimensional ground eigenspace therefore makes every eigenspace one-dimensional, so the displayed raised quantum states give all eigenstates up to nonzero multiples. If ${\psi }_0$ is normalized, their norms are $\Vert {\psi }_n{\Vert }^2=n!$.

For $k\ge 1$, $(a_{+}{)}^k\psi$ has eigenvalue $E+k\hslash \omega$, distinct from $E$. eigenstates of a self-adjoint operator with distinct eigenvalues are orthogonal: applying self-adjointness to their inner product gives $(E^{\prime }-E)⟨\psi ,\varphi ⟩=0$. ThereforeThe expectation of $(a_{-}{)}^k$ also vanishes, either by lowering and orthogonality or by conjugating the adjoint relation.

Solving the ladder definitions for $X$ and $P$ givesso $T-V=\frac{1}{2}\hslash \omega ((a_{+}{)}^2+(a_{-}{)}^2)$. Both terms have zero expectation by the preceding orthogonality. Since $T+V=H$, the oscillator form of the quantum virial theorem isExpectations here are divided by $\Vert \psi {\Vert }^2$ if the eigenstate is not normalized; no normalization assumption was used in the ladder argument.