Source author record

S. Simula

S. Simula 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

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

83 published item(s)

preprint2022arXiv

$\vert V_{cb} \vert$ and $R(D^{(*)})$ using lattice QCD and unitarity

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

preprint2022arXiv

Exclusive B-meson semileptonic decays from unitarity and lattice QCD

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

preprint2022arXiv

Exclusive semileptonic $B \to π\ell ν_\ell$ and $B_s \to K \ell ν_\ell$ decays through unitarity and lattice QCD

The Cabibbo-Kobayashi-Maskawa (CKM) matrix element $\vert V_{ub}\vert$ is obtained from exclusive semileptonic $B \to π\ell ν_\ell$ and $B_s \to K \ell ν_\ell$ decays adopting the unitarity-based dispersion matrix approach for the determination of the hadronic form factors (FFs) in the whole kinematical range. We use lattice computations of the relevant susceptibilities and of the FFs in the large-$q^2$ regime in order to derive their behavior in the low-$q^2$ region without assuming any specific momentum dependence and without constraining their shape using experimental data. Then, we address the extraction of $\vert V_{ub}\vert$ from the experimental data, obtaining $\vert V_{ub}\vert = (3.62 \pm 0.47) \cdot 10^{-3}$ from $B \to π$ and $\vert V_{ub}\vert = (3.77 \pm 0.48) \cdot 10^{-3}$ from $B_s \to K$, which after averaging yield $\vert V_{ub}\vert = (3.69 \pm 0.34) \cdot 10^{-3}$. These results are compatible with the most recent inclusive value $\vert V_{ub} \vert_{incl} = 4.13\,(26) \cdot 10^{-3}$ at the 1$σ$ level. We also present purely theoretical estimates of the ratio of the $τ/μ$ decay rates $R^{τ/μ}_{π(K)}$, the normalized forward-backward asymmetry $\bar{\mathcal{A}}_{FB}^{\ell,π(K)}$ and the normalized lepton polarization asymmetry $\bar{\mathcal{A}}_{polar}^{\ell,π(K)}$.

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

preprint2020arXiv

FLAG Review 2019

We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to π$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_π$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ 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. We review the status of lattice determinations of the strong coupling constant $α_s$. Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal.

preprint2020arXiv

Ratios of the hadronic contributions to the lepton $g-2$ from Lattice QCD+QED simulations

The ratios among the leading-order (LO) hadronic vacuum polarization (HVP) contributions to the anomalous magnetic moments of electron, muon and tau-lepton, $a_{\ell=e,μτ}^{HVP,LO}$, are computed using lattice QCD+QED simulations. The results include the effects at order $O(α_{em}^2)$ as well as the electromagnetic and strong isospin-breaking corrections at orders $O(α_{em}^3)$ and $O(α_{em}^2(m_u-m_d))$, respectively, where $(m_u-m_d)$ is the $u$- and $d$-quark mass difference. We employ the gauge configurations generated by the Extended Twisted Mass Collaboration with $N_f=2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range 210 - 450 MeV. We show that in the case of the electron-muon ratio the hadronic uncertainties in the numerator and in the denominator largely cancel out, while in the cases of the electron-tau and muon-tau ratios such a cancellation does not occur. For the electron-muon ratio we get $R_{e/μ} \equiv (m_μ/m_e)^2 (a_e^{HVP,LO} / a_μ^{HVP,LO}) = 1.1456~(83)$ with an uncertainty of $\simeq 0.7 \%$. Our result, which represents an accurate Standard Model (SM) prediction, agrees very well with the estimate obtained using the results of dispersive analyses of the experimental $e^+ e^- \to$ hadrons data. Instead, it differs by $\simeq 2.7$ standard deviations from the value expected from present electron and muon (g - 2) experiments after subtraction of the current estimates of the QED, electro-weak, hadronic light-by-light and higher-order HVP contributions, namely $R_{e/μ} = 0.575~(213)$. An improvement of the precision of both the experiment and the QED contribution to the electron (g - 2) by a factor of $\simeq 2$ could be sufficient to reach a tension with our SM value of the ratio $R_{e/μ}$ at a significance level of $\simeq 5$ standard deviations.

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.

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

QCD sum-rule results for heavy-light meson decay constants and comparison with lattice QCD

Updated predictions for the decay constants of the D, Ds, B and Bs mesons obtained from Borel QCD sum rules for heavy-light currents are presented and compared with the recent lattice averages performed by the Flavor Lattice Averaging Group. An excellent agreement is obtained in the charm sector, while some tension is observed in the bottom sector. Moreover, available lattice and QCD sum-rule calculations of the decay constants of the vector D*, Ds*, B* and Bs* mesons are compared. Again some tension in the bottom sector is observed.

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

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

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

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

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.

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

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

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

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

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.

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

Measurement of the Deuteron Structure Function F2 in the Resonance Region and Evaluation of Its Moments

Inclusive electron scattering off the deuteron has been measured to extract the deuteron structure function F2 with the CEBAF Large Acceptance Spectrometer (CLAS) at the Thomas Jefferson National Accelerator Facility. The measurement covers the entire resonance region from the quasi-elastic peak up to the invariant mass of the final-state hadronic system W~2.7 GeV with four-momentum transfers Q2 from 0.4 to 6 (GeV/c)^2. These data are complementary to previous measurements of the proton structure function F2 and cover a similar two-dimensional region of Q2 and Bjorken variable x. Determination of the deuteron F2 over a large x interval including the quasi-elastic peak as a function of Q2, together with the other world data, permit a direct evaluation of the structure function moments for the first time. By fitting the Q2 evolution of these moments with an OPE-based twist expansion we have obtained a separation of the leading twist and higher twist terms. The observed Q2 behaviour of the higher twist contribution suggests a partial cancellation of different higher twists entering into the expansion with opposite signs. This cancellation, found also in the proton moments, is a manifestation of the "duality" phenomenon in the F2 structure function.

preprint2006arXiv

Neutron structure function moments at leading twist

The experimental data on F2 structure functions of the proton and deuteron were used to construct their moments. In particular, recent measurements performed with CLAS detector at Jefferson Lab allowed to extend our knowledge of structure functions in the large-x region. The phenomenological analysis of these experimental moments in terms of the Operator Product Expansion permitted to separate the leading and higher twist contributions. Applying nuclear corrections to extracted deuteron moments we obtained the contribution of the neutron. Combining leading twist moments of the neutron and proton we found d/u ratio at x->1 approaching 0, although 1/5 value could not be excluded. The twist expansion analysis suggests that the contamination of higher twists influences the extraction of the d/u ratio at x->1 even at Q2-scale as large as 12 (GeV/c)^2.

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

Leading twist moments of the neutron structure function F2n

We perform a global analysis of neutron F2n structure function data, obtained by combining proton and deuteron measurements over a large range of kinematics. From these data the lowest moments (n <= 10) of the leading twist neutron F2n structure function are extracted. Particular attention is paid to nuclear effects in the deuteron, which become increasingly important for the higher moments. Our results for the nonsinglet, isovector (p - n) combination of the leading twist moments are used to test recent lattice simulations. We also determine the lowest few moments of the higher twist contributions, and find these to be approximately isospin independent, suggesting the possible dominance of ud correlations over uu and dd in the nucleon.

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.

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.

preprint2003arXiv

A Kinematically Complete Measurement of the Proton Structure Function F2 in the Resonance Region and Evaluation of Its Moments

We measured the inclusive electron-proton cross section in the nucleon resonance region (W < 2.5 GeV) at momentum transfers Q**2 below 4.5 (GeV/c)**2 with the CLAS detector. The large acceptance of CLAS allowed for the first time the measurement of the cross section in a large, contiguous two-dimensional range of Q**2 and x, making it possible to perform an integration of the data at fixed Q**2 over the whole significant x-interval. From these data we extracted the structure function F2 and, by including other world data, we studied the Q**2 evolution of its moments, Mn(Q**2), in order to estimate higher twist contributions. The small statistical and systematic uncertainties of the CLAS data allow a precise extraction of the higher twists and demand significant improvements in theoretical predictions for a meaningful comparison with new experimental results.

preprint2003arXiv

Analysis of the Axial Anomaly on the Lattice with O(a)-improved Wilson Action

Flavor singlet and non-singlet axial Ward identities are investigated using the Wilson formulation of lattice QCD with Clover O(a)-improvement, which breaks explicitly chiral symmetry. The matching at one-loop order of all the relevant renormalization constants with the continuum \bar{MS} scheme is presented. Our calculations include: 1) the contributions arising from the Clover term of the action; 2) the complete mixing of the gluon operator G \tilde{G} with the divergence of the singlet axial current; 3) the use of both local and extended definitions of the fermionic bilinear operators. A definition of the gluon operator G \tilde{G} on the lattice outside the chiral limit is proposed. Our definition takes into account the possible power-divergent mixing with the pseudoscalar density, generated by the breaking of chiral symmetry. A non-perturbative procedure for the evaluation of such mixing constant is worked out. Finally, the renormalization properties of the composite insertion of the topological charge operator relevant for the lattice calculation of the neutron electric dipole moment, induced by the strong CP-violating term of the QCD Lagrangian, are discussed.

preprint2003arXiv

New extraction of alpha_s(MZ) from proton DIS data

An exploratory study for a new determination of the strong coupling constant alpha_s(MZ) from existing world data on the proton structure function F2p in the Q**2-range from ~ 5 to ~ 120 (GeV/c)**2 is presented. The main features of our approach are: 1) the use of low-order Nachtmann moments evaluated with a direct contribution from data larger than 70% of the total; 2) the inclusion of high-order perturbative effects through the soft gluon resummation technique; 3) a direct control over higher-twist effects; and 4) the independence from any specific choice of the x-shape of the input parton distributions. At next-to-leading order we get alpha_s(MZ) = 0.1209 +/- 0.0010 (stat.) +/- 0.0015 (syst.) with a significant dependence upon the order of the moment used. Including soft gluon effects we obtain alpha_s(MZ) = 0.1188 +/- 0.0010 (stat.) +/- 0.0014 (syst.) with a remarkable better stability against the order of the moment. Our findings are compared with recent DIS determinations of alpha_s(MZ) and directions for future improvements are discussed.

preprint2003arXiv

Possible evidence of extended objects inside the proton

Recent experimental determinations of the Nachtmann moments of the inelastic structure function of the proton F2p(x, Q**2), obtained at Jefferson Lab, are analyzed for values of the squared four-momentum transfer Q**2 ranging from ~ 0.1 to ~ 2 (GeV/c)**2. It is shown that such inelastic proton data exhibit a new type of scaling behavior and that the resulting scaling function can be interpreted as a constituent form factor consistent with the elastic nucleon data. These findings suggest that at low momentum transfer the inclusive proton structure function originates mainly from the elastic coupling with extended objects inside the proton. We obtain a constituent size of ~ 0.2 - 0.3 fm.

preprint1999arXiv

A connection between inclusive semileptonic decays of bound and free heavy quarks

A relativistic constituent quark model, formulated on the light-front, is used to derive a new parton approximation for the inclusive semileptonic decay width of the B-meson. A simple connection between the decay rate of a free heavy-quark and the one of a heavy-quark bound in a meson or in a baryon is established. The main features of the new approach are the treatment of the b-quark as an on-mass-shell particle and the inclusion of the effects arising from the b-quark transverse motion in the B-meson. In a way conceptually similar to the deep-inelastic scattering case, the B-meson inclusive width is expressed as the integral of the free b-quark partial width multiplied by a bound-state factor related to the b-quark distribution function in the B-meson. The non-perturbative meson structure is described through various quark-model wave functions, constructed via the Hamiltonian light-front formalism using as input both relativized and non-relativistic potential models. A link between spectroscopic quark models and the B-meson decay physics is obtained in this way. Our predictions for the B -> X_c l nu_l and B -> X_u l nu_l decays are used to extract the CKM parameters |V_cb| and |V_ub| from available inclusive data. After averaging over the various quark models adopted and including leading-order perturbative QCD corrections, we obtain |V_cb| = (43.0 +/- 0.7_exp +/- 1.8_th) 10^-3 and |V_ub| = (3.83 +/- 0.48_exp +/- 0.14_th) 10^-3, implying |V_ub / V_cb| = 0.089 +/- 0.011_exp +/- 0.005_th, in nice agreement with existing predictions.

preprint1999arXiv

Bloom-Gilman duality of the nucleon structure function and the elastic peak contribution

The occurrence of the Bloom-Gilman duality in the nucleon structure function is investigated by analyzing the Q**2-behavior of low-order moments, both including and excluding the contribution arising from the nucleon elastic peak. The Natchmann definition of the moments has been adopted in order to cancel out target-mass effects. It is shown that the onset of the Bloom-Gilman duality occurs around Q**2 ~ 2 (GeV/c)**2 if only the inelastic part of the nucleon structure function is considered, whereas the inclusion of the nucleon elastic peak contribution leads to remarkable violations of the Bloom-Gilman duality.

preprint1998arXiv

Light-Baryon Spectroscopy and the Electromagnetic Form Factors in the Quark Model

The momentum distributions of the constituent quarks inside the nucleon and the prominent electroproduced nucleon resonances are investigated in the two most sophisticated, available quark potential models, based respectively on the assumption of the valence + gluon dominance and on the exchange of the pseudoscalar Goldstone-bosons arising from the spontaneous breaking of chiral symmetry. It is shown that both models predict a large, similar content of high-momentum components, due to the short-range part of the interquark interaction, which affect the behaviour of both elastic and transition electromagnetic form factors at large values of the momentum transfer. The electromagnetic form factors are calculated within a relativistic approach formulated on the light-front, adopting a one-body current with constituent quark form factors. The results suggest that soft, non-perturbative effects can play a relevant role for explaining the existing data on elastic as well as transition form factors (at least) for Q**2 ~ 10 - 20 (GeV/c)**2.

preprint1995arXiv

Electromagnetic form factors of the $ρ$ meson in a light-front constituent quark model

The electromagnetic form factors of the $ρ$ meson are evaluated adopting a relativistic constituent quark model based on the light-front formalism, and using a meson wave function with the high-momentum tail generated by the one-gluon-exchange interaction. The breakdown of the rotational covariance for the one-body component of the current operator is investigated and the sensitivity of the ratio of the $ρ$-meson form factors to the pion (charge) form factor to the spin-dependent component of the effective $q \bar{q}$ interaction is illustrated.

preprint1995arXiv

Electromagnetic Structure of Mesons in a Light-Front Constituent Quark Model

The elastic and transition electromagnetic form factors of mesons are investigated within a light-front constituent quark model. The general formulae for the space-like matrix elements of the one-body component of the electromagnetic current, including both Dirac and Pauli form factors of the constituent quarks, are presented. Our results for the pion charge form factor and the $πρ$ transition form factor are reported and compared with those obtained within various relativistic approaches in a range of values of the four-momentum transfer accessible to $CEBAF$.

preprint1995arXiv

Nucleon and pion electromagnetic structure and constituent quark form factors

We have evaluated nucleon and pion electromagnetic form factors using for the first time eigenfunctions of a light-front mass operator, reproducing a large set of energy levels, and constituent quark form factors in the one-body current. The comparison with the experimental data yields valuable information on the electromagnetic structure of light constituent quarks and on the different role played by their Dirac and Pauli form factors.

preprint1995arXiv

Nucleon-Nucleon Correlations and Inclusive Electron Scattering off Few-Nucleon Systems at $x>1$

Inclusive electron scattering off few-nucleon systems is investigated at $x>1$ and high momentum transfer, including the contributions from quasi-elastic and deep inelastic scattering processes. It is shown that at $x>1$ the inclusive cross section is dominated by the process of virtual photon absorption on a pair of correlated nucleons both in case of quasi-elastic and deep inelastic scattering events. The sensitivity of the nuclear response to the effects arising from the possible presence of multiquark cluster configurations at short internucleon separations as well as from possible medium-dependent modifications of the nucleon structure function, is illustrated.

preprint1995arXiv

Nucleon-Nucleon Correlations and Multiquark Cluster Effects in Deep Inelastic Electron Scattering off Few-Nucleon Systems at $x>1$

Inclusive $A(e,e')X$ and semi-inclusive $A(e,e'N)X$ deep inelastic electron scattering processes off few-nucleon systems are investigated at $x > 1$, showing some of the relevant features of the cross section which are sensitive to the effects arising from nucleon-nucleon correlations and possible exotic multiquark cluster configurations at short internucleon separations.