Thermodynamic integration methods perturbation

Considerable use has been made of the thermodynamic perturbation and thermodynamic integration methods in biochemical modelling, calculating the relative Gibbs energies of binding of inhibitors of biological macromolecules (e.g. proteins) with the aid of suitable thermodynamic cycles. Some applications to materials are described by Alfe et al. [11]. 

The evaluation of the free energy is essential to quantitatively treat a chemical process in condensed phase. In this section, we review methods of free-energy calculation within the context of classical statistical mechanics. We start with the standard free-energy perturbation and thermodynamic integration methods. We then introduce the method of distribution functions in solution. The method of energy representation is described in its classical form in this section, and is combined with the QM/MM methodology in the next section. 

Free-Energy Perturbation and Thermodynamic Integration Methods... 

The standard and often used methods to circumvent the difficulty associated with the form of Eq. (17-30) are the free-energy perturbation and thermodynamic integration methods [13,42,43], These methods are generally applicable to free-energy... 

Simulations, Time-dependent Methods and Solvation Models 16.1 Simulation Methods 16.1.1 Free Energy Methods 16.1.2 Thermodynamic Perturbation Methods 16.1.3 Thermodynamic Integration Methods 16.2 Time-dependent Methods ill 373 380 380 381 383 ... 

The relative free energy AA can now be computed by either thermodynamic integration or perturbation methods. 

What has been developed within the last 20 years is the computation of thermodynamic properties including free energy and entropy [12, 13, 14]. But the ground work for free energy perturbation was done by Valleau and Torrie in 1977 [15], for particle insertion by Widom in 1963 and 1982 [16, 17] and for umbrella sampling by Torrie and Valleau in 1974 and 1977 [18, 19]. These methods were primarily developed for use with Monte Carlo simulations continuous thermodynamic integration in MD was first described in 1986 [20]. 

As noted above, it is very difficult to calculate entropic quantities with any reasonable accmacy within a finite simulation time. It is, however, possible to calculate differences in such quantities. Of special importance is the Gibbs free energy, as it is the natoal thermodynamical quantity under normal experimental conditions (constant temperature and pressme. Table 16.1), but we will illustrate the principle with the Helmholtz free energy instead. As indicated in eq. (16.1) the fundamental problem is the same. There are two commonly used methods for calculating differences in free energy Thermodynamic Perturbation and Thermodynamic Integration. 

The next three chapters deal with the most widely used classes of methods free energy perturbation (FEP) [3], methods based on probability distributions and histograms, and thermodynamic integration (TI) [1, 2], These chapters represent a mix of traditional material that has already been well covered, as well as the description of new techniques that have been developed only recendy. The common thread followed here is that different methods share the same underlying principles. Chapter 5 is dedicated to a relatively new class of methods, based on calculating free energies from nonequilibrium dynamics. In Chap. 6, we discuss an important topic that has not received, so far, sufficient attention - the analysis of errors in free energy calculations, especially those based on perturbative and nonequilibrium approaches. 

---