[molpro-user] the error in the do loop

qyf_jh qyf_jh at 126.com
Fri Apr 23 03:12:59 BST 2010


dear molpro users:
I try to use the "do loop" to get the potential energy surface, but it puzzled me.
when I set the i=1,120, I will get the right outcome. But I set the i=1,150.it will be end with the error:
 ? Error
 ? MXDEL
 ? The problem occurs in AUGI
 
Maybe, the C-P bond has broken when the i=130, C-P bond length=2.6, and the molecular symmtry has been changed!
 
Who can do me a faver? Please!
PS:here is my input file and output file!
ang
C1
C2,C1,1.283
P,C2,r,C1,180
end
r=1.30
hf
do i=1,150
casscf
ecas(i)=energy
rasc(i)=r
emrci(i)=energy
emrda(i)=energd
r=r+0.01
end do
 
The error:
 ?WARNING, ROTATION 2.2 -  3.2  D1E= 0.136D-08  D2E= 0.191D-08 ELIMINATED DUE TO SMALL DERIVATIVES
   1   35    7    0   -2309.72914750   -2309.72916903   -0.00002154    0.00733804 0.00000319 0.00000000  0.45D-02      0.27
 ? Error
 ? MXDEL
 ? The problem occurs in AUGI
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20100423/f16f7302/attachment.html>


More information about the Molpro-user mailing list