Past exam of the mathematics course of the University of Cambridge / 2016 / iii / Paper 301 (Quantum Field Theory)

Use units $\hslash =c=1$, the Minkowski metric ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$, and the Fourier transform convention $f(x)=\int d^{4}p{e}^{-ip\cdot x}\tilde{f}(p)/(2\pi {)}^4$. Here gamma matrices satisfy $\{{\gamma }^a,{\gamma }^b\}=2{\eta }^{ab}{I}_4$, and $\bar{\psi }={\psi }^{†}{\gamma }^0$ is the Dirac adjoint. A single Dirac field describes both electrons and Positrons; these are the particle and antiparticle sectors of that field.

The free Maxwell Lagrangian and Dirac action, with a linear covariant gauge condition for the photon, giveBefore gauge fixing, the Maxwell Lagrangian has zero modes in field theory along $A_a\mapsto A_a+{\partial }_a\lambda$, so its quadratic operator cannot be inverted on all potentials. In the linear Lorenz gauge, the Faddeev-Popov determinant is $\det \square$. It is independent of $A$ and may be absorbed into the normalization; the corresponding Faddeev-Popov ghost fields have no interacting vertices in this Abelian linear gauge. A residual gauge symmetry is removed by the specified boundary conditions.

The free generating functional isThe Dirac field variables and their sources are independent Grassmann fields in this path integral. They anticommute; replacing them by ordinary commuting fields would give the wrong statistics and the wrong functional determinant. The bosonic Gaussian functional integral contributes an inverse square root of a determinant, and the Grassmann Gaussian integral contributes a determinant. If $K_A$ is the photon quadratic operator and $K_D=i{\gamma }^a{\partial }_a-m$, completing the square givesIn the last expression $D_F=i{K}_A^{-1}$ and $S_F=i{K}_D^{-1}$ are the photon propagator and Dirac propagator, with Feynman i-epsilon prescription. The products include the appropriate spacetime integrals and index contractions. Functional derivatives with respect to $J$, and consistently ordered left or right Grassmann derivatives with respect to the fermionic sources, generate the time ordering of the corresponding fields. In Feynman gauge, $\alpha =1$, the momentum-space two-point functions areThe Feynman i-epsilon prescription specifies vacuum boundary conditions rather than an arbitrary inverse of the differential operator.

The covariant photon propagator uses four potential components. Its operator formalism counterpart is Gupta-Bleuler quantization: impose $({\partial }_a{A}^a{)}^{(+)}|\mathrm{phys}⟩=0$ and take the Gupta-Bleuler null-state quotient. This leaves a positive physical state space with two transverse photon polarization vectors. The temporal and longitudinal oscillator components occur in intermediate covariant expressions; the Ward identity removes their dependence from physical amplitudes. The free Dirac field uses the canonical anticommutation relations, producing the same fermionic signs as its Grassmann fields in the path integral.

The operator-path-integral equivalence can be seen directly with a regulator. Divide time into small intervals and insert complete sets of field-coordinate states for bosons, and resolutions in fermionic coherent states for fermions. The bosonic matrix elements produce the phase-space factor $\exp i\sum [\pi \Delta \varphi -H\Delta t]$; integrating out the quadratic canonical momentum produces the bosonic action. The fermionic coherent state overlaps produce the first-order term $i{\psi }^{†}\dot{\psi }$ and the Berezin integral measure. Multiplying the short-time kernels recovers the path integral. Projecting the remote endpoints onto the Fock vacuum with an infinitesimal damping selects the same Feynman propagators as the operator formalism. Field insertions become time-ordered products under this construction. Conversely, their quadratic generating functional obeys the canonical free-field equations and has precisely the oscillator two-point functions, so its higher free correlators agree by the Wick theorem. This establishes the equivalence for the regulated free theory and order by order in the perturbation series.

To couple the matter field electromagnetically, promote its global phase symmetry to the local transformationThe gauge covariant derivative $D_a={\partial }_a+ie{A}_a$ obeys $D_a^{\prime }{\psi }^{\prime }=e^{-ie\lambda }{D}_a\psi$. Replacing ${\partial }_a$ by $D_a$ in the Dirac action therefore gives the invariant matter densityThe electromagnetic field tensor is unchanged by the local transformation. Thus the unfixed quantum electrodynamics action is gauge-invariant. The Dirac current $j^a=\bar{\psi }{\gamma }^a\psi$ is conserved by the matter equations, and the electron and Positron excitations carry opposite charges. The added gauge fixing density selects a representative and is not itself invariant under arbitrary local transformations; it does not change gauge-invariant observables. At a free fermion vertex,For external on-shell Dirac spinors the corresponding current contraction vanishes. The quantum extension is the Ward identity, which makes physical amplitudes insensitive to adding a multiple of the photon momentum to its polarization vector. A gauge-compatible regularization preserves this vector-current identity.

With ${\mathcal{L}}_{\mathrm{int}}=-e\bar{\psi }{\gamma }^a\psi A_a$, expand $\exp (i\int {\mathcal{L}}_{\mathrm{int}})$ in powers of $e$. A term of order $n$ contains $n$ spacetime integrations and $1/n!$. Applying the Wick theorem pairs the free fields: an $A$-$A$ Wick contraction supplies a photon propagator, and a $\psi$-$\bar{\psi }$ Wick contraction supplies an oriented Dirac propagator. Each insertion supplies an interaction vertex. The permutations of contractions cancel the expansion factorials except for the Feynman-diagram symmetry factor. Interchanging Grassmann fields produces the fermionic sign, including a minus sign for every closed fermion loop. This is how Feynman diagrams arise from expectation values, rather than an extra dynamical assumption.

The normalized vacuum generating functional removes components with no external insertions. In the operator formalism, for an interacting-vacuum expectation value of an inserted product $\mathcal{O}$, the same cancellation appears asVacuum projection and the Feynman i-epsilon prescription are implicit. The linked-cluster theorem exponentiates all connected vacuum bubbles into the same factor in numerator and denominator, so it cancels. Normalize by $Z[0]$ to remove every vacuum component. This cancellation of vacuum bubbles still leaves products of disconnected diagrams that each contain external insertions. If only connected correlators are wanted, differentiate the connected generating functional $W=-i\log [Z[J,\bar{\eta },\eta ]/Z[0]]$.

The resulting momentum-space QED Feynman rules for the bare theory can be stated in Feynman gauge as follows:

The illustrated interaction vertex has incoming fermion momentum $p$, incoming photon momentum $k$, and outgoing fermion momentum $p^{\prime }=p+k$. Its solid-line arrows indicate fermion flow. The photon wavy line has no fermion-flow arrow. The propagators, the vertex $-ie{\gamma }^a$, momentum conservation, and the fermionic signs determine the perturbative amplitudes.

For this question use the mostly-plus Dirac convention: the Minkowski metric is ${\eta }_{ab}=\mathrm{diag}(-1,1,1,1)$ andThis convention matches the printed plane-wave phase and final identity. It is related to the usual mostly-minus Dirac equation by reversing the metric and taking the negatives of the usual gamma matrices. In particular, the resulting Dirac action and its Dirac adjoint describe the same physical massive field. The Dirac gamma matrices are four complex $4\times 4$ matrices representing the spacetime Clifford algebra with quadratic form $-\eta$. The irreducible complex representation has dimension four. A convenient explicit choice is the negative of the standard Dirac representation of the gamma matrices:where ${\sigma }^i$ are the Pauli matrices. Their multiplication law verifies the displayed anticommutator.

In this representation the gamma matrix adjoint and transpose identities areandThe invariant way to express the latter pattern uses the charge-conjugation matrix:Individual transpose signs depend on the basis. More generally, a similarity transformation ${\gamma }^{\prime a}=M{\gamma }^a{M}^{-1}$ changes the Hermitizing matrix to $H^{\prime }=M^{-†}{\gamma }^0{M}^{-1}$ and the charge-conjugation matrix to $C^{\prime }=MC{M}^T$. Then ${\gamma }^{\prime a†}=H^{\prime }{\gamma }^{\prime a}{H}^{\prime -1}$ and $C^{\prime -1}{\gamma }^{\prime a}{C}^{\prime }=-{\gamma }^{\prime aT}$. Thus the simple formula with ${\gamma }^{\prime 0}$ itself presumes a compatible Hermitian basis, rather than an arbitrary nonunitary similarity transformation.

Applying $(i{\gamma }^a{\partial }_a-m)$ to the Dirac equation gives $(\square -m^2)\Psi =0$. Its mass shell is $p^2=-m^2$, so the frequencies are $p^0=\pm \sqrt{{\mathbf{p}}^2+m^2}$. The Dirac spinor transforms in the four-component Spinor representation of the Lorentz group. Under spatial rotations, the two upper and the two lower components each transform as a two-component spin-$1/2$ representation: the spin angular momentum matrices are $\mathrm{diag}({\sigma }^i,{\sigma }^i)/2$. At rest the positive-energy equation selects the upper two components, giving two independent spin polarizations, and the negative-frequency equation selects the lower two.

In the quantum theory a mode expansion isWith the mostly-plus Dirac convention, $e^{ip\cdot x}=e^{-iEt+i\mathbf{p}\cdot \mathbf{x}}$ is positive frequency. The negative-frequency coefficient obeys $(/p+m)v_s(p)=0$. The fermionic annihilation operators $b_s$ and $d_s$ satisfy the canonical anticommutation relations, with their respective fermionic creation operators. The $b_s^{†}$ excitations are particles; the $d_s^{†}$ excitations are antiparticles with the same positive energy, mass and spin-$1/2$ but opposite charge. After normal ordering, the Hamiltonian operator contains positive multiples of $b_s^{†}{b}_s+d_s^{†}{d}_s$. Reinterpreting the negative-frequency part as antiparticle creation supplies a spectrum bounded below rather than a physical tower of negative-energy particles.

For the printed $e^{ip\cdot x}$ wave, ${\partial }_a\Psi =i{p}_a\Psi$. Substitution givesThe spin label $s$ indexes the two states of a spin-$1/2$ particle, rather than varying the particle's total spin. For real on-shell $p$, Hermitian conjugation and ${\gamma }^{a†}={\gamma }^0{\gamma }^a{\gamma }^0$ giveThese are right and left null-vector equations for the same on-shell matrix.

To obtain the Gordon identity, take both external Dirac spinors to have the same real mass $m$. Their two equations implyDefine ${\gamma }^{ab}=[{\gamma }^a,{\gamma }^b]/2$. The Clifford algebra relation yieldsThereforeMultiplying the previous null-vector equation by $-1$ proves the required formula exactly:For $m\ne 0$ it can be solved for the vector-current matrix element, separating a momentum term from the antisymmetric Dirac spinor term. The identity before division also holds at $m=0$. The signs depend jointly on the metric, Clifford relation, Dirac mass term and plane-wave phase. In the mostly-minus convention of Questions 1 and 3, the printed phase instead gives $(/p+m)u=0$; the corresponding identity uses $(p_b-p_b^{\prime })$ in the antisymmetric term. Mixing that convention with the formula proved here would produce an apparent sign error.

Use ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =1$. The free real scalar field has actionThe canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Hamiltonian operator is obtained by the Legendre transform of the density:Canonical quantization promotes $\varphi$ and $\pi$ to Hermitian operator-valued fields and imposes the equal-time canonical commutation relationsTheir Heisenberg equation of motion gives $\dot{\varphi }=\pi$, $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, hence the Klein-Gordon equation $(\square +m^2)\varphi =0$.

Let $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$ and $p\cdot x=E_{\mathbf{p}}t-\mathbf{p}\cdot \mathbf{x}$. The Hermitian field mode expansion isThe scalar field oscillator inversion extractsand its adjoint obtained by Hermitian conjugation extracts $a^{†}$. Substitute these expressions into the equal-time canonical commutation relation. The two mixed field-momentum terms giveFor $[a(\mathbf{p}),a(\mathbf{q})]$, the corresponding coefficient is a difference of those square roots and multiplies ${\delta }^3(\mathbf{p}+\mathbf{q})$; it vanishes because $E_{\mathbf{p}}=E_{-\mathbf{p}}$. Taking the adjoint gives $[a^{†},a^{†}]=0$. Thus these are bosonic annihilation operators and creation operators. The vacuum satisfies $a(\mathbf{p})|0⟩=0$, and normal ordering gives $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}(\mathbf{p})a(\mathbf{p})/(2\pi {)}^3$, after removing the constant zero-point energy. A one-particle excitation has energy $E_{\mathbf{p}}$ and spin zero.

The Feynman propagator $G_F(x-y)=⟨0|T\varphi (x)\varphi (y)|0⟩$ is the vacuum expectation value of the time-ordered product of two field insertions. For $x^0>y^0$, the field at $y$ creates a one-particle excitation from the vacuum and the field at $x$ annihilates it; the opposite time ordering reverses the roles. It is a propagation amplitude and correlation function of vacuum fluctuations, rather than a transition probability. From the oscillator expansion, only the $a$-$a^{†}$ contraction survives, so with $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$,Here $\theta$ is the Heaviside step function, and changing $\mathbf{p}$ to $-\mathbf{p}$ in the second term gives the last line.

For the Fourier transform convention ${\tilde{G}}_F(\omega ,\mathbf{p})=\int dt{d}^{3}r{e}^{i\omega t-i\mathbf{p}\cdot \mathbf{r}}{G}_F(t,\mathbf{r})$, insert a positive damping factor $e^{-ϵ|t|}$ and integrate the positive and negative half-lines separately:The last equality is a distribution identity; the infinitesimals in the partial fractions need not have the same finite magnitude as the infinitesimal in the combined denominator. The Feynman i-epsilon prescription means that the positive-frequency pole is just below the real energy axis and the negative-frequency pole is just above it. Closing the contour integral below for $t>0$, and above for $t<0$, reproduces the oscillator result by the residue theorem. This prescription fixes which homogeneous solutions are added to the Green function. As an independent normalization check, $e^{-iE|t|}/(2E)$ has derivative jump $-i$ at $t=0$, soIt is the expectation-value convention for $G_F$, including the numerator $i$, that determines this source normalization.

The vacuum is a centered free Gaussian state. The Wick theorem expresses the four-field time-ordered product as the sum of all pair Wick contractions plus terms containing a normal-ordered product. The latter terms have zero vacuum expectation value. There are exactly three complete pairings, with no fermionic signs for this bosonic field. Thus the free scalar four-point function isThese formulas describe canonical quantization of a real scalar field; in particular, a real scalar field uses a single oscillator family rather than independent charged-particle and antiparticle families.

Return to the mostly-minus Minkowski metric of Question 3, and write $f={\partial }_a{A}^a$. The field $\xi (x)$ is assumed real and nonzero wherever its inverse occurs. Under the usual Abelian gauge transformation $A_a\mapsto A_a+{\partial }_a\lambda$, with $\xi$ inert, the electromagnetic field tensor is invariant but $f\mapsto f+\square \lambda$. The change in the added density isThis is not generally a total derivative. The action is not invariant under arbitrary gauge transformations. It retains residual Lorenz gauge symmetry for transformations satisfying $\square \lambda =0$, with the same boundary conditions imposed before and after the transformation. The action itself is undefined at $\xi =0$; that value can only be considered as a limiting gauge.

To vary the electromagnetic four-potential, use $\delta F_{ab}={\partial }_a\delta A_b-{\partial }_b\delta A_a$ and $\delta f={\partial }_a\delta A^a$. The integration by parts of both terms giveswith surface terms removed by the variational boundary conditions. HenceFor variable $\xi$, the derivative must act on $1/\xi$ as well as on $f$:Since $\xi$ enters algebraically, its Euler-Lagrange equation isFor a real field this gives $f=0$. Thus a dynamical gauge-fixing parameter imposes a constraint equation in field theory; it is not an ordinary propagating scalar. On configurations satisfying both equations, ${\partial }_a{F}^{ab}=0$, ${\partial }_a{A}^a=0$, and therefore $\square A^b=0$. There is no independent kinetic equation that determines $\xi$.

For the momentum-space equation at a prescribed constant $\xi \ne 0$, use $A_a(x)=\int d^{4}p{e}^{-ip\cdot x}{A}_a(p)/(2\pi {)}^4$. Then ${\partial }_a\mapsto -i{p}_a$, and the linear $A$ equation isEquivalently, $-p^2{A}^a+(1+1/\xi )p^a(p\cdot A)=0$. If the separately varied $\xi$ equation is also imposed on a real classical solution, then $p\cdot A=0$ and its nonzero modes lie on the massless mass shell. For constructing a full-field Green function, however, invert the fixed-background quadratic operator before imposing that on-shell constraint. Holding $\xi$ fixed and integrating over $A$ is a Gaussian functional integral; integrating over the original variable $\xi$ as well is a different constrained problem.

For $p^2\ne 0$, the Lorentzian versions of the transverse projector of a vector field and longitudinal projector of a vector field areThey satisfy $P_T^2=P_T$, $P_L^2=P_L$, $P_T{P}_L=0$ and $P_T+P_L=I$. The mixed-index operator isIts inverse follows by inverting these two scalar eigenvalues. Define $G_{ab}$ by $K^{ac}{G}_{cb}={\delta }_b^a$. The positive-sign covariant gauge-fixing inverse isThis is the algebraic inverse of the kinetic operator. A vacuum photon propagator instead has the factor $i{G}_{ab}$ and a Feynman i-epsilon prescription for its simple and double poles. If the name $G$ is used for the vacuum two-point function, include that factor $i$ consistently in its defining source equation. A retarded Green function would use different boundary conditions at the same poles. The original plus sign corresponds to the usual covariant gauge parameter $\alpha =-\xi$; in particular, $\xi =-1$ gives Feynman gauge and $G_{ab}=-{\eta }_{ab}/p^2$.

Contracting gives the longitudinal gauge propagator contractionThus the contraction is not identically zero as a function of four-momentum for any $\xi \ne 0$, and it is a nonzero vector at every nonzero non-null four-momentum. A particular component may vanish when $p_b=0$. At $p^2=0$, the displayed inverse has poles and must be interpreted using the chosen Green function prescription, rather than as a pointwise finite matrix. In the limiting covariant Landau gauge, $\xi \to 0$, the longitudinal part vanishes. A covariant photon Green function may have a longitudinal component even though physical photon polarizations are transverse.