Using Linear Combinations Of Original ZMAT Coordinates
A type of calculation that might arise is one in which it is desirable to use combinations of the ZMAT file internal coordinates. For example, suppose that you wanted to optimize the structure of planar vinylidene (CCH2) subject to the constraint that the two CCH bond angles (those equilvalent in the C2v equilibrium structure differ by 10 degrees. While it might be possible to obtain such a result by doing several different calculations, there are features in CFOUR that allow this to be done easily. An area of application in which such calculations are commonly needed is that associated with developing potential energy surfaces for nuclear motion calculations. The new NITROGEN program, developed by Changala, is the choice of the CFOUR developers for most variational nuclear motion calculations, and this code often makes use of such linear combinations of coordinates.
To specify an alternative set of internal coordinates to be used in an optimizations which are linear transformations of those given in the ZMAT file, include the string
in your ZMAT file (the name is motivated by NITROGEN). The syntax used to specify the coordinates and the associated optimization is best communicated by example. For the vinylidene calculation mentioned above, the body of the usual ZMAT specification might be
The %nitroq field looks like this: