Sampling along reaction coordinates with the Wang-Landau method

The multiple range random walk algorithm recently proposed by Wang and Landau [2001, Phys. Rev. Lett., 86, 2050] is adapted to the computation of free energy profiles for molecular systems along reaction coordinates. More generally, we show how to extract partial averages in various statistical ensembles without invoking simulations with constraints, biasing potentials or unknown parameters. The method is illustrated on a model 10-dimensional potential energy surface, for which analytical results are obtained. It is then applied to the potential of mean force associated with the dihedral angle of the butane molecule in the gas phase and in carbon tetrachloride solvent. Finally, isomerization in a small rocksalt cluster, (NaF)4, is investigated in the microcanonical ensemble, and the results are compared to those of parallel tempering Monte Carlo.

[1]  Watanabe,et al.  Direct dynamical calculation of entropy and free energy by adiabatic switching. , 1990, Physical review letters.

[2]  F. Calvo,et al.  Monte-Carlo simulations of rotating clusters , 1998 .

[3]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[4]  Mihaly Mezei,et al.  Convergence acceleration in Monte Carlo computer simulation on water and aqueous solutions , 1983 .

[5]  Density-of-states Monte Carlo method for simulation of fluids , 2002, cond-mat/0201470.

[6]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[7]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[8]  Joost VandeVondele,et al.  Efficient multidimensional free energy calculations for ab initio molecular dynamics using classical bias potentials , 2000 .

[9]  G. Ciccotti,et al.  Constrained reaction coordinate dynamics for the simulation of rare events , 1989 .

[10]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[11]  Juan J. de Pablo,et al.  A biased Monte Carlo technique for calculation of the density of states of polymer films , 2002 .

[12]  Giovanni Ciccotti,et al.  Sampling of molecular conformations by molecular dynamics techniques , 1993 .

[13]  G. Ciccotti,et al.  Activation energies by molecular dynamics with constraints , 1991 .

[14]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[15]  Christoph Dellago,et al.  Transition path sampling: throwing ropes over mountains in the dark , 2000 .

[16]  Harold A. Scheraga,et al.  Preferential sampling near solutes in monte carlo calculations on dilute solutions , 1977 .

[17]  H. C. Andersen,et al.  Molecular dynamics study of melting and freezing of small Lennard-Jones clusters , 1987 .

[18]  K. J. Resch,et al.  Comment on “Manipulating the frequency-entangled states by an acoustic-optical modulator” , 2001 .

[19]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[20]  M. Tuckerman,et al.  An Adiabatic Molecular Dynamics Method for the Calculation of Free Energy Profiles , 2002 .

[21]  D. L. Freeman,et al.  Phase changes in 38-atom Lennard-Jones clusters. II. A parallel tempering study of equilibrium and dynamic properties in the molecular dynamics and microcanonical ensembles , 2000, physics/0003072.

[22]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[23]  Sidney Yip,et al.  Optimized free-energy evaluation using a single reversible-scaling simulation , 1999 .

[24]  Charles L. Brooks,et al.  The thermodynamics of solvophobic effects: A molecular‐dynamics study of n‐butane in carbon tetrachloride and water , 1990 .

[25]  Juan J. de Pablo,et al.  Monte Carlo simulation of proteins through a random walk in energy space , 2002 .

[26]  M. Tuckerman,et al.  On the use of the adiabatic molecular dynamics technique in the calculation of free energy profiles , 2002 .