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 an updated determination of the values of $\vert V_{cb} \vert$, $R(D^*)$ and $\vert V_{ub} \vert/\vert V_{cb} \vert$ based on the new data on semileptonic $B \to D^* \ell \nu_\ell$ decays by the Belle and Belle-II Collaborations and on the recent theoretical progress in the calculation of the form factors relevant for semileptonic $B \to D^* \ell \nu_\ell$ and $B_s \to K \ell \nu_\ell$ decays. In particular we present results derived by applying either the Dispersive Matrix (DM) method of Refs. [1-6] or the more standard Boyd-Grinstein-Lebed (BGL) [7] approach to the most recent values of the form factors determined in lattice QCD. Using all the available lattice results for the form factors from the DM method we get the theoretical value $R^{\rm th}(D^*) = 0.262 \pm 0.009$ and we extract from a bin-per-bin analysis of the experimental data the value $\vert V_{cb} \vert = (39.92 \pm 0.64) \cdot10^{-3}$. Our result for $R(D^*)$ is consistent with the latest experimental world average $R^{\rm exp}(D^*) = 0.284 \pm 0.012$ [8]} at the $\simeq 1.5\,\sigma$ level. Our value for $\vert V_{cb} \vert$ is compatible with the latest inclusive determinations $\vert V_{cb} \vert^{\rm incl} = (41.97 \pm 0.48) \cdot 10^{-3}$ [9] and $\vert V_{cb} \vert^{\rm incl} = (41.69\pm 0.63) \cdot 10^{-3}$ [10] within $\simeq 2.6$ and $\simeq 2.0$ standard deviations, respectively. From a reappraisal of the calculations of $\vert V_{ub} \vert / \vert V_{cb} \vert$, we also obtain $\vert V_{ub} \vert / \vert V_{cb} \vert = 0.087\pm 0.009$ in good agreement with the result $\vert V_{ub} \vert / \vert V_{cb} \vert = 0.0844\pm 0.0056$ from the latest FLAG review [11]., Comment: 28 pages, 9 figures, 6 tables. Recent update of the HPQCD results for the semileptonic $B \to D^*$ form factors included in the lattice inputs for the bin-per-bin analyses. Figs. 1-2,4,6-8 and Tables 1-4 updated. Main conclusions unchanged. Covariance matrices for various theoretical observables added as appendix for applications to phenomenological studies

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 the results of the application of the Dispersion Matrix approach to semileptonic heavy-to-heavy and heavy-to-light $B$-meson decays. This method allows to determine the hadronic form factors in a non-perturbative and model-independent way. Starting from the available lattice results at large values of the momentum transfer, we obtain the behaviour of the form factors in their whole kinematical range without introducing any parameterization of their momentum dependence. We will focus on the determination of the Cabibbo-Kobayashi-Maskawa matrix elements $\vert V_{cb} \vert$ and $\vert V_{ub} \vert$ through the analysis of $B \to D^{(*)} \ell \nu$, $B_s \to D_s^{(*)} \ell \nu$, $B \to \pi \ell \nu$ and $B_s \to K \ell \nu$ decays. New theoretical determinations of the Lepton Flavour Universality ratios relevant for these transitions will be also presented, by focusing in particular on the $R(D_{(s)}^{(*)})$ ratios., Comment: 9 pages, 3 figures. Proceedings for the 39th International Symposium on Lattice Field Theory (Lattice 2022). arXiv admin note: substantial text overlap with arXiv:2209.15413, arXiv:2205.13952; text overlap with arXiv:2211.07236

Agglomeration-based strategies are important both within adaptive refinement algorithms and to construct scalable multilevel algebraic solvers. In order to automatically perform agglomeration of polygonal grids, we propose the use of Machine Learning (ML) strategies, that can naturally exploit geometrical information about the mesh in order to preserve the grid quality, enhancing performance of numerical methods and reducing the overall computational cost. In particular, we employ the k-means clustering algorithm and Graph Neural Networks (GNNs) to partition the connectivity graph of a computational mesh. Moreover, GNNs have high online inference speed and the advantage to process naturally and simultaneously both the graph structure of mesh and the geometrical information, such as the areas of the elements or their barycentric coordinates. These techniques are compared with METIS, a standard algorithm for graph partitioning, which is meant to process only the graph information of the mesh. We demonstrate that performance in terms of quality metrics is enhanced for ML strategies. Such models also show a good degree of generalization when applied to more complex geometries, such as brain MRI scans, and the capability of preserving the quality of the grid. The effectiveness of these strategies is demonstrated also when applied to MultiGrid (MG) solvers in a Polygonal Discontinuous Galerkin (PolyDG) framework. In the considered experiments, GNNs show overall the best performance in terms of inference speed, accuracy and flexibility of the approach.

We present the results of the application of the Dispersion Matrix approach to exclusive semileptonic $B$-meson decays. This method allows to determine the hadronic form factors in a non-perturbative and completely model-independent way. Starting from lattice results available at large values of the momentum transfer, the behaviour of the form factors in their whole kinematical range is obtained without introducing any parameterization of their momentum dependence. We will focus on the determination of the Cabibbo-Kobayashi-Maskawa matrix elements $\vert V_{cb} \vert$ and $\vert V_{ub} \vert$ through the analysis of $B_{(s)} \to D_{(s)}^{(*)} \ell \nu$ and $B_{(s)} \to \pi(K) \ell \nu$ decays. New theoretical determinations of the Lepton Flavour Universality ratios relevant for these transitions will be also presented., Comment: 8 pages, 5 figures, 2 tables; contribution to QCD@Work - International Workshop on QCD - Theory and Experiment, 27 - 30 June 2022, Lecce (Italy). arXiv admin note: substantial text overlap with arXiv:2205.13952

In addition to the well-known $B \to D^{(*)} \ell \nu_\ell$ decays, semileptonic $B_s \to D_s^{(*)} \ell \nu_\ell$ processes offer the possibility to determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $\vert V_{cb}\vert$. We implement the Dispersive Matrix (DM) approach to describe the hadronic Form Factors (FFs) for the $B_s \to D_s^{(*)}$ transition in the whole kinematical range, starting from recent Lattice QCD computations at large values of the 4-momentum transfer. We extract $\vert V_{cb} \vert$ from the experimental data, obtaining $\vert V_{cb} \vert \cdot 10^3 = (41.7 \pm 1.9)$ from $B_s \to D_s \ell \nu_\ell$ and $\vert V_{cb} \vert \cdot 10^3 =(40.7 \pm 2.4)$ from $B_s \to D_s^* \ell \nu_\ell$ decays. After averaging with the values of $\vert V_{cb} \vert$ obtained from the $B \to D^{(*)}$ channels [1,2] we get $\vert V_{cb} \vert \cdot 10^3 =(41.2 \pm 0.8)$, which is compatible with the most recent inclusive estimate $\vert V_{cb} \vert_{\rm{incl}} \cdot 10^3 = 42.16 \pm 0.50$ [3] at the $1 \sigma$ level. In addition we test the Lepton Flavour Universality (LFU) by computing the $\tau / \ell$ ratios of the total decay rates (where $\ell$ is a light lepton), obtaining $R(D_s) = 0.298\,(5)$ and $R(D_s^*)= 0.250\,(6)$. We also address the issue of the $SU(3)_F$ symmetry breaking by comparing the hadronic FFs entering the semileptonic $B \to D^{(*)}$ and $B_s \to D_s^{(*)}$ channels., Comment: 31 pages, 8 figures, 3 tables. Revised and reformatted version to appear in PRD. Results and conclusions unchanged

We examine the semileptonic $B \to D^{(*)} \ell \nu_\ell$ and $B \to \pi \ell \nu_\ell$ decays adopting the unitarity-based Dispersive Matrix (DM) method, which allows to determine the shape of the relevant hadronic form factors (FFs) in their whole kinematical range, using only lattice QCD results available at large values of the 4-momentum transfer without making any assumption on their momentum dependence. Moreover, the experimental data are not used to constrain the shape of the FFs, but only to obtain our final exclusive determination of $\vert V_{cb} \vert$ and $\vert V_{ub} \vert$, namely: $\vert V_{cb} \vert \cdot 10^3 = 41.1 \pm 1.0$ and $\vert V_{ub} \vert \cdot 10^3 = 3.88 \pm 0.32$, which are consistent with the latest inclusive determinations at the $1 \sigma$ level or better. Our calculation of the FFs allows to obtain pure theoretical estimates of the $\tau / \mu$ ratios of differential decay rates, $R(D) = 0.296 \pm 0.008$ and $R(D^*) = 0.275 \pm 0.008$, which turn out to be compatible with the experimental world averages within $\simeq 1.4$ standard deviations., Comment: 6 pages, 2 figures. Contribution to the CKM Workshop 2021, 22-26 November 2021, Melbourne, Australia

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

The Cabibbo-Kobayashi-Maskawa (CKM) matrix element $\vert V_{ub}\vert$ is obtained from exclusive semileptonic $B \to \pi \ell \nu_\ell$ and $B_s \to K \ell \nu_\ell$ decays adopting the unitarity-based dispersion matrix approach for the determination of the hadronic form factors (FFs) in the whole kinematical range. We use lattice computations of the relevant susceptibilities and of the FFs in the large-$q^2$ regime in order to derive their behavior in the low-$q^2$ region without assuming any specific momentum dependence and without constraining their shape using experimental data. Then, we address the extraction of $\vert V_{ub}\vert$ from the experimental data, obtaining $\vert V_{ub}\vert = (3.62 \pm 0.47) \cdot 10^{-3}$ from $B \to \pi$ and $\vert V_{ub}\vert = (3.77 \pm 0.48) \cdot 10^{-3}$ from $B_s \to K$, which after averaging yield $\vert V_{ub}\vert = (3.69 \pm 0.34) \cdot 10^{-3}$. These results are compatible with the most recent inclusive value $\vert V_{ub} \vert_{incl} = 4.13\,(26) \cdot 10^{-3}$ at the 1$\sigma$ level. We also present purely theoretical estimates of the ratio of the $\tau/\mu$ decay rates $R^{\tau/\mu}_{\pi(K)}$, the normalized forward-backward asymmetry $\bar{\mathcal{A}}_{FB}^{\ell,\pi(K)}$ and the normalized lepton polarization asymmetry $\bar{\mathcal{A}}_{polar}^{\ell,\pi(K)}$., Comment: 46 pages, 9 figures, 9 tables. Version improved by the referee's suggestions. Conclusions unchanged. To appear in JHEP. arXiv admin note: text overlap with arXiv:2105.07851

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 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 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 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

In this work we apply the Dispersive Matrix (DM) method of Refs. [1,2] to the lattice computations of the Form Factors (FFs) entering the semileptonic $B \to D^* \ell \nu_\ell$ decays, recently produced by the FNAL/MILC Collaborations [3] at small, but non-vanishing values of the recoil variable ($w-1$). Thanks to the DM method we obtain the FFs in the whole kinematical range accessible to the decay in a completely model-independent and non-perturbative way, implementing exactly both unitarity and kinematical constraints. Using our theoretical bands of the FFs we extract $\vert V_{cb} \vert$ from the experimental data and compute the theoretical value of $R(D^*)$. Our final result for $\vert V_{cb} \vert$ reads $\vert V_{cb} \vert = (41.3 \pm 1.7) \cdot 10^{-3}$, compatible with the most recent inclusive estimate at the $0.5\sigma$ level. Moreover, we obtain the pure theoretical value $R(D^*) = 0.275 \pm 0.008$, which is compatible with the experimental world average at the $\sim 1.3 \sigma$ level., Comment: 20 pages, 5 figures, 2 tables. Revised version including appendices describing the DM method, the application of the unitarity filters and the implementation of the kinematical constraints. Results and conclusions unchanged. Matches published version in EPJC

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

The Cabibbo-Kobayashi-Maskawa (CKM) matrix element $|V_{cb}|$ is extracted from exclusive semileptonic $B \to D^{(*)}$ decays adopting a novel unitarity-based approach which allows to determine in a full non-perturbative way the relevant hadronic form factors (FFs) in the whole kinematical range. By using existing lattice computations of the $B \to D^{(*)}$ FFs at small recoil from FNAL/MILC and JLQCD Collaborations, we show that it is possible to extrapolate their behavior also at large recoil without assuming any specific momentum dependence and without constraining their shape using experimental data. Thus, we address the extraction of $|V_{cb}|$ from the experimental data on the semileptonic $B \to D^{(*)} \ell \nu_\ell$ decays, obtaining $\vert V_{cb}\vert = (41.0 \pm 1.2 ) \cdot 10^{-3}$ from $B \to D$ using as input the final FNAL/MILC lattice data for the FFs and $|V_{cb}| = (40.4 \pm 1.8 ) \cdot 10^{-3}$ from $B \to D^*$ using the preliminary JLQCD lattice data. Our result from $B \to D$ is consistent within $\sim 1$ standard deviation with the most recent inclusive determination $|V_{cb}|_{incl} = (42.00 \pm 0.65) \cdot 10^{-3}$. The resulting uncertainty is comparable with those obtained in literature using experimental data to constrain the shape of the FFs. Our result from $B \to D^*$, though consistent with $|V_{cb}|_{incl} $, is still based on preliminary lattice data for the FFs and its uncertainty is greater than the ones obtained in literature using experimental data to constrain the shape of the FFs. We investigate also the issue of Lepton Flavor Universality thanks to new theoretical estimates of the ratios $R(D^{(*)})$, namely $R(D) = 0.296(8)$ using final FNAL/MILC lattice results, and $R(D^{*}) = 0.261(20)$ using preliminary JLQCD and FNAL/MILC lattice data. Our findings differ by $\sim 1.4\sigma$ from the latest experimental determinations., Comment: 48 pages, 15 figures, 6 tables. Matches the published version in PRD

In this work we present the first non-perturbative determination of the hadronic susceptibilities that constrain the form factors entering the semileptonic $B \to D^{(*)} \ell \nu_\ell $ transitions due to unitarity and analyticity. The susceptibilities are obtained by evaluating moments of suitable two-point correlation functions obtained on the lattice. Making use of the gauge ensembles produced by the Extended Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) and with pion masses in the range $\simeq 210 - 450$ MeV, we evaluate the longitudinal and transverse susceptibilities of the vector and axial-vector polarization functions at the physical pion point and in the continuum and infinite volume limits. The ETMC ratio method is adopted to reach the physical $b$-quark mass $m_b^{phys}$. At zero momentum transfer for the $b \to c$ transition we get $\chi_{0^+}(m_b^{phys}) = 7.58\,(59) \cdot 10^{-3}$, $\chi_{1^-}(m_b^{phys}) = 6.72\,(41) \cdot 10^{-4}$ GeV$^{-2}$, $\chi_{0^-}(m_b^{phys}) = 2.58\,(17) \cdot 10^{-2}$ and $\chi_{1^+}(m_b^{phys}) = 4.69\,(30) \cdot 10^{-4}$ GeV$^{-2}$ for the scalar, vector, pseudoscalar and axial susceptibilities, respectively. In the case of the vector and pseudoscalar channels the one-particle contributions due to $B_c^*$- and $B_c$-mesons are evaluated and subtracted to improve the bounds, obtaining: $\chi_{1^-}^{sub}(m_b^{phys}) = 5.84\,(44) \cdot 10^{-4}$ GeV$^{-2}$ and $\chi_{0^-}^{sub}(m_b^{phys}) = 2.19\,(19) \cdot 10^{-2}$., Comment: 48 pages, 20 figures, 9 tables. Version to appear in PRD. arXiv admin note: text overlap with arXiv:2105.02497

In this work we discuss in detail the non-perturbative determination of the momentum dependence of the form factors entering in semileptonic decays using unitarity and analyticity constraints. The method contains several new elements with respect to previous proposals and allows to extract, using suitable two-point functions computed non-perturbatively, the form factors at low momentum transfer $q^2$ from those computed explicitly on the lattice at large $q^2$, without any assumption about their $q^2$-dependence. The approach will be very useful for exclusive semileptonic $B$-meson decays, where the direct calculation of the form factors at low $q^2$ is particularly difficult due to large statistical fluctuations and discretisation effects. As a testing ground we apply our approach to the semileptonic $D \to K \ell \nu_\ell$ decay, where we can compare the results of the unitarity approach to the explicit direct lattice calculation of the form factors in the full $q^2$-range. We show that the method is very effective and that it allows to compute the form factors with rather good precision., Comment: 66 pages, 13 figures and 10 tables. Minor changes matched to the version accepted by PRD

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 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 novel strategy to renormalize lattice operators in QCD+QED, including first order QED corrections to the non-perturbative evaluation of QCD renormalization constants. Our procedure takes systematically into account the mixed non-factorizable QCD+QED effects which were neglected in previous calculations, thus significantly reducing the systematic uncertainty on renormalization corrections. The procedure is presented here in the RI'-MOM scheme, but it can be applied to other schemes (e.g. RI-SMOM) with appropriate changes. We discuss the application of this strategy to the calculation of the leading isospin breaking corrections to the leptonic decay rates $\Gamma(\pi_{\mu 2})$ and $\Gamma(K_{\mu 2})$, evaluated for the first time on the lattice. The precision in the matching to the $W$-regularization scheme is improved to $\mathcal{O}(\alpha_{em}\alpha_s(M_W))$ with respect to previous calculations. Finally, we show the updated precise result obtained for the Cabibbo-Kobayashi-Maskawa matrix element $|V_{us}|$., Comment: Presented at the 37th International Symposium on Lattice Field Theory (Lattice 2019), 16-22 June 2019, Wuhan, China

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

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.

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

We present a lattice calculation of the leading-order electromagnetic and strong isospin-breaking corrections to the hadronic vacuum polarization (HVP) contribution to the anomalous magnetic moment of the muon. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses between $\simeq 210$ and $\simeq 450$ MeV. The results are obtained adopting the RM123 approach in the quenched-QED approximation, which neglects the charges of the sea quarks. Quark disconnected diagrams are not included. After the extrapolations to the physical pion mass and to the continuum and infinite-volume limits the contributions of the light, strange and charm quarks are respectively equal to $\delta a_\mu^{\rm HVP}(ud) = 7.1 ~ (2.5) \cdot 10^{-10}$, $\delta a_\mu^{\rm HVP}(s) = -0.0053 ~ (33) \cdot 10^{-10}$ and $\delta a_\mu^{\rm HVP}(c) = 0.0182 ~ (36) \cdot 10^{-10}$. At leading order in $\alpha_{em}$ and $(m_d - m_u) / \Lambda_{QCD}$ we obtain $\delta a_\mu^{\rm HVP}(udsc) = 7.1 ~ (2.9) \cdot 10^{-10}$, which is currently the most accurate determination of the isospin-breaking corrections to $a_\mu^{\rm HVP}$., Comment: 23 pages, 7 figures, 5 tables. Version to appear in PRD. A bug in the update of the strange and charm contributions is removed and an extended discussion on the identification of the ground-state is included. arXiv admin note: text overlap with arXiv:1808.00887, arXiv:1707.03019

We present the results of the first lattice QCD calculation of the $K \to \pi$ matrix elements of the chromomagnetic operator $O_{CM} = g\, \bar s\, \sigma_{\mu\nu} G_{\mu\nu} d$, which appears in the effective Hamiltonian describing $\Delta S = 1$ transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been determined non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing. Our result for the B-parameter of the chromomagnetic operator at the physical pion and kaon point is $B_{CMO}^{K \pi} = 0.273 ~ (70)$, while in the SU(3) chiral limit we obtain $B_{CMO} = 0.072 ~ (22)$. Our findings are significantly smaller than the model-dependent estimate $B_{CMO} \sim 1 - 4$, currently used in phenomenological analyses, and improve the uncertainty on this important phenomenological quantity., Comment: 20 pages, 4 figures, 2 table. Refined SU(3) ChPT analysis with no changes in the final result. Version to appear in PRD

The leading-order electromagnetic and strong isospin-breaking corrections to the ratio of $K_{\mu 2}$ and $\pi_{\mu 2}$ decay rates are evaluated for the first time on the lattice, following a method recently proposed. The lattice results are obtained using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks. Systematics effects are evaluated and the impact of the quenched QED approximation is estimated. Our result for the correction to the tree-level $K_{\mu 2} / \pi_{\mu 2}$ decay ratio is $-1.22\,(16) \%$ to be compared to the estimate $-1.12\,(21) \%$ based on Chiral Perturbation Theory and adopted by the Particle Data Group., Comment: 5 pages, 6 figures; extended supplemental material with 1 table and 1 figure, results unchanged

We present a study of the isospin-breaking (IB) corrections to pseudoscalar (PS) meson masses using the gauge configurations produced by the ETM Collaboration with $N_f=2+1+1$ dynamical quarks at three lattice spacings varying from 0.089 to 0.062 fm. Our method is based on a combined expansion of the path integral in powers of the small parameters $(\widehat{m}_d - \widehat{m}_u)/\Lambda_{QCD}$ and $\alpha_{em}$, where $\widehat{m}_f$ is the renormalized quark mass and $\alpha_{em}$ the renormalized fine structure constant. We obtain results for the pion, kaon and $D$-meson mass splitting; for the Dashen's theorem violation parameters $\epsilon_\gamma(\overline{\mathrm{MS}}, 2~\mbox{GeV})$, $\epsilon_{\pi^0}$, $\epsilon_{K^0}(\overline{\mathrm{MS}}, 2~\mbox{GeV})$; for the light quark masses $(\widehat{m}_d - \widehat{m}_u)(\overline{\mathrm{MS}}, 2~\mbox{GeV})$, $(\widehat{m}_u / \widehat{m}_d)(\overline{\mathrm{MS}}, 2~\mbox{GeV})$; for the flavour symmetry breaking parameters $R(\overline{\mathrm{MS}}, 2~\mbox{GeV})$ and $Q(\overline{\mathrm{MS}}, 2~\mbox{GeV})$ and for the strong IB effects on the kaon decay constants., Comment: 8 pages, 2 figures, talk delivered at the XXXV International Symposium on Lattice Field Theory, Granada (Spain), 18-24 June 2017

We present a lattice calculation of the Hadronic Vacuum Polarization (HVP) contribution of the strange and charm quarks to the anomalous magnetic moment of the muon including leading-order electromagnetic (e.m.) corrections. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The strange and charm quark masses are tuned at their physical values. Neglecting disconnected diagrams and after the extrapolations to the physical pion mass and to the continuum limit we obtain: $a_\mu^s(\alpha_{em}^2) = (53.1 \pm 2.5) \cdot 10^{-10}$, $a_\mu^s(\alpha_{em}^3) = (-0.018 \pm 0.011) \cdot 10^{-10}$ and $a_\mu^c(\alpha_{em}^2) = (14.75 \pm 0.56) \cdot 10^{-10}$, $a_\mu^c(\alpha_{em}^3) = (-0.030 \pm 0.013) \cdot 10^{-10}$ for the strange and charm contributions, respectively., Comment: 8 pages, 3 figures, talk delivered at the XXXV Int'l Symposium on Lattice Field Theory (LAT17), Granada (Spain), June 18-24, 2017

We present a lattice calculation of the Hadronic Vacuum Polarization (HVP) contribution of the strange and charm quarks to the anomalous magnetic moment of the muon including leading-order electromagnetic corrections. We employ the gauge configurations generated by the European Twisted Mass Collaboration (ETMC) with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The strange and charm quark masses are tuned at their physical values. Neglecting disconnected diagrams and after the extrapolations to the physical pion mass and to the continuum limit we obtain: $a_\mu^s(\alpha_{em}^2) = (53.1 \pm 2.5) \cdot 10^{-10}$, $a_\mu^s(\alpha_{em}^3) = (-0.018 \pm 0.011) \cdot 10^{-10}$ and $a_\mu^c(\alpha_{em}^2) = (14.75 \pm 0.56) \cdot 10^{-10}$, $a_\mu^c(\alpha_{em}^3) = (-0.030 \pm 0.013) \cdot 10^{-10}$ for the strange and charm contributions, respectively., Comment: 34 pages, 10 figures, 5 tables; version to appear in JHEP

We present a lattice computation of the isospin-breaking corrections to pseudoscalar meson masses using the gauge configurations produced by the European Twisted Mass collaboration with $N_f = 2 + 1 + 1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082$ and $0.089$ fm) with pion masses in the range $M_\pi \simeq 210 - 450$ MeV. The strange and charm quark masses are tuned at their physical values. We adopt the RM123 method based on the combined expansion of the path integral in powers of the $d$- and $u$-quark mass difference ($\widehat{m}_d - \widehat{m}_u$) and of the electromagnetic coupling $\alpha_{em}$. Within the quenched QED approximation, which neglects the effects of the sea-quark charges, and after the extrapolations to the physical pion mass and to the continuum and infinite volume limits, we provide results for the pion, kaon and (for the first time) charmed-meson mass splittings, for the prescription-dependent parameters $\epsilon_{\pi^0}$, $\epsilon_\gamma(\overline{MS}, 2~\mbox{GeV})$, $\epsilon_{K^0}(\overline{MS}, 2~\mbox{GeV})$, related to the violations of the Dashen's theorem, and for the light quark mass difference $(\widehat{m}_d - \widehat{m}_u)(\overline{MS}, 2~\mbox{GeV})$., Comment: 47 pages, 20 figures, 4 tables; comments on QED and QCD splitting prescriptions added; version to appear in PRD

In Carrasco et al. we have recently proposed a method to calculate $O(e^2)$ electromagnetic corrections to leptonic decay widths of pseudoscalar mesons. The method is based on the observation that the infrared divergent contributions (that appear at intermediate stages of the calculation and that cancel in physical quantities thanks to the Bloch-Nordsieck mechanism) are universal, i.e. depend on the charge and the mass of the meson but not on its internal structure. In this talk we perform a detailed analysis of the finite-volume effects associated with our method. In particular we show that also the leading $1/L$ finite-volume effects are universal and perform an analytical calculation of the finite-volume leptonic decay rate for a point-like meson., Comment: Corrected a typo affecting the 1/L^3 term. The typo has been discovered from the authors of arXiv:2109.05002

We demonstrate that the leading and next-to-leading finite-volume effects in the evaluation of leptonic decay widths of pseudoscalar mesons at $O(\alpha)$ are universal, i.e. they are independent of the structure of the meson. This is analogous to a similar result for the spectrum but with some fundamental differences, most notably the presence of infrared divergences in decay amplitudes. The leading non-universal, structure-dependent terms are of $O(1/L^2)$ (compared to the $O(1/L^3)$ leading non-universal corrections in the spectrum). We calculate the universal finite-volume effects, which requires an extension of previously developed techniques to include a dependence on an external three-momentum (in our case, the momentum of the final state lepton). The result can be included in the strategy proposed in Ref.\,\cite{Carrasco:2015xwa} for using lattice simulations to compute the decay widths at $O(\alpha)$, with the remaining finite-volume effects starting at order $O(1/L^2)$. The methods developed in this paper can be generalised to other decay processes, most notably to semileptonic decays, and hence open the possibility of a new era in precision flavour physics.

Electromagnetic effects in the leptonic decay rates $\pi^+ \to \mu^+ \nu$ and $K^+ \to \mu^+ \nu$ are evaluated for the first time on the lattice. Following a method recently proposed in Ref. [1] the emission of virtual photons at leading order in the electromagnetic coupling is evaluated on the lattice and the infrared divergence computed for a point-like meson at finite lattice volume is subtracted. The physical decay rate is then obtained by adding the emission of real and virtual photons regularised with a photon mass. Using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks the feasibility of our approach is demonstrated. Preliminary results for the electromagnetic corrections to charged (neutral) pion and kaon masses as well as to the leptonic decay rates of charged pions and kaons are presented., Comment: 7 pages, 4 figures, 1 table, proceedings of the 34th annual International Symposium on Lattice Field Theory, 24-30 July 2016, University of Southampton (UK)

We present our study of the renormalization of the chromomagnetic operator,O(CM), which appears in the effective Hamiltonian describing Delta S = 1 transitions in and beyond the Standard Model. We have computed, perturbatively to one-loop, the relevant Green's functions with two (quark-quark) and three (quark-quark-gluon) external fields, at nonzero quark masses, using both the lattice and dimensional regularizations. The perturbative computation on the lattice is carried out using the maximally twisted-mass action for the fermions, while for the gluons we employed the Symanzik improved gauge action for different sets of values of the Symanzik coefficients. We have identified all the operators which can possibly mix with O(CM), including lower dimensional and non gauge invariant operators, and we have calculated those elements of the mixing matrix which are relevant for the renormalization of O(CM). We have also performed numerical lattice calculations to determine non-perturbatively the mixings of the chromomagnetic operator with lower dimensional operators, through proper renormalization conditions. For the first time the 1/a**2-divergent mixing of the chromomagnetic operator with the scalar density has been determined non-perturbatively with high precision. Moreover, the 1/a-divergent mixing with the pseudoscalar density, due to the breaking of parity within the twisted-mass regularization of QCD, has been calculated non-perturbatively and found to be smaller than its one-loop perturbative estimate. The QCD simulations have been carried out using the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2 + 1 + 1 dynamical quarks, which include in the sea, besides two light mass degenerate quarks, also the strange and charm quarks with masses close to their physical values., Comment: 28 pages, 13 figures, 7 tables

In this paper, for the first time a method is proposed to compute electromagnetic effects in hadronic processes using lattice simulations. The method can be applied, for example, to the leptonic and semileptonic decays of light or heavy pseudoscalar mesons. For these quantities the presence of infrared divergences in intermediate stages of the calculation makes the procedure much more complicated than is the case for the hadronic spectrum, for which calculations already exist. In order to compute the physical widths, diagrams with virtual photons must be combined with those corresponding to the emission of real photons. Only in this way do the infrared divergences cancel as first understood by Bloch and Nordsieck in 1937. We present a detailed analysis of the method for the leptonic decays of a pseudoscalar meson. The implementation of our method, although challenging, is within reach of the present lattice technology., Comment: 38 pages, 14 figures

The Chromomagnetic operator (CMO) mixes with a large number of operators under renormalization. We identify which operators can mix with the CMO, at the quantum level. Even in dimensional regularization (DR), which has the simplest mixing pattern, the CMO mixes with a total of 9 other operators, forming a basis of dimension-five, Lorentz scalar operators with the same flavor content as the CMO. Among them, there are also gauge noninvariant operators; these are BRST invariant and vanish by the equations of motion, as required by renormalization theory. On the other hand using a lattice regularization further operators with $d \leq 5$ will mix; choosing the lattice action in a manner as to preserve certain discrete symmetries, a minimul set of 3 additional operators (all with $d<5$) will appear. In order to compute all relevant mixing coefficients, we calculate the quark-antiquark (2-pt) and the quark-antiquark-gluon (3-pt) Green's functions of the CMO at nonzero quark masses. These calculations were performed in the continuum (dimensional regularization) and on the lattice using the maximally twisted mass fermion action and the Symanzik improved gluon action. In parallel, non-perturbative measurements of the $K-\pi$ matrix element are being performed in simulations with 4 dynamical ($N_f = 2+1+1$) twisted mass fermions and the Iwasaki improved gluon action., Comment: 7 pages, 1 figure, 3 tables, LATTICE2014 proceedings

We present preliminary results of the first lattice QCD calculation of the K -> pi matrix elements of the chromomagnetic operator O_{CM}=g sbar sigma_{munu} G_{munu} d, which appears in the effective Hamiltonian describing Delta S=1 transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one-loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been computed non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with N_f=2+1+1 dynamical quarks at three values of the lattice spacing. Our preliminary result for the B-parameter of the chromomagnetic operator is B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~1-4 currently used in phenomenological analyses., Comment: 7 pages, 3 figures, LATTICE2014 proceedings