Related papers: Comparing the R algorithm and RHMC for staggered f…
Arrays of Rydberg atoms are a powerful platform to realize strongly-interacting quantum many-body systems. A common Rydberg Hamiltonian is free of the sign problem, meaning that its equilibrium properties are amenable to efficient…
Polynomial approximations to the inverse of the fermion matrix are used to filter the dynamics of the upper energy scales in HMC simulations. The use of a multiple time-scale integration scheme allows the filtered pseudofermions to be…
Complete spectra of the staggered Dirac operator $\Dirac$ are determined in four-dimensional $SU(2)$ gauge fields with and without dynamical fermions. An attempt is made to relate the performance of multigrid and conjugate gradient…
Three topics concerning fermion simulation algorithms are discussed: 1.) A performance comparison of the multiboson technique to simulate dynamical fermions and the Kramers equation algorithm, 2.) the question of reversibility in the Hybrid…
We introduce a numerical algorithm to stochastically sample the dual fermion perturbation series around the dynamical mean field theory, generating all topologies of two-particle interaction vertices. We show results in the weak and strong…
A Dirac choice for the averaging kernel $C$ is implemented numerically. This improved kernel will be needed in gauge covariant multigrid computations for propagators of staggered fermions. Results for $C$ and the variational coarse grid…
We compare the performances of the exact one-flavor algorithm (EOFA) and the rational hybrid Monte Carlo algorithm (RHMC), for dynamical simulations of lattice QCD with domain-wall fermion.
An exact, nonlocal, finite step-size algorithm for Monte Carlo simulation of theories with dynamical fermions is proposed. The algorithm is based on obtaining the new configuration U' from the old one U by solving the equation $ M(U') \eta…
We study discretization effects of the Wilson and staggered Dirac operator with $N_{\rm c}>2$ using chiral random matrix theory (chRMT). We obtain analytical results for the joint probability density of Wilson-chRMT in terms of a…
We apply the Hybrid Monte Carlo method to the simulation of overlap fermions. We give the fermionic force for the molecular dynamics update. We present early results on a small dynamical chiral ensemble.
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…
Continuous-time determinantal algorithm is proposed for the quantum Monte Carlo simulation of the interacting fermions. The scheme does not invoke Hubbard-Stratonovich transformation. The fermionic action is divided into two parts. One of…
Ridge regression (RR) is an important machine learning technique which introduces a regularization hyperparameter $\alpha$ to ordinary multiple linear regression for analyzing data suffering from multicollinearity. In this paper, we present…
Diagrammatic Monte Carlo (DiagMC) is a numeric technique that allows one to calculate quantities specified in terms of diagrammatic expansions, the latter being a standard tool of many-body quantum statistics. The sign problem that is…
In the simplified setting of the Schwinger model we present a systematic study on the simulation of dynamical fermions by global accept/reject steps that take into account the fermion determinant. A family of exact algorithms is developed,…
We propose improved estimators to compute the reweighting factors which are needed for lattice QCD calculations that rely on twisted-mass reweighting for the light quark contribution and the Rational Hybrid Monte Carlo (RHMC) algorithm for…
We report on simulations of QCD with many flavors of degenerate quarks, the DBW2 gauge action and naive staggered fermions, using the rational hybrid Monte Carlo algorithm. We primarily focus on eight degenerate quark flavors where a…
We propose a modification of the Hybrid-Monte-Carlo algorithm that allows for a larger step-size of the integration scheme at constant acceptance rate. The key ingredient is that the pseudo-fermion action is split into two parts. We test…
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…
Hamiltonian Monte Carlo (HMC) improves the computational efficiency of the Metropolis algorithm by reducing its random walk behavior. Riemannian Manifold HMC (RMHMC) further improves HMC's performance by exploiting the geometric properties…