19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

10 MOLECULAR GEOMETRY 64<br />

geometry specification is given, but usually this does not need to be done, as a geometry specification<br />

where the first line is a single integer will be recognized as XYZ format, as will the case<br />

of the first line consisting of a chemical symbol followed by three cartesian coordinates. The<br />

geometry block should then contain the cartesian coordinates in Minnesota Computer Centre,<br />

Inc. XYZ format. Variable names, and in general expressions that are linear in all dependent<br />

variables, may be used as well as fixed numerical values. Non-linear expressions should not be<br />

used, because the geometry optimization module is unable to differentiate them.<br />

The XYZ file format consists of two header lines, the first of which contains the number of<br />

atoms, and the second of which is a title. The remaining lines each specify the coordinates of<br />

one atom, with the chemical symbol in the first field, and the x, y, z coordinates following. A<br />

sequence number may be appended to the chemical symbol; it is then interpreted as the atomic<br />

group number, which can be used when different basis sets are wanted for different atoms of the<br />

same kind. The basis set is then specified for this group number rather than the atomic symbol.<br />

As a further extension, the first two header lines can be omitted.<br />

Note that for XYZ input the default is not to reorient the molecule. Orientation can be forced,<br />

however, by the MASS or CHARGE options on the ORIENT directive.<br />

geomtyp=xyz<br />

geometry={<br />

3 ! number of atoms<br />

This is an example of geometry input for water with an XYZ file<br />

O ,0.0000000000,0.0000000000,-0.1302052882<br />

H ,1.4891244004,0.0000000000, 1.0332262019<br />

H,-1.4891244004,0.0000000000, 1.0332262019<br />

}<br />

hf<br />

http://www.molpro.net/info/current/examples/h2o_xyzinput.com<br />

The XYZ format is specified within the documentation distributed with MSCI’s XMol package.<br />

Note that <strong>MOLPRO</strong> has the facility to write XYZ files with the PUT command (see section 10.3).<br />

10.2 Symmetry specification<br />

If standard Z-matrix input is used, <strong>MOLPRO</strong> determines the symmetry automatically by default.<br />

However, sometimes it is necessary to use a lower symmetry or a different orientation<br />

than obtained by the default, and this can be achieved by explicit specification of the symmetry<br />

elements to be used, as described below.<br />

Generating symmetry elements, which uniquely specify the point group, can be specified on the<br />

SYMMETRY directive. This must be given before the geometry block. Each symmetry directive<br />

only affects the subsequent geometry block; after a geometry block has been processed, the<br />

defaults are restored. Note that the specification of symmetry elements inside the geometry<br />

block is no longer allowed.<br />

The dimension of the point group is 2**(number of fields given). Each field consists of one or<br />

more of X, Y, or Z (with no intervening spaces) which specify which coordinate axes change<br />

sign under the corresponding generating symmetry operation. It is usually wise to choose z to<br />

be the unique axis where appropriate (essential for C 2 and C 2h ). In that case, the possibilities<br />

are:<br />

(null card)<br />

Z<br />

C 1 (i.e., no point group symmetry)<br />

C s

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!