Related papers: GIBS: A grand-canonical Monte Carlo simulation pro…
Molecular simulations in an open environment, involving ion exchange, are necessary to study various systems, from biosystems to confined electrolytes. However, grand-canonical simulations are often computationally demanding in condensed…
Monte Carlo simulation provides a powerful tool for understanding and exploring thermodynamic phase equilibria in many-particle interacting systems. Among the most physically intuitive simulation methods is Gibbs ensemble Monte Carlo…
Gibbs measures, such as Coulomb gases, are popular in modelling systems of interacting particles. Recently, we proposed to use Gibbs measures as randomized numerical integration algorithms with respect to a target measure $\pi$ on $\mathbb…
We present an open-source code for the simulation of electron and ion transport for user-defined gas mixtures with static uniform electric and magnetic fields. The program provides microscopic interaction simulation and is interfaced with…
A Grand-canonical Monte-Carlo simulation method extended to simulate a mixture of salts is presented. Due to charge neutrality requirement of electrolyte solutions, ions must be added to or removed from the system in groups. This leads to…
``Glauber'' models are used to calculate geometric quantities in the initial state of heavy ion collisions, such as impact parameter, number of participating nucleons and initial eccentricity. The four RHIC experiments have different…
Artifacts arise in the simulations of electrolytes using periodic boundary conditions (PBC). We show the origin of these artifacts are the periodic image charges and the constraint of charge neutrality inside the simulation box, both of…
In this work we present an efficient implementation of Canonical Monte Carlo simulation for Coulomb many body systems on graphics processing units (GPU). Our method takes advantage of the GPU Single Instruction, Multiple Data (SIMD)…
A physical model based on a Monte-Carlo approach is proposed to calculate the ionization dynam- ics of warm dense matters (WDM) within particle-in-cell simulations, and where the impact (col- lision) ionization (CI), electron-ion…
The ion production in an ECRIS is modelled using a particle-in-cell Monte-Carlo collision code in a three-dimensional geometry. Only the heavy particles (ions and atoms) are tracked, with the electron density determined from the requirement…
Using the specific model of a system of like charged ions confined between two planar like charged surfaces, we compare the predictions for the energy and density profile of four simulation methods available to treat the long range Coulomb…
This work describes a new 1D hybrid approach for modeling atmospheric pressure discharges featuring complex chemistry. In this approach electrons are described fully kinetically using Particle-In-Cell/Monte-Carlo (PIC/MCC) scheme, whereas…
We present an extended version of GLISSANDO, a Monte-Carlo generator for Glauber-like models of the initial stage of relativistic heavy-ion collisions. The increased functionality of the code incorporates a parametrization of shape of…
This chapter surveys the most standard Monte Carlo methods available for simulating from a posterior distribution associated with a mixture and conducts some experiments about the robustness of the Gibbs sampler in high dimensional Gaussian…
Fluid phase behavior of charge-stabilized colloidal suspensions is explored by applying a new variant of the Gibbs ensemble Monte Carlo simulation method to a coarse-grained one-component model with implicit microions and solvent. The…
The Glauber model is extensively applied to heavy ion collision for describing a number of interaction processes over a wide range of energies from near the Coulomb barrier to higher energies. The model gives the nucleus-nucleus interaction…
We have designed an improved algorithm that enables us to simulate the chemistry of cold dense interstellar clouds with a full gas-grain reaction network. The chemistry is treated by a unified microscopic-macroscopic Monte Carlo approach…
Ultracold atomic systems have been of great research interest in the past, with more recent attention being paid to systems of mixed species. In this work we carry out non-perturbative Path Integral Monte Carlo (PIMC) simulations of N…
The results of analytical approximations and extensive calculations based on a path integral Monte Carlo (PIMC) scheme are presented. A new (direct) PIMC method allows for a correct determination of thermodynamic properties such as energy…
The collisions between high-energy ions and solids can result in significant physical and chemical changes to the material. These effects are potentially important for better understanding the chemistry of interstellar and planetary bodies,…