[molpro-user] MCSCF and save,ref= card

Peter REINHARDT reinh at lct.jussieu.fr
Thu Mar 10 09:58:56 GMT 2005


Dear MOLPRO users and maintainers,
while trying to calculate MCSCF wavefunctions with 26 electrons in 18 
active orbitals I encounter the following strange behaviour of MOLPRO:

when asking for a simple MCSCF with 

! MCSCF septet with many open shells, symmetry B1u
!
multi;
closed, 7,2,4,0,5,2,2,0;
wf,70,5,6;

all is going well with about 250 seconds/iteration. The job is finished 
within half an hour and 4 iterations in D2H.

However, when asking for 
! MCSCF septet with many open shells, symmetry B1u
!
multi;
closed, 7,2,4,0,5,2,2,0;
save,ref=2150;
wf,70,5,6;

for saving the configurations on a permanent file for a restart,
the job takes about 15 hours (!) for the first iteration, doing what?
There are these 'new' lines in the output 

 Maximum number of shells:             5
 Maximum number of spin couplings:   637

 Storage for configuration lists12277771, coupling coefficients7418564, address lists*******

and final configurations (and total energies) are different:

with save,ref=2150 (after 7 iterations):

 CI vector
 ---------
 22+0 22 2+ 2+ 22+0 2 2+ +      0.9533840
 22+0 22 2+ 2+ 22+2 2 0+ +     -0.0748168
 22++ 22 -+ 2+ 22++ 2 -+ +     -0.0571668

 TOTAL ENERGIES                     -2912.42558696

without save,ref=2150:

 CI vector
 ---------
 2220 2+ 2+ 2+ 2220 + 2+ +      0.9471155
 2220 2+ 2+ 2+ 2222 + 0+ +     -0.0733663
 2220 2+ 2+ 2+ 2022 + 2+ +     -0.0533393

 TOTAL ENERGIES                     -2912.43052852

Cpu times are (on a 64bit Opteron 2.4 GHz)


 PROGRAMS   *       MULTI     MULTI   RHF-SCF       INT
 CPU TIMES  *      2149.72      6.71      0.32      2.60


 PROGRAMS   *       MULTI     MULTI   RHF-SCF       INT
 CPU TIMES  *     50168.32      6.28      0.34      2.60


Any hint to this behaviour/problem(?) is highly welcome!

Sincerely yours,
   Peter Reinhardt

--
      v  v  v  v  v  v  v  v  v  v  v  v  v  v  v  v  v  v  v  v  v
     >  Peter Reinhardt              |                             <
     >  Lab. Chimie Theorique        |                             <
     >  Universite Paris VI          |                             <
     >  4 place Jussieu, case 137    |  Tel.: +33(0)1 44 27 9657   <
     >  F -- 75252 Paris, France     |  Fax.:               4117   <
     >          email: Peter.Reinhardt at lct.jussieu.fr              <
     >          Web: http://www.lct.jussieu.fr/reinhardt           <
      ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^





More information about the Molpro-user mailing list