Related papers: The Tunneling Hybrid Monte-Carlo algorithm

We apply the topology conserving gauge action proposed by Luescher to the four-dimensional lattice QCD simulation in the quenched approximation. With this gauge action the topological charge is stabilized along the hybrid Monte Carlo…

We apply the UV-filtering preconditioner, previously used to improve the Multi-Boson algorithm, to the Polynomial Hybrid Monte Carlo (UV-PHMC) algorithm. The performance test for the algorithm is given for the plaquette gauge action and the…

We present a set of related Hybrid Monte Carlo methods to simulate an arbitrary number of dynamical overlap fermions. Each fermion is represented by a chiral pseudo-fermion field. The new algorithm reduces critical slowing down in the…

We perform dynamical QCD simulations with $n_f=2$ overlap fermions by hybrid Monte-Carlo method on $6^4$ to $8^3\times 16$ lattices. We study the problem of topological sector changing. A new method is proposed which works without…

The probability density of low-lying eigenvectors of the hermitian Wilson-Dirac operator is examined. Comparisons in position and size between eigenvectors, topological charge and action density are made. We do this for standard Monte-Carlo…

In this paper we propose new algorithm to reduce autocorrelation in Markov chain Monte-Carlo algorithms for euclidean field theories on the lattice. Our proposing algorithm is the Hybrid Monte-Carlo algorithm (HMC) with restricted Boltzmann…

We construct positive-definite pseudofermion actions for one fermion flavor in lattice field theory, for Wilson and domain-wall fermions respectively. The positive definiteness of these actions ensures that they can be simulated with the…

We study the spectrum of the hermitian Wilson Dirac operator in the epsilon-regime of QCD in the quenched approximation and compare it to predictions from Wilson Random Matrix Theory. Using the distributions of single eigenvalues in the…

We introduce modifications to Monte Carlo simulations of the Feynman path integral that improve sampling of localised interactions. The new algorithms generate trajectories in simple background potentials designed to concentrate them about…

The probability of accepting a candidate move in the hybrid Monte Carlo algorithm can be increased by considering a transition to be between windows of several states at the beginning and end of the trajectory, with a state within the…

UKQCD's dynamical fermion project uses the Generalised Hybrid Monte-Carlo (GHMC) algorithm to generate QCD gauge configurations for a non-perturbatively O(a) improved Wilson action with two degenerate sea-quark flavours. We describe our…

It has been suggested to project out a number of low-lying eigenvalues of the four-dimensional Wilson--Dirac operator that generates the transfer matrix of domain-wall fermions in order to improve simulations with domain-wall fermions. We…

In this work, we introduce three algorithmic improvements to reduce the cost and improve the scaling of orbital space variational Monte Carlo (VMC). First, we show that by appropriately screening the one- and two-electron integrals of the…

The efficiency of Hamiltonian Monte Carlo (HMC) can suffer when sampling a distribution with a wide range of length scales, because the small step sizes needed for stability in high-curvature regions are inefficient elsewhere. To address…

We present the preliminary tests on two modifications of the Hybrid Monte Carlo (HMC) algorithm. Both algorithms are designed to travel much farther in the Hamiltonian phase space for each trajectory and reduce the autocorrelations among…

Random Walk Metropolis Hastings (RWMH) algorithm, is quite inefficient in high dimensions because of its abysmally slow acceptance rate. The slow acceptance rate results from the fact that RWMH separately updates each coordinate of the…

We present an exact pseudofermion action for hybrid Monte Carlo simulation (HMC) of one-flavor domain-wall fermion (DWF), with the effective 4-dimensional Dirac operator equal to the optimal rational approximation of the overlap-Dirac…

We describe an adaptive multigrid algorithm for solving inverses of the domain-wall fermion operator. Our multigrid algorithm uses an adaptive projection of near-null vectors of the domain-wall operator onto coarser four-dimensional…

At stronger gauge-field couplings, the domain wall fermion (DWF) residual mass, a measure of chiral symmetry breaking, grows rapidly. This measure is largely due to near zero fermion eigenmodes of logarithm of the 4D transfer matrix along…

This paper surveys in detail the relations between numerical integration and the Hamiltonian (or hybrid) Monte Carlo method (HMC). Since the computational cost of HMC mainly lies in the numerical integrations, these should be performed as…