Calculation of free energy through successive umbrella sampling

被引:238
作者
Virnau, P [1 ]
Müller, M
机构
[1] MIT, Dept Phys, Cambridge, MA 02139 USA
[2] Univ Mainz, Inst Phys, D-55099 Mainz, Germany
关键词
D O I
10.1063/1.1739216
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
We consider an implementation of umbrella sampling in which the pertinent range of states is subdivided into small windows that are sampled consecutively and linked together. This allows us to simulate without a weight function or to extrapolate the results to the neighboring window in order to estimate a weight function. Additionally, we present a detailed error analysis in which we demonstrate that the error in umbrella sampling is controlled and, in the absence of sampling difficulties, independent of the window sizes. In this case, the efficiency of our implementation is comparable to a multicanonical simulation with a very good weight function, which in our scheme does not need to be known ahead of time. The analysis also allows us to detect sampling difficulties such as correlations between adjacent windows and provides a test of equilibration. We exemplify the scheme by simulating the liquid-vapor coexistence in a Lennard-Jones system. (C) 2004 American Institute of Physics.
引用
收藏
页码:10925 / 10930
页数:6
相关论文
共 19 条
[1]   Multicanonical recursions [J].
Berg, BA .
JOURNAL OF STATISTICAL PHYSICS, 1996, 82 (1-2) :323-342
[2]   MULTICANONICAL ENSEMBLE - A NEW APPROACH TO SIMULATE 1ST-ORDER PHASE-TRANSITIONS [J].
BERG, BA ;
NEUHAUS, T .
PHYSICAL REVIEW LETTERS, 1992, 68 (01) :9-12
[3]   FINITE-SIZE EFFECTS AT ASYMMETRIC 1ST-ORDER PHASE-TRANSITIONS [J].
BORGS, C ;
KOTECKY, R .
PHYSICAL REVIEW LETTERS, 1992, 68 (11) :1734-1737
[4]   OPTIMIZED MONTE-CARLO DATA-ANALYSIS [J].
FERRENBERG, AM ;
SWENDSEN, RH .
PHYSICAL REVIEW LETTERS, 1989, 63 (12) :1195-1198
[5]   NEW MONTE-CARLO TECHNIQUE FOR STUDYING PHASE-TRANSITIONS [J].
FERRENBERG, AM ;
SWENDSEN, RH .
PHYSICAL REVIEW LETTERS, 1988, 61 (23) :2635-2638
[6]   Global optimization by energy landscape paving [J].
Hansmann, UHE ;
Wille, LT .
PHYSICAL REVIEW LETTERS, 2002, 88 (06) :68105/1-68105/4
[7]   HISTOGRAM MONTE-CARLO RENORMALIZATION-GROUP METHOD FOR PHASE-TRANSITION MODELS WITHOUT CRITICAL SLOWING DOWN [J].
HU, CK .
PHYSICAL REVIEW LETTERS, 1992, 69 (19) :2739-2742
[8]   The evaporation/condensation transition of liquid droplets [J].
MacDowell, LG ;
Virnau, P ;
Müller, M ;
Binder, K .
JOURNAL OF CHEMICAL PHYSICS, 2004, 120 (11) :5293-5308
[9]   SIMULATED TEMPERING - A NEW MONTE-CARLO SCHEME [J].
MARINARI, E ;
PARISI, G .
EUROPHYSICS LETTERS, 1992, 19 (06) :451-458
[10]   2D crystal shapes, droplet condensation, and exponential slowing down in simulations of first-order phase transitions [J].
Neuhaus, T ;
Hager, JS .
JOURNAL OF STATISTICAL PHYSICS, 2003, 113 (1-2) :47-83