Related papers: Polynomial Filtered HMC -- an algorithm for lattic…
Accurately calculating the mass of flavor-singlet meson states from numerical lattice simulations is an important milestone for lattice QCD. Careful measurement of the full pseudoscalar flavor-singlet propagator is also a crucial step in…
In quantum chromodynamics (QCD) at nonzero chemical potential, the eigenvalues of the Dirac operator are scattered in the complex plane. Can the fluctuation properties of the Dirac spectrum be described by universal predictions of…
Smearing the gauge links of dynamical configurations removes small scale unphysical vacuum fluctuations and thus improves the chiral properties of lattice fermions. Recently we proposed the hypercubic smearing (HYP) that improves the flavor…
We report on first results of an ongoing effort to simulate lattice QCD with two degenerate flavours of quarks by means of the twisted mass formulation tuned to maximal twist. By utilising recent improvements of the HMC algorithm,…
Hamiltonian Monte Carlo (HMC) has become routinely used for sampling from posterior distributions. Its extension Riemann manifold HMC (RMHMC) modifies the proposal kernel through distortion of local distances by a Riemannian metric. The…
The R algorithm is widely used for simulating two flavours of dynamical staggered fermions. We give a simple proof that the algorithm converges to the desired probability distribution to within O(dt^2) errors, but show that the relevant…
This paper introduces a novel ansatz-based technique for solution of the Hubbard model over two length scales. Short range correlations are treated exactly using a dynamical cluster approximation QMC simulation, while longer-length-scale…
Recently, the Hamilton Monte Carlo (HMC) has become widespread as one of the more reliable approaches to efficient sample generation processes. However, HMC is difficult to sample in a multimodal posterior distribution because the HMC chain…
The extended Hubbard model on a two-dimensional lattice captures key physical phenomena, but is challenging to simulate due to the presence of long-range interactions. In this work, we present an efficient quantum algorithm for simulating…
One of the tasks in color image processing and computer vision is to recover clean data from partial observations corrupted by noise. To this end, robust quaternion matrix completion (QMC) has recently attracted more attention and shown its…
Lattice Quantum Chromodynamics (QCD) is an approach used by theoretical physicists to model the strong nuclear force. This works at the sub-nuclear scale to bind quarks together into hadrons including the proton and neutron. One of the long…
Least-squares optimized polynomials are discussed which are needed in the two-step multi-bosonic algorithm for Monte Carlo simulations of quantum field theories with fermions. A recurrence scheme for the calculation of necessary…
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…
In recent papers a microscopic model for the SM Higgs mechanism has been proposed, and an idea how to determine the 24 quark and lepton masses of all 3 generations has emerged in that framework. This idea is worked out in detail here by…
Lattice simulations of Quantum Chromodynamics (QCD) enable one to calculate the low-energy properties of the strong interaction among quarks and gluons based on the first principle. The most time-consuming part of the numerical simulations…
Polynomial filtering can provide a highly effective means of computing all eigenvalues of a real symmetric (or complex Hermitian) matrix that are located in a given interval, anywhere in the spectrum. This paper describes a technique for…
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not…
We introduce a numerical method for generating the approximating polynomials used in fermionic calculations with smeared link actions. We investigate the stability of the algorithm and determine the optimal weight function and the optimal…
We propose a versatile and efficient algorithmic framework for optimizing fermion-to-qubit mappings by generalizing the idea of randomized block coordinate descent. Our greedy approach, termed Randomized Subsystem Descent, iteratively…
We have presented some practical consequences on the molecular-dynamics simulations arising from the numerical algorithm published recently in paper Int. J. Mod. Phys. C 16, 413 (2005). The algorithm is not a finite-difference method and…