Critical slowing down and error analysis in lattice QCD simulations
arXiv:1009.5228 · doi:10.1016/j.nuclphysb.2010.11.020
Abstract
We study the critical slowing down towards the continuum limit of lattice QCD simulations with Hybrid Monte Carlo type algorithms. In particular for the squared topological charge we find it to be very severe with an effective dynamical critical exponent of about 5 in pure gauge theory. We also consider Wilson loops which we can demonstrate to decouple from the modes which slow down the topological charge. Quenched observables are studied and a comparison to simulations of full QCD is made. In order to deal with the slow modes in the simulation, we propose a method to incorporate the information from slow observables into the error analysis of physical observables and arrive at safer error estimates.
33 pages including 15 figures and 4 tables
References in corpus (4)
- Perturbative analysis of the gradient flow in non-abelian gauge theories
- Accelerating Dynamical Fermion Computations using the Rational Hybrid Monte Carlo (RHMC) Algorithm with Multiple Pseudofermion Fields
- Exploring the HMC trajectory-length dependence of autocorrelation times in lattice QCD
- The Monte Carlo method in quantum field theory
Cited by in corpus (6)
- Probing the chiral regime of Nf=2 QCD with mixed actions
- Comparison of the mass preconditioned HMC and the DD-HMC algorithm for two-flavour QCD
- Non-renormalizability of the HMC algorithm
- Algorithms for lattice QCD: progress and challenges
- Wilson fermions at fine lattice spacings: scale setting, pion form factors and (g-2)_mu
- Autocorrelations in Hybrid Monte Carlo Simulations