The Yukawa potential

A force, in the modern telling, is not a field of influence reaching silently across empty space — it is a conversation, carried by a particle passed back and forth between the bodies that feel it. Two electrons repel because one hurls a photon at the other; two nucleons stick together because they trade pions. This picture, worked out inside time-ordered perturbation theory, replaces the classical potential with the exchange of a quantum, and it answers a question the classical potential cannot even ask: why do some forces reach across a galaxy while others die inside a nucleus? The answer is the mass of the exchanged quantum. Follow the process a+b→c+da+b\to c+d through the two orderings in which the mediator XX can be emitted and reabsorbed, sum the amplitudes as quantum mechanics demands, and the energy-violating intermediate terms recombine into a single covariant object, the propagator 1/(q2−mX2)1/(q^2-m_X^2). The exchanged XX never sits on its mass shell — momentum is conserved at each vertex, energy is not — so it is virtual, a bookkeeping quantum that exists only for the instant of the exchange. Fourier-transform that propagator back to position space and out drops the Yukawa potential V(r)∝e−mXr/rV(r)\propto e^{-m_Xr}/r: a Coulomb 1/r1/r tail, but throttled by an exponential whose scale is the Compton wavelength 1/mX1/m_X of the mediator. That one exponential is the hinge of the whole story. Set mX=0m_X=0 and it vanishes, restoring the infinite-range Coulomb law of the massless photon. Give the mediator a mass and the force is guillotined at r∼1/mXr\sim1/m_X. Yukawa read this backwards in 1935: the nuclear force reaches about a femtometre, so its carrier must weigh a couple of hundred electron masses. The pion, found twelve years later, weighed exactly that.

Key takeaways
  • A force is the exchange of a virtual quantum; summing the two time-ordered diagrams for a+b→c+da+b\to c+d gives the covariant amplitude Mfi=gagbq2−mX2\mathcal{M}_{fi}=\dfrac{g_ag_b}{q^2-m_X^2}.
  • The exchanged mediator is off mass-shell: momentum is conserved at each vertex, energy is not, so q2≠mX2q^2\neq m_X^2 — it is virtual, permitted only fleetingly.
  • Fourier-transforming the propagator 1/(q2+mX2)1/(\mathbf q^2+m_X^2) into position space yields the Yukawa potential V(r)=−g24πe−mXrrV(r)=-\dfrac{g^2}{4\pi}\dfrac{e^{-m_Xr}}{r}.
  • The mediator mass sets the range R=1/mX=ℏ/mXcR=1/m_X=\hbar/m_Xc: a massive mediator gives a short-range force, and taking mX→0m_X\to0 kills the exponential and restores the infinite-range Coulomb potential of the massless photon.
  • Yukawa inverted the logic in 1935: a ∼1.4 fm\sim1.4\ \mathrm{fm} nuclear range demands a mediator of ∼140 MeV\sim140\ \mathrm{MeV} — the pion, discovered in 1947.
  • Key idea: a force is the exchange of a virtual quantum, and the mass of that quantum fixes the range of the force, because the propagator 1/(q2−mX2)1/(q^2-m_X^2) Fourier-transforms into a screened e−mXr/re^{-m_Xr}/r potential rather than a bare 1/r1/r one.
The core insight

A force is the exchange of a virtual quantum, and the mass of that quantum fixes the range of the force: the propagator 1/(q2−mX2)1/(q^2-m_X^2) Fourier-transforms into a screened e−mXr/re^{-m_Xr}/r Yukawa potential rather than a bare 1/r1/r Coulomb one, so a massless mediator gives an infinite-range force and a massive one a force cut off at the Compton wavelength 1/mX1/m_X.

Where this lives today

The Yukawa idea is not a historical relic; it is the load-bearing template for how every force in the Standard Model is computed. In nuclear physics the long-range tail of the nucleon–nucleon interaction is still the one-pion-exchange Yukawa potential, and modern chiral effective field theory builds the full force as a controlled sum of pion and heavier-meson exchanges — one-boson-exchange potentials that fit scattering data to high precision, with QCD understood as the deeper layer beneath. In electroweak physics the massive propagators of the WW and ZZ are Yukawa potentials of ∼10−3 fm\sim10^{-3}\ \mathrm{fm} range, and their low-energy contraction to a contact interaction is Fermi theory, with the mediator mass fixing GFG_F. The very same mathematics reappears far from particle physics: screened Coulomb (Debye–Hückel) potentials in plasmas and electrolytes, and Thomas–Fermi screening of impurities in metals, are Yukawa potentials e−r/λ/re^{-r/\lambda}/r with the screening length λ\lambda playing the role of 1/mX1/m_X. And every internal line in every Feynman diagram — the photon propagator, the gluon propagator, the Higgs propagator — is the momentum-space Yukawa object 1/(q2−m2)1/(q^2-m^2). A separate modern usage borrows Yukawa's name for the Yukawa couplings between the Higgs field and fermions, which generate fermion masses after electroweak symmetry breaking; this is historically and conceptually distinct from Yukawa's meson-mediated nuclear-force mechanism.

The physical picture

The exchange picture lives in two frames at once. In the amplitude frame it is a sum over histories: the mediator XX can be launched by aa and caught by bb, or launched by bb (as its antiparticle) and caught by aa, and quantum mechanics adds these two time-ordered stories. Neither story alone is Lorentz-invariant, and in each the exchanged quantum flagrantly violates energy conservation; only their sum is covariant, and in the sum the two energy-denominators fuse into the single propagator 1/(q2−mX2)1/(q^2-m_X^2). In the potential frame the same propagator, read as scattering off a static source, is a screened Coulomb law whose screening length is the mediator's Compton wavelength. The first figure shows the two time orderings combining into one covariant propagator; the second shows how the mediator's mass sets the range of the force.

Figure
The process a + b → c + d proceeds through two time orderings — a emits the mediator X which b absorbs, or b emits anti-X which a absorbs. Momentum is conserved at every vertex but energy is not, so X is virtual (off-shell). Summing both orderings recombines them into one covariant propagator 1/(q2−mX2)1/(q^2 - m_X^2)
Summing the two time orderings gives one covariant propagator 1/(q² − m_X²) time ac bd X a emits X, b absorbs + ac bd anti-X b emits anti-X, a absorbs = ac bd q 1/(q² − m_X²)

The exchanged quantum in each ordering carries a definite three-momentum but the "wrong" energy — that mismatch Ej≠EiE_j\neq E_i is exactly what makes it virtual, and it is why no detector ever catches an exchange particle in flight. What survives into the covariant amplitude is not the two illegal intermediate states but their sum, the propagator, whose denominator q2−mX2q^2-m_X^2 is a Lorentz scalar. When that amplitude is fed back into the Born formula for scattering off a static source, its Fourier transform is a potential — and the mediator's mass, buried in the denominator, resurfaces as an exponential screening factor. The second figure is the payoff: the Yukawa potential of a real pion beside the Coulomb potential it would become if the pion were massless.

Figure
The Yukawa potential of a pion (m = 139.6 MeV) versus the Coulomb potential. The mass term cuts the force off beyond the range ℏ/m_π c = 197.3/139.6 = 1.41 fm — the reach of the nuclear force — while the massless Coulomb tail (m = 0) runs on to infinity as 1/r
The mass sets the range: massless 1/r versus massive e^(−m_X r)/r separation r (fm) |V(r)| (arb. units) 0.5 1.41 3 5 range = ℏ/m_π c = 1.41 fm (pion Compton wavelength 1/m_π) m = 0: V ∝ 1/r (Coulomb, infinite range) m_π = 139.6 MeV: V ∝ e^(−m_π r)/r

The core equations

Everything in this article can be run in either direction along one chain: field equation → potential, or amplitude → propagator → potential. We work in natural units ℏ=c=1\hbar=c=1, where mass, energy, momentum, and inverse length all share a single unit. Start with the field a static point source of a massive scalar sets up: it obeys the time-independent Klein–Gordon equation with a delta-function source of strength gg.

(∇2−mX2) ϕ(r)=−g δ3(r)(1)(\nabla^2 - m_X^2)\,\phi(\mathbf r) = -g\,\delta^3(\mathbf r) \qquad (1)

where: ϕ\phi — the static scalar field (potential per unit test charge); mXm_X — the mediator mass; gg — the coupling (source strength) at the vertex; δ3(r)\delta^3(\mathbf r) — a point source at the origin. Reading it: this is Poisson's equation ∇2ϕ=−gδ3\nabla^2\phi=-g\delta^3 with one extra term, −mX2ϕ-m_X^2\phi. That single mass term is the whole difference between an infinite-range and a short-range force. Solving equation (1) gives the Yukawa potential — the massive cousin of the Coulomb potential. Multiplying by the coupling gg of a second particle placed in the field gives the interaction energy between them.

ϕ(r)=g4πe−mXrr  ⟹  V(r)=−g24πe−mXrr(2)\phi(r)=\frac{g}{4\pi}\frac{e^{-m_Xr}}{r}\;\Longrightarrow\; V(r)=-\frac{g^2}{4\pi}\frac{e^{-m_Xr}}{r} \qquad (2)

where: r=∣r∣r=|\mathbf r| — separation; g2/4πg^2/4\pi — the dimensionless coupling strength (the scalar analogue of α\alpha); the sign is negative because scalar exchange is attractive between like sources. Reading it: the familiar 1/r1/r is multiplied by the screening factor e−mXre^{-m_Xr}, which crushes the potential to zero once rr exceeds 1/mX1/m_X. The same physics arrives from the scattering side. Following a+b→c+da+b\to c+d through both time orderings of the exchange and summing the amplitudes gives a manifestly covariant matrix element — the propagator sandwiched between two vertex couplings.

Mfi=ga gbq2−mX2,qμ=paμ−pcμ(3)\mathcal{M}_{fi}=\frac{g_a\,g_b}{q^2-m_X^2},\qquad q^\mu=p_a^\mu-p_c^\mu \qquad (3)

where: ga,gbg_a,g_b — the couplings at the two vertices; qμq^\mu — the four-momentum transferred by XX; q2=(Ea−Ec)2−(pa−pc)2q^2=(E_a-E_c)^2-(\mathbf p_a-\mathbf p_c)^2 — its invariant square. Reading it: 1/(q2−mX2)1/(q^2-m_X^2) is the propagator, the amplitude for the virtual XX to travel from one vertex to the other. In elastic scattering the transfer is spacelike, q2=−q2<0q^2=-\mathbf q^2<0, so M=−gagb/(q2+mX2)\mathcal{M}=-g_ag_b/(\mathbf q^2+m_X^2) never hits the pole. The bridge between equations (2) and (3) is the Born approximation: the static-limit amplitude is minus the Fourier transform of the potential, so inverting the transform of the propagator reproduces the Yukawa form exactly.

V(r)=− ⁣∫ ⁣d3q(2π)3 g2 eiq⋅rq2+mX2=−g24πe−mXrr(4)V(\mathbf r)=-\!\int\!\frac{d^3q}{(2\pi)^3}\,\frac{g^2\,e^{i\mathbf q\cdot\mathbf r}}{\mathbf q^2+m_X^2}=-\frac{g^2}{4\pi}\frac{e^{-m_Xr}}{r} \qquad (4)

where: the integrand is the momentum-space propagator times the coupling g2g^2; the transform is the standard three-dimensional integral ∫d3q(2π)3eiq⋅rq2+m2=e−mr4πr\int\frac{d^3q}{(2\pi)^3}\frac{e^{i\mathbf q\cdot\mathbf r}}{\mathbf q^2+m^2}=\frac{e^{-mr}}{4\pi r}. Reading it: momentum space and position space hold the same information — a pole at q2=−mX2\mathbf q^2=-m_X^2 is an exponential of scale 1/mX1/m_X in rr. Finally, the two numbers that summarize a Yukawa force: its range and its strength. The range is fixed by the mediator mass alone; the strength by the dimensionless coupling.

R=1mX=ℏmXc,αX≡g24π,V(r)→mX→0−αXr(5)R=\frac{1}{m_X}=\frac{\hbar}{m_Xc},\qquad \alpha_X\equiv\frac{g^2}{4\pi},\qquad V(r)\xrightarrow[m_X\to0]{}-\frac{\alpha_X}{r} \qquad (5)

where: RR — the range (the reduced Compton wavelength of XX); αX\alpha_X — the coupling constant (for the photon this is the fine-structure constant α≈1/137\alpha\approx1/137). Reading it: mass and range are reciprocals, and the massless limit smoothly restores the pure Coulomb law — the Yukawa potential contains Coulomb as its mX→0m_X\to0 special case.

Derivation ASumming the two time orderings into a covariant propagator

The goal is to show, without hand-waving, how the propagator 1/(q2−mX2)1/(q^2-m_X^2) emerges from ordinary second-order perturbation theory applied to particle exchange. In time-ordered perturbation theory each diagram contains an off-shell intermediate state whose energy does not satisfy the on-shell relation, and the individual terms are frame-dependent; their sum, however, recombines into the Lorentz-covariant propagator. We track the two energy denominators through their recombination and watch a factor of 2EX2E_X cancel to leave the invariant result. This is the calculation that first makes "a force is an exchanged virtual particle" a quantitative statement rather than a slogan.

Step 1. Write the second-order amplitude Tfi=∑jVfjVji/(Ei−Ej)T_{fi}=\sum_j V_{fj}V_{ji}/(E_i-E_j), where VV are the (non-invariant) transition matrix elements and jj runs over intermediate states. Each vertex contributes a coupling gg and, in the standard relativistic normalization, a factor (2Ek)−1/2(2E_k)^{-1/2} for every particle kk attached to it. Why: this is textbook Rayleigh–Schrödinger perturbation theory; the only new ingredient is the relativistic normalization that will later strip off cleanly into the invariant amplitude.

Step 2. Enumerate the two orderings. In ordering (a), aa emits XX (intermediate state c+X+bc+X+b) which bb later absorbs; energy conservation is broken by Ei−Ej=Ea−Ec−EXE_i-E_j=E_a-E_c-E_X. In ordering (b), bb emits anti-XX (intermediate state a+d+Xa+d+X) which aa absorbs; here Ei−Ej=Eb−Ed−EXE_i-E_j=E_b-E_d-E_X. In both, three-momentum is conserved at every vertex, so XX carries q=pa−pc\mathbf q=\mathbf p_a-\mathbf p_c and EX=+q2+mX2E_X=+\sqrt{\mathbf q^2+m_X^2}. Why: energy non-conservation at intermediate stages is the precise, technical meaning of "virtual" — the mediator has a real momentum but the wrong energy for its mass.

Step 3. Add the two denominators. Overall energy conservation Ea+Eb=Ec+EdE_a+E_b=E_c+E_d gives Eb−Ed=Ec−EaE_b-E_d=E_c-E_a, so writing Δ≡Ea−Ec\Delta\equiv E_a-E_c the two orderings contribute 1Δ−EX\frac{1}{\Delta-E_X} and 1−Δ−EX\frac{1}{-\Delta-E_X}. Their sum is 1Δ−EX−1Δ+EX=2EXΔ2−EX2\frac{1}{\Delta-E_X}-\frac{1}{\Delta+E_X}=\frac{2E_X}{\Delta^2-E_X^2}. Why: the two frame-dependent, energy-violating pieces are individually ugly, but their sum has a clean denominator Δ2−EX2\Delta^2-E_X^2 — the first sign that a covariant quantity is hiding inside.

Step 4. Cancel the 2EX2E_X. Each vertex carried a (2EX)−1/2(2E_X)^{-1/2}, so the exchanged particle contributed an overall 1/2EX1/2E_X; multiplying by the 2EX2E_X from Step 3 leaves no factor of EXE_X at all. What remains is Tfi∝1/(Δ2−EX2)T_{fi}\propto 1/(\Delta^2-E_X^2). Why: the cancellation is the technical miracle of the derivation — the arbitrary intermediate energy EXE_X disappears, and with it every trace of the frame in which we ordered the diagrams.

Step 5. Recognize the invariant. Since EX2=q2+mX2E_X^2=\mathbf q^2+m_X^2 and Δ=Ea−Ec\Delta=E_a-E_c, the denominator is Δ2−EX2=(Ea−Ec)2−q2−mX2=q2−mX2\Delta^2-E_X^2=(E_a-E_c)^2-\mathbf q^2-m_X^2=q^2-m_X^2 with qμ=paμ−pcμq^\mu=p_a^\mu-p_c^\mu. Stripping the external normalization factors defines the Lorentz-invariant amplitude Mfi=gagb/(q2−mX2)\mathcal{M}_{fi}=g_ag_b/(q^2-m_X^2). Why: q2q^2 and mX2m_X^2 are both Lorentz scalars, so the summed amplitude is frame-independent — the propagator is the covariant residue of two non-covariant histories.

Go deeper · Derivation B · 2 minFourier-transforming the propagator into the Yukawa potential
Derivation BFourier-transforming the propagator into the Yukawa potential

The aim is to turn the momentum-space propagator of equation (3) into the position-space potential of equation (2), making explicit why a mediator mass in the denominator becomes an exponential screening length in space. The vehicle is the Born approximation, which equates a scattering amplitude with the Fourier transform of the potential doing the scattering; running it backwards recovers the potential from the amplitude, and the entire calculation reduces to one standard integral whose exponential answer is the origin of "short-range force". This is the step that converts a Feynman diagram into a force law you could put into a Schrödinger equation.

Step 1. In the non-relativistic static limit, set Ea=EcE_a=E_c so q2=−q2q^2=-\mathbf q^2 and the amplitude (3) becomes M(q)=−g2/(q2+mX2)\mathcal{M}(\mathbf q)=-g^2/(\mathbf q^2+m_X^2) (taking ga=gb=gg_a=g_b=g). The Born approximation states that scattering off a static potential V(r)V(\mathbf r) has amplitude M(q)=−V~(q)\mathcal{M}(\mathbf q)=-\tilde V(\mathbf q), the Fourier transform of VV. Why: this is the dictionary between the two languages — a Feynman amplitude in momentum space and a potential in position space are Fourier partners, so identifying them lets us read one off the other.

Step 2. Identify V~(q)=g2/(q2+mX2)\tilde V(\mathbf q)=g^2/(\mathbf q^2+m_X^2) and invert the transform: V(r)=∫d3q(2π)3 V~(q) eiq⋅r=∫d3q(2π)3g2 eiq⋅rq2+mX2V(\mathbf r)=\int\frac{d^3q}{(2\pi)^3}\,\tilde V(\mathbf q)\,e^{i\mathbf q\cdot\mathbf r}=\int\frac{d^3q}{(2\pi)^3}\frac{g^2\,e^{i\mathbf q\cdot\mathbf r}}{\mathbf q^2+m_X^2}. Why: the potential is nothing but the propagator viewed in position space; every property of the force is already encoded in the analytic structure of 1/(q2+mX2)1/(\mathbf q^2+m_X^2).

Step 3. Do the angular integral. Aligning r\mathbf r with the polar axis, ∫dΩ eiqrcos⁡θ=4π sin⁡(qr)/(qr)\int d\Omega\,e^{iqr\cos\theta}=4\pi\,\sin(qr)/(qr), so V(r)=g22π2r∫0∞q sin⁡(qr)q2+mX2 dqV(r)=\frac{g^2}{2\pi^2 r}\int_0^\infty\frac{q\,\sin(qr)}{q^2+m_X^2}\,dq. Why: the potential is spherically symmetric because the propagator depends only on ∣q∣|\mathbf q|, and the angular integral collapses the vector transform to a single radial one.

Step 4. Evaluate the radial integral, a standard contour result: ∫0∞q sin⁡(qr)q2+mX2 dq=π2e−mXr\int_0^\infty\frac{q\,\sin(qr)}{q^2+m_X^2}\,dq=\frac{\pi}{2}e^{-m_Xr}. Hence V(r)=g22π2r⋅π2e−mXr=g24πe−mXrrV(r)=\frac{g^2}{2\pi^2r}\cdot\frac{\pi}{2}e^{-m_Xr}=\frac{g^2}{4\pi}\frac{e^{-m_Xr}}{r} (attractive: −g2/4π e−mXr/r-g^2/4\pi\,e^{-m_Xr}/r). Why: the pole of the propagator at q=±imXq=\pm im_X sits off the real axis by exactly mXm_X, and closing the contour turns that imaginary pole into the real exponential e−mXre^{-m_Xr} — the mass literally is the inverse screening length.

Step 5. Check the massless limit. As mX→0m_X\to0 the exponent vanishes and V(r)→g2/4πrV(r)\to g^2/4\pi r, the Coulomb potential; a massless mediator has an infinite screening length and hence an infinite range. Why: this both validates the calculation against the known photon case and shows the deep unity — Coulomb is not a separate law but the mX=0m_X=0 boundary of the Yukawa family.

Go deeper · Derivation C · 1 minSolving the static field equation directly and reading off the range
Derivation CSolving the static field equation directly and reading off the range

The purpose of this shorter derivation is to confirm the Yukawa potential straight from the field equation (1), independent of any scattering argument, and to extract the physical range as the distance over which the field falls by a factor ee. Verifying the same e−mXr/re^{-m_Xr}/r from a completely different starting point — a classical field sourced by a point charge — is exactly the kind of cross-check that turns a formula into a law, and it isolates the range as a clean, measurable number.

Step 1. Away from the origin (r≠0r\neq0) the source term drops and equation (1) reads (∇2−mX2)ϕ=0(\nabla^2-m_X^2)\phi=0. Using ∇2f(r)=1rd2dr2(rf)\nabla^2 f(r)=\frac1r\frac{d^2}{dr^2}(rf), substitute the trial form ϕ=Ae−mXr/r\phi=Ae^{-m_Xr}/r: then rϕ=Ae−mXrr\phi=Ae^{-m_Xr} and d2dr2(rϕ)=mX2Ae−mXr\frac{d^2}{dr^2}(r\phi)=m_X^2Ae^{-m_Xr}, so ∇2ϕ=mX2ϕ\nabla^2\phi=m_X^2\phi and the equation is satisfied identically. Why: the trial function is engineered so that the Laplacian returns mX2m_X^2 times itself, exactly matching the mass term — the exponential is the unique spherically symmetric decaying solution.

Step 2. Fix the constant AA from the source. Near r=0r=0 the exponential →1\to1 and ϕ→A/r\phi\to A/r, whose Laplacian is the Coulomb singularity ∇2(A/r)=−4πA δ3(r)\nabla^2(A/r)=-4\pi A\,\delta^3(\mathbf r). Matching to the source −gδ3(r)-g\delta^3(\mathbf r) gives A=g/4πA=g/4\pi, so ϕ(r)=g4πe−mXrr\phi(r)=\frac{g}{4\pi}\frac{e^{-m_Xr}}{r}. Why: the strength of a 1/r1/r field at the origin is set entirely by the enclosed source, so the normalization is inherited unchanged from the Coulomb problem — only the tail is modified.

Step 3. Read off the range. The field falls to 1/e1/e of its (screening-stripped) value over the distance where mXr=1m_Xr=1, i.e. R=1/mX=ℏ/mXcR=1/m_X=\hbar/m_Xc. Restoring units, R=ℏc/(mXc2)R=\hbar c/(m_Xc^2), so a mediator of mass mXc2=200 MeVm_Xc^2=200\ \mathrm{MeV} gives R=197.3/200≈1 fmR=197.3/200\approx1\ \mathrm{fm}. Why: this reciprocal relation between mass and range is the entire predictive content of the Yukawa idea — measure a force's reach and you have weighed its carrier.

Limiting cases

CaseResultMeaning
Massless mediator, mX→0m_X\to0V(r)→−g24πrV(r)\to-\dfrac{g^2}{4\pi r}pure Coulomb 1/r1/r; infinite range — the photon, the graviton
Long range, r≫1/mXr\gg1/m_XV∼e−mXr→0V\sim e^{-m_Xr}\to0force cut off exponentially beyond the Compton wavelength
Low momentum transfer, q2≪mX2\mathbf q^2\ll m_X^2M→−g2/mX2\mathcal{M}\to-g^2/m_X^2 (constant)point/contact interaction — Fermi's four-fermion theory
High momentum transfer, q2≫mX2\mathbf q^2\gg m_X^2M→−g2/q2\mathcal{M}\to-g^2/\mathbf q^2mass irrelevant; Coulomb/Rutherford 1/q21/\mathbf q^2 scattering
On mass-shell, q2→mX2q^2\to m_X^2propagator pole 1/(q2−mX2)→∞1/(q^2-m_X^2)\to\inftymediator becomes real and can be produced (a resonance)

Worked examples

1 · Weighing the pion from the range of the nuclear force (fully worked)

The long-range part of the nuclear force between nucleons extends over roughly 1.4 fm and is dominated by one-pion exchange. Modelling this long-range interaction with a Yukawa potential V∝e−mπr/rV\propto e^{-m_\pi r}/r, (a) what mediator mass does that range imply? (b) Invert with the measured pion mass mπ=139.6 MeV/c2m_\pi=139.6\ \mathrm{MeV}/c^2 to predict the range. Use ℏc=197.3 MeV⋅fm\hbar c=197.3\,{\rm MeV}\cdot{\rm fm}.

Step 1. From equation (5) the range and mediator mass are reciprocals, R=ℏc/(mXc2)R=\hbar c/(m_Xc^2), so a mass is a range read backwards: mXc2=ℏc/Rm_Xc^2=\hbar c/R. The mass term in the field equation is nothing but the inverse range of the force it carries.

Step 2. Insert the observed range R≈1.4 fmR\approx1.4\ \mathrm{fm}: mXc2=197.3 MeV⋅fm1.4 fm=141 MeVm_Xc^2=\dfrac{197.3\ \mathrm{MeV\cdot fm}}{1.4\ \mathrm{fm}}=\mathbf{141\ MeV}. This is essentially Yukawa's 1935 estimate — a particle a few hundred times heavier than the electron, intermediate in mass between electron and proton, which is why he called it a meson ("middle-weight").

Step 3. Invert with the modern pion mass to check consistency: R=ℏcmπc2=197.3 MeV⋅fm139.6 MeV=1.41 fmR=\dfrac{\hbar c}{m_\pi c^2}=\dfrac{197.3\ \mathrm{MeV\cdot fm}}{139.6\ \mathrm{MeV}}=\mathbf{1.41\ fm}. The measured pion reproduces the nuclear range on the nose.

Step 4. Interpret. Yukawa read a particle's mass off a force's reach a full twelve years before that particle was seen (the pion was found in cosmic rays in 1947). A Yukawa interaction mediated by a particle of mass mπm_\pi has characteristic range ℏ/(mπc)\hbar/(m_\pi c); a 1.4 fm range therefore implies mπ≈140 MeVm_\pi\approx140\,{\rm MeV}, remarkably close to the measured pion mass. The Yukawa potential turned a qualitative force into a quantitative prediction of a new particle.

2 · From the W-boson propagator to Fermi's constant

The charged-current weak interaction is mediated by the WW boson, with mass MW=80.4 GeVM_W=80.4\,{\rm GeV} and coupling gW≈0.65g_W\approx0.65. (a) Estimate its range. (b) Show that at low momentum transfer the massive propagator reduces to a contact interaction, and recover Fermi's constant GFG_F.

Step 1. Range from equation (5): R=ℏcMWc2=197.3 MeV⋅fm80 400 MeV=2.45×10−3 fm=2.5×10−18 mR=\dfrac{\hbar c}{M_Wc^2}=\dfrac{197.3\ \mathrm{MeV\cdot fm}}{80\,400\ \mathrm{MeV}}=2.45\times10^{-3}\ \mathrm{fm}=\mathbf{2.5\times10^{-18}\ m}. The weak force reaches roughly a thousandth of a proton radius — effectively a point contact, which is why β\beta decay looked for decades like a zero-range interaction.

Step 2. Take the low-qq limit of the propagator. For everyday weak processes q2≪MW2\mathbf q^2\ll M_W^2, so from the limiting-cases table M=gW2q2+MW2→gW2MW2\mathcal{M}=\dfrac{g_W^2}{\mathbf q^2+M_W^2}\to\dfrac{g_W^2}{M_W^2} — a constant, independent of q\mathbf q. A momentum-independent amplitude is the Fourier transform of a delta function: the exchange has shrunk to a contact interaction. The heavy mediator hides the propagator's structure and mimics Fermi's original four-fermion vertex.

Step 3. Match to Fermi theory, whose contact coupling is GF/2=gW2/(8MW2)G_F/\sqrt2=g_W^2/(8M_W^2). Evaluate: gW28MW2=(0.65)28 (80.4 GeV)2=0.42258×6464 GeV2=8.2×10−6 GeV−2\dfrac{g_W^2}{8M_W^2}=\dfrac{(0.65)^2}{8\,(80.4\ \mathrm{GeV})^2}=\dfrac{0.4225}{8\times6464\ \mathrm{GeV^2}}=8.2\times10^{-6}\ \mathrm{GeV^{-2}}, so GF=2×8.2×10−6=1.16×10−5 GeV−2G_F=\sqrt2\times8.2\times10^{-6}=\mathbf{1.16\times10^{-5}\ GeV^{-2}}.

Step 4. Compare with experiment: the measured Fermi constant is GF=1.166×10−5 GeV−2G_F=1.166\times10^{-5}\ \mathrm{GeV^{-2}}. The Yukawa propagator of an 80 GeV80\ \mathrm{GeV} mediator reproduces Fermi's constant to three significant figures — the weak force is not intrinsically feeble, it is merely short-ranged, its apparent weakness at low energy being the 1/MW21/M_W^2 suppression of a very heavy propagator.

Physics of the portrait

The Yukawa potential is where a single equation predicted a single particle — a rare, clean instance of theory running ahead of experiment, and of a young physicist reading a new particle's mass straight off the reach of a force.

Origins · 1 minHideki Yukawa and how this was discovered
HY
Hideki Yukawa
1907–1981 · meson theory and the nuclear force

In 1935 Hideki Yukawa, a twenty-seven-year-old at Osaka University, faced a puzzle that had defeated everyone: what holds the nucleus together against the electrostatic repulsion of its protons, and why does that force vanish just outside the nucleus? Electromagnetism, mediated by the massless photon, reaches to infinity — but the nuclear force plainly does not. Yukawa's leap was to ask what kind of mediator would give a short-range force, and to answer with the field equation: a mediator with mass. Working from the relativistic wave equation for a massive field, he found the potential e−mr/re^{-m r}/r and saw that its range is the mediator's Compton wavelength ℏ/mc\hbar/mc. Turning it around, the observed nuclear range of about a femtometre demanded a quantum some two hundred times the electron mass — heavier than an electron, lighter than a proton, a genuinely new kind of particle he named the meson. The prediction was met first with skepticism and then with a false alarm: the muon, discovered in cosmic rays in 1936 with roughly the right mass, was mistaken for Yukawa's particle until it proved to ignore the nuclear force entirely. The true carrier, the pion, was finally found in 1947 by Powell's group in photographic emulsions exposed on mountaintops, at 140 MeV — precisely Yukawa's estimate. He received the 1949 Nobel Prize in Physics, Japan's first. His deeper legacy outlived the specific model: the propagator 1/(q2−m2)1/(q^2-m^2) he introduced is the internal line of every Feynman diagram, and "the mass of the mediator sets the range of the force" is now one of the load-bearing sentences of the Standard Model.

History
1935
— Yukawa publishes the meson theory of nuclear forces, predicting a mediator of about 200 electron masses
1936
— Anderson and Neddermeyer find a penetrating cosmic-ray particle, later identified as the muon
1947
— Powell, Lattes, and Occhialini discover the pion at 140 MeV, the genuine Yukawa meson
1948
— Feynman's spacetime diagrams recast exchange forces as covariant propagators
1949
— Yukawa receives the Nobel Prize in Physics, the first Japanese laureate

Going further

Texts. Thomson, Modern Particle Physics, §5.1–5.1.1 (time-ordered perturbation theory and the propagator) and §5.2 · Griffiths, Introduction to Elementary Particles, 2nd ed., Ch. 6 (the Feynman calculus) and §7.1 · Halzen and Martin, Quarks and Leptons, §4.5–4.6 (the propagator and the Yukawa potential) · Peskin and Schroeder, An Introduction to Quantum Field Theory, §4.7 (the Yukawa potential from the Born approximation) · Schwartz, Quantum Field Theory and the Standard Model, §3.4 · Aitchison and Hey, Gauge Theories in Particle Physics, Vol. 1, Ch. 5. The following key results tie the article together. The two time-ordered denominators recombine as 1/(Δ−EX)−1/(Δ+EX)=2EX/(Δ2−EX2)1/(\Delta-E_X)-1/(\Delta+E_X)=2E_X/(\Delta^2-E_X^2), with the 2EX2E_X cancelling against the vertex normalisation to give M=gagb/(q2−mX2)\mathcal{M}=g_ag_b/(q^2-m_X^2). The Fourier integral used throughout is ∫d3q(2π)3eiq⋅rq2+m2=e−mr4πr\int\frac{d^3q}{(2\pi)^3}\frac{e^{i\mathbf q\cdot\mathbf r}}{\mathbf q^2+m^2}=\frac{e^{-mr}}{4\pi r}, and the field-equation identity (∇2−mX2)(e−mXr/r)=−4πδ3(r)(\nabla^2-m_X^2)(e^{-m_Xr}/r)=-4\pi\delta^3(\mathbf r) confirms the static solution directly. Numerically: the pion range is ℏc/mπc2=197.3/139.6=1.41 fm\hbar c/m_\pi c^2=197.3/139.6=1.41\ \mathrm{fm}, and the weak-force contraction GF/2=gW2/8MW2G_F/\sqrt{2}=g_W^2/8M_W^2 yields GF=1.166×10−5 GeV−2G_F=1.166\times10^{-5}\ \mathrm{GeV^{-2}} for MW=80.4 GeVM_W=80.4\ \mathrm{GeV} and gW=0.65g_W=0.65. The values mπ=139.6 MeVm_\pi=139.6\ \mathrm{MeV}, MW=80.4 GeVM_W=80.4\ \mathrm{GeV}, and ℏc=197.3 MeV⋅fm\hbar c=197.3\ \mathrm{MeV\cdot fm} are standard tabulated quantities.

Prerequisites. Fermi's Golden Rule · Perturbation Theory and Particle Exchange · Natural Units · Relativistic Kinematics. Leads to. Feynman Diagrams and Virtual Particles · Introduction to QED · The Photon Propagator · The W-Boson Propagator · Rutherford and Mott Scattering.

Ask the tutor: Ask me to walk through how the same 1/(q2−m2)1/(q^2-m^2) propagator gives Rutherford scattering for a massless photon and a contact four-fermion interaction for a heavy W — one formula, two famous force laws.

Why does giving the mediator a mass make the force short-ranged

Because the mediator mass is the inverse of the force's range. In momentum space the exchange amplitude is 1/(q2+mX2)1/(\mathbf q^2+m_X^2), whose pole sits a distance mXm_X off the real axis; Fourier-transforming to position space turns that displaced pole into the factor e−mXre^{-m_Xr}, an exponential that falls by 1/e1/e over the distance R=1/mX=ℏ/mXcR=1/m_X=\hbar/m_Xc, the mediator's Compton wavelength. A massless mediator has its pole on the axis (mX=0m_X=0), the exponential becomes 1, and the force reaches to infinity as 1/r1/r — Coulomb. So "massive" and "short-range" are two descriptions of the same pole location, and the heavier the carrier, the shorter its reach.

Is the exchanged virtual particle real, and does it violate energy conservation

It is not a real particle, and no net energy is violated. In each time-ordered diagram the exchanged quantum has a definite momentum but the "wrong" energy for its mass, so intermediate energy is unbalanced (Ej≠EiE_j\neq E_i) — that imbalance is precisely what "virtual" means. But the initial and final states have equal total energy, so overall energy is conserved; the mismatch lives only in the transient internal state. In the covariant sum the internal line is off-shell (q2≠mX2q^2\neq m_X^2) and is a mathematical propagator, an integration variable, not something that could be caught in a detector between the vertices. The uncertainty-principle "borrowing" picture is a mnemonic that gets the range right for the wrong reason.

Why is scalar (Yukawa) exchange attractive while photon exchange repels like charges

The sign of the force is set by the spin of the mediator. Exchange of a spin-0 scalar produces a universally attractive potential between like sources — every pair of nucleons attracts, which is what binds a nucleus. Exchange of a spin-1 vector (the photon) carries an extra relative minus sign from the vector coupling, so like charges repel and opposites attract, as in electrostatics. Exchange of a spin-2 tensor (the graviton) is attractive again, and universally so, which is why gravity only pulls. The Yukawa envelope e−mr/re^{-m r}/r is common to all three; only the mediator's spin decides the sign in front of it.

How is the Yukawa potential related to the Klein-Gordon equation

They are the static solution and the field equation. The Klein–Gordon equation (□+m2)ϕ=0(\Box+m^2)\phi=0 is the full relativistic wave equation of a free massive scalar field; adding a static point source turns it into (∇2−m2)ϕ=−gδ3(r)(\nabla^2-m^2)\phi=-g\delta^3(\mathbf r), and that equation's solution is the Yukawa potential ϕ=(g/4π)e−mr/r\phi=(g/4\pi)e^{-mr}/r. So the Yukawa potential is what a massive Klein–Gordon field looks like in the neighbourhood of a stationary charge, exactly as the Coulomb potential is what a massless field (obeying Laplace's equation) looks like around a point charge. The mass term in the wave equation becomes the exponential in the potential.

Why was Yukawa's meson confused with the muon

Because the muon, discovered in cosmic rays in 1936, happened to have about the right mass — roughly 100 MeV, in Yukawa's predicted band — so it was natural to identify it as the nuclear-force carrier. The trouble emerged over the next decade: the muon passed through matter far too easily, barely interacting with nuclei, whereas Yukawa's meson had to couple strongly to nucleons to bind them. The 1947 discovery of the pion — which does interact strongly and decays into the muon — resolved it. The lesson: matching a mass is not enough; the mediator must also have the right coupling. The muon turned out to be a heavy electron, unrelated to the nuclear force.

If the strong force is Yukawa exchange, why do we need QCD

Because one-pion exchange is only the long-range tail of a deeper theory. The Yukawa potential is an excellent effective description of the nucleon–nucleon force at separations beyond about a femtometre, where a single pion is exchanged; at shorter distances heavier mesons and multi-pion exchanges pile in, and inside a nucleon the relevant degrees of freedom are quarks and gluons governed by quantum chromodynamics. QCD is the fundamental theory; the Yukawa potential is what it reduces to at low energy and long range, in the same way that a molecular van der Waals force is a residue of the underlying electromagnetism. Chiral effective field theory makes this precise, deriving the Yukawa tail systematically from QCD's symmetries.

How does the coupling g²/4π compare to the fine-structure constant, and why is the strong force strong

The dimensionless coupling of a Yukawa force is αX=g2/4π\alpha_X=g^2/4\pi, the direct analogue of the electromagnetic α=e2/4π≈1/137≈0.0073\alpha=e^2/4\pi\approx1/137\approx0.0073. For the pion–nucleon interaction the measured value is g2/4π≈14g^2/4\pi\approx14 — nearly two thousand times larger than α\alpha. That is the quantitative meaning of "strong": at each vertex the amplitude is proportional to gg, so a coupling of order 10 rather than 0.1 makes strong-interaction rates overwhelm electromagnetic ones, and it also means perturbation theory (an expansion in g2/4πg^2/4\pi) converges poorly for the strong force — one of the reasons QCD required non-perturbative tools like the lattice. The range and the strength are independent knobs: the pion sets the range through its mass, the coupling sets the strength through g2/4πg^2/4\pi.

Before this

Where this leads