Related papers: On the length of one-dimensional reactive paths
The replica exchange method is a powerful tool for overcoming slow relaxation in molecular simulations, but its efficiency depends strongly on the choice of the number and interval of replicas and their exchange probabilities. Here, we…
We present a path integral Monte Carlo study of the global superfluid fraction and local superfluid density in cylindrically-symmetric reservoirs of liquid $^{4}$He separated by nanoaperture arrays. The superfluid response to both…
We consider the first exit point distribution from a bounded domain $\Omega$ of the stochastic process $(X_t)_{t\ge 0}$ solution to the overdamped Langevin dynamics $$d X_t = -\nabla f(X_t) d t + \sqrt{h} \ d B_t$$ starting from the…
We address the possibility of performing numerical Monte Carlo simulations for the thermodynamics of quantum dissipative systems. Dissipation is considered within the Caldeira-Leggett formulation, which describes the system in the…
Markov chain Monte Carlo methods such as Gibbs sampling and simple forms of the Metropolis algorithm typically move about the distribution being sampled via a random walk. For the complex, high-dimensional distributions commonly encountered…
The author's presentation of multilevel Monte Carlo path simulation at the MCQMC 2006 conference stimulated a lot of research into multilevel Monte Carlo methods. This paper reviews the progress since then, emphasising the simplicity,…
The direct integration of the harmonic oscillator path integral obscures the fundamental structure of its discrete, imaginary time propagator (density matrix). This work, by first proving an operator identity for contracting two free…
Monte Carlo simulations are widely used to simulate complex molecular systems, but standard approaches suffer from metastability. Lately, the use of non-local proposal updates in a collective-variable (CV) space has been proposed in several…
When a Brownian particle, initially being in the potential well, overcomes the barrier and moves to the absorptive border, it still has a chance to be scattered back to the well by thermal fluctuations. We study this phenomenon carefully…
We analyse a multilevel Monte Carlo method for the approximation of distribution functions of univariate random variables. Since, by assumption, the target distribution is not known explicitly, approximations have to be used. We provide an…
Transition of a system between two states is an important but difficult problem in natural science. In this article we study the transition problem in the framework of transition path ensemble. Using the overdamped Langevin method, we…
A simple Monte Carlo procedure is described for simulating the multiple scattering and absorption of electrons with the incident energy in the range 1-50 keV moving through a slab of uniformly distributed material of given atomic number,…
We have developed a path integral Monte Carlo method for simulating helium films and apply it to the second layer of helium adsorbed on graphite. We use helium-helium and helium-graphite interactions that are found from potentials which…
Algorithms based on discretizing Langevin diffusion are popular tools for sampling from high-dimensional distributions. We develop novel connections between such Monte Carlo algorithms, the theory of Wasserstein gradient flow, and the…
A novel method for simulating the statistical mechanics of molecular systems in which both nuclear and electronic degrees of freedom are treated quantum mechanically is presented. The scheme combines a path integral description of the…
The interactions between a group of components are commonly studied in several areas of science (social science, biology, material science, complex dynamical systems, among others) using the methods of thermodynamics and statistical…
Motivated by the study of rare events for a typical genetic switching model in systems biology, in this paper we aim to establish the general two-scale large deviations for chemical reaction systems. We build a formal approach to explicitly…
We present a method to sample Markov-chain trajectories constrained to both the initial and final conditions, which we term Markov bridges. The trajectories are conditioned to end in a specific state at a given time. We derive the master…
We present a path integral formulation of Darcy's equation in one dimension with random permeability described by a correlated multi-variate lognormal distribution. This path integral is evaluated with the Markov chain Monte Carlo method to…
Path integral Monte Carlo approach is used to study the coupled quantum dynamics of the electron and nuclei in hydrogen molecule ion. The coupling effects are demonstrated by comparing differences in adiabatic Born--Oppenheimer and…