Related papers: Diagrammatic Monte Carlo algorithm for the resonan…
A general-purpose, self-adapting Monte Carlo (MC) algorithm implemented in the program {\tt Foam} is described. The high efficiency of the MC, that is small maximum weight or variance of the MC weight is achieved by means of dividing the…
We present a numerical method and computer code to calculate the radiative transfer and excitation of molecular lines. Formulating the Monte Carlo method from the viewpoint of cells rather than photons allows us to separate local and…
The reptation Monte Carlo algorithm is a simple, physically motivated and efficient method for equilibrating semi-dilute solutions of linear polymers. Here we propose two simple generalizations for the analogue {\it Amoeba} algorithm for…
We describe a Monte Carlo procedure for the simulation of dynamically triangulate random surfaces with a boundary (topology of a disk). The algorithm keeps the total number of triangles fixed, while the length of the boundary is allowed to…
We consider population-imbalanced two-component Fermi gases under external harmonic confinement interacting through short-range two-body potentials with diverging s-wave scattering length. Using the fixed-node diffusion Monte Carlo method,…
The dynamical responses of Ising metamagnet (layered antiferromagnet) in the presence of a sinusoidally oscillating magnetic field are studied by Monte Carlo simulation. The time average staggered magnetisation plays the role of dynamic…
This manuscript describes an adjoint/reverse Monte Carlo method to calculate the flux of charged plasma particles to the wall of e.g. a tokamak. Two applications are described: a fusion product activation probe and a neutral beam injection…
We present a simulation algorithm for dynamical fermions that combines the multiboson technique with the Hybrid Monte Carlo algorithm. We find that the algorithm gives a substantial gain over the standard methods in practical simulations.…
We address the calculation of dynamical correlation functions for many fermion systems at zero temperature, using the auxiliary-field quantum Monte Carlo method. The two-dimensional Hubbard hamiltonian is used as a model system. Although…
The behavior of a Lattice Monte Carlo algorithm (if it is designed correctly) must approach that of the continuum system that it is designed to simulate as the time step and the mesh step tend to zero. However, we show for an algorithm for…
The Hubbard model and extended Hubbard model on the honeycomb lattice can be seen as prototype models of single layer graphene placed in a high dielectric constant environment that screens the Coulomb interaction. Taking advantage of the…
The short-range modes of the fermionic determinant can be absorbed in the gauge action using the loop expansion. The coefficients of this expansion and the zeroes of the polynomial approximating the remainder can be optimized by a simple,…
A recent reformulation [1] of the problem of Coulomb gases in the presence of a dynamical dielectric medium showed that finite temperature simulations of such systems can be accomplished on the basis of completely local Hamiltonians on a…
The Monte Carlo method is a thriving and mathematically beautiful numerical technique used extensively, nowadays, to deal with many demanding problems in diverse fields. Here, we present an iterative Monte Carlo algorithm to work out very…
We explore the possibilities of using the fermionic functional renormalization group to compute the phase diagram of systems with competing instabilities. In order to overcome the ubiquituous divergences encountered in RG flows, we propose…
We elaborate on the methodology to simulate bulk systems in the absence of time-reversal symmetry by the phase-fixed path-integral Monte Carlo method under (possibly twisted) periodic boundary conditions. Such systems include…
We demonstrate that summation of connected diagrams to high order starting from a BCS hamiltonian is a viable generic unbiased approach for strongly correlated fermions in superconducting or superfluid phases. For the 3D attractive Hubbard…
We develop a novel Monte Carlo algorithm for the vector consisting of the supremum, the time at which the supremum is attained and the position at a given (constant) time of an exponentially tempered L\'evy process. The algorithm, based on…
We study the Fermi gas at unitarity and at T=0 by assuming that, at high polarizations, it is a normal Fermi liquid composed of weakly interacting quasiparticles associated with the minority spin atoms. With a quantum Monte Carlo approach…
We adapt the Quantum Monte Carlo method to the cascaded formalism of quantum optics, allowing us to simulate the emission of photons of known energy. Statistical processing of the photon clicks thus collected agrees with the theory of…