19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

19 THE MCSCF PROGRAM MULTI 135<br />

transforming the final orbitals to natural orbitals (the first order density matrix is diagonalized),<br />

to pseudo-canonical orbitals (an effective Fock-operator is diagonalized), or of localizing the<br />

orbitals.<br />

19.5.1 Defining the starting guess<br />

START,record,[options];<br />

record: dump record containing starting orbitals. As usual, record has the form irec.ifil, where<br />

irec is the record number (e.g., 2140), and ifil the file number (usually 2). The options can be<br />

used to select orbitals of a specific type; for details, see section 4.11.<br />

If this card is missing, the program tries to find suitable starting orbitals as follows:<br />

First:<br />

Second:<br />

Third:<br />

Try to read orbitals from the record specified on the ORBITAL<br />

card (or the corresponding default, see ORBITAL). All files are<br />

searched.<br />

Try to find orbitals from the most recent MCSCF calculation.<br />

All files are searched.<br />

Try to find orbitals from the most recent SCF calculation. All<br />

files are searched.<br />

If no orbitals are found, a starting orbital guess is generated.<br />

It is often useful to employ MCSCF orbitals from a neighbouring geometry as starting guess<br />

(this will happen automatically if orbitals are found, see the above defaults). Note, however,<br />

that frozen-core orbitals should always be taken from an SCF or MCSCF calculation at the<br />

present geometry and must be specified separately on the FROZEN card. Otherwise the program<br />

is likely to stop with error “non-orthogonal core orbitals”. The program remembers where to<br />

take the core orbitals from if these have been specified on a FROZEN card in a previous MCSCF<br />

calculation.<br />

19.5.2 Rotating pairs of initial orbitals<br />

ROTATE,orb1.sym,orb2.sym,angle<br />

Performs a 2×2 rotation of the initial orbitals orb1 and orb2 in symmetry sym by angle degrees.<br />

With angle=0 the orbitals are exchanged. ROTATE is meaningful only after the START card.<br />

See MERGE for other possibilities to manipulate orbitals.<br />

19.5.3 Saving the final orbitals<br />

ORBITAL,record.file<br />

The orbitals are dumped to record record.file. Default for record is 2140 and file=2. This default<br />

record number is incremented by one for each subsequent MCSCF calculation in the same job<br />

(see section 4.11). Therefore, if several different MCSCF calculations at several geometries<br />

are performed in one job, each MCSCF will normally start with appropriate orbitals even if no<br />

ORBITAL or START card is present.<br />

The ORBITAL card can be omitted if a NATORB, CANORB or LOCORB card is present, since<br />

orb can also be specified on these cards (the same defaults for orb as above apply in these cases).

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

Saved successfully!

Ooh no, something went wrong!