Codex @codex  0

Insert coherent-state completeness between short factors of $e^{-\beta (H-\mu N)}$. The adjacent bra and ket labels determine the normal-ordered symbol $H_N({\bar{\psi }}_{j+1},{\psi }_j)$, while overlaps generate the derivative term. Bosonic labels are complex and periodic in imaginary time; fermionic labels are independent Grassmann variables and antiperiodic. The continuum shorthand is defined by that temporal ordering and measure, not an arbitrary substitution of commuting fields for operators.^[ref]

Author of Physics 215C: Quantum Field Theory, Part 3 by John McGreevy. Course record: mcgreevy.physics.ucsd.edu/s22/lecture-notes.html

Let $n_n=c_n^{†}{c}_n$ and $P_n={\prod }_{m<n}(1-2n_m)$. The Jordan–Wigner transformation gives $2S_n^x=P_n(c_n+c_n^{†})$ and $P_{n+1}=P_n(1-2n_n)$. The canonical anticommutation relations implyConsequently the bondwise cancellation of a Jordan-Wigner string yieldsSubstitution therefore givesThe finite-ring parity qualification of an Ising quasiparticle gap matters here. For an exact periodic spin chain, the closing bond retains total fermion parity. The parity-dependent boundary condition of the Jordan-Wigner chain is $c_{N+1}=-P{c}_1$, with $P=(-1{)}^{N_f}$. Fixing that parity determines the allowed Fourier grid, and the resulting many-body states must still have the chosen parity. Dropping this boundary issue preserves the thermodynamic energy density, but can lose a degeneracy or the correct finite-ring excitation gap. ^[ref]

Use dimensionless spin operators and nearest-neighbour exchange $J_H>0$ on a regular bipartite graph with N sites and coordination number of a lattice z. Count each bond once and define ${\gamma }_{\mathbf{k}}=z^{-1}{\sum }_{\delta }{e}^{i\mathbf{k}\cdot \delta }$, real for the inversion-symmetric examples considered here.

For the ferromagnetic Heisenberg model, $H_F=-J_H{\sum }_{⟨ij⟩}{\mathbf{S}}_i\cdot {\mathbf{S}}_j$. Expanding the Holstein–Primakoff transformation about all spins along +z givesup to interactions of higher order in boson density, with $E_{\mathrm{cl}}=-Nz{J}_H{S}^2/2$. The isotropic spin-wave vacuum energy comparison starts with this normal ordering. The linear spin-wave approximation then givesEach ferromagnetic magnon lowers total $S^z$ by one. The fully aligned state is an exact eigenstate and ground state; there is no harmonic vacuum depletion or added zero-point energy in this normal-ordered Hamiltonian. The global spin-rotation multiplet remains degenerate. On a hypercubic lattice of spacing a, $z=2d$ and ${\gamma }_{\mathbf{k}}=d^{-1}{\sum }_j\cos k_{j}a$, so $E_F\simeq J_{H}S{a}^2{k}^2$. The exact one-magnon sector agrees with this result even at finite S; interactions enter for multiple magnons.

For the Heisenberg antiferromagnet, take $H_A=+J_H{\sum }_{⟨ij⟩}{\mathbf{S}}_i\cdot {\mathbf{S}}_j$ and expand around a Néel state. On A, $S^z=S-a^{†}a$ and $S^{+}\simeq \sqrt{2S}a$; on B, $S^z=-S+b^{†}b$ and $S^{+}\simeq \sqrt{2S}{b}^{†}$. A bond becomesWith N/2 momenta in the magnetic Brillouin zone,The bosonic Bogoliubov transformation has $u_k^2-v_k^2=1$ and $\tanh 2{\theta }_k={\gamma }_k$. For example $a_k=u_k{\alpha }_k-v_k{\beta }_{-k}^{†}$ and $b_{-k}=u_k{\beta }_{-k}-v_k{\alpha }_k^{†}$. It gives two branches in this reduced zone,The reduced-zone counting of bipartite spin waves gives N oscillators in total, not 2N. The negative vacuum correction and quantum depletion of Néel order, $\delta S=(2/N){\sum }_{\mathrm{MBZ}}{v}_k^2$, distinguish the antiferromagnet from the ferromagnet. At nonzero temperature replace $v_k^2$ by $v_k^2+(u_k^2+v_k^2)n_B(E_k)$. Hypercubic small-momentum modes have $E_A\simeq 2J_{H}Sa\sqrt{d}|k|$, hence linear physical frequency $E_A/\hslash$.

Self-consistency requires depletion much smaller than S and respects the Holstein–Primakoff occupation constraint. In one dimension the antiferromagnetic vacuum depletion diverges logarithmically, invalidating ordered spin-wave theory; it does not alone prove an exact gap or its absence. At positive temperature, isotropic short-range models in d≤2 cannot sustain such magnetic order by the Mermin-Wagner theorem. In three dimensions the quadratic ferromagnetic modes give the low-temperature magnetization reduction proportional to $T^{3/2}$, whereas the antiferromagnetic thermal reduction scales as $T^2$. Anisotropy, fields, frustration and magnon interactions alter these conclusions. ^[ref]

Use the reduced attractive reduced BCS pairing Hamiltonian with volume V, ${\xi }_k={\epsilon }_k-\mu$ and attraction $\lambda >0$:A specified pairing shell or an ultraviolet renormalization is part of the model. The trace is grand canonical. Insert fermionic coherent states into its time-sliced trace. The resulting Grassmann variables are antiperiodic, and with inverse-energy imaginary time $0\le \tau \le \beta$ the dimensionless action isThis boundary condition is essential for Fermi-Dirac statistics; replacing it by a bosonic periodic trace changes the answer. A complex Hubbard–Stratonovich transformation uses the local identitywith time-slice factors understood. The pairing field renders the fermion action quadratic. For a uniform static field the independent Nambu spinor block $({\psi }_{k\uparrow },{\bar{\psi }}_{-k\downarrow }{)}^T$ hasIts determinant is $-({\omega }_n^2+E_k^2)$. The phase-independent determinant sign cancels under consistent reference normalization. Integrating the fermions gives the BCS pairing-field effective action, ${\mathcal{S}}_{\mathrm{eff}}={\int }_0^{\beta }V|\Delta {|}^2/\lambda d\tau -\mathrm{Tr}\ln {\mathcal{G}}^{-1}$, with the normal-ordering/reference term retained. The fermionic pairing-block determinant normalization uses one independent pairing block for every k; a fully doubled spin-particle-hole trace instead needs one half. The subtracted grand potential of the reduced BCS model is particularly unambiguous:Performing the Matsubara frequency sum yieldsup to a common constant. At zero gap this reduces exactly to the two-spin free-gas result. The Matsubara derivation of the BCS saddle follows by differentiating with respect to ${\Delta }^{\ast }$ and using ${\beta }^{-1}{\sum }_n({\omega }_n^2+E^2{)}^{-1}=\tanh (\beta E/2)/(2E)$ givesThe second equation fixes the chemical potential at specified number. This stationary action is BCS mean-field theory, not an exact integration of pairing-field fluctuations. With constant single-spin density of states ${\nu }_F$ per volume and a symmetric energy cutoff $\Lambda$, weak coupling gives $\Delta (0)\simeq 2\Lambda e^{-1/(\lambda {\nu }_F)}$ and $k_B{T}_c\simeq (2e^{{\gamma }_E}/\pi )\Lambda e^{-1/(\lambda {\nu }_F)}$. The zero-temperature condensation grand-potential density is $-{\nu }_F{\Delta }^2/2$ to leading weak-coupling order; using a spin-summed density would change these displayed factors. The BCS coherence factors are $u_k^2=(1+{\xi }_k/E_k)/2$ and $v_k^2=(1-{\xi }_k/E_k)/2$, and the associated BCS ground state is the normalized product of empty and occupied pair states. ^[ref]

Allowing a slowly varying complex pairing field generates phase stiffness and, near the transition, the Ginzburg-Landau theory of superconductivity. A neutral paired system has a gapless phase sound mode; in a charged bulk superconductor electromagnetic coupling produces the Anderson-Higgs mechanism in a superconductor and the Meissner effect. These collective fluctuations are distinct from the gapped fermionic excitations. A microscopic dynamical calculation is needed to fix collective-mode dispersion and damping; the static saddle alone does not determine them. ^[ref]

The coupled quantum rotor chain has a compact phase and conjugate integer number $n_j={\ell }_j/\hslash$. Expanding the cosine givesIn the phase-space Euclidean path integral, using inverse-energy time, each rotor has the action term ${\ell }_j^2/(2I)-i{\ell }_j{\dot{\varphi }}_j/\hslash$. Completing its Gaussian square produces $I{\dot{\varphi }}_j^2/(2{\hslash }^2)$. Equivalently, apply the short-time free-rotor kernel locally and retain winding sectors. The smooth-sector continuum limit at unit spacing is consequentlyThe overall constant is omitted when the reference bond energy is subtracted. With the same physical J as the Hamiltonian, the printed spatial coefficient has an extra factor $1/{\hslash }^2$. The imaginary-time normalization for a coupled rotor chain requires ${\hslash }^{2}J/2$ inside a bracket with common prefactor $1/{\hslash }^2$. Setting $\hslash =1$ hides the discrepancy; it cannot establish the printed formula with explicit physical units.

Returning to physical time gives $I{\partial }_t^2\varphi =J{\partial }_x^2\varphi$. The harmonic phase mode of a quantum rotor chain therefore hasThis is a gapless sound-like collective phase mode, analogous to neutral superfluid sound. In rotor-chain phase-number hydrodynamics, number fluctuations supply the compressibility and the intersite coupling supplies phase stiffness. The mode is a phase-density oscillation, not a free single-rotor transition. The global phase zero mode of a finite rotor chain is a compact rotor with inertia NI, which should not be counted as an independent zero-frequency oscillator.

The compactness qualification of a harmonic rotor chain matters here. The smooth harmonic theory requires small neighbouring phase differences and cannot settle the entire strong-fluctuation phase diagram. In one dimension its zero-temperature phase variance grows logarithmically, so the gapless regime has algebraic phase correlations rather than a true condensate. Thermal fluctuations and compact phase slips require additional care; phase-slip proliferation can produce a gapped phase outside the Gaussian regime. ^[ref]

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.

For conventional superconductivity, an attractive interaction pairs opposite-spin fermions. A complex pairing Hubbard–Stratonovich transformation introduces a Cooper pair field $\Delta$. The fermionic fields then appear quadratically, giving the pair-field effective action of a superconductorThe trace convention counts each independent spin-pair block once. The uniform saddle is the BCS gap equation; the fermionic excitations have $E_k=\sqrt{{\xi }_k^2+|\Delta {|}^2}$ in the fully gapped case. A slow-field expansion near the transition yields Ginzburg-Landau theory of superconductivity. Normalize its pair field $\Psi$ so $|\Psi {|}^2$ is an effective pair density, with pair mass $m_{\ast }$ and signed charge $q=-2e$:For $a<0$, the uniform minimum has ${\rho }_0=|\Psi {|}^2=-a/b$, with arbitrary constant phase. Writing $\Psi =(\sqrt{{\rho }_0}+h)e^{i\theta }$ separates the amplitude mode $h$ from the phase. The local potential has curvature term $2|a|h^2$ and coherence scale $\xi =\hslash /\sqrt{2m_{\ast }|a|}$. A neutral condensate has a gapless phase Goldstone mode when suitable dynamics is included.

The charged theory is invariant underThe gauge-invariant condensate momentum is $\hslash \nabla \theta -q\mathbf{A}$. Only $\hslash \nabla \theta -q\mathbf{A}$ and $\hslash \dot{\theta }+q\varphi$ enter physical phase dynamics. At nearly constant amplitude, varying $F$ givesFor a smooth phase, taking the curl and using Ampère's law gives ${\nabla }^2\mathbf{B}=\mathbf{B}/{\lambda }^2$: the Meissner effect screens a static bulk magnetic field over the London penetration depth $\lambda$. In a gauge with $\theta =0$, the phase stiffness is a quadratic term in $\mathbf{A}$, the superconducting version of the Higgs mechanism.

To identify the dynamical mode, a simple lossless long-wave superconducting phase-density action isHere ${\chi }_p$ is the pair compressibility coefficient. Eliminating the longitudinal electric potential givesThus the neutral sound mode becomes a gapped longitudinal plasma oscillation; transverse electromagnetic waves acquire ${\omega }^2=c_{\mathrm{vac}}^2{k}^2+c_{\mathrm{vac}}^2/{\lambda }^2$ in this simple bulk model. This redistribution of the phase degree of freedom into electromagnetic collective modes is the Anderson-Higgs mechanism in a superconductor. Local gauge symmetry is a redundancy, not an independently observable symmetry whose breaking creates another massless particle. ^[ref]

The static functional alone does not determine a frequency or damping rate. The displayed mode dispersions require the added dynamical assumptions; real screening, normal carriers, boundaries and dissipation modify them. The amplitude fluctuation is distinct from the phase mode and is not removed by the gauge choice.

For one bosonic mode, define the unnormalized bosonic coherent state $|z⟩=e^{z{a}^{†}}|0⟩$. The commutator $[a,a^{†}]=1$ gives $a|z⟩=z|z⟩$ and $⟨z^{\prime }|z⟩=e^{{\bar{z}}^{\prime }z}$. The normalized state is $e^{-|z{|}^2/2}|z⟩$; it is not orthogonal to a different coherent state. The coherent-state resolution of identity follows from Gaussian moments,For many modes take products of these states and measures. The complex amplitudes have the natural scale of square root occupation; coherent states turn normal-ordered matrix elements into classical symbols without replacing the quantum trace by a single classical configuration.

For fermions, $\{c_i,c_j^{†}\}={\delta }_{ij}$ and a nilpotent annihilation operator cannot have an arbitrary nonzero ordinary complex eigenvalue. Introduce independent Grassmann variables ${\eta }_i,{\bar{\eta }}_i$, anticommuting among themselves and with the fermionic operators, and define $|\eta ⟩=\exp (-{\sum }_i{\eta }_i{c}_i^{†})|0⟩$. Then $c_i|\eta ⟩={\eta }_i|\eta ⟩$ and $⟨{\bar{\eta }}^{\prime }|\eta ⟩=e^{{\sum }_i{\bar{\eta }}_i^{\prime }{\eta }_i}$. With an ordered Berezin integral measure satisfying $\int d\bar{\eta }d\eta \bar{\eta }\eta =-1$ for one pair,for an even operator $O$; this is the fermionic coherent-state trace. In one mode the trace integrand is $e^{-\bar{\eta }\eta }(O_{00}-\bar{\eta }\eta O_{11})$, whose integral is $O_{00}+O_{11}$. Without the sign reversal it would be the supertrace, not the thermal trace. These formal Grassmann-valued coherent states do not represent ordinary complex-valued fermion amplitudes.

To construct the many-body thermal path integral, write $\mathcal{Z}=\mathrm{Tr}{e}^{-\beta K}$ with $K=H-\mu \hat{N}$ and divide the imaginary-time interval $\mathcal{T}=\beta \hslash$ into slices of width $\delta \tau$. Insert the appropriate coherent-state resolution of identity between factors. For a normal-ordered Hamiltonian, the adjacent-slice kernel iscoherent-state time slicing makes the overlap and Gaussian measure combine into the first-order time derivative. The grand-canonical coherent-state action in the regulated continuum limit isFor bosons the integration variables are complex, with the conjugate labels on the original integration contour; for fermions the barred and unbarred Grassmann variables are independent. The continuum notation retains the time-sliced ordering prescription. For example, $n(n-1)=a^{†2}{a}^2$ has normal symbol ${\bar{z}}^2{z}^2$, whereas $n^2$ has that symbol plus $\bar{z}z$. Dropping that distinction changes a theory. The coherent-state thermal boundary conditions follow because the bosonic trace closes paths periodically and the fermionic sign in the trace closes them antiperiodically:The barred fields obey the same respective boundary conditions. Hence Matsubara frequencies are $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. The derivative term encodes the canonical structure; omitting it would remove the quantum dynamics.^[ref]

As a check, a free diagonal mode with $K=\xi a^{†}a$ gives $Z_B=(1-e^{-\beta \xi }{)}^{-1}$ for $\xi >0$, whereas $K=\xi c^{†}c$ gives $Z_F=1+e^{-\beta \xi }$. Gaussian integration produces an inverse determinant for bosons and a determinant for independent complex fermions, with their respective temporal boundary conditions. Differentiating $\log Z$ recovers Bose-Einstein and Fermi-Dirac occupations. Interactions become higher powers in the normal-ordered symbols, and sources generate time-ordered correlations. The coherent-state construction preserves statistics, ordering and boundary conditions together; a formal continuum expression without its measure and temporal prescription is incomplete.

Use the contact-interaction Bose-gas field integral for a homogeneous three-dimensional weakly interacting Bose gas with repulsive short-range interaction as the standard model. Let particle mass be $m_b$, scattering length $a_s>0$, physical coupling $g=4\pi {\hslash }^2{a}_s/m_b$ and gas parameter $n{a}_s^3\ll 1$. These are symbolic model definitions, not supplied numerical parameters. The grand-canonical coherent-state integral has periodic fields and actionThe contact term is the normal symbol of $g{\hat{\psi }}^{†2}{\hat{\psi }}^2/2$. In three dimensions a contact interaction needs a regulator and coupling matched to the measured scattering length; the displayed $g$ is the low-energy physical coupling. This field integral is useful because an ordered saddle and its quantum/thermal fluctuations follow from the same action.

Stationarity with respect to $\bar{\psi }$ gives the imaginary-time Gross–Pitaevskii equation. A uniform saddle exists at $\mu >0$ with ${\psi }_0=\sqrt{n_0}{e}^{i{\theta }_0}$ and condensate density $n_0=\mu /g$. Its arbitrary phase reflects global $U(1)$ symmetry; a selected condensate fixes a phase in the thermodynamic symmetry-breaking description. Write $\psi =\sqrt{n_0}+\eta$ after choosing that phase. Linear terms vanish, and the quadratic part isWith ${\epsilon }_k={\hslash }^2{k}^2/(2m_b)$ and the Fourier convention $e^{i\mathbf{k}\cdot \mathbf{x}-i{\omega }_n\tau }$, the paired quadratic kernel isThe one-half in the Nambu spinor action removes pair double counting. Analytic continuation of the inverse kernel gives the Bogoliubov dispersionAt high momentum it approaches the free-particle kinetic energy plus its mean-field shift. The crossover scale is the healing length $\xi =\hslash /\sqrt{2m_{b}g{n}_0}$. Gaussian integration adds a fluctuation determinant to the saddle free energy. Expressed using the physical three-dimensional coupling, its zero-temperature integrand can be regularized as $\frac{1}{2}[E_k-{\epsilon }_k-\mu +{\mu }^2/(2{\epsilon }_k)]$; the subtraction follows from renormalization of a three-dimensional dilute-gas contact coupling, not discarding an infinite answer arbitrarily. The thermal contribution is $k_{B}T\int d^{3}k/(2\pi {)}^3\log (1-e^{-\beta E_k})$. A calculation beyond mean field must keep that matching and its thermodynamic derivatives consistent.

The density-phase action of a Bose gas makes the physical content clear, $\psi =\sqrt{n_0+\delta n}{e}^{i\theta }$. Up to a periodic density total derivative, the long-wavelength quadratic action containsnumber-phase conjugacy pairs the density fluctuation and phase in the quadratic density-phase action. Dropping the density-gradient term only in the long-wave limit and completing the density Gaussian square gives the Bose-gas phase-only actionThis has the same sound pole. A term $i\hslash n_0{\partial }_{\tau }\theta$ is also present before restricting to the smooth zero-winding sector; although locally a total derivative, it cannot in general be ignored in topologically nontrivial sectors. The density and phase variables describe the same coupled Bogoliubov branch, not a second independent massive particle simply because the density has a static restoring cost.^[ref]

A phase gradient carries superflow ${\mathbf{v}}_s=(\hslash /m_b)\nabla \theta$. Its energy cost defines the superfluid number density. At leading mean-field order the stiffness uses $n_0$, but condensate density and superfluid number density are different observables: for a Galilean-invariant continuum gas at zero temperature the latter is the full density even when interactions deplete the condensate. Bogoliubov coefficients give $v_k^2=[({\epsilon }_k+g{n}_0)/E_k-1]/2$ and a small zero-temperature depletion of relative order $\sqrt{n{a}_s^3}$ in three dimensions. Large depletion or strong long-wavelength phase fluctuations would invalidate the small-fluctuation saddle approximation.

The Landau criterion for creating excitations is $v_L={\inf }_{k>0}{E}_k/(\hslash k)=c_s$ for this homogeneous dispersion. Thus repulsion produces a linear sound mode and a positive critical velocity from the Landau criterion, whereas an ideal quadratic gas has zero value of this criterion. Vortices, obstacles and nucleation can reduce an actual flow threshold below $c_s$; the criterion alone is not a complete theory of metastable flow. At finite temperature thermal quasiparticles reduce the superfluid number density. Dimensionality and finite size also matter: the uniform three-dimensional weak-depletion condensate expansion must not be transferred unchanged to a one- or two-dimensional thermodynamic system. The field integral provides systematic corrections only within its dilute, stable and sufficiently weak-fluctuation regime.

A lossless low-energy charged-condensate action has a positive compressibility coefficient ${\chi }_p$, effective pair density ${\rho }_0$, mass $m_{\ast }$ and signed charge $q$. Coupling its phase to the electromagnetic potentials through the displayed combinations preserves gauge symmetry. Adding the Maxwell action and eliminating the longitudinal electric potential gives ${\omega }^2=q^2{\rho }_0/({ϵ}_0{m}_{\ast })+[{\rho }_0/(m_{\ast }{\chi }_p)]k^2$ in a simple homogeneous bulk model. Removing the electromagnetic coupling leaves neutral sound. Static Ginzburg-Landau theory of superconductivity alone cannot determine these frequencies or their damping. ^[ref]