• 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.
  • 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.
  • 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)$.
  • 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.
  • 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.
  • 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.
  • 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).
  • 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.
  • 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).
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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.
  • 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].
  • 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.
  • 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)$.
  • 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 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)$.
  • 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.