Past exam of the quantum Condensed Matter Field Theory course of the University of Cambridge / 2015 / Paper 1 / 2

The Holstein–Primakoff transformation converts a spin algebra into bosonic operators while retaining the finite physical spin Hilbert space. With $[a_n,a_m^{†}]={\delta }_{nm}$, writeThe Holstein–Primakoff occupation constraint is $0\le n_n\le 2S$. Acting on $|n⟩$, the raising matrix element is $\sqrt{n(2S-n+1)}$ and the lowering element is $\sqrt{(n+1)(2S-n)}$. Their squared difference is $2(S-n)$, reproducing $[S^{+},S^{-}]=2S^z$. Operator ordering inside the square roots is essential. For low boson occupation relative to $S$, expand them to obtain $S^{+}\simeq \sqrt{2S}a$ and $S^{-}\simeq \sqrt{2S}{a}^{†}$. Keeping quadratic terms gives the linear spin-wave approximation; higher terms describe magnon interactions. Large $S$ helps control the expansion, but a low-dimensional infrared divergence can invalidate it even at large finite $S$.

For the ferromagnet choose $H_F=-J{\sum }_n{\mathbf{S}}_n\cdot {\mathbf{S}}_{n+1}$, $J>0$, with unit lattice spacing and $\hslash =1$. Expand around the all-up state. A longitudinal bond contributes $-J{S}^2+JS(n_n+n_{n+1})$, and the transverse exchange contributes $-JS(a_n^{†}{a}_{n+1}+a_{n+1}^{†}{a}_n)$. After a Fourier transform,These are ferromagnetic magnons. The polarized vacuum is an exact ground state of the full ferromagnetic Hamiltonian: every bond is in its maximum total-spin sector, minimizing that bond energy. It therefore has no quantum zero-point reduction of its magnetization. The one-magnon dispersion is also exact, although interacting many-magnon states need not be free bosons. Rotating the polarization generates the maximum-total-spin ground-state multiplet. Selecting a definite direction in a thermodynamic symmetry-breaking description breaks $SU(2)$ to $U(1)$; the two transverse spin components form a canonical pair, producing a single quadratic type-B Goldstone boson. A small aligning field gaps it by the Zeeman energy.

At nonzero temperature the one-dimensional ferromagnet is different: its thermal magnon number contains $\int dk[e^{\beta JS{k}^2}-1{]}^{-1}$, with infrared behaviour $\int dkT/k^2$. The divergent magnetization reduction signals the absence of finite-temperature long-range order in an infinite short-range isotropic chain, consistent with the Mermin-Wagner theorem. It does not disallow the polarized zero-temperature ground states just constructed.^[ref]

For the antiferromagnet use $H_{AF}=+J{\sum }_n{\mathbf{S}}_n\cdot {\mathbf{S}}_{n+1}$, still $J>0$, on an even bipartite chain. The classical reference state alternates up and down. Rotate one sublattice through $\pi$ so that both local quantization axes point along their reference spins. In unrotated physical components, to leading order $S_A^z=S-a^{†}a$, $S_B^z=-S+b^{†}b$, $S_A^{+}\simeq \sqrt{2S}a$ and $S_B^{+}\simeq \sqrt{2S}{b}^{†}$. Choosing Fourier phases using the actual site positions givesThe reduced Brillouin zone has $M=N/2$ momenta and length $\pi$. The pair creation and annihilation terms express the fact that an antiparallel reference product state is not an eigenstate. A bosonic Bogoliubov transformation with $u_k=\cosh {\theta }_k$, $v_k=\sinh {\theta }_k$, $\tanh 2{\theta }_k=\cos k$ diagonalizes each two-mode block. With the sign convention $a_k=u_k{\alpha }_k-v_k{\beta }_{-k}^{†}$ and $b_{-k}=u_k{\beta }_{-k}-v_k{\alpha }_k^{†}$, the result isThus the harmonic antiferromagnetic modes have linear, rather than quadratic, long-wave dispersion. In an unfolded Brillouin-zone description the same result is represented by zeros at both zero and the antiferromagnetic wavevector; these are not extra independent degrees of freedom beyond the two-sublattice accounting.

The Bogoliubov vacuum has no $\alpha$ or $\beta$ excitations but contains original spin bosons, with $v_k^2=\frac{1}{2}(|\sin k{|}^{-1}-1)$. Its staggered magnetization is $m_s=S-M^{-1}{\sum }_{k\in \mathrm{RBZ}}{v}_k^2$. This quantum depletion of Néel order diverges logarithmically at small $k$ in one dimension. Linear spin-wave theory therefore does not establish a stable ordered antiferromagnetic ground state for this chain. Its nominal gapless dispersion is a result of expanding around an invalid ordered reference, not a general proof of the exact excitation spectrum. For the nearest-neighbour isotropic chain, integer-spin systems have a Haldane gap, whereas half-integer-spin systems are gapless; the spin-one-half chain in particular has fractional spin excitations rather than a simple dilute gas of ordered-state magnons. This distinction requires physics beyond the harmonic expansion, including the topological term in the continuum theory.^[ref]

For one coordinate with Hamiltonian $H={\hat{p}}^2/(2m)+V(\hat{x})$, the propagator is $K(x_f,t_f;x_i,t_i)=⟨x_f|e^{-iHT/\hslash }|x_i⟩$, where $T=t_f-t_i>0$. Divide $T$ into $M$ intervals $\delta t$, use the Trotter product formula to separate kinetic and potential evolution, and insert complete position and momentum bases at each slice. The momentum integral is Gaussian:Multiplying the short-time kernels gives the time-sliced configuration-space path integral,with $x_0=x_i$, $x_M=x_f$. Continuum notation abbreviates this limit as ${\int }_{x_i}^{x_f}\mathcal{D}x{e}^{iS[x]/\hslash }$, where $S=\int dt[m{\dot{x}}^2/2-V(x)]$. The measure includes the displayed normalization; it is not an unweighted count of paths. A small damping or analytic-continuation prescription defines the oscillatory integrals. The paths contribute complex amplitudes, so alternative histories interfere before probabilities are calculated. This is the central difference from a classical statistical mixture of trajectories.^[ref]

When actions vary rapidly relative to $\hslash$, stationary phase selects classical paths. Put $x=x_{\mathrm{cl}}+\eta$, with $\eta$ zero at the endpoints. Vanishing first variation gives $m{\ddot{x}}_{\mathrm{cl}}+V^{\prime }(x_{\mathrm{cl}})=0$. The quadratic term isThe resulting Gaussian fluctuation determinant is the semiclassical propagator prefactor. In one dimension, away from caustics,The branch convention makes $\nu =0$ for a direct free-particle path before any caustic. The Maslov index tracks the phase changes as fluctuations pass through a focal singularity. In several dimensions replace the endpoint derivative by the Van Vleck determinant. The approximation sums all relevant classical paths; real-time action need only be stationary, not a minimum. It fails in this simple form at a caustic or coalescing saddles, where a uniform approximation is needed. Quadratic actions have no higher fluctuation terms, making their Gaussian evaluation exact. For example $S_{\mathrm{cl}}=m(x_f-x_i{)}^2/(2T)$ gives the exact free-particle kernel $\sqrt{m/(2\pi i\hslash T)}{e}^{i{S}_{\mathrm{cl}}/\hslash }$. Tunnelling involves classically forbidden regions and can be described by imaginary-time saddles rather than real classical trajectories.

Quantum statistical mechanics follows by Wick rotation $t=-i\tau$. The imaginary-time evolution operator is $e^{-\tau H/\hslash }$, with Euclidean actionTaking a thermal trace identifies and integrates the endpoint positions:This thermal trace as a closed Euclidean path integral explains the periodic imaginary-time direction. A density-matrix element instead retains open endpoints. For many identical particles the trace also sums over endpoint permutations, with a plus sign for bosons and permutation parity for fermions. In coherent-state language this becomes periodic bosonic fields and antiperiodic Grassmann fields, with corresponding Matsubara frequencies.

The Euclidean weight resembles a classical Boltzmann weight. A time-sliced local quantum field in $d$ spatial dimensions becomes a classical-looking field model in $d+1$ dimensions, with anisotropic couplings and extra extent $\beta \hslash$. Thus zero-temperature quantum criticality can be treated using an infinitely extended imaginary-time direction, while nonzero temperature imposes a finite thickness. The dynamical scaling of time need not equal the spatial scaling: the dynamical critical exponent specifies their relationship. This quantum-to-classical imaginary-time mapping requires qualification when weights have signs or imaginary terms; such a functional integral is not automatically an ordinary positive classical probability distribution. Even with a positive weight, the mapping concerns equilibrium statistics, not literal dissipative real-time trajectories. For a free particle the Euclidean kernel is a diffusion kernel with diffusion coefficient $\hslash /(2m)$, illustrating the mathematical correspondence.

The classical high-temperature limit is a further approximation, rather than simply another name for imaginary time. If quantum level spacings and exchange effects are negligible, a one-particle partition function in $d$ dimensions becomesFor indistinguishable particles there is a factor $1/N!$ in the dilute classical gas, and exchange is negligible when the number density times the thermal de Broglie wavelength to power $d$ is small. For a harmonic oscillator, the exact quantum result $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}$ tends to $1/(\beta \hslash \omega )$ as $\beta \hslash \omega \to 0$, recovering the classical phase-space integral. stationary phase explains classical mechanics; the short imaginary-time, weak-exchange limit explains classical statistical mechanics. These are related but distinct limits of the same quantum formulation.

Consider a conventional spin-singlet, weak-coupling, three-dimensional superconductor. Set $\hslash =k_B=1$ in the microscopic action and restore electromagnetic units below. An attraction $\lambda >0$ in an energy shell about the Fermi surface gives a Grassmann actionThe attractive sign permits pairing even when the single-particle Fermi gas is otherwise stable. A complex Hubbard–Stratonovich transformation introduces a pair field $\Delta$ proportional at the saddle to $\lambda ⟨{\psi }_{\downarrow }{\psi }_{\uparrow }⟩$. The action is then quadratic in fermions. Integrating them out gives the BCS pairing-field effective actionHere the trace uses one independent Nambu spinor pairing block; a fully doubled particle-hole representation instead carries its corresponding one-half factor. At a uniform saddle the Bogoliubov--de Gennes Hamiltonian can be written asThe negative branch describes the paired vacuum; excitations are Bogoliubov quasiparticles mixing particles and holes. Extremizing the effective action givesThe physical cutoff or a renormalized interaction is part of the microscopic model; an unregulated contact integral must not be used indiscriminately. With constant density of states $N(0)$ per spin and shell half-width ${\Omega }_c$, the zero-temperature weak-coupling solution is ${\Delta }_0\simeq 2{\Omega }_c{e}^{-1/[N(0)\lambda ]}$. The exponentially small scale follows from the logarithmic pairing susceptibility of states near the Fermi surface. For this ideal weak-coupling model, $2{\Delta }_0/(k_B{T}_c)\simeq 3.53$. The microscopic theory relates the gap, pairing and thermodynamic transition.^[ref]

Close to a continuous transition, for fields varying slowly compared with microscopic scales, expand the effective action in powers and gradients of $\Delta$. Its static part is the Ginzburg-Landau theory of superconductivity:The normalization of $\Delta$ fixes the coefficients; here $a(T)$ changes sign at $T_c$. The uniform minimum has $|{\Delta }_0{|}^2=-a/b$ for $a<0$, condensation free-energy density $-a^2/(2b)$, and characteristic coherence-length scale $\xi =\sqrt{K/|a|}$. Above the transition the minimum is at zero. The quartic stabilization and gradient cost are precisely the ingredients of Landau theory for a complex order parameter. This controlled near-transition expansion supplies a bridge between microscopic BCS theory and macroscopic magnetic response; its coefficients and validity depend on the material and fluctuation regime.^[ref]

Write $\Delta =(|{\Delta }_0|+\eta )e^{i\varphi }$. The radial fluctuation has positive local restoring curvature, whereas the phase has no local potential in a neutral theory with global $U(1)$ symmetry. Freezing the amplitude givesThe neutral case supports a gapless collective phase mode. In the charged case the electromagnetic field couples to this phase rigidity. Under $\mathbf{A}\mapsto \mathbf{A}+\nabla \chi$ and $\varphi \mapsto \varphi +q\chi /\hslash$, the displayed combination is invariant. Locally choosing a gauge with constant phase leaves a quadratic term in $\mathbf{A}$. Varying the free energy gives the transverse London current and magnetic screening scale,Combining this relation with static Maxwell equations gives exponential field penetration, the Meissner effect. Single-valuedness of the pair field around a loop with negligible current gives flux in units ${\Phi }_0=2\pi \hslash /|q|=\pi \hslash /e$. These are gauge-invariant physical consequences of the condensate stiffness.

For the dynamics, the neutral phase and density mode is reorganized by the Coulomb interaction and the dynamical electromagnetic field. In a bulk three-dimensional system the long-wave longitudinal mode has a plasma gap, rather than remaining the neutral sound mode. This is the Anderson-Higgs mechanism in a superconductor: a phase coordinate becomes part of the electromagnetic collective response instead of surviving as an independent massless physical particle. Local gauge transformations are redundancies of the description, so a literal breaking of a local redundancy is not an observable. Phase coherence, current response and magnetic screening provide the physical statements.^[ref]

A static Landau functional alone cannot determine that full dynamical spectrum. A two-dimensional film, for example, need not have the same long-wave plasma dispersion as a bulk charged system. Nor does a positive radial curvature automatically guarantee a sharp Higgs amplitude resonance: coupling to quasiparticles and other modes can damp it. The chosen mean-field phase is also a thermodynamic symmetry-breaking description; an isolated finite system can be formulated with definite particle number. The microscopic pairing saddle, the near-transition Landau expansion and the electromagnetic Anderson-Higgs response describe successive scales of the same superconducting state, with distinct assumptions at each step.