4 Inria - Institut National de Recherche en Informatique et en Automatique (Domaine de Voluceau Rocquencourt - BP 105 78153 Le Chesnay Cedex - France)
4 Inria - Institut National de Recherche en Informatique et en Automatique (Domaine de Voluceau Rocquencourt - BP 105 78153 Le Chesnay Cedex - France)
5 CUHK - The Chinese University of Hong Kong [Hong Kong] (Shatin, NT, Hong Kong SAR - Chine) "> CUHK - The Chinese University of Hong Kong [Hong Kong]
Résumé
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 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 ∈ 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 ∈ 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 L1-error as the time parameter (iteration) goes to infinity and this is at most of O(exp{-κl}/N ), for κ > 0. In addition a non-asymptotic in time, and time uniform L1-bound is proved which is O(l/ √ 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, which is utilized in a great variety of scientific and engineering fields.
Origine | Fichiers produits par l'(les) auteur(s) |
---|
Michel Caffarel : Connectez-vous pour contacter le contributeur
https://hal.science/hal-04705748
Soumis le : lundi 23 septembre 2024-12:00:16
Dernière modification le : lundi 4 novembre 2024-14:48:47
Dates et versions
- HAL Id : hal-04705748 , version 1
- ARXIV : 2407.15494
- DOI : 10.48550/arXiv.2407.15494
Citer
Elsiddig Awadelkarim, Michel Caffarel, Pierre Del Moral, Ajay Jasra. On the Particle Approximation of Lagged Feynman-Kac Formulae. 2024. ⟨hal-04705748⟩
Collections
45
Consultations
9
Téléchargements