Related papers: Low-lying Wilson Dirac operator eigenvector mixing…
We investigate the QCD phase diagram in the strong-coupling lattice QCD with fluctuation and $1/g^2$ effects by using the auxiliary field Monte-Carlo simulations. The complex phase of the Fermion determinant at finite chemical potential is…
As a prerequisite to dynamical fermion simulations a detailed study of optimal parameters and scaling behavior is conducted for the quenched Schr\"odinger functional at fixed renormalized coupling. We compare standard hybrid overrelaxation…
We describe a controllable and unbiased strong-coupling diagrammatic Monte Carlo technique that is applicable to a wide range of fermionic systems and spin models. Unlike previous strong coupling methods that generally rely on the…
The improvement of simulations of QCD with dynamical Wilson fermions by combining the Hybrid Monte Carlo algorithm with parallel tempering is studied. As an indicator for decorrelation the topological charge is used.
The system-level dynamics of multivalent biomolecular interactions can be simulated using a rule-based kinetic Monte Carlo method in which a rejection sampling strategy is used to generate reaction events. This method becomes inefficient…
Molecular dynamics algorithms are subject to some amount of error dependent on the size of the time step that is used. This error can be corrected by periodically updating the system with a Metropolis criteria, where the integration step is…
The low-lying eigenvalues of a (sparse) hermitian matrix can be computed with controlled numerical errors by a conjugate gradient (CG) method. This CG algorithm is accelerated by alternating it with exact diagonalisations in the subspace…
We perform hybrid Monte-Carlo simulation of $N_f=2+1+1 $ lattice QCD with domain-wall quarks at the physical point. The simulation is carried out on the $ L^3 \times T = 64^3 \times 64 $ lattice with lattice spacing $a \sim 0.064 $ fm ($ L…
We define a sparse hermitian lattice Dirac matrix, $H$, coupling $2n+1$ Dirac fermions. When $2n$ fermions are integrated out the induced action for the last fermion is a rational approximation to the hermitian overlap Dirac operator. We…
We show how the integrators used for the molecular dynamics step of the Hybrid Monte Carlo algorithm can be further improved. These integrators not only approximately conserve some Hamiltonian $H$ but conserve exactly a nearby shadow…
The low energy physics of both graphene and surface states of three-dimensional topological insulators is described by gapless Dirac fermions with linear dispersion. In this work, we predict the emergence of a "heavy" Dirac fermion in a…
We introduce a factorization of the fermion determinant in lattice QCD with Wilson-type fermions that leads to a bosonic action which is local in the block fields. The interaction among gauge fields on distant blocks is mediated by…
The study of ultracold optically trapped atoms has opened new vistas in the physics of correlated quantum systems. Much attention has now turned to mixtures of bosonic and fermionic atoms. A central puzzle is the disagreement between the…
We study interacting Majorana fermions in two dimensions as a low-energy effective model of a vortex lattice in two-dimensional time-reversal-invariant topological superconductors. For that purpose, we implement ab-initio quantum Monte…
Using exact continuous quantum Monte Carlo techniques, we study the zero and finite temperature properties of a system of harmonically trapped one dimensional spin 1/2 fermions with short range interactions. Motivated by experimental…
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…
We present a polynomial Hybrid Monte Carlo (PHMC) algorithm as an exact simulation algorithm with dynamical Kogut-Susskind fermions. The algorithm uses a Hermitian polynomial approximation for the fractional power of the KS fermion matrix.…
In application of the complex Langevin method to QCD at high density and low temperature, the singular-drift problem occurs due to the appearance of near-zero eigenvalues of the Dirac operator. In order to avoid this problem, we proposed to…
Large N gauge theories with adjoint matter can be numerically studied using lattice techniques. Eguchi-Kawai reductions holds for this theory and one can reduce the lattice model to a single site. Hybrid Monte Carlo algorithm can be used to…
We discuss the lattice formulation of gauge theories with fermions in arbitrary representations of the color group, and present in detail the implementation of the HMC/RHMC algorithm for simulating dynamical fermions. We discuss the…