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.

16 THE SCF PROGRAM 101<br />

16.1.8 Printing options<br />

PRINT|ORBPRINT=value<br />

DEBUG=value<br />

Number of virtual orbitals to be printed. If value=0, the occupied<br />

orbitals are printed.<br />

Option for debug print.<br />

16.2 Defining the wavefunction<br />

The number of electrons and the total symmetry of the wavefunction are specified on the WF<br />

card:<br />

WF,elec,sym,spin<br />

where<br />

elec<br />

sym<br />

spin<br />

is the number of electrons<br />

is the number of the irreducible representation<br />

defines the spin symmetry, spin = 2 ∗ S (singlet=0, doublet=1,<br />

triplet=2 etc.)<br />

Note that these values take sensible defaults if any or all are not specified (see section 10.2).<br />

16.2.1 Defining the number of occupied orbitals in each symmetry<br />

OCC,n 1 ,n 2 ,...,n 8<br />

To avoid convergence problems in cases with high symmetry, this card should be included whenever<br />

the occupation pattern is known in advance. n i is the number of occupied orbitals in the<br />

irreducible representation i. The total number of orbitals must be equal to (elec+spin)/2 (see WF<br />

card).<br />

16.2.2 Specifying closed-shell orbitals<br />

CLOSED,n 1 ,n 2 ,...,n 8<br />

This optional card can be used in open-shell calculations to specify the number of closed-shell<br />

orbitals in each symmetry. This makes possible to force specific states in the absence of an<br />

OPEN card.<br />

16.2.3 Specifying open-shell orbitals<br />

OPEN,orb 1 .sym 1 ,orb 2 .sym 2 ,...,orb n .sym n<br />

This optional card can be used to specify the singly occupied orbitals. The number of singly<br />

occupied orbitals must be equal to spin, and their symmetry product must be equal to sym (see<br />

WF card). If the OPEN card is not present, the open shell orbitals are selected automatically.<br />

The algorithm tries to find the ground state, but it might happen that a wrong state is obtained<br />

if there are several possibilities for distributing the open shell electrons among the available<br />

orbitals. This can also be avoided using the CLOSED card. If orb i .sym is negative, this orbital<br />

will be occupied with negative spin (only allowed in UHF).

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

Saved successfully!

Ooh no, something went wrong!