Partial-Global Stochastic Metropolis Update for Dynamical Smeared Link Fermions

TL;DR

This paper introduces an efficient partial-global stochastic Metropolis algorithm for dynamical simulations with HYP smeared staggered fermions, demonstrating improved autocorrelation times especially with two flavors, enabling larger volume simulations.

Contribution

The paper presents a novel implementation of the partial-global stochastic Metropolis algorithm with fermion matrix reduction and determinant breakup for improved efficiency in lattice QCD simulations.

Findings

Efficient autocorrelation times with two flavors.

Feasibility of larger volume simulations.

Parameter choices optimize performance.

Abstract

We performed dynamical simulations with HYP smeared staggered fermions using the recently proposed partial-global stochastic Metropolis algorithm with fermion matrix reduction and determinant breakup improvements. In this paper we discuss our choice of the action parameters and study the autocorrelation time both with four and two fermionic flavors at different quark mass values on approximately 10 fm^4 lattices. We find that the update is especially efficient with two flavors making simulations on larger volumes feasible.