378 articles – 175 references  [version française]
HAL: hal-00477302, version 1

Short view  Export this paper
Long-time convergence of an Adaptive Biasing Force method: the bi-channel case
Lelievre T., Minoukadeh K.
http://hal.archives-ouvertes.fr/hal-00477302
Available versions:
Preprint, Working Paper, ...
Mathematics/Analysis of PDEs
Long-time convergence of an Adaptive Biasing Force method: the bi-channel case
Tony Lelievre () 1, Kimiya Minoukadeh () 1
1:  Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS)
http://cermics.enpc.fr/
Ecole des Ponts ParisTech
6 et 8 avenue Blaise Pascal Cité Descartes - Champs sur Marne 77455 Marne la Vallée Cedex 2
France
We present convergence results for an adaptive algorithm to compute free energies, namely the adaptive biasing force (ABF) method. The free energy is the effective potential associated to a so-called reaction coordinate (RC). Computing free energy differences remains an important challenge in molecular dynamics due to the presence of meta-stable regions in the potential energy surface. The ABF method uses an on-the-fly estimate of the free energy to bias dynamics and overcome metastability. Using entropy arguments and logarithmic Sobolev inequalities, previous results have shown that the rate of convergence of the ABF method is limited by the metastable features of the canonical measures conditioned to being at fixed values of the RC. In this paper, we present an improvement on the existing results, in the presence of such metastabilities, which is a generic case encountered in practice. More precisely, we study the so-called bi-channel case, where two channels along the RC direction exist between an initial and final state, the channels being separated from each other by a region of very low probability. With hypotheses made on `channel-dependent' conditional measures, we show on a bi-channel model that we introduce, that the convergence of the ABF method is in fact not limited by metastabilities in directions orthogonal to the RC under two crucial assumptions: (i) exchange between the two channels is possible for some values of the RC and (ii) the free energy is a good bias in each channel.
English

free energy computation – nonlinear partial differential equation – Fokker-Planck equation – logarithmic Sobolev inequality

Project Id Méthodes Géométriques et échantillonnage: Application à la simulation moléculaire, MEGAS, ANR-09-BLAN-0216-01

Attached file list to this document: 
TEX
Bichannel.tex(93.4 KB)
bichannel-sketch2.eps(232.7 KB)
PDF
Bichannel.pdf(404.1 KB)
PS
Bichannel.ps(1.1 MB)