Quark model calculation of $η\to l^+ l^-$ to all orders in the bound state relative momentum

Abstract

The electromagnetic box diagram for the leptonic decays of pseudoscalar mesons in the quark model is evaluated to all orders in ${\bf p} / m_q$, where ${\bf p}$ is the relative three-momentum of the quark-antiquark pair and $m_q$ is the quark mass. We compute $B_P \equiv Γ(η\to l^+ l^-) / Γ(η\to γγ)$ using a popular nonrelativistic (NR) harmonic oscillator wave function, and with a relativistic momentum space wave function that we derive from the MIT bag model. We also compare with a calculation in the limit of extreme NR binding due to Bergström. Numerical calculations of $B_P$ using these three parameterizations of the wave function agree to within a few percent over a wide kinematical range. We find that the quark model leads in a natural way to a negligible value for the ratio of dispersive to absorptive parts of the electromagnetic amplitude for $η\to μ^+ μ^-$ (unitary bound). However we find substantial deviations from the unitary bound in other kinematical regions, such as $η,π^0 \to e^+ e^-$. These quark models yield $B(η\to μ^+μ^-) \approx 4.3 \times 10^{-6}$, within errors of the recent SATURNE measurement of $5.1 \pm 0.8 \times 10^{-6}$, $B(η\to e^+ e^-) \approx 6.3 \times 10^{-9}$, and $B(π^0 \to e^+ e^-) \approx 1.0 \times 10^{-7}$.