We prove the consistency of an adaptive importance sampling strategy based on biasing the potential energy function V of a diffusion process ; for the sake of simplicity, periodic boundary conditions are assumed, so that lives on the flat d-dimensional torus. The goal is to sample its invariant distribution . The bias , where is the new (random and time-dependent) potential function, acts only on some coordinates of the system, and is designed to flatten the corresponding empirical occupation measure of the diffusion X in the large-time regime. The diffusion process writes , where the bias is function of the key quantity : a probability occupation measure which depends on the past of the process, i.e. on . We are thus dealing with a self-interacting diffusion. In this note, we prove that when t goes to infinity, almost surely converges to μ. Moreover, the approach is justified by the convergence of the bias to a limit that has an interpretation in terms of a free energy. The main argument is a change of variables, which formally validates the consistency of the approach. The convergence is then rigorously proven adapting the ODE method from stochastic approximation.
Nous prouvons la consistance d'une méthode adaptative d'échantillonnage préférentiel, basée sur le biaisage du potentiel d'un processus de diffusion ; pour simplifier, des conditions au bord périodiques sont appliquées, si bien que est à valeurs dans le tore plat d-dimensionnel. L'objectif est d'échantillonner sa mesure invariante . Le biais , où est le nouveau potentiel (aléatoire et dépendant du temps), n'agit que sur certaines coordonnées du système, et est construit de façon à rendre la mesure d'occupation empirique correspondante uniforme en temps long. Le processus de diffusion s'écrit , où le biais est fonction de la quantité : une mesure d'occupation qui dépend du passé du processus, i.e. de . Nous étudions ainsi un processus de diffusion en auto-interaction. Dans cette note, nous prouvons la convergence presque sûre de vers μ quand t tend vers l'infini. De plus, l'approche est justifiée par la convergence du biais, vers une limite qui s'interprète en termes d'énergie libre. L'argument principal de la preuve est un changement de variables, qui formellement valide la consistance de l'approche. La convergence est alors prouvée rigoureusement en adaptant la méthode de l'EDO venant de l'approximation stochastique.
Accepted:
Published online:
Michel Benaïm 1; Charles-Édouard Bréhier 2
@article{CRMATH_2016__354_8_842_0, author = {Michel Bena{\"\i}m and Charles-\'Edouard Br\'ehier}, title = {Convergence of adaptive biasing potential methods for diffusions}, journal = {Comptes Rendus. Math\'ematique}, pages = {842--846}, publisher = {Elsevier}, volume = {354}, number = {8}, year = {2016}, doi = {10.1016/j.crma.2016.05.011}, language = {en}, }
Michel Benaïm; Charles-Édouard Bréhier. Convergence of adaptive biasing potential methods for diffusions. Comptes Rendus. Mathématique, Volume 354 (2016) no. 8, pp. 842-846. doi : 10.1016/j.crma.2016.05.011. https://comptes-rendus.academie-sciences.fr/mathematique/articles/10.1016/j.crma.2016.05.011/
[1] Well-tempered metadynamics: a smoothly converging and tunable free-energy method, Phys. Rev. Lett., Volume 100 (2008) no. 2
[2] M. Benaïm, C.-É. Bréhier, Convergence of adaptive biasing potential methods for diffusion processes, in preparation.
[3] Self-interacting diffusions, Probab. Theory Relat. Fields, Volume 122 (2002) no. 1, pp. 1-41
[4] Free energy calculations: an efficient adaptive biasing potential method, J. Phys. Chem. B, Volume 114 (2010), pp. 5823-5830
[5] Molecular Dynamics, Interdisciplinary Applied Mathematics, vol. 39, Springer, Cham, 2015 (With deterministic and stochastic numerical methods)
[6] Free Energy Computations, Imperial College Press, London, 2010 (A mathematical perspective)
[7] Self-healing umbrella sampling: a non-equilibrium approach for quantitative free energy calculations, J. Phys. Chem. B, Volume 110 (2006) no. 29, pp. 14011-14013
Cited by Sources:
Comments - Policy