skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Lattice QCD simulations with light dynamical quarks

Abstract

I report recent results from full QCD simulations by CP-PACS and JLQCD collaborations.

Authors:
 [1];  [2]
  1. Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8571 (Japan)
  2. (United States)
Publication Date:
OSTI Identifier:
21056806
Resource Type:
Journal Article
Resource Relation:
Journal Name: AIP Conference Proceedings; Journal Volume: 892; Journal Issue: 1; Conference: QCHS7: 7. conference on quark confinement and the hadron spectrum, Ponta Delgada, Acores (Portugal), 2-7 Sep 2006; Other Information: DOI: 10.1063/1.2714339; (c) 2007 American Institute of Physics; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; COMPUTERIZED SIMULATION; LATTICE FIELD THEORY; NUMERICAL ANALYSIS; PERTURBED ANGULAR CORRELATION; QUANTUM CHROMODYNAMICS; QUARKS

Citation Formats

Aoki, Sinya, and Riken BNL Research Center, Physics 510A, BNL, Upton, NY 11973. Lattice QCD simulations with light dynamical quarks. United States: N. p., 2007. Web. doi:10.1063/1.2714339.
Aoki, Sinya, & Riken BNL Research Center, Physics 510A, BNL, Upton, NY 11973. Lattice QCD simulations with light dynamical quarks. United States. doi:10.1063/1.2714339.
Aoki, Sinya, and Riken BNL Research Center, Physics 510A, BNL, Upton, NY 11973. Tue . "Lattice QCD simulations with light dynamical quarks". United States. doi:10.1063/1.2714339.
@article{osti_21056806,
title = {Lattice QCD simulations with light dynamical quarks},
author = {Aoki, Sinya and Riken BNL Research Center, Physics 510A, BNL, Upton, NY 11973},
abstractNote = {I report recent results from full QCD simulations by CP-PACS and JLQCD collaborations.},
doi = {10.1063/1.2714339},
journal = {AIP Conference Proceedings},
number = 1,
volume = 892,
place = {United States},
year = {Tue Feb 27 00:00:00 EST 2007},
month = {Tue Feb 27 00:00:00 EST 2007}
}
  • We present the first results for lattice QCD at finite temperature T and chemical potential {mu} with four flavors of Wilson quarks. The calculations are performed using the imaginary chemical potential method at {kappa}=0,0.001,0.15,0.165,0.17 and 0.25, where {kappa} is the hopping parameter, related to the bare quark mass m and lattice spacing a by {kappa}=1/(2ma+8). Such a method allows us to do large scale Monte Carlo simulations at imaginary chemical potential {mu}=i{mu}{sub I}. By analytic continuation of the data with {mu}{sub I}<{pi}T/3 to real values of the chemical potential, we expect at each {kappa} set-membership sign [0,{kappa}{sub chiral}], a phasemore » transition line on the ({mu},T) plane, in a region relevant to the search for quark-gluon plasma in heavy-ion collision experiments. The transition is first order at small or large quark mass, and becomes a crossover at intermediate quark mass.« less
  • We calculate the B-meson decay constants f B, f Bs, and their ratio in unquenched lattice QCD using domain-wall light quarks and relativistic b-quarks. We use gauge-field ensembles generated by the RBC and UKQCD collaborations using the domain-wall fermion action and Iwasaki gauge action with three flavors of light dynamical quarks. We analyze data at two lattice spacings of a ≈ 0.11, 0.086 fm with unitary pion masses as light as M π ≈ 290 MeV; this enables us to control the extrapolation to the physical light-quark masses and continuum. For the b-quarks we use the anisotropic clover action withmore » the relativistic heavy-quark interpretation, such that discretization errors from the heavy-quark action are of the same size as from the light-quark sector. We renormalize the lattice heavy-light axial-vector current using a mostly nonperturbative method in which we compute the bulk of the matching factor nonperturbatively, with a small correction, that is close to unity, in lattice perturbation theory. We also improve the lattice heavy-light current through O(α sa). We extrapolate our results to the physical light-quark masses and continuum using SU(2) heavy-meson chiral perturbation theory, and provide a complete systematic error budget. We obtain f B0 = 196.2(15.7) MeV, f B+ = 195.4(15.8) MeV, f Bs = 235.4(12.2) MeV, f Bs/f B0 = 1.193(59), and f Bs/f B+ = 1.220(82), where the errors are statistical and total systematic added in quadrature. In addition, these results are in good agreement with other published results and provide an important independent cross check of other three-flavor determinations of B-meson decay constants using staggered light quarks.« less
  • We demonstrate a method for calculating the neutral B-meson decay constants and mixing matrix elements in unquenched lattice QCD with domain-wall light quarks and static b-quarks. Our computation is performed on the '2+1' flavor gauge configurations generated by the RBC and UKQCD Collaborations with a lattice spacing of a{approx_equal}0.11 fm (a{sup -1}=1.729 GeV) and a lattice spatial volume of approximately (1.8 fm){sup 3}. We simulate at three different light sea quark masses with pion masses down to approximately 430 MeV, and extrapolate to the physical quark masses using a phenomenologically-motivated fit function based on next-to-leading order heavy-light meson SU(2) chiralmore » perturbation theory. For the b-quarks, we use an improved formulation of the Eichten-Hill action with static link-smearing to increase the signal-to-noise ratio. We also improve the heavy-light axial current used to compute the B-meson decay constant to O({alpha}{sub s}pa) using one-loop lattice perturbation theory. We present initial results for the SU(3)-breaking ratios f{sub B{sub s}}/f{sub B{sub d}} and {xi}=f{sub B{sub s{radical}}}(B{sub B{sub s}})/f{sub B{sub d{radical}}}(B{sub B{sub d}}), thereby demonstrating the viability of the method. For the ratio of decay constants, we find f{sub B{sub s}}/f{sub B{sub d}}=1.15(12) and for the ratio of mixing matrix elements, we find {xi}=1.13(12), where in both cases the errors reflect the combined statistical and systematic uncertainties, including an estimate of the size of neglected O(1/m{sub b}) effects.« less
  • We calculate the form factors for B → πℓν and B s → Kℓν decay in dynamical lattice quantum chromodynamics (QCD) using domain-wall light quarks and relativistic b-quarks. We use the (2+1)-flavor gauge-field ensembles generated by the RBC and UKQCD collaborations with the domain-wall fermion action and Iwasaki gauge action. For the b-quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze data at two lattice spacings of a ≈ 0.11, 0.086 fm with unitary pion masses as light as M π ≈ 290 MeV. We simultaneously extrapolate our numerical results to the physical light-quark massesmore » and to the continuum and interpolate in the pion/kaon energy using SU(2) “hard-pion” chiral perturbation theory for heavy-light meson form factors. We provide complete systematic error budgets for the vector and scalar form factors f + (q 2) and f 0(q 2) for both B → πℓν and B s → Kℓν at three momenta that span the q 2 range accessible in our numerical simulations. Next we extrapolate these results to q 2 = 0 using a model-independent z-parametrization based on analyticity and unitarity. We present our final results for f +(q 2) and f 0(q 2)as the coefficients of the series in z and the matrix of correlations between them; this provides a parametrization of the form factors valid over the entire allowed kinematic range. Our results agree with other three-flavor lattice-QCD determinations using staggered light quarks, and have comparable precision, thereby providing important independent cross-checks. Both B → πℓν and B s → Kℓν decays enable determinations of the Cabibbo-Kobayashi-Maskawa matrix element |V ub|. Furthermore, we perform a combined z-fit of our numerical B → πℓν form-factor data with the experimental measurements of the branching fraction from BABAR and Belle leaving the relative normalization as a free parameter; we obtain |V ub| = 3.61(32)×10 -3, where the error includes statistical and all systematic uncertainties. The same approach can be applied to the decay B s → Kℓν to provide an alternative determination of |V ub| once the process has been measured experimentally. In anticipation of future experimental measurements, we make predictions for B → πℓν and B s → Kℓν differential branching fractions and forward-backward asymmetries in the Standard Model.« less