Search

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 in isospin symmetric QCD. We employ the gauge configurations generated by the Extended Twisted Mass Collaboration (ETMC) with $N_f = 2 + 1 + 1$ flavors of Wilson-clover twisted-mass quarks at five lattice spacings and at values of the quark mass parameters that are close and/or include the isospin symmetric QCD point of interest. After computing the small corrections necessary to precisely match this point, and carrying out an extrapolation to the continuum limit based on the data at lattice spacings $a \simeq 0.049, 0.057, 0.068, 0.080$~fm and spatial lattice sizes up to $L \simeq 7.6$~fm, we obtain $a_\mu^{\rm HVP}(s) = (53.57 \pm 0.63) \times 10^{-10}$ and $a_\mu^{\rm HVP}(c) = (14.56 \pm 0.13) \times 10^{-10}$, for the quark-connected strange and charm contributions, respectively. Our findings agree well with the corresponding results by other lattice groups., Comment: 33 pages, 18 figures, 3 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 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 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 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 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 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 calculation of the Mellin moments of the transverse quark spin densities in the nucleon using lattice QCD. The densities are extracted from the unpolarized and transversity generalized form factors extrapolated to the continuum limit using three $N_f=2+1+1$ twisted mass fermion gauge ensembles simulated with physical quark masses and spanning three lattice spacings. The first moment of transversely polarized quarks in an unpolarized nucleon shows an interesting distortion, which can be traced back to the sharp falloff of the transversity generalized form factor $\bar{B}_{Tn0}(t)$. The isovector tensor anomalous magnetic moment is determined to be $\kappa_T=1.051(94)$, which confirms a negative and large Boer-Mulders function, $h_1^{\perp}$, in the nucleon.

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 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 a lattice QCD computation and preliminary results for the leptonic decay constants of the pseudoscalar mesons $K$, $D$ and $D_{s}$ in the isosymmetric QCD limit. The computation is based on simulations of $N_f = 2+1+1$ dynamical quarks performed by the Extended Twisted Mass Collaboration (ETMC), where the light, strange and charm quark masses are all tuned at their physical values. We also present preliminary unitarity checks for the first and second rows of the Cabibbo-Kobayashi-Maskawa matrix., Comment: 9 pages, 3 figures, 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@MIT

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 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 a determination of the ratio of kaon and pion leptonic decay constants in isosymmetric QCD (isoQCD), $f_K / f_\pi$, making use 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, including configurations close to the physical point for all dynamical flavors. The simulations are carried out at three values of the lattice spacing ranging from $\sim 0.068$ to $\sim 0.092$ fm with linear lattice size up to $L \sim 5.5$~fm. The scale is set by the PDG value of the pion decay constant, $f_\pi^{isoQCD} = 130.4~(2)$ MeV, at the isoQCD pion point, $M_\pi^{isoQCD} = 135.0~(2)$ MeV, obtaining for the gradient-flow (GF) scales the values $w_0 = 0.17383~(63)$ fm, $\sqrt{t_0} = 0.14436~(61)$ fm and $t_0 / w_0 = 0.11969~(62)$ fm. The data are analyzed within the framework of SU(2) Chiral Perturbation Theory (ChPT) without resorting to the use of renormalized quark masses. At the isoQCD kaon point $M_K^{isoQCD} = 494.2~(4)$ MeV we get $(f_K / f_\pi)^{isoQCD} = 1.1995~(44)$, where the error includes both statistical and systematic uncertainties. Implications for the Cabibbo-Kobayashi-Maskawa (CKM) matrix element $|V_{us}|$ and for the first-row CKM unitarity are discussed., Comment: 68 pages, 14 figures, 12 tables. Version to appear in PRD

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 results for the moments of nucleon isovector vector and axial generalised parton distribution functions computed within lattice QCD. Three ensembles of maximally twisted mass clover-improved fermions simulated with a physical value of the pion mass are analyzed. Two of these ensembles are generated using two degenerate light quarks. A third ensemble is used having, in addition to the light quarks, strange and charm quarks in the sea. A careful analysis of the convergence to the ground state is carried out that is shown to be essential for extracting the correct nucleon matrix elements. This allows a controlled determination of the unpolarised, helicity and tensor second Mellin moments. The vector and axial-vector generalised form factors are also computed as a function of the momentum transfer square up to about 1 GeV$^2$. The three ensembles allow us to check for unquenching effects and to assess lattice finite volume effects., Comment: Published version. Typos corrected, references and figures added. 22 pages, 17 figures

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

We present an investigation of the electromagnetic pion form factor, $F_\pi(Q^2)$, at small values of the four-momentum transfer $Q^2$ ($\lesssim 0.25$ GeV$^2$), based on the gauge configurations generated by European Twisted Mass Collaboration with $N_f = 2$ twisted-mass quarks at maximal twist including a clover term. Momentum is injected using non-periodic boundary conditions and the calculations are carried out at a fixed lattice spacing ($a \simeq 0.09$ fm) and with pion masses equal to its physical value, 240 MeV and 340 MeV. Our data are successfully analyzed using Chiral Perturbation Theory at next-to-leading order in the light-quark mass. For each pion mass two different lattice volumes are used to take care of finite size effects. Our final result for the squared charge radius is $\langle r^2 \rangle_\pi = 0.443~(29)$ fm$^2$, where the error includes several sources of systematic errors except the uncertainty related to discretization effects. The corresponding value of the SU(2) chiral low-energy constant $\overline{\ell}_6$ is equal to $\overline{\ell}_6 = 16.2 ~ (1.0)$., Comment: 23 pages, 9 figures, 7 tables

We present results on the light, strange and charm nucleon scalar and tensor charges from lattice QCD, using simulations with $N_f=2$ flavors of twisted mass Clover-improved fermions with a physical value of the pion mass. Both connected and disconnected contributions are included, enabling us to extract the isoscalar, strange and charm charges for the first time directly at the physical point. Furthermore, the renormalization is computed non-perturbatively for both isovector and isoscalar quantities. We investigate excited state effects by analyzing several sink-source time separations and by employing a set of methods to probe ground state dominance. Our final results for the scalar charges are $g_S^u = 5.20(42)(15)(12)$, $g_S^d = 4.27(26)(15)(12)$, $g_S^s=0.33(7)(1)(4)$, $g_S^c=0.062(13)(3)(5)$ and for the tensor charges $g_T^u = 0.782(16)(2)(13)$, $g_T^d = -0.219(10)(2)(13)$, $g_T^s=-0.00319(69)(2)(22)$, $g_T^c=-0.00263(269)(2)(37)$ in the $\overline{\rm MS}$ scheme at 2~GeV. The first error is statistical, the second is the systematic error due to the renormalization and the third the systematic arising from possible contamination due to the excited states., Comment: 20 pages and 13 figures

We present results for the isospin-0 $\pi\pi$ s-wave scattering length calculated in twisted mass lattice QCD. We use three $N_f = 2$ ensembles with unitary pion mass at its physical value, 240~MeV and 330~MeV respectively. We also use a large set of $N_f = 2 + 1 +1$ ensembles with unitary pion masses varying in the range of 230~MeV - 510~MeV at three different values of the lattice spacing. A mixed action approach with the Osterwalder-Seiler action in the valence sector is adopted to circumvent the complications arising from isospin symmetry breaking of the twisted mass quark action. Due to the relatively large lattice artefacts in the $N_f = 2 + 1 +1$ ensembles, we do not present the scattering lengths for these ensembles. Instead, taking the advantage of the many different pion masses of these ensembles, we qualitatively discuss the pion mass dependence of the scattering properties of this channel based on the results from the $N_f = 2 + 1 +1$ ensembles. The scattering length is computed for the $N_f = 2$ ensembles and the chiral extrapolation is performed. At the physical pion mass, our result $M_\pi a^\mathrm{I=0}_0 = 0.198(9)(6)$ agrees reasonably well with various experimental measurements and theoretical predictions., Comment: 7 pages, 2 figures, talk presented at the 34th annual International Symposium on Lattice Field Theory (Lattice2016), 24-30 July 2016, University of Southampton, UK

We present results for the isospin-0 $\pi\pi$ s-wave scattering length calculated with Osterwalder-Seiler valence quarks on Wilson twisted mass gauge configurations. We use three $N_f = 2$ ensembles with unitary (valence) pion mass at its physical value (250$\sim$MeV), at 240$\sim$MeV (320$\sim$MeV) and at 330$\sim$MeV (400$\sim$MeV), respectively. By using the stochastic Laplacian Heaviside quark smearing method, all quark propagation diagrams contributing to the isospin-0 $\pi\pi$ correlation function are computed with sufficient precision. The chiral extrapolation is performed to obtain the scattering length at the physical pion mass. Our result $M_\pi a^\mathrm{I=0}_0 = 0.198(9)(6)$ agrees reasonably well with various experimental measurements and theoretical predictions. Since we only use one lattice spacing, certain systematics uncertainties, especially those arising from unitary breaking, are not controlled in our result., Comment: 21 pages, 5 figures, 6 tables

We present precise lattice computations for the b-quark mass, the quark mass ratios mb/mc and mb/ms as well as the leptonic B-decay constants. We employ gauge configurations with four dynamical quark flavors, up/down, strange and charm, at three values of the lattice spacing (a ~ 0.06 - 0.09 fm) and for pion masses as low as 210 MeV. Interpolation in the heavy quark mass to the bottom quark point is performed using ratios of physical quantities computed at nearby quark masses exploiting the fact that these ratios are exactly known in the static quark mass limit. Our results are also extrapolated to the physical pion mass and to the continuum limit and read: mb(MSbar, mb) = 4.26(10) GeV, mb/mc = 4.42(8), mb/ms = 51.4(1.4), fBs = 229(5) MeV, fB = 193(6) MeV, fBs/fB = 1.184(25) and (fBs/fB)/(fK/fpi) = 0.997(17)., Comment: Version to appear in PRD. Added comments to simulation setup and error budget discussion. 1+20 pages, 9 figures

We present physics results from simulations of QCD using $N_f = 2$ dynamical Wilson twisted mass fermions at the physical value of the pion mass. These simulations were enabled by the addition of the clover term to the twisted mass quark action. We show evidence that compared to previous simulations without this term, the pion mass splitting due to isospin breaking is almost completely eliminated. Using this new action, we compute the masses and decay constants of pseudoscalar mesons involving the dynamical up and down as well as valence strange and charm quarks at one value of the lattice spacing, $a \approx 0.09$ fm. Further, we determine renormalized quark masses as well as their scale-independent ratios, in excellent agreement with other lattice determinations in the continuum limit. In the baryon sector, we show that the nucleon mass is compatible with its physical value and that the masses of the $\Delta$ baryons do not show any sign of isospin breaking. Finally, we compute the electron, muon and tau lepton anomalous magnetic moments and show the results to be consistent with extrapolations of older ETMC data to the continuum and physical pion mass limits. We mostly find remarkably good agreement with phenomenology, even though we cannot take the continuum and thermodynamic limits., Comment: 45 pages, 15 figures

We present results on the nucleon scalar, axial and tensor charges as well as on the momentum fraction, and the helicity and transversity moments. The pion momentum fraction is also presented. The computation of these key observables is carried out using lattice QCD simulations at a physical value of the pion mass. The evaluation is based on gauge configurations generated with two degenerate sea quarks of twisted mass fermions with a clover term. We investigate excited states contributions with the nucleon quantum numbers by analyzing three sink-source time separations. We find that, for the scalar charge, excited states contribute significantly and to a less degree to the nucleon momentum fraction and helicity moment. Our analysis yields a value for the nucleon axial charge agrees with the experimental value and we predict a value of 1.027(62) in the $\overline{\text{MS}}$ scheme at 2 GeV for the isovector nucleon tensor charge directly at the physical point. The pion momentum fraction is found to be $\langle x\rangle_{u-d}^{\pi^\pm}=0.214(15)(^{+12}_{-9})$ in the $\overline{\rm MS}$ at 2 GeV., Comment: The experimental value of the helicity plotted in Fig.17 is corrected. ETM Collaboration. 25 pages, 19 figures. Version accepted for PRD (volume 92, eid 114513)

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

We present unquenched lattice QCD results for the matrix elements of four-fermion operators relevant to the description of the neutral K and D mixing in the Standard Model and its extensions. We have employed simulations with Nf = 2 + 1 + 1 dynamical sea quarks at three values of the lattice spacings in the interval 0.06 - 0.09 fm and pseudoscalar meson masses in the range 210 - 450 MeV. Our results are extrapolated to the continuum limit and to the physical pion mass. Renormalization constants have been determined non-perturbatively in the RI-MOM scheme. In particular, for the Kaon bag-parameter, which is relevant for the \overline{K}^0-K^0 mixing in the Standard Model, we obtain B_K^{RGI} = 0.717(24)., Comment: Added comments to error budget discussion; fig.19 corrected. Version to appear in PRD

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

We present a lattice QCD calculation of the pseudoscalar decay constants fK, fD and fDs performed 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 values in the real world. The simulations are based on a unitary setup for the two light mass-degenerate quarks and on a mixed action approach for the strange and charm quarks. We use data simulated at three different values of the lattice spacing in the range 0.06 - 0.09 fm and at pion masses in the range 210 - 450 MeV. Our main results are: fK+ / fpi+ = 1.184 (16), fK+ = 154.4 (2.0) MeV, which incorporate the leading strong isospin breaking correction due to the up- and down-quark mass difference, and fK = 155.0 (1.9) MeV, fD = 207.4 (3.8) MeV, fDs = 247.2 (4.1) MeV, fDs / fD = 1.192 (22) and (fDs / fD) / (fK / fpi) = 1.003 (14) obtained in the isospin symmetric limit of QCD. Combined with the experimental measurements of the leptonic decay rates of kaon, pion, D- and Ds-mesons our results lead to the following determination of the CKM matrix elements: |Vus| = 0.2269 (29), |Vcd| = 0.2221 (67) and |Vcs| = 1.014 (24). Using the latest value of |Vud| from superallowed nuclear beta decays the unitarity of the first row of the CKM matrix is fulfilled at the permille level., Comment: 20 pp., 4 figures; revised version to appear in PRD; improved calculation of IB effects for fK+; minor changes in the final values. arXiv admin note: text overlap with arXiv:1403.4504

In this contribution, results from $N_f=2$ lattice QCD simulations at one lattice spacing using twisted mass fermions with a clover term at the physical pion mass are presented. The mass splitting between charged and neutral pions (including the disconnected contribution) is shown to be around $20(20)~\mathrm{MeV}$. Further, a first measurement using the clover twisted mass action of the average momentum fraction of the pion is given. Finally, an analysis of pseudoscalar meson masses and decay constants is presented involving linear interpolations in strange and charm quark masses. Matching to meson mass ratios allows the calculation of quark mass ratios: $\mu_s/\mu_l=27.63(13)$, $\mu_c/\mu_l=339.6(2.2)$ and $\mu_c/\mu_s=12.29(10)$. From this mass matching the quantities $f_K=153.9(7.5)~\mathrm{MeV}$, $f_D=219(11)~\mathrm{MeV}$, $f_{D_s}=255(12)~\mathrm{MeV}$ and $M_{D_s}=1894(93)~\mathrm{MeV}$ are determined without the application of finite volume or discretization artefact corrections and with errors dominated by a preliminary estimate of the lattice spacing.

We discuss a lattice QCD computation of the $B$-meson decay constants by the ETM collaboration where suitable ratios allow to reach the bottom quark sector by combining simulations around the charm-quark mass with an exactly known static limit. The different steps involved in this ratio method are discussed together with an account of the assessment of various systematic effects. A comparison of results from simulations with two and four flavour dynamical quarks is presented., Comment: 8 pages, 4 figures. Proceedings of the 8th International Workshop on the CKM Unitarity Triangle (CKM 2014), Vienna, Austria, September 8-12, 2014

We present precision lattice calculations of the pseudoscalar decay constants of the charmed sector as well as determinations of the bottom quark mass and its ratio to the charm quark mass. We employ Nf=2+1+1 dynamical quark gauge configurations generated by the European Twisted Mass Collaboration, using data at three values of the lattice spacing and pion masses as low as 210 MeV. Strange and charm sea quark masses are close to their physical values., Comment: 7 pages, Contribution to the Proceedings of the 37th International Conference on High Energy Physics (ICHEP2014), Valencia, Spain, July 2-9, 2014

We present the first unquenched lattice QCD results for the bag parameters controlling the short distance contribution to D meson oscillations in the Standard Model and beyond. We have used the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2 dynamical quarks, at four lattice spacings and light meson masses in the range 280-500 MeV. Renormalization is carried out non-perturbatively with the RI-MOM method. The bag-parameter results have been used to constrain New Physics effects in D-Dbar mixing, to put a lower bound to the generic New Physics scale and to constrain off-diagonal squark mass terms for TeV-scale Supersymmetry., Comment: 18 pages, 4 figures, 10 tables, updated to match the version accepted for publication in Physical Review D. This article draws form arXiv:1308.1851 which is on B-Bbar mixing instead of D-Dbar mixing

We present a lattice QCD calculation of the up, down, strange and charm quark masses performed 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. The simulations are based on a unitary setup for the two light quarks and on a mixed action approach for the strange and charm quarks. The analysis uses data at three values of the lattice spacing and pion masses in the range 210 - 450 MeV, allowing for accurate continuum limit and controlled chiral extrapolation. The quark mass renormalization is carried out non-perturbatively using the RI-MOM method. The results for the quark masses converted to the bar{MS} scheme are: mud(2 GeV) = 3.70(17) MeV, ms(2 GeV) = 99.6(4.3) MeV and mc(mc) = 1.348(46) GeV. We obtain also the quark mass ratios ms/mud = 26.66(32) and mc/ms = 11.62(16). By studying the mass splitting between the neutral and charged kaons and using available lattice results for the electromagnetic contributions, we evaluate mu/md = 0.470(56), leading to mu = 2.36(24) MeV and md = 5.03(26) MeV., Comment: 61 pages, 22 figures, 21 tables. Final version to appear in NPB

Taking inspiration from lattice QCD data, we argue that a finite non-perturbative contribution to the quark mass is generated as a consequence of the dynamical phenomenon of spontaneous chiral symmetry breaking, in turn triggered by the explicitly breaking of chiral symmetry induced by the critical Wilson term in the action. In pure lattice QCD this mass term cannot be separated from the unavoidably associated linearly divergent contribution. However, if QCD is enlarged to a theory where also a scalar field is present, coupled to an SU(2) doublet of fermions via a Yukawa and a Wilson-like term, then in the phase where the scalar field takes a non-vanishing expectation value, a dynamically generated and "naturally" light fermion mass (numerically unrelated to the expectation value of the scalar field) is conjectured to emerge at a critical value of the Yukawa coupling where the symmetry of the model is maximally enhanced. Masses dynamically generated in this way display a natural hierarchy according to which the stronger is the strongest of the interactions the fermion is subjected to the larger is its mass., Comment: Added more information in Fig 1. Added Fig. 10. Added an extra Appendix. Restructured a few sentences according to referee suggestions. Corrected a few misprints. All results unchanged. Now 50 pages and 10 Figures

We study matrix elements of the "chromomagnetic" operator on the lattice. This operator is contained in the strangeness-changing effective Hamiltonian which describes electroweak effects in the Standard Model and beyond. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with other operators of equal and lower dimensionality, including also non gauge invariant quantities; it is thus quite a challenge to extract from lattice simulations a clear signal for the hadronic matrix elements of this operator. We compute all relevant mixing coefficients to one loop in lattice perturbation theory; this necessitates calculating both 2-point (quark-antiquark) and 3-point (gluon-quark-antiquark) Green's functions at nonzero quark masses. We use the twisted mass lattice formulation, with Symanzik improved gluon action. For a comprehensive presentation of our results, along with detailed explanations and a more complete list of references, we refer to our forthcoming publication [1]., Comment: 7 pages, 1 figure. Talk presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

In this contribution, a first look at simulations using maximally twisted mass Wilson fermions at the physical point is presented. A lattice action including clover and twisted mass terms is presented and the Monte Carlo histories of one run with two mass-degenerate flavours at a single lattice spacing are shown. Measurements from the light and heavy-light pseudoscalar sectors are compared to previous $N_f = 2$ results and their phenomenological values. Finally, the strategy for extending simulations to $N_f = 2 + 1 + 1$ is outlined., Comment: presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a lattice QCD calculation of the pseudoscalar decay constants $f_K$, $f_D$ and $f_{D_s}$ performed by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical fermions. We simulated at three different values of the lattice spacing, the smallest being approximately $0.06fm$, and with pion masses as small as $210$MeV. Our main results are: $f_{K^+}/f_{\pi^+}=1.183(17)$, $f_{K^+}=154.4(2.1)$MeV, $f_{D_s}=242.1(8.3)$MeV, $f_D=201.9(8.0)$MeV, $f_{D_s}/f_D=1.199(25)$ and $(f_{D_s}/f_D) / (f_K/f_\pi) = 1.005(15)$., Comment: 7 pages, 5 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a lattice QCD determination of the $b$-quark mass and of the $f_{B_s}$ and $f_B$ decay constants performed with $N_f = 2 + 1 + 1$ twisted mass Wilson fermions. We have used simulations at three values of the lattice spacing generated by ETMC with pion masses ranging from 210 to 440 MeV. To obtain physical quantities we performed a combined chiral and continuum limit and an extrapolation in the heavy quark mass from the charm to the $b$-quark region using suitable ratios calculated at nearby heavy-quark masses having an exactly known static limit. Our results are: $m_b(m_b) = 4.29 (13)$ GeV, $f_B = 196 (9)$ MeV, $f_{B_s} = 235 (9)$ MeV, $f_{B_s} / f_B = 1.201 (25)$, $(f_{B_s}/f_B)/(f_K/f_\pi) = 1.007 (16)$ and $(f_{B_s}/f_B)/(f_{D_s}/f_D) = 1.008 (13)$., Comment: 7 pages, 7 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a lattice QCD determination of the average up-down, strange and charm quark masses based on simulations performed by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical fermions. We simulated at three different values of the lattice spacing, the smallest being approximately $0.06fm$, and with pion masses as small as $210 \text{MeV}$. Our results are: $m_{ud}(2\text{GeV})=3.70(17)\text{MeV}$, $m_s(2\text{GeV})=99.2(3.9)\text{MeV}$, $m_c(m_c)=1.350(49)\text{GeV}$, $m_s/m_{ud}=26.64(30)$ and $m_c/m_s=11.65(12)$.

We present an accurate lattice QCD computation of the b-quark mass, the B and Bs decay constants, the B-mixing bag-parameters for the full four-fermion operator basis, as well as estimates for \xi and f_{Bq}\sqrt{B_q} extrapolated to the continuum limit and the physical pion mass. We have used Nf = 2 dynamical quark gauge configurations at four values of the lattice spacing generated by ETMC. Extrapolation in the heavy quark mass from the charm to the bottom quark region has been carried out using ratios of physical quantities computed at nearby quark masses, having an exactly known infinite mass limit., Comment: 7 pages, 4 figures, presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present a lattice QCD computation of the b-quark mass, the B and B_s decay constants, the B-mixing bag parameters for the full four-fermion operator basis as well as determinations for \xi and f_{Bq}\sqrt{B_i^{(q)}} extrapolated to the continuum limit and to the physical pion mass. We used N_f = 2 twisted mass Wilson fermions at four values of the lattice spacing with pion masses ranging from 280 to 500 MeV. Extrapolation in the heavy quark mass from the charm to the bottom quark region has been carried out on ratios of physical quantities computed at nearby quark masses, exploiting the fact that they have an exactly known infinite mass limit. Our results are m_b(m_b, \overline{\rm{MS}})=4.29(12) GeV, f_{Bs}=228(8) MeV, f_{B}=189(8) MeV and f_{Bs}/f_B=1.206(24). Moreover with our results for the bag-parameters we find \xi=1.225(31), B_1^{(s)}/B_1^{(d)}=1.01(2), f_{Bd}\sqrt{\hat{B}_{1}^{(d)}} = 216(10) MeV and f_{Bs}\sqrt{\hat{B}_{1}^{(s)}} = 262(10) MeV. We also computed the bag parameters for the complete basis of the four-fermion operators which are required in beyond the SM theories. By using these results for the bag parameters we are able to provide a refined Unitarity Triangle analysis in the presence of New Physics, improving the bounds coming from B_{(s)}-\bar B_{(s)} mixing., Comment: version accepted for publication in JHEP, extended discussion on error budget, lattice results unchanged, updated upper bound for C1 in Tables 5 and 6, references added; 55 pages, 25 figures, 14 tables

We present a method to evaluate on the lattice the leading isospin breaking effects due to both the small mass difference between the up and down quarks and the QED interaction. Our proposal is applicable in principle to any QCD+QED gauge invariant hadronic observable which can be computed on the lattice. It is based on the expansion of the path-integral in powers of the small parameters (m_d - m_u)/Lambda_{QCD} and alpha_{em}, where m_f is the renormalized quark mass and alpha_{em} the renormalized fine structure constant. In this paper we discuss in detail the general strategy of the method and the conventional, although arbitrary, separation of QCD from QED isospin breaking corrections. We obtain results for the pion mass splitting, M_{pi+}^2-M_{pi0}^2= 1.44(13)(16) x 10^3 MeV^2, for the Dashen's theorem breaking parameter epsilon_{gamma}= 0.79(18)(18), for the light quark masses, [m_d - m_u](MSbar,2 GeV)= 2.39(8)(17) MeV, [m_u / m_d](MSbar,2 GeV)= 0.50(2)(3) and for the flavour symmetry breaking parameters R and Q. We also update our previous results for the QCD isospin breaking corrections to the Kl2 decay rate and for the QCD contribution to the neutron-proton mass splitting., Comment: 39 pages, 17 figures

We present a precise lattice QCD determination of the b-quark mass, of the B and Bs decay constants and first results for the B-meson bag parameters. For our computation we employ the so-called ratio method and our results benefit from the use of improved interpolating operators for the B-mesons. QCD calculations are performed with Nf = 2 dynamical light-quarks at four values of the lattice spacing and the results are extrapolated to the continuum limit. The preliminary results are mb(mb) = 4.35(12) GeV for the MSbar b-quark mass, fBs = 234(6) MeV and fB = 197(10) MeV for the B-meson decay constants, BBs(mb) = 0.90(5) and BB(mb) = 0.87(5) for the B-meson bag parameters., Comment: 6 pages, 3 figures. Proceedings of the 36th International Conference on High Energy Physics - ICHEP 2012; July 4-11 2012; Melbourne, Australia

We present a precise lattice QCD determination of the b-quark mass, of the B and Bs decay constants and first preliminary results for the B-mesons bag parameter. Simulations are performed with Nf = 2 Wilson twisted mass fermions at four values of the lattice spacing and the results are extrapolated to the continuum limit. Our calculation benefits from the use of improved interpolating operators for the B-mesons and employs the so-called ratio method. The latter allows a controlled interpolation at the b-quark mass between the relativistic data around and above the charm quark mass and the exactly known static limit., Comment: 7 pages, 4 figures, 1 table. Proceedings of the 30th International Symposium on Lattice Field Theory - Lattice 2012; June 24-29, 2012; Cairns, Australia