Simulation of activation free energies in molecular systems

Abstract

A method is presented for determining activation free energies in complex molecular systems. The method relies on knowledge of the minimum energy path and bases the activation free energy calculation on moving along this path from a minimum to a saddle point. Use is made of a local reaction coordinate which describes the advance of the reaction in each segment of the minimum energy path. The activation free energy is formulated as a sum of two terms. The first is due to the change in the local reaction coordinate between the endpoints of each segment of the path. The second is due to the change in direction of the minimum energy path between consecutive segments. Both contributions can be obtained by molecular dynamics simulations with a constraint on the local reaction coordinate. The method is illustrated by applying it to a model potential and to the C7eq to C7ax transition in the alanine dipeptide. It is found that the term due to the change of direction in the reaction path can make a substantial contribution to the activation free energy.