Related papers: A Finite Baryon Density Algorithm
Ab-initio studies of strongly interacting bosonic and fermionic systems is greatly facilitated by efficient Monte Carlo algorithms. This article emphasizes this requirement, and outlines the ideas behind the construction of the cluster…
We present a new lattice Monte Carlo approach developed for studying large numbers of strongly interacting nonrelativistic fermions, and apply it to a dilute gas of unitary fermions confined to a harmonic trap. Our lattice action is highly…
We develop and test a spectral-density analysis method, based on the introduction of smeared energy kernels, to extract physical information from two-point correlation functions computed numerically in lattice field theory. We apply it to a…
The performance of the Monte Carlo sampling methods relies on the crucial choice of a proposal density. The notion of optimality is fundamental to design suitable adaptive procedures of the proposal density within Monte Carlo schemes. This…
Typical fermion algorithms require the computation (or sampling) of the fermion determinant. We focus instead on cluster algorithms which do not involve the determinant and involve a more physically relevant sampling of the configuration…
We study improved variants of L\"uscher's algorithm, including the cases of staggered fermions and finite density.
Three possibilities to speed up the Hybrid Monte Carlo algorithm are investigated. Changing the step-size adaptively brings no practical gain. On the other hand, substantial improvements result from using an approximate Hamiltonian or a…
We present a new quantum Monte Carlo algorithm suitable for generically complex problems, such as systems coupled to external magnetic fields or anyons in two spatial dimensions. We find that the choice of gauge plays a nontrivial role, and…
QCD at finite density presents specific challenges to lattice gauge theory. Nonetheless, a region of the QCD phase diagram up to moderately large baryon chemical potentials has been successfully explored on the lattice and new results and…
We present a quantum algorithm to compute the logarithm of the determinant of the fermion matrix, assuming access to a classical lattice gauge field configuration. The algorithm uses the quantum eigenvalue transform, and quantum mean…
We study fluctuations in the canonical ensemble, where the net baryon number is exactly conserved. The focus is on cumulants and factorial cumulants linked to the baryon and antibaryon multiplicities and their sum or difference in full…
Bayesian models have many desirable properties, most notable is their ability to generalize from limited data and to properly estimate the uncertainty in their predictions. However, these benefits come at a steep computational cost as…
We present a massively parallel quantum Monte Carlo based implementation of real-space dynamical mean-field theory for general inhomogeneous correlated fermionic lattice systems. As a first application, we study magnetic order in a binary…
We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the…
In this manuscript, we consider a finite nonparametric mixture model with non-independent marginal density functions. Dependence between the marginal densities is modeled using a copula device. Until recently, no deterministic algorithms…
We demonstrate the applicability of a recently proposed multiscale thermalization algorithm to two-color quantum chromodynamics (QCD) with two mass-degenerate fermion flavors. The algorithm involves refining an ensemble of gauge…
The state estimation problem for nonlinear systems with stochastic uncertainties can be formulated in the Bayesian framework, where the objective is to replace the state completely by its probability density function. Without the…
We derive a lattice approximation for a class of equilibrium quantum statistics describing the behaviour of any combination and number of bosonic and fermionic particles with any sufficiently binding potential. We then develop an intuitive…
Charge neutralisation procedure based on a Markov chain Monte Carlo algorithm is applied to a system of generated hadrons. The algorithm changes the charge of a randomly picked particle by shifting it within its isomultiplet. For baryons…
We present extensive new \emph{ab initio} path integral Monte Carlo (PIMC) results for a variety of structural properties of warm dense hydrogen and beryllium. To deal with the fermion sign problem -- an exponential computational bottleneck…