Related papers: Numerical heating in particle-in-cell simulations …
The discrete time path integral Monte Carlo (PIMC) with a one-particle density matrix approximation is applied to study the quantum phase transition in the coupled double-well chain. To improve the convergence properties, the exact action…
We develop a parallel rejection algorithm to tackle the problem of low acceptance in Monte Carlo methods, and apply it to the simulation of the hopping conduction in Coulomb glasses using Graphics Processing Units, for which we also…
Two first-principles simulation techniques, path integral Monte Carlo (PIMC) and density functional molecular dynamics (DFT-MD), are applied to study hot, dense helium in the density-temperature range of 0.387 - 5.35 g/cc and 500 K -…
Monte Carlo simulation code has been developed and tested for studying the passage of charged particle beams and radiation through the crystalline matter at energies from tens of MeV up to hundreds of GeV. The developed Monte Carlo code…
We describe the Monte Carlo (MC) simulation package of the Borexino detector and discuss the agreement of its output with data. The Borexino MC 'ab initio' simulates the energy loss of particles in all detector components and generates the…
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 introduce an extension of the particle-in-cell (PIC) method that captures the Landau collisional effects in the Vlasov-Maxwell-Landau equations. The method arises from a regularisation of the variational formulation of the Landau…
An efficient Path Integral Monte Carlo procedure is proposed to simulate the behavior of quantum many-body dissipative systems described within the framework of the influence functional. Thermodynamic observables are obtained by Monte Carlo…
A numerical method to implement a linearized Coulomb collision operator in the two-weight $\delta f$ Monte Carlo method for multi-ion-species neoclassical transport simulation is developed. The conservation properties and the adjointness…
A grand canonical Monte Carlo (MC) algorithm is presented for studying the lattice gas model (LGM) of multiple protein sequence alignment, which coherently combines long-range interactions and variable-length insertions. MC simulations are…
Particle-in-Cell Monte Carlo simulations on large-scale systems play a fundamental role in understanding the complexities of plasma dynamics in fusion devices. Efficient handling and analysis of vast datasets are essential for advancing…
Ideally, binary-collision algorithms conserve kinetic momentum and energy. In practice, the finite size of collision cells and the finite difference in the particle locations affect the conservation properties. In the present work, we…
The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact $N$-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC…
Modeling the dynamics of a quantum system connected to the environment is critical for advancing our understanding of complex quantum processes, as most quantum processes in nature are affected by an environment. Modeling a macroscopic…
Monte Carlo (MC) simulations are extensively used for various purposes in modern high-energy physics (HEP) experiments. Precision measurements of established Standard Model processes or searches for new physics often require the collection…
Comptonization is the process in which photon spectrum changes due to multiple Compton scatterings in the electronic plasma. It plays an important role in the spectral formation of astrophysical X-ray and gamma-ray sources. There are…
The results of numerical simulation using a modified Monte Carlo method with a heat bath algorithm for the pseudospin model of cuprates are presented. The temperature phase diagrams are constructed for various degrees of doping and for…
The computer revolution has been driven by a sustained increase of computational speed of approximately one order of magnitude (a factor of ten) every five years since about 1950. In natural sciences this has led to a continuous increase of…
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 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…