Related papers: Parallel tempering in full QCD with Wilson fermion…
We study efficiency of higher order integrator schemes for the hybrid Monte Carlo (HMC) algorithm. Numerical tests are performed for Quantum Chromo Dynamics (QCD) with two flavors of Wilson fermions. We compare 2nd, 4th and 6th order…
Lattice QCD should allow a derivation of the $\Delta I=1/2$ rule from first principles, but numerical calculations to date have been plagued by a variety of problems. After a brief review of these problems, we present several new methods…
We present the first results for lattice QCD at finite temperature $T$ and chemical potential $\mu$ with four flavors of Wilson quarks. The calculations are performed using the imaginary chemical potential method at $\kappa=0$, 0.001, 0.15,…
Developing efficient MCMC algorithms is indispensable in Bayesian inference. In parallel tempering, multiple interacting MCMC chains run to more efficiently explore the state space and improve performance. The multiple chains advance…
We employ the Complex Langevin method for simulation of complex-valued actions. First, we show how to test for convergence of the method by explicitely computing boundary terms and demonstrate this in a model. Then we investigate the…
We present a parallelizable SSOR preconditioning scheme for Krylov subspace iterative solvers which proves to be efficient in lattice QCD applications involving Wilson fermions. Our preconditioner is based on a locally lexicographic…
We study the opportunities for parallelism for the simulation of surface reactions. We introduce the concept of a partition and we give new simulation methods based on Cellular Automaton using partitions. We elaborate on the advantages and…
We show that it is possible to improve the chiral behaviour and the approach to the continuum limit of correlation functions in lattice QCD with Wilson fermions by taking arithmetic averages of correlators computed in theories regularized…
The investigation of the decorrelation efficiency of the HMC algorithm with respect to vacuum topology is a prerequisite for trustworthy full QCD simulations, in particular for the computation of topology sensitive quantities. We…
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 present two improvements to our previous dynamical overlap HMC algorithm. We introduce a new method of differentiating the eigenvectors of the Kernel operator, which removes an instability in the fermionic force. Secondly, by simulating…
Recent years have seen a rise in the application of machine learning techniques to aid the simulation of hard-to-sample systems that cannot be studied using traditional methods. Despite the introduction of many different architectures and…
The variational method is used widely for determining eigenstates of the QCD hamiltonian for actions with a conventional transfer matrix, e.g., actions with improved Wilson fermions. An alternative lattice fermion formalism, staggered…
We propose and use a novel, hybrid Monte Carlo algorithm that combines configurational bias particle swaps with parallel tempering. We use this new method to simulate a standard model of a glass forming binary mixture above and below the…
We discuss multi-dimensional generalizations of multicanonical algorithm, simulated tempering, and replica-exchange method. We generalize the original potential energy function $E_0$ by adding any physical quantity $V$ of interest as a new…
We discuss the adaptation of the Hybrid Monte Carlo algorithm to overlap fermions. We derive a method which can be used to account for the delta function in the fermionic force caused by the differential of the sign function. We discuss the…
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on…
Smearing the gauge links of dynamical configurations removes small scale unphysical vacuum fluctuations und thus improves the chiral properties of lattice fermions. We present a new algorithm for the simulation of dynamical fermions coupled…
The standard kinetic Monte Carlo algorithm is an extremely efficient method to carry out serial simulations of dynamical processes such as thin-film growth. However, in some cases it is necessary to study systems over extended time and…
JLQCD collaboration recently started the $N_f=3$ QCD simulations with the $O(a)$-improved Wilson fermion action employing an exact fermion algorithm developed for odd number of quark flavors. It is found that this theory has an unexpected…