Related papers: Schur complement solver for Quantum Monte-Carlo si…
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 the splitting of the pseudo-fermion action into two parts. We test…
In this paper, a new progressive mesh algorithm is introduced in order to perform fast physical simulations by the use of a lattice Boltzmann method (LBM) on a single-node multi-GPU architecture. This algorithm is able to mesh automatically…
We investigate the QCD phase diagram in the strong coupling limit by using a newly developed auxiliary field Monte-Carlo (AFMC) method. Starting from an effective action in the leading order of the 1/g^2 and 1/d expansion with one species…
Perturbative expansions of several small Wilson loops are computed through next-to-next-to-leading order in unquenched lattice QCD, from Monte Carlo simulations at weak couplings. This approach provides a much simpler alternative to…
We discuss designer Hamiltonians---lattice models tailored to be free from sign problems ("de-signed") when simulated with quantum Monte Carlo methods but which still host complex many-body states and quantum phase transitions of interest…
We present the first GPU-based conjugate gradient (CG) solver for lattice QCD with domain-wall fermions (DWF). It is well-known that CG is the most time-consuming part in the Hybrid Monte Carlo simulation of unquenched lattice QCD, which…
We investigate the quench of half-filled 1D and 2D fermionic Hubbard models to models without Coulomb interaction. Since the time propagation is gaussian we can use a variety of time-dependent quantum Monte Carlo methods to tackle this…
We consider an efficient scheme to simulate fermionic Hubbard models with nonlocal density-density interactions in two dimensions, based on bond-centered auxiliary-field quantum Monte Carlo. The simulations are shown to be sign-problem free…
We present the results of an effort to accelerate a Rational Hybrid Monte Carlo (RHMC) program for lattice quantum chromodynamics (QCD) simulation for 2 flavours of staggered fermions on multiple Kepler K20X GPUs distributed on different…
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…
Fourier solvers have become efficient tools to establish structure-property relations in heterogeneous materials. Introduced as an alternative to the Finite Element (FE) method, they are based on fixed-point solutions of the…
In this work, results are presented of Hybrid-Monte-Carlo simulations of the tight-binding Hamiltonian of graphene, coupled to an instantaneous long-range two-body potential which is modeled by a Hubbard-Stratonovich auxiliary field. We…
Stimulated by the successful descriptions of strongly correlated electron systems by fractionalized fermions, correspondence between interacting fermions and non-interacting multi-component fermions is formulated in examples of the Hubbard…
We study an exactly solvable quantum field theory (QFT) model describing interacting fermions in 2+1 dimensions. This model is motivated by physical arguments suggesting that it provides an effective description of spinless fermions on a…
Lattice Quantum Chromodynamics (Lattice QCD) is a quantum field theory on a finite discretized space-time box so as to numerically compute the dynamics of quarks and gluons to explore the nature of subatomic world. Solving the equation of…
In order to develop fast inversion algorithms we have used overlap solvers in two dimensions. Lattice QED theory with U(1) group symmetry in two dimensional space-times dimensions has always been a testing ground for algorithms. By the…
The absence of negative sign problem in quantum Monte Carlo simulations of spin and fermion systems has different origins. World-line based algorithms for spins require positivity of matrix elements whereas auxiliary field approaches for…
We present a method of simulating lattice QCD at nonzero chemical potential in the chiral limit. By adding a weak four-fermi interaction to the standard staggered fermion SU(3) QCD action, we produce an algorithm in which the limit of…
This work introduces CLIP, a CUDA-accelerated phase-field lattice Boltzmann framework for simulating immiscible two-phase flows with high density and viscosity ratios in both two- and three-dimensional domains. By leveraging GPU…
The solution for non-linear, complex partial differential Equations (PDEs) is achieved through numerical approximations, which yield a linear system of equations. This approach is prevalent in Computational Fluid Dynamics (CFD), but it…