Related papers: Efficient calculation of imaginary time displaced …
The eigenvalue-function pair of the 3D Schr\"odinger equation can be efficiently computed by use of high order, imaginary time propagators. Due to the diffusion character of the kinetic energy operator in imaginary time, algorithms…
We outline how auxiliary-field quantum Monte Carlo (AFQMC) can leverage graphical processing units (GPUs) to accelerate the simulation of solid state sytems. By exploiting conservation of crystal momentum in the one- and two-electron…
Numerical simulations of models and theories that describe complex systems such as spin glasses are becoming increasingly important. Beyond fundamental research, these computational methods also find practical applications in fields like…
We propose a method for computing n-time correlation functions of arbitrary spinorial, fermionic, and bosonic operators, consisting of an efficient quantum algorithm that encodes these correlations in an initially added ancillary qubit for…
Theoretical analysis typically involves imaginary-time correlation functions. Inferring real-time dynamical response functions from this information is notoriously difficult. However, as we articulate here, it is straightforward to compute…
We consider Monte Carlo simulations of classical spin models of statistical mechanics using the massively parallel architecture provided by graphics processing units (GPUs). We discuss simulations of models with discrete and continuous…
A number of applications require the computation of the trace of a matrix that is implicitly available through a function. A common example of a function is the inverse of a large, sparse matrix, which is the focus of this paper. When the…
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…
A generalization to the quantum case of a recently introduced algorithm (Y. Tomita and Y. Okabe, Phys. Rev. Lett. {\bf 86}, 572 (2001)) for the determination of the critical temperature of classical spin models is proposed. We describe a…
Monte Carlo simulations are methods for simulating statistical systems. The aim is to generate a representative ensemble of configurations to access thermodynamical quantities without the need to solve the system analytically or to perform…
We analyze a new Monte Carlo method which uses transition matrix in the space of energy. This method gives an efficient reweighting technique. The associated artificial dynamics is a constrained random walk in energy, producing the result…
In this work we compare numerically exact Quantum Monte Carlo (QMC) calculations and Green function theory (GFT) calculations of thin ferromagnetic films including second order anisotropies. Thereby we concentrate on easy plane systems,…
This article was motivated by the desire to improve Markov chain Monte Carlo methods for spatial survival models in which the locations of individuals in space are known. For a dataset comprising information on n individuals, standard…
Whenever real particle production occurs in quantum field theory, the imaginary part of the Hadamard Elementary function $G^{(1)}$ is non-vanishing. A method is presented whereby the imaginary part of $G^{(1)}$ may be calculated for a…
A general method for computing kinetic isotope effects is described. The method uses the quantum-instanton approximation and is based on the thermodynamic integration with respect to the mass of the isotopes and on the path-integral…
The time dependent quantum Monte Carlo method for fermions is introduced and applied for calculation of entanglement of electrons in one-dimensional quantum dots with several spin-polarized and spin-compensated electron configurations. The…
We develop, implement and test a set of algorithms for estimating N-point correlation functions from pixelized sky maps. These algorithms are slow, in the sense that they do not break the O(N_pix^N) barrier, and yet, they are fast enough…
A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick…
We propose a variant of the Simulated Annealing method for optimization in the multivariate analysis of differentiable functions. The method uses global actualizations via the Hybrid Monte Carlo algorithm in their generalized version for…
The directed-loop quantum Monte Carlo method is generalized to the case of retarded interactions. Using the path integral, fermion-boson or spin-boson models are mapped to actions with retarded interactions by analytically integrating out…