Related papers: Boundary processing for Monte Carlo Simulations in…
We present a new implementation of the Monte Carlo method to simulate the evolution of star clusters. The major improvement with respect to the previously developed codes is the treatment of the external tidal field taking into account for…
We review a recently devised Monte Carlo simulation method for the direct study of quasi-stationary properties of stochastic processes with an absorbing state. The method is used to determine the static correlation function and the…
Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However,…
The importance-sampling Monte Carlo algorithm appears to be the universally optimal solution to the problem of sampling the state space of statistical mechanical systems according to the relative importance of configurations for the…
A numerical method using implicit surface representations is proposed to solve the linearized Poisson-Boltzmann equations that arise in mathematical models for the electrostatics of molecules in solvent. The proposed method used an implicit…
A model for two-dimensional colloids confined laterally by "structured boundaries" (i.e., ones that impose a periodicity along the slit) is studied by Monte Carlo simulations. When the distance D between the confining walls is reduced at…
Kinetic equations model distributions of particles in position-velocity phase space. Often, one is interested in studying the long-time behavior of particles in high-collisional regimes in which an approximate (advection)-diffusion model…
Tempering is a popular tool in Bayesian computation, being used to transform a posterior distribution $p_1$ into a reference distribution $p_0$ that is more easily approximated. Several algorithms exist that start by approximating $p_0$ and…
Diffusion models enable the synthesis of highly accurate samples from complex distributions and have become foundational in generative modeling. Recently, they have demonstrated significant potential for solving Bayesian inverse problems by…
We propose numerical schemes that enable the application of particle methods for advection problems in general bounded domains. These schemes combine particle fields with Cartesian tensor product splines and a fictitious domain approach.…
Fractal Interpolation has been proposed in the literature as an efficient way to construct closure models for the numerical solution of coarse-grained Navier-Stokes equations. It is based on synthetically generating a scale-invariant…
While boundary plasmas in present-day tokamaks generally fall in a fluid regime, neutral species near the boundary often require kinetic models due to long mean-free-paths compared to characteristic spatial scales in the region. Monte-Carlo…
We address the computation of ground-state properties of chemical systems and realistic materials within the auxiliary-field quantum Monte Carlo method. The phase constraint to control the fermion phase problem requires the random walks in…
Generative diffusions are a powerful class of Monte Carlo samplers that leverage bridging Markov processes to approximate complex, high-dimensional distributions, such as those found in image processing and language models. Despite their…
Radial basis functions are a common mathematical tool used to construct a smooth interpolating function from a set of data points. A spatial prior based on thin-plate spline radial basis functions can be easily implemented resulting in a…
A propagator-based approach is investigated for Monte-Carlo (MC) modeling of neutral particles transport in fusion boundary plasmas. The propagator is essentially a Green function for the neutral kinetic equation, which depends on the…
We perform numerical studies including Monte Carlo simulations of high rotational symmetry random tilings. For computational convenience, our tilings obey fixed boundary conditions in regular polygons. Such tilings are put in correspondence…
We demonstrate the feasibility of performing sufficient configurational sampling of disordered oxides directly from first principles without resorting to the use of fitted models such as cluster expansion. This is achieved by harnessing the…
In this paper we design a neural interpolation operator to improve the boundary data for regional weather models, which is a challenging problem as we are required to map multi-scale dynamics between grid resolutions. In particular, we…
We present a method of deriving two boundary conditions at a thin membrane for diffusion from experimental data. This method can be really useful in complex membrane systems in which we do not know mechanisms of processes occurring within…