Geometry optimization for molecules

For such applications of classical optimization theory, the data on energy and gradients are so computationally expensive that only the most efficient optimization methods can be considered, no matter how elaborate. The number of quantum chemical wave function calculations must absolutely be minimized for overall efficiency. The computational cost of an update algorithm is always negligible in this context. Data from successive iterative steps should be saved, then used to reduce the total number of steps. Any algorithm dependent on line searches in the parameter hyperspace should be avoided. 

An equilibrium state is defined for generalized coordinates such that the total energy E( q ) is minimized. The energy gradients p = are forces with reversed sign. Coordinate and gradient displacements from m successive iterations are saved as m x n column matrices, Aq and Ap, respectively. The Hessian matrix is Ftj =, , such that Newton s extrapolation formula is Aq = —Gp°, where G = F 1. 

In quasi-Newton methods, a parametrized estimate of F or G is used initially, then updated at each iterative step. Saved values of q, p can be used in standard algorithms such as BFGS, described above. 

---

G. E. Chudinov and D. V. Napolov, Chem. Phys. Lett., 201, 250 (1993). A New Method of Geometry Optimization for Molecules in Solution in the Framework of the Born-Kirkwood-Onsager Approach. 

Among the various types of G derivatives that can be computed with the PCM program, a special importance is taken by the derivatives with respect to nuclear coordinates. Without analytical expressions of the energy gradient it is not possible to perform studies of geometry optimizations for molecules of medium size, and without analitycal formulation of the Hessian both the characterization of the saddle points and the determination of vibrational frequencies become quite expensive. 

Eckert F, Pulay P and Werner H-J 1997 Ab initio geometry optimization for large molecules J. Comput. Chem. 18 1473... 

In principle, an accurate calculation of these barriers should require full geometry optimization for the planar molecule. This is no easy task, however, and the result is very... 

With all those assumptions and limitations in mind we can analyze the number of matrix-vector products necessary to perform the geometry optimization for the three model molecules, reported in Table 1.3. 

YETI [198] is a force field designed for the accurate representation of nonbonded interactions. It is most often used for modeling interactions between biomolecules and small substrate molecules. The molecular geometry optimization for the component molecules is not previewed so that it has been obtained from some other source, such as AMBER. Then YETI is used to model the docking. 

For the first time, analytic gradients have been devised and implemented for the DKH approach by Nasluzov and Rosch (1996), which is a prerequisite for efficient geometry optimizations of molecules with more than just a few atoms. 

The analysis of the excited state properties of Ml as a function of the polarity factor k s), was performed with the geometry optimized for the ground state and adopting, for the sake of comparison, both the solvaton set corresponding to the net charges of the unsolvated molecule (calc. A) and the self-... 

What is correct in the Grotthuss mechanism is the stepwise transfer mechanism. If the mechanism had been one of ordinary diffusion of a given (classical) proton through water, there would have been all kinds of barriers. The important thing is that the proton that transfers is not the same proton. In Figure 9.11, a cluster of water molecules is geometry optimized. The molecule at the center is a hydronium molecule. The distance to one of the other molecules is very short (2.14 A) and almost certainly barrierless. From this molecule, there is a chain of hydrogen bonds suitable for successive PT. 

The structures and stabilities of the two NH4 isomers have been studied theoretically with various quantum-chemical methods. Geometry optimizations for the H NH3 species by the CEPA [9], MP2 and SDQ-CI [10], MP2 [11, 12], SD-CI [13], and MBPT(2) [14] methods all agreed on a structure of Cg symmetry in which the H" ion is bound nearly linearly to one of the N-H bonds of the slightly distorted NH3 molecule with a small tilt (maximal 15 ) toward the other hydrogen atoms. The relatively large H-H distance (around 2 A) Indicates a bond arising from dipole-dipole attraction. Values between 25 and 63 kJ/mol... 

---