Related papers: On the length of one-dimensional reactive paths
We give an introduction to the calculation of path integrals on a lattice, with the quantum harmonic oscillator as an example. In addition to providing an explicit computational setup and corresponding pseudocode, we pay particular…
Traditional stochastic modeling of reactive systems limits the domain of applicability of the associated path thermodynamics to systems involving a single elementary reaction at the origin of each observed change in composition. An…
We present a numerical and partially analytical study of classical particles obeying a Langevin equation that describes diffusion on a surface modeled by a two dimensional potential. The potential may be either periodic or random. Depending…
Using the Monte Carlo simulation method for bosonic reaction-diffusion systems introduced recently [S.-C. Park, Phys. Rev. E {\bf 72}, 036111 (2005)], one dimensional bosonic models are studied and compared to the corresponding Langevin…
Langevin algorithms are popular Markov chain Monte Carlo methods that are often used to solve high-dimensional large-scale sampling problems in machine learning. The most classical Langevin Monte Carlo algorithm is based on the overdamped…
In plasmas, distribution functions often demonstrate long anisotropic tails or otherwise significant deviations from local Maxwellians. The tails, especially if they are pulled out from the bulk, pose a serious challenge for numerical…
In numerical studies of diffusive dynamics, two different action functionals are often used to specify the probability distribution of trajectories, one of which requiring the evaluation of the second derivative of the potential in addition…
Simulating percolation and critical phenomena of labelled species inside films composed of single-component linear homogeneous macromolecules using molecular Monte Carlo method in 3 dimensions, we study dependence of these conducting…
Recent work incorporating geometric ideas in Markov chain Monte Carlo is reviewed in order to highlight these advances and their possible application in a range of domains beyond Statistics. A full exposition of Markov chains and their use…
The random walk problem is studied in two and three dimensions in the presence of a random distribution of static traps. An efficient Monte Carlo method, based on a mapping onto a polymer model, is used to measure the survival probability…
Based on the path integral representation of the partition function of a many body system with separable two body interaction we propose a systematic extension of the perturbed static path approximation (PSPA) to lower temperatures.…
We demonstrate that the multicanonical approach is not restricted to Monte Carlo simulations, but can also be applied to simulation techniques such as molecular dynamics, Langevin, and hybrid Monte Carlo algorithms. The effectiveness of the…
A theoretical investigation of quantum-transport phenomena in mesoscopic systems is presented. In particular, a generalization to ``open systems'' of the well-known semiconductor Bloch equations is proposed. The presence of spatial boundary…
We present some considerations about the parallel implementations of the kinetic (Monte Carlo) version of the Ising model. In some cases the equilibrium distribution of the parallel version does not present the symmetry breaking phenomenon…
Consider a randomly-oriented two dimensional Manhattan lattice where each horizontal line and each vertical line is assigned, once and for all, a random direction by flipping independent and identically distributed coins. A deterministic…
We study a sequential Monte Carlo algorithm to sample from the Gibbs measure with a non-convex energy function at a low temperature. We use the practical and popular geometric annealing schedule, and use a Langevin diffusion at each…
The time evolution of the thermally activated decay rates is considered. This evolution is of particular importance for the recent nanoscale experiments discussed in the literature, where the potential barrier is relatively low (or the…
An efficient Path Integral Monte Carlo procedure is proposed to simulate the behavior of quantum many-body dissipative systems described within the framework of the influence functional. Thermodynamic observables are obtained by Monte Carlo…
We investigate Monte Carlo simulation strategies for determining the effective ("depletion") potential between a pair of hard spheres immersed in a dense sea of much smaller hard spheres. Two routes to the depletion potential are…
Tau-leaping is a popular discretization method for generating approximate paths of continuous time, discrete space, Markov chains, notably for biochemical reaction systems. To compute expected values in this context, an appropriate…