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 Berry-Esseen theorem for Feynman-Kac and interacting particle models

On the concentration properties of Interacting particle processes

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

Berry–Esseen Bounds for Kernel Estimates of Stationaryprocesses

Weighted Sums and Berry-Esseen type estimates in Free Probability Theory

Eigenfunction martingale estimators for interacting particle systems and their mean field limit

Improved nonuniform Berry--Esseen-type bounds

A multivariate Berry--Esseen theorem for time-dependent expanding dynamical systems

Long run convergence of discrete-time interacting particle systems of the McKean-Vlasov type

On the Particle Approximation of Lagged Feynman-Kac Formulae

Non-uniform Berry-Esseen theorem and Edgeworth expansions with applications to transport distances for weakly dependent random variables

Stein's method and a cubic mean-field model

Exponential Convergence for the Fredrickson–Andersen One-Spin Facilitated Model

Berry-Esseen estimates for regenerative processes under weak moment assumptions

A Berry-Esseen bound of order $ 1/\sqrt{n} $ for martingales

On a particle approximation to the Dean-Kawasaki type equation with logarithmic interactions

Kinetically constrained models out of equilibrium

Non-uniform Berry--Esseen bounds for Gaussian, Poisson and Rademacher processes

Nonstationary generalized TASEP in KPZ and jamming regimes

Non-extinction of a Fleming-Viot particle model

An exact representation of the fermion dynamics in terms of Poisson processes and its connection with Monte Carlo algorithms