Abstract:In this paper we examine the numerical approximation of the limiting invariant measure associated with Feynman-Kac formulae. These are expressed in a discrete time formulation and are associated with a Markov chain and a potential function. The typical application considered here is the computation of eigenvalues associated with non-negative operators as found, for example, in physics or particle simulation of rare-events. We focus on a novel \emph{lagged} approximation of this invariant measure, based upon the introduction of a ratio of time-averaged Feynman-Kac marginals associated with a positive operator iterated $l \in\mathbb{N}$ times; a lagged Feynman-Kac formula. This estimator and its approximation using Diffusion Monte Carlo (DMC) have been extensively employed in the physics literature. In short, DMC is an iterative algorithm involving $N\in\mathbb{N}$ particles or walkers simulated in parallel, that undergo sampling and resampling operations. In this work, it is shown that for the DMC approximation of the lagged Feynman-Kac formula, one has an almost sure characterization of the $\mathbb{L}_1$-error as the time parameter (iteration) goes to infinity and this is at most of $\mathcal{O}(\exp\{-\kappa l\}/N)$, for $\kappa>0$. In addition a non-asymptotic in time, and time uniform $\mathbb{L}_1-$bound is proved which is $\mathcal{O}(l/\sqrt{N})$. We also prove a novel central limit theorem to give a characterization of the exact asymptotic in time variance. This analysis demonstrates that the strategy used in physics, namely, to run DMC with $N$ and $l$ small and, for long time enough, is mathematically justified. Our results also suggest how one should choose $N$ and $l$ in practice. We emphasize that these results are not restricted to physical applications; they have broad relevance to the general problem of particle simulation of the Feynman-Kac formula.

A duality formula and a particle Gibbs sampler for continuous time Feynman-Kac measures on path spaces

Coag-Frag duality for a class of stable Poisson-Kingman mixtures

On the Particle Approximation of Lagged Feynman-Kac Formulae

Feynman-Kac equation revisited

Compound Poisson distributions for random dynamical systems using probabilistic approximations

The particle approximation of quasi-stationary distributions. Part~I: concentration bounds in the uniform case

Feynman-Kac formula for stochastic hybrid systems

A Berry-Esseen theorem for Feynman-Kac and interacting particle models

On the path integral simulation of space-time fractional Schroedinger equation with time independent potentials

Long-time behavior for discretization schemes of Fokker-Planck equations via couplings

Subordination principle and Feynman-Kac formulae for generalized time-fractional evolution equations

Brownian motion on disconnected sets, basic hypergeometric functions, and some continued fractions of Ramanujan

Birkhoff sum convergence of Fréchet observables to stable laws for Gibbs-Markov systems and applications

Gibbs measures on Brownian currents

On particle Gibbs sampling

Quantitative approximate independence for continuous mean field Gibbs measures

Ergodicity of Langevin Dynamics and its Discretizations for Non-smooth Potentials

The Bouncy Particle Sampler: A Nonreversible Rejection-Free Markov Chain Monte Carlo Method

Compound Poisson particle approximation for McKean-Vlasov SDEs

Parametrized Families of Gibbs Measures and their Statistical Inference

Propagation of chaos and residual dependence in Gibbs measures on finite sets