[molpro-user] Convergence Problem in MCSCF Calculations

Sunil sunil.phys at gmail.com
Mon Jan 26 05:46:37 GMT 2009


Dear Sir,

I am facing a convergence problem when using Molpro to calculate potential
energy surfaces of CO2 molecule. The input of the program is as follows:


*** CO2 PE curve                        ! Title: PE Curve for CO2
memory,32,m
maxiter=100
geomtyp=zmat
geometry
ang                                     ! ang is used if we intend to
provide distance in Angstrom; default for Z-matrix is au
C
O1 C rco1(i)
O2 C rco2(i) O1 oco
end

r1(1) = 0.400000 Ang
r2(1) = 0.400001 Ang
dr = 0.1 Ang
do i = 1,19
r1(i+1) = r1(i)
r2(i+1) = r2(i) + dr
enddo

oco = 180 Deg

basis = vqz                             ! Basis set specification

do i = 1,20

rco1(i) = r1(i)
rco2(i) = r2(i)

hf;
enhf(i)=energy

mp2;
enmp2(i)=energy

multi;
emcscf(i) = energy

enddo

table rco1,rco2,enhf,enmp2,emcscf               !emscf          !,ecas,emrci
head R1,R2,SCF,MP2,MCSCF
save co2_mcscf.tab
Title Results for CO2 (Angle = 180 Deg)
!sort 3,1,2
---

It produces the output but without convergence for the MCSCF calculations.

Thank you in advance for any kind of help.

regards
Sunil
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20090126/cb1a6f52/attachment.html>


More information about the Molpro-user mailing list