<div>Dear Users,</div><div><br></div><div>I have been trying to calculate the potential energy surface of three-times-ionized CO2 molecule using the program given below. The program does not converge even in the HF level when the variable "i" is above 10 or so. Can you please help me in deciding what modification I can do to the program so as to get convergence? Thanks in advance.</div>
<div><br></div><div>-sunil</div><div>---------------------------------------------------------------------------------------</div><div><br></div><div>*** CO2 PE curve                        ! Title: PE Curve for CO2</div>
<div>print basis,orbitals,civector<br></div><div>geomtyp=zmat</div><div>geometry</div><div>ang                                  </div><div>C</div><div>O1 C rco1(k)</div><div>O2 C rco2(k) O1 oco</div><div>end</div><div><br>
</div><div>r1(1) = 0.7500000 Ang</div><div>r2(1) = 0.7500001 Ang</div><div>dr = 0.05 Ang</div><div><br></div><div>do i = 1,50</div><div>r1(i+1) = r1(i) + dr</div><div>r2(i+1) = r2(i) + dr</div><div>enddo</div><div><br></div>
<div>oco = 180 Degree</div><div><br></div><div>basis = 6-311G++(2DF,2PD)               ! Basis set specification</div><div>k = 0</div><div>do i = 1,#r1</div><div>do j = 1,i</div><div>k=k+1</div><div>rco1(k) = r1(i)</div><div>
rco2(k) = r2(j)</div><div><br></div><div>hf;</div><div>occ 7 2 1 0</div><div>wf 19 2 1</div><div>open 2.2</div><div>enhf(k)=energy</div><div><br></div><div>enddo</div><div>enddo</div><div><br></div><div>table rco1,rco2,enhf                   </div>
<div>head R1,R2,SCF                          </div><div>save co2.tab</div><div>Title Results for CO2 (Angle = 180 Deg)</div><div><br></div><div>---</div><div><br></div>