Related papers: Direct simulation Monte Carlo schemes for Coulomb …
A gyrokinetic Coulomb collision operator is derived, which is particularly useful to describe the plasma dynamics at the periphery region of magnetic confinement fusion devices. The derived operator is able to describe collisions occurring…
The linearized collision operator of the Boltzmann equation can in a natural way be written as a sum of a positive multiplication operator, the collision frequency, and an integral operator. Compactness of the integral operator for…
In this study, we consider a numerical implementation of the nonlinear Rosenbluth-Trubnikov collision operator for particle simulations in plasma physics in the framework of the finite element method (FEM). The relevant particle evolution…
We present novel Monte Carlo methods for treating the interacting shell model that allow exact calculations much larger than those heretofore possible. The two-body interaction is linearized by an auxiliary field; Monte Carlo evaluation of…
Recently, a Quantum Monte Carlo method alternative to the Path Integral Monte Carlo method was developed for the numerical solution of the N-boson problem; it is based on the stochastic evolution of classical fields. Here we apply it to…
Quantum nanosystems involve the coupled dynamics of fermions or bosons across multiple scales in space and time. Examples include quantum dots, superconducting or magnetic nanoparticles, molecular wires, and graphene nanoribbons. The number…
This chapter is devoted to the computation of equilibrium (thermodynamic) properties of quantum systems. In particular, we will be interested in the situation where the interaction between particles is so strong that it cannot be treated as…
This article describes methods for the deterministic simulation of the collisional Boltzmann equation. It presumes that the transport and collision parts of the equation are to be simulated separately in the time domain. Time stepping…
We created an efficient algorithm suitable for graphics processing units (GPUs) to perform Monte Carlo simulations of a subset of reaction-diffusion models. The algorithm uses techniques that are specific to GPU programming, and combines…
We consider M systems (each an electron in a long square cylinder) uniformly arranged on a ring and with Coulomb interactions. Exact straightforward numerical time-dependent perturbation calculation of a single N-level ($\lesssim 7$)…
A Monte Carlo method for the collisional guiding-center Fokker-Planck kinetic equation is derived to include the effects of background magnetic-field nonuniformity. It is shown that, in the limit of a homogeneous magnetic field, the…
We present a promising coarse-graining strategy for linking micro- and mesoscales of soft matter systems. The approach is based on effective pairwise interaction potentials obtained from detailed atomistic molecular dynamics (MD)…
Polarizable particle systems, including charged colloids, polarizable ions, biomolecular assemblies, and soft nanomaterials, can exhibit contact electrostatic interactions that depart strongly from Coulomb behavior when dielectric mismatch…
The study of ultracold optically trapped atoms has opened new vistas in the physics of correlated quantum systems. Much attention has now turned to mixtures of bosonic and fermionic atoms. A central puzzle is the disagreement between the…
This paper introduces a class of Monte Carlo algorithms which are based upon the simulation of a Markov process whose quasi-stationary distribution coincides with a distribution of interest. This differs fundamentally from, say, current…
We study the dynamical evolution of globular clusters containing populations of primordial binaries, using our newly updated Monte Carlo cluster evolution code with the inclusion of direct integration of binary scattering interactions. We…
We propose a simplified version of local molecular field (LMF) theory to treat Coulomb interactions in simulations of ionic fluids. LMF theory relies on splitting the Coulomb potential into a short-ranged part that combines with other…
We develop a random batch Ewald (RBE) method for molecular dynamics simulations of particle systems with long-range Coulomb interactions, which achieves an $O(N)$ complexity in each step of simulating the $N$-body systems. The RBE method is…
Monte Carlo (MC) methods are important computational tools for molecular structure optimizations and predictions. When solvent effects are explicitly considered, MC methods become very expensive due to the large degree of freedom associated…
In this paper we present a dynamical Monte Carlo algorithm which is applicable to systems satisfying a clustering condition: during the dynamical evolution the system is mostly trapped in deep local minima (as happens in glasses, pinning…