[molpro-user] maxiter puzzle in multi

Haitao Ma mht at mrdlab.icas.ac.cn
Thu Apr 21 02:44:13 BST 2005


Dear Molpro-user,

    I got a puzzle as following, which had been posted here last year but 
no reply. So I re-sent it here to get more suggestion.
    *****from  CERMM, Concordia University*********
    We are experiencing trouble getting our CASSCF calculations to
converge. We know that the selected active space is appropriate for
our system (O4), and we employ equilibrium geometry for different
isomers.
We usually get the calculation to converge more
easily by doing successive calculations with increasing active space
and required number of states.
In such calculation,

     casscf;
       occ,   20;
       closed,8;
       wf,32,1,0; state,1;
       wf,32,1,2; state,7;
       wf,32,1,4; state,1;

     casscf;
       occ,   20;
       core,8;
       wf,32,1,0; state,5;
       wf,32,1,2; state,7;
       wf,32,1,4; state,1;

     casscf;
       maxiter,100;
       IPRINT,MICRO;            !print information on convergence problem
       occ,   20;
       closed,8;

the first casscf achieves convergence, but the next two won't, and
both print something like:

 ITER. MIC  NCI  NEG     ENERGY(VAR)     ENERGY(PROJ)   ENERGY CHANGE     
GRAD(0)  GRAD(ORB)   GRAD(CI)    
STEP       TIME

   1   40   29    0    -299.46358617    -299.46327296    0.00031321    
0.00605306 0.00000000 0.28965734 
0.16D-01    907.83

 *** NO CONVERGENCE ***

 We faile to understand what parameters will determine the number of
 interations, and what should be done to increse the number of
 iterations. As you can see on below the second calculation, the maxiter 
did
 not seem to help.

 CONVERGENCE PROGRESS
 --------------------

 ITERATION                 GRADIENT                       ENERGY                    
CHANGE IN ENERGY
     0                    0.001496911537             -299.463625399664
     1                    0.000000000000             -299.436598472010                
0.027026927654   
(PROJECTED)

 ** WVFN ****  START OF ITERATION  2 AT TIME  1240.64

 Entering internal optimisation at time  1254.91

            Time  Orb. grad.   Orb. change   CI gradient    CI change         
Energy      DIIS

   1  1  1298.62  0.00007543    0.00000000    1.00000000    0.00000000   
-299.43659916     0
   1  2  1298.66  0.00000020    0.00060323    1.00000000    0.00000000   
-299.43659919     0
 Too many P-space configurations for threshold      0.40
 Actual threshold used:      0.34
   2  1  1406.90  0.00086038    0.00060336    0.00391926    0.00084597   
-299.40184999     0
   2  2  1406.94  0.00000689    0.00873427    0.00391926    0.00000000   
-299.40185670     0
   3  1  1406.99  0.00000001    0.00874368    0.00000000    0.00000000   
-299.40185670     0

 CPU FOR TRAOPM:      1.33 ( 490.7 MFLOPS)  NOP_CORE=  40 NOP_FILE=   0 
MINWK=   612061  NWUSED=  6372060 
NWL= 89999500

 End of internal optimisation at time    1408.33

 ENERGY =   -299.4018566977  ENERGY CHANGE =   0.0617687020     
0.0347417743 (PROJECTED)  GRADIENT =  
0.72E-02  STEP =   0.25E-01  CONV = 1

 *** NO CONVERGENCE ***


Any insights appreciated.




More information about the Molpro-user mailing list