NOTE: This tutorial was intended to be used with Gromacs version 3.1.4 and may not be compatible with newer versions. The tutorial will be revised in the near future. Note, a brief tutorial on Free Energy calculation with Gromacs 3.3.x can be found at the DillGroup Wiki.

Free energy calculation methods

Experimentally, a difference in free energy is determined either from the relative probability of finding the system in a given state, or from the reversible work required to transform the system from one initial state into another. Computationally, both approaches can be used. However, counting how often the system is in one state, for instance how often a complex is formed when simulating a mixture of two compounds, is extremely inefficient. In practice, a perturbation is applied to force the transition from one state to another. Then, statistical procedures are used to calculate the work done on the system by the perturbation.

The free energy is a state function, which means that the free energy difference is only depending on the initial and end state, no matter what path is taken to go from one to the other. As a consequence, you can choose any nonphysical path you can think of to perform your calculations as long as you can relate them through thermodynamic cycles to the physical process you are interested in.

The most widely used methods to calculate free energy differences are:

Free energy calculations can be performed using either Monte-Carlo or Molecular Dynamics simulations, but during this practical only Molecular Dynamics will be used. As always with this type of methods, the extent of the sampling reached during the simulation will be one of the two primarily limiting factors concerning the accuracy of the calculations. The other one being the underlying model, or force field, used to describe the system.

The following exercises show some of the practical applications that can be made of the free energy calculations.

  1. solvation free energy:

  2. direct determination of the hydration free energy of toluene.
  3. relative solvation free energy:

  4. determination of the relative hydration free energy of p-cresol with respect to toluene.


<<- Back


Proceed ->>
This page was last modified on May 12th 2003