Source author record

R. Frezzotti

R. Frezzotti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

53works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

53 published item(s)

preprint2022arXiv

Lattice calculation of the pion mass difference $M_{π^{+}}-M_{π^{0}}$ at order $\mathcal{O}(α_{em})$

We present a lattice calculation of the charged/neutral pion mass difference $M_{π^{+}}-M_{π^{0}}$ at order $\mathcal{O}(α_{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_π \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 $α_{em}$. Making use of the recently proposed RTM scheme, we evaluate the full $\mathcal{O}(α_{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_{π^{+}}-M_{π^{0}}= 4.622~(95)~{\rm MeV}$ which is in good agreement with the experimental result $[ M_{π^{+}} - M_{π^{0}} ]^{exp.} = 4.5936(5)~{\rm MeV}$.

preprint2022arXiv

Virtual Photon Emission in Leptonic Decays of Pseudoscalar Mesons

We present a preliminary non-perturbative lattice calculation of the form factors entering the processes $K\to \ell\,ν_\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.

preprint2020arXiv

First lattice calculation of radiative leptonic decay rates of pseudoscalar mesons

We present a non-perturbative lattice calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar ν_\ell γ$, 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 ν_\ell$ due to the exchange of a virtual photon, this allows accurate predictions at $O(α_{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(α_{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.

preprint2020arXiv

Moments of nucleon generalized parton distributions from lattice QCD simulations at physical pion mass

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.

preprint2016arXiv

Mass of the b-quark and B-decay constants from Nf=2+1+1 twisted-mass Lattice QCD

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

preprint2016arXiv

Nucleon and pion structure with lattice QCD simulations at physical value of the pion mass

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}^{π^\pm}=0.214(15)(^{+12}_{-9})$ in the $\overline{\rm MS}$ at 2 GeV.

preprint2015arXiv

A non-perturbative mechanism for elementary particle mass generation

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.

preprint2015arXiv

Leptonic decay constants fK, fD and fDs with Nf = 2+1+1 twisted-mass lattice QCD

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.

preprint2014arXiv

B-physics from Nf=2 tmQCD: the Standard Model and beyond

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 ξ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 ξ=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.

preprint2014arXiv

D-Dbar Mixing in the Standard Model and Beyond from Nf=2 Twisted Mass QCD

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.

preprint2014arXiv

Heavy flavour precision physics from Nf = 2 + 1 + 1 lattice simulations

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.

preprint2014arXiv

K -> pi matrix elements of the chromagnetic operator on the lattice

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.

preprint2014arXiv

Lattice QCD Study of $B$-meson Decay Constants from ETMC

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.

preprint2014arXiv

Perturbative and non-perturbative renormalization results of the Chromomagnetic Operator on the Lattice

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-π$ matrix element are being performed in simulations with 4 dynamical ($N_f = 2+1+1$) twisted mass fermions and the Iwasaki improved gluon action.

preprint2014arXiv

Progress in Simulations with Twisted Mass Fermions at the Physical Point

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: $μ_s/μ_l=27.63(13)$, $μ_c/μ_l=339.6(2.2)$ and $μ_c/μ_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.

preprint2014arXiv

Up, down, strange and charm quark masses with Nf = 2+1+1 twisted mass lattice QCD

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.

preprint2013arXiv

A $N_f = 2 + 1 + 1$ "twisted" determination of the $b$-quark mass, $f_{B}$ and $f_{B_s}$

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_π) = 1.007 (16)$ and $(f_{B_s}/f_B)/(f_{D_s}/f_D) = 1.008 (13)$.

preprint2013arXiv

A determination of the average up-down, strange and charm quark masses from $N_f=2+1+1$

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

preprint2013arXiv

A first look at maximally twisted mass lattice QCD calculations at the physical point

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.

preprint2013arXiv

B-physics computations from Nf=2 tmQCD

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

preprint2013arXiv

Leading isospin breaking effects on the lattice

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.

preprint2013arXiv

Pseudoscalar decay constants $f_K/f_π$, $f_D$ and $f_{D_s}$ with $N_f = 2 + 1 + 1$ ETMC configurations

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_{π^+}=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_π) = 1.005(15)$.

preprint2013arXiv

The chromomagnetic operator on the lattice

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

preprint2012arXiv

B-physics from lattice QCD...with a twist

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.

preprint2012arXiv

B-physics from the ratio method with Wilson twisted mass fermions

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.

preprint2012arXiv

Isospin breaking effects due to the up-down mass difference in Lattice QCD

We present a new method to evaluate with high precision isospin breaking effects due to the small mass difference between the up and down quarks using lattice QCD. Our proposal is applicable in principle to any hadronic observable which can be computed on the lattice. It is based on the expansion of the path-integral in powers of the small parameter md-mu. In this paper, we apply this method to compute the leading isospin breaking effects for several physical quantities of interest: the kaon meson masses, the kaon decay constant, the form factors of semileptonic Kl3 decays and the neutron-proton mass splitting.

preprint2012arXiv

Lattice QCD calculation of strong isospin breaking effects

We present a new method to evaluate with high precision the isospin breaking effects due to the mass difference between the up and down quarks using lattice QCD. Our proposal is applicable in principle to any hadronic observable which can be computed on the lattice. It is based on the expansion of the path-integral in powers of the small parameter $m_d - m_u$. In this talk we discuss how to apply this method to compute the leading isospin breaking effects for several physical quantities of interest: the kaon masses, the kaon decay constants and the neutron-proton mass splitting.

preprint2012arXiv

Lattice QCD determination of m_b, f_B and f_Bs with twisted mass Wilson fermions

We present a lattice QCD determination of the b quark mass and of the B and B_s decay constants, performed with N_f=2 twisted mass Wilson fermions, by simulating at four values of the lattice spacing. In order to study the b quark on the lattice, two methods are adopted in the present work, respectively based on suitable ratios with exactly known static limit and on the interpolation between relativistic data, evaluated in the charm mass region, and the static point, obtained by simulating the HQET on the lattice. The two methods provide results in good agreement. For the b quark mass in the MSbar scheme and for the decay constants we obtain m_b(m_b)=4.29(14) GeV, f_B=195(12) MeV, f_Bs=232(10) MeV and f_Bs/f_B=1.19(5). As a byproduct of the analysis we also obtain the results for the f_D and f_Ds decay constants: f_D=212(8) MeV, f_Ds=248(6) MeV and f_Ds/f_D=1.17(5).

preprint2012arXiv

Neutral meson oscillations in the Standard Model and beyond from Nf=2 Twisted Mass Lattice QCD

We present the ETMC results for the bag parameters describing the neutral kaon mixing in the Standard Model and beyond and preliminary results for the bag parameters controlling the short distance contributions in the D^0-\bar{D}^0 oscillations. We also present preliminary results for the B_{Bd}, B_{Bs}, B_{Bs}/B_{Bd} and ξ-parameter controlling B^0_-\bar{B}^0 oscillations in the Standard Model employing the so-called ratio method. Using Nf=2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks we achieve both O(a)-improvement and continuum like renormalization pattern. Simulations are performed at three-values of the lattice spacing and several values of quark masses in the light, strange, charm region and above charm up to ~2.5m_c. Our results are extrapolated to the continuum limit and extrapolated/interpolated to the physical quark masses.

preprint2011arXiv

Average up/down, strange and charm quark masses with Nf=2 twisted mass lattice QCD

We present a high precision lattice calculation of the average up/down, strange and charm quark masses performed with Nf=2 twisted mass Wilson fermions. The analysis includes data at four values of the lattice spacing and pion masses as low as ~270 MeV, allowing for accurate continuum limit and chiral extrapolation. The strange and charm masses are extracted by using several methods, based on different observables: the kaon and the eta_s meson for the strange quark and the D, D_s and eta_c mesons for the charm. The quark mass renormalization is carried out non-perturbatively using the RI-MOM method. The results for the quark masses in the MSbar scheme read: m_ud(2 GeV)= 3.6(2) MeV, m_s(2 GeV)=95(6) MeV and m_c(m_c)=1.28(4) GeV. We also obtain the ratios m_s/m_ud=27.3(9) and m_c/m_s=12.0(3).

preprint2011arXiv

BK-parameter from Nf = 2 twisted mass lattice QCD

We present an unquenched $N_f=2$ lattice computation of the $B_{K}$ parameter which controls $K^0-\bar K^0$ oscillations. A partially quenched setup is employed with two maximally twisted dynamical (sea) light Wilson quarks, and valence quarks of both the maximally twisted and the Osterwalder--Seiler variety. Suitable combinations of these two kinds of valence quarks lead to a lattice definition of the $B_{K}$ parameter which is both multiplicatively renormalizable and O($a$) improved. Employing the non-perturbative RI-MOM scheme, in the continuum limit and at the physical value of the pion mass we get $B^{\rm RGI}_K=0.729\pm 0.030$, a number well in line with the existing quenched and unquenched determinations.

preprint2011arXiv

K^0-\bar{K}^0 mixing in the Standard Model from Nf=2+1+1 Twisted Mass Lattice QCD

We present preliminary results at β = 1.95 (a = 0.077 fm) on the first unquenched N_f=2+1+1 lattice computation of the B_K parameter which controls the neutral kaon oscillations in the Standard Model. Using N_f=2+1+1 maximally twisted sea quarks and Osterwalder-Seiler valence quarks we achieve O(a) improvement and a continuum-like renormalization pattern for the four-fermion operator. Our results are extrapolated/interpolated to the physical light/strange quark mass but not yet to the continuum limit. The computation of the relevant renormalization constants is performed non perturbatively in the RI'-MOM scheme using dedicated simulations with N_f=4 degenerate sea quark flavours produced by the ETM collaboration. We get B_K^{RGI} (a = 0.077) = 0.747(18), which when compared to our previous unquenched N_f=2 determination and most of the existing results, suggests a rather weak B_K^{RGI} dependence on the number of dynamical flavours. We are at the moment analysing lattice data at two additional β values which will allow us to perform an extrapolation to the continuum limit.

preprint2011arXiv

Renormalisation constants of quark bilinears in lattice QCD with four dynamical Wilson quarks

We present preliminary results of the non-perturbative computation of the RI-MOM renormalisation constants in a mass-independent scheme for the action with Iwasaki glue and four dynamical Wilson quarks employed by ETMC. Our project requires dedicated gauge ensembles with four degenerate sea quark flavours at three lattice spacings and at several values of the standard and twisted quark mass parameters. The RI-MOM renormalisation constants are obtained from appropriate O(a) improved estimators extrapolated to the chiral limit.

preprint2011arXiv

Renormalization constants for Wilson fermion lattice QCD with four dynamical flavours

We report on an ongoing non-perturbative computation of RI-MOM scheme renormalization constants for the lattice action with four dynamical flavours currently in use by ETMC. For this goal dedicated simulations with four degenerate sea quark flavours are performed at several values of the standard and twisted quark mass parameters. We discuss a method for removing possible O(a) artifacts at all momenta and extrapolating renormalization constant estimators to the chiral limit. We give preliminary results at one lattice spacing.

preprint2010arXiv

A proposal for B-physics on current lattices

A method to extract B-physics parameters (b-quark mass and fB, fBs decay constants) from currently available lattice data is presented and tested. The approach is based on the idea of constructing appropriate ratios of heavy-light meson masses and decay constants, respectively, possessing a precisely known static limit, and evaluating them at various pairs of heavy quark masses around the charm. Via a smooth interpolation in the heavy quark mass from the easily accessible charm region to the asymptotic point, B-physics parameters are computed with a few percent (statistical + systematic) error using recently produced Nf=2 maximally twisted Wilson fermions data.

preprint2010arXiv

K^0-\bar{K}^0 Mixing Beyond the SM from Nf=2 tmQCD

We present preliminary results on the of neutral kaon oscillations in extensions of the Standard Model. Using Nf=2 maximally twisted sea quarks and Osterwalder-Seiler valence quarks, we achieve both O(a)-improvement and continuum-like renormalization pattern for the relevant four-fermion operators. We perform simulations at three values of the lattice spacing and extrapolate/interpolate our results to the continuum limit and physical light/strange quark mass. The calculation of the renormalization constants of the complete operator basis is performed non- perturbatively in the RI-MOM scheme.

preprint2010arXiv

Non-perturbative renormalization of quark bilinear operators with Nf=2 (tmQCD) Wilson fermions and the tree-level improved gauge action

We present results for the renormalization constants of bilinear quark operators obtained by using the tree-level Symanzik improved gauge action and the Nf=2 twisted mass fermion action at maximal twist, which guarantees automatic O(a)-improvement. Our results are also relevant for the corresponding standard (un-twisted) Wilson fermionic action since the two actions only differ, in the massless limit, by a chiral rotation of the quark fields. The scale-independent renormalization constants ZV, ZA and the ratio ZP/ZS have been computed using the RI-MOM approach, as well as other alternative methods. For ZA and ZP/ZS, the latter are based on both standard twisted mass and Osterwalder-Seiler fermions, while for ZV a Ward Identity has been used. The quark field renormalization constant Zq and the scale dependent renormalization constants ZS, ZP and ZT are determined in the RI-MOM scheme. Leading discretization effects of O(g^2 a^2), evaluated in one-loop perturbation theory, are explicitly subtracted from the RI-MOM estimates.

preprint2010arXiv

Quark masses with Nf=2 twisted mass lattice QCD

We present the results of the recent high precision lattice calculation of the average up/down, strange and charm quark masses performed by ETMC with Nf=2 twisted mass Wilson fermions. The analysis includes data at four values of the lattice spacing and pion masses as low as ~270 MeV, allowing for accurate continuum limit and chiral extrapolation. The strange and charm masses are extracted by using several methods, based on different observables: the kaon and the eta_s meson for the strange quark and the D, D_s and eta_c mesons for the charm. The quark mass renormalization is carried out non-perturbatively using the RI-MOM method. The results for the quark masses in the MSbar scheme read: m_ud(2 GeV)= 3.6(2) MeV, m_s(2 GeV)=95(6) MeV and m_c(m_c)=1.28(4) GeV. We have also obtained the ratios m_s/m_ud=27.3(9) and m_c/m_s=12.0(3). Moreover, we provide the updated result for the bottom quark mass, m_b(m_b)=4.3(2) GeV, obtained using the method presented in 0909.3187 [hep-lat].

preprint2009arXiv

f_B and f_Bs with maximally twisted Wilson fermions

We present a lattice QCD calculation of the heavy-light decay constants f_B and f_Bs performed with Nf=2 maximally twisted Wilson fermions, at four values of the lattice spacing. The decay constants have been also computed in the static limit and the results are used to interpolate the observables between the charm and the infinite-mass sectors, thus obtaining the value of the decay constants at the physical b quark mass. Our preliminary results are f_B=191(14) MeV, f_Bs=243(14) MeV, f_Bs/f_B=1.27(5). They are in good agreement with those obtained with a novel approach, recently proposed by our Collaboration (ETMC), based on the use of suitable ratios having an exactly known static limit.

preprint2009arXiv

Kaon oscillations in the Standard Model and Beyond using Nf=2 dynamical quarks

We compute non-perturbatively the B-parameters of the complete basis of four-fermion operators needed to study the Kaon oscillations in the SM and in its supersymmetric extension. We perform numerical simulations with two dynamical maximally twisted sea quarks at three values of the lattice spacing on configurations generated by the ETMC. Unwanted operator mixings and O(a) discretization effects are removed by discretizing the valence quarks with a suitable Osterwalder-Seiler variant of the Twisted Mass action. Operators are renormalized non-perturbatively in the RI/MOM scheme. Our preliminary result for BK(RGI) is 0.73(3)(3).

preprint2009arXiv

Light Meson Physics from Maximally Twisted Mass Lattice QCD

We present a comprehensive investigation of light meson physics using maximally twisted mass fermions for two mass-degenerate quark flavours. By employing four values of the lattice spacing, spatial lattice extents ranging from 2.0 fm to 2.5 fm and pseudo scalar masses in the range 280 MeV to 650 MeV we control the major systematic effects of our calculation. This enables us to confront our data with chiral perturbation theory and extract low energy constants of the effective chiral Lagrangian and derived quantities, such as the light quark mass, with high precision.

preprint2009arXiv

O(a^2) cutoff effects in lattice Wilson fermion simulations

In this paper we propose to interpret the large discretization artifacts affecting the neutral pion mass in maximally twisted lattice QCD simulations as O(a^2) effects whose magnitude is roughly proportional to the modulus square of the (continuum) matrix element of the pseudoscalar density operator between vacuum and one-pion state. The numerical size of this quantity is determined by the dynamical mechanism of spontaneous chiral symmetry breaking and turns out to be substantially larger than its natural magnitude set by the value of Lambda_QCD.

preprint2009arXiv

Pseudoscalar decay constants of kaon and D-mesons from Nf=2 twisted mass Lattice QCD

We present the results of a lattice QCD calculation of the pseudoscalar meson decay constants fpi, fK, fD and fDs, performed with Nf=2 dynamical fermions. The simulation is carried out with the tree-level improved Symanzik gauge action and with the twisted mass fermionic action at maximal twist. We have considered for the final analysis three values of the lattice spacing, a~0.10 fm, 0.09 fm and 0.07 fm, with pion masses down to mpi~270 MeV. Our results for the light meson decay constants are fK=158.1(2.4) MeV and fK/fpi=1.210(18). From the latter ratio, by using the experimental determination of Gamma(K-->mu nu_mu (gamma))/ Gamma(pi--> mu nu_mu (gamma)) and the average value of |Vud| from nuclear beta decays, we obtain |Vus|=0.2222(34), in good agreement with the determination from semileptonic Kl3 decays and the unitarity constraint. For the D and Ds meson decay constants we obtain fD=197(9) MeV, fDs=244(8) MeV and fDs/fD=1.24(3). Our result for fD is in good agreement with the CLEO experimental measurement. For fDs our determination is smaller than the PDG 2008 experimental average but in agreement with a recent improved measurement by CLEO at the 1.4 sigma level.

preprint2008arXiv

Dynamical Twisted Mass Fermions with Light Quarks: Simulation and Analysis Details

In a recent paper [hep-lat/0701012] we presented precise lattice QCD results of our European Twisted Mass Collaboration (ETMC). They were obtained by employing two mass-degenerate flavours of twisted mass fermions at maximal twist. In the present paper we give details on our simulations and the computation of physical observables. In particular, we discuss the problem of tuning to maximal twist, the techniques we have used to compute correlators and error estimates. In addition, we provide more information on the algorithm used, the autocorrelation times and scale determination, the evaluation of disconnected contributions and the description of our data by means of chiral perturbation theory formulae.

preprint2007arXiv

Dynamical Twisted Mass Fermions with Light Quarks

We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.

preprint2007arXiv

Light quark masses and pseudoscalar decay constants from Nf=2 Lattice QCD with twisted mass fermions

We present the results of a lattice QCD calculation of the average up-down and strange quark masses and of the light meson pseudoscalar decay constants with Nf=2 dynamical fermions. The simulation is carried out at a single value of the lattice spacing with the twisted mass fermionic action at maximal twist, which guarantees automatic O(a)-improvement of the physical quantities. Quark masses are renormalized by implementing the non-perturbative RI-MOM renormalization procedure. Our results for the light quark masses are m_ud^{msbar}(2 GeV)= 3.85 +- 0.12 +- 0.40 MeV, m_s^{msbar}(2 GeV) = 105 +- 3 +- 9 MeV and m_s/m_ud = 27.3 +- 0.3 +- 1.2. We also obtain fK = 161.7 +- 1.2 +- 3.1 MeV and the ratio fK/fpi=1.227 +- 0.009 +- 0.024. From this ratio, by using the experimental determination of Gamma(K-> mu nu (gamma))/Gamma(pi -> mu nu (gamma)) and the average value of |Vud| from nuclear beta decays, we obtain |Vus|=0.2192(5)(45), in agreement with the determination from Kl3 decays and the unitarity constraint.

preprint2004arXiv

Chirally improving Wilson fermions II. Four-quark operators

In this paper we discuss how the peculiar properties of twisted lattice QCD at maximal twist can be employed to set up a consistent computational scheme in which, despite the explicit breaking of chiral symmetry induced by the presence of the Wilson and mass terms in the action, it is possible to completely bypass the problem of wrong chirality and parity mixings in the computation of the CP-conserving matrix elements of the $ΔS=1,2$ effective weak Hamiltonian and at the same time have a positive determinant for non-degenerate quarks as well as full O($a$) improvement in on-shell quantities with no need of improving the lattice action and the operators.

preprint2001arXiv

Cutoff effects in twisted mass lattice QCD

We present a first numerical study of lattice QCD with O(a) improved Wilson quarks and a chirally twisted mass term. Renormalized correlation functions are derived from the Schroedinger functional and evaluated in an intermediate space-time volume of size 0.75^3 x 1.5 fm^4. In the quenched approximation precise results are then obtained with a moderate computational effort, allowing for a detailed study of the continuum approach. The latter is discussed in terms of observables which converge to meson masses and decay constants in the limit of large space-time volume. In the O(a) improved theory we find residual cutoff effects to be at the level of a few percent for lattice spacings of about 0.1 fm.

preprint2001arXiv

Quenched twisted mass QCD at small quark masses and in large volume

As a test of quenched lattice twisted mass QCD, we compute the non-perturbatively O($a$) improved pseudoscalar and vector meson masses and the pseudoscalar decay constant down to $M_{\rm PS}/M_{\rm V} = 0.467(13)$ at $β=6$ in large volume. We check the absence of exceptional configurations and -- by further data at $β=6.2$ -- the size of scaling violations. The CPU time cost for reaching a given accuracy is close to that with ordinary Wilson quarks at $M_{\rm PS}/M_{\rm V} \simeq 0.6$ and grows smoothly as $M_{\rm PS}/M_{\rm V}$ decreases.

preprint2001arXiv

Status of APEmille

This paper presents the status of the APEmille project, which is essentially completed, as far as machine development and construction is concerned. Several large installations of APEmille are in use for physics production runs leading to many new results presented at this conference. This paper briefly summarizes the APEmille architecture, reviews the status of the installations and presents some performance figures for physics codes.