<div dir="ltr">Dear users,<div>                I am computing potential energy surface of CCH+ using Molpro. I have observed that there is a difference in energy values for the same geometry when it was computed in a loop and when it is computed separately as a single point.</div><div><br></div><div>r(C-H) = 2.4 a.u , r(C-C) = 21.9 a.u, C-C-H angle = 179 degree : -75.74572157 hartree (Single point)</div><div>                                                                                             -75.78211095 hartree (in a loop)</div><div><br></div><div>My input for loop is</div><div><br></div><div><div> ***,CCH+</div><div>memory,70,m</div><div>gthresh,energy=0.32d-6</div><div><br></div><div>r=1.0d0</div><div>rc=2.4d0</div><div>theta = 179 deg</div><div><br></div><div>geometry={</div><div>C1;</div><div>C2,C1,r;</div><div>H,C2,rc,C1,theta}</div><div><br></div><div>basis=aug-cc-pVQZ</div><div><br></div><div>do i=1,56</div><div><br></div><div>{hf;wf,12,2,2}<br></div><div><br></div><div>{multi;</div><div>maxiter,40;</div><div>pspace,0.3;</div><div>wf,12,2,2}            </div><div><br></div><div>{ci;</div><div>option,maxit=90                                                                 </div><div>option,maxiti=9000;</div><div>!option,nstati=2;</div><div>option,cluster=2;                                                                 </div><div>wf,12,2,2}                                                               </div><div><br></div><div>e(i)=energy                                                                 </div><div>ed(i)=energd2</div><div><br></div><div>if(i.lt.3) then</div><div>r=r+0.2d0</div><div>elseif(i.ge.3).and.(i.lt.24) then</div><div>r=r+0.1d0</div><div>elseif(i.ge.24).and.(i.lt.36) then</div><div>r=r+0.2d0</div><div>elseif(i.ge.36).and.(i.lt.44) then</div><div>r=r+0.5d0</div><div>else</div><div>r=r+1.0d0</div><div>endif</div><div><br></div><div>{table,sr,capr,e,ed</div><div>save,179d.dat,new</div><div>}</div><div><br></div><div>enddo</div></div><div><br></div><div>Is there something I am doing wrong? Any help is greatly appreciated. Thanks in advance.</div><div><br></div><div>Bhargava,</div><div>Dept. of Chemistry,</div><div>IIT Madras.</div></div>