Skip to content

Yukawa theory couples a scalar field to a fermion bilinear. It shows how spin-zero exchange generates forces and how a scalar can decay into a particle–antiparticle pair. Its simplicity makes it useful for learning amplitudes and renormalization without the gauge constraints of QED.

For one real scalar ϕ\phi and one Dirac field ψ\psi in 3+13+1 dimensions, take

L=12(∂ϕ)2−V(ϕ)+ψˉ(iγμ∂μ−m)ψ−yϕψˉψ.\begin{aligned} \mathcal L={}&\frac12(\partial\phi)^2-V(\phi) +\bar\psi(i\gamma^\mu\partial_\mu-m)\psi\\ &-y\phi\bar\psi\psi . \end{aligned}

Expand about a chosen stable vacuum with ⟨ϕ⟩=0\langle\phi\rangle=0 and V′′(0)=M2>0V''(0)=M^2>0. The coupling yy is real and dimensionless. A bounded potential and the counterterms allowed by the actual symmetries complete the model; specifying only the Yukawa vertex is insufficient at loop level. A fermion mass and a scalar Yukawa term do not generally preserve ϕ↦−ϕ\phi\mapsto-\phi by itself.

When M>2mM>2m, the scalar can decay into ψψˉ\psi\bar\psi. For a scalar at rest, the tree-level amplitude is proportional to yuˉ(p)v(p′)y\bar u(p)v(p'). Summing the two outgoing spin labels gives

∑spins∣M∣2=y2tr⁡[(p ⁣ ⁣ ⁣/+m)(p′ ⁣ ⁣ ⁣/−m)]=2y2(M2−4m2).\sum_{\mathrm{spins}}\lvert\mathcal M\rvert^2 =y^2\operatorname{tr}\bigl[(p\!\!\!/+m)(p'\!\!\!/-m)\bigr] =2y^2(M^2-4m^2).

The two-body phase-space factor supplies another square root, yielding

Γϕ→ψψˉ=y2M8π(1−4m2M2)3/2.\Gamma_{\phi\to\psi\bar\psi} =\frac{y^2M}{8\pi} \left(1-\frac{4m^2}{M^2}\right)^{3/2}.

This is for one Dirac species, no additional color multiplicity, and leading order in yy. Below threshold the on-shell two-body decay closes. The 3/23/2 power is characteristic of this scalar bilinear; a pseudoscalar coupling has a different threshold law. The spin sums and phase space follow Srednicki 2007, §§11 and 45–48.

If the scalar is heavy and external momenta are small compared with MM, its equation of motion gives ϕ≃−yψˉψ/M2\phi\simeq-y\bar\psi\psi/M^2. Substitution produces the leading local term

ΔLeff=y22M2(ψˉψ)2,\Delta\mathcal L_{\mathrm{eff}} =\frac{y^2}{2M^2}(\bar\psi\psi)^2,

with derivative corrections suppressed by further powers of momentum over MM. This is an effective interaction, not an exact replacement near the scalar pole.

Pole masses, decay rates, and scattering cross sections are the useful observables. Weak-coupling calculations require small loop corrections and controlled logarithms. A scalar expectation value shifts the fermion mass, but neither that shift nor the tree decay calculation establishes the model’s full vacuum stability.

Fermion and Yukawa Tree Amplitudes gives the detailed amplitude conventions. Compare Dirac Field for the free limit and Gross–Neveu Model for a lower-dimensional four-fermion theory with different ultraviolet behavior.

  • Srednicki, Mark. Quantum Field Theory. Cambridge University Press, 2007, §§11, 45–48, and 51–52. DOI.