No Arabic abstract
We propose a new strategy to evaluate the partition function of lattice QCD with Wilson gauge action coupled to staggered fermions, based on a strong coupling expansion in the inverse bare gauge coupling $beta= 2N/g^{2}$. Our method makes use of the recently developed formalism to evaluate the ${rm SU}(N)$ $1-$link integrals and consists in an exact rewriting of the partition function in terms of a set of additional dual degrees of freedom which we call Decoupling Operator Indices (DOI). The method is not limited to any particular number of dimensions or gauge group ${rm U}(N)$, ${rm SU}(N)$. In terms of the DOI the system takes the form of a Tensor Network which can be simulated using Worm-like algorithms. Higher order $beta$-corrections to strong coupling lattice QCD can be, in principle, systematically evaluated, helping to answer the question whether the finite density sign problem remains mild when plaquette contributions are included. Issues related to the complexity of the description and strategies for the stochastic evaluation of the partition function are discussed.
Our knowledge about the QCD phase diagram at finite baryon chemical potential $mu_{B}$ is limited by the well known sign problem. The path integral measure, in the standard determinantal approach, becomes complex at finite $mu_{B}$ so that standard Monte Carlo techniques cannot be directly applied. As the sign problem is representation dependent, by a suitable choice of the fundamental degrees of freedom that parameterize the partition function, it can get mild enough so that reweighting techniques can be used. A successful formulation, capable to tame the sign problem, is known since decades in the limiting case $betato 0$, where performing the gauge integration first, gives rise to a dual formulation in terms of color singlets (MDP formulation). Going beyond the strong coupling limit represents a serious challenge as the gauge integrals involved in the computation are only partially known analytically and become strongly coupled for $beta>0$. We will present explict formulae for all the integral relevant for ${rm SU}(N)$ gauge theories discretised `a la Wilson, and will discuss how they can be used to obtain a positive dual formulation, valid for all $beta$, for pure Yang Mills theory.
The nucleon axial form factor is a dominant contribution to errors in neutrino oscillation studies. Lattice QCD calculations can help control theory errors by providing first-principles information on nucleon form factors. In these proceedings, we present preliminary results on a blinded calculation of $g_A$ and the axial form factor using HISQ staggered baryons with 2+1+1 flavors of sea quarks. Calculations are done using physical light quark masses and are absolutely normalized. We discuss fitting form factor data with the model-independent $z$ expansion parametrization.
Recent results from lattice QCD simulations provide a realistic picture, based upon first principles, of~$Upsilon$ physics. We combine these results with the experimentally measured mass of the $Upsilon$~meson to obtain an accurate and reliable value for the $b$-quarks pole mass. We use two different methods, each of which yields a mass consistent with $M_b = 5.0(2)$~GeV. This corresponds to a bare mass of $M_b^0 = 4.0(1)$~GeV in our lattice theory and an $msbar$~mass of $M_b^msbar(M_b)=4.0(1)$~GeV. We discuss the implications of this result for the $c$-quark mass. ******************************************************************************* THIS IS THE VERSION WHICH WILL BE PUBLISHED IN PRL. SUBSTANTIAL MATERIAL HAS BEEN ADDED, INCLUDING RESULTS WITH DYNAMICAL FERMIONS AND A CALCULATION OF THE MSBAR MASS. *******************************************************************************
The simulation of lattice QCD on massively parallel computers stimulated the development of scalable algorithms for the solution of sparse linear systems. We tackle the problem of the Wilson-Dirac operator inversion by combining a Schwarz alternating procedure (SAP) in multiplicative form with a flexible variant of the GMRES-DR algorithm. We show that restarted GMRES is not able to converge when the system is poorly conditioned. By adding deflation in the form of the FGMRES-DR algorithm, an important fraction of the information produced by the iterates is kept between successive restarts leading to convergence in cases in which FGMRES stagnates.
We present results from our simulations of quantum chromodynamics (QCD) with four flavors of quarks: u, d, s, and c. These simulations are performed with a one-loop Symanzik improved gauge action, and the highly improved staggered quark (HISQ) action. We are generating gauge configurations with four values of the lattice spacing ranging from 0.06 fm to 0.15 fm, and three values of the light quark mass, including the value for which the Goldstone pion mass is equal to the physical pion mass. We discuss simulation algorithms, scale setting, taste symmetry breaking, and the autocorrelations of various quantities. We also present results for the topological susceptibility which demonstrate the improvement of the HISQ configurations relative to those generated earlier with the asqtad improved staggered action.