Related papers: A boundary-based net-exchange Monte Carlo method f…
The behavior of a Lattice Monte Carlo algorithm (if it is designed correctly) must approach that of the continuum system that it is designed to simulate as the time step and the mesh step tend to zero. However, we show for an algorithm for…
In image processing, solving inverse problems is the task of finding plausible reconstructions of an image that was corrupted by some (usually known) degradation operator. Commonly, this process is done using a generative image model that…
We describe a numerical model for the interaction of light with large raindrops using realistic nonspherical drop shapes. We apply geometrical optics and a Monte Carlo technique to perform ray traces through the drops. We solve the problem…
A new algorithm for Monte Carlo calculation of the double exchange model is studied. The algorithm is commonly applicable to wide classes of strongly correlated electron systems which involve itinerant electrons coupled with…
We describe a Monte Carlo procedure which allows sampling of the disjoint configuration spaces associated with crystalline and fluid phases, within a single simulation. The method utilises biased sampling techniques to enhance the…
Successful computer studies of glass-forming materials need to overcome both the natural tendency to structural ordering and the dramatic increase of relaxation times at low temperatures. We present a comprehensive analysis of eleven…
We introduce an efficient numerical implementation of a Markov Chain Monte Carlo method to sample a probability distribution on a manifold (introduced theoretically in Zappa, Holmes-Cerfon, Goodman (2018)), where the manifold is defined by…
We derive an efficient method for the insertion of structured particles in grand canonical Monte Carlo simulations of adsorption in very confining geometries. We extend this method to path integral simulations and use it to calculate the…
We present Arepo-MCRT, a novel Monte Carlo radiative transfer (MCRT) radiation-hydrodynamics (RHD) solver for the unstructured moving-mesh code Arepo. Our method is designed for general multiple scattering problems in both optically thin…
A new Monte-Carlo algorithm for calculating time-dependent radiative-transfer under the assumption of LTE is presented. Unlike flux-limited diffusion the method is polychromatic, includes scattering, and is able to treat the optically thick…
The standard Kernel Quadrature method for numerical integration with random point sets (also called Bayesian Monte Carlo) is known to converge in root mean square error at a rate determined by the ratio $s/d$, where $s$ and $d$ encode the…
One of the most demanding calculations is to generate random samples from a specified probability distribution (usually with an unknown normalizing prefactor) in a high-dimensional configuration space. One often has to resort to using a…
We demonstrate that the effects of matter upon neutrino propagation may be recast as the scattering of the initial neutrino wavefunction. Exchanging the differential, Schrodinger equation for an integral equation for the scattering matrix S…
In this work we present a theoretical study on the propagation of light in heterogeneous systems with fluctuating optical properties. To understand the consequences of the fluctuations we perform numerical calculations with uniform and non…
We propose a modification, based on the RESTART (repetitive simulation trials after reaching thresholds) and DPR (dynamics probability redistribution) rare event simulation algorithms, of the standard diffusion Monte Carlo (DMC) algorithm.…
A new code for evaluation of light absorption and scattering by interstellar dust grains is presented. The radiative transfer problem is solved using ray-tracing algorithm in a self-consistent and highly efficient way. The code demonstrates…
Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also…
Monte Carlo techniques based on indivisible energy packets are described for computing light curves and spectra for 3-D supernovae. The radiative transfer is time-dependent and includes all effects of O(v/c). Monte Carlo quantization is…
Distributed learning methods have gained substantial momentum in recent years, with communication overhead often emerging as a critical bottleneck. Gradient compression techniques alleviate communication costs but involve an inherent…
Inverse medium scattering is an ill-posed, nonlinear wave-based imaging problem arising in medical imaging, remote sensing, and non-destructive testing. Machine learning (ML) methods offer increased inference speed and flexibility in…