Lattice QCD Calculations Research Articles

Strong-Isospin-Breaking Correction to the Muon Anomalous Magnetic Moment from Lattice QCD at the Physical Point.

All lattice-QCD calculations of the hadronic-vacuum-polarization contribution to the muon's anomalous magnetic moment to date have been performed with degenerate up- and down-quark masses. Here we calculate directly the strong-isospin-breaking correction to a_{μ}^{HVP} for the first time with physical values of m_{u} and m_{d} and dynamical u, d, s, and c quarks, thereby removing this important source of systematic uncertainty. We obtain a relative shift to be applied to lattice-QCD results obtained with degenerate light-quark masses of δa_{μ}^{HVP,m_{u}≠m_{d}}=+1.5(7)%, in agreement with estimates from phenomenology.

K→π matrix elements of the chromomagnetic operator on the lattice

We present the results of the first lattice QCD calculation of the $K \to \pi$ matrix elements of the chromomagnetic operator $O_{CM} = g\, \bar s\, \sigma_{\mu\nu} G_{\mu\nu} d$, which appears in the effective Hamiltonian describing $\Delta S = 1$ transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been determined non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with $N_f = 2+1+1$ dynamical quarks at three values of the lattice spacing. Our result for the B-parameter of the chromomagnetic operator at the physical pion and kaon point is $B_{CMO}^{K \pi} = 0.273 ~ (70)$, while in the SU(3) chiral limit we obtain $B_{CMO} = 0.072 ~ (22)$. Our findings are significantly smaller than the model-dependent estimate $B_{CMO} \sim 1 - 4$, currently used in phenomenological analyses, and improve the uncertainty on this important phenomenological quantity.

Analysis of angular observables of $\Lambda_b \to \Lambda (\to p\pi)\mu^{+}\mu^{-}$ decay in the standard and $Z^{\prime}$ models

In 2015, the LHCb collaboration has measured $\frac{d{\mathcal{B}}}{dq^2}$, the lepton- and hadron-side forward-backward asymmetries, denoted by $A^\ell_{FB}$ and $A^{\Lambda}_{FB}$, respectively in the range $15 < q^2(=s) < 20$ GeV$^2$ with 3 fb$^{-1}$ of data. Motivated by these measurements, we perform an analysis of $q^2$ dependent $\Lambda_b \to \Lambda (\to p \pi ) \mu^+\mu^-$ angular observables at large- and low-recoil in the SM and in a family non-universal $Z^{\prime}$ model. In the present study we use the recently performed high-precision lattice QCD calculations of the form factors that have well controlled uncertainties especially in $15 < s < 20$ GeV$^2$ bin. Using the full four-folded angular distribution of $\Lambda_b \to \Lambda (\to p \pi ) \mu^+\mu^-$ decay, firstly we calculate the values of these measured quantitites in the SM and compare their numerical values with the measurements in appropriate bins of $s$. In case of the possible discrepancy between the SM prediction and measurements, we try to see if these can be accommodated though the extra neutral $Z^{\prime}$ boson. In addition, the fraction of longitudinal polarization of the dimuon $F_{L}$ is measured to be $0.61^{+0.11}_{-0.14}\pm 0.03$ in $15 < s < 20$ GeV$^2$ at the LHCb. We find that in this bin the value found in the $Z^{\prime}$ model is close to the observed values. After comparing the results of these observables, we have proposed a number of other observables whose values are calculated in different bins of $s$ in the SM and $Z^{\prime}$ model. We illustrate that the experimental observations of these observables in several bins of $s$ can help to test the predictions of the SM and unravel NP contributions arises due to $Z^{\prime}$ model in these decays.

New methods for B meson decay constants and form factors from lattice NRQCD

We determine the normalisation of scalar and pseudoscalar current operators made from non-relativistic $b$ quarks and Highly Improved Staggered light quarks in lattice Quantum Chromodynamics (QCD) through $\mathcal{O}(\alpha_s)$ and $\Lambda_{\text{QCD}}/m_b$. We use matrix elements of these operators to extract $B$ meson decay constants and form factors, then compare to those obtained using the standard vector and axial-vector operators. This provides a test of systematic errors in the lattice QCD determination of the $B$ meson decay constants and form factors. We provide a new value for the $B$ and $B_s$ meson decay constants from lattice QCD calculations on ensembles that include $u$, $d$, $s$ and $c$ quarks in the sea and those which have the $u/d$ quark mass going down to its physical value. Our results are $f_B=0.196(6)$ GeV, $f_{B_s}=0.236(7)$ GeV and $f_{B_s}/f_B =1.207(7)$, agreeing well with earlier results using the temporal axial current. By combining with these previous results, we provide updated values of $f_B=0.190(4)$ GeV, $f_{B_s}=0.229(5)$ GeV and $f_{B_s}/f_B = 1.206(5)$.

Baryon chiral perturbation theory combined with the 1/Nc expansion in SU(3): Framework

Baryon chiral perturbation theory combined with the $1/{N}_{c}$ expansion is implemented for three flavors. Baryon masses, vector charges and axial vector couplings are studied to one-loop and organized according to the $\ensuremath{\xi}$-expansion, in which the $1/{N}_{c}$ and the low-energy power countings are linked according to $1/{N}_{c}=\mathcal{O}(\ensuremath{\xi})=\mathcal{O}(p)$. The renormalization to $\mathcal{O}({\ensuremath{\xi}}^{3})$ necessary for the mentioned observables is provided, along with applications to the baryon masses and axial couplings as obtained in lattice QCD calculations.

Domain wall fermion QCD with the exact one flavor algorithm

Lattice QCD calculations including the effects of one or more non-degenerate sea quark flavors are conventionally performed using the Rational Hybrid Monte Carlo (RHMC) algorithm, which computes the square root of the determinant of $\mathscr{D}^{\dagger} \mathscr{D}$, where $\mathscr{D}$ is the Dirac operator. The special case of two degenerate quark flavors with the same mass is described directly by the determinant of $\mathscr{D}^{\dagger} \mathscr{D}$ --- in particular, no square root is necessary --- enabling a variety of algorithmic developments, which have driven down the cost of simulating the light (up and down) quarks in the isospin-symmetric limit of equal masses. As a result, the relative cost of single quark flavors --- such as the strange or charm --- computed with RHMC has become more expensive. This problem is even more severe in the context of our measurements of the $\Delta I = 1/2$ $K \rightarrow \pi \pi$ matrix elements on lattice ensembles with $G$-parity boundary conditions, since $G$-parity is associated with a doubling of the number of quark flavors described by $\mathscr{D}$, and thus RHMC is needed for the isospin-symmetric light quarks as well. In this paper we report on our implementation of the exact one flavor algorithm (EOFA) introduced by the TWQCD collaboration for simulations including single flavors of domain wall quarks. We have developed a new preconditioner for the EOFA Dirac equation, which both reduces the cost of solving the Dirac equation and allows us to re-use the bulk of our existing high-performance code. Coupling these improvements with careful tuning of our integrator, the time per accepted trajectory in the production of our 2+1 flavor $G$-parity ensembles with physical pion and kaon masses has been decreased by a factor of 4.2.

Lattice QCD calculation of the B(s)→D(s)*ℓν form factors at zero recoil and implications for |Vcb|

We present results of a lattice QCD calculation of $B\to D^*$ and $B_s\to D_s^*$ axial vector matrix elements with both states at rest. These zero recoil matrix elements provide the normalization necessary to infer a value for the CKM matrix element $|V_{cb}|$ from experimental measurements of $\bar{B}^0\to D^{*+}\ell^-\bar{\nu}$ and $\bar{B}^0_s\to D_s^{*+}\ell^-\bar{\nu}$ decay. Results are derived from correlation functions computed with highly improved staggered quarks (HISQ) for light, strange, and charm quark propagators, and nonrelativistic QCD for the bottom quark propagator. The calculation of correlation functions employs MILC Collaboration ensembles over a range of three lattice spacings. These gauge field configurations include sea quark effects of charm, strange, and equal-mass up and down quarks. We use ensembles with physically light up and down quarks, as well as heavier values. Our main results are $\mathcal{F}^{B\to D^*}(1)= 0.895\pm 0.010_{\mathrm{stat}}\pm{{0.024}_{\mathrm{sys}}}$ and $\mathcal{F}^{B_s\to D_s^*}(1)= 0.883\pm 0.010_{\mathrm{stat}}\pm{0.028_{\mathrm{sys}}}$. We discuss the consequences for $|V_{cb}|$ in light of recent investigations into the extrapolation of experimental data to zero recoil.

Variance reduction and cluster decomposition

It is a common problem in lattice QCD calculation of the mass of the hadron with an annihilation channel that the signal falls off in time while the noise remains constant. In addition, the disconnected insertion calculation of the three-point function and the calculation of the neutron electric dipole moment with the $\theta$ term suffer from a noise problem due to the $\sqrt{V}$ fluctuation. We identify these problems to have the same origin and the $\sqrt{V}$ problem can be overcome by utilizing the cluster decomposition principle. We demonstrate this by considering the calculations of the glueball mass, the strangeness content in the nucleon, and the CP violation angle in the nucleon due to the $\theta$ term. It is found that for lattices with physical sizes of 4.5 - 5.5 fm, the statistical errors of these quantities can be reduced by a factor of 3 to 4. The systematic errors can be estimated from the Akaike information criterion. For the strangeness content, we find that the systematic error is of the same size as that of the statistical one when the cluster decomposition principle is utilized. This results in a 2 to 3 times reduction in the overall error.

Data-driven analysis for the temperature and momentum dependence of the heavy-quark diffusion coefficient in relativistic heavy-ion collisions

By applying a Bayesian model-to-data analysis, we estimate the temperature and momentum dependence of the heavy quark diffusion coefficient in an improved Langevin framework. The posterior range of the diffusion coefficient is obtained by performing a Markov chain Monte Carlo random walk and calibrating on the experimental data of $D$-meson $R_{\mathrm{AA}}$ and $v_2$ in three different collision systems at RHIC and the LHC: AuAu collisions at 200 GeV, PbPb collisions at 2.76 and 5.02 TeV. The spatial diffusion coefficient is found to be consistent with lattice QCD calculations and comparable with other models' estimation. We demonstrate the capability of our improved Langevin model to simultaneously describe the $R_{\mathrm{AA}}$ and $v_2$ at both RHIC and the LHC energies, as well as the higher order flow coefficient such as $D$-meson $v_3$. We show that by applying a Bayesian analysis, we are able to quantitatively and systematically study the heavy flavor dynamics in heavy-ion collisions.

Octet baryons in large magnetic fields

Magnetic properties of octet baryons are investigated within the framework of chiral perturbation theory. Utilizing a power counting for large magnetic fields, the Landau levels of charged mesons are treated exactly giving rise to baryon energies that depend non-analytically on the strength of the magnetic field. In the small-field limit, baryon magnetic moments and polarizabilities emerge from the calculated energies. We argue that the magnetic polarizabilities of hyperons provide a testing ground for potentially large contributions from decuplet pole diagrams. In external magnetic fields, such contributions manifest themselves through decuplet-octet mixing, for which possible results are compared in a few scenarios. These scenarios can be tested with lattice QCD calculations of the octet baryon energies in magnetic fields.

Exploring Partonic Structure of Hadrons Using abinitio Lattice QCD Calculations.

Following our previous proposal, we construct a class of good "lattice cross sections" (LCSs), from which we can study the partonic structure of hadrons from abinitio lattice QCD calculations. These good LCSs, on the one hand, can be calculated directly in lattice QCD, and on the other hand, can be factorized into parton distribution functions (PDFs) with calculable coefficients, in the same way as QCD factorization for factorizable hadronic cross sections. PDFs could be extracted from QCD global analysis of the lattice QCD generated data of LCSs. We also show that the proposed functions for lattice QCD calculation of PDFs in the literature are special cases of these good LCSs.

The gluon structure of hadrons and nuclei from lattice QCD

I discuss recent lattice QCD studies of the gluon structure of hadrons and light nuclei. After very briefly highlighting new determinations of the gluon contributions to the nucleon’s momentum and spin, presented by several collaborations over the last year, I describe first calculations of gluon generalised form factors. The generalised transversity gluon distributions are of particular interest since they are purely gluonic; they do not mix with quark distributions at leading twist. In light nuclei they moreover provide a clean signature of non-nucleonic gluon degrees of freedom, and I present the first evidence for such effects, based on lattice QCD calculations. The planned Electron-Ion Collider, designed to access gluon structure quantities, will have the capability to test this prediction, and measure a range of gluon observables including generalised gluon distributions and transverse momentum dependent gluon distributions, within the next decade.

Heavy-light mesons in a contact interaction

We present a unified formalism for the analysis of mesons. It is provided by a symmetry preserving Schwinger-Dyson-Bethe-Salpeter-Equation (SDBSE) treatment of a vector-vector contact interaction. The contact interaction (CI) model offers a simple-to-implement alternative to perform exploratory studies of QCD, within the SDBSE framework. Within the limitation of this model, we calculate the spectrum of D and B mesons. Our results are in agreement within 2% when compare with experimental data, lattice QCD and other SDBSE calculations involving sophisticated interaction kernels.

HQE parameters from unquenched lattice data on pseudoscalar and vector heavy-light meson masses

We present a new lattice determination of some of the parameters appearing both in the Operator Product Expansion (OPE) analysis of the inclusive semileptonic B-meson decays and in the Heavy Quark Expansion (HQE) of the pseudoscalar (PS) and vector (V) heavy-light meson masses. We perform a lattice QCD (LQCD) computation of PS and V heavy-light meson masses for heavy-quark masses mh in the range from [see formula in PDF] to ≃ [see formula in PDF]. We employed the Nf = 2 + 1 + 1 gauge configurations of the European Twisted Mass Collaboration (ETMC) at three values of the lattice spacing a ≃ (0.062,0.082,0.089) fm with pion masses in the range Mπ ≃ (210 - 450) MeV. The heavy-quark mass is simulated directly on the lattice up to ≃ [see formula in PDF]. The interpolation to the physical [see formula in PDF] is performed using the ETMC ratio method and adopting the kinetic mass scheme. We obtain [see formula in PDF] (1 GeV) = 4.61(20) GeV (m̅b(m̅b) = 4.26(18) GeV in the MS scheme). The lattice data are analyzed in terms of the HQE and the matrix elements of dimension-4 and dimension-5 operators are extracted with good precision, namely: [see formula in PDF] and [see formula in PDF]. The data also allow for an estimate of the dimension-6 operator matrix elements.

Thermo-magnetic dependent couplings for quark matter

In this work, we improve our model by adding the strange quark and the Polyakov loop. We then build a thermo-magnetic dependence for the G and K couplings of the SU(3) PNJL model by fitting lattice QCD calculations for the average and the difference of u and d quark condensates under a strong magnetic field.

Nucleon axial coupling from Lattice QCD

We present state-of-the-art results from a lattice QCD calculation of the nucleon axial coupling, gA, using Möbius Domain-Wall fermions solved on the dynamical Nf = 2 + 1 + 1 HISQ ensembles after they are smeared using the gradient-flow algorithm. Relevant three-point correlation functions are calculated using a method inspired by the Feynman-Hellmann theorem, and demonstrate significant improvement in signal for fixed stochastic samples. The calculation is performed at five pion masses of mπ ~ {400, 350, 310, 220, 130} MeV, three lattice spacings of a ~ {0.15, 0.12, 0.09} fm, and we do a dedicated volume study with mπL ~ {3.22, 4.29, 5.36}. Control over all relevant sources of systematic uncertainty are demonstrated and quantified. We achieve a preliminary value of gA = 1.285(17), with a relative uncertainty of 1.33%.

Lattice QCD Application Development within the US DOE Exascale Computing Project

In October, 2016, the US Department of Energy launched the Exascale Computing Project, which aims to deploy exascale computing resources for science and engineering in the early 2020’s. The project brings together application teams, software developers, and hardware vendors in order to realize this goal. Lattice QCD is one of the applications. Members of the US lattice gauge theory community with significant collaborators abroad are developing algorithms and software for exascale lattice QCD calculations. We give a short description of the project, our activities, and our plans.

Neutron Electric Dipole Moment on the Lattice

For the neutron to have an electric dipole moment (EDM), the theory of nature must have T, or equivalently CP, violation. Neutron EDM is a very good probe of novel CP violation in beyond the standard model physics. To leverage the connection between measured neutron EDM and novel mechanism of CP violation, one requires the calculation of matrix elements for CP violating operators, for which lattice QCD provides a first principle method. In this paper, we review the status of recent lattice QCD calculations of the contributions of the QCD Θ-term, the quark EDM term, and the quark chromo-EDM term to the neutron EDM.

D → Klv semileptonic decay using lattice QCD with HISQ at physical pion masses

he quark flavor sector of the Standard Model is a fertile ground to look for new physics effects through a unitarity test of the Cabbibo-Kobayashi-Maskawa (CKM) matrix. We present a lattice QCD calculation of the scalar and the vector form factors (over a large q2 region including q2 = 0) associated with the D→ Klv semi-leptonic decay. This calculation will then allow us to determine the central CKM matrix element, Vcs in the Standard Model, by comparing the lattice QCD results for the form factors and the experimental decay rate. This form factor calculation has been performed on the Nf = 2 + 1 + 1 MILC HISQ ensembles with the physical light quark masses.

Calm Multi-Baryon Operators

There are many outstanding problems in nuclear physics which require input and guidance from lattice QCD calculations of few baryons systems. However, these calculations suffer from an exponentially bad signal-to-noise problem which has prevented a controlled extrapolation to the physical point. The variational method has been applied very successfully to two-meson systems, allowing for the extraction of the two-meson states very early in Euclidean time through the use of improved single hadron operators. The sheer numerical cost of using the same techniques in two-baryon systems has so far been prohibitive. We present an alternate strategy which offers some of the same advantages as the variational method while being significantly less numerically expensive. We first use the Matrix Prony method to form an optimal linear combination of single baryon interpolating fields generated from the same source and different sink interpolating fields. Very early in Euclidean time this optimal linear combination is numerically free of excited state contamination, so we coin it a calm baryon. This calm baryon operator is then used in the construction of the two-baryon correlation functions.To test this method, we perform calculations on the WM/JLab iso-clover gauge configurations at the SU(3) flavor symmetric point with mπ~ 800 MeV — the same configurations we have previously used for the calculation of two-nucleon correlation functions. We observe the calm baryon significantly removes the excited state contamination from the two-nucleon correlation function to as early a time as the single-nucleon is improved, provided non-local (displaced nucleon) sources are used. For the local two-nucleon correlation function (where both nucleons are created from the same space-time location) there is still improvement, but there is significant excited state contamination in the region the single calm baryon displays no excited state contamination.