research29 min studyLesson 635 in Physics

Rare flavor decays and effective Hamiltonians

Organize loop-suppressed flavor amplitudes with Wilson coefficients, GIM cancellation, and matched hadronic matrix elements.

On this page

What you’ll understand

  • State an effective-Hamiltonian normalization and operator dimensions.
  • Derive the cancellation of a flavor-independent loop contribution.
  • Identify coefficient degeneracies and hadronic limitations in rate interpretation.

A rare decay can probe short-distance interactions whose particles are too heavy to produce directly. Sensitivity comes from comparing a small Standard Model amplitude with well-controlled possible modifications. Effective Hamiltonians organize that comparison by separating local operators from their coefficients and hadronic matrix elements. Rarity alone does not guarantee theoretical cleanliness. We will derive a flavor cancellation, fix one coefficient convention, and examine why a precisely measured rate can still admit several microscopic interpretations. Natural units are used throughout.

Integrate out heavy fields without discarding their effects

At energies well below the W mass, a process can be expanded in local operators built from the remaining fields. Their Wilson coefficients encode heavy-particle effects, while matrix elements describe the chosen external hadrons. For an illustrative bottom-to-strange transition choose Heff=−(4GF/2)VtbVts∗∑iCiOi+h.c.\mathcal H_{\mathrm{eff}}=-(4G_F/\sqrt2)V_{tb}V_{ts}^*\sum_i C_iO_i+\mathrm{h.c.} as an operator-basis normalization. The displayed factorization is conventional; CKM ratios for additional contributions may remain inside coefficients. A number quoted without its operator definitions and extracted prefactors is not a portable observable.

For example, define the electromagnetic dipole operator using the running bottom mass, a chirality projector, and the photon field strength. Each quark field has mass dimension three halves, the field strength dimension two, and the inserted mass dimension one, giving operator dimension six. The Fermi constant has dimension minus two, so the Hamiltonian density has dimension four with dimensionless coefficients. Including the mass or loop factor in a different place changes the numerical coefficient convention without changing the underlying decay amplitude.

Heff=−4GF2VtbVts∗∑iCi(μ)Oi(μ)+h.c.,O7=e mb(μ)16π2sˉσμνPRb Fμν,PR=1+γ52\mathcal H_{\mathrm{eff}}=-\frac{4G_F}{\sqrt2}V_{tb}V_{ts}^*\sum_i C_i(\mu)O_i(\mu)+\mathrm{h.c.},\qquad O_7=\frac{e\,m_b(\mu)}{16\pi^2}\bar s\sigma^{\mu\nu}P_R b\,F_{\mu\nu},\qquad P_R=\frac{1+\gamma_5}{2}
This declares the dipole normalization; sigma-mu-nu is i[gamma-mu,gamma-nu]/2. GF has inverse-energy-squared units, O7 has energy dimension six, and C7 is dimensionless. It is not the full operator basis for a measured decay.↔ Scroll wide equations horizontally.

Unitarity removes flavor-blind loop pieces

Standard Model flavor-changing neutral currents have no tree-level neutral gauge coupling between distinct same-charge mass eigenstates. A loop transition can instead contain ∑i=u,c,tVibVis∗f(mi2/mW2)\sum_{i=u,c,t}V_{ib}V_{is}^*f(m_i^2/m_W^2). For different external flavors, column orthogonality makes the CKM sum zero. Any part of f independent of the internal flavor therefore cancels. If all internal masses were equal, the entire displayed contribution would vanish. This GIM mechanism explains a structural suppression, not a universal statement that every rare process is numerically dominated by the same internal quark.

Subtracting an arbitrary common value of the loop function exposes the cancellation explicitly and often clarifies which mass differences matter. Heavy-particle behavior then depends on the specific operator and symmetry structure; it cannot be inferred solely from dimensional counting. A full prediction also requires matching at a high scale and renormalization-group evolution to a hadronic scale. Operators generally mix under this evolution. Keeping one coefficient while neglecting its coupled matrix elements can spoil the cancellation of artificial scale dependence and create a misleading apparent effect.

∑iVibVis∗=0,∑iVibVis∗f(xi)=∑iVibVis∗[f(xi)−f(x0)],A∝∑iCi(μ)⟨f∣Oi(μ)∣B⟩\sum_iV_{ib}V_{is}^*=0,\qquad \sum_iV_{ib}V_{is}^*f(x_i)=\sum_iV_{ib}V_{is}^*[f(x_i)-f(x_0)],\qquad \mathcal A\propto\sum_i C_i(\mu)\langle f|O_i(\mu)|B\rangle
The first identity assumes a unitary three-family CKM matrix and b distinct from s. The arbitrary reference x0 cancels from the exact sum. Consistent coefficient and matrix-element scale and scheme dependence cancel in the physical amplitude.↔ Scroll wide equations horizontally.

Separate an informative constraint from a unique explanation

The rate is a squared coherent amplitude, so changing the sign or phase of one contribution may leave a rate unchanged. Angular observables, CP asymmetries, other decay modes, and spectral shapes can help resolve these degeneracies. Long-distance hadronic effects may imitate some short-distance coefficient shifts, particularly when resonances or nonlocal weak-current insertions contribute. Quantifying that ambiguity is part of the prediction. A clean mathematical coefficient fit becomes a new-interaction claim only after experimental and theoretical uncertainties are propagated consistently.

A heavy scale can sometimes be bounded through a coefficient proportional to coupling squared divided by mass squared, but an unknown coupling prevents a unique mass inference. Correlated operators, interference, and matching assumptions further affect that translation. Our hand example deliberately holds all matrix elements and phase space fixed and retains one real coefficient. It illustrates sensitivity and a sign degeneracy, not a numerical bound on a particular particle model. This distinction keeps an effective description useful without assigning it information that the chosen observable never supplied.

A coefficient shift can lower a rate

In a one-operator toy decay with unchanged matrix element and phase space, take real CSM=−0.30C_{\mathrm{SM}}=-0.30 and ΔC=+0.06\Delta C=+0.06. Find the rate ratio to the Standard Model toy rate and a different real shift giving the same ratio.

  1. The new coefficient is C=−0.30+0.06=−0.24C=-0.30+0.06=-0.24. Add amplitudes before squaring; adding the squared shifts would discard interference.

  2. The ratio is ∣−0.24∣2/∣−0.30∣2=(4/5)2=16/25=0.64|-0.24|^2/|-0.30|^2=(4/5)^2=16/25=0.64. A positive shift has lowered the magnitude because it partially cancels a negative reference coefficient.

  3. The opposite final coefficient, +0.24, has the same squared magnitude. It requires ΔC=0.24−(−0.30)=+0.54\Delta C=0.24-(-0.30)=+0.54, exhibiting a rate-only sign degeneracy.

The result

The toy rate becomes 64 percent of its reference value; shifts +0.06 and +0.54 give that same ratio through opposite final coefficient signs.

Explore the model

Interactive lab

Interactive lab · waves

A phase affects the squared coherent sum

Reset and set BOTH Phase difference = 90° and Amplitude ratio = 1; predict sum amplitude √2. For a fresh comparison, Reset and set Phase difference = 0° and Amplitude ratio = 1; predict sum amplitude 2. The corresponding squared-amplitude factors are two and four.

A phase affects the squared coherent sumTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe lab adds two equal-frequency sine waves. It supports coherent-amplitude reasoning but computes no Wilson coefficients, loop functions, hadronic uncertainties, or rare-decay branching fractions.

Check your understanding

Why does a flavor-independent loop-function term cancel in a unitary three-family b-to-s amplitude?

Try it yourself

Use common amplitude units with a short-distance contribution 2 and a long-distance contribution i. An illustrative new interaction changes only the short-distance contribution from 2 to 1. Find the old and new squared amplitudes and their ratio.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The relative phase is ninety degrees, so compute the modulus of each complete complex sum.

Compare with the solution

Initially the amplitude is 2+i2+i, whose squared magnitude is five. After the shift it is 1+i1+i, whose squared magnitude is two. The new-to-old rate ratio is therefore 2/52/5 with fixed phase space. Squaring only the short-distance piece would instead produce one quarter and miss the unchanged long-distance contribution. A realistic interpretation must justify both that contribution's magnitude and phase and the assumption that the proposed interaction leaves it unchanged.

Take it with you

  • Wilson coefficients are meaningful only with their normalization, scale, scheme, and operator basis.
  • GIM cancellation removes flavor-independent loop terms through CKM unitarity.
  • One rare rate can constrain amplitudes while leaving sign, phase, and hadronic ambiguities.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.