Radiative leptonic decays on the lattice

Radiative leptonic decays of pseudoscalar mesons probe both the weak interaction and the hadronic structure in useful ways. Adding a sufficiently energetic photon to the final state can actually increase the branching fraction [2], as it removes the helicity suppression. Perhaps the most interesting example is $B^{-}\to\ell^{-}\bar{\nu}\,\gamma$, shown in Fig. 1 (left). For large $E_{\gamma}^{(0)}$, this process is the cleanest probe of the first inverse moment of the $B$-meson light-cone distribution amplitude, $1/\lambda_{B}=\int_{0}^{\infty}\frac{\Phi_{B^{+}}(\omega)}{\omega}\>d\omega$, an important input in QCD-factorization predictions for nonleptonic $B$ decays that is presently poorly determined [3, 4, 5, 6, 7, 8, 9]. A recent search for this decay by Belle gave an upper limit $\mathcal{B}(B^{-}\to\ell^{-}\bar{\nu}\gamma,\>E_{\gamma}^{(0)}>1\>{\rm GeV})<3.0\times 10^{-6}$, close to the Standard-Model expectation [10]. Lattice QCD results for the $B^{-}\to\ell^{-}\bar{\nu}\gamma$ form factors could be used to constrain $\lambda_{B}$. Also very interesting are the flavor-changing neutral-current decays $B^{0}\to\ell^{+}\ell^{-}\gamma$ and $B_{s}\to\ell^{+}\ell^{-}\gamma$ (shown in Fig. 1, right). While the purely leptonic decays are sensitive to $C_{10,S,P}-C^{\prime}_{10,S,P}$ only, the radiative leptonic decays probe all Wilson coefficients in the weak effective Hamiltonian, including $C_{9}$, in which global fits of experimental results for other $b\to s\ell^{+}\ell^{-}$ decays indicate a deviation from the Standard Model that violates lepton flavor universality (LFU) [1]. Since the radiative leptonic decays are not helicity-suppressed, they are well-suited for testing LFU with light leptons [11, 12]. For the charmed-meson radiative leptonic decays $D^{+}\to e^{+}\nu\gamma$ and $D_{s}^{+}\to e^{+}\nu\gamma$, the BESIII collaboration has reported upper limits on the branching fractions with $E_{\gamma}^{(0)}>10\>{\rm MeV}$ of $3.0\times 10^{-5}$ and $1.3\times 10^{-4}$, respectively [13, 14]. Finally, in contrast to the heavy-meson decays, there are already precise measurements of the differential branching fractions of $K^{-}\to e^{-}\bar{\nu}\gamma$, $K^{-}\to\mu^{-}\bar{\nu}\gamma$, $\pi^{-}\to e^{-}\bar{\nu}\gamma$, and $\pi^{-}\to\mu^{-}\bar{\nu}\gamma$, as reviewed in Ref. [15]. These decay modes can therefore be used to test the lattice QCD methods.

In the following, we show how radiative leptonic decays can be calculated on a Euclidean lattice, and we present early numerical results. One of us previously reported on this project at the Lattice 2018 conference [16]. At Lattice 2019, radiative leptonic decays were also discussed by G. Martinelli [17].

To define the form factors for charged-current radiative leptonic decays of pseudoscalar mesons, we use the notation for $B^{-}\to\ell^{-}\bar{\nu}\,\gamma$. The quark electromagnetic and weak currents are given by $J_{\mu}=\sum_{q}e_{q}\,\bar{q}\gamma_{\mu}q$ and $J^{\rm weak}_{\nu}=\bar{u}\gamma_{\nu}(1-\gamma_{5})b$. The decay amplitude depends on the hadronic tensor, which is defined as

in Minkowski space. Throughout this work, we asssume that the photon is real, i.e., $p_{\gamma}^{2}=0$. The hadronic tensor can be decomposed as [7]

where p_{B}=m_{B}{\textsl{v}} and the $(p_{\gamma})_{\mu}\text{-terms}$ will disappear when contracting with the photon polarization vector. The form factors $F_{V}$ and $F_{A}$ are functions of the photon energy in the $B$-meson rest frame, E_{\gamma}^{(0)}=p_{\gamma}\cdot{\textsl{v}}=(m_{B}^{2}-q^{2})/(2m_{B}). Also appearing in Eq. (2) is the $B$-meson decay constant $f_{B}$.

To prepare for the discussion in the next section, it is useful to write down the spectral representation of $T_{\mu\nu}$ in Minkowski space for the two different time orderings of the currents. By inserting complete sets of energy/momentum eigenstates and performing the time integrals, we find

(in infinite volume, the sums over $n$ and $m$ include integrals over the continuous spectrum of multi-particle states).

In this section, we show that $T_{\mu\nu}$ can be extracted from the Euclidean three-point function

where $\phi_{B}\sim\bar{u}\gamma_{5}b$ is an interpolating field for the $B$ meson, and $t$, $t_{B}$ now denote the Euclidean time. We define the integrals

with a finite integration range $T$. Here we take $t_{B}$ to be large and negative (with $t_{B}<-T$), such that ground-state saturation is achieved for the $B$ meson. Inserting again complete sets of energy/momentum eigenstates, we find, for the first time ordering,

The sum over states in Eq. (7) differs from the sum in Eq. (3) by the factor in the last line. However, the exponential $e^{-(E_{\gamma}+E_{n,(\mathbf{p}_{B}-\mathbf{p}_{\gamma})}-E_{B})T}$ will vanish for large $T$ if $E_{\gamma}+E_{n,(\mathbf{p}_{B}-\mathbf{p}_{\gamma})}>E_{B}$. Because the states $|n(\mathbf{p}_{B}-\mathbf{p}_{\gamma})\rangle$ have the same quark-flavor quantum numbers as the $B$ meson, we have $E_{n,(\mathbf{p}_{B}-\mathbf{p}_{\gamma})}\geq E_{B,(\mathbf{p}_{B}-\mathbf{p}_{\gamma})}=\sqrt{m_{B}^{2}+(\mathbf{p}_{B}-\mathbf{p}_{\gamma})^{2}}$. Thus, we need $\displaystyle\sqrt{\mathbf{p}_{\gamma}^{2}}+\sqrt{m_{B}^{2}+(\mathbf{p}_{B}-\mathbf{p}_{\gamma})^{2}}>\sqrt{m_{B}^{2}+\mathbf{p}_{B}^{2}}$. This is in fact always true if $\mathbf{p}_{\gamma}\neq 0$.

For the other time ordering, we find

The unwanted exponential $e^{(E_{\gamma}-E_{m,\mathbf{p}_{\gamma}})T}$ in the last line goes to zero for large $T$ if $E_{m,\mathbf{p}_{\gamma}}>E_{\gamma}$. Because the states $|m(\mathbf{p}_{\gamma})\rangle$ are hadronic and have nonzero masses, their energies are larger than the energy of a photon with the same spatial momentum, showing that this condition is also always satisfied. In summary, for $\mathbf{p}_{\gamma}\neq 0$,

where $I_{\mu\nu}$ is the integral from $-T$ to $T$. The energy $E_{B}$ and the overlap factor $\langle B(\mathbf{p}_{B})|\phi_{B}^{\dagger}(0)|0\rangle$ can be obtained from the two-point function $\int d^{3}x\>\>e^{-i\mathbf{p}_{B}\cdot\mathbf{x}}\langle\phi_{B}(t,\mathbf{x})\>\phi_{B}^{\dagger}(0)\rangle$.

Note that similar nonlocal matrix elements appear in processes with two photons, whose lattice calculation has been discussed, for example, in Refs. [18, 19, 20].

In this section, we present some early numerical results for the $D_{s}^{+}\to\ell^{+}\nu\gamma$ and $K^{-}\to\ell^{-}\bar{\nu}\gamma$ form factors. These results are from only 25 configurations of the “24I” RBC/UKQCD ensemble [22] with $2+1$ flavors of domain-wall fermions and the Iwasaki gauge action, with $a^{-1}=1.785(5)$ GeV and $m_{\pi}=340(1)$ MeV. For the light and strange valence quarks, we use the same domain-wall action as in Ref. [22]. The valence charm quark is implemented with a Möbius domain-wall action with stout-smeared gauge links ($N=3$, $\rho=0.1$), $L_{5}/a=12$, $aM_{5}=1.0$, $am_{f}=0.6$ [23], which approximately corresponds to the physical charm-quark mass. We use local currents with “mostly nonperturbative” renormalization. Gaussian smearing is performed for the lighter quark in the meson interpolating field. We start with a $\mathbb{Z}_{2}$ random-wall source at the time slice of the weak current (denoted as time “0” here) and perform sequential inversions through the meson interpolating field; disconnected diagrams are presently neglected. All-mode averaging [24] with 16 sloppy and 1 exact samples per configuration is employed; the 16 sloppy samples correspond to 16 different starting time slices. Our initial calculations used $\mathbf{p}_{K/D_{s}}=0$ and $\mathbf{p}^{2}_{\gamma}\in\{1,2,3,4,5\}\left(\frac{2\pi}{L}\right)^{2}$.

Figure 2 shows examples of the $D_{s}^{+}\to\ell^{+}\nu\gamma$ three-point functions. Multiplying by $e^{E_{\gamma}t}$ and summing over $t$ gives $T_{\mu\nu}$ for sufficiently large summation range $T$. The form factors $F_{V}$ and $F_{A}$ extracted from $T_{\mu\nu}$ (at the lowest photon momentum) are shown as a function of $T$ in Fig. 3. The results plateau at approximately $T/a=8$. We also extracted the meson decay constants from the $v^{\mu}v^{\nu}$ term in $T_{\mu\nu}$. As can be seen in Fig. 4, the results agree with the known values, which is a valuable test of our calculation. Finally, Fig. 5 shows the form factors $F_{V}$ and $F_{A}$ as a function of the photon energy. Note that, with our current choice of momenta, all of the photon energies are above the physical region for $K^{-}\to\ell^{-}\bar{\nu}\gamma$. The results for $F_{A}$ are dominated by the point-like contribution equal to $-e_{\ell}f_{K/D_{s}}/E_{\gamma}^{(0)}$.

We have shown that the form factors describing radiative leptonic decays can be calculated on the lattice; even though they involve a nonlocal matrix element, the use of imaginary time poses no difficulty in this case. The early results shown here for $D_{s}^{+}\to\ell^{+}\nu\gamma$ and $K^{-}\to\ell^{-}\bar{\nu}\gamma$ cover photon energies from approximately 0.5 to 1 GeV. For $K^{-}\to\ell^{-}\bar{\nu}\gamma$ we need to reach lower photon energies to compare with experiment; this can be achieved by using moving frames (i.e., nonzero $\mathbf{p}_{K}$) and/or a larger volume. To study the $B_{(s)}$ radiative leptonic decays with the domain-wall action for the heavy quark, we will need to extrapolate in the mass. We are also considering calculations directly at the physical $b$-quark mass using the “relativistic heavy-quark action” [25], but, because this action is only on-shell improved, additional steps are likely needed to remove unphysical behavior occurring when the electromagnetic and weak currents get close to each other.

Acknowledgments: We thank the RBC and UKQCD Collaborations for providing the gauge-field configurations. C.K. and S.M. are supported by the US DOE, Office of Science, Office of HEP under award number DE-SC0009913. S.M. is also supported by the RIKEN BNL Research Center. A.S. and C.L. are supported in part by US DOE Contract No. DESC0012704(BNL). During a part of this work, C.L. was also supported by a DOE Office of Science Early Career Award. This work used resources at TACC that are part of XSEDE, supported by NSF grant number ACI-1548562.