Characterization of anharmonicities on complex potential energy surfaces: Perturbation theory and simulation

We have systematically investigated the effect of anharmonicity on the equilibrium properties of systems with a complex potential energy surface. Anharmonicities are modeled by the temperature dependence of the harmonic frequencies {νi} near a stationary point of the PES. The low-temperature behavior is described by a simple thermal expansion ν(i)(β)=ν0(i)[1−α1(i)/β+α2(i)/2β2+⋯], where the coefficients {αj(i)} are obtained from perturbation theory. Using a simple diagrammatic representation, we give the complete expressions for the first two coefficients α1 and α2 in terms of derivatives of the potential. This approach is illustrated for the example of a bulk Lennard-Jones system of 32 particles, in both the solid and the liquid states. We also determine the anharmonic frequencies from reversible-scaling Monte Carlo simulations, which appear particularly well suited to this problem. As an example, we have studied a model biopolymer that exhibits significant first and second order anharmonicities. To show ...

[1]  M. A. Hoef,et al.  Free energy of the Lennard-Jones solid , 2000 .

[2]  Thomas A. Weber,et al.  Hidden structure in liquids , 1982 .

[3]  E. R. Cowley,et al.  Summation of free-energy diagrams of an anharmonic crystal and equation of state for a Lennard-Jones solid , 1998 .

[4]  Philippe Choquard,et al.  The anharmonic crystal , 1967 .

[5]  D Thirumalai,et al.  The nature of folded states of globular proteins , 1992, Biopolymers.

[6]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

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

[8]  Liu,et al.  Variational path-integral theory of thermal properties of solids. , 1991, Physical review. B, Condensed matter.

[9]  Tartaglia,et al.  Role of unstable directions in the equilibrium and aging dynamics of supercooled liquids , 2000, Physical review letters.

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

[11]  R. Mattuck,et al.  A guide to Feynman diagrams in the many-body problem , 1976 .

[12]  David J. Wales,et al.  Coexistence in small inert gas clusters , 1993 .

[13]  Lee,et al.  New Monte Carlo algorithm: Entropic sampling. , 1993, Physical review letters.

[14]  R. Berry,et al.  Realistic master equation modeling of relaxation on complete potential energy surfaces: Partition function models and equilibrium results , 1998 .

[15]  K. Hinsen,et al.  Harmonicity in slow protein dynamics , 2000 .

[16]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[17]  Jonathan P. K. Doye,et al.  Calculation of thermodynamic properties of small Lennard‐Jones clusters incorporating anharmonicity , 1995 .

[18]  J. D. Doll,et al.  A heat capacity estimator for Fourier path integral simulations , 2000 .

[19]  Andersen,et al.  Testing mode-coupling theory for a supercooled binary Lennard-Jones mixture I: The van Hove correlation function. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[20]  H. Stanley,et al.  Thermodynamic and structural aspects of the potential energy surface of simulated water. , 2000, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  J. Straub,et al.  Approximate solution of the classical Liouville equation using Gaussian phase packet dynamics: Application to enhanced equilibrium averaging and global optimization , 1993 .

[22]  C. Brooks,et al.  Thermodynamics of protein folding: A statistical mechanical study of a small all-β protein , 1997 .

[23]  Jonathan P. K. Doye,et al.  Quantum partition functions from classical distributions: Application to rare-gas clusters , 2001 .

[24]  D. Wales,et al.  Energy landscapes of some model glass formers , 2001 .

[25]  F. Sciortino,et al.  Aging as dynamics in configuration space , 1999, cond-mat/9905090.

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

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

[28]  Liquid limits: glass transition and liquid-gas spinodal boundaries of metastable liquids. , 2000, Physical review letters.

[29]  Srikanth Sastry,et al.  The relationship between fragility, configurational entropy and the potential energy landscape of glass-forming liquids , 2000, Nature.

[30]  G. Voth,et al.  Modeling physical systems by effective harmonic oscillators: The optimized quadratic approximation , 1995 .

[31]  Potential energy landscape of a model glass former: thermodynamics, anharmonicities, and finite size effects. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[32]  Srikanth Sastry,et al.  Evaluation of the configurational entropy of a model liquid from computer simulations , 2000, cond-mat/0005225.