Related papers: Random batch particle methods for the homogeneous …
The multiple range random walk algorithm recently proposed by Wang and Landau [Phys. Rev. Lett. 86, 2050 (2001)] is adapted to the computation of free energy profiles for molecular systems along reaction coordinates. More generally, we show…
Hamiltonian Monte Carlo (HMC) algorithms which combine numerical approximation of Hamiltonian dynamics on finite intervals with stochastic refreshment and Metropolis correction are popular sampling schemes, but it is known that they may…
Particle-in-cell methods with stochastic collision models are commonly used to simulate collisional plasma dynamics, with applications ranging from hypersonic flight to semiconductor manufacturing. Code verification of such methods is…
To minimise systematic errors in Monte Carlo simulations of charged particles, long range electrostatic interactions have to be calculated accurately and efficiently. Standard approaches, such as Ewald summation or the naive application of…
In the current noisy intermediate scale quantum era of quantum computation, available hardware is severely limited by both qubit count and noise levels, precluding the application of many current hybrid quantum-classical algorithms to…
In this paper, we develop an ensemble-based time-stepping algorithm to efficiently find numerical solutions to a group of linear, second-order parabolic partial differential equations (PDEs). Particularly, the PDE models in the group could…
We give an explicit stochastic Hamiltonian model of discontinuous unitary evolution for quantum spontaneous jumps like in a system of atoms in quantum optics, or in a system of quantum particles that interacts singularly with "bubbles"…
We propose a stochastic branching particle-based method for solving nonlinear non-conservative advection-diffusion-reaction equations. The method splits the evolution into an advection-diffusion step, based on a linearized Kolmogorov…
We present a rigorous derivation for off-lattice implementations of the so-called "random-walk" algorithm recently introduced by Wang and Landau [PRL 86, 2050 (2001)]. Originally developed for discrete systems, the algorithm samples…
For a large class of inhomogeneous interacting particle systems (IPS) on a lattice we develop a rigorous method for mapping them onto homogeneous IPS. Our novel approach provides a direct way of obtaining the statistical properties of such…
A multiscale numerical method is proposed for the solution of semi-linear elliptic stochastic partial differential equations with localized uncertainties and non-linearities, the uncertainties being modeled by a set of random parameters. It…
In general, adding a stochastic perturbation to a differential equation possessing an invariant manifold destroys the invariance as far as the It\^o formalism is used. In this article, we propose an invariantization method for perturbations…
We consider the (numerically motivated) Nanbu stochastic particle system associated to the spatially homogeneous Boltzmann equation for true hard potentials. We establish a rate of propagation of chaos of the particle system to the unique…
We present a numerical framework for the simulation of collisional plasma dynamics, based on a coupling between Direct Simulation Monte Carlo (DSMC) and Particle-in-Cell (PIC) methods for the Vlasov-Maxwell-Landau system. The approach…
We propose a splitting Hamiltonian Monte Carlo (SHMC) algorithm, which can be computationally efficient when combined with the random mini-batch strategy. By splitting the potential energy into numerically nonstiff and stiff parts, one…
A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is…
We are interested in developing a numerical method for capturing stationary sheaths, that a plasma forms in contact with a metallic wall. This work is based on a bi-species (ion/electron) Vlasov-Amp{\`e}re model proposed in [3]. The main…
This paper deals with the problem of simulating dense dispersed systems composed by large numbers of particles undergoing ballistic aggregation. The most classical approaches for dealing with such problems are represented by the so-called…
We present an adaptive spectral method for solving the Landau/Fokker-Planck equation for electron-ion systems. The heart of the algorithm is an expansion in Laguerre polynomials, which has several advantages, including automatic…
An iterative algorithm is presented for solving the RPA equations of linear response. The method optimally computes the energy-weighted moments of the strength function, allowing one to match the computational effort to the intrinsic…