Finite difference thermodynamic integration

A combination of thermodynamic perturbation, from which the slope of the free energy change with respect to X is obtained, and thermodynamic integration of the extrapolated values of these slopes has also been used. This approach has been called finite difference thermodynamic integration (FDTI). 

M. Mezei, J. Chem. Phys., 86, 7084 (1987). The Finite Difference Thermodynamic Integration, Tested on Calculating the Hydration Free Energy Difference between Acetone and Dimethylamine in Water. 

The finite difference thermodynamic integration combines the thermodynamic perturbation and integration methods. The value of the integrands at the integration points is determined by thermodynamic perturbation. In certain cases this approach has been shown to lead to faster convergence of the free energy difference calculation. 

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. 

Internal energy and a number of other thermodynamic variables (defined later) are state functions and are, therefore, properties of the system. Since state functions can be expressed mathematically as functions of thermodynamic coordinates such as temperature and pressure, their values can always be identified with points on a graph. The differential of a state function is spoken of as an infinitesimal change in the property. The integration of such a differential results in a finite difference between two values of the property. For example,... 

No matter how defined, restricting our A a-like solutions to be the lowest-energy solutions, FONs never occur. Instead, for the case of near-degeneracy, such solutions will involve orbitals that are linear combinations of symmetry-distinct orbitals with complex coefficients. With such extensions we can use integral occupation numbers and employ finite differences to evaluate derivatives in Parr s exact HK-based thermodynamic theory of the electron gas. The HK variational problem... 

Numerical methods include those based on finite difference calculus. They are ideally suited for tabulated experimental data such as one finds in thermodynamic tables. They also include methods of solving simultaneous linear equations, curve fitting, numerical solution of ordinary and partial differential equations and matrix operations. In this appendix, numerical interpolation, integration, and differentiation are considered. Information about the other topics is available in monographs by Hornbeck [2] and Lanczos [3]. 

The two-phase integrated finite difference flow code TOUGH2 was selected to perform the hydro-logical calculations for this I BEX analysis. TOUGH2 is a general-purpose program for the simulation of multi-dimensional, multiphase fluid flow in porous and fractured media. Fluid advec-tion is calculated by a multiphase extension of Darcy s law, and diffusion is included. Local thermodynamic equilibrium is assumed in all... 

Equation (4) is the fundamental equation of thermodynamic integration (TI). When the derivative of the potential energy with respect to X is known analytically, equation (4) can be applied via a numerical quadrature for a series of X values. In practice, equation (4) has mostly served as the basis of the slow growth (SG) procedure used with MD simulations. Specifically, X is changed incrementally over the full M time steps of the simulation and the average in equation (4) is approximated by a finite difference between the time steps to yield... 

The overall simulation of high-temperature corrosion processes under near-service conditions requires both a thermodynamic model to predict phase stabilities for given conditions and a mathematical description of the process kinetics, i.e. solid state diffusion. Such a simulation has been developed by integrating the thermodynamic program library, ChemApp, into a numerical finite-difference diffusion calculation, InCorr, to treat internal oxidation and nitridation of Ni-base alloys [10]. This simulation was intended to serve as a basis for an advanced computer model for internal oxidation and sulfidation of low-alloy boiler steels. 

Within the context of materials science, and particularly the area of tilloys, the realization that the materials we deal with are invariably of finite size has two different connotations. Finite-size materials that conform to and can be described by the laws of thermodynamics, and those that do not. It is indeed a remarkable fact that nature manages to integrate... 

---