[molpro-user] rs3 opt

Mintz, Benjamin J. bm0037 at gmail.com
Mon Dec 20 20:14:45 GMT 2010


Dear Users,

I am currently trying to perform a CASPT3 optimization, and I am running into a problem.  Specifically, I am trying to optimize the first excited state for my molecule.  The optimization converges, but the output looks like the optimization is utilizing the RS2 energies and not the RS3 energies.  I pasted part of the output below, which lists the RSPT2 and RSPT3 energies and the energies from the optimization.  I also pasted part of my input.  Why is the optimization not optimizing the RS3 energies?  Any help would be appreciate.  Let me know if I need to provide any more information. 

 {hf;wf,35,4,1;orbprint,4}
 {multi;occ,9,3,7,2;frozen,3,0,2,0;closed,8,0,6,0;wf,35,4,1;orbprint,4
 select
 con,2,2,2,0,0,1,0
 con,0,2,2,0,2,1,0
 con,2,2,0,0,0,2,1
 con,0,2,0,0,2,2,1
 }
 {rs3;state,1,2}
 {optg,numerical;
 inactive,d1;
 inactive,d2;
 inactive,d3;
 inactive,d4}


 

 Reference energy                    -192.00677670
 Nuclear energy                       130.22929885
 Kinetic energy                       192.49019540
 One electron energy                 -515.10328409
 Two electron energy                  192.28496840
 Virial quotient                       -1.00051338
 Correlation energy                    -0.58224014
 !RSPT2 STATE 2.4 Energy             -192.589016841281

 !RSPT2 STATE 2.4 Dipole moment         0.00000000     0.00000000     0.56683718
 Dipole moment /Debye                   0.00000000     0.00000000     1.44066205

 !RSPT expec            <2.4|H|2.4>  -192.490833418072

 Correlation energy                    -0.60081552
 !RSPT3 STATE 1.4 Energy             -192.607592225598

.
.
.

 END OF GEOMETRY OPTIMIZATION.    TOTAL CPU:     53836.6 SEC

 ITER.   ENERGY(OLD)    ENERGY(NEW)      DE          GRADMAX     GRADNORM    GRADRMS     STEPMAX     STEPLEN     STEPRMS
   1  -192.58901684  -192.59615035    -0.00713351  0.12323135  0.12793186  0.01421465  0.08500958  0.09811296  0.01090144
   2  -192.59615035  -192.59670839    -0.00055804  0.02893302  0.03046462  0.00338496  0.02788328  0.03482585  0.00386954
   3  -192.59670839  -192.59673641    -0.00002801  0.00508719  0.00606747  0.00067416  0.00630698  0.00935194  0.00103910
   4  -192.59673641  -192.59673813    -0.00000173  0.00063479  0.00121156  0.00013462  0.00214935  0.00300464  0.00033385
   5  -192.59673813  -192.59673824    -0.00000011  0.00021685  0.00036503  0.00004056  0.00049652  0.00066890  0.00007432



Thanks,
Benjamin J. Mintz
Research Associate, Postdoctoral
Oak Ridge Leadership Computing Facility (OLCF)
Oak Ridge National Laboratory
Phone: (865) 574-6639
Email: mintzbj at ornl.gov
Gmail: bm0037 at gmail.com




More information about the Molpro-user mailing list