[molpro-user] Can you average over ionised/attached states using SA-CASSCF in Molpro?

[Seth Olsen]

I thought I had done this before in the past, but I can’t find it, so I’ll just ask you: is it possible to average over states with different electron number using Molpro’s SA-CASSCF?  For example, the mcscf deck was written with the intention of averaging the ground state of a 90-electron molecule with it’s electron-attached state:

{mcscf
occ,49
closed,43
config,csf
wf,90,1,0
state,1
wf,91,1,1
state,1
natorb,2141.2;dm
}

It ends up failing with the message:

1PROGRAM * MULTI (Direct Multiconfiguration SCF)       Authors: P.J. Knowles, H.-J. Werner (1984)     S.T. Elbert (1988)


 Number of closed-shell orbitals:  43 (  43 )
 Number of active  orbitals:        6 (   6 )
 Number of external orbitals:      95 (  95 )

 State symmetry 1

 Number of electrons:     4    Spin symmetry=Singlet   Space symmetry=1
 Number of states:        1
 Number of CSFs:        105

 State symmetry 2

 Number of electrons:     5    Spin symmetry=Doublet   Space symmetry=1
 Number of states:        1
 Number of CSFs:        210

 ? Error
 ? illegal nelec
 ? The problem occurs in cicon

 GLOBAL ERROR fehler on processor   0



Can this be done?  I thought to ask before wrestling with it further…

Cheers,
Seth
===========================
Seth Olsen, PhD.
Honorary Fellow
School of Mathematics & Physics
The University of Queensland
QLD 4072  Australia
Ph: +61 7 3365 2816
===========================
A PGP public key for this address has been uploaded to the key servers.



===========================
Seth Olsen, PhD.
Honorary Fellow
School of Mathematics & Physics
The University of Queensland
QLD 4072  Australia
Ph: +61 7 3365 2816
===========================
A PGP public key for this address has been uploaded to the key servers.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20170828/b857eee8/attachment-0001.html>