Related papers: Efficient Methods for Handling Long-Range Forces i…
In a previous paper, we developed a new particle-in-cell method for the Vlasov-Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a…
In this paper the normal collision of spherical particles is investigated. The particle interaction is modelled in a macroscopic way using the Hertzian contact force with additional linear damping. The goal of the work is to develop an…
Continuum solvation models are becoming increasingly relevant in condensed matter simulations, allowing to characterize materials interfaces in the presence of wet electrified environments at a reduced computational cost with respect to all…
We present a computational algorithm for computing short range forces between particles. The algorithm has two distinguishing features. First, it is optimized for multi-processor computers, and will use as many processors as are available.…
We use a quantum Monte Carlo method to investigate various classes of 2D spin models with long-range interactions at low temperatures. In particular, we study a dipolar XXZ model with U(1) symmetry that appears as a hard-core boson limit of…
The Direct Simulation Monte Carlo (DSMC) method is widely employed for simulating rarefied nonequilibrium gas flows. With advances in aerospace engineering and micro/nano-scale technologies, gas flows exhibit the coexistence of rarefied and…
This paper describes an algorithm for selecting parameter values (e.g. temperature values) at which to measure equilibrium properties with Parallel Tempering Monte Carlo simulation. Simple approaches to choosing parameter values can lead to…
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)…
Kinetic plasma processes, such as magnetic reconnection, collisionless shocks, and turbulence, are fundamental to the dynamics of astrophysical and laboratory plasmas. Simulating these processes often requires particle-in-cell (PIC)…
Continuous-time random disturbances from the renewable generation pose a significant impact on power system dynamic behavior. In evaluating this impact, the disturbances must be considered as continuous-time random processes instead of…
We derived a number of numerical methods to treat biomolecular systems with multiple time scales. Based on the splitting of the operators associated with the slow-varying and fast-varying forces, new multiple time-stepping (MTS) methods are…
The approximate computation of all gravitational forces between $N$ interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than $\mathcal{O}(N)$ operations. FMM groups…
A general synthetic iterative scheme is proposed to solve the Enskog equation within a Monte Carlo framework. The method demonstrates rapid convergence by reducing intermediate Monte Carlo evolution and preserves the asymptotic-preserving…
The simulation of large nonlinear dynamical systems, including systems generated by discretization of hyperbolic partial differential equations, can be computationally demanding. Such systems are important in both fluid and kinetic…
Unbiased stochastic sampling of the one- and two-body reduced density matrices is achieved in full configuration interaction quantum Monte Carlo with the introduction of a second, "replica" ensemble of walkers, whose population evolves in…
Gauge theories establish the standard model of particle physics, and lattice gauge theory (LGT) calculations employing Markov Chain Monte Carlo (MCMC) methods have been pivotal in our understanding of fundamental interactions. The present…
This paper gives an overview of the particle transport theory essentials, the basics of particle-material interaction simulation, physical quantities needed to simulate particle transport and interactions in materials, Monte Carlo…
A $\theta$ term, which couples to topological charge, is added to the two-dimensional lattice CP^3 model and U(1) gauge theory. Monte Carlo simulations are performed and compared to strong-coupling character expansions. In certain…
We propose two new Bayesian smoothing methods for general state-space models with unknown parameters. The first approach is based on the particle learning and smoothing algorithm, but with an adjustment in the backward resampling weights.…
We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo…