Related papers: The RHMC Algorithm for 2 Flavours of Dynamical Sta…
Nested integration of the form $\int f\left(\int g(\bs{y},\bs{x})\di{}\bs{x}\right)\di{}\bs{y}$, characterized by an outer integral connected to an inner integral through a nonlinear function $f$, is a challenging problem in various fields,…
A general formalism was introduced for reorganization of QCD evolution equations and derivation of hierarchical solution to DGLAP equation. The hierarchical solution separates two types of parton emissions: the flavour changing emissions…
We present an update on the MILC Collaboration's light hadron spectrum calculation with two flavors of dynamical, staggered quarks. We present extrapolations of the nucleon to rho mass ratio to the continuum limit for fixed values of the pi…
We compute the condensate in QCD with two flavors of dynamical fermions using numerical simulation. The simulations use overlap fermions, and the condensate is extracted by fitting the distribution of low lying eigenvalues of the Dirac…
Diagrammatic Monte Carlo (DiagMC) is a numeric technique that allows one to calculate quantities specified in terms of diagrammatic expansions, the latter being a standard tool of many-body quantum statistics. The sign problem that is…
For thermodynamics studies it is desirable to simulate two degenerate flavors and retain at least a remnant of the chiral symmetry. Staggered fermions can achieve this at the cost of rooting the determinant. Rooting can be avoided using…
We report on recent progress in employing the Highly Improved Staggered Quark (HISQ) action introduced by the HPQCD/UKQCD collaboration in simulations with dynamical fermions. The HISQ action is an order $a^2$ Symanzik-improved action with…
We introduce a new Monte Carlo method for pure gauge theories. It is not intended for use with dynamical fermions. It belongs to the class of Local Hybrid Monte Carlo (LHMC) algorithms, which make use of the locality of the action by…
To obtain a light mode in two-dimensional staggered fermions, we introduce four new local operators keeping the rotational invariance for a staggered Dirac operator. To split masses of tastes, three cases are considered. The mass matrix and…
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…
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 present a random matrix theory (RMT) for the staggered lattice QCD Dirac operator. The staggered RMT is equivalent to the zero-momentum limit of the staggered chiral Lagrangian and includes all taste breaking terms at their leading…
Many recent Markov chain Monte Carlo (MCMC) samplers leverage continuous dynamics to define a transition kernel that efficiently explores a target distribution. In tandem, a focus has been on devising scalable variants that subsample the…
Modern implementations of Hamiltonian Monte Carlo and related MCMC algorithms support sampling of probability functions that embed numerical root-finding algorithms, thereby allowing fitting of statistical models involving analytically…
Hamiltonian Monte Carlo (HMC) has emerged as a powerful Markov Chain Monte Carlo (MCMC) method to sample from complex continuous distributions. However, a fundamental limitation of HMC is that it can not be applied to distributions with…
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…
We present preliminary results from UKQCD simulations at light quark masses using two flavours of non-pertubatively improved Wilson fermions. We report on the performance of the standard HMC algorithm at these quark masses where m_pi/m_rho…
The dual-fermion approach provides a formally exact prescription for calculating properties of a correlated electron system in terms of a diagrammatic expansion around dynamical mean-field theory (DMFT). Most practical implementations,…
We investigate the combination of a two-level sampling algorithm with distillation techniques to compute disconnected fermionic correlation functions. The method relies on a factorization of the quark propagator into domain-local…
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…