research29 min studyLesson 633 in Physics

Neutral meson mixing and time evolution

Evolve an unstable two-flavor meson system with declared eigenstate signs and separate survival probability from conditional flavor oscillation.

On this page

What you’ll understand

  • Construct flavor evolution from two decaying propagation eigenstates.
  • Keep the mass-eigenstate and g-minus signs consistent.
  • Calculate instantaneous and time-integrated mixing probabilities.

A neutral meson and its antiparticle can communicate through weak interactions and common decay channels. Their propagation eigenstates are consequently superpositions of flavor states with different masses and widths. The reduced two-state system is unstable: probability leaves it when the meson decays. Treating its evolution as an ordinary normalized qubit would discard essential physics. We will derive the signs directly from our eigenstate definitions and then make a controlled equal-width approximation for hand calculations.

An effective Hamiltonian includes loss to decay channels

On the meson-antimeson subspace use Heff=M−iΓ/2H_{\mathrm{eff}}=M-i\Gamma/2, with M and Gamma Hermitian matrices. This is the usual exponential-decay approximation after eliminating the decay products, not an assertion that fundamental quantum mechanics is nonunitary. CPT invariance makes the two diagonal entries equal. Off-diagonal dispersive and absorptive terms describe virtual conversion and shared on-shell decay channels. The eigenvalues can be written mass minus i times width over two. Natural units are used until numerical time units are introduced explicitly.

Define ∣ML⟩=p∣M⟩+q∣Mˉ⟩|M_L\rangle=p|M\rangle+q|\bar M\rangle and ∣MH⟩=p∣M⟩−q∣Mˉ⟩|M_H\rangle=p|M\rangle-q|\bar M\rangle, normalized by ∣p∣2+∣q∣2=1|p|^2+|q|^2=1. Here H labels the larger mass and L the smaller, so Δm=mH−mL>0\Delta m=m_H-m_L>0. The ratio q over p depends on flavor-state phases, whereas its magnitude does not. When CP is violated in mixing these propagation eigenstates need not be orthogonal in the reduced description. The general formalism handles that feature; our probability example additionally sets the ratio's magnitude to one.

Heff=M−i2Γ,eL,H(t)=e−imL,Ht−ΓL,Ht/2,g±(t)=eL(t)±eH(t)2H_{\mathrm{eff}}=M-\frac{i}{2}\Gamma,\qquad e_{L,H}(t)=e^{-im_{L,H}t-\Gamma_{L,H}t/2},\qquad g_\pm(t)=\frac{e_L(t)\pm e_H(t)}2
Time is proper time of the propagating meson. In natural units masses and widths have energy units; with ps as time units, use mass frequency delta m divided by hbar and decay rate Gamma divided by hbar.↔ Scroll wide equations horizontally.

Invert the eigenstate definitions before evolving

Adding and subtracting the eigenstates expresses the initial flavor M as their equal-amplitude combination divided by twice p. Evolve each component and transform back. Our definition of g minus is e L minus e H over two, so the converted component carries a plus sign multiplying (q/p)g−(q/p)g_-. Some references define the opposite g minus and compensate by placing a minus sign in the flavor evolution. Either convention works; mixing one convention's definition with the other's evolution changes interference signs and is incorrect.

For equal widths and ∣q/p∣=1|q/p|=1, the two flavor probabilities share an exponential survival factor. The unmixed component has cosine squared and the mixed component sine squared of half the mass phase. Their sum is the survival probability, not one. Dividing by that sum gives probabilities conditioned on the meson still being present. This distinction matters experimentally because flavor-tagged decay rates include both the oscillation and the decay density. A large converted fraction among surviving mesons need not mean a large number of surviving events.

∣M(t)⟩=g+(t)∣M⟩+qpg−(t)∣Mˉ⟩,PM→Mˉ(t)=e−γtsin⁡2ωt2,PM→M(t)=e−γtcos⁡2ωt2|M(t)\rangle=g_+(t)|M\rangle+\frac qp g_-(t)|\bar M\rangle,\qquad P_{M\to\bar M}(t)=e^{-\gamma t}\sin^2\frac{\omega t}{2},\qquad P_{M\to M}(t)=e^{-\gamma t}\cos^2\frac{\omega t}{2}
The last two formulas assume equal widths and no CP violation in mixing. Here omega = delta m/hbar and gamma = common width/hbar are frequencies in the same inverse-time units as the chosen t.↔ Scroll wide equations horizontally.

Integrating decay time changes the observable

For flavor-specific decays with equal tagging efficiencies and the same partial-width normalization, integrate the mixed decay density over time. In the equal-width, CP-conserving limit its fraction among all such decays is χ=x2/[2(1+x2)]\chi=x^2/[2(1+x^2)], where x is the mass-oscillation frequency divided by the decay rate. This follows by integrating gamma times the two probabilities; the sum integrates to one. Very slow mixing gives a small fraction, while rapid unresolved mixing tends to one half rather than complete conversion.

Mixing itself is not evidence for CP violation. A real, symmetric effective mixing interaction can already produce mass splitting and oscillation. CP tests compare particle and antiparticle behavior or their interference with decay amplitudes, which the next lesson develops. Unequal widths, detector time acceptance, imperfect initial tagging, and experimental resolution modify observable time distributions. The simple integrated fraction should therefore be used only after stating its assumptions. A laboratory decay time must also be converted to proper time when the meson is moving relativistically.

Survival and flavor after one picosecond

Assume equal widths, ∣q/p∣=1|q/p|=1, decay rate γ=ln⁡2 ps−1\gamma=\ln2\,\mathrm{ps^{-1}}, and mass frequency ω=π/2 ps−1\omega=\pi/2\,\mathrm{ps^{-1}}. Starting with M, find survival, mixed, and unmixed probabilities at proper time 1 ps.

  1. The survival factor is e−γt=e−ln⁡2=1/2e^{-\gamma t}=e^{-\ln2}=1/2. Half the original ensemble has decayed before this proper time.

  2. The half phase is ωt/2=π/4\omega t/2=\pi/4, so sine squared and cosine squared both equal one half. Multiply each by the survival factor rather than renormalizing the initial ensemble.

  3. The mixed and unmixed probabilities are each one quarter. Their sum is one half, and the conditional mixed fraction among survivors is (1/4)/(1/2)=1/2(1/4)/(1/2)=1/2.

The result

Survival is 1/2; mixed and unmixed probabilities are each 1/4 relative to the initial ensemble, with a conditional mixed fraction of 1/2.

Explore the model

Interactive lab

Interactive lab · quantum

Normalized weights omit meson decay

Reset and set BOTH Mixing angle θ = 60° and Relative phase φ = 0°; predict P(1) = 0.25. For a fresh comparison, Reset and set θ = 120° and φ = 0°; predict P(1) = 0.75. In each run Play changes relative phase but leaves these basis probabilities fixed.

Normalized weights omit meson decayA normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Current readings: Probability |0⟩: 71.1%; Probability |1⟩: 28.9%; Bloch (x, y, z): 0.78, 0.45, 0.42. All parameters are adjustable with the labeled controls below.xyz|0⟩|1⟩Basis probabilitiesP(0)71.1%P(1)28.9%
Bloch vector / P(0)P(1)Phase latitude
Probability |0⟩
71.1%
Probability |1⟩
28.9%
Bloch (x, y, z)
0.78, 0.45, 0.42
t = 0 s

A normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. 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 forThis pure-state instrument preserves total probability and illustrates normalized two-state weights. It has no decay width, proper-time mixing Hamiltonian, or meson survival channel; its playback is not the displayed meson evolution.

Check your understanding

In the equal-width CP-conserving example, what is P(M to M) plus P(M to anti-M) before conditioning on survival?

Try it yourself

In the same equal-width approximation, a meson has mass frequency ω=4 ps−1\omega=4\,\mathrm{ps^{-1}} and decay rate γ=2 ps−1\gamma=2\,\mathrm{ps^{-1}}. Find the time-integrated mixed fraction for flavor-specific decays, and compare it with the rapid-mixing limit.

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

Give me a hint

Use x=ω/γx=\omega/\gamma in the integrated fraction, not either dimensional frequency by itself.

Compare with the solution

The dimensionless ratio is x = 2. Therefore χ=4/[2(1+4)]=2/5\chi=4/[2(1+4)]=2/5. This is a fraction of the time-integrated flavor-specific decay sample under equal partial-width and efficiency assumptions, not the instantaneous survival probability. As x becomes arbitrarily large the fraction approaches one half. The value two fifths is below that limit because an appreciable subset decays before completing a large oscillation phase.

Take it with you

  • An unstable meson subspace evolves with a non-Hermitian effective Hamiltonian.
  • Consistent eigenstate and g-minus definitions determine interference signs.
  • Mixing, survival, and CP violation are distinct observables.

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.