Minimum energy coordinates components

To coordinate components, the generalized flows and the thermodynamic forces can be used to define the trajectories of the evolution of nonequilibriun systems in time. A trajectory specifies the curve represented by the flow and force components as a function of time in the flow-force space. A useful trajectory can be found and analyzed by a variation principle. In thermodynamics, the variation principles lead to the least energy dissipation and minimum entropy generation at steady states. According to the most general evolutionary criterion, open chemical reaction systems are dissipative, and evolve toward an asymptotic state in time. 

A careful distinction must be drawn between transition states and intermediates. As noted in Chapter 4, an intermediate occupies a potential energy minimum along the reaction coordinate. Additional activation, whether by an intramolecular process (distortion, rearrangement, dissociation) or by a bimolecular reaction with another component, is needed to enable the intermediate to react further it may then return to the starting materials or advance to product. One can divert an intermediate from its normal course by the addition of another reagent. This substance, referred to as a trap or scavenger, can be added prior to the start of the reaction or (if the lifetime allows) once the first-formed intermediate has built up. Such experiments are the trapping experiments referred to in Chapters 4 and 5. 

We begin by defining a set of N atoms and writing their Cartesian coordinates as a single vector with 3N components, r = (r, ..., ro,N). If locating the atoms at ro is a local minimum in the energy of the atoms, then it is convenient to define new coordinates x = r ro. The Taylor expansion of the atom s energy about the minimum at ro is, to second order,... 

Table 33.1 Components of the SAPT interaction energy of the water dimer for four intermolecular separations R, computed in the [5s3p3d2flglh/3s3pldlf ] basis taken from Ref. [103], containing the 3s2pldlf set of midbond functions. The angular coordinates are close to those of the global minimum configuration. All energies in kcal/mol...

I is the identity matrix. The six first derivatives of the energy with respect to the strain components e, measure the forces acting on the unit cell. When combined with the atomic coordinates we get a matrix with 3N - - 6 dimensions. At a minimum not only should there be no force on any of the atoms but the forces on the unit cell should also be zero. Application of a standard iterative minimisation procedure such as the Davidon-Fletcher-Powell method will optimise all these degrees of freedom to give a strain-free final structure. In such procedures a reasonably accurate estimate of the initial inverse Hessian matrix is usually required to ensure that the changes in the atomic positions and in the cell dimensions are matched. 

---