Main /
## Use Of Internal And Symmetry-adapted Internal Coordinates[[Page link?]]By default, force constants obtained in a vibrational frequency calculation are available in both the Cartesian representation and also in terms of the internal coordinates specified in the ZMAT file. On some occasions, it may be desirable to have force constants expressed in other ways, such as in symmetry coordinates. Because of the requirement that a "proper" ZMAT geometry specification must include symmetrically complete sets of internal coordinates, it is sometimes impossible to construct a ZMAT file with the same coordinates that you would like to use for calculating force constants. For example, most documented harmonic force fields for benzene include f To facilitate comparisons of ab initio force constants with those in the literature, the user may specify an additional set or sets of coordinates. These will not be used in the course of the geometry optimization (those specified in the Z-matrix are always used), and are read only at the end of a vibrational frequency calculation. They are then used to construct matrices needed to transform first and second derivatives from the Cartesian representation to the coordinates that are specified. The guidelines for defining alternative coordinate systems are as follows. If the coordinates of interest are simple internal coordinates (bond stretching, bond angle, dihedral angle or angle between a bond and a plane), then only one file is needed. This file is called INTCOORD and has the format described below: Each internal coordinate is defined by two consecutive lines in the file. The first line contains a CHARACTER*1 variable that specifies the coordinate as an internuclear distance (R); the angle ABC defined by vectors from atom B to atom A and from atom B to atom C (A); the dihedral angle ABCD between the projections of the vectors CD and BA onto the plane perpendicular to the BC internuclear vector (D); and the angle between the plane formed by the BC and BD bonds and the bond between atoms A and B (W). The second line is an integer string and specifies the atoms involved in the particular coordinate. The numbering scheme for the atoms corresponds precisely to that in the ZMAT file, except that dummy atoms are not counted. Thus, atom number 6 is the sixth real atom that appears in the ZMAT file. For R, the two integers correspond to the atoms defining the distance; for angle ABC, the integers correspond to atoms A, B and C; for dihedral angle ABCD, the integers correspond to atoms A,B,C and D; and for the "Wilson" angle, the four atoms correspond to atoms C,D,A and B (the order is very important here!). As an example, suppose the H Each composite coordinate is defined by three consecutive lines in the file. The first contains the number of primitive (i.e. internal) coordinates that contribute to the composite coordinate. The second contains a string of integers that defines which of the INTCOORD file coordinates are used; the indices correspond with the order of definition in the INTCOORD file. The third line specifies the coefficients that are attached to each primitive coordinate in the composite coordinate. These are automatically normalized by the program, so it is not possible to use unnormalized coordinates at the present time. An example which helps to clarify the SYMCOORD file protocol is given below, which defines the usual symmetry coordinates for water [S A comment for internal coordinate aficionados: If only the INTCOORD file is supplied, the program is capable of handling redundancies in the coordinate system. It uses the eigenvectors of the |

Page last modified on May 17, 2013, at 09:26 PM

CFOUR is partially supported by the U.S. National Science Foundation.