Often, no further input is necessary. By default, the number of electrons is equal to the nuclear charge, the wavefunction is assumed to be totally symmetric (symmetry 1), and the spin multiplicity is 1 (singlet) for an even number of electrons and 2 (doublet) otherwise. The Aufbau principle is used to determine the occupation numbers in each symmetry. Normally, this works well in closed-shell cases, but sometimes wrong occupations are obtained or the wavefunction alternates between different orbital spaces. In such cases, the OCC directive must be used to force convergence to the desired state. The default behaviour can be modified either by options on the command line, or by directives.
In open-shell cases, we recommend to use the WF, OCC, CLOSED, or OPEN cards to define the wavefunction uniquely. Other commands frequently used are START and ORBITAL (or SAVE) to modify the default records for starting and optimized orbitals, respectively. The SHIFT option or directive allows to modify the level shift in the RHF program, and EXPEC to calculate expectation values of one-electron operators (see section 6.13).
molpro@molpro.net