Molecule Types Vibrational Structure Force Field Analysis <Prev Next>

Internal Coordinates

The internal coordinates object allows 3N-6 (or more) internal valence coordinates, R, to be set up which will be used to construct the matrix, B, that transforms between the 3N Cartesian displacements, Δx and R:

R = B Δx

Redundant valence coordinates can be specified (i.e. more than 3N-6), though this imposes some constraints on the possible calculations. The internal coordinates window, shown below for NH3 in a planar configuration, allows the display and editing of the valence coordinates and (optionally) the force field in terms of valence coordinates. To create an internal coordinates object, right click on the Electronic State object object and select "Add New..., Internal Coords". To bring up the window right click on the internal coordinates object and select "View...".


The Controls

The plot at the bottom shows the current internal coordinate, i.e. the selected row in the top left grid. The arrows indicate how far each atom moves, and the atoms involved are linked by thick gray lines. The plot updates semi-automatically; use "Operate, Check" to force an update.

G

PlotPlaneGroup The "Plot Plane" control allows the selection of the plane in which to plot the molecule.
MagnificationSpinThe magnification selects the scale factor to use for the arrows.
StatusBar The status bar shows the latest mouse position; draw a box with the mouse to measure distance and angles.
OperateButton
The operate button brings up a menu of possible operations:
Check
 		Check the contents of the top grids, and update the plot to match
Print
 		Check the contents of the top grids, and display the details of the vibrational mode calculation in the Log Window.
Show VIBCA file
 		Writes a data file for the VIBCA program to the Log Window. See http://www.ifpan.edu.pl/~kisiel/vibr/vibr.htm; note that not all PGOPHER files can be converted.
Calculate f Matrix from l Print the f matrix in the Log Window calculated from the l matrix elements and frequencies given in the vibrational mode objects. This is only possible if there are no redundant internal co-ordinates. If this is not the case, the calculation must be done via symmetry coordinates and the F matrix.
Set f Matrix from l Clear the f matrix grid (top right) and set it from the l matrix elements and frequencies given in the vibrational mode objects. This is only possible if there are no redundant internal co-ordinates. If this is not the case, the calculation must be done via symmetry coordinates and the F matrix.
Calculate f matrix from F matrix Print the f matrix calculated from the F matrix using the following equation:
f = UT F U
Set f matrix from F matrix Clear the f matrix grid (top right) and set it by calculation from the F matrix using the following equation:
f = UT F U
Sort f Matrix by Value Sort the f matrix grid (top right) in descending order of the first column.
Sort f Matrix by Mode Sort the f matrix grid (top right) by mode numbers, with diagonal elements first.

The Force Constant Matrix

FGrid
The grid at the top right allows the force field to be given in terms of valence coordinates:

V = ½ RT f R = ½ Σij fijRiRj

The non zero elements of the f matrix, fij are specified in the grid, with the force constant in the first column, and i and j in the next two columns. (The fourth column is to allow future expansion to cubic and higher terms.). The units of (stretching) force constants are  force constants (for stretches) in millidyne/Å = aJ/Å2 (1 millidyne/Å = 1 aJ/Å2 = 100 J/m2 = 100 N/m).

The force constants can be given as simple numerical values, but can also be given as general expressions involving symbolic values. This is typically more useful, as it allows symmetry constraints to be added naturally. For the NH3 example shown here, note the use of three different force constants, kstretch, kbend and kumb for each of the three different types of internal coordinate. The force constants can be given as general expressions if required; the common usage of scaling ab initio force constants could be done by expressing the constants as kscale*kstretch, and so on, with kscale initialized to 1.When reading the force constant matrix, unknown variables are automatically created as parameters to the current internal coordinates object, and they can then be edited in the Constants Window. If dealing with more than one isotopologue, these variables are better placed in a variables object, (with Global set to true) so the same value can be used in all species.

Note that the force constant grid will expand automatically as atoms are added. The grid can be left blank if the potential energy is to be expressed in terms of symmetry coordinates.

Internal Valence Coordinates

Grid
The grid at the top right allows a set of valence internal coordinates, R, to be specified. The first column dictates the valence coordinate type; a drop down box shows the available types and allows the type to be changed. See Valence Coordinates for more detail on the available types. Each coordinate requires the specification of 2-4 atoms, here labeled k - m, depending on the type of the coordinate.

Internal Coordinates Settings

Active Set true to use this internal coordinates object; there must be no more than one active internal coordinates object at a time.
InPrincipalAxis Toggles selected calculations between the principal axis system and the axis system used in the z matrix specification.