Characterization of anharmonicities on complex potential energy surfaes: 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 {v_i} near a stationary point of the PES. The low-temperature behavior is described by a simple thermal expansion v~(i)(#beta#)=v_0~(i)[1-#alpha#_1~(i)/#beta# +#alpha#_2~(i)/2#beta#~2+…], where the coefficients {a_j~(i)} near obtained from perturbation theory. Using a simple diagrammatic representation, we give the complete expressions for the first two coefficients #alpah#_1 and #alpha#_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 we4ll suited to this problem. As an example, we have studied a model biopolymer that exhibits significant first and second order anharmonicities. To show the importance of treating anharmonicities properly, we have calculated the caloric curve (heat capacity) of the quantum Ne_(13) cluster in both the classical and quantum regimes. For this calculation we have used a superposition approximation and exact anharmonic classical corrections to second order in perturbation theory. When every vibrational mode of each inherent structure is treated separately, we find god agreement between out results and previous quantum Monte Carlo calculations.