Related papers: A Binary Collision Method for Screened Coulomb Col…
An expression for the effective Coulomb logarithm in one-component-plasma is proposed, which allows to extend the applicability of the classical formula for the self-diffusion coefficient to the strongly coupled regime. The proposed…
When a Monte Carlo algorithm is used to evaluate a physical observable A, it is possible to slightly modify the algorithm so that it evaluates simultaneously A and the derivatives $\partial$ $\varsigma$ A of A with respect to each…
We present a lattice Monte Carlo algorithm based on the one originally proposed by Maggs and Rossetto for simulating electrostatic interactions in inhomogeneous dielectric media. The original algorithm is known to produce attractive…
In this paper we develop a direct simulation Monte Carlo (DSMC) method for simulating highly nonequilibrium dynamics of nearly degenerate ultra-cold gases. We show that our method can simulate the high-energy collision of two thermal clouds…
Quasi-monochromatic, high energy and highly polarized $\gamma$-ray beam sources based on Compton scattering of laser photons (LCS) on relativistic electrons have developed for the last few decades as established instruments for nuclear…
Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop a novel approach to reduce this…
We develop path-integral Monte Carlo simulations for a parabolic two-dimensional (2D) quantum dot containing $N$ interacting electrons in the presence of Dresselhaus and/or Rashba spin-orbit couplings. Our method solves in a natural way the…
In this work we propose a Hybrid method with Deviational Particles (HDP) for a plasma modeled by the inhomogeneous Vlasov-Poisson-Landau system. We split the distribution into a Maxwellian part evolved by a grid based fluid solver and a…
We present a method which extends Monte Carlo studies to situations that require a large dynamic range in particle number. The underlying idea is that, in order to calculate the collisional evolution of a system, some particle interactions…
Generalising the work of Lenard and Bernstein, we introduce a new, fully relativistic model to describe collisional plasmas. Like the Fokker-Planck operator, this equation represents velocity diffusion and conserves particle number.…
A model is presented for the solution of electrokinetic phenomena of colloidal suspensions in fluid mixtures. We solve the discrete Boltzmann equation with a BGK collision operator using the lattice Boltzmann method to simulate binary fluid…
We apply the diagrammatic Monte Carlo approach to three-dimensional Fermi-polaron systems with mass-imbalance, where an impurity interacts resonantly with a noninteracting Fermi sea whose atoms have a different mass. This method allows to…
We carry out a sign-problem-free quantum Monte Carlo calculation of a bilayer model with a repulsive intra-layer Hubbard interaction and a ferromagnetic inter-layer interaction. The latter breaks the global $SU(2)$ spin rotational symmetry…
We use a Monte Carlo approach to study hadron azimuthal angular correlations in high energy proton-proton and central nucleus-nucleus collisions at the BNL Relativistic Heavy Ion Collider (RHIC) energies at mid-rapidity. We build a hadron…
We report on Monte Carlo studies of the hard X-ray polarimeter X-Calibur. The polarimeter will be used in the focal plane of a grazing incidence hard X-ray telescope. It combines a low-Z Compton scatterer with a high-Z Cadmium Zinc…
The differential formulation of the Landau-Fokker-Planck collision integral is developed for the case of relativistic electromagnetic interactions.
Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically…
We present a Monte Carlo code dedicated to the simulation of bremsstrahlung arising in collisions of polarized electrons with thin target foils. The program consists of an electron transport algorithm taking into account elastic…
Kinetic descriptions of runaway electrons (RE) are usually based on Fokker-Planck models that determine the probability distribution function (PDF) of RE in 2-dimensional momentum space. Despite of the simplification involved, the…
We use path-integral Monte Carlo (PIMC) to study the effects of adding a long-range repulsive Coulomb interaction to the usual Van der Waals interaction between two atoms of a submonolayer quantum film such as helium on graphite or a pure…