manual   quickstart   instguide   update   basis

Next: 7.4 MCSCF with selected Up: 7 MCSCF and CASSCF Previous: 7.2 Restricted Active Space   Contents

## 7.3 State-averaged MCSCF

In order to compute excited states it is usually best to optimize the energy average for all states under consideration. This avoids root-flipping problems during the optimization process and yields a single set of compromise orbitals for all states.

The number of states to be optimized in a given symmetry is specified on a state directive, which must follow directly after the wf directive, e.g.,

wf,16,1,0;state,2 !optimize two states of symmetry 1

It is also possible to optimize states of different symmetries together. In this case several wf / state directives can follow each other, e.g.,

wf,16,1,0;state,2 !optimize two states of symmetry 1
wf,16,2,0;state,1 !optimize one states of symmetry 2

etc. Optionally also the weights for each state can be specified, e.g.

wf,16,1,0;state,2;weight,0.2,0.8   !optimize two states of symmetry 1
!first state has weight 0.2,
!second state weight 0.8


By default, the weights of all states are identical, which is normally the most sensible choice. The following example shows a state-averaged calculation for , in which the valence states (, ) are treated together.

***,O2
print,basis,orbitals
geometry={                 !geometry specification, using z-matrix
o1
o2,o1,r
}

r=2.2 bohr                 !bond distance

basis=vtz                  !triple-zeta basis set

{hf                        !invoke RHF program
wf,16,4,2                  !define wavefunction: 16 electrons, symmetry 4, triplet
occ,3,1,1,,2,1,1           !number of occupied orbitals in each symmetry
open,1.6,1.7}              !define open shell orbitals

{casscf                    !invoke CASSCF program
wf,16,4,2                  !triplet Sigma-
wf,16,4,0                  !singlet delta (xy)
wf,16,1,0}                 !singlet delta (xx - yy)


Note that averaging of states with different spin-multiplicity, as in the present examples, is possible only for CASSCF, but not for more restricted RASSCF or MCSCF, wavefunctions.

Next: 7.4 MCSCF with selected Up: 7 MCSCF and CASSCF Previous: 7.2 Restricted Active Space   Contents

manual   quickstart   instguide   update   basis

molpro@molpro.net 2013-06-18