Related papers: Comparing the R algorithm and RHMC for staggered f…
We present an exact local bosonic algorithm for the simulation of dynamical fermions in lattice QCD. We show that this algorithm is a viable alternative to the Hybrid Monte Carlo algorithm.
Correlated fermions are of high interest in condensed matter (Fermi liquids, Wigner molecules), cold atomic gases and dense plasmas. Here we propose a novel approach to path integral Monte Carlo (PIMC) simulations of strongly degenerate…
We present a mathematical framework for constructing and analyzing parallel algorithms for lattice Kinetic Monte Carlo (KMC) simulations. The resulting algorithms have the capacity to simulate a wide range of spatio-temporal scales in…
Riemannian manifold Hamiltonian Monte Carlo (RMHMC) is a sampling algorithm that seeks to adapt proposals to the local geometry of the posterior distribution. The specific form of the Hamiltonian used in RMHMC necessitates {\it…
We present a technique that enables the evaluation of perturbative expansions based on one-loop-renormalized vertices up to large expansion orders. Specifically, we show how to compute large-order corrections to the random phase…
We investigate exact symmetries of a staggered fermion in D dimensions. The Dirac operator is reformulated by SO(2D) Clifford algebra. The chiral symmetry, rotational invariance and parity symmetries are clarified in any dimension. Local…
Markov chain Monte Carlo algorithms are used to simulate from complex statistical distributions by way of a local exploration of these distributions. This local feature avoids heavy requests on understanding the nature of the target, but it…
The randomized Kaczmarz (RK) method is an iterative method for approximating the least-squares solution of large linear systems of equations. The standard RK method uses sequential updates, making parallel computation difficult. Here, we…
We investigate the validity of the square rooting procedure of the staggered determinant in the context of the Schwinger model. We find some evidence that at fixed physical quark mass the square root of the staggered determinant becomes…
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 report on our study of the Riemannian Manifold HMC (RMHMC) algorithm with the mass term for the gauge momenta replaced by rational functions of the gauge covariant Laplace operator. A comparison of HMC and RMHMC on a 2+1+1 flavor…
The Direct Simulation Monte Carlo (DSMC) method is widely employed for simulating rarefied nonequilibrium gas flows. With advances in aerospace engineering and micro/nano-scale technologies, gas flows exhibit the coexistence of rarefied and…
Round Robin (RR) scheduling algorithm is a preemptive scheduling algorithm. It is designed especially for time sharing Operating System (OS). In RR scheduling algorithm the CPU switches between the processes when the static Time Quantum…
We discuss the lattice formulation of gauge theories with fermions in arbitrary representations of the color group, and present the implementation of the RHMC algorithm for simulating dynamical Wilson fermions. A first dataset is presented…
We investigate the extension of the Prokof'ev-Svistunov worm algorithm to Wilson lattice fermions in an external scalar field. We effectively simulate by Monte Carlo the graphs contributing to the hopping expansion of the two-point function…
We describe a new Hybrid Monte Carlo (HMC) algorithm for dynamical overlap fermions, which improves the rate of topological index changes by adding an additional (intensive) term to the action for the molecular dynamics part of the…
The Rugged Metropolis (RM) algorithm is a biased updating scheme, which aims at directly hitting the most likely configurations in a rugged free energy landscape. Details of the one-variable (RM$_1$) implementation of this algorithm are…
We describe an HMC algorithm for dynamical overlap fermions which makes use of their good chiral properties. We test the algorithm in the Schwinger model. Topological sectors are readily changed even in the massless case.
The choice of appropriate interaction models is among the major disadvantages of conventional methods such as molecular dynamics and Monte Carlo simulations. On the other hand, the so-called reverse Monte Carlo (RMC) method, based on…
We discuss R package SQUAREM for accelerating iterative algorithms which exhibit slow, monotone convergence. These include the well-known expectation-maximization algorithm, majorize-minimize (MM), and other EM-like algorithms such as…