Energy quadratic approximation

There are several reasons that Newton-Raphson minimization is rarely used in mac-romolecular studies. First, the highly nonquadratic macromolecular energy surface, which is characterized by a multitude of local minima, is unsuitable for the Newton-Raphson method. In such cases it is inefficient, at times even pathological, in behavior. It is, however, sometimes used to complete the minimization of a structure that was already minimized by another method. In such cases it is assumed that the starting point is close enough to the real minimum to justify the quadratic approximation. Second, the need to recalculate the Hessian matrix at every iteration makes this algorithm computationally expensive. Third, it is necessary to invert the second derivative matrix at every step, a difficult task for large systems. 

The superscript (2) marks the quadratic approximation of the electrostatic energy. 

The concept of dipole hardness permit to explore the relation between polarizability and reactivity from first principles. The physical idea is that an atom is more reactive if it is less stable relative to a perturbation (here the external electric field). The atomic stability is measured by the amount of energy we need to induce a dipole. For very small dipoles, this energy is quadratic (first term in Equation 24.19). There is no linear term in Equation 24.19 because the energy is minimum relative to the dipole in the ground state (variational principle). The curvature hi of E(p) is a first measure of the stability and is equal exactly to the inverse of the polarizability. Within the quadratic approximation of E(p), one deduces that a low polarizable atom is expected to be more stable or less reactive as it does in practice. But if the dipole is larger, it might be useful to consider the next perturbation order ... 

Their conclusion tended to the negative on this. They also concluded that quadratic approximations for energy variations were not applicable. 

Fig. 3.2 Changes in bond lengths or in bond angles result in changes in the energy of a molecule. Such changes are handled by the Es,re,Ch and Ebend terms in the molecular mechanics forcefield. The energy is approximately a quadratic function of the change in bond length or angle...

Let us first define the external MEC in M, consisting of m atoms. Consider the global equilibrium of M in contact with a hypothetical electron reservoir (r) fi0 = fj1 where fi= fi, the chemical potential of r. Let z = N — N° = d/V denotes the vector of a hypothetical AIM electron-population displacements from their equilibrium values N°. Since d/V = - d/Vr, the assumed equilibrium removes the first-order contribution to the associated change due to z in the energy, = M + , of the combined (closed) system (Mir) moreover, taking into account the infinitely soft character of a macroscopic reservoir, the only contribution to the energy change in the quadratic approximation is ... 

Thus, the corresponding expression for the energy change in the quadratic approximation becomes ... 

The last terms f(bare) fs(bare) j in Eq. (7) are the bare contributions to the surface excess free energy. Since we consider here only symmetric walls, we assume the same functional form for both walls at z=0 and at z=D, assuming the quadratic approximation [11] for simplicity,... 

From the frontier orbital energies an approximated absolute hardness is also obtained. In fact, under a finite difference approximation and a quadratic dependence of the energy on the number of electrons, absolute hardness is defined as... 

Using intersection-adapted coordinates, the quadratic approximation, in other words the local harmonic approximation, of the adiabatic energy difference for a finite displacement around Qo reads thus... 

The component of displacement of the F ion is denoted Uia- The rationale for our strategy is cast in geometric terms in fig. 5.2 which shows fhe quadratic approximation to a fully nonlinear energy surface in the neighborhood of the minimum for an idealized two-dimensional configuration space. 

Fig. 5.2. Nonlinear potential energy V(xi, X2) and corresponding quadratic approximation to that potential. The cutaway in the upper left hand corner shows one quadrant of the two potentials to illustrate their correspondence in the limit of small displacements.

Vibrational spectroscopies are particularly useful for the analysis of the adsorbed layers on metallic particles. Among them, infrared spectroscopy is of widespread use and provides a powerful tool in the study of metal-based catalysts under reaction conditions. Under the approximation of vibrational and rotational coordinate separation, the vibrational wavefunction by is a function of the internal coordinates (Qk) and is a solution of the vibrational hamiltonian. Assuming a quadratic approximation of the potential energy in terms of the internal coordinates, then ... 

In a pivotal development. Miller, Handy and Adams [12] derived the classical Hamiltonian for a simple potential based on the MEP. The idea of the reaction path Hamiltonian is, conceptually, to consider the potential as a trough or as a stream bed along with 3N-7 harmonic walls that are free to close in or widen out as one proceeds along the trough. The potential energy surface is approximated as the potential energy of the MEP Vo(s) plus a quadratic approximation to the energy in directions perpendicular to the MEP,... 

Nichols J A, Taylor H, Schmidt P P and Simons J 1990 Walking on potential energy surfaces J. Chem. Phys. 92 340-6 Simons J and Nichols J 1990 Strategies for walking on potential energy surfaces using local quadratic approximations... 

This equation can be solved in terms of the eigenvectors of F(rQ) to advance one step forward from 5( o) on the MEP. Much larger steps can be taken using this local quadratic approximation method [41,44]. Even more substantial increases in the step size have been achieved by accounting approximately for third derivatives of the energy along the path (although only ab initio second derivatives are actually computed [45]). 

---