<DIV>dear molpro users:</DIV>
<DIV>I try to use the "do loop" to get the potential energy surface, but it puzzled me.</DIV>
<DIV>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:</DIV>
<DIV> ? Error<BR> ? MXDEL<BR> ? The problem occurs in AUGI</DIV>
<DIV> </DIV>
<DIV>Maybe, the C-P bond has broken when the i=130, C-P bond length=2.6, and the molecular symmtry has been changed!</DIV>
<DIV> </DIV>
<DIV>Who can do me a faver? Please!</DIV>
<DIV>PS:here is my input file and output file!</DIV>
<DIV>ang<BR>C1<BR>C2,C1,1.283<BR>P,C2,r,C1,180<BR>end</DIV>
<DIV>r=1.30<BR>hf</DIV>
<DIV>do i=1,150</DIV>
<DIV>casscf<BR>ecas(i)=energy<BR>rasc(i)=r<BR>emrci(i)=energy<BR>emrda(i)=energd</DIV>
<DIV>r=r+0.01<BR>end do</DIV>
<DIV> </DIV>
<DIV>The error:</DIV>
<DIV> ?WARNING, ROTATION 2.2 -  3.2  D1E= 0.136D-08  D2E= 0.191D-08 ELIMINATED DUE TO SMALL DERIVATIVES<BR>   1   35    7    0   -2309.72914750   -2309.72916903   -0.00002154    0.00733804 0.00000319 0.00000000  0.45D-02      0.27<BR> ? Error<BR> ? MXDEL<BR> ? The problem occurs in AUGI</DIV><br><br><span title="neteasefooter"/></span>