14.1 Analysis of Normal Modes in Terms of Internal Coordinates

A note in advance: The analysis of normal modes can (at nearly no computational cost) always be redone as long as you keep a copy of the file hessian.

A general prerequisite for this option is that you have defined a set of non-redundant coordinates for all 3N-6 (3N-5) degrees of freedom of your molecule. To make sure that this is the case, you should switch off redundant coordinates (currently, this is only possible by manually removing the data group $redundant and also removing the entry redundant on in $optimize). Run define to generate non-redundant coordinates by using the iaut command in the internal coordinate menu (or by creating them manually via idef). We recommend to use the irem command first to delete all previous definitions of internal coordinates. See Section 4 for further details. If the molecules point group is not C1, define will set some of the coordinate to status d (display) or i (ignore). Use the ic command to change all coordinates to k. You can also achieve this by editing in the $intdef data-group manually.

The analysis in internal coordinates is switched on by adding a line in the data-group $drvopt that has the following syntax:

analysis [only] intcoord [print print-level]

Keywords in square brackets are optional. If only is added, the program assumes that the file hessian exists and runs only the analysis part of aoforce. The program will give the following output (controlled by the print level given in parenthesis):

Note that for large molecules or complicated topologies the B-matrix (that is used to transform from Cartesian coordinates into internal coordinates and vice versa) may become singular. In this case only the normal modes in the internal coordinate basis can be listed.