Coordinates internal

In linear XYZ (c), bent XY2 (b), and pyramidal XY, (c) molecules, the number of internal coordinates is the same as the number of normal vibrations. In a nonplanar X2Y2 molecule (d) such as H2O2, the number of internal coordinates is the same as the number of vibrations if the twisting angle around the central bond (At) is considered. In a tetrahedral XY4 molecule (c), however, the number of internal coordinates exceeds the number of normal vibrations by one. This is due to the fact that the six angle coordinates around the central atom are not independent of each other, that is, they must satisfy the relation 

This is called a redundant condition. In planar XY3 molecule (/), the number of internal coordinates is seven when the coordinate A6, which represents the deviation from planarity, is considered. Since the number of vibrations is six, one redundant condition such as 

The redundant conditions are more complex in ring compounds. For example, the number of internal coordinates in a triangular X3 molecule (h) exceeds the number of vibrations by three. One of these redundant conditions (A species) is 

The other two redundant conditions (JE species) involve bond stretching and angle deformation coordinates such as 

The procedure for finding the number of normal vibrations in each species was described in Sec. 1-7. This procedure is, however, considerably simplified if internal coordinates are used. Again consider a pyramidal XY3 molecule. Using the internal coordinates shown in Fig. 1-1 Ic, we can write the representation for the C3 operation as 

In principle, is a known function of p and therefore also of Spjjj, through the equilibrium conditions 

Note that, in stead of the original F, which is singular or 

As the arbitrary initial conformation x we may choose the conformation we know from minimisation, x (p)  

Approximating P computed at equilibrium with P at the initial conformation, just as before, we get 

This means that both the gradient and the Hessian matrix are calculated at the equilibrium conformation as found with unchanged energy parameters, but now with one parameter changed at a time. The set of 

---

Reaction is assumed to have occurred if a particular internal coordinate q, such as a bond length, attains a... 

The transfonnation matrix L is obtained from a nonnal-mode analysis perfonned in internal coordmates [59, ]. Thus, as the evolution of the nonnal-mode coordinates versus time is evaluated from equation (A3.12.49), displacements in the internal coordinates and a value for q are found from equation (A3.12.50). The variation in q with time results from a superposition of the nonnal modes. At a particular time, the... 

Electronic spectra are almost always treated within the framework of the Bom-Oppenlieimer approxunation [8] which states that the total wavefiinction of a molecule can be expressed as a product of electronic, vibrational, and rotational wavefiinctions (plus, of course, the translation of the centre of mass which can always be treated separately from the internal coordinates). The physical reason for the separation is that the nuclei are much heavier than the electrons and move much more slowly, so the electron cloud nonnally follows the instantaneous position of the nuclei quite well. The integral of equation (BE 1.1) is over all internal coordinates, both electronic and nuclear. Integration over the rotational wavefiinctions gives rotational selection rules which detemiine the fine structure and band shapes of electronic transitions in gaseous molecules. Rotational selection rules will be discussed below. For molecules in condensed phases the rotational motion is suppressed and replaced by oscillatory and diflfiisional motions. 

These electronic energies dependence on the positions of the atomic centres cause them to be referred to as electronic energy surfaces such as that depicted below in figure B3.T1 for a diatomic molecule. For nonlinear polyatomic molecules having atoms, the energy surfaces depend on 3N - 6 internal coordinates and thus can be very difficult to visualize. In figure B3.T2, a slice tln-oiigh such a surface is shown as a fimction of two of the 3N - 6 internal coordinates. 

In this way the optimization can be cast m temis of the original coordinate set, including the redundancies. Exactly the same transfomiations between Cartesian and internal coordinate quantities hold as for the non-redundant case (see the next section), but with the generalized inverse replacing the regular inverse. 

The redundant optimization scheme [53] can be applied to natural internal coordinates, which are sometimes redundant for polycyclic and cage compounds. It can also be applied directly to the underlying primitives. 

Table B3.5.1 Number of cycles to converge for geometry optimizations of some typical organic molecules usmg Cartesian, Z-matrix and delocalized internal coordinates. ...

Thus the transfonnation matrix for the gradient is the inverse transpose of that for the coordinates. In the case of transfonnation from Cartesian displacement coordmates (Ax) to internal coordinates (Aq), the transfonnation is singular becanse the internal coordinates do not specify the six translational and rotational degrees of freedom. One conld angment the internal coordinate set by the latter bnt a simpler approach is to rise the generalized inverse [58]... 

Wlien working with any coordinate system other than Cartesians, it is necessary to transfonn finite displacements between Cartesian and internal coordinates. Transfomiation from Cartesians to internals is seldom a problem as the latter are usually geometrically defined. However, to transfonn a geometry displacement from internal coordinates to Cartesians usually requires the solution of a system of coupled nonlinear equations. These can be solved by iterating the first-order step [47]... 

By combining the Lagrange multiplier method with the highly efficient delocalized internal coordinates, a very powerfiil algoritlun for constrained optimization has been developed [ ]. Given that delocalized internal coordinates are potentially linear combinations of all possible primitive stretches, bends and torsions in the system, cf Z-matrix coordinates which are individual primitives, it would seem very difficult to impose any constraints at all however, as... 

Baker J and Chan F 1996 The location of transition states a comparison of Cartesian, Z-matrix, and natural internal coordinates J. Comput. Chem. 17 888... 

Baker J, KInghorn D and Pulay P 1999 Geometry optimization In delocalized Internal coordinates An efficient quadratically scaling algorithm for large molecules J. Chem. Phys. 110 4986... 

FogarasI G, Zhou X, Taylor P W and Pulay P 1992 The calculation of ab initio molecular geometries efficient optimization by natural Internal coordinates and empirical correction by offset forces J. Am. 

Pye C C and Poirier R A 1998 Graphical approach for defining natural Internal coordinates J. Comput. Chem. 19 504... 

Pulay P and FogarasI G 1992 Geometry optimization In redundant Internal coordinates J. Chem. Phys. 96 2856... 

Peng C, Ayala P Y, Schlegel H B and Frisch M J 1996 Using redundant Internal coordinates to optimize equilibrium geometries and transition states J. Comput. Chem. 17 49... 

Baker J, KessI A and Delley B 1996 The generation and use of delocalized Internal coordinates In geometry optimization J. Chem. Phys. 105 192... 

The three internal coordinates aie expressed as combinations of squares of the interparticle distances ... 

The full dynamical treatment of electrons and nuclei together in a laboratory system of coordinates is computationally intensive and difficult. However, the availability of multiprocessor computers and detailed attention to the development of efficient software, such as ENDyne, which can be maintained and debugged continually when new features are added, make END a viable alternative among methods for the study of molecular processes. Eurthemiore, when the application of END is compared to the total effort of accurate determination of relevant potential energy surfaces and nonadiabatic coupling terms, faithful analytical fitting and interpolation of the common pointwise representation of surfaces and coupling terms, and the solution of the coupled dynamical equations in a suitable internal coordinates, the computational effort of END is competitive. 

To obtain the Hamiltonian at zeroth-order of approximation, it is necessary not only to exclude the kinetic energy of the nuclei, but also to assume that the nuclear internal coordinates are frozen at R = Ro, where Ro is a certain reference nucleai configuration, for example, the absolute minimum or the conical intersection. Thus, as an initial basis, the states t / (r,s) = t / (r,s Ro) are the eigenfunctions of the Hamiltonian s, R ). Accordingly, instead of Eq. (3), one has... 

Atoms not explicitly included in the trajectory must be generated. The position at which an atom may be placed is in some sense arbitrary, the approach being analogous to the insertion of a test particle. Chemically meaningful end states may be generated by placing atoms based on internal coordinates. It is required, however, that an atom be sampled in the same relative location in every configuration. An isolated molecule can, for example, be inserted into... 

---