Search

We determine, by means of lattice QCD calculations, the local form factors describing the $B_{s}\to \mu^{+}\mu^{-}\gamma$ decay. For this analysis we make use of the gauge configurations produced by the ETM Collaboration with $N_{f}=2+1+1$ flavour of Wilson-Clover twisted-mass fermions at maximal twist. To obtain the $B_{s}$ meson form-factors, we perform simulations for several heavy-strange meson masses $m_{H_{s}}$ in the range $m_{H_{s}} \in [ m_{D_{s}}, 2 m_{D_{s}} ]$, and extrapolate to the physical $B_{s}$ meson point $m_{B_{s}}\simeq 5.367~{\rm GeV}$ making use of the HQET scaling laws. We cover the region of large di-muon invariant masses $\sqrt{q^{2}} > 4.16\,{\rm GeV}$, and use our results to determine the branching fraction for $B_{s}\to \mu^{+}\mu^{-}\gamma$, which has been recently measured by LHCb in the region $\sqrt{q^{2}} > 4.9\,{\rm GeV}$. The largest contribution to the uncertainty in the partial branching fractions at values of $\sqrt{q^{2}} < 4.8\,{\rm GeV}$ is now due to resonance and other long-distance effects, including those from "charming penguins", which we estimate by summing over the contributions from the $J_P=1^-$ charmonium resonances., Comment: 42 page, 25 figures, 6 tables

We present a lattice determination of the inclusive decay rate of the process $\tau\mapsto X_{us} \nu_\tau$ in which the $\tau$ lepton decays into a generic hadronic state $X_{us}$ with $\bar u s$ flavour quantum numbers. Our results have been obtained in $n_f=2+1+1$ iso-symmetric QCD with full non-perturbative accuracy, without any OPE approximation and, except for the presently missing long-distance isospin-breaking corrections, include a solid estimate of all sources of theoretical uncertainties. This has been possible by using the Hansen-Lupo-Tantalo method [1] that we have already successfully applied in [2] to compute the inclusive decay rate of the process $\tau\mapsto X_{ud} \nu_\tau$ in the $\bar u d$ flavour channel. By combining our first-principles theoretical results with the presently-available experimental data we extract the CKM matrix element $\vert V_{us}\vert$, the Cabibbo angle, with a $0.9$\% accuracy, dominated by the experimental error., Comment: 9 pages, 5 figures, 1 table. Version accepted by PRL, expanded technical discussion moved to an appendix, results unchanged

We compute the structure-dependent axial and vector form factors for the radiative leptonic decays $D_s\to \ell\nu_\ell\gamma$, where $\ell$ is a charged lepton, as functions of the energy of the photon in the rest frame of the $D_s$ meson. The computation is performed using gauge-field configurations with 2+1+1 sea-quark flavours generated by the European Twisted Mass Collaboration and the results have been extrapolated to the continuum limit. For the vector form factor we observe a very significant partial cancellation between the contributions from the emission of the photon from the strange quark and that from the charm quark. The results for the form factors are used to test the reliability of various Anz\"atze based on single-pole dominance and its extensions, and we present a simple parametrization of the form factors which fits our data very well and which can be used in future phenomenological analyses. Using the form factors we compute the differential decay rate and the branching ratio for the process $D_s\to e\nu_e\gamma$ as a function of the lower cut-off on the photon energy. With a cut-off of 10 MeV for example, we find a branching ratio of Br$(E_\gamma>10\,\mathrm{MeV})=4.4(3)\times 10^{-6}$ which, unlike some model calculations, is consistent with the upper bound from the BESIII experiment Br$(E_\gamma>10\,\mathrm{MeV})<1.3\times 10^{-4}$ at 90% confidence level. Even for photon energies as low as 10 MeV, the decay $D_s\to e\nu_e\gamma$ is dominated by the structure-dependent contribution to the amplitude (unlike the decays with $\ell=\mu$ or $\tau$), confirming its value in searches for hypothetical new physics as well as in determining the Cabibbo-Kobayashi-Maskawa (CKM) parameters at $O(\alpha_\mathrm{em})$, where $\alpha_{\mathrm{em}}$ is the fine-structure constant., Comment: 31 pages, 14 figures

We compute the structure-dependent axial and vector form factors for the radiative leptonic decays $D_s\to \ell\nu_\ell\gamma$, where $\ell$ is a charged lepton, as functions of the energy of the photon in the rest frame of the $D_s$ meson. The computation is performed using gauge-field configurations with 2+1+1 sea-quark flavours generated by the European Twisted Mass Collaboration and the results have been extrapolated to the continuum limit. For the vector form factor we observe a very significant partial cancellation between the contributions from the emission of the photon from the strange quark and that from the charm quark. The results for the form factors are used to test the reliability of various Anz\"atze based on single-pole dominance and its extensions, and we present a simple parametrization of the form factors which fits our data very well and which can be used in future phenomenological analyses. Using the form factors we compute the differential decay rate and the branching ratio for the process $D_s\to e\nu_e\gamma$ as a function of the lower cut-off on the photon energy. With a cut-off of 10 MeV for example, we find a branching ratio of Br$(E_\gamma>10\,\mathrm{MeV})=4.4(3)\times 10^{-6}$ which, unlike some model calculations, is consistent with the upper bound from the BESIII experiment Br$(E_\gamma>10\,\mathrm{MeV})<1.3\times 10^{-4}$ at 90% confidence level. Even for photon energies as low as 10 MeV, the decay $D_s\to e\nu_e\gamma$ is dominated by the structure-dependent contribution to the amplitude (unlike the decays with $\ell=\mu$ or $\tau$), confirming its value in searches for hypothetical new physics as well as in determining the Cabibbo-Kobayashi-Maskawa (CKM) parameters at $O(\alpha_\mathrm{em})$, where $\alpha_{\mathrm{em}}$ is the fine-structure constant., Comment: 31 pages, 14 figures

We compute the structure-dependent axial and vector form factors for the radiative leptonic decays $D_s\to \ell\nu_\ell\gamma$, where $\ell$ is a charged lepton, as functions of the energy of the photon in the rest frame of the $D_s$ meson. The computation is performed using gauge-field configurations with 2+1+1 sea-quark flavours generated by the European Twisted Mass Collaboration and the results have been extrapolated to the continuum limit. For the vector form factor we observe a very significant partial cancellation between the contributions from the emission of the photon from the strange quark and that from the charm quark. The results for the form factors are used to test the reliability of various Anz\"atze based on single-pole dominance and its extensions, and we present a simple parametrization of the form factors which fits our data very well and which can be used in future phenomenological analyses. Using the form factors we compute the differential decay rate and the branching ratio for the process $D_s\to e\nu_e\gamma$ as a function of the lower cut-off on the photon energy. With a cut-off of 10 MeV for example, we find a branching ratio of Br$(E_\gamma>10\,\mathrm{MeV})=4.4(3)\times 10^{-6}$ which, unlike some model calculations, is consistent with the upper bound from the BESIII experiment Br$(E_\gamma>10\,\mathrm{MeV})<1.3\times 10^{-4}$ at 90% confidence level. Even for photon energies as low as 10 MeV, the decay $D_s\to e\nu_e\gamma$ is dominated by the structure-dependent contribution to the amplitude (unlike the decays with $\ell=\mu$ or $\tau$), confirming its value in searches for hypothetical new physics as well as in determining the Cabibbo-Kobayashi-Maskawa (CKM) parameters at $O(\alpha_\mathrm{em})$, where $\alpha_{\mathrm{em}}$ is the fine-structure constant., Comment: 31 pages, 14 figures

We present a novel method to determine on the lattice both the real and imaginary parts of complex electroweak amplitudes involving two external currents and a single hadron or the QCD vacuum in the external states. The method is based on the spectral representation of the relevant time-dependent correlation functions and, by extending the range of applicability of other recent proposals built on the same techniques, overcomes the difficulties related to the analytic continuation from Minkowskian to Euclidean time, arising when intermediate states with energies smaller than the external states contribute to the amplitude. In its simplest form, the method relies on the standard $i \varepsilon$ prescription to regularize the Feynman integrals and at finite $\varepsilon$ it requires to verify the condition $1/L \ll \varepsilon \ll \Delta(E)$, where $L$ is the spatial extent of the lattice and, for any given energy $E$, $\Delta(E)$ represents the typical size of the interval around $E$ in which the hadronic amplitude is significantly varying. In order to illustrate the effectiveness of this approach in a realistic case, we apply the method to evaluate non-perturbatively the hadronic amplitude contributing to the radiative leptonic decay $D_s \to\ell\nu_\ell\,\gamma^*$, working for simplicity with a single lattice ensemble at fixed volume and lattice spacing., Comment: 22 pages, 8 figures

We study the radiative leptonic decays $P\to\ell\nu_\ell\,\ell^{\prime\,+}\ell^{\prime\,-}$, where $P$ is a pseudoscalar meson and $\ell$ and $\ell^\prime$ are charged leptons. In such decays the emitted photon is off-shell and, in addition to the "point-like" contribution in which the virtual photon is emitted either from the lepton or the meson treated as a point-like particle, four structure-dependent (SD) form factors contribute to the amplitude. We present a strategy for the extraction of the SD form factors and implement it in an exploratory lattice computation of the decay rates for the four channels of kaon decays ($\ell,\ell^\prime=e,\mu$). It is the SD form factors which describe the interaction between the virtual photon and the internal hadronic structure of the decaying meson, and in our procedure we separate the SD and point-like contributions to the amplitudes. We demonstrate that the form factors can be extracted with good precision and, in spite of the unphysical quark masses used in our simulation ($m_\pi\simeq 320\,$MeV and $m_K\simeq 530\,$MeV), the results for the decay rates are in reasonable semiquantitative agreement with experimental data (for the channels where these exist). Following this preparatory work, the emphasis of our future work will be on obtaining results at physical quark masses and on the control of the systematic uncertainties associated with discretisation and finite-volume errors., Comment: 30 pages, 47 figures

We present a lattice calculation of the charged/neutral pion mass difference $M_{\pi^{+}}-M_{\pi^{0}}$ at order $\mathcal{O}(\alpha_{em})$ using the gauge configurations produced by the Extended Twisted Mass Collaboration with $N_{f}=2+1+1$ dynamical quark flavours at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089~{\rm fm}$) and pion masses in the range $M_{\pi} \simeq 250-500~{\rm MeV}$. We employ the RM123 method and expand the path-integral around the isospin symmetric point at leading order in the electromagnetic coupling $\alpha_{em}$. Making use of the recently proposed RTM scheme, we evaluate the full $\mathcal{O}(\alpha_{em})$ contribution, with the inclusion of the disconnected diagram. At the physical point, after performing the continuum and infinite volume extrapolation, we obtain the value $M_{\pi^{+}}-M_{\pi^{0}}= 4.622~(95)~{\rm MeV}$ which is in good agreement with the experimental result $[ M_{\pi^{+}} - M_{\pi^{0}} ]^{exp.} = 4.5936(5)~{\rm MeV}$., Comment: 13 pages, 5 figures, 1 table

We present a preliminary non-perturbative lattice calculation of the form factors entering the processes $K\to \ell\,\nu_\ell\,\ell'^+\,\ell'^-$ and of the corresponding branching ratios. These form factors describe the interaction between the mediating virtual photon and the internal hadronic structure of the meson. By separating them from the point-like contribution to the matrix element we are able to isolate and reconstruct the structure-dependent contribution to the decay width. Our numerical analysis employs only one gauge ensemble and so it is affected by systematic uncertainties due to the missing continuum and physical point extrapolation. Despite this, we already find a reasonable agreement with the experimental data and with the next to leading order Chiral Perturbation Theory predictions. The method is general and can be applied to any pseudoscalar meson, though for heavier mesons the possibility of internal lighter states becomes problematic and still needs a proper study. A non-perturbative, model-independent lattice evaluation of these processes would allow further progress in the theoretical predictions of SM hadronic quantities and in the search of New Physics., Comment: 11 pages, 16 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{SD}$ and $a_{\mu}^W$, defined by the RBC/UKQCD Collaboration [1]. We employ a subset of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks, which are close to the physical point for the masses of all the dynamical flavors. The simulations are carried out at three values of the lattice spacing ranging from $\simeq 0.057$ to $\simeq 0.080$ fm with linear lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{SD}({\rm ETMC}) = 69.33\,(29) \cdot 10^{-10}$, which is consistent with the recent dispersive value $a_\mu^{SD}(e^+ e^-) = 68.4\,(5) \cdot 10^{-10}$ [2] within $\simeq 1.6 \, \sigma$. In the case of the intermediate window we get the value $a_\mu^W({\rm ETMC}) = 235.0\,(1.1) \cdot 10^{-10}$, which is consistent with the result $a_\mu^W({\rm BMW}) = 236.7\,(1.4) \cdot 10^{-10}$ [3] by the BMW collaboration as well as with the recent determination by the CLS/Mainz group of $a_\mu^W({\rm CLS}) = 237.30\,(1.46) \cdot 10^{-10}$ [4] at the $\sim 1.0 - 1.3 \, \sigma$ level. However, it is larger than the dispersive result $a_\mu^W(e^+ e^-) = 229.4\,(1.4) \cdot 10^{-10}$ [2] by $\simeq 3.1 \, \sigma$. The tension increases to $\simeq 4.2 \, \sigma$ if we average our ETMC result with the BMW and the CLS/Mainz ones. Our accurate lattice results in the short and intermediate windows hint at possible deviations of the $e^+ e^-$ cross section data with respect to Standard Model (SM) predictions distributed somewhere in the low (and possibly intermediate) energy regions, but not in the high energy region., Comment: 82 pages, 33 figures, 14 tables

We present a lattice calculation of the charged/neutral pion mass difference $M_{\pi^{+}}-M_{\pi^{0}}$ at order $\mathcal{O}(\alpha_{em})$ using the gauge configurations produced by the Extended Twisted Mass Collaboration with $N_{f}=2+1+1$ dynamical quark flavours at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089~{\rm fm}$) and pion masses in the range $M_{\pi} \simeq 250-500~{\rm MeV}$. We employ the RM123 method and expand the path-integral around the isospin symmetric point at leading order in the electromagnetic coupling $\alpha_{em}$. Making use of the recently proposed RTM scheme, we evaluate the full $\mathcal{O}(\alpha_{em})$ contribution, with the inclusion of the disconnected diagram. At the physical point, after performing the continuum and infinite volume extrapolation, we obtain the value $M_{\pi^{+}}-M_{\pi^{0}}= 4.622~(95)~{\rm MeV}$ which is in good agreement with the experimental result $[ M_{\pi^{+}} - M_{\pi^{0}} ]^{exp.} = 4.5936(5)~{\rm MeV}$., Comment: 13 pages, 5 figures, 1 table

We study the radiative leptonic decays $P\to\ell\nu_\ell\,\ell^{\prime\,+}\ell^{\prime\,-}$, where $P$ is a pseudoscalar meson and $\ell$ and $\ell^\prime$ are charged leptons. In such decays the emitted photon is off-shell and, in addition to the "point-like" contribution in which the virtual photon is emitted either from the lepton or the meson treated as a point-like particle, four structure-dependent (SD) form factors contribute to the amplitude. We present a strategy for the extraction of the SD form factors and implement it in an exploratory lattice computation of the decay rates for the four channels of kaon decays ($\ell,\ell^\prime=e,\mu$). It is the SD form factors which describe the interaction between the virtual photon and the internal hadronic structure of the decaying meson, and in our procedure we separate the SD and point-like contributions to the amplitudes. We demonstrate that the form factors can be extracted with good precision and, in spite of the unphysical quark masses used in our simulation ($m_\pi\simeq 320\,$MeV and $m_K\simeq 530\,$MeV), the results for the decay rates are in reasonable semiquantitative agreement with experimental data (for the channels where these exist). Following this preparatory work, the emphasis of our future work will be on obtaining results at physical quark masses and on the control of the systematic uncertainties associated with discretisation and finite-volume errors., Comment: 30 pages, 47 figures

We present a preliminary non-perturbative lattice calculation of the form factors entering the processes $K\to \ell\,\nu_\ell\,\ell'^+\,\ell'^-$ and of the corresponding branching ratios. These form factors describe the interaction between the mediating virtual photon and the internal hadronic structure of the meson. By separating them from the point-like contribution to the matrix element we are able to isolate and reconstruct the structure-dependent contribution to the decay width. Our numerical analysis employs only one gauge ensemble and so it is affected by systematic uncertainties due to the missing continuum and physical point extrapolation. Despite this, we already find a reasonable agreement with the experimental data and with the next to leading order Chiral Perturbation Theory predictions. The method is general and can be applied to any pseudoscalar meson, though for heavier mesons the possibility of internal lighter states becomes problematic and still needs a proper study. A non-perturbative, model-independent lattice evaluation of these processes would allow further progress in the theoretical predictions of SM hadronic quantities and in the search of New Physics., Comment: 11 pages, 16 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{SD}$ and $a_{\mu}^W$, defined by the RBC/UKQCD Collaboration [1]. We employ a subset of the gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks, which are close to the physical point for the masses of all the dynamical flavors. The simulations are carried out at three values of the lattice spacing ranging from $\simeq 0.057$ to $\simeq 0.080$ fm with linear lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{SD}({\rm ETMC}) = 69.33\,(29) \cdot 10^{-10}$, which is consistent with the recent dispersive value $a_\mu^{SD}(e^+ e^-) = 68.4\,(5) \cdot 10^{-10}$ [2] within $\simeq 1.6 \, \sigma$. In the case of the intermediate window we get the value $a_\mu^W({\rm ETMC}) = 235.0\,(1.1) \cdot 10^{-10}$, which is consistent with the result $a_\mu^W({\rm BMW}) = 236.7\,(1.4) \cdot 10^{-10}$ [3] by the BMW collaboration as well as with the recent determination by the CLS/Mainz group of $a_\mu^W({\rm CLS}) = 237.30\,(1.46) \cdot 10^{-10}$ [4] at the $\sim 1.0 - 1.3 \, \sigma$ level. However, it is larger than the dispersive result $a_\mu^W(e^+ e^-) = 229.4\,(1.4) \cdot 10^{-10}$ [2] by $\simeq 3.1 \, \sigma$. The tension increases to $\simeq 4.2 \, \sigma$ if we average our ETMC result with the BMW and the CLS/Mainz ones. Our accurate lattice results in the short and intermediate windows hint at possible deviations of the $e^+ e^-$ cross section data with respect to Standard Model (SM) predictions distributed somewhere in the low (and possibly intermediate) energy regions, but not in the high energy region., Comment: 82 pages, 33 figures, 14 tables

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavours of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavours tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD} = 69.27\,(34) \cdot 10^{-10}$, in agreement with the dispersive determination based on experimental $e^+ e^-$ data. For the intermediate window we get instead $a_\mu^{\rm W} = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with recent determinations by other lattice collaborations, but disagrees with the dispersive determination at the level of $3.6\,\sigma$., Comment: 10 pages, 4 figures, 1 table, LATTICE 2022

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavours of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavours tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD} = 69.27\,(34) \cdot 10^{-10}$, in agreement with the dispersive determination based on experimental $e^+ e^-$ data. For the intermediate window we get instead $a_\mu^{\rm W} = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with recent determinations by other lattice collaborations, but disagrees with the dispersive determination at the level of $3.6\,\sigma$., Comment: 10 pages, 4 figures, 1 table, LATTICE 2022

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavours of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavours tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD} = 69.27\,(34) \cdot 10^{-10}$, in agreement with the dispersive determination based on experimental $e^+ e^-$ data. For the intermediate window we get instead $a_\mu^{\rm W} = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with recent determinations by other lattice collaborations, but disagrees with the dispersive determination at the level of $3.6\,\sigma$., Comment: 10 pages, 4 figures, 1 table, LATTICE 2022

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$, defined by the RBC/UKQCD Collaboration [1]. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavors tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD}({\rm ETMC}) = 69.27\,(34) \cdot 10^{-10}$, which is consistent with the recent dispersive value of $a_\mu^{\rm SD}(e^+ e^-) = 68.4\,(5) \cdot 10^{-10}$ [2]. In the case of the intermediate window we get the value $a_\mu^{\rm W}({\rm ETMC}) = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with the result $a_\mu^{\rm W}({\rm BMW}) = 236.7\,(1.4) \cdot 10^{-10}$ [3] by the BMW collaboration as well as with the recent determination by the CLS/Mainz group of $a_\mu^{\rm W}({\rm CLS}) = 237.30\,(1.46) \cdot 10^{-10}$ [4]. However, it is larger than the dispersive result of $a_\mu^{\rm W}(e^+ e^-) = 229.4\,(1.4) \cdot 10^{-10}$ [2] by approximately $3.6$ standard deviations. The tension increases to approximately $4.5$ standard deviations if we average our ETMC result with those by BMW and CLS/Mainz. Our accurate lattice results in the short and intermediate windows point to a possible deviation of the $e^+ e^-$ cross section data with respect to Standard Model predictions in the low and intermediate energy regions, but not in the high energy region., Comment: 78 pages, 22 figures, 14 tables. Analysis improved with more data and fits, presentation reorganized, more material in appendices, fixed typos

We present a preliminary non-perturbative lattice calculation of the form factors entering the processes $K\to \ell\,\nu_\ell\,\ell'^+\,\ell'^-$ and of the corresponding branching ratios. These form factors describe the interaction between the mediating virtual photon and the internal hadronic structure of the meson. By separating them from the point-like contribution to the matrix element we are able to isolate and reconstruct the structure-dependent contribution to the decay width. Our numerical analysis employs only one gauge ensemble and so it is affected by systematic uncertainties due to the missing continuum and physical point extrapolation. Despite this, we already find a reasonable agreement with the experimental data and with the next to leading order Chiral Perturbation Theory predictions. The method is general and can be applied to any pseudoscalar meson, though for heavier mesons the possibility of internal lighter states becomes problematic and still needs a proper study. A non-perturbative, model-independent lattice evaluation of these processes would allow further progress in the theoretical predictions of SM hadronic quantities and in the search of New Physics., Comment: 11 pages, 16 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present a lattice determination of the leading-order hadronic vacuum polarization (HVP) contribution to the muon anomalous magnetic moment, $a_{\mu}^{\rm HVP}$, in the so-called short and intermediate time-distance windows, $a_{\mu}^{\rm SD}$ and $a_{\mu}^{\rm W}$, defined by the RBC/UKQCD Collaboration [1]. We employ gauge ensembles produced by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks with masses of all the dynamical quark flavors tuned close to their physical values. The simulations are carried out at three values of the lattice spacing equal to $\simeq 0.057, 0.068$ and $0.080$ fm with spatial lattice sizes up to $L \simeq 7.6$~fm. For the short distance window we obtain $a_\mu^{\rm SD}({\rm ETMC}) = 69.27\,(34) \cdot 10^{-10}$, which is consistent with the recent dispersive value of $a_\mu^{\rm SD}(e^+ e^-) = 68.4\,(5) \cdot 10^{-10}$ [2]. In the case of the intermediate window we get the value $a_\mu^{\rm W}({\rm ETMC}) = 236.3\,(1.3) \cdot 10^{-10}$, which is consistent with the result $a_\mu^{\rm W}({\rm BMW}) = 236.7\,(1.4) \cdot 10^{-10}$ [3] by the BMW collaboration as well as with the recent determination by the CLS/Mainz group of $a_\mu^{\rm W}({\rm CLS}) = 237.30\,(1.46) \cdot 10^{-10}$ [4]. However, it is larger than the dispersive result of $a_\mu^{\rm W}(e^+ e^-) = 229.4\,(1.4) \cdot 10^{-10}$ [2] by approximately $3.6$ standard deviations. The tension increases to approximately $4.5$ standard deviations if we average our ETMC result with those by BMW and CLS/Mainz. Our accurate lattice results in the short and intermediate windows point to a possible deviation of the $e^+ e^-$ cross section data with respect to Standard Model predictions in the low and intermediate energy regions, but not in the high energy region., Comment: 78 pages, 22 figures, 14 tables. Analysis improved with more data and fits, presentation reorganized, more material in appendices, fixed typos

We study the radiative leptonic decays $P\to\ell\nu_\ell\,\ell^{\prime\,+}\ell^{\prime\,-}$, where $P$ is a pseudoscalar meson and $\ell$ and $\ell^\prime$ are charged leptons. In such decays the emitted photon is off-shell and, in addition to the "point-like" contribution in which the virtual photon is emitted either from the lepton or the meson treated as a point-like particle, four structure-dependent (SD) form factors contribute to the amplitude. We present a strategy for the extraction of the SD form factors and implement it in an exploratory lattice computation of the decay rates for the four channels of kaon decays ($\ell,\ell^\prime=e,\mu$). It is the SD form factors which describe the interaction between the virtual photon and the internal hadronic structure of the decaying meson, and in our procedure we separate the SD and point-like contributions to the amplitudes. We demonstrate that the form factors can be extracted with good precision and, in spite of the unphysical quark masses used in our simulation ($m_\pi\simeq 320\,$MeV and $m_K\simeq 530\,$MeV), the results for the decay rates are in reasonable semiquantitative agreement with experimental data (for the channels where these exist). Following this preparatory work, the emphasis of our future work will be on obtaining results at physical quark masses and on the control of the systematic uncertainties associated with discretisation and finite-volume errors., Comment: 30 pages, 47 figures

We present a lattice calculation of the charged/neutral pion mass difference $M_{\pi^{+}}-M_{\pi^{0}}$ at order $\mathcal{O}(\alpha_{em})$ using the gauge configurations produced by the Extended Twisted Mass Collaboration with $N_{f}=2+1+1$ dynamical quark flavours at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089~{\rm fm}$) and pion masses in the range $M_{\pi} \simeq 250-500~{\rm MeV}$. We employ the RM123 method and expand the path-integral around the isospin symmetric point at leading order in the electromagnetic coupling $\alpha_{em}$. Making use of the recently proposed RTM scheme, we evaluate the full $\mathcal{O}(\alpha_{em})$ contribution, with the inclusion of the disconnected diagram. At the physical point, after performing the continuum and infinite volume extrapolation, we obtain the value $M_{\pi^{+}}-M_{\pi^{0}}= 4.622~(95)~{\rm MeV}$ which is in good agreement with the experimental result $[ M_{\pi^{+}} - M_{\pi^{0}} ]^{exp.} = 4.5936(5)~{\rm MeV}$., Comment: 13 pages, 5 figures, 1 table

We present our preliminary results concerning the charged/neutral pion mass difference $M_{\pi^{+}} - M_{\pi^{0}}$ at order $\mathcal{O}(\alpha_{em})$ in the QED interactions, and for $M_{\pi^{+}} - M_{\pi^{0}}$ at order $\mathcal{O}\left( (m_{d}-m_{u})^{2}\right)$ in the strong isospin-breaking term. The latter contribution provides a determination of the $\rm{SU}(2)$ chiral perturbation theory low-energy constant $\ell_{7}$, whose present estimate is affected by a rather large uncertainty. The disconnected contributions appearing in the diagrammatic expansion of $M_{\pi^{+}} - M_{\pi^{0}}$, being very noisy, are notoriously difficult to evaluate and have been neglected in our previous calculations. By making use of twisted mass Lattice QCD simulations and adopting the RM123 method, we will show that taking profit from our recently proposed rotated twisted-mass (RTM) scheme, tailored to improve the signal on these kinds of observables, it is possible to evaluate the disconnected diagrams with good precision. For the QED induced pion mass difference, we obtain, after performing the extrapolation towards the continuum and thermodynamic limit and at the physical point, the preliminary value $M_{\pi^{+}}-M_{\pi^{0}} = 4.622~(95)~{\rm MeV}$, that is in good agreement with the experimental result. For the determination of the low-energy constant $\ell_{7}$, our result $\ell_{7} = 2.5~(1.4)\times 10^{-3}$, which is limited so far to a single lattice spacing, is in agreement and improves phenomenological estimates., Comment: 12 pages, 6 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We evaluate by means of lattice QCD calculations the low-energy constant $\ell_{7}$ which parametrizes strong isospin effects at NLO in $\rm{SU}(2)$ chiral perturbation theory. Among all low-energy constants at NLO, $\ell_{7}$ is the one known less precisely, and its uncertainty is currently larger than $50\%$. Our strategy is based on the RM123 approach in which the lattice path-integral is expanded in powers of the isospin breaking parameter $\Delta m= (m_{d}-m_{u})/2$. In order to evaluate the relevant lattice correlators we make use of the recently proposed rotated twisted-mass (RTM) scheme. Within the RM123 approach, it is possible to cleanly extract the value of $\ell_{7}$ from either the pion mass splitting $M_{\pi^{+}}-M_{\pi^{0}}$ induced by strong isospin breaking at order $\mathcal{O}\left((\Delta m)^{2}\right)$ (mass method), or from the coupling of the neutral pion $\pi^{0}$ to the isoscalar operator $\left(\bar{u}\gamma_{5}u + \bar{d}\gamma_{5} d\right)/\sqrt{2}$ at order $\mathcal{O}(\Delta m)$ (matrix element method). In this pilot study we limit the analysis to a single ensemble generated by the Extended Twisted Mass Collaboration (ETMC) with $N_{f}=2+1+1$ dynamical quark flavours, which corresponds to a lattice spacing $a\simeq 0.095~{\rm fm}$ and to a pion mass $M_{\pi}\simeq 260~{\rm MeV}$. We find that the matrix element method outperforms the mass method in terms of resulting statistical accuracy. Our determination, $\ell_{7} = 2.5(1.4)\times 10^{-3}$, is in agreement and improves previous calculations., Comment: 17 pages, 3 figures

We present results for the light, strange and charm quark masses using $N_f=2+1+1$ twisted mass fermion ensembles at three values of the lattice spacing, including two ensembles simulated with the physical value of the pion mass. The analysis is done both in the meson and baryon sectors. The difference in the mean values found in the two sectors is included as part of the systematic error. The presentation is based on the work of Ref. [1], where more details can be found., Comment: Talk presented at the 38th International Symposium on Lattice Field Theory, Lattice 2021, Zoom/Gather@Massachusetts Institute of Technology, 26-30 Jul. 2021; 8 pages, 5 figures

We propose a scheme of lattice twisted-mass fermion regularization which is particularly convenient for application to isospin breaking (IB) QCD and QED calculations, based in particular on the so called RM123 approach, in which the IB terms of the action are treated as a perturbation. The main, practical advantage of this scheme is that it allows the calculation of IB effects on some mesonic observables, like e.g. the pi+ - pi0 mass splitting, using lattice correlation functions in which the quark and antiquark fields in the meson are regularized with opposite values of the Wilson parameter r. These correlation functions are found to be affected by much smaller statistical fluctuations, with respect to the analogous functions in which quark and antiquark fields are regularized with the same value of r. Two numerical application of this scheme, that we call "rotated twisted-mass", within pure QCD and QCD+QED respectively, are also provided for illustration., Comment: 17 pages, 2 figures. v2: minor corrections, references added. Matches the published version

We present a calculation of the up, down, strange and charm quark masses performed within the lattice QCD framework. We use the twisted mass fermion action and carry out simulations that include in the sea two light mass-degenerate quarks, as well as the strange and charm quarks. In the analysis we use gauge ensembles simulated at three values of the lattice spacing and with light quarks that correspond to pion masses in the range from 350 MeV to the physical value, while the strange and charm quark masses are tuned approximately to their physical values. We use several quantities to set the scale in order to check for finite lattice spacing effects and in the continuum limit we get compatible results. The quark mass renormalization is carried out non-perturbatively using the RI'-MOM method converted into the $\overline{\rm MS}$ scheme. For the determination of the quark masses we use physical observables from both the meson and the baryon sectors, obtaining $m_{ud} = 3.636(66)(^{+60}_{-57})$~MeV and $m_s = 98.7(2.4)(^{+4.0}_{-3.2})$~MeV in the $\overline{\rm MS}(2\,{\rm GeV})$ scheme and $m_c = 1036(17)(^{+15}_{-8})$~MeV in the $\overline{\rm MS}(3\,{\rm GeV})$ scheme, where the first errors are statistical and the second ones are combinations of systematic errors. For the quark mass ratios we get $m_s / m_{ud} = 27.17(32)(^{+56}_{-38})$ and $m_c / m_s = 11.48(12)(^{+25}_{-19})$., Comment: 37 pages, 23 figures, 24 tables. One reference added

We present the calculation of the short-distance power corrections to the CP-violation parameter $\varepsilon_K$ coming from dimension-8 operators in the $\Delta S=2$ effective Hamiltonian. A first estimate of this contribution, obtained for large-$N_c$ and in the chiral limit, was provided in ref. [arXiv:hep-ph/0406094]. Here we evaluate and include the $\mathcal{O}(m_K^2/m_c^2)$ and $\mathcal{O}(\Lambda_{QCD}^2/m_c^2)$ contributions that, a priori, could induce $\mathcal{O}(1)$ corrections to previous estimates, as $m_K$ is numerically of order $\Lambda_{QCD}$. Our computation shows that there are several dimension-8 operators besides the one considered before. Their effect on $\varepsilon_K$, however, accidentally cancels out to a large extent, leaving the final correction at the level of 1%., Comment: 23 pages, 2 figures. v2: added references and CKM parameter values. v3: matches published version

We present our preliminary results concerning the charged/neutral pion mass difference $M_{\pi^{+}} - M_{\pi^{0}}$ at order $\mathcal{O}(\alpha_{em})$ in the QED interactions, and for $M_{\pi^{+}} - M_{\pi^{0}}$ at order $\mathcal{O}\left( (m_{d}-m_{u})^{2}\right)$ in the strong isospin-breaking term. The latter contribution provides a determination of the $\rm{SU}(2)$ chiral perturbation theory low-energy constant $\ell_{7}$, whose present estimate is affected by a rather large uncertainty. The disconnected contributions appearing in the diagrammatic expansion of $M_{\pi^{+}} - M_{\pi^{0}}$, being very noisy, are notoriously difficult to evaluate and have been neglected in our previous calculations. By making use of twisted mass Lattice QCD simulations and adopting the RM123 method, we will show that taking profit from our recently proposed rotated twisted-mass (RTM) scheme, tailored to improve the signal on these kinds of observables, it is possible to evaluate the disconnected diagrams with good precision. For the QED induced pion mass difference, we obtain, after performing the extrapolation towards the continuum and thermodynamic limit and at the physical point, the preliminary value $M_{\pi^{+}}-M_{\pi^{0}} = 4.622~(95)~{\rm MeV}$, that is in good agreement with the experimental result. For the determination of the low-energy constant $\ell_{7}$, our result $\ell_{7} = 2.5~(1.4)\times 10^{-3}$, which is limited so far to a single lattice spacing, is in agreement and improves phenomenological estimates., Comment: 12 pages, 6 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We present results for the light, strange and charm quark masses using $N_f=2+1+1$ twisted mass fermion ensembles at three values of the lattice spacing, including two ensembles simulated with the physical value of the pion mass. The analysis is done both in the meson and baryon sectors. The difference in the mean values found in the two sectors is included as part of the systematic error. The presentation is based on the work of Ref. [1], where more details can be found., Comment: Talk presented at the 38th International Symposium on Lattice Field Theory, Lattice 2021, Zoom/Gather@Massachusetts Institute of Technology, 26-30 Jul. 2021; 8 pages, 5 figures

We evaluate by means of lattice QCD calculations the low-energy constant $\ell_{7}$ which parametrizes strong isospin effects at NLO in $\rm{SU}(2)$ chiral perturbation theory. Among all low-energy constants at NLO, $\ell_{7}$ is the one known less precisely, and its uncertainty is currently larger than $50\%$. Our strategy is based on the RM123 approach in which the lattice path-integral is expanded in powers of the isospin breaking parameter $\Delta m= (m_{d}-m_{u})/2$. In order to evaluate the relevant lattice correlators we make use of the recently proposed rotated twisted-mass (RTM) scheme. Within the RM123 approach, it is possible to cleanly extract the value of $\ell_{7}$ from either the pion mass splitting $M_{\pi^{+}}-M_{\pi^{0}}$ induced by strong isospin breaking at order $\mathcal{O}\left((\Delta m)^{2}\right)$ (mass method), or from the coupling of the neutral pion $\pi^{0}$ to the isoscalar operator $\left(\bar{u}\gamma_{5}u + \bar{d}\gamma_{5} d\right)/\sqrt{2}$ at order $\mathcal{O}(\Delta m)$ (matrix element method). In this pilot study we limit the analysis to a single ensemble generated by the Extended Twisted Mass Collaboration (ETMC) with $N_{f}=2+1+1$ dynamical quark flavours, which corresponds to a lattice spacing $a\simeq 0.095~{\rm fm}$ and to a pion mass $M_{\pi}\simeq 260~{\rm MeV}$. We find that the matrix element method outperforms the mass method in terms of resulting statistical accuracy. Our determination, $\ell_{7} = 2.5(1.4)\times 10^{-3}$, is in agreement and improves previous calculations., Comment: 17 pages, 3 figures

We present a calculation of the up, down, strange and charm quark masses performed within the lattice QCD framework. We use the twisted mass fermion action and carry out simulations that include in the sea two light mass-degenerate quarks, as well as the strange and charm quarks. In the analysis we use gauge ensembles simulated at three values of the lattice spacing and with light quarks that correspond to pion masses in the range from 350 MeV to the physical value, while the strange and charm quark masses are tuned approximately to their physical values. We use several quantities to set the scale in order to check for finite lattice spacing effects and in the continuum limit we get compatible results. The quark mass renormalization is carried out non-perturbatively using the RI'-MOM method converted into the $\overline{\rm MS}$ scheme. For the determination of the quark masses we use physical observables from both the meson and the baryon sectors, obtaining $m_{ud} = 3.636(66)(^{+60}_{-57})$~MeV and $m_s = 98.7(2.4)(^{+4.0}_{-3.2})$~MeV in the $\overline{\rm MS}(2\,{\rm GeV})$ scheme and $m_c = 1036(17)(^{+15}_{-8})$~MeV in the $\overline{\rm MS}(3\,{\rm GeV})$ scheme, where the first errors are statistical and the second ones are combinations of systematic errors. For the quark mass ratios we get $m_s / m_{ud} = 27.17(32)(^{+56}_{-38})$ and $m_c / m_s = 11.48(12)(^{+25}_{-19})$., Comment: 37 pages, 23 figures, 24 tables. One reference added

We propose a scheme of lattice twisted-mass fermion regularization which is particularly convenient for application to isospin breaking (IB) QCD and QED calculations, based in particular on the so called RM123 approach, in which the IB terms of the action are treated as a perturbation. The main, practical advantage of this scheme is that it allows the calculation of IB effects on some mesonic observables, like e.g. the pi+ - pi0 mass splitting, using lattice correlation functions in which the quark and antiquark fields in the meson are regularized with opposite values of the Wilson parameter r. These correlation functions are found to be affected by much smaller statistical fluctuations, with respect to the analogous functions in which quark and antiquark fields are regularized with the same value of r. Two numerical application of this scheme, that we call "rotated twisted-mass", within pure QCD and QCD+QED respectively, are also provided for illustration., Comment: 17 pages, 2 figures. v2: minor corrections, references added. Matches the published version

We propose a scheme of lattice twisted-mass fermion regularization which is particularly convenient for application to isospin breaking (IB) QCD and QED calculations, based in particular on the so called RM123 approach, in which the IB terms of the action are treated as a perturbation. The main, practical advantage of this scheme is that it allows the calculation of IB effects on some mesonic observables, like e.g. the pi+ - pi0 mass splitting, using lattice correlation functions in which the quark and antiquark fields in the meson are regularized with opposite values of the Wilson parameter r. These correlation functions are found to be affected by much smaller statistical fluctuations, with respect to the analogous functions in which quark and antiquark fields are regularized with the same value of r. Two numerical application of this scheme, that we call "rotated twisted-mass", within pure QCD and QCD+QED respectively, are also provided for illustration., Comment: 17 pages, 2 figures. v2: minor corrections, references added. Matches the published version

We present results for the light, strange and charm quark masses using $N_f=2+1+1$ twisted mass fermion ensembles at three values of the lattice spacing, including two ensembles simulated with the physical value of the pion mass. The analysis is done both in the meson and baryon sectors. The difference in the mean values found in the two sectors is included as part of the systematic error. The presentation is based on the work of Ref. [1], where more details can be found., Comment: Talk presented at the 38th International Symposium on Lattice Field Theory, Lattice 2021, Zoom/Gather@Massachusetts Institute of Technology, 26-30 Jul. 2021; 8 pages, 5 figures

We present our preliminary results concerning the charged/neutral pion mass difference $M_{\pi^{+}} - M_{\pi^{0}}$ at order $\mathcal{O}(\alpha_{em})$ in the QED interactions, and for $M_{\pi^{+}} - M_{\pi^{0}}$ at order $\mathcal{O}\left( (m_{d}-m_{u})^{2}\right)$ in the strong isospin-breaking term. The latter contribution provides a determination of the $\rm{SU}(2)$ chiral perturbation theory low-energy constant $\ell_{7}$, whose present estimate is affected by a rather large uncertainty. The disconnected contributions appearing in the diagrammatic expansion of $M_{\pi^{+}} - M_{\pi^{0}}$, being very noisy, are notoriously difficult to evaluate and have been neglected in our previous calculations. By making use of twisted mass Lattice QCD simulations and adopting the RM123 method, we will show that taking profit from our recently proposed rotated twisted-mass (RTM) scheme, tailored to improve the signal on these kinds of observables, it is possible to evaluate the disconnected diagrams with good precision. For the QED induced pion mass difference, we obtain, after performing the extrapolation towards the continuum and thermodynamic limit and at the physical point, the preliminary value $M_{\pi^{+}}-M_{\pi^{0}} = 4.622~(95)~{\rm MeV}$, that is in good agreement with the experimental result. For the determination of the low-energy constant $\ell_{7}$, our result $\ell_{7} = 2.5~(1.4)\times 10^{-3}$, which is limited so far to a single lattice spacing, is in agreement and improves phenomenological estimates., Comment: 12 pages, 6 figures. Talk presented at the 38th International Symposium on Lattice Field Theory (Lattice 2021), 26-30 July, 2021, Zoom/Gather@Massachusetts Institute of Technology

We evaluate by means of lattice QCD calculations the low-energy constant $\ell_{7}$ which parametrizes strong isospin effects at NLO in $\rm{SU}(2)$ chiral perturbation theory. Among all low-energy constants at NLO, $\ell_{7}$ is the one known less precisely, and its uncertainty is currently larger than $50\%$. Our strategy is based on the RM123 approach in which the lattice path-integral is expanded in powers of the isospin breaking parameter $\Delta m= (m_{d}-m_{u})/2$. In order to evaluate the relevant lattice correlators we make use of the recently proposed rotated twisted-mass (RTM) scheme. Within the RM123 approach, it is possible to cleanly extract the value of $\ell_{7}$ from either the pion mass splitting $M_{\pi^{+}}-M_{\pi^{0}}$ induced by strong isospin breaking at order $\mathcal{O}\left((\Delta m)^{2}\right)$ (mass method), or from the coupling of the neutral pion $\pi^{0}$ to the isoscalar operator $\left(\bar{u}\gamma_{5}u + \bar{d}\gamma_{5} d\right)/\sqrt{2}$ at order $\mathcal{O}(\Delta m)$ (matrix element method). In this pilot study we limit the analysis to a single ensemble generated by the Extended Twisted Mass Collaboration (ETMC) with $N_{f}=2+1+1$ dynamical quark flavours, which corresponds to a lattice spacing $a\simeq 0.095~{\rm fm}$ and to a pion mass $M_{\pi}\simeq 260~{\rm MeV}$. We find that the matrix element method outperforms the mass method in terms of resulting statistical accuracy. Our determination, $\ell_{7} = 2.5(1.4)\times 10^{-3}$, is in agreement and improves previous calculations., Comment: 17 pages, 3 figures

We present a calculation of the up, down, strange and charm quark masses performed within the lattice QCD framework. We use the twisted mass fermion action and carry out simulations that include in the sea two light mass-degenerate quarks, as well as the strange and charm quarks. In the analysis we use gauge ensembles simulated at three values of the lattice spacing and with light quarks that correspond to pion masses in the range from 350 MeV to the physical value, while the strange and charm quark masses are tuned approximately to their physical values. We use several quantities to set the scale in order to check for finite lattice spacing effects and in the continuum limit we get compatible results. The quark mass renormalization is carried out non-perturbatively using the RI'-MOM method converted into the $\overline{\rm MS}$ scheme. For the determination of the quark masses we use physical observables from both the meson and the baryon sectors, obtaining $m_{ud} = 3.636(66)(^{+60}_{-57})$~MeV and $m_s = 98.7(2.4)(^{+4.0}_{-3.2})$~MeV in the $\overline{\rm MS}(2\,{\rm GeV})$ scheme and $m_c = 1036(17)(^{+15}_{-8})$~MeV in the $\overline{\rm MS}(3\,{\rm GeV})$ scheme, where the first errors are statistical and the second ones are combinations of systematic errors. For the quark mass ratios we get $m_s / m_{ud} = 27.17(32)(^{+56}_{-38})$ and $m_c / m_s = 11.48(12)(^{+25}_{-19})$., Comment: 37 pages, 23 figures, 24 tables. One reference added

In the original version of this manuscript, an error was introduced on pp352. '2.7nb:1.6nb' has been corrected to '2.4nb:1.3nb' in the current online and printed version. doi:10.1093/ptep/ptz106.

We present a comparison of existing experimental data for the radiative leptonic decays $P\to\ell\nu_\ell\gamma$, where $P=K$ or $\pi$ and $\ell=e$ or $\mu$, from the KLOE, PIBETA, E787, ISTRA+ and OKA collaborations with theoretical predictions based on the recent non-perturbative determinations of the structure-dependent vector and axial-vector form factors, $F_V$ and $F_A$ respectively. These were obtained using lattice QCD+QED simulations at order $O(\alpha_{\mathrm{em}})$ in the electromagnetic coupling. We find good agreement with the KLOE data on $K\to e\nu_e\gamma$ decays from which the form factor $F^+=F_V+F_A$ can be determined. For $K\to\mu\nu_\mu\gamma$ decays we observe differences of up to 3\,-\,4 standard deviations at large photon energies between the theoretical predictions and the data from the E787, ISTRA+ and OKA experiments and similar discrepancies in some kinematical regions with the PIBETA experiment on radiative pion decays. A global study of all the kaon-decay data within the Standard Model results in a poor fit, largely because at large photon energies the KLOE and E787 data cannot be reproduced simultaneously in terms of the same form factor $F^+$. The discrepancy between the theoretical and experimental values of the form factor $F^-=F_V-F_A$ is even more pronounced. These observations motivate future improvements of both the theoretical and experimental determinations of the structure-dependent form factors $F^+$ and $F^-$, as well as further theoretical investigations of models of "new physics" which might for example, include possible flavor changing interactions beyond $V - A$ and/or non-universal corrections to the lepton couplings., Comment: 31 pages, 8 figures, 11 tables. Version to appear in Phys. Rev. D

We present a non-perturbative lattice calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar \nu_\ell \gamma$, where $P$ is a pseudoscalar meson and $\ell$ is a charged lepton. Together with the non-perturbative determination of the corrections to the processes $P\to \ell \bar \nu_\ell$ due to the exchange of a virtual photon, this allows accurate predictions at $O(\alpha_{em})$ to be made for leptonic decay rates for pseudoscalar mesons ranging from the pion to the $D_s$ meson. We are able to separate unambiguously and non-pertubatively the point-like contribution, from the structure-dependent, infrared-safe, terms in the amplitude. The fully non-perturbative $O(a)$ improved calculation of the inclusive leptonic decay rates will lead to the determination of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements also at $O(\alpha_{em})$. Prospects for a precise evaluation of leptonic decay rates with emission of a hard photon are also very interesting, especially for the decays of heavy $D$ and $B$ mesons for which currently only model-dependent predictions are available to compare with existing experimental data., Comment: 46 pages, 12 figures

We present a comparison of existing experimental data for the radiative leptonic decays $P\to\ell\nu_\ell\gamma$, where $P=K$ or $\pi$ and $\ell=e$ or $\mu$, from the KLOE, PIBETA, E787, ISTRA+ and OKA collaborations with theoretical predictions based on the recent non-perturbative determinations of the structure-dependent vector and axial-vector form factors, $F_V$ and $F_A$ respectively. These were obtained using lattice QCD+QED simulations at order $O(\alpha_{\mathrm{em}})$ in the electromagnetic coupling. We find good agreement with the KLOE data on $K\to e\nu_e\gamma$ decays from which the form factor $F^+=F_V+F_A$ can be determined. For $K\to\mu\nu_\mu\gamma$ decays we observe differences of up to 3\,-\,4 standard deviations at large photon energies between the theoretical predictions and the data from the E787, ISTRA+ and OKA experiments and similar discrepancies in some kinematical regions with the PIBETA experiment on radiative pion decays. A global study of all the kaon-decay data within the Standard Model results in a poor fit, largely because at large photon energies the KLOE and E787 data cannot be reproduced simultaneously in terms of the same form factor $F^+$. The discrepancy between the theoretical and experimental values of the form factor $F^-=F_V-F_A$ is even more pronounced. These observations motivate future improvements of both the theoretical and experimental determinations of the structure-dependent form factors $F^+$ and $F^-$, as well as further theoretical investigations of models of "new physics" which might for example, include possible flavor changing interactions beyond $V - A$ and/or non-universal corrections to the lepton couplings., Comment: 31 pages, 8 figures, 11 tables. Version to appear in Phys. Rev. D

We present a non-perturbative lattice calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar \nu_\ell \gamma$, where $P$ is a pseudoscalar meson and $\ell$ is a charged lepton. Together with the non-perturbative determination of the corrections to the processes $P\to \ell \bar \nu_\ell$ due to the exchange of a virtual photon, this allows accurate predictions at $O(\alpha_{em})$ to be made for leptonic decay rates for pseudoscalar mesons ranging from the pion to the $D_s$ meson. We are able to separate unambiguously and non-pertubatively the point-like contribution, from the structure-dependent, infrared-safe, terms in the amplitude. The fully non-perturbative $O(a)$ improved calculation of the inclusive leptonic decay rates will lead to the determination of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements also at $O(\alpha_{em})$. Prospects for a precise evaluation of leptonic decay rates with emission of a hard photon are also very interesting, especially for the decays of heavy $D$ and $B$ mesons for which currently only model-dependent predictions are available to compare with existing experimental data., Comment: 46 pages, 12 figures

We present a comparison of existing experimental data for the radiative leptonic decays $P\to\ell\nu_\ell\gamma$, where $P=K$ or $\pi$ and $\ell=e$ or $\mu$, from the KLOE, PIBETA, E787, ISTRA+ and OKA collaborations with theoretical predictions based on the recent non-perturbative determinations of the structure-dependent vector and axial-vector form factors, $F_V$ and $F_A$ respectively. These were obtained using lattice QCD+QED simulations at order $O(\alpha_{\mathrm{em}})$ in the electromagnetic coupling. We find good agreement with the KLOE data on $K\to e\nu_e\gamma$ decays from which the form factor $F^+=F_V+F_A$ can be determined. For $K\to\mu\nu_\mu\gamma$ decays we observe differences of up to 3\,-\,4 standard deviations at large photon energies between the theoretical predictions and the data from the E787, ISTRA+ and OKA experiments and similar discrepancies in some kinematical regions with the PIBETA experiment on radiative pion decays. A global study of all the kaon-decay data within the Standard Model results in a poor fit, largely because at large photon energies the KLOE and E787 data cannot be reproduced simultaneously in terms of the same form factor $F^+$. The discrepancy between the theoretical and experimental values of the form factor $F^-=F_V-F_A$ is even more pronounced. These observations motivate future improvements of both the theoretical and experimental determinations of the structure-dependent form factors $F^+$ and $F^-$, as well as further theoretical investigations of models of "new physics" which might for example, include possible flavor changing interactions beyond $V - A$ and/or non-universal corrections to the lepton couplings., Comment: 31 pages, 8 figures, 11 tables. Version to appear in Phys. Rev. D

We present a non-perturbative lattice calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar \nu_\ell \gamma$, where $P$ is a pseudoscalar meson and $\ell$ is a charged lepton. Together with the non-perturbative determination of the corrections to the processes $P\to \ell \bar \nu_\ell$ due to the exchange of a virtual photon, this allows accurate predictions at $O(\alpha_{em})$ to be made for leptonic decay rates for pseudoscalar mesons ranging from the pion to the $D_s$ meson. We are able to separate unambiguously and non-pertubatively the point-like contribution, from the structure-dependent, infrared-safe, terms in the amplitude. The fully non-perturbative $O(a)$ improved calculation of the inclusive leptonic decay rates will lead to the determination of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements also at $O(\alpha_{em})$. Prospects for a precise evaluation of leptonic decay rates with emission of a hard photon are also very interesting, especially for the decays of heavy $D$ and $B$ mesons for which currently only model-dependent predictions are available to compare with existing experimental data., Comment: 46 pages, 12 figures

The Belle II detector will provide a major step forward in precision heavy flavor physics, quarkonium and exotic states, searches for dark sectors, and many other areas. The sensitivity to a large number of key observables can be improved by about an order of magnitude compared to the current measurements, and up to two orders in very clean search measurements. This increase in statistical precision arises not only due to the increased luminosity, but also from improved detector efficiency and precision for many channels. Many of the most interesting observables tend to have very small theoretical uncertainties that will therefore not limit the physics reach. This book has presented many new ideas for measurements, both to elucidate the nature of current anomalies seen in flavor, and to search for new phenomena in a plethora of observables that will become accessible with the Belle II dataset. The simulation used for the studiesinthis book was state ofthe artat the time, though weare learning a lot more about the experiment during the commissioning period. The detector is in operation, and working spectacularly well.

The Belle II detector will provide a major step forward in precision heavy flavor physics, quarkonium and exotic states, searches for dark sectors, and many other areas. The sensitivity to a large number of key observables can be improved by about an order of magnitude compared to the current measurements, and up to two orders in very clean search measurements. This increase in statistical precision arises not only due to the increased luminosity, but also from improved detector efficiency and precision for many channels. Many of the most interesting observables tend to have very small theoretical uncertainties that will therefore not limit the physics reach. This book has presented many new ideas for measurements, both to elucidate the nature of current anomalies seen in flavor, and to search for new phenomena in a plethora of observables that will become accessible with the Belle II dataset. The simulation used for the studiesinthis book was state ofthe artat the time, though weare learning a lot more about the experiment during the commissioning period. The detector is in operation, and working spectacularly well.

We present a non-perturbative calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar \nu_\ell \gamma$, where $P$ is a pseudoscalar meson and $\ell$ is a charged lepton. Together with the non-perturbative determination of the virtual photon corrections to the processes $P\to \ell \bar \nu_\ell$, this will allow accurate predictions to be made at $O(\alpha_{em})$ for leptonic decay rates for pseudoscalar mesons ranging from the pion to the $B$ meson. We are able to separate unambiguously the point-like contribution, the square of which leads to the infrared divergence in the decay rate, from the structure dependent, infrared-safe, terms in the amplitude. The fully non-perturbative, $O(a)$ improved calculation of the inclusive leptonic decay rates will lead to significantly improved precision in the determination of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Precise predictions for the emission of a hard photon are also very interesting, especially for the decays of heavy $D$ and $B$ mesons for which currently only model-dependent predictions are available to compare with existing experimental data.

We discuss the theoretical framework required for the computation of radiative corrections to semileptonic decay rates in lattice simulations, and in particular to those for $K_{\ell3}$ decays. This is an extension of the framework we have developed and successfully implemented for leptonic decays. New issues which arise for semileptonic decays, include the presence of unphysical terms which grow exponentially with the time separation between the insertion of the weak Hamiltonian and the sink for the final-state meson-lepton pair. Such terms must be identified and subtracted. We discuss the cancellation of infrared divergences and show that, with the QED$_\mathrm{\,L}$ treatment of the zero mode in the photon propagator, the $O(1/L)$ finite-volume corrections are "universal". These corrections however, depend not only on the semileptonic form factors $f^\pm(q^2)$ but also on their derivatives $df^\pm/dq^2$. (Here $q$ is the momentum transfer between the initial and final state mesons.) We explain the perturbative calculation which would need to be performed to subtract the $O(1/L)$ finite-volume effects., Comment: Presented at the 37th International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, China

The leading electromagnetic (e.m.) and strong isospin-breaking corrections to the $\pi^+ \to \mu^+ \nu[\gamma]$ and $K^+ \to \mu^+ \nu[\gamma]$ leptonic decay rates are evaluated for the first time on the lattice. The results are obtained using gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks. The relative leading-order e.m.~and strong isospin-breaking corrections to the decay rates are 1.53(19)\% for $\pi_{\mu 2}$ decays and 0.24(10)\% for $K_{\mu 2}$ decays. Using the experimental values of the $\pi_{\mu 2}$ and $K_{\mu 2}$ decay rates and updated lattice QCD results for the pion and kaon decay constants in isosymmetric QCD, we find that the Cabibbo-Kobayashi-Maskawa matrix element $ | V_{us}| = 0.22538(46)$, reducing by a factor of about $1.8$ the corresponding uncertainty in the Particle Data Group review. Our calculation of $|V_{us}|$ allows also an accurate determination of the first-row CKM unitarity relation $| V_{ud}|^2 + | V_{us}|^2 + | V_{ub}|^2 = 0.99988(46)$. Theoretical developments in this paper include a detailed discussion of how QCD can be defined in the full QCD+QED theory and an improved renormalisation procedure in which the bare lattice operators are renormalised non-perturbatively into the (modified) Regularization Independent Momentum subtraction scheme and subsequently matched perturbatively at $O(\alpha_{em}\alpha_s(M_W))$ into the W-regularisation scheme appropriate for these calculations., Comment: 63 pages, 10 figures and 2 tables. Version matches the published paper