Source author record

V. Lubicz

V. Lubicz 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

96works
3topics
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

96 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

Power corrections to the CP-violation parameter $\varepsilon_K$

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

preprint2022arXiv

Virtual Photon Emission in Leptonic Decays of Charged Pseudoscalar Mesons

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

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.

preprint2016arXiv

Electromagnetic corrections to the leptonic decay rates of charged pseudoscalar mesons: lattice results

Electromagnetic effects in the leptonic decay rates $π^+ \to μ^+ ν$ and $K^+ \to μ^+ ν$ are evaluated for the first time on the lattice. Following a method recently proposed in Ref. [1] the emission of virtual photons at leading order in the electromagnetic coupling is evaluated on the lattice and the infrared divergence computed for a point-like meson at finite lattice volume is subtracted. The physical decay rate is then obtained by adding the emission of real and virtual photons regularised with a photon mass. Using the gauge ensembles produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ dynamical quarks the feasibility of our approach is demonstrated. Preliminary results for the electromagnetic corrections to charged (neutral) pion and kaon masses as well as to the leptonic decay rates of charged pions and kaons are presented.

preprint2016arXiv

Hypercubic Effects in semileptonic $D \to π$ decays on the lattice

We present a lattice determination of the vector and scalar form factors of the semileptonic $D \to π\ell ν$ decays, which are relevant for the extraction of the CKM matrix element $\lvert V_{cd} \rvert$ from experimental data. Our analysis is based on the gauge configurations produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ flavors of dynamical quarks. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV. Quark momenta are injected on the lattice using non-periodic boundary conditions. The matrix elements of both vector and scalar currents are determined for a plenty of kinematical conditions in which parent and child mesons are either moving or at rest. Lorentz symmetry breaking due to hypercubic effects is clearly observed in the data and included in the decomposition of the current matrix elements in terms of additional form factors. Our preliminary estimate for the vector form factor at zero 4-momentum transfer is $f_+^{D \to π}(0) = 0.631\,(40)$, which can be compared with the latest FLAG average $f_+^{D \to π}(0) = 0.666\,(29)$ available only at $N_f = 2 + 1$.

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

Review of lattice results concerning low-energy particle physics

We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.

preprint2015arXiv

D -> πl νand D -> K l νform factors with Nf=2+1+1 Twisted fermions

We present a lattice QCD determination of the vector and scalar form factors of the semileptonic decays D->πl νand D -> K l νwhich are relevant for the extraction of the CKM matrix elements |Vcd| and |Vcs| from experimental data. Our analysis is based on the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2+1+1 dynamical fermions. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV. Our preliminary estimates for the vector form factor at zero 4-momentum transfer are f+(D -> π)(0) = 0.610 (23) and f+(D -> K)(0) = 0.747 (22), where the uncertainties are only statistical. By combining our results with the experimental values of f+(D -> π)(0) |Vcd| and f+(D -> K)(0) |Vcs| we obtain |Vcd| = 0.2336 (93) and |Vcs| = 0.975 (30), which together with the PDG determination of |Vcb| are in agreement with the unitarity constraint of the Standard Model.

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.

preprint2015arXiv

Momentum dependence of kaon semileptonic form factors with Nf=2+1+1 Twisted Mass fermions

We present a lattice QCD determination of the vector and scalar form factors of the kaon semileptonic decay $K \to π\ell ν$, which is relevant for the determination of the CKM matrix element $|V_{us}|$ from experimental data. Our results are based on the gauge configurations produced by the European Twisted Mass Collaboration with Nf = 2+1+1 dynamical fermions. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV. Our estimate for the vector form factor at zero 4-momentum transfer is $f_+(0) = 0.9709 (46)$, where the uncertainty is both statistical and systematic. By combining our result with the latest experimental value of $f_+(0)|V_{us}|$ we obtain $|V_{us}| = 0.2230 (11)$, which satisfies the unitarity constraint of the Standard Model at the permille level using the updated determination of $|V_{ud}|$ coming from superallowed nuclear $β$ decays. We present also the momentum dependence of the vector and scalar form factors in the whole range of values of the squared 4-momentum transfer measured in $K_{\ell 3}$ decays, obtaining a good agreement with the experimental data.

preprint2015arXiv

QED Corrections to Hadronic Processes in Lattice QCD

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

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

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

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

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

Matrix elements of the electromagnetic operator between kaon and pion states

We compute the matrix elements of the electromagnetic (EM) operator between kaon and pion states, using lattice QCD with maximally twisted-mass fermions and two flavors of dynamical quarks (Nf = 2). The EM operator is renormalized non-perturbatively in the RI'/MOM scheme and our simulations cover pion masses as light as 260 MeV and three values of the lattice spacing, ranging from ~ 0.07 up to ~ 0.1 fm. At the physical point our preliminary result for the K -> pion tensor form factor at zero-momentum transfer is fT[Kπ](0) = 0.42(2_stat), which differs significantly from the old quenched result fT[Kπ](0) = 0.78(6) obtained by the SPQcdR Collaboration with pion masses above 500 MeV. We investigate the source of this difference and conclude that it is mainly related to the chiral extrapolation of the quenched data. For the case of the tensor charge of the pion we obtain the preliminary value fT[ππ](0) = 0.200(14_stat), which can be compared with the result fT[ππ](0) = 0.216(34) obtained at Nf = 2 by the QCDSF Collaboration using higher pion masses.

preprint2011arXiv

Matrix elements of the electromagnetic operator between kaon and pion states

We compute the matrix elements of the electromagnetic (EM) operator (\bar{s} F{μν} σ{μν} d) between kaon and pion states, using lattice QCD with maximally twisted-mass fermions and two flavors of dynamical quarks (Nf = 2). The operator is renormalized non-perturbatively in the RI'/MOM scheme and our simulations cover pion masses as light as 270 MeV and three values of the lattice spacing from ~ 0.07 up to ~ 0.1 fm. At the physical point our result for the corresponding tensor form factor at zero-momentum transfer is fT{Kπ}(0) = 0.417 (14_stat) (5_syst), where the systematic error does not include the effect of quenching the strange and charm quarks. Our result differs significantly from the old quenched result fT{Kπ}(0) = 0.78 (6) obtained by the SPQcdR Collaboration with pion masses above 500 MeV. We investigate the source of this difference and conclude that it is mainly related to the chiral extrapolation. We also study the tensor charge of the pion and obtain the value fT{ππ}(0) = 0.195 (8_stat) (6_syst) in good agreement with, but more accurate than the result fT{ππ}(0) = 0.216 (34) obtained by the QCDSF Collaboration using higher pion masses.

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

An Improved Standard Model Prediction Of BR(B -> tau nu) And Its Implications For New Physics

The recently measured B -> tau nu branching ratio allows to test the Standard Model by probing virtual effects of new heavy particles, such as a charged Higgs boson. The accuracy of the test is currently limited by the experimental error on BR(B -> tau nu) and by the uncertainty on the parameters fB and |Vub|. The redundancy of the Unitarity Triangle fit allows to reduce the error on these parameters and thus to perform a more precise test of the Standard Model. Using the current experimental inputs, we obtain BR(B -> tau nu)_SM = (0.84 +- 0.11)x10^{-4}, to be compared with BR(B -> tau nu)_exp = (1.73 +- 0.34)x10^{-4}. The Standard Model prediction can be modified by New Physics effects in the decay amplitude as well as in the Unitarity Triangle fit. We discuss how to disentangle the two possible contributions in the case of minimal flavour violation at large tan beta and generic loop-mediated New Physics. We also consider two specific models with minimal flavour violation: the Type-II Two Higgs Doublet Model and the Minimal Supersymmetric Standard Model.

preprint2010arXiv

Flavor Physics in the Quark Sector

One of the major challenges of particle physics has been to gain an in-depth understanding of the role of quark flavor and measurements and theoretical interpretations of their results have advanced tremendously: apart from masses and quantum numbers of flavor particles, there now exist detailed measurements of the characteristics of their interactions allowing stringent tests of Standard Model predictions. Among the most interesting phenomena of flavor physics is the violation of the CP symmetry that has been subtle and difficult to explore. Till early 1990s observations of CP violation were confined to neutral $K$ mesons, but since then a large number of CP-violating processes have been studied in detail in neutral $B$ mesons. In parallel, measurements of the couplings of the heavy quarks and the dynamics for their decays in large samples of $K, D$, and $B$ mesons have been greatly improved in accuracy and the results are being used as probes in the search for deviations from the Standard Model. In the near future, there will be a transition from the current to a new generation of experiments, thus a review of the status of quark flavor physics is timely. This report summarizes the results of the current generation of experiments that is about to be completed and it confronts these results with the theoretical understanding of the field.

preprint2010arXiv

Improved analysis of the scalar and vector form factors of kaon semileptonic decays with N_f = 2 twisted-mass fermions

We investigate the vector and scalar form factors relevant for K_{\ell 3} semileptonic decays using maximally twisted-mass fermions with two flavors of dynamical quarks (N_f = 2). The simulations cover pion masses as light as 260 MeV and four values of the lattice spacing, ranging from ~0.05 up to ~0.1 fm, which allow to compute directly, for the first time, the continuum limit for the vector form factor at zero-momentum transfer, f_+(0). The preliminary result is f_+(0) = 0.9544(68), where the error is statistical only. We also extrapolate both form factors to the physical point and study their momentum dependence. Our results are in good agreement with those obtained from a dispersion analyses of the experimental data. Together with the form factors, we analyze the ratio of the leptonic decay constants f_K / f_π, by imposing the constraint coming from the Callan-Treiman theorem, obtaining at the physical point f_K / f_π= 1.190(8). Combining our results for f_+(0) and f_K / f_πwith the experimental measurements of the leptonic and semilpetonic decay rates, and using the determination of |V_{ud}| from nuclear beta decays, we determine the values of the Cabibbo angle |V_{us}| from both K_{\ell 3} and K_{\ell 2} decays, obtaining |V_{us}|^{K_{\ell 3}} = 0.2266(17) and $|V_{us}|^{K_{\ell 2}} = 0.2258(16).

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

Pseudoscalar meson decay constants f_K, f_D and f_Ds, from Nf=2 twisted mass Lattice QCD

We present the results of a lattice QCD calculation of the pseudoscalar meson decay constants f_K, f_D and f_Ds, performed with N_f=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. With respect to our previous study (0709.4574 [hep-lat]), here we have analysed data at three values of the lattice spacing (a=0.10 fm, 0.09 fm, 0.07 fm) and performed the continuum limit, and we have included at a=0.09 fm data with a lighter quark mass (m_pi = 260 MeV) and a larger volume (L = 2.7 fm), thus having at each lattice spacing L >= 2.4 fm and m_pi*L >= 3.6. Our result for the kaon decay constant is f_K=(157.5 +- 0.8|_{stat.} +- 3.3|_{syst.}) MeV and for the ratio f_K/f_pi=1.205 +- 0.006|_{stat.} +- 0.025|_{syst.}, in good agreement with the other N_f=2 and N_f=2+1 lattice calculations. For the D and D_s meson decay constants we obtain f_D=(205 +- 7|_{stat.} +- 7|_{syst.}) MeV, in good agreement with the CLEO-c experimental measurement and with other recent N_f=2 and N_f=2+1 lattice calculations, and f_{Ds}=(248 +- 3|_{stat.} +- 8|_{syst.}) MeV that, instead, is 2.3 sigma below the CLEO-c/BABAR experimental average, confirming the present tension between lattice calculations and experimental measurements.

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

preprint2010arXiv

Update of the Unitarity Triangle Analysis

We present the status of the Unitarity Triangle Analysis (UTA), within the Standard Model (SM) and beyond, with experimental and theoretical inputs updated for the ICHEP 2010 conference. Within the SM, we find that the general consistency among all the constraints leaves space only to some tension (between the UTA prediction and the experimental measurement) in BR(B -> tau nu), sin(2 beta) and epsilon_K. In the UTA beyond the SM, we allow for New Physics (NP) effects in (Delta F)=2 processes. The hint of NP at the 2.9 sigma level in the B_s-\bar B_s mixing turns out to be confirmed by the present update, which includes the new D0 result on the dimuon charge asymmetry but not the new CDF measurement of phi_s, being the likelihood not yet released.

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

First Evidence of New Physics in b <--> s Transitions

We combine all the available experimental information on Bs mixing, including the very recent tagged analyses of Bs to J/Psi phi by the CDF and D0 collaborations. We find that the phase of the Bs mixing amplitude deviates more than 3 sigma from the Standard Model prediction. While no single measurement has a 3 sigma significance yet, all the constraints show a remarkable agreement with the combined result. This is a first evidence of physics beyond the Standard Model. This result disfavours New Physics models with Minimal Flavour Violation with the same significance.

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) corrections to the one-loop propagator and bilinears of clover fermions with Symanzik improved gluons

We calculate corrections to the fermion propagator and to the Green's functions of all fermion bilinear operators of the form $\barΨΓΨ$, to one-loop in perturbation theory. We employ the Wilson/clover action for fermions and a family of Symanzik improved actions for gluons. The novel aspect of our calculations is that they are carried out to second order in the lattice spacing, $O(a^2)$. Consequently, they have addressed a number of new issues, most notably the appearance of loop integrands with strong IR divergences (convergent only beyond 6 dimensions). Such integrands are not present in $O(a^1)$ improvement calculations; there, IR divergent terms are seen to have the same structure as in the $O(a^0)$ case, by virtue of parity under integration, and they can thus be handled by well-known techniques. We explain how to correctly extract the full $O(a^2)$ dependence; in fact, our method is generalizable to any order in $a$. The $O(a^2)$ corrections to the quark propagator and Green's functions computed in this paper are useful to improve the nonperturbative RI-MOM determination of renormalization constants for quark bilinear operators. Our results depend on a large number of parameters: coupling constant, number of colors, lattice spacing, external momentum, clover parameter, Symanzik coefficients, gauge parameter. To make these results most easily accessible to the reader, we have included them in the distribution package of this paper, as an ASCII file named: Oa2results.m ; the file is best perused as Mathematica input.

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.

preprint2009arXiv

Status of the Unitarity Triangle Analysis

We present an update of the Unitarity Triangle (UT) analysis, within the Standard Model (SM) and beyond. Within the SM the main novelties are the inclusion in epsilon_K of the contributions of xi and phi_epsilon \neq π/4 pointed out by A.J. Buras and D. Guadagnoli, and an accurate prediction of BR(B -> tau nu), by using the indirect determination of |V_ub| from the UT fit, which can be compared to the present experimental result. In the generalization of the UT analysis to investigate New Physics (NP) effects, the estimate of xi is more delicate and only the effect of phi_epsilon \neq π/4 has been included. We confirm an hint of NP in the B_s-\bar B_s mixing at the 2.9 sigma level, which makes a comparison with new experimental data certainly desired.

preprint2009arXiv

Vector and scalar form factors for K- and D-meson semileptonic decays from twisted mass fermions with Nf = 2

We present lattice results for the form factors relevant in the K -> pion and D -> pion semileptonic decays, obtained from simulations with two flavors of dynamical twisted-mass fermions and pion masses as light as 260 MeV. For K -> pion decays we discuss the estimates of the main sources of systematic uncertainties, including the quenching of the strange quark, leading to our final result f+(0) = 0.9560 (57) (62). Combined with the latest experimental data, our value of f+(0) implies for the CKM matrix element |Vus| the value 0.2267 (5) (20) consistent with the first-row CKM unitarity. For D -> pion decays the application of Heavy Meson Chiral Perturbation Theory allows to extrapolate our results for both the scalar and the vector form factors at the physical point with quite good accuracy, obtaining a nice agreement with the experimental data. In particular at zero-momentum transfer we obtain f+(0) = 0.64 (5).

preprint2008arXiv

B, D and K decays

With the advent of the LHC, we will be able to probe New Physics (NP) up to energy scales almost one order of magnitude larger than it has been possible with present accelerator facilities. While direct detection of new particles will be the main avenue to establish the presence of NP at the LHC, indirect searches will provide precious complementary information, since most probably it will not be possible to measure the full spectrum of new particles and their couplings through direct production. In particular, precision measurements and computations in the realm of flavour physics are expected to play a key role in constraining the unknown parameters of the Lagrangian of any NP model emerging from direct searches at the LHC. The aim of Working Group 2 was twofold: on one hand, to provide a coherent, up-to-date picture of the status of flavour physics before the start of the LHC; on the other hand, to initiate activities on the path towards integrating information on NP from high-pT and flavour data.

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.

preprint2008arXiv

Model-independent constraints on Delta F=2 operators and the scale of New Physics

We update the constraints on new-physics contributions to Delta F=2 processes from the generalized unitarity triangle analysis, including the most recent experimental developments. Based on these constraints, we derive upper bounds on the coefficients of the most general Delta F=2 effective Hamiltonian. These upper bounds can be translated into lower bounds on the scale of new physics that contributes to these low-energy effective interactions. We point out that, due to the enhancement in the renormalization group evolution and in the matrix elements, the coefficients of non-standard operators are much more constrained than the coefficient of the operator present in the Standard Model. Therefore, the scale of new physics in models that generate new Delta F=2 operators, such as next-to-minimal flavour violation, has to be much higher than the scale of minimal flavour violation, and it most probably lies beyond the reach of direct searches at the LHC.

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.

preprint2006arXiv

Exploring twisted mass Lattice QCD with the Clover term

It has been shown that in the twisted mass formulation of Lattice QCD at maximal twist large cutoff effects are generated when the quark mass becomes of O(a Lambda_QCD^2). In general, these effects can be suppressed in two ways: either by choosing the critical quark mass in an "optimal way", or by adding the Clover term to the twisted action. We investigate the second option by performing a quenched lattice QCD simulation with twisted Clover fermions and pion masses as low as 280 MeV. We show that the Clover term is indeed efficient in reducing the large cutoff effects. In particular, the so-called bending phenomenon observed in the determination of the pion decay constant is cured in this way. In addition, by using maximally twisted Clover fermions, we provide a non-perturbative determination of the vector current renormalization constant Z_V as well as of the non-perturbatively renormalized light quark masses. Finally, we calculate the connected contribution to the charged-neutral pseudoscalar meson mass splitting, finding that the introduction of the Clover term in the twisted action is also beneficial, in the quenched approximation, in reducing cutoff effects related to the isospin symmetry breaking at finite lattice spacing.

preprint2006arXiv

Next-to-leading order strong interaction corrections to the Delta F=2 effective Hamiltonian in the MSSM

We compute the next-to-leading order strong interaction corrections to gluino-mediated Delta F=2 box diagrams in the Minimal Supersymmetric Standard Model. These corrections are given by two loop diagrams which we have calculated in three different regularization schemes in the mass insertion approximation. We obtain the next-to-leading order Wilson coefficients of the Delta F=2 effective Hamiltonian relevant for neutral meson mixings. We find that the matching scale uncertainty is largely reduced at the next-to-leading order, typically from about 10-15% to few percent.

preprint2006arXiv

The Unitarity Triangle Fit in the Standard Model and Hadronic Parameters from Lattice QCD: A Reappraisal after the Measurements of Delta m_s and BR(B to tau nu)

The recent measurements of the B_s meson mixing amplitude by CDF and of the leptonic branching fraction BR(B to tau nu) by Belle call for an upgraded analysis of the Unitarity Triangle in the Standard Model. Besides improving the previous constraints on the parameters of the CKM matrix, these new measurements, combined with the recent determinations of the angles alpha, beta and gamma from non-leptonic decays, allow, in the Standard Model, a quite accurate extraction of the values of the hadronic matrix elements relevant for K-Kbar and B_{s,d}-B_{s,d}bar mixing and of the leptonic decay constant f_B. These values, obtained ``experimentally'', can then be compared with the theoretical predictions, mainly from lattice QCD. In this paper we upgrade the UT fit, we determine from the data the kaon B-parameter B_Khat, the B^0 mixing amplitude parameters f_Bs B^{1/2}_Bs and xi, the decay constant f_B, and make a comparison of the obtained values with lattice predictions. We also discuss the different determinations of V_ub and show that current data do not favour the value measured in inclusive decays.

preprint2006arXiv

The UTfit Collaboration Report on the Status of the Unitarity Triangle beyond the Standard Model I. Model-independent Analysis and Minimal Flavour Violation

Starting from a (new physics independent) tree level determination of rhobar and etabar, we perform the Unitarity Triangle analysis in general extensions of the Standard Model with arbitrary new physics contributions to loop-mediated processes. Using a simple parameterization, we determine the allowed ranges of non-standard contributions to |Delta F|=2 processes. Remarkably, the recent measurements from B factories allow us to determine with good precision the shape of the Unitarity Triangle even in the presence of new physics, and to derive stringent constraints on non-standard contributions to |Delta F|=2 processes. Since the present experimental constraints favour models with Minimal Flavour Violation, we present the determination of the Universal Unitarity Triangle that can be defined in this class of extensions of the Standard Model. Finally, we perform a combined fit of the Unitarity Triangle and of new physics contributions in Minimal Flavour Violation, reaching a sensitivity to a new physics scale of about 5 TeV. We also extrapolate all these analyses into a "year 2010" scenario for experimental and theoretical inputs in the flavour sector. All the results presented in this paper are also available at the URL http://www.utfit.org, where they are continuously updated.

preprint2005arXiv

Chiral extrapolation of hyperon vector form factors

We present a new study of SU(3)-breaking corrections in hyperon vector form factors relevant for the extraction of Vus. A lattice quenched simulation has been performed, showing that it is possible to reach the required precision to extract SU(3)-breaking corrections in the regime of simulated masses. In order to perform the chiral extrapolation we calculated the chiral corrections to the vector form factor in HBChPT. Besides the one-loop O(p^2) contribution, we included also the subleading O(p^3) and O(1/M_B) corrections that, due to the Ademollo-Gatto theorem, are free from the contamination of unknown low energy constants. The results complete and correct previous calculations, and show that subleading corrections cannot be neglected. We also studied decuplet contributions within HBChPT and show that, in this case, the chiral expansion breaks down, rising doubts on the consistency of the theory.

preprint2005arXiv

Non-perturbatively renormalised light quark masses from a lattice simulation with N_f=2

We present results for the light quark masses obtained from a lattice QCD simulation with N_f=2 degenerate Wilson dynamical quark flavours. The sea quark masses of our lattice, of spacing a ~ 0.06 fm, are relatively heavy, i.e., they cover the range corresponding to 0.60 <~ M_P/M_V <~ 0.75. After implementing the non-perturbative RI-MOM method to renormalise quark masses, we obtain m_{ud}^{MS}(2 GeV)=4.3 +- 0.4^{+1.1}_{-0} MeV, and m_s^{MS}(2 GeV)=101 +- 8^{+25}_{-0} MeV, which are about 15% larger than they would be if renormalised perturbatively. In addition, we show that the above results are compatible with those obtained in a quenched simulation with a similar lattice.

preprint2005arXiv

Non-perturbatively Renormalized Light Quark Masses with Two Dynamical Fermions

We present the results of a partially quenched lattice QCD calculation of light quark masses with $N_f=2$ degenerate dynamical flavors. Numerical simulations are carried out using the plaquette gauge action and the Wilson quark action at $β= 5.8$ ($a^{-1} \simeq 3.2\gev$). The spatial extension of the $24^3 \times 48$ lattice is about 1.5 fm. Configurations have been generated at four values of the sea quark masses, for which the ratio of pseudoscalar over vector meson masses is in the range $M_P/M_V \simeq 0.60 ÷0.75$. An important feature of the present study is the use of non-perturbative renormalization, performed with the $\ri$ method. The effects of dynamical sea quarks in the determination of light quark masses have been investigated by performing a quenched calculation on a similar lattice. Our results for the average up-down and strange quark masses are $m_{ud}^{MS} (2\gev) = 4.3(4) (^{+1.1}_{-0.4})\mev$ and $m_s^{MS}(2\gev)=101(8)(^{+25}_{-9}) \mev$. These values are larger than those obtained by evaluating the quark mass renormalization constants with one-loop (boosted) perturbation theory. Our results for the light quark masses are compatible with those obtained in the quenched simulation. No significant sea quark effects are seen, at the values of sea quark masses used in the present study.

preprint2005arXiv

Operator product expansion and quark condensate from Lattice QCD in coordinate space

We present a Lattice QCD determination of the chiral quark condensate based on a new method. We extract the quark condensate from the operator product expansion of the quark propagator at short euclidean distances, where it represents the leading contribution in the chiral limit. From this study we obtain <\bar q q>^ms(2 GeV)=-(265+-5+-22 MeV)^3$, in good agreement with determinations of this quantity based on different approaches. The simulation is performed by using the O(a)-improved Wilson action at beta=6.45 on a volume 32^3\times70 in the quenched approximation.

preprint2005arXiv

The 2004 UTfit Collaboration Report on the Status of the Unitarity Triangle in the Standard Model

Using the latest determinations of several theoretical and experimental parameters, we update the Unitarity Triangle analysis in the Standard Model. The basic experimental constraints come from the measurements of |V_ub/V_cb|, Delta M_d, the lower limit on Delta M_s, epsilon_k, and the measurement of the phase of the B_d - anti B_d mixing amplitude through the time-dependent CP asymmetry in B^0 to J/psi K^0 decays. In addition, we consider the direct determination of alpha, gamma, 2 beta + gamma and cos(2 beta) from the measurements of new CP-violating quantities, recently performed at the B factories. We also discuss the opportunities offered by improving the precision of the various physical quantities entering in the determination of the Unitarity Triangle parameters. The results and the plots presented in this paper can also be found at http://www.utfit.org, where they are continuously updated with the newest experimental and theoretical results.

preprint2004arXiv

Light Hadron Spectrum, Renormalization Constants and Light Quark Masses with Two Dynamical Fermions

The results of a preliminary partially quenched (N_f=2) study of the light hadron spectrum, renormalization constants and light quark masses are presented. Numerical simulations are carried out with the LL-SSOR preconditioned Hybrid Monte Carlo with two degenerate dynamical fermions, using the plaquette gauge action and the Wilson quark action at beta = 5.8. Finite volume effects have been investigated employing two lattice volumes: 16^3 x 48 and 24^3 x 48. Configurations have been generated at four values of the sea quark mass corresponding to M_{PS}/M_V ~ 0.6 - 0.8.

preprint2004arXiv

The role of Lattice QCD in flavor physics

Understanding flavor physics is one of the most important tasks of particle physics today, which is motivating an extraordinary experimental and theoretical investigational effort. Important progress in this field has already been achieved in the last few years, with Lattice QCD calculations playing an essential role in this effort. I will describe some lattice contributions to the studies of flavor physics by focusing particularly on the determination of the CKM matrix and on the study of CP violation in the Standard Model.

preprint2003arXiv

Lifetime Differences and CP Violation Parameters of Neutral B Mesons at the Next-to-Leading Order in QCD

We compute the next-to-leading order QCD corrections to the off-diagonal elements of the decay-width matrix Gamma entering the neutral B-meson oscillations. From this calculation the width differences DeltaGamma and the CP violation parameters (q/p) of Bd and Bs mesons are estimated, including the complete O}(alpha_s) QCD corrections and the 1/mb contributions. For the width difference DeltaGamma_s we agree with previous results. By using the lattice determinations of the relevant hadronic matrix elements we obtain the theoretical predictions DeltaGamma_d/Gamma_d=(2.42+-0.59) x 10^(-3) and DeltaGamma_s/Gamma_s=(7.4+-2.4) x 10^(-2). For the CP violation parameters, we find |(q/p)_d|-1=(2.96+-0.67) x 10^(-4) and |(q/p)_s|-1=(-1.28+-0.27) x 10^(-5). These predictions are compatible with the experimental measurements which, however, suffer at present from large uncertainties.

preprint2003arXiv

Matrix elements of DI=3/2 K-->pi pi decays

We present a numerical computation of matrix elements of DI=3/2 K-->pi pi decays by using Wilson fermions. In order to extrapolate to the physical point we work at unphysical kinematics and we resort to Chiral Perturbation Theory at the next-to-leading order. In particular we explain the case of the electroweak penguins O_{7,8} which can contribute significantly in the theoretical prediction of epsilon'/epsilon. The study is done at beta=6.0 on a 24^3x64 lattice.

preprint2003arXiv

Unitarity Triangle Analysis in the Standard Model and Sensitivity to New Physics

By using the most recent determinations of the several theoretical and experimental input parameters, we update the Unitarity Triangle analysis in the Standard Model and discuss the sensitivity to New Physics effects. We investigate the interest of measuring with a better precision the various physical quantities entering the Unitarity Triangle analysis and study in a model independent way whether, despite the undoubted success of the CKM mechanism in the Standard Model, the Unitarity Triangle analysis still allows the presence of New Physics.

preprint2002arXiv

Kaon Weak Matrix Elements with Wilson Fermions

We present results of several numerical studies with Wilson fermions relevant for kaon physics. We compute the B_K parameter by using two different methods and extrapolate to the continuum limit. Our preliminary result is B_K(2 GeV)=0.66(7). Delta I=3/2 K->pi pi matrix elements are obtained by using the next-to-leading order expressions derived in chiral perturbation theory in which the low energy constants are determined by the lattice results computed at unphysical kinematics. From the simulation at beta=6.0 our (preliminary) results read: <pi pi|O_7(2 GeV)|K>_{I=2}=0.14(1)(1) GeV^3 and <pi pi|O_8(2GeV)|K>_{I=2}=0.69(6)(6) GeV^3.

preprint2002arXiv

Next-to-Leading Order QCD Corrections to Spectator Effects in Lifetimes of Beauty Hadrons

Theoretical predictions of beauty hadron lifetimes, based on the heavy quark expansion up to and including order 1/mb^2, do not to reproduce the experimental measurements of the lifetime ratios tau(B+)/tau(Bd) and tau(Lambdab)/tau(Bd). Large corrections to these predictions come from phase-space enhanced 1/mb^3 contributions, i.e. hard spectator effects. In this paper we calculate the next-to-leading order QCD corrections to the Wilson coefficients of the local operators appearing at O(1/mb^3). We find that these corrections improve the agreement with the experimental data. The lifetime ratio of charged to neutral B mesons, tau(B+)/tau(Bd), turns out to be in very good agreement with the corresponding measurement, whereas for tau(Bs)/tau(Bd) and tau(Lambdab)/tau(Bd) there is a residual difference at the 1-sigma level. We discuss, however, why the theoretical predictions are less accurate in the latter cases.

preprint2002arXiv

Non-perturbative renormalization constants and light quark masses

We present the results of an extensive non-perturbative calculation of the renormalization constants of bilinear quark operators for the non-perturbatively O(a)-improved Wilson action. The results are obtained at four values of the lattice coupling, by using the RI/MOM and the Ward identities methods. A new non-perturbative renormalization technique, which is based on the study of the lattice correlation functions at short distance in x-space, is also numerically investigated. We then use our non-perturbative determination of the quark mass renormalization constants to compute the values of the strange and the average up/down quark masses. After performing an extrapolation to the continuum limit, we obtain ms^msbar(2 GeV) = (106 +- 2 +- 8) MeV and ml^msbar(2 GeV)=(4.4 +- 0.1 +- 0.4) MeV.

preprint2001arXiv

2000 CKM-Triangle Analysis A Critical Review with Updated Experimental Inputs and Theoretical Parameters

Within the Standard Model, a review of the current determination of the sides and angles of the CKM unitarity triangle is presented, using experimental constraints from the measurements of |ε_K|, |V_{ub}/V_{cb}|, Δm_d and from the limit on Δm_s, available in September 2000. Results from the experimental search for {B}^0_s-\bar{B}^0_s oscillations are introduced in the present analysis using the likelihood. Special attention is devoted to the determination of the theoretical uncertainties. The purpose of the analysis is to infer regions where the parameters of interest lie with given probabilities. The BaBar "95 %, C.L. scanning" method is also commented.

preprint2001arXiv

Extraction of K-->pipi Matrix Elements with Wilson Fermions

We present the status of a lattice calculation for the K-->pipi matrix elements of the (delta S=1) effective weak Hamiltonian, directly with two pion in the final state. We study the energy shift of two pion in a finite volume both in the I=0 and I=2 channels. We explain a method to avoid the Goldstone pole contamination in the computation of renormalization constants for (delta I=3/2) operators. Finally we show some preliminary results for the matrix elements of (delta I=1/2) operators. Our quenched simulation is done at beta=6.0, with Wilson fermions, on a (24^3 X 64) lattice.

preprint2001arXiv

Matrix elements of (delta S=2) operators with Wilson fermions

We test the recent proposal of using the Ward identities to compute the K0-K0bar mixing amplitude with Wilson fermions, without the problem of spurious lattice subtractions. From our simulations, we observe no difference between the results obtained with and without subtractions. In addition, from the standard study of the complete set of (delta S=2) operators, we quote the following (preliminary) results (in the MS(NDR) scheme): Bk(2 GeV)=0.70(10), < O7^{3/2}>_{K->pi pi} = 0.10(2)(1) GeV^3, < O8^{3/2}>_{K->pi pi} = 0.49(6)(0) GeV^3.

preprint2000arXiv

A Theoretical Prediction of the Bs-Meson Lifetime Difference

We present the results of a quenched lattice calculation of the operator matrix elements relevant for predicting the Bs width difference. Our main result is (ΔΓ_Bs/Γ_Bs)= (4.7 +/- 1.5 +/- 1.6) 10^(-2), obtained from the ratio of matrix elements, R(m_b)=<\bar B_s^0|Q_S|B_s^0>/<\bar B_s^0|Q_L|B_s^0>=-0.93(3)^(+0.00)_(-0.01). R(m_b) was evaluated from the two relevant B-parameters, B_S^{MSbar}(m_b)=0.86(2)^(+0.02)_(-0.03) and B_Bs^{MSbar}(m_b) = 0.91(3)^(+0.00)_(-0.06), which we computed in our simulation.

preprint2000arXiv

First Lattice Calculation of the Electromagnetic Operator Amplitude <pi0|Q+|K0>

We present the first lattice calculation of the matrix element of the electromagnetic operator <pi0|Q+|K0>, where Q+ = (Q_d e/16 pi^2)* (\bar s_L sigma{mu,nu} F{mu,nu} d_R + \bar s_R sigma{mu,nu} F{mu,nu} d_L). This matrix element plays an important role, since it contributes to enhance the CP violating part of the K_L -> pi0 e+ e- amplitude in supersymmetric extensions of the Standard Model.

preprint1999arXiv

Quark Masses and Renormalization Constants from Quark Propagator and 3-point Functions

We have computed the light and strange quark masses and the renormalization constants of the quark bilinear operators, by studying the large-p^2 behaviour of the lattice quark propagator and 3-point functions. The calculation is non-perturbatively improved, at O(a), in the chiral limit. The method used to compute the quark masses has never been applied so far, and it does not require an explicit determination of the quark mass renormalization constant.

preprint1998arXiv

Improved Renormalization of Lattice Operators: A Critical Reappraisal

We systematically examine various proposals which aim at increasing the accuracy in the determination of the renormalization of two-fermion lattice operators. We concentrate on three finite quantities which are particularly suitable for our study: the renormalization constants of the vector and axial currents and the ratio of the renormalization constants of the scalar and pseudoscalar densities. We calculate these quantities in boosted perturbation theory, with several running boosted couplings, at the "optimal" scale q*. We find that the results of boosted perturbation theory are usually (but not always) in better agreement with non-perturbative determinations of the renormalization constants than those obtained with standard perturbation theory. The finite renormalization constants of two-fermion lattice operators are also obtained non-perturbatively, using Ward Identities, both with the Wilson and the tree-level Clover improved actions, at fixed cutoff ($β$=6.4 and 6.0 respectively). In order to amplify finite cutoff effects, the quark masses (in lattice units) are varied in a large interval 0<am<1. We find that discretization effects are always large with the Wilson action, despite our relatively small value of the lattice spacing ($a^{-1} \simeq 3.7$ GeV). With the Clover action discretization errors are significantly reduced at small quark mass, even though our lattice spacing is larger ($a^{-1} \simeq 2$ GeV). However, these errors remain substantial in the heavy quark region. We have implemented a proposal for reducing O(am) effects, which consists in matching the lattice quantities to their continuum counterparts in the free theory. We find that this approach still leaves appreciable, mass dependent, discretization effects.

preprint1998arXiv

Non-perturbatively Improved Heavy-Light Mesons: Masses and Decay Constants

We present a study of the heavy-light spectrum and of the D- and B-meson decay constants. The results wer e obtained in the quenched approximation, by using the non-perturbatively improved Clover lattice action at beta=6.2, with a sample of 100 configurations, on a 24^3 x 64 lattice. After a careful analysis of th e systematic errors present in the extraction of the physical results, by assuming quite conservative discretization errors, we find f_Ds=231 +/- 12^{+6}_{-1} MeV, f_D = 211 +/- 14^{+0}_{-12} MeV, f_Ds/f_D=1.10(2), f_Bs = 204 +/- 16^{+28}_{-0} MeV, f_B = 179 +/- 18^{+26}_{-9} MeV, f_Bs/f_B=1.14(3)^{+0}_{-1}. Our results, which have smaller discretization errors than many previous estimates at fixed value of the lattice spacing $a$, support a large value of f_B in the quenched approximation.

preprint1998arXiv

Non-perturbatively Renormalized Light-Quark Masses with the Alpha Action

We have computed the light quark masses using the O(a^2) improved Alpha action, in the quenched approximation. The renormalized masses have been obtained non-perturbatively. By eliminating the systematic error coming from the truncation of the perturbative series, our procedure removes the discrepancies, observed in previous calculations, between the results obtained using the vector and the axial-vector Ward identities. It also gives values of the quark masses larger than those obtained by computing the renormalization constants using (boosted) perturbation theory. Our main results, in the RI (MOM) scheme and at a renormalization scale μ=2 GeV, are m^{RI}_s= 138(15) MeV and m^{RI}_l= 5.6(5) MeV, where m^{RI}_s is the mass of the strange quark and m^{RI}_l=(m^{RI}_u+m^{RI}_d)/2 the average mass of the up-down quarks. From these results, which have been obtained non-perturbatively, by using continuum perturbation theory we derive the \bar{MS} masses, at the same scale, and the renormalization group invariant (m^{RGI}) masses. We find m^{NLO \bar{MS}}_s= 121(13)$ MeV and m^{NLO\bar{MS}}_l= 4.9(4) MeV at the next-to-leading order; m^{N^2LO \bar{MS}}_s= 111(12) MeV, m^{N^2LO \bar{MS}}_l= 4.5(4) MeV, m_s^{RGI}= 177(19) MeV and m^{RGI}_l= 7.2(6) MeV at the next-to-next-to-leading order.

preprint1995arXiv

DEPENDENCE OF THE CURRENT RENORMALISATION CONSTANTS ON THE QUARK MASS

We study the behaviour of the vector and axial current renormalisation constants $Z_V$ and $Z_A$ as a function of the quark mass, $m_q$. We show that sizeable $O(am_q)$ and $O(g_0^2 a m_q)$ systematic effects are present in the Wilson and Clover cases respectively. We find that the prescription of Kronfeld, Lepage and Mackenzie for correcting these artefacts is not always successful.