Related papers: Improvements in charged lepton and photon propagat…
Monte Carlo simulations play a crucial role in all stages of particle collider experiments. There has been a long-term trend in HEP of both increasing collision energies and the luminosity. As a result, the requirements for MC simulations…
Simulating samples from arbitrary probability distributions is a major research program of statistical computing. Recent work has shown promise in an old idea, that sampling from a discrete distribution can be accomplished by perturbing and…
Classical molecular dynamics simulation is performed mostly using the established velocity Verlet integrator or other symplectic propagation schemes. In this work, an alternative formulation of numerical propagators for classical molecular…
We propose an algorithm for molecular dynamics or Monte Carlo simulations that uses an interpolation procedure to estimate potential energy values from energies and gradients evaluated previously at points of a simplicial mesh. We chose an…
Electron microscopy is a powerful tool for studying the properties of materials down to their atomic structure. In many cases, the quantitative interpretation of images requires simulations based on atomistic structure models. These…
For the particle detectors based on noble gases or liquids, it is essential to understand the transport dynamic and the properties of the electrons. We report the development of a tool for electron transport in noble gases He, Ne, Ar, Kr,…
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 study some aspects of a Monte Carlo method invented by Maggs and Rossetto for simulating systems of charged particles. It has the feature that the discretized electric field is updated locally when charges move. Results of simulations of…
We present an experimental demonstration of boson sampling as a hardware accelerator for Monte Carlo integration. Our approach leverages importance sampling to factorize an integrand into a distribution that can be sampled using quantum…
AI accelerators have proliferated in data centers in recent years and are now almost ubiquitous. In addition, their computational power and, most importantly, their energy efficiency are up to orders of magnitude higher than that of…
A Monte Carlo model of radiative transfer in multi-phase dusty media is applied to the situation of stars and clumpy dust in a sphere or a disk. The distribution of escaping and absorbed photons are shown for various filling factors and…
Motivated by applications to Bayesian inference for statistical models with orthogonal matrix parameters, we present $\textit{polar expansion},$ a general approach to Monte Carlo simulation from probability distributions on the Stiefel…
Policy-guided Monte Carlo is an adaptive method to simulate classical interacting systems. It adjusts the proposal distribution of the Metropolis-Hastings algorithm to maximize the sampling efficiency, using a formalism inspired by…
Many sophisticated computer models have been developed to understand the behaviour of particle accelerators. Even these complex models often do not describe the measured data. Interactions of the beam with external fields, other particles…
A significant number of high power proton beams are available or will go online in the near future. This provides exciting opportunities for new fixed target experiments and the search for new physics in particular. In this note we will…
Cellular automata are a class of computational models based on simple rules and algorithms that can simulate a wide range of complex phenomena. However, when using conventional computers, these 'simple' rules are only encapsulated at the…
We develop a comprehensive approach to simulate the deformation of mirrors and lenses due to thermal and mechanical stresses that couples efficiently to photon-based optics simulations. This expands upon previous work where we demonstrated…
The Transfer Matrix Method (TMM) is a widely used technique for modeling linear propagation of electromagnetic waves through stratified layered media. However, since its extension to inhomogeneous and nonlinear systems is not…
A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed…
A simple reweighting scheme is proposed for Monte Carlo simulations of interacting particle systems, permitting one to study various parameter values in a single study, and improving efficiency by an order of magnitude. Unlike earlier…