Related papers: CATCH 1.4 User's Guide
Recent developments in parallel Markov chain Monte Carlo (MCMC) algorithms allow us to run thousands of chains almost as quickly as a single chain, using hardware accelerators such as GPUs. While each chain still needs to forget its initial…
The shell-model Monte Carlo (SMMC) technique transforms the traditional nuclear shell-model problem into a path-integral over auxiliary fields. We describe below the method and its applications to four physics issues: calculations of sdpf-…
Crystallization is a key step in macromolecular structure determination by crystallography. While a robust theoretical treatment of the process is available, due to the complexity of the system, the experimental process is still largely one…
Jet quenching is one of the major discoveries of the heavy-ion program at RHIC. While there is a wealth of data from RHIC that will soon be supplemented with measurements at the LHC, on the theoretical side the situation is less clear. A…
Simulation of the so-called monatomic and diatomic chains by the computer program CHAINPLOT is described. The simulation provides useful teaching material for undergraduate condensed matter physics lecture courses, and has revealed that for…
Crystals with small thickness along the beam exhibit top performance for steering particle beams through planar channeling. For such crystals, the effect of nuclear dechanneling plays an important role because it affects their efficiency.…
Sequential Monte Carlo methods, also known as particle methods, are a popular set of techniques for approximating high-dimensional probability distributions and their normalizing constants. These methods have found numerous applications in…
The concept of measurement in classical scattering is interpreted as an overlap of a particle packet with some area in phase space that describes the detector. Considering that usually we record the passage of particles at some point in…
We present a new Monte Carlo algorithm for studying site or bond percolation on any lattice. The algorithm allows us to calculate quantities such as the cluster size distribution or spanning probability over the entire range of site or bond…
The Dirichlet process (DP) is a fundamental mathematical tool for Bayesian nonparametric modeling, and is widely used in tasks such as density estimation, natural language processing, and time series modeling. Although MCMC inference…
We introduce the Quantization Monte Carlo method to solve thermal radiative transport equations with possibly several collision regimes, ranging from few collisions to massive number of collisions per time unit. For each particle in a given…
Triple-GEM detectors are a well known technology in high energy physics. In order to have a complete understanding of their behavior, in parallel with on beam testing, a Monte Carlo code has to be developed to simulate their response to the…
Sequential Monte Carlo (SMC), or particle filtering, is widely used in nonlinear state-space systems, but its performance often suffers from poorly approximated proposal and state-transition distributions. This work introduces a…
A new method, based on the simulated annealing algorithm and aimed at the inverse problem in the analysis of intergalactic (interstellar) complex spectra of hydrogen and metal lines, is presented. We consider the process of line formation…
We propose a hierarchy of multi-level kinetic Monte Carlo methods for sampling high-dimensional, stochastic lattice particle dynamics with complex interactions. The method is based on the efficient coupling of different spatial resolution…
Polymer-assisted ion transport underpins both energy storage technologies and emerging neuromorphic computing devices. Efficient modeling of ion migration is essential for understanding the performance of batteries and memristors, but it…
Planar channeling of 855 MeV electrons and positrons in straight and bent tungsten (110) crystal is simulated by means of the \MBNExplorer software package. The results of simulations for a broad range of bending radii are analyzed in terms…
We discuss the improvement in the accuracy of a Monte Carlo integration that can be obtained by optimization of the `a-priori weights' of the various channels. These channels may be either the strata in a stratified-sampling approach, or…
Monte Carlo simulations have been performed for 855 MeV and 6.3 GeV electrons channeling in silicon single crystals at circular bent (111) planes. The aim was to identify critical experimental parameters which effect the volume-deflection…
Monte Carlo sampling is a powerful toolbox of algorithmic techniques widely used for a number of applications wherein some noisy quantity, or summary statistic thereof, is sought to be estimated. In this paper, we survey the literature for…