Quadratic equation, energy surfaces

The method of moments of coupled-cluster equations (MMCC) is extended to potential energy surfaces involving multiple bond breaking by developing the quasi-variational (QV) and quadratic (Q) variants of the MMCC theory. The QVMMCC and QMMCC methods are related to the extended CC (ECC) theory, in which products involving cluster operators and their deexcitation counterparts mimic the effects of higher-order clusters. The test calculations for N2 show that the QMMCC and ECC methods can provide spectacular improvements in the description of multiple bond breaking by the standard CC approaches. 

In the case of potentials having a lower permutation symmetry, the coordinates (i = 1,2,3) may still be used, but other terms besides those of equation (15) are obviously allowed for the representation. For example, for AB2-type systems, there are two linear terms, two quadratic terms, and one cubic term with C2 symmetry which may be used for the representation of the potential-energy surface, namely29... 

As we noted above, the kinetic energy is positive definite. Furthermore, it is quadratic in the momenta. As a consequence, we can reduce the search for points of stationary flow in phase space to one of finding the stationary points of the potential energy surface. To see how this comes about, consider the Hamilton s equations for the three velocities... 

The common feature observed in both DFT and GCMC simulations is that these results overpredict the amount adsorbed in the reduced pressure region greater than about 0.2. This seems to indicate that the fluid—fluid interaction energy is overestimated in the presence of a solid surface, and therefore the usual assumption of pairwise additivity of fluid-fluid and solid-fluid potential energies is questionable. One way of resolving this issue is the application of the following quadratic equation for the potential of one molecule [83] ... 

In the case where we only take into account the quadratic part of the normal form, or, equivalently, if we linearize Hamilton s equations, we have Kcnf I, h, , Jd) = XI + Y1=2 ( kjk and the energy surface Kcnf(0, ]i,. .., ]d) = E encloses a simplex in (/2,..., Jd) whose volume leads to the well-known result [45]... 

Exercise C.2 Consider case 3 of Exercise C.lb. Evaluate the gradient and the Hessian at c = 3.3, y = 1.8 and solve the Newton-Raphson equation of Exercise C.la for xi, yi). Comment on this. In practical cases, the quantum mechanical potential energy surface is not quadratic, especially for regions as far away from the neighborhood of the minimum as (jc, y) in this example. 

This PSP formalism can be used for both concentrated as well as infinitely dilute systems, for small molecules as well as for high polymers and copolymers. It may, in particular, be used for the surface characterization of polymer and other solid surfaces. In fact, it may be used in order to obtain the corresponding surface energy components, Ya. and y, of the solid surfaces [42]. The widely used van der Waals-Lifshitz surface tension component is nearly the sum of Yh- and while the contact angle 0 of a liquid L on a polymer surface P is given by the quadratic equation [42] ... 

Standard numerical methods such as second-order Runge-Kutta could be used, but a more effective approach is to expand the potential energy surface in equation (2) to second order and integrate the resulting expression from Xi to 3c,+i. This yields the local quadratic approximation (LQA) of Page and Mclver " which is an explicit second-order method. 

The adiabatic potential-energy surfaces are obtained by diagonalizing the Hamiltonian of equations (27) in the fixed-nuclei limit (Tn = 0). For vanishing quadratic coupling (g = 0), this yields... 

The analysis of Section 2.4 may now be taken over in its entirety. If the initial point is near a stationary point on the energy surface, and the surface is locally quadratic (at least in good approximation), then the stationary point may be reached by choosing the parameters so that (2.4.20) is satisfied. More commonly, when the parameters are real, it is sufficient to solve the reduced equation (2.4.21), which may be written... 

Once the tensor components have been evaluated, for some given point on the energy surface, the nearest stationary point (in quadratic approximation) can be found by solving the usual equations (cf. (2.4.20)). When the matrices are blocked according to orbital and coefficient indices, these equations may be written... 

Newton s method is based on a local quadratic model of the energy surface. The orbital rotations generated by this method therefore behave incorrectly for large rotations. In particular, the Newton equations are not periodic in the orbital-rotation parameters as we would expect fix)m a consideration of the global behaviour of the energy function (10.1.22). The one-electron approximation of the SCF method, by contrast, is correct only to zero order in the rotations, but - provided the effective Hamiltonian (i.e. the Fock operator) is a reasonable one - it exhibits the correct global behaviour. In particular, it is periodic in the orbital rotations. 

---