A way to include anharmonicity fully in free energy predictions is provided by thermodynamic integration (TI). These set of techniques are based on performing a series of molecular dynamics (MD) simulations where an external parameter (real or artificial) is varied and free energy differences are obtained by integration of a related quantity.