Electron-positron annihilation into a muon pair Created 2026-10-05 Updated 2026-10-06
At tree level, the Feynman diagram has an electron-positron current and a muon-antimuon current joined by one virtual photon. In the massless limit, the spin average and final-state spin sum give , and the integrated relativistic scattering cross-section is .
When the full Lorentz-invariant phase-space measure treats identical final particles as labeled, each physical final configuration is counted times. The unlabelled relativistic scattering cross-section therefore includes . For two identical scalars, this divides the full-angle result by two; equivalently, integrate over a region containing only one representative of each exchanged pair.
With relativistically normalized external states, is the reduced S-matrix element after removing the four-momentum conservation distribution. A relativistic scattering cross-section combines with the Lorentz-invariant phase-space measure and the invariant flux factor. Identical-leg factorials in Feynman vertices and the identical final-state symmetry factor in the event integral have different roles.
In the static charge density model, define the normalized Fourier transform
Here is the spatial momentum transfer divided by the reduced Planck constant: has units of inverse length. Physical momentum transfer is . We keep and the speed of light explicit; compute , expressing the numerator in and in , to obtain in inverse femtometres. The spatial is distinct from the invariant four-momentum transfer used in relativistic scattering cross-sections.
For spherical symmetry, take the polar axis along and use spherical coordinates. With the cosine of the polar angle,
Every radial shell contributes an amplitude. At all shells add with the same sign; at larger their contributions can cancel. This cancellation is the origin of a diffraction minimum in the nuclear charge form factor, rather than an absence of electric charge in the nucleus.
The static Fourier transform interpretation is appropriate for the nuclear model and kinematic regime specified here. A fully relativistic composite system requires frame and current conventions; its electric form factor is not universally the transform of a frame-independent three-dimensional charge density.
Put and neglect electron and final-quark masses compared with . First fix a final color charge; its multiplicity will be restored for the inclusive rate. The photon-exchange scattering amplitude is, up to an overall phase,
Averaging over the two spins of each incident particle and summing final spins gives
By the gamma matrix trace identities, the first trace is , and likewise for the second. Contracting the two tensors cancels the terms proportional to , giving
In the centre-of-mass frame all four energies are . The first pair of scalar products equals each, and the second pair equals each. Hence .
For clarity, the printed flux hint has a factor-of-four inconsistency: when , the prefactor multiplying the spin-averaged amplitude is , not . Alternatively, is correct if . Use this massless two-body scattering flux normalization, keeping the spin average separate.
To obtain the two-body Lorentz-invariant phase space directly, integrate the momentum-conserving spatial delta function, leaving . The remaining radial integral is , so . The relativistic scattering cross-section is consequently
Since , the integrated electron-positron annihilation into a quark pair result, for this fixed color, is
The first two printed formulas are thus color-stripped. Distinct final colors are orthogonal states and must be summed, producing . At leading order each accessible quark flavor contributes independently. Away from thresholds and narrow resonances, quark-hadron duality relates the inclusive parton rate to the observed inclusive electron-positron annihilation into hadrons, giving
This formula concerns the photon channel, assumes is large compared with the relevant quark masses and the strong-interaction scale, and does not include Z boson exchange.
Beyond leading order, an exactly exclusive pair of colored quarks is not an observable asymptotic state because of confinement and hadronization. Moreover its virtual corrections have infrared divergences; real soft or collinear gluon emission cannot be excluded if a finite rate is wanted. The inclusive sum includes both the corrected two-parton channel and channels such as , with their degenerate soft and collinear limits. Real and virtual infrared divergences then cancel for an infrared-safe observable. A specified jet observable can also be meaningful if it defines an infrared-safe observable; a bare two-quark final state is not.
A convenient inclusive calculation uses the optical theorem. For the electromagnetic hadronic current , define
The hadronic electromagnetic-current spectral density is obtained from its discontinuity. In this convention the hadronic R ratio is , where at this order. Cutting the forward amplitude sums all allowed hadronic intermediate states, so perturbative Quantum chromodynamics applied to this inclusive quantity implements the required sum of real and virtual contributions.
Finally, write , , , and . The given truncated correction is
Differentiate with respect to , using the stated one-loop beta function and . The derivative of cancels the explicit logarithmic derivative at order . Keeping the residual derivative of the truncated expression gives
Thus the principle of minimal sensitivity, applied to this expression with nonzero coupling and , selects the one-loop stationary scale for the hadronic annihilation correction
This is a prescription based on the residual order- scale dependence of an order- approximation. At strictly order , the derivative already vanishes for every renormalization scale; the prescription does not assert exact scale dependence of the observable. It is useful only while the selected running coupling remains perturbative. If or , the stated derivative imposes no unique scale.
Use the Feynman vertex for the scalar Yukawa coupling, the internal scalar propagator , incoming factors and outgoing Dirac adjoints . Define as the coefficient in the connected S-matrix element after removing its overall momentum-conserving Dirac delta function. A direct pairing contributes
The crossed pairing interchanges the identical outgoing fermions and therefore has a relative fermionic sign of minus. With and , the result is
It changes sign under exchange of the two outgoing momenta and spin labels, as required for identical fermions. No fermion spin sum or identical-particle phase-space factor belongs in this amplitude; those enter a relativistic scattering cross-section calculation.
Work in natural units , at temperatures far below the -boson mass but of order or above the electron mass. The Fermi interaction then supplies a local four-fermion amplitude. The Fermi constant has dimension energy. In scattering of a thermal lepton from a nonrelativistic nucleon, the nucleon current matrix element has size , while the lepton current has size . Thus . The two-body relativistic scattering cross-section contains the kinematic factor ; when the final and initial lepton energies are comparable, their momentum ratio is of order one. Consequently . This estimates temperatures above the electron mass and mass-splitting thresholds; those scales only affect the dimensionless rate function near freeze-out.
The light leptons in the thermal bath have number density and speed of order unity. Their conversion rate per nucleon is therefore
For example, and change the neutron-to-proton ratio. The bath density here is the lepton density, not the much smaller baryon density. Dimensionless weak couplings, phase-space integrals and the several reaction channels set a numerical coefficient. Electron masses and the neutron-proton splitting also become important for a precision rate near cosmological weak freeze-out.
For an unpolarized initial fermion, average over its two spin states and sum over the unobserved final spin:
The initial scalar has only one spin state. One can evaluate this sum directly from normalized Dirac spinors, or use consistent fermion spin sums to express it as a trace. It is the squared sum of both tree scattering amplitudes, not the sum of their separate squares.
The relativistic scattering cross-section is obtained by integrating the Lorentz-invariant phase-space measure and dividing by the invariant flux factor:
There is no identical-final-particle factor, because the outgoing scalar and fermion are distinct. All energies are positive, with and .
Equivalently, in the centre-of-momentum frame set . Integrating the energy delta function in the relativistic two-body phase space gives
For this elastic process , so integrate over the full solid angle. This supplies the requested prescription without evaluating the angular integral.
There are exactly two tree Feynman diagrams: the incident scalar can be absorbed before the final scalar is emitted, or after it. The internal fermion four-momenta are respectively and . A scalar exchange diagram would require an absent scalar self-interaction, so there is no additional tree channel.
Define the S-matrix convention . The two tree scattering amplitudes are
Using the Dirac propagator from the preceding Feynman rules,
The two terms add with the same relative sign: neither diagram exchanges identical external fermions. Their interference must be retained when squaring the complete scattering amplitude. In the phase convention above, , which gives an equivalent simplified numerator. An overall phase depends on the S-matrix convention and does not change a relativistic scattering cross-section.
Keep the mostly-plus Minkowski metric and the Fourier transform . Write , with . To fix the otherwise unspecified phase of and the Dirac adjoint, take
where are ordinary mostly-minus gamma matrices. Then , , and the fermionic time-derivative term is . In these conventions the interaction with real is Hermitian: in mostly-minus notation it is . If one instead calls the square-one chirality matrix, its coefficient must be to represent the same interaction. These phase choices leave physical relativistic scattering cross-sections unchanged.
Expanding yields the following Feynman rules with relativistically normalized external states:
Strip the overall four-momentum conservation delta function when defining the scattering amplitude. There are no further bare interaction vertices, no gauge fixing and no Faddeev-Popov ghost fields in this theory. Renormalized higher-order calculations add the required counterterms; these are additional to the rules of the displayed classical Lagrangian density.
Substitute the preceding spin average into the supplied differential scattering cross-section. Since ,
The angular integral is . Consequently the relativistic scattering cross-section for the massless electron-positron annihilation into a muon pair is
The outgoing particles are distinct, so there is no identical-particle factor. Comparing with the requested monomial form gives
The integrated answer depends only on the incoming invariant ; the angular variables have been integrated out.
In four spacetime dimensions with , a relativistic scattering cross-section has mass dimension , since it is an area. The action is dimensionless, so the Lagrangian density has dimension four. The scalar kinetic term then gives , and the cubic term gives .
Each term of the four-point tree amplitude in phi cubed theory has dimension . Hence . In the final integral,
while the numerical factors are dimensionless. Therefore
Including the identical-particle factor changes only a dimensionless constant and preserves this check.
Continuum amplitudes. The leading quantum electrodynamics diagram is electron-positron annihilation into a virtual photon, followed by conversion to a muon pair. Replacing the outgoing muons with a quark-antiquark pair gives the quark-level hadronic diagram. The first two panels below show these processes.
Figure 1.
Original Feynman diagrams for electromagnetic pair production, vector-quarkonium annihilation, and the open-bottom decay topology
. Details.
The Feynman rules give, up to signs and factors of ,
The two vertices supply , and the photon propagator supplies . Each external current scales as , so the complete amplitude is of order , and its squared magnitude is of order . The two-body relativistic scattering cross-section supplies the remaining from flux and phase space. Integrating the angular distribution yields . Keeping only the squared propagator while ignoring the external currents would incorrectly predict .
A quark of charge contributes times the massless muon rate. With open and , the hadronic R ratio is
This is the requested massless leading-QED continuum approximation. Strong corrections, physical hadron thresholds, resonance contributions and Z exchange are separate effects; the pure photon formula is not a complete high-energy electroweak cross section throughout the entire stated range.
Efficiency normalization. Read the original PDF plot as a continuum of about and a maximum around at . The supplied conversion constant gives
Hence near the peak,
If multiplies the displayed cross sections, the required scale is
Here is the implied detection efficiency under the common-scale assumption. Using continuum readings – gives approximately –, illustrating the graph-reading precision. The nearby open-bottom threshold also limits the precision of the ideal continuum approximation.
Resonance quantum numbers and widths. Direct formation through a single photon requires
These are the quantum numbers of the electromagnetic current. This concerns single-photon formation; associated production or two-photon reactions need not obey this restriction.
For a Breit-Wigner resonance cross section, subtract the continuum before finding the full width at half maximum. The plotted excess is about , so its half-height occurs at a total displayed cross section near , not . The crossings are roughly and , with a separation around –. Thus
For spin-one formation from unpolarized electrons and positrons, and . With , the peak resonance cross section and branching fractions relation becomes
The corrected excess is . Taking as a representative graph reading gives
Thus the plot supports
The width and electronic partial width are estimates from the printed graph under the prescribed Breit-Wigner and efficiency assumptions; these are not precision world-average values. A scale correction changes the peak height and inferred partial width, but not the peak position or width.
Why the lower Upsilon states are narrow. In electromagnetic production and decay of vector quarkonium, every state can annihilate as . The electronic width probes the wavefunction at the origin, . The radial overlaps differ, but the same short-distance annihilation mechanism operates, so electronic widths remain of comparable order.
The hadronic mechanism changes sharply at . The first three Upsilon mesons lie below it and cannot produce an on-shell B-meson pair. Their leading purely strong annihilation channel is , followed by hadronization: a color-singlet pair cannot annihilate to one gluon, and two gluons cannot match its negative charge-conjugation parity. Three-gluon annihilation carries three powers of in the rate and requires quark-antiquark overlap. Electromagnetic annihilation to light quarks and, for excited states, transitions to lower bottomonium plus pions or photons also contribute.
The fourth state lies above the pair threshold and can instead undergo through light-pair creation, keeping the b and anti-b in separate outgoing mesons. The final panel depicts this quark-line topology; it represents a nonperturbative strong decay, not a free on-shell b quark radiating a real gluon. For two pseudoscalar B mesons the pair is in a P wave to give . Although close to threshold, this open strong channel has a much larger rate than the annihilation channels. Consequently a total width of tens of MeV can coexist with an electronic width below a keV; below threshold the total widths are smaller by hundreds or more, while their electronic widths remain of similar order.
Hyperfine estimate. For two spin-half constituents,
The meson hyperfine mass splitting therefore raises the vector mass above the pseudoscalar mass. To infer its coefficient from only one measured mass, adopt the intended simple constituent ansatz that the spin-independent baseline is . Then
For the spin-one partner,
The estimate is the same for and with the supplied equal light-quark masses. The splitting is below , so the strong decay is closed. An electromagnetic transition is open and dominates over weak decay:
It preserves electric charge, giving the corresponding charged or neutral B meson; the photon energy is approximately , slightly reduced by recoil. Strictly, the statement that the mass merely includes a spin-spin term does not fix : an independent spin-independent binding energy would leave the baseline unknown. The numerical prediction above explicitly uses the additive constituent-mass approximation.