Free energy reconstruction from steered dynamics without post-processing

被引:21
作者
Athenes, Manuel [1 ,2 ]
Marinica, Mihai-Cosmin [1 ]
机构
[1] CEA Saclay, Serv Rech Met Phys, Dept Mat Nucl, F-91191 Gif Sur Yvette, France
[2] LLNL, Condensed Matter & Mat Div, Phys & Life Sci Directorate, Livermore, CA 94551 USA
关键词
Free-energy calculations; Statistical thermodynamics; Computer chemistry; Molecular simulation; Monte Carlo method; MOLECULAR-DYNAMICS; ENSEMBLE AVERAGES; SIMULATIONS; INFORMATION; TRANSITION; LANDSCAPE; ALGORITHM; DIFFUSION; EFFICIENT; PATHWAYS;
D O I
10.1016/j.jcp.2010.06.003
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Various methods achieving importance sampling in ensembles of nonequilibrium trajectories enable one to estimate free energy differences and, by maximum-likelihood post-processing, to reconstruct free energy landscapes. Here, based on Bayes theorem, we propose a more direct method in which a posterior likelihood function is used both to construct the steered dynamics and to infer the contribution to equilibrium of all the sampled states. The method is implemented with two steering schedules. First, using non-autonomous steering, we calculate the migration barrier of the vacancy in Fe-a. Second, using an autonomous scheduling related to metadynamics and equivalent to temperature-accelerated molecular dynamics, we accurately reconstruct the two-dimensional free energy landscape of the 38-atom Lennard-Jones cluster as a function of an orientational bond-order parameter and energy, down to the solid-solid structural transition temperature of the cluster and without maximum-likelihood post-processing. (C) 2010 Elsevier Inc. All rights reserved.
引用
收藏
页码:7129 / 7146
页数:18
相关论文
共 70 条
[1]   Large-scale conformational sampling of proteins using temperature-accelerated molecular dynamics [J].
Abrams, Cameron F. ;
Vanden-Eijnden, Eric .
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, 2010, 107 (11) :4961-4966
[2]   Development of an interatomic potential for phosphorus impurities in α-iron [J].
Ackland, GJ ;
Mendelev, MI ;
Srolovitz, DJ ;
Han, S ;
Barashev, AV .
JOURNAL OF PHYSICS-CONDENSED MATTER, 2004, 16 (27) :S2629-S2642
[3]   Free energy landscape from path-sampling:: application to the structural transition in LJ38 [J].
Adjanor, G. ;
Athenes, M. ;
Calvo, F. .
EUROPEAN PHYSICAL JOURNAL B, 2006, 53 (01) :47-60
[4]   Gibbs free-energy estimates from direct path-sampling computations -: art. no. 234104 [J].
Adjanor, G ;
Athènes, M .
JOURNAL OF CHEMICAL PHYSICS, 2005, 123 (23)
[5]  
Adjanor G, 2008, AIP CONF PROC, V999, P186, DOI 10.1063/1.2918104
[6]   Web ensemble averages for retrieving relevant information from rejected Monte Carlo moves [J].
Athenes, M. .
EUROPEAN PHYSICAL JOURNAL B, 2007, 58 (01) :83-95
[7]   Parallel Monte Carlo simulations using a residence weight algorithm -: art. no. 016701 [J].
Athènes, M .
PHYSICAL REVIEW E, 2002, 66 (01) :1-016701
[8]   A path-sampling scheme for computing thermodynamic properties of a many-body system in a generalized ensemble [J].
Athènes, M .
EUROPEAN PHYSICAL JOURNAL B, 2004, 38 (04) :651-663
[9]   Identification of novel diffusion cycles in B2 ordered phases by Monte Carlo simulation [J].
Athenes, M ;
Bellon, P ;
Martin, G .
PHILOSOPHICAL MAGAZINE A-PHYSICS OF CONDENSED MATTER STRUCTURE DEFECTS AND MECHANICAL PROPERTIES, 1997, 76 (03) :565-585
[10]   Multiple-Replica Exchange with Information Retrieval [J].
Athenes, Manuel ;
Calvo, Florent .
CHEMPHYSCHEM, 2008, 9 (16) :2332-2339