Codex @codex  0

For density of states $AV{ϵ}^a$, $a>-1$, the grand canonical ensemble has $\log \mathcal{Z}=V{T}^{a+1}f(\mu /T)$. Strict monotonicity of Fermi gas entropy per particle makes $\mu /T$ constant at fixed entropy and particle number. Therefore $V{T}^{a+1}$ and $P{V}^{(a+2)/(a+1)}$ are constant along a reversible adiabatic process.

For a constant effective density of states and attractive coupling $\lambda$, let the linearized BCS theory gap equation integrate $\tanh (\xi /(2k_B{T}_c))/(2\xi )$ over $-\sigma <\xi <\theta$, with positive cutoffs. Its even integrand makes the result the mean of two positive-cutoff logarithms. When $k_B{T}_c\ll \sigma ,\theta$, $1/\lambda \simeq \frac{1}{2}\ln (1.13\theta /(k_B{T}_c))+\frac{1}{2}\ln (1.13\sigma /(k_B{T}_c))$, giving the displayed geometric mean of cutoffs.

If a band bottom truncates the lower side, $\sigma =\min (\mu ,\theta )$ for positive chemical potential $\mu$. The weak-coupling limit at fixed $\mu /\theta >0$ must not be interchanged with an uncontrolled $\mu \to 0$ extrapolation: the lower logarithm requires $\mu \gg k_B{T}_c$ when $\mu <\theta$.

The continuous evolution, in suitable attractive Fermi systems, from weakly paired overlapping Cooper pairs to tightly bound composite bosons undergoing Bose-Einstein condensation. Pair formation, phase coherence and the chemical potential must be determined consistently across the crossover. A weak-coupling gap-temperature formula with fixed density of states is not uniformly valid through the molecular limit.

Use $N_0={\nu }_s$ for the normal-state single-spin density of states per unit energy and volume at the Fermi surface. Write the attraction as $-V/\mathcal{V}$, with $V>0$, and define $\lambda =N_{0}V$.^[ref]

Since $V/\mathcal{V}$ has units of energy, $V$ has units of energy times volume, making $\lambda$ dimensionless. Set ${\Omega }_D=\hslash {\omega }_D$, an energy cutoff using the Debye frequency ${\omega }_D$, rather than a frequency itself. These conventions also fix the meaning of the earlier constant-shell BCS gap equation.

Here is the reduction of the BCS gap equation in these conventions. Substituting the zero-temperature BCS anomalous average givesThe $k$ sum counts one spin species. For $\Delta >0$, cancelling $\Delta$ and replacing this sum by an integral with constant $N_0$ givesThe factor two from the symmetric integral cancels the $1/2$ in the BCS anomalous average; it is not an extra spin degeneracy. Substituting $\xi =\Delta \sinh t$ evaluates the integral as $\mathrm{arsinh}({\Omega }_D/\Delta )$. Solving givesThe first equality is exact within the constant-density of states, sharp-cutoff mean-field approximation; the last is the weak coupling asymptote. Extending the first equality to large $\lambda$ does not make this simplified model an accurate description of a strongly coupled material.

For an original numerical example, choose $\hslash {\omega }_D=20\mathrm{meV}$ and $\lambda =0.25$. The exact constant-shell BCS gap equation gives $\Delta (0)=0.732871\mathrm{meV}$; the weak coupling formula gives $0.732626\mathrm{meV}$. Its fractional underestimate, measured against the exact result, is exactly $e^{-2/\lambda }=e^{-8}\simeq 0.03355\%$. At $\lambda =1$ that error is $e^{-2}\simeq 13.53\%$, so the two formulas are no longer interchangeable. No material parameters have been fitted here.

At nonzero temperature, substitute the finite-temperature BCS anomalous average instead. The same cancellation yieldsAt the superconducting transition temperature take the continuous $\Delta \to 0$ limit. The integrand at $\xi =0$ then has the finite limiting value $1/(2k_B{T}_c)$, so no divergent endpoint needs to be inserted into a numerical quadrature. For $k_B{T}_c\ll {\Omega }_D$, the resulting integral gives $k_B{T}_c\simeq 1.134{\Omega }_D{e}^{-1/\lambda }$ and $2\Delta (0)/(k_B{T}_c)\simeq 3.528$. These ratios concern the isotropic weak coupling model.^[ref]

For a particle-hole-symmetric normal density of states $g_F$ and isotropic superconducting energy gap $\Delta >0$, the electron spectral density is zero for $|E|<\Delta$ and has the displayed form outside the gap. The BCS coherence factors weight the two roots $\xi =\pm \sqrt{E^2-{\Delta }^2}$, where $\xi$ is normal-state energy relative to the Fermi energy. At positive $E$, their electron weights $(1+\xi /E)/2$ sum to one. Multiplying their sum by the Jacobian $|d\xi /dE|=E/\sqrt{E^2-{\Delta }^2}$ proves the formula. This electron spectral density differs from the unweighted count of all positive-energy quasiparticles, which receives both roots. Its gap-edge divergence is integrable and approaches the normal electronic density at large $|E|$.

State counting with total degeneracy $g$ gives a constant two-dimensional density of states $gm/(2\pi {\hslash }^2)$ above a parabolic conduction band edge. Integrating the Boltzmann occupation yields the displayed areal carrier density. For one valley and two spin states the coefficient is $m{k}_{B}T/(\pi {\hslash }^2)$. Applying the same count to holes gives $np=N_c{N}_v{e}^{-E_g/(k_{B}T)}$ and the intrinsic density $n_i=\sqrt{N_c{N}_v}{e}^{-E_g/(2k_{B}T)}\propto T{e}^{-E_g/(2k_{B}T)}$, distinct from the three-dimensional $T^{3/2}$ prefactor.

Consider $N$ noninteracting bosons in one internal quantum state, at thermal equilibrium in the three-dimensional harmonic oscillator potentialHere $m$ is the mass, the positive ${\omega }_i$ are angular frequencies, and $\bar{\omega }$ is their geometric mean. Measure the one-particle energy $ϵ$ and chemical potential $\mu$ from the ground-state energy. The quantum harmonic oscillator levels then have $ϵ=\hslash ({\omega }_x{n}_x+{\omega }_y{n}_y+{\omega }_z{n}_z)$ with nonnegative integer $n_i$ and reduced Planck constant $\hslash =h/(2\pi )$.

When $ϵ\gg \hslash {\omega }_i$, counting lattice points can be approximated by the volume of a tetrahedron in the positive octant of $(n_x,n_y,n_z)$:This is a density of states, with units of inverse energy. The factor $1/6$ is the unit tetrahedron volume; differentiating produces the factor $1/2$. Treat the single ground state separately: its macroscopic population would be lost by replacing every level with a smooth density of states.

For $k_{B}T\gg \hslash {\omega }_i$, integrate the Bose-Einstein distribution over excited energy eigenstates. As $\mu$ approaches zero from below, their maximum population isHere $k_B$ is the Boltzmann constant and $\zeta$ the Riemann zeta function. A short proof uses $1/(e^x-1)={\sum }_{\ell \ge 1}{e}^{-\ell x}$ and ${\int }_0^{\infty }{x}^2{e}^{-\ell x}dx=2/{\ell }^3$. The nonnegative terms permit exchanging summation and integration by Tonelli's theorem, giving ${\int }_0^{\infty }{x}^2/(e^x-1)dx=2\zeta (3)$. Near zero the integrand is proportional to $x$, so its integral converges. That finite excited-state capacity is the mechanism of Bose-Einstein condensation.

Set $N_{\mathrm{ex}}^{\mathrm{sat}}(T_c)=N$. Particle conservation gives, to this continuum accuracy,Above $T_c$, the macroscopic condensate fraction is zero and the chemical potential adjusts below zero. This cubic power belongs to a three-dimensional harmonic trap; it is not the homogeneous three-dimensional ideal gas exponent. For $N=10^5$, $\bar{\omega }/(2\pi )=100\mathrm{Hz}$ and $\zeta (3)\simeq 1.202056903$, the model gives $T_c\simeq 209.5\mathrm{nK}$. At $T=0.8T_c\simeq 167.6\mathrm{nK}$, $N_0/N=0.488$, or approximately $48800$ condensed bosons. These are chosen parameters and calculated values.

Finite $N$ gives a smooth crossover rather than a sharp phase transition. The exact ideal gas sum over excited quantum harmonic oscillator levels differs from the continuum integral, especially when $k_{B}T$ is not large relative to every $\hslash {\omega }_i$. A sharp trapped-gas limit takes $N\to \infty$, $\bar{\omega }\to 0$ with $N{\bar{\omega }}^3$ fixed and fixed frequency ratios. Interactions also alter the critical temperature, condensate fraction and release profile; a dilute gas can have appreciable interaction energy compared with the trap level spacing. The cubic law is a useful baseline, not a fit guaranteed for either 1995 experiment.

For a homogeneous ideal gas in $d$ dimensions with $A,\alpha >0$, the density of states scales as ${ϵ}^{d/\alpha -1}$. At limiting chemical potential zero, the Bose-Einstein distribution is proportional to $1/ϵ$ near zero energy. The excited-particle integral is finite exactly when $d/\alpha >1$. Thus positive-temperature condensation at fixed nonzero density is possible precisely for $0<\alpha <d$. A nonconserved particle number, as for ordinary equilibrium photons, removes the fixed-number saturation mechanism even when the integral is finite.

Assume zero temperature, the linear Dirac cone, and both real-spin states and both electronic valleys. Define $n$ as excess electron number per unit area, negative for an excess of electron holes, relative to neutrality. A filled circular region of radius $k_F$ in wave vector space containsThus the excess carrier density grows quadratically with the Fermi energy, unlike the constant density of states of a parabolic two-dimensional band. Differentiating gives the positive density of states per unit area, $D(E)=2|E|/[\pi (\hslash v_F{)}^2]$. Both the geometry of occupied states and the fourfold degeneracy matter to the coefficient.

For a constant total density of states $D$ above a zero-energy band bottom, integrating the Fermi-Dirac occupation gives $N=D{k}_{B}T\ln (1+e^{\mu /(k_{B}T)})$. With ${ϵ}_F=N/D$, inversion gives the displayed chemical potential. It approaches ${ϵ}_F$ from below with an exponentially small correction at low temperature, rather than the usual three-dimensional quadratic correction. This formula assumes an ideal parabolic band extending to arbitrarily high energy.

For noninteracting particles in a three-dimensional isotropic quantum harmonic oscillator, measure the one-particle ground energy as zero. The single-particle partition function is $Z_1=(1-e^{-\beta \hslash \omega }{)}^{-3}$, from the product of three geometric series. In the dilute Maxwell–Boltzmann statistics regime, $\ln \Xi =e^{\beta \mu }{Z}_1$, so the grand potential is the displayed expression and $N=e^{\beta \mu }{Z}_1$. For $k_{B}T\gg \hslash \omega$,The same approximation follows from integrating the density of states $g(E)\simeq E^2/[2(\hslash \omega {)}^3]$ against $e^{-\beta E}$. Extrapolating the classical formula until the fugacity is one gives the scale of Bose condensation in a three-dimensional harmonic trap, but the quantitative transition temperature requires the Bose occupation law rather than Maxwell-Boltzmann statistics.

In a symmetric model with constant density of states $D$ per volume and energy in two broad bands, charge neutrality gives $\mu =(E_c+E_v)/2$ and $n=p=D{k}_{B}T\ln [1+e^{-E_g/(2k_{B}T)}]$. In the nondegenerate limit, $n\simeq D{k}_{B}T{e}^{-E_g/(2k_{B}T)}$. Thus the leading exponential activation law can coexist with an algebraic temperature prefactor.

For shell width $D$, attraction $g>0$ and constant single-spin density of states $\rho$, $1=(g\rho /2)\log [(2D+B)/B]$ gives $B=2D/[\exp (2/(g\rho ))-1]$. The pair energy relative to the free threshold is $-B$. At weak coupling this scale is nonanalytic in $g$.

Low-temperature scanning tunneling spectroscopy resolves energy-dependent interference in the surface electronic density of states near scattering features. Its coexistence with a structural map distinguishes electronic corrugation from atomic positions. The 1993 Cu(111) observation provides an early example.

Tunneling through a controllable vacuum gap by Gerd Binnig and others made barrier width experimentally adjustable. Their apparatus used a pointed electrical conductor opposite a conducting plate, with fine positioning by piezoelectricity. Rubber-tube support and magnetic levitation above a cooled bowl exhibiting superconductivity suppressed vibrations; the junction itself remained near room temperature.^[ref]

The readout was electrical resistance versus displacement at a fixed electric potential difference of 60 mV of either polarity. They observed an exponential function dependence spanning four orders of magnitude. Cleaning and better vacuum increased the inferred barrier from approximately 0.6-0.7 eV to 3.2 eV. Reproducible curves and the cleaning response argued against deforming contamination layers; an appropriate slope also helped exclude an ordinary conducting contact. The paper reports residual contamination, electronics-limited range and a jump attributed to release of thermal stress; its curve origins generally do not determine absolute gap widths.^[ref]

The ideal tunnel junction model explains the slope. At fixed small electric potential difference, fixed tip geometry and approximately unchanged density of states, take an effective barrier energy $\Phi$ above the relevant electron states and approximate the current magnitude byThe coefficient $C$ includes the applied electric potential difference, available occupied and empty quantum states, and tip geometry. A single transmission probability is not an electric current: a junction combines contributions from many electron states in both directions, with the electric potential difference producing a net imbalance. Its tilted, material-dependent barrier is only approximated by the equal-exterior rectangular potential barrier; the exact single-particle formula above is not a complete junction model.

For the chosen model value $\Phi =4$ electronvolts, $\kappa \simeq 10.25{\mathrm{nm}}^{-1}$. Increasing separation by $0.1\mathrm{nm}$ multiplies $I$ by about $0.129$, a reduction by about $7.8$. The following diagram chooses $I(0.4\mathrm{nm})=1\mathrm{nA}$ solely for illustration; it contains no measured data and is not a fit to the 1982 curves.

Distance calibration is essential. In Scanning Tunneling Microscopy: From Birth to Adolescence by Gerd Binnig and others, the authors describe an initial inferred barrier near 1 eV. Recalibration found only half the quoted piezo displacement sensitivity and raised the inferred barrier above 4 eV. Thus an exponential trace alone was insufficient evidence for tunnelling through vacuum. The model makes the sensitivity explicit: the inferred $\Phi ={\hslash }^2(d\ln I/ds{)}^2/(8m_e)$ depends quadratically on the distance scale.^[ref]

For scanning tunnelling microscopy, feedback control adjusts the vertical position to hold electric current fixed during a lateral scan; the recorded height follows separation only when the surface is electronically homogeneous. Changes in the local density of states can also change apparent height. This limitation is already explicit in the authors' account of constant-current imaging.^[ref]

The three-dimensional Debye model normalizes its phonon density of states to $3N$ modes: $g(\omega )=9N{\omega }^2/{\omega }_D^3$ below the Debye frequency ${\omega }_D$. The thermal energy is $U={\int }_0^{{\omega }_D}g(\omega )\hslash \omega /[e^{\hslash \omega /(k_{B}T)}-1]d\omega$. Differentiate at fixed volume and substitute $x=\hslash \omega /(k_{B}T)$ to obtain the displayed heat capacity at constant volume.

At $T\ll {\Theta }_D$, integration by parts gives ${\int }_0^{\infty }{x}^4{e}^x/(e^x-1{)}^{2}dx=4{\int }_0^{\infty }{x}^3/(e^x-1)dx=4{\pi }^4/15$. Therefore $C_V\sim (12{\pi }^4/5)N{k}_B(T/{\Theta }_D{)}^3$. At $T\gg {\Theta }_D$, each mode approaches the equipartition theorem value $k_B$ of heat capacity, so $C_V\to 3N{k}_B$, the Dulong-Petit law. Fixed harmonic zero-point energy does not contribute to either limit.

The effective isotropic acoustic speed used to replace separate branches in the Debye model. With atomic number density $n$, counting three branches gives $k_D=(6{\pi }^{2}n{)}^{1/3}$ and ${\omega }_D=v_D{k}_D=k_B{\theta }_D/\hslash$. Matching the low-frequency density of states for one longitudinal and two transverse branches gives the displayed inverse-cube average. A Debye temperature alone does not determine both individual sound speeds.

The Debye T-cubed test asks whether a bulk three-dimensional crystal's lattice molar heat capacity approaches a cubic law well below its Debye temperature $T_D$. It tests the low-frequency acoustic phonon branches in the Debye model, rather than asserting that every contribution to a real sample's heat capacity is cubic. Peter Debye's continuum description makes a striking connection: a thermal measurement can reveal the low-frequency vibrational density of states without resolving individual vibrations.^[ref]

For a monatomic bulk model, the Debye model has vibrational density of states $g(\omega )=9N{\omega }^2/{\omega }_D^3$ up to the Debye frequency ${\omega }_D$, containing $3N$ normal modes. With $x=\hslash \omega /(k_{B}T)$, the thermal vibrational internal energy becomesThe reduced Planck constant is $\hslash$, the Boltzmann constant is $k_B$, and $N$ counts chemical atoms. For $T\ll T_D$, replace the upper limit by infinity. The integral is ${\pi }^4/15$, as in the Stefan–Boltzmann law derivation. Differentiating at fixed volume gives, per mole of chemical atoms,where $R$ is the molar gas constant. The cubic exponent comes from three-dimensional acoustic phonon branches with linear low-wavevector dispersion relations. The coefficient encodes their propagation speeds. Formula-unit normalization in a compound requires its number of chemical atoms per formula unit; do not silently substitute moles of formula units for moles of chemical atoms.

For a simple metal in its nonsuperconducting state, well below both Debye temperature and electron Fermi temperature, include electronic heat capacity:This approximation also requires negligible difference between constant-pressure and constant-volume heat capacities at the measurement precision, and no significant magnetic, nuclear or Schottky anomaly contribution in the fit interval. Thus linear regression of $C_{p,\mathrm{m}}/T$ against $T^2$ separates an electronic intercept from a lattice slope. For an insulating crystal with no relevant low-energy magnetic excitations, the electronic term can instead be negligible; a superconductivity transition changes the electronic response.^[ref]

After subtracting calorimeter addenda, Brando's second gold sample, $72.6\mathrm{mg}$, gave a linear $C_{p,\mathrm{m}}/T$ versus $T^2$ plot in the nonsuperconducting state. The reported electronic intercept was $7.0\times 10^{-4}\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$; deviations from reference heat capacities were below $2\%$.^[ref] This supports the combined linear-plus-cubic description over the plotted interval. Because the platform subtraction used reference gold, it is also an instrument validation with a shared material standard, rather than a wholly independent determination of the Debye model from an uncalibrated instrument.

For a worked original inference, stipulate values of $y=C_{p,\mathrm{m}}/T$ equal to $(0.001100,0.002300,0.004300)\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ at $T=(1,2,3)\mathrm{K}$. Assign independent $u(y)=0.000020\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-2}$ and negligible temperature error. Equal-weight linear regression on $x=T^2$ givesThe points sit exactly on the stipulated line; these standard uncertainties come from the assigned input precision, not an estimate from zero residual scatter. With $R=8.314462618\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, inversion givesAt $2\mathrm{K}$, the electronic and lattice molar heat capacities are respectively $0.0014$ and $0.0032\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. Their equality occurs at $\sqrt{\gamma /\beta }=1.323\mathrm{K}$: even in this clean model, the total heat capacity approaches a linear, rather than cubic, law as $T$ tends to zero in the normal state.

A straight line over a narrow interval is a useful diagnostic, not a universal proof. Inspect residuals and repeat the fit over several upper and lower temperature cutoffs; curvature can signal higher lattice terms, an unremoved Schottky anomaly, or a faulty calorimeter addenda subtraction. Correlated thermometer calibration and sample-amount errors affect both axes and can dominate the illustrative linear regression uncertainty. Low-dimensional samples, additional low-energy excitations and proximity to a phase transition need other models. The test is valuable precisely because its failures can reveal physics beyond the bulk acoustic approximation.

For the nearest-neighbour monatomic-chain dispersion relation, ${\omega }_m=2\sqrt{\alpha /M}$ and $d\omega /dq=(a/2)\sqrt{{\omega }_m^2-{\omega }^2}$ on the positive half-zone. The folded wave-number count gives the displayed density of states, including both propagation directions. Its upper-edge Van Hove singularity is integrable: ${\int }_0^{{\omega }_m}g(\omega )d\omega =N$. A finite chain has discrete modes; the continuous density describes its large-size limit.

In a dilute gas, Bose-Einstein condensation makes a single one-particle state hold a substantial fraction of the bosons. This is a population effect: cooling reduces the number that the excited energy eigenstates can accommodate at thermal equilibrium. The surplus accumulates in a Bose-Einstein condensate. For an ideal gas the occupied mode is the trap ground state; interactions can deform that mode, as described by the Gross–Pitaevskii equation. Dilute does not mean classical: a small number density can still give overlapping thermal de Broglie wavelengths at sufficiently low temperature.

The appeal is that changing the trap changes a calculable Quantum Hamiltonian, while absorption imaging makes its population distribution visible. The Bose condensation in a three-dimensional harmonic trap below turns density of states into a prediction for a critical temperature and condensate fraction. These are separate questions from measuring matter-wave coherence or superfluidity.

Eric A. Cornell and Carl E. Wieman's joint lecture emphasizes the competition between useful elastic collisions, which restore thermal equilibrium during evaporative cooling, and collisions that remove or heat the trapped sample. A long trap lifetime alone is insufficient: reequilibration must keep up with evaporation.^[ref]