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

Take a uniform nearest-neighbour harmonic chain with displacement $u_j$, conjugate momentum $p_j$, mass $m_0$, spring stiffness $\kappa$, lattice spacing $a$ and periodic length $\ell =Na$. These specify a standard model for the requested notes, not additional numerical exam data. Its Lagrangian and Hamiltonian areVariation gives $m_0{\ddot{u}}_j=\kappa (u_{j+1}-2u_j+u_{j-1})$. For a mode $e^{i(kja-\omega t)}$, periodicity requires $k=2\pi n/\ell$ andThe lattice makes $k$ periodic modulo $2\pi /a$, so choose one Brillouin zone. At long wavelength ${\omega }_k\simeq v|k|$, where $v=a\sqrt{\kappa /m_0}$ is the sound speed. The continuum limit of a harmonic chain expands $u_{j+1}-u_j\simeq a{\partial }_{x}u$ and replaces ${\sum }_j$ by $a^{-1}\int dx$. Defining mass density $\rho =m_0/a$ and one-dimensional elastic coefficient $K=\kappa a$ givesThe momentum density is $\Pi =\rho {\partial }_{t}u$, and the Hamiltonian density is ${\Pi }^2/(2\rho )+K({\partial }_{x}u{)}^2/2$. Translational invariance of $u$ makes total momentum conserved. The continuum approximation describes acoustic wavelengths large compared with $a$; it is not a justification for extending the lattice dispersion and zero-point sums to arbitrarily large wavevector.^[ref]

canonical quantization promotes the coordinates and momenta to operators with $[u_j,p_l]=i\hslash {\delta }_{jl}$ and the other equal-time commutators zero. In the continuum, $[u(x),\Pi (x^{\prime })]=i\hslash {\delta }_{\mathrm{per}}(x-x^{\prime })$. The nonzero normal modes can be represented by independent phonon operators obeying $[b_k,b_q^{†}]={\delta }_{kq}$. A useful lattice expansion iswhere $X=N^{-1}{\sum }_j{u}_j$ is the centre-of-mass displacement. The corresponding momentum isIncluding the centre-of-mass term restores the full site commutator: the nonzero-mode sum alone gives $i\hslash ({\delta }_{jl}-1/N)$. Substitution diagonalizes the Hamiltonian,A phonon is one quantum of a collective normal mode, carrying energy $\hslash {\omega }_k$ and lattice momentum $\hslash k$ modulo a reciprocal-lattice vector. Reality of $u_j$ relates its Fourier coordinates by $u_k^{†}=u_{-k}$; it does not remove the independent oppositely travelling phonon modes. Using one complete Brillouin zone avoids double-counting real standing-wave coordinates.^[ref]

The free centre of mass has no restoring force and must not be treated as an oscillator of zero frequency. Fixing it, weakly pinning it, or specifying a finite centre-of-mass volume defines the corresponding zero-mode factor. With a fixed centre of mass, the internal thermal partition function and occupations areThe internal energy is ${\sum }_{k\ne 0}\hslash {\omega }_k(⟨n_k⟩+1/2)$; its high-temperature heat capacity approaches $(N-1)k_B$. For a large one-dimensional chain at temperatures below the lattice cutoff but above its finite-size mode spacing,This one-dimensional acoustic phonon heat capacity shows that quantization turns classical normal waves into bosonic excitations with a linear low-temperature heat capacity in one dimension. The Euclidean field action $S_E=\frac{1}{2}{\int }_0^{\beta \hslash }d\tau dx[\rho ({\partial }_{\tau }u{)}^2+K({\partial }_{x}u{)}^2]$ is Gaussian, with periodic imaginary-time boundary conditions; its Fourier modes reproduce these oscillators. Absolute displacement correlations require the zero-mode convention, and the lattice regulates short-distance fluctuations.

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.