DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors

Abstract

In this study, we demonstrate the applicability of a recently proposed multiscale thermalization algorithm to two-color quantum chromodynamics (QCD) with two mass-degenerate fermion flavors. The algorithm involves refining an ensemble of gauge configurations that had been generated using a renormalization group (RG) matched coarse action, thereby producing a fine ensemble that is close to the thermalized distribution of a target fine action; the refined ensemble is subsequently rethermalized using conventional algorithms. Although the generalization of this algorithm from pure Yang-Mills theory to QCD with dynamical fermions is straightforward, we find that in the latter case, the method is susceptible to numerical instabilities during the initial stages of rethermalization when using the hybrid Monte Carlo algorithm. We find that these instabilities arise from large fermion forces in the evolution, which are attributed to an accumulation of spurious near-zero modes of the Dirac operator. We propose a simple strategy for curing this problem, and demonstrate that rapid thermalization--as probed by a variety of gluonic and fermionic operators--is possible with the use of this solution. Also, we study the sensitivity of rethermalization rates to the RG matching of the coarse and fine actions, and identify effective matching conditions based on a variety ofmore » measured scales.« less

Authors:
 [1];  [1]
  1. Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Center for Theoretical Physics
Publication Date:
Research Org.:
Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), High Energy Physics (HEP)
OSTI Identifier:
1424917
Alternate Identifier(s):
OSTI ID: 1334216
Grant/Contract Number:  
SC0010495; SC0011090
Resource Type:
Accepted Manuscript
Journal Name:
Physical Review D
Additional Journal Information:
Journal Volume: 94; Journal Issue: 11; Journal ID: ISSN 2470-0010
Publisher:
American Physical Society (APS)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS

Citation Formats

Detmold, William, and Endres, Michael G. Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors. United States: N. p., 2016. Web. doi:10.1103/PhysRevD.94.114502.
Detmold, William, & Endres, Michael G. Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors. United States. https://doi.org/10.1103/PhysRevD.94.114502
Detmold, William, and Endres, Michael G. Fri . "Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors". United States. https://doi.org/10.1103/PhysRevD.94.114502. https://www.osti.gov/servlets/purl/1424917.
@article{osti_1424917,
title = {Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors},
author = {Detmold, William and Endres, Michael G.},
abstractNote = {In this study, we demonstrate the applicability of a recently proposed multiscale thermalization algorithm to two-color quantum chromodynamics (QCD) with two mass-degenerate fermion flavors. The algorithm involves refining an ensemble of gauge configurations that had been generated using a renormalization group (RG) matched coarse action, thereby producing a fine ensemble that is close to the thermalized distribution of a target fine action; the refined ensemble is subsequently rethermalized using conventional algorithms. Although the generalization of this algorithm from pure Yang-Mills theory to QCD with dynamical fermions is straightforward, we find that in the latter case, the method is susceptible to numerical instabilities during the initial stages of rethermalization when using the hybrid Monte Carlo algorithm. We find that these instabilities arise from large fermion forces in the evolution, which are attributed to an accumulation of spurious near-zero modes of the Dirac operator. We propose a simple strategy for curing this problem, and demonstrate that rapid thermalization--as probed by a variety of gluonic and fermionic operators--is possible with the use of this solution. Also, we study the sensitivity of rethermalization rates to the RG matching of the coarse and fine actions, and identify effective matching conditions based on a variety of measured scales.},
doi = {10.1103/PhysRevD.94.114502},
journal = {Physical Review D},
number = 11,
volume = 94,
place = {United States},
year = {Fri Dec 02 00:00:00 EST 2016},
month = {Fri Dec 02 00:00:00 EST 2016}
}

Journal Article:

Citation Metrics:
Cited by: 5 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Lattice gauge fields, principal bundles and the calculation of topological charge
journal, December 1986

  • Phillips, Anthony; Stone, David
  • Communications in Mathematical Physics, Vol. 103, Issue 4
  • DOI: 10.1007/BF01211167

Multigrid Monte Carlo Method for Lattice Field Theories
journal, March 1986


Multigrid monte carlo in the sine gordon model
journal, April 1994


Chiral symmetry breaking in confining theories
journal, June 1980


Hybrid Monte Carlo
journal, September 1987


Multigrid preconditioning for the overlap operator in lattice QCD
journal, May 2015


Multigrid Monte Carlo algorithms for SU(2) lattice gauge theory: Two versus four dimensions
journal, December 1994


Analytic smearing of SU ( 3 ) link variables in lattice QCD
journal, March 2004


Multi-grid Monte Carlo (II). Two-dimensional XY model
journal, May 1991


Normalization of currents in lattice QCD
journal, April 1984


Dark nuclei. II. Nuclear spectroscopy in two-color QCD
journal, December 2014


Topology of lattice gauge fields
journal, August 1982

  • Lüscher, M.
  • Communications in Mathematical Physics, Vol. 85, Issue 1
  • DOI: 10.1007/BF02029132

Again on SU(3) glueball mass
journal, January 1985


Hamiltonian evolution for the hybrid Monte Carlo algorithm
journal, August 1992


The connection between local operators on the lattice and in the continuum and its relation to meson decay constants
journal, April 1983


Multi-grid Monte Carlo (III). Two-dimensional O(4)-symmetric non-linear σ-model
journal, August 1992


Multiscale Monte Carlo equilibration: Pure Yang-Mills theory
journal, December 2015


Speeding up the hybrid Monte Carlo algorithm for dynamical fermions
journal, October 2001


Critical slowing down and error analysis in lattice QCD simulations
journal, April 2011


Confinement of quarks
journal, October 1974


Adaptive Multigrid Algorithm for the Lattice Wilson-Dirac Operator
journal, November 2010


Metadynamics surfing on topology barriers: the CP N−1 case
journal, July 2016

  • Laio, A.; Martinelli, G.; Sanfilippo, F.
  • Journal of High Energy Physics, Vol. 2016, Issue 7
  • DOI: 10.1007/JHEP07(2016)089

Study of chiral symmetry in quenched QCD using the overlap Dirac operator
journal, April 1999


Multicanonical multigrid Monte Carlo method and effective autocorrelation time
journal, April 1994


Lattice regularization of gauge theories without loss of chiral symmetry
journal, April 1995


Topology of the SU (2) vacuum: a lattice study using improved cooling
journal, August 1997

  • de Forcrand, Philippe; García Pérez, Margarita; Stamatescu, Ion-Olimpiu
  • Nuclear Physics B, Vol. 499, Issue 1-2
  • DOI: 10.1016/S0550-3213(97)00275-7

Adaptive Aggregation-Based Domain Decomposition Multigrid for the Lattice Wilson--Dirac Operator
journal, January 2014

  • Frommer, A.; Kahl, K.; Krieg, S.
  • SIAM Journal on Scientific Computing, Vol. 36, Issue 4
  • DOI: 10.1137/130919507

Domain decomposition, multilevel integration, and exponential noise reduction in lattice QCD
journal, May 2016


Locality and exponential error reduction in numerical lattice gauge theory
journal, September 2001


PRIMME: preconditioned iterative multimethod eigensolver—methods and software description
journal, April 2010

  • Stathopoulos, Andreas; McCombs, James R.
  • ACM Transactions on Mathematical Software, Vol. 37, Issue 2
  • DOI: 10.1145/1731022.1731031

Kaon B -parameter from quenched domain-wall QCD
journal, May 2006


Finite renormalisations of currents in lattice gauge theories
journal, March 1983


Multilevel algorithm for flow observables in gauge theories
journal, April 2016


Properties and uses of the Wilson flow in lattice QCD
journal, August 2010


The Chroma Software System for Lattice QCD
journal, March 2005


Glueball masses and string tension in lattice QCD
journal, June 1987


Works referencing / citing this record:

Flow-based generative models for Markov chain Monte Carlo in lattice field theory
journal, August 2019


Flow-based generative models for Markov chain Monte Carlo in lattice field theory
text, January 2019