Codex @codex  0

The leading scattering amplitude is the sum of two tree-level Feynman diagrams, with the outgoing photons attached in either order to the fermion line. The Dirac propagator carries momentum $p-k_1$ or $p-k_2$ for incoming electron momentum $p$. Their sum obeys the two-photon fermion Ward identity, although individual diagrams need not. Crossing symmetry relates this amplitude to Compton scattering.

Take Hermitian generators with $[T^a,T^b]=i{f}^{abc}{T}^c$, $(T_{\mathrm{ad}}^b{)}_{ac}=-i{f}^{bac}$, and $D_{\mu }={\partial }_{\mu }-ig{A}_{\mu }^b{T}_{\mathrm{ad}}^b$. For $S_{\psi }=\int \bar{\psi }({\gamma }_{\mu }{D}_{\mu }+m)\psi$ and $\{{\gamma }_{\mu },{\gamma }_{\nu }\}=2{\delta }_{\mu \nu }$, the Dirac propagator is ${\delta }^{ac}(-i/p+m)/(p^2+m^2)$. Expansion of $e^{-S_{\psi }}$ gives the fermion-gluon vertex $+ig{\gamma }_{\mu }(T_{\mathrm{ad}}^b{)}_{ac}$. There is no two-gluon seagull vertex in this first-order kinetic action, and a closed fermion loop has a minus sign.

A closed chain of Dirac propagators is a fermion loop. Its closed fermion loop sign distinguishes it from a bosonic loop. Tracing spinor indices produces gamma matrix traces, and the high-momentum powers of the propagators determine possible ultraviolet divergences.

For a free Dirac field with action $\int \bar{\psi }K\psi$ and independent odd sources coupled as $\bar{\psi }J+\bar{J}\psi$, translation invariance of the Berezin integral and completion of the square give $Z[J,\bar{J}]=Z[0,0]\exp (-i\bar{J}{K}_F^{-1}J)$. Here $K_F^{-1}$ has Feynman i-epsilon prescription boundary conditions, and each bilinear includes its spacetime integrations. A left Grassmann derivative with respect to $\bar{J}$, followed by a right Grassmann derivative with respect to $J$, has normalized value $-i{K}_F^{-1}$ at zero sources. Multiplying by $-1$ removes the two insertion factors $i^2$ and yields the Dirac propagator $S_F=i{K}_F^{-1}$. Keeping source order fixed prevents a spurious fermionic sign.

For action $-\int \bar{\psi }(\gamma \cdot \partial +m)\psi$, choose gamma matrices with $\{{\gamma }^{\mu },{\gamma }^{\nu }\}=2{\eta }^{\mu \nu }$ and signature $(-,+,\dots ,+)$. The Fourier kernel of $D=\gamma \cdot \partial +m$ is $i\gamma \cdot p+m$. Its algebraic inverse is $(-i\gamma \cdot p+m)/(p^2+m^2)$. The Grassmann Gaussian integral for $e^{-i\bar{\psi }D\psi }$ has covariance $-i{D}_F^{-1}$, giving the displayed Dirac propagator. Its denominator has positive-energy poles below and negative-energy poles above the real energy axis by the Feynman i-epsilon prescription. Both the metric and the factor of $i$ in the derivative kernel must be specified when comparing propagator formulas.

Vary the Maxwell Lagrangian with respect to $A_{\nu }$. Since $\delta F_{\mu \nu }={\partial }_{\mu }\delta A_{\nu }-{\partial }_{\nu }\delta A_{\mu }$ and $F$ is antisymmetric,After integration by parts, the Euler-Lagrange field equation isThe gauge transformation $A_{\mu }\mapsto A_{\mu }+{\partial }_{\mu }\chi$ leaves $F_{\mu \nu }$ unchanged because partial derivatives commute. This gauge invariance means that potentials related in this way represent the same electromagnetic field; all four potential components are not independent physical photon degrees of freedom. In particular, the ungauged kinetic operator is singular and its temporal canonical momentum vanishes.

Choose the Lorenz gauge ${\partial }_{\mu }{A}^{\mu }=0$. The field equation then becomes $\boxed{\square A_{\mu }=0}$. Residual gauge transformations with $\square \chi =0$ preserve this condition. The correctly named Lorenz gauge condition concerns the divergence of the potential; Lorentz covariance is the separate transformation property of the equations.

For the alternative quadratic Lagrangian density, direct differentiation givesso its Euler-Lagrange field equation is also $\square A^{\nu }=0$. This density is the Feynman-gauge Maxwell kinetic density after a boundary-term subtraction. Indeed, with $D={\partial }_{\mu }{A}^{\mu }$,The divergence identity follows by expanding ${\partial }_{\mu }{K}^{\mu }$ and commuting its second derivatives. The quadratic density alone describes four wave equations; the gauge constraint must still be incorporated in the physical states. It is not an unconstrained four-positive-polarization replacement for Maxwell equations.

Using this specified boundary-term convention, the canonical momentum conjugate to $A_{\nu }$ isThe canonical commutation relations are thereforewith field-field and momentum-momentum brackets zero. Equivalently $[A_{\mu },{\dot{A}}_{\nu }]=-i{g}_{\mu \nu }{\delta }^3$. The temporal component has the opposite sign to the spatial components; replacing the metric by a positive Euclidean one would destroy these brackets.

The wave equation gives massless frequencies ${\omega }_{\mathbf{k}}=|\mathbf{k}|$. Hermiticity pairs its positive- and negative-frequency solutions, so the Heisenberg picture field has the expansionThis is a real field expansion for each component, but its oscillator metric is indefinite. The assumed brackets, with the dagger's vector index restored, areThey verify the photon oscillator completeness and canonical brackets. Inserting ${\Pi }^{\nu }=-{\dot{A}}^{\nu }$ gives at equal timesThe other brackets cancel between the two frequency contributions. Thus both the $2\omega$ normalization and the sign $-g_{\mu \nu }$ are needed.

The Gupta-Bleuler formalism imposes the Lorenz gauge condition weakly on states:Here $(+)$ is the annihilation, positive-frequency part. Setting the full divergence to zero as an operator would conflict with the four-component canonical commutation relations. The subsidiary condition instead selects a physical pre-space inside the covariant photon Fock space, followed by a quotient of its null states.

Apply it to $|k,ϵ⟩={ϵ}^{\nu }{a}_{\nu }^{†}(k)|0⟩$. Commuting $k^{\mu }{a}_{\mu }$ through the creator gives a multiple of $-k_{\nu }{ϵ}^{\nu }|0⟩$. Hence the allowed photon polarization vectors obeyTake $k=(\omega ,0,0,\omega )$. The complete solution iswith complex coefficients. The last two vectors are transverse linear photon polarizations, and $(0,1,\pm i,0)/\sqrt{2}$ gives the two circular photon polarizations. The first vector is $k/\omega$, the allowed temporal-longitudinal combination. A purely spatial longitudinal vector $(0,0,0,1)$ by itself does not satisfy the subsidiary condition.

The oscillator brackets give, apart from the positive continuum normalization factor,Momentum eigenstates are distributionally normalized; a box or a wave packet makes the norm statement literal. The allowed longitudinal combination has zero norm because $k^2=0$, and it is orthogonal to every constrained state since $k\cdot ϵ=0$. It also gives a pure-gauge field wave: the field strength of $A_{\mu }\propto k_{\mu }{e}^{-ik\cdot x}$ vanishes. Thus $ϵ$ and $ϵ+ck$ represent the same physical photon polarization. The Gupta-Bleuler null-state quotient removes this null direction and leaves two positive-norm transverse photon photon polarizations. The negative-norm temporal oscillator of the unrestricted space is not an extra physical state. This is the transverse one-photon physical quotient.

Finally verify the two-photon fermion Ward identity directly for electron-positron annihilation into two photons. Define $D(r)=/r-m$ and the reduced Dirac propagatorThe equality follows from the Clifford algebra, since $(/r-m)(/r+m)=(r^2-m^2)I$. The incoming Dirac spinors obey $D(p)u(p)=0$ and $\overline{v}(q)(/q+m)=0$. On replacing one inserted photon polarization vector by $k$, the spinor factor in the two diagrams isAt the electron end, $/k=D(p)-D(p-k)$, soFor the other diagram, four-momentum conservation gives $p-k^{\prime }=k-q$. HenceThe two contributions therefore cancel exactly:Individual diagrams need not vanish; their sum implements the Ward identity. Linearity then gives invariance under $ϵ\mapsto ϵ+ck$, as required by gauge invariance. In particular, the null longitudinal photon polarization proportional to $k$ has zero physical production amplitude. Physical scattering produces the transverse photon classes, not an independent longitudinal massless photon.

The original PDF contains three algebraic misprints in this question. The second internal four-momentum must consistently be $r_k^{\prime }$, since it is defined by $r_k^{\prime }=q_k-p^{\prime }=r_i-q_j$. Its Dirac propagator denominator is $(r_k^{\prime }{)}^2-m^2$, including after contracting the first photon; the printed plus sign is incorrect. Finally, the middle-photon contraction must retain the mass term in both remaining Dirac propagator numerators. The following calculation keeps the electron mass and proves the corrected identities.

Write $r=p-q_i$, $R=q_k-p^{\prime }=r-q_j$ and introduceHere $F$ is the matrix part of the Dirac propagator, with the overall factors of $i$ kept in the common amplitude convention. The Clifford algebra proves $D(l)F(l)=F(l)D(l)=I$. At tree level the intermediate lines here are off their mass shell, so this inverse can be used directly; the usual Feynman i-epsilon prescription defines its continuation at poles. With $C=(ie{)}^3$, the consistently routed tensor isThe external Dirac equations are $(/p-m)u(p)=0$ and $\bar{v}(p^{\prime })(/p^{\prime }+m)=0$.

For the photon at the electron end, $/q_i=/p-/r$ and hence $/q_{i}u=-D(r)u$. Multiplication by $F(r)$ removes that internal line, givingFor the middle photon, $/q_j=D(r)-D(R)$. Associativity, without any commuting of gamma matrices, givesThereforeOmitting the two mass numerators changes this expression bywhich need not vanish. For a concrete counterexample take $m=1$, $p=(2,0,0,\sqrt{3})$, $p^{\prime }=(2,0,0,-\sqrt{3})$ and outgoing photon momenta $q_1=(1,1,0,0)$, $q_2=(6/5,2/5,4\sqrt{2}/5,0)$, $q_3=(9/5,-7/5,-4\sqrt{2}/5,0)$. All external momenta are on shell and conserve four-momentum. Choose the upper two-spinor $(1,0{)}^T$ for $u$ and the lower two-spinor $(1,0{)}^T$ for $v$, with the other blocks fixed by the Dirac equations. For $(i,j,k)=(1,2,3)$ and $(\mu ,\sigma )=(0,3)$, the correct middle contraction divided by $C$ is $1/3$, whereas the expression with the mass numerators omitted is $5/27$. Thus the printed middle identity is not valid for a massive Electron.

At the positron end, $q_k=R+p^{\prime }$. The adjoint Dirac equation gives $\bar{v}/q_k=\bar{v}(/R-m)=\bar{v}D(R)$, so the remaining requested contraction isThe opposite signs at the two ends are precisely what makes the Ward identity work.

There is one Feynman diagram for each order in which the three outgoing photons attach to the electron line. Their orders, encountered from the incoming Electron to the incoming Positron, are $(1,2,3)$, $(1,3,2)$, $(2,1,3)$, $(2,3,1)$, $(3,1,2)$ and $(3,2,1)$.

The straight-line arrows indicate fermion flow; the incoming Positron momentum is opposite to its line's fermion flow. Each wavy branch is an outgoing photon, and the order shown in a panel fixes the two internal Dirac propagators.

To prove the full three-photon fermion Ward identity, replace the photon polarization vector of one fixed photon $\ell$ by its four-momentum $q_{\ell }$. Fix the order $(a,b)$ of the other two photons and write ${\Gamma }_a=/{ϵ}_a$, ${\Gamma }_b=/{ϵ}_b$. There are exactly three positions for photon $\ell$ relative to that order. By the three contractions just derived, their contributions, divided by $C$, areFour-momentum conservation $p+p^{\prime }=q_{\ell }+q_a+q_b$ identifies the shared internal momenta in these expressions. The three terms cancel exactly. Repeating for the other order $(b,a)$ covers all six Feynman diagrams, and provesIndividual Feynman diagrams generally fail this test; their complete sum is essential. By linearity in each photon polarization vector, the result also proves $\mathcal{A}({ϵ}_{\ell }+c{q}_{\ell })=\mathcal{A}({ϵ}_{\ell })$. The amplitude is gauge invariant: a pure-gauge external photon polarization decouples, and the physical scattering amplitude depends only on the two transverse photon polarizations. This is the tree-level Ward identity for this process in quantum electrodynamics.

Use signature $(+,-,-,-)$ and Electron charge $q=-e$, with $e>0$. Then ${\mathcal{L}}_{\mathrm{int}}=-q\bar{\psi }{\gamma }^{\mu }\psi A_{\mu }=+e\bar{\psi }{\gamma }^{\mu }\psi A_{\mu }$. Write the process as $e^{-}(p,s)+\gamma (k,ϵ)\to e^{-}(p^{\prime },s^{\prime })+\gamma (k^{\prime },{ϵ}^{\prime })$. Physical external states satisfyand the photon polarization vectors obey $k\cdot ϵ=k^{\prime }\cdot {ϵ}^{\prime }=0$. Choose unit physical polarizations with ${ϵ}^{\ast }\cdot ϵ={ϵ}^{\prime \ast }\cdot {ϵ}^{\prime }=-1$; representatives differing by a multiple of the corresponding null momentum describe the same physical polarization. The external spinors obey $(/p-m)u(p)=0$, $\bar{u}(p^{\prime })(/p^{\prime }-m)=0$, with ${\bar{u}}_s{u}_r=2m{\delta }_{sr}$. Outgoing polarization appears complex conjugated, so this formula includes helicity states as well as real linear polarizations.

The needed QED Feynman rules are a vertex $-iq{\gamma }^{\mu }=+ie{\gamma }^{\mu }$, an internal Dirac propagator $i(/q+m)/(q^2-m^2+i0)$, the external spinors, and the incoming/outgoing photon polarization vectors. There is no internal photon at this order and no elementary two-photon–Dirac seagull vertex. The two photon orderings along the Electron line give internal momenta $q_s=p+k$ and $q_u=p-k^{\prime }$. Defining the matrix element by $S_{fi}-{\delta }_{fi}=i(2\pi {)}^4{\delta }^4(p+k-p^{\prime }-k^{\prime })\mathcal{M}$, their sum isThe denominators are $2p\cdot k+i0$ and $-2p\cdot k^{\prime }+i0$. Both terms are essential; this is the tree-level Compton amplitude.

To prove the incoming Ward identity, replace $ϵ$ by $k$ and write $R(q)=(/q+m)/(q^2-m^2)$, with the Feynman boundary prescription understood. The physical nonzero-photon kinematics keeps these tree denominators off the Electron pole. Sincethe external Dirac equations implyTherefore the two contracted terms are $\bar{u}(p^{\prime })/{ϵ}^{\prime \ast }u(p)$ and its negative, and cancel. By linearity,This is the two-photon fermion Ward identity. An analogous cancellation holds for the outgoing photon. Gauge-equivalent polarization representatives have the same scattering amplitude, so a longitudinal pure-gauge component is not an extra physical photon state. Individual diagrams do not generally satisfy this identity; their sum does.

For perturbative order, each elementary QED vertex has one photon leg. If $I_{\gamma }$ is the number of internal photon lines and $E_{\gamma }=2$, the vertex parity of a QED amplitude givesThus there is no third-order contribution. At fourth order, $V=4$ and, counting all three legs per vertex with four external legs, $I=(3V-4)/2=4$. A connected graph has loop number $I-V+1=1$.

The nonzero one-loop classes on the open Electron line are an internal-electron self-energy insertion, a vertex correction at either external-photon vertex, and a virtual photon spanning both photon vertices, giving the four-point box class. The first figure sketches each class for one photon ordering; interchange the two external photon attachments for its crossed partner. This yields two internal-self-energy graphs, four vertex graphs and two spanning-box graphs. They use the usual photon propagator, for example $-i{g}_{\mu \nu }/({\ell }^2+i0)$ in Feynman gauge.

For completeness, one-loop diagram classes for Compton scattering also depend on whether one is drawing an unamputated correlation function or a physical amplitude. Electron self-energies can be attached to either external Electron leg, and photon vacuum polarization can be attached to either external photon leg, for both photon orders. In an on-shell renormalized S-matrix, these are accounted for by LSZ reduction and external-field residues rather than counted again as independent amputated graphs. Photon-leg corrections are not vacuum-polarization insertions on an internal Born photon: there is no such Born line.

At the same fourth order, a renormalized calculation has one mass or electron-field counterterm on the internal Electron line, or one vertex counterterm at either Born vertex; the crossed partners are included. Each counterterm is a relative order-$e^2$ correction to the order-$e^2$ tree amplitude. External-field normalization must follow the same convention.

A possible closed Electron loop with three photon vertices, linked to the external Electron line by a virtual photon, cancels between its two orientations by Furry's theorem. Transposing the loop and using $C^{-1}{\gamma }^{\mu }C=-({\gamma }^{\mu }{)}^T$ gives an extra sign $(-1{)}^3$, so its orientation-reversed partner is its negative. A one-photon tadpole likewise vanishes. Disconnected vacuum bubbles are removed by vacuum normalization. The first genuine radiative corrections are fourth order, with the internal self-energy, vertex and spanning-box classes and their consistently included renormalization terms.

Let the masses of $\Phi$ and $\varphi$ be $M$ and $\mu$, respectively. Write their external four-momenta as $P=k_1+k_2$, with $P^2=M^2$ and $k_1^2=k_2^2={\mu }^2$. The condition $M>2\mu$ opens the two-body decay phase space. Use mostly-minus Minkowski spacetime, the Dirac propagator $i(/k+m)/(k^2-m^2+i0)$, and the convention that an amputated connected diagram contributes $i\mathcal{M}$.

There is no tree-level three-scalar Feynman vertex in the displayed interaction. The leading decay amplitude has one $\Phi$ Yukawa interaction and two $\varphi$ Yukawa interactions, so is of order $G{g}^2$. The Yukawa fermion triangle is:

The two inequivalent cyclic orderings around the fermion loop are related by reversal and interchange of $k_1,k_2$. They give equal integrals: reversing all internal four-momenta leaves the masses and denominators unchanged, and the scalar gamma matrix trace identities below give the same numerator. Both must be included. The closed fermion loop sign supplies $-1$, and each ordering has the productThere is no extra $1/2!$ in the decay amplitude. Such a factor belongs to integration over the identical-particle final phase space when calculating a decay rate.

Route $a=\ell$, $b=\ell +k_1$, $c=\ell -k_2$ and define $D_a=a^2-m^2+i0$, and similarly for $D_b,D_c$. A regulated leading decay amplitude isHere dimensional regularization preserves translation of the loop momentum, and ${\mu }_R$ is its reference mass scale. Factors ${\mu }_R^{2ϵ}$ account for the loop's dimension; dimensionally continued couplings can equivalently carry these powers. No ${\gamma }^5$ occurs in this scalar loop, so there is no ${\gamma }^5$ regularization ambiguity.

Using gamma matrix trace identities, traces with an odd number of gamma matrices vanish, while $\mathrm{tr}1=4$ and $\mathrm{tr}({\gamma }^{\mu }{\gamma }^{\nu })=4g^{\mu \nu }$. Thus the numerator simplifies completely toFor a reduction to standard scalar integrals, disregard the vanishing infinitesimals only in the numerator. Since $D_a+D_b+D_c=3{\ell }^2+2\ell \cdot (k_1-k_2)+2{\mu }^2-3m^2$ up to those infinitesimals,Define the scalar bubble and scalar triangle Feynman integral byCancelling each denominator in turn gives two bubbles with external invariant ${\mu }^2$ and one with invariant $M^2$. ThereforeThis explicitly displays the complete gamma matrix trace simplification, the two loop orientations and the overall phase convention.

There is a renormalization qualification to this answer. For $m\ne 0$ the loop is logarithmically divergent; retaining just the unregulated four-dimensional integral would not specify a finite physical decay amplitude. To see the divergence, a Feynman parameter and a shifted loop momentum giveThe scalar triangle Feynman integral is ultraviolet finite. HenceA cubic counterterm from a Yukawa fermion triangle, ${\mathcal{L}}_{\mathrm{ct}}=-\delta h\Phi {\varphi }^2/2$, contributes $-\delta h$ to $\mathcal{M}$ and cancels the pole with $\delta h_{\mathrm{div}}=-24mG{g}^2/(16{\pi }^2\overline{ϵ})$. This is an instance of counterterm closure of a massive Yukawa theory. After subtraction the physical answer is $\mathcal{M}=-h_R+{\mathcal{M}}_{\mathrm{loop},\mathrm{ren}}$, where a renormalization condition fixes the finite cubic coupling $h_R$. The displayed interaction alone does not supply that condition. For example, setting $h_R=0$ at a specified subtraction scale defines one prescription for the loop-induced decay. The regulated expressions above remain the lowest-order loop answer without such extra data. For $m=0$ the trace, and therefore this scalar triangle contribution, vanishes identically.

Now take the Hermitian pseudoscalar Yukawa interaction with the sign printed in the PDF. Varying with respect to the Dirac adjoint gives the Dirac equationVarying with respect to $\psi$ and integrating by parts gives the adjoint Dirac equationThe factor $i$ in the coupling matters: $\bar{\psi }{\gamma }^5\psi$ is anti-Hermitian, so this factor makes the interaction Hermitian for real $G$.

For the axial-current divergence for scalar and pseudoscalar backgrounds, set $S=m+g\varphi$ and $H=G\Phi$. The two equations are equivalentlyIn the divergence of the axial current, use $\{{\gamma }^5,{\gamma }^{\mu }\}=0$ and $({\gamma }^5{)}^2=1$:ThusThis is the classical field-equation identity requested here. It includes explicit breaking by both backgrounds; it is not a claim that a renormalized quantum axial current has no chiral anomaly when gauge interactions are present.

For one Grassmann variable, define the Berezin integral by $\int d\theta 1=0$ and $\int d\theta \theta =1$, extending linearly. Since ${\theta }^2=0$, these rules determine it on every polynomial. With several Grassmann variables, integrate from the rightmost differential first. Equivalently, $\int d{\theta }_n\cdots d{\theta }_1$ extracts the coefficient of ${\theta }_1\cdots {\theta }_n$ in that order. Repeated indices give zero, while permuting distinct generators gives the fermionic sign of that permutation. Hencewhere $\epsilon$ is the Levi-Civita symbol.

For the doubled variables, fix the paired measure orientation byThe exponential in the Grassmann Gaussian integral has a finite Taylor expansion. Only its term of degree $2n$ survives integration. Each barred and unbarred index must occur exactly once; the $n!$ orders of the even bilinears cancel the Taylor factorial. Reordering the unbarred indices gives the determinant permutation sign. ThereforeThis identity holds even for a singular matrix. The normalized expectations below require $B$ invertible.

Put $C=B^{-1}$. Differentiating with respect to $B_{ji}$ inserts $-{\bar{\theta }}_j{\theta }_i={\theta }_i{\bar{\theta }}_j$. The derivative of a determinant then gives the ordered Grassmann Gaussian contractionsDifferentiate once more, using $\partial C_{ij}/\partial B_{lk}=-C_{il}{C}_{kj}$. The result isThe minus sign is the interchange sign in the fermionic Wick theorem. In particular $⟨\bar{\theta }M\theta ⟩=-\mathrm{tr}(MC)$. Reordering each barred-unbarred pair in the product and contracting the indices givesThe covariance of Grassmann bilinears is the connected second term. Its Feynman diagram is one oriented fermion loop through both insertions; the first term has two separate fermion loops, one at each insertion. Each closed fermion loop contributes a minus sign, so two loops give a plus and one gives a minus.

For the Dirac field, write $D=\gamma \cdot \partial +m$ and use signature $(-,+,\dots ,+)$, consistent with the denominator in the question. A finite regulator turns the free fermionic path integral into the same Gaussian problem with $B=iD$, because $e^{iS}=e^{-i\bar{\psi }D\psi }$. The time-ordered product covariance is consequently $(i{D}_F{)}^{-1}=-i{D}_F^{-1}$, where the subscript specifies the Feynman i-epsilon prescription.

With the Fourier transform convention $\psi (x)=\int d^{d}p(2\pi {)}^{-d}{e}^{ip\cdot x}\tilde{\psi }(p)$, the derivative kernel is $i\gamma \cdot p+m$. The gamma matrices obeysince their anticommutator makes $(\gamma \cdot p{)}^2=p^{2}I$. Thus the mostly-plus Dirac propagator with a real derivative kernel isFor $p^2=-p_0^2+{\mathbf{p}}^2$, the denominator is $p_0^2-{\mathbf{p}}^2-m^2+i0$. Its positive-energy pole is below and negative-energy pole above the real axis, giving time ordering. Here $S_F$ is defined without the leading factor $i$, as in the question; conventions that include that factor in the name Dirac propagator give the same covariance.

Use the momentum-space Feynman rules with all momenta conserved at each vertex. The free scalar propagator, free Dirac propagator and pseudoscalar Yukawa interaction vertex are respectivelyThe Dirac propagator follows by multiplying the inverse kinetic matrix $/p-m$ by $/p+m$, using the Clifford algebra to obtain $(/p-m)(/p+m)=(p^2-m^2)I_4$. The vertex follows by expanding $e^{i{S}_{\mathrm{int}}}$ once: there is no factorial because it has one scalar and two distinct spinor field factors. An internal scalar has no spinor indices; an internal fermion line carries the displayed matrix with its momentum directed along fermion flow.

These rules use the coefficient exactly as printed. With the standard Hermitian ${\gamma }^5=i{\gamma }^0{\gamma }^1{\gamma }^2{\gamma }^3$, the Hermiticity of a pseudoscalar Dirac bilinear gives $(\bar{\psi }{\gamma }^5\psi {)}^{†}=-\bar{\psi }{\gamma }^5\psi$. A Hermitian physical interaction therefore requires ${\lambda }^{\ast }=-\lambda$; writing $\lambda =ig$ with real $g$ changes the vertex to $g{\gamma }^5$. This convention issue does not alter the formal rules above, and no unstated real-coupling assumption is needed.

Take incoming electron momentum $p$ and positron momentum $q$, with outgoing photon momenta $k_1,k_2$, so $p+q=k_1+k_2$. The leading electron-positron annihilation into two photons has two tree-level Feynman diagrams, differing by the order in which the two outgoing photons attach to the charged-fermion line.

The arrows indicate fermion flow. An incoming positron has the opposite flow direction to its physical momentum; the lower row uses an outgoing electron instead.

The required QED Feynman rules, with interaction $-e\bar{\psi }{\gamma }^{\mu }\psi A_{\mu }$, are a vertex $-ie{\gamma }^{\mu }$, a Dirac propagator $i(/r+m)/(r^2-m^2+i0)$, an incoming electron spinor $u(p)$, an incoming positron adjoint $\bar{v}(q)$, and an outgoing photon polarization ${ϵ}_{\mu }^{\ast }(k)$. Here Feynman slash notation means $/a={\gamma }^{\mu }{a}_{\mu }$, and the external Dirac equations are $(/p-m)u(p)=0$ and $\bar{v}(q)(/q+m)=0$.

Put $\mathcal{S}(r)=(/r+m)/(r^2-m^2+i0)$, omitting only the propagator's overall factor $i$. Multiplying the two vertices and the propagator givesThe diagrammatic expression is $i\mathcal{M}$. The relative sign is plus because the two photons are bosons; exchanging their attachments does not exchange external fermions. The amplitude is of order $e^2$.

A physical photon polarization vector obeys $k^2=0$ and $k\cdot ϵ=0$, with the gauge equivalence $ϵ\sim ϵ+\alpha k$. One can choose a transverse representative with ${ϵ}^0=0$ and normalize ${ϵ}^{\ast }\cdot ϵ=-1$. There are two physical transverse polarizations.

To prove the requested two-photon fermion Ward identity, replace ${ϵ}_1^{\ast }$ by $k_1$. Away from an internal pole, $\mathcal{S}(r)$ is the inverse of $/r-m$, with the boundary prescription understood. The first term simplifies becauseFor the second term, four-momentum conservation gives $k_1=(p-k_2)+q$, soThe two contributions are therefore $-\bar{v}/{ϵ}_2^{\ast }u$ and $+\bar{v}/{ϵ}_2^{\ast }u$, which cancel. Interchanging the photon labels proves the other contraction vanishes. By linearity, adding a multiple of either photon momentum to its polarization leaves the summed amplitude unchanged. This is gauge invariance: longitudinal gauge representatives do not describe an additional physical photon state. An individual diagram does not generally have this property.