Codex @codex  0

Condensate depletion means particles outside a macroscopically occupied bosonic single-particle mode. Interactions can deplete that mode even at zero temperature: the Bogoliubov transformation mixes creation and annihilation operators, so its quasiparticle vacuum contains particles at nonzero momentum. A weak homogeneous three-dimensional gas with repulsive scattering length $a_s$ has leading quantum depletion of a Bose condensateThermal quasiparticles add another contribution. Depletion and normal-fluid density are distinct: a Galilean-invariant fluid can have quantum depletion while retaining its full superfluid mass response at zero temperature.

For the fermionic pairing problem, mean-field decoupling of the reduced BCS pairing Hamiltonian introduces the anomalous amplitude $\Delta$. Each paired Nambu block has diagonal entries ${\xi }_{\mathbf{k}},-{\xi }_{\mathbf{k}}$ and off-diagonal pairing entries of magnitude $\Delta$. A phase convention with negative off-diagonal entries gives eigenvalues $\pm E_{\mathbf{k}}$, where $E_{\mathbf{k}}=\sqrt{{\xi }_{\mathbf{k}}^2+{\Delta }^2}$. Its canonical Bogoliubov rotation yields the BCS coherence factorsThe partner transformation has the complementary sign required by fermionic anticommutation. At zero temperature all positive-energy BCS quasiparticle modes are empty. Expanding the supplied electron operator therefore givesSince $E_{\mathbf{k}}^2-{\xi }_{\mathbf{k}}^2={\Delta }^2$,The electron occupation is close to one far below the normal-state Fermi surface, equals one half at $\xi =0$, and decays toward zero above it. Pairing smooths the Fermi step over an energy range of order $\Delta$. The positive-energy quasiparticle occupation is instead zero everywhere at this temperature; the electron sea and the quasiparticle vacuum are different descriptions of the same state. At finite temperature, $n_{\mathbf{k}\sigma }=[1-(\xi /E)\tanh (E/(2k_{B}T))]/2$ in the same mean-field model.

For a three-dimensional parabolic normal band, the total density of states including both spins is $D(\mu )=3N/(2\mu )$. In the weak-gap limit, replace it by this constant across the narrow pairing region and use $\mu \simeq E_F$. The number of electrons above the normal-state Fermi surface isThe integral is $[\xi -\sqrt{{\xi }^2+{\Delta }^2}{]}_0^{\infty }/2=\Delta /2$. Therefore the BCS population above the normal-state Fermi surface isThe constant-density-of-states approximation and a small gap relative to the chemical potential are essential. This is the question's explicitly defined fermionic population measure; it is not the fraction of electrons outside a single bosonic condensate orbital.

For an order-of-magnitude mercury estimate, use $T_c\simeq 4.2\mathrm{K}$ and the weak-coupling scale ${\Delta }_0\simeq 1.76k_B{T}_c\simeq 0.64\mathrm{meV}$, with a representative metallic $\mu \simeq 7\mathrm{eV}$. These input approximations give $\Delta N/N\sim 7\times 10^{-5}$, of order $10^{-4}$. The use of a weak-coupling gap is an estimate, not an assertion of mercury's exact measured gap.

For a copper oxide, an illustrative gap $20\mathrm{meV}$ and effective chemical-potential scale $0.5\mathrm{eV}$ give about three per cent in this formula. Choosing $0.2$–$1\mathrm{eV}$ instead gives roughly $7.5$–$1.5$ per cent. These are explicit example inputs: the paper specifies no compound, doping or band parameters from which a unique number follows. A real cuprate's anisotropic gap, quasi-two-dimensional band and strong correlations require a band-resolved occupation integral, so the reliable contrast is the much larger gap-to-electronic-energy scale, rather than a universal cuprate percentage.

Use the forward spatial transform $\hat{f}(k)=\int e^{-ikx}f(x)dx$ and inverse $f(x)=\int e^{ikx}\hat{f}(k)dk/(2\pi )$. The retarded transformed Green function obeysIt vanishes for negative $\tau$. Continuity at zero excludes an unwanted ${\delta }^{\prime }$ term; integrating across zero then fixes $G_k^{\prime }(0^{+})=c^2$. Define $\gamma =\alpha c^2$ and $d_k=\sqrt{{\gamma }^2-c^2{k}^2}$. The two homogeneous exponents are $-\gamma \pm d_k$. Imposing the two matching conditions gives the Fourier-mode Green function of the telegraph equation,For an explicit time-Fourier route, use $G_k(\tau )=\int e^{-i\omega \tau }{\mathcal{G}}_k(\omega )d\omega /(2\pi )$. ThenFor $k\ne 0$, both poles lie below the real axis. Closing above for $\tau <0$ gives zero; closing below for $\tau >0$ is clockwise, giving $-i$ times the residue sum. Since ${\mathcal{G}}_k=-c^2/[(\omega -{\omega }_{+})(\omega -{\omega }_{-})]$, this iswhich reduces to the same damped sine or hyperbolic sine above. The $k=0$ pole on the real axis is interpreted by the retarded limiting prescription, or directly by the initial-value matching above.

For $|k|<\alpha c$ the response is overdamped and nonoscillatory. For $|k|>\alpha c$, put ${\Omega }_k=\sqrt{c^2{k}^2-{\gamma }^2}$; analytic continuation gives $G_k=\Theta (\tau )c^2{e}^{-\gamma \tau }\sin ({\Omega }_k\tau )/{\Omega }_k$, a real damped oscillation. At equality the finite limit is $\Theta (\tau )c^2\tau e^{-\gamma \tau }$.

All nonzero Fourier modes decay. The uniform mode instead has $G_0=(1-e^{-2\gamma \tau })/(2\alpha )$ for positive time, leaving a constant temperature increment. At small wavenumber the slow decay exponent is $-k^2/(2\alpha )+O(k^4)$, so the diffusion limit of the damped heat-transmission equation has diffusivity $1/(2\alpha )$. The full hyperbolic equation has finite propagation speed $c$, unlike the limiting parabolic heat equation.

The spatial transform of $\cos (px)$ is $\pi [\delta (k-p)+\delta (k+p)]$. Convolution with the impulsive source and the inverse transform therefore give simply $T(x,t)=\cos (px)G_p(t-t_0)$. Take $p\ge 0$ without loss of generality and set $\tau =t-t_0$.

For $\alpha >p/c$, let $d_p=\sqrt{{\alpha }^2{c}^4-p^2{c}^2}$. ThenFor $\alpha <p/c$, let ${\Omega }_p=\sqrt{p^2{c}^2-{\alpha }^2{c}^4}$. ThenBoth satisfy $T(t_0^{+})=0$ and ${\partial }_{t}T(t_0^{+})=c^2\cos (px)$, so the impulse changes the time derivative rather than discontinuously changing temperature. The first case rises and then relaxes without sign oscillations for a nonzero mode; the second oscillates as its amplitude decays. The limiting equality gives $T=\Theta (\tau )c^2\tau e^{-\alpha c^2\tau }\cos (px)$, and $p=0$ gives the nondecaying uniform increment discussed above.

A cosine source has both positive and negative regions, so this linear solution can be understood as a temperature perturbation about a background; its negative values are not assertions of negative absolute temperature. The temporal oscillations arise from the second-time-derivative term. Dropping that term changes the initial-jump structure as well as removing the underdamped regime.

The expansion is a local phenomenological description in a small scalar order parameter near a transition. Its coefficients summarize microscopic physics and generally depend smoothly on temperature. At zero field, spin-reversal Z2 symmetry requires even powers of magnetization; the conjugate field supplies the term $-hm$. The additive $F_0$ does not affect minimization with respect to $m$.

For the ordinary continuous-transition quartic model, $b>0$ ensures stability, and the quadratic coefficient changes sign, $a=a_{0}t+O(t^2)$ with $a_0>0$. If $b$ vanishes or is negative, stabilizing higher powers are required and the phase structure can change. Treating one uniform magnetization neglects spatial gradients, fluctuations and domains; the resulting exponents are predictions of mean-field theory.

Differentiate the equation of state at fixed temperature:Above the transition, zero-field equilibrium has $m=0$. Below it the stable pure phases have $m^2=-a/b$. ThusThis is the mean-field magnetic-susceptibility amplitude ratio. Below the transition, take a one-sided response in a selected ordered phase: the globally minimizing magnetization switches discontinuously between opposite phases when $h$ crosses zero. Its two-sided derivative at precisely zero field is therefore not the finite pure-phase magnetic susceptibility used in the ratio.

For the isotropic zero-field classical Heisenberg model of a ferromagnet, the analytic vector Landau expansion uses the Cartesian components of the vector order parameter near $\mathbf{m}=0$. rotational symmetry makes its scalar terms functions of $\mathbf{m}\cdot \mathbf{m}$. Global spin inversion symmetry also requires $f(\mathbf{m})=f(-\mathbf{m})$. Through fourth order this givesA positive quartic coefficient stabilizes the unconstrained truncated free energy. The disordered phase has $b>0$ and its minimum at $m=0$; the aligned phase has $b<0$ and a minimum at nonzero $m$. At the continuous transition $b=0$. The attractive alignment interaction supplies an energetic tendency toward negative quadratic curvature; orientational entropy favors the disordered phase at high temperature, so the sign of $b$ is a free-energy question rather than a consequence of the pair energy alone.

The analyticity assumption matters: the variable printed as $m$ is a nonnegative magnitude, so reversing $\mathbf{m}$ does not send that scalar to $-m$. Terms proportional to $|\mathbf{m}|$ or $|\mathbf{m}{|}^3$ would be rotationally invariant but nonanalytic at the vector origin. Standard analytic Landau theory excludes them. Without that assumption, the stated pair symmetry alone would not exclude odd powers of the magnitude.^[ref]

Differentiation gives $f^{\prime }(m)=2m(b+2d{m}^2)$. For $b>0$ the origin is the global minimum; for $b<0$, $m=\sqrt{-b/(2d)}$ has $f^{\prime \prime }=-4b>0$ and $f=-b^2/(4d)<f(0)$. Therefore, close enough to the transition for the quartic expansion to apply,The order parameter approaches zero continuously with mean-field critical exponent $\beta =1/2$: the transition is continuous. In the aligned phase a particular direction is chosen. The vector-spin rotation symmetry is spontaneously broken from $SO(3)$ to $SO(2)$, or from $O(3)$ to $O(2)$ including reflections. The manifold of vector order parameter directions is a sphere, so this is not just a two-state sign choice.^[ref]

The plot uses illustrative $d=1/2$ and temperature measured in units of $T_c$. A unit-spin average must satisfy $m\le 1$; if the quartic model is used literally far below the transition, its constrained minimum is $\min \{1,\sqrt{(T_c-T)/(2d{T}_c)}\}$. In practice the near-transition expansion is not a controlled prediction of deep low-temperature saturation. Also, a finite zero-field statistical ensemble that averages all directions has $⟨\mathbf{m}⟩=0$ even below $T_c$; the curve describes the magnitude in a selected ordered state, obtained by appropriate boundary conditions or the thermodynamic limit before removing a small magnetic field.

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.

Set $q=e^{-\beta \hslash \omega }$. The oscillator populations are $(1-q)q^n$, and differentiating the geometric series gives $\bar{n}=q/(1-q)$. The canonical oscillator coordinate variance follows from $2\bar{n}+1=(1+q)/(1-q)=\coth (\beta \hslash \omega /2)$, andIt approaches $\hslash /(2m\omega )$ at zero temperature and $k_{B}T/(m{\omega }^2)$ at high temperature, agreeing respectively with ground-state fluctuations and classical equipartition.

For the canonical quantum Gibbs state,If $H|n⟩=E_n|n⟩$, its spectral form is $\rho =Z^{-1}{\sum }_n{e}^{-\beta E_n}|n⟩⟨n|$, with a sum over any degeneracy and an integral where a continuous spectrum needs one. It is Hermitian, positive and has unit trace, and it commutes with $H$. The formula presumes a well-defined finite partition function. At zero temperature a nondegenerate stable system approaches its ground-state projector; a degenerate ground space requires specifying the ensemble or preparation.

Specify the magnetic coupling before fixing its sign: let the magnetic moment be $\mathbf{\mu }=\gamma \mathbf{S}$, where $\gamma$ is a signed gyromagnetic ratio. The Zeeman effect gives $H_Z=-\gamma \mathbf{B}\cdot \mathbf{S}$. For spin one-half, $\mathbf{S}=\hslash \mathbf{\sigma }/2$, so with $B=|\mathbf{B}|$ and $u=\beta \gamma \hslash B/2$,This is the quantum Gibbs state. Its eigenvalues are $(1\pm \tanh u)/2$, so its purity is $(1+{\tanh }^{2}u)/2$. For nonzero level splitting it approaches a pure state as $|u|\to \infty$: positive inverse temperature selects the lower energy, while negative inverse temperature selects the upper one. It is a maximally mixed state at $u=0$, including infinite temperature or zero magnetic splitting. At zero field the Gibbs state remains $I/2$; lowering temperature alone does not select a unique pure state from that degeneracy.

Negative temperature is meaningful for an internally equilibrated spin sector with a bounded energy spectrum and population inversion. It must be sufficiently isolated from unbounded motional degrees of freedom: the oscillator's unbounded spectrum, for instance, makes its partition sum diverge at negative inverse temperature. The two-level spin partition sum itself stays finite for either sign.

For $\mathbf{B}=B{\mathbf{e}}_z$, the quantum Gibbs state is diagonal in the $S_z$ basis. Its two probabilities are proportional to $e^{\pm \beta \gamma \hslash B/2}$. ThusAt small inverse temperature, $⟨S_z⟩\simeq \beta \gamma {\hslash }^{2}B/4$. As a function of $u=\beta \gamma \hslash B/2$, the normalized spin component is the smooth odd curve $\tanh u$: it crosses zero at $u=0$ and saturates at $\pm 1$. As a function of $\beta$ its direction is reversed if $\gamma B<0$. For $\gamma B>0$, positive temperature favors $S_z=+\hslash /2$, while negative temperature gives the inverted population.

The fixed-particle-number sector is spanned by $|1⟩=c_1^{†}|0⟩$ and $|2⟩=c_2^{†}|0⟩$. In this basis,The relevant quantum Gibbs state is normalized within this two-dimensional sector, not over all particle numbers. With $\beta =(k_{B}T{)}^{-1}$ and $x=\beta t$, the exponential of a Pauli vector givesThereforeIts eigenstates are $|\pm ⟩=(|1⟩\pm |2⟩)/\sqrt{2}$, with energies $\mp t$ and probabilities $p_{\pm }=e^{\pm x}/(2\cosh x)$. These are nonnegative and sum to one, so the matrix is Hermitian and has trace one. Its purity of a density operator isFor $t\ne 0$, the zero-temperature limit is the pure ground-state projector,with purity one. As $T\to \infty$, $x\to 0$ and $\rho \to I/2$, with purity $1/2$. At every strictly positive finite temperature and finite $t$, both eigenvalues are nonzero, so the state is mixed. If $t=0$, the two levels are degenerate and the Gibbs state is $I/2$ at all temperatures, including its zero-temperature limit.

Using dimensionless Von Neumann entropy, $S=-\mathrm{Tr}(\rho \ln \rho )$, givesFor nonzero $t$, $S\to 0$ as $T\to 0$, while $S\to \ln 2$ as $T\to \infty$. More precisely, the high-temperature expansion is $S=\ln 2-x^2/2+O(x^4)$, and the low-temperature form is $S\sim (1+2|x|)e^{-2|x|}$. The thermodynamic entropy includes an overall factor $k_B$.

Take the two magnetic moments along $z$ to be $\pm \mu$, so $H=-\mu B_z{\sigma }_z$ and the two energies are $\mp \mu B_z$. The quantum Gibbs state at inverse temperature $\beta =1/(k_{B}T)$ isIts Bloch vector is $(0,0,\tanh \beta \mu B_z)$. At $B_z=0$, or in the high-temperature limit, it approaches $I/2$; at $T\to 0^{+}$ with nonzero $\mu B_z$, it approaches the lower-energy pure state. Here $\mu$ is the magnitude or signed coefficient multiplying ${\sigma }_z$ in the magnetic Hamiltonian, rather than a moment with an additional unstated factor of two.

For each pair of adjacent spins, define the transfer matrix for a classical spin chain by $T_{s{s}^{\prime }}=e^{-\beta E_0(s,s^{\prime })}$. Multiplying its entries gives the Boltzmann factor of a whole configuration. The periodic boundary condition identifies the last index with the first, soThe assumed interchange symmetry of $E_0$ makes $T$ a real symmetric matrix, hence it has an orthonormal eigenbasis. At finite positive temperature, every matrix entry is strictly positive, so the Perron–Frobenius theorem supplies a unique positive largest eigenvalue ${\lambda }_{+}$, larger in modulus than every other one. Therefore $L^{-1}\log Z_L\to \log {\lambda }_{+}$.

Local observables become matrix insertions. For example, with $S_{s{s}^{\prime }}=s{\delta }_{s{s}^{\prime }}$,Thus partition functions and correlation functions reduce to a fixed-dimensional eigenvalue problem, rather than a sum over exponentially many microstates.

Let $S=\mathrm{diag}(1,-1)$ be the base-spin insertion matrix. In the eigenbasis from solution, its diagonal matrix elements vanish: $⟨+|S|+⟩=⟨-|S|-⟩=0$. Hencenot only at large $L$ but for every finite ring at zero field.

This vanishing is not merely an averaging artifact hiding a finite-temperature ordered phase. The two-point correlation function is $⟨{\sigma }_1{\sigma }_{1+r}⟩=({\lambda }_{-}^r{\lambda }_{+}^{L-r}+{\lambda }_{+}^r{\lambda }_{-}^{L-r})/Z_L$. In the thermodynamic limit it becomes $({\lambda }_{-}/{\lambda }_{+}{)}^r$, with finite correlation length ${\xi }^{-1}=\log ({\lambda }_{+}/{\lambda }_{-})$ for every $T>0$ and $J>0$. Thermal domain walls destroy long-range alignment in this one-dimensional short-range system. At zero temperature, the two fully aligned ground states are degenerate; that limiting case must be distinguished from any fixed positive temperature.

Insert $D=\mathrm{diag}(1,-1)$ at the first undecorated site. In the eigenbasis of the transfer matrix, $D$ exchanges $v_{+}$ and $v_{-}$, and both diagonal matrix elements vanish. Explicitly,Thus its large-$L$ limit is zero at every fixed $T>0$. The result is exact even before taking that limit. Equivalently, simultaneous reversal of all spins pairs configurations of equal energy and opposite ${\sigma }_1$. Positivity of the transfer matrix at fixed positive temperature prevents a degenerate dominant eigenvalue and spontaneous magnetization in this one-dimensional zero-field chain. Selecting one aligned state at strictly zero temperature is a different order of limits.

Insert the decorating-spin observable into the one bond that contains it. The decorated-spin occupation insertion uses the matrix:The partition function still contains $L$ decorating-spin sums; one of them has been replaced by this weighted sum. HenceTaking the dominant eigenvector givesThis is the probability that the decorating site is in either nonzero state. It tends to $2/3$ at high temperature, when its three states are equally likely, and to one at low temperature, when the nonzero aligned states gain bond energy. Unlike the magnetization, this even observable is not forced to vanish by spin reversal.

For low temperature, factor out the largest Boltzmann factor from the dominant eigenvalue:Thus the Helmholtz free energy per site has the expansionDifferentiate at fixed $J$, remembering that $K=J/(k_{B}T)$. The entropy per site isThe nonzero residual entropy comes from an exponentially large ground-state degeneracy, generated by geometrical frustration. A triangle with antiferromagnetic energy $J[st+(s+t)\tau ]$ has minimum energy $-J$: if $s=t$, the apex must have the opposite sign; if $s=-t$, either apex sign is permitted. The two all-equal configurations instead have energy $3J$. Thus exactly one of a minimum-energy triangle's bonds is frustrated. The bridge can simultaneously minimize its energy at $-J$ by requiring $t=-s^{\prime }$, where $s^{\prime }$ is the next cell's first base spin.

For fixed $s$ and $s^{\prime }$, the number of allowed apex choices is therefore two when $s=s^{\prime }$ and one when $s=-s^{\prime }$. The zero-temperature counting transfer matrix on successive first base spins isThe eigenvalues of $G$ are $3$ and $1$. Dividing $k_B\ln g_M$ by the $3M$ sites and taking the thermodynamic limit gives $k_B\ln 3/3$, independently confirming the derivative calculation. This does not assert a unique ordered ground state: the extensive choice of frustrated bonds is precisely what leaves entropy at zero temperature. At high temperature, all Ising spins become independent and $s\to k_B\ln 2$ per site.

In the $s=+1,-1$ basis the transfer matrix and its eigenvalues areFor finite positive temperature, $a,b>0$ and ${\lambda }_{+}>|{\lambda }_{-}|$. Consequently $Z={\lambda }_{+}^L+{\lambda }_{-}^L$ and $L^{-1}\log Z\to \log {\lambda }_{+}$ in the thermodynamic limit. There are $2L$ spins, so the Helmholtz free energy per spin isThe factor $1/2$ is essential: one transfer matrix step counts two spins. For completeness, if the undefined $r$ is retained, set $a_r=A_{-1,-1}^{(r)}$, $d_r=A_{r+1,r+1}^{(r)}$ and $b_r=A_{-1,r+1}^{(r)}$. The same thermodynamic limit gives $f_r=-(2\beta {)}^{-1}\log [(a_r+d_r+\sqrt{(a_r-d_r{)}^2+4b_r^2})/2]$ whenever these are two distinct spin values.