Linear dependence in geometry optimization
Dave Moore
dtmoore at rijnh.nl
Fri Feb 7 14:03:13 GMT 2003
Dear all,
I am having a problem with a fairly simple geometry optimization (mp2/VDZ
for dimethyl ether) in MOLPRO 2002.3 with a recently updated patch-set (I
believe it to be completely current). No matter what symmetry I specify
(C2V, Cs, or C1), this optimization always crashed with an error message
like:
?ERROR IN ORTHS: LINEAR DEPENDENCY DETECTED FOR VECTOR 37.1
OVERLAP= -0.4437E-15
after 1 or 2 steps. I have gone through the input and output several times
looking for evidence of a careless mistake on my part, and I can't seem to
find it. Is this a problem with MOLPRO or am I just overlooking something
silly? I would appreciate any insights ... the input file is appended to
the end of this email. Thanks in advance for any help.
Dave Moore
FOM Rijnhuizen
Edisonbaan 14
3439 MN, Nieuwegein
The Netherlands
+31+30 609-6826
FAX: 603-1204
****** BEGIN MOLPRO INPUT ******
***
mem,60,M
basis,vdz
gthresh, energy=1.0e-10, orbital=1.0e-10, optgrad=4.0e-5
gexpec, dm, sm, qm
file,2,methyl_ether_monomer.C2v_mp2.vdz.wf
geometry = {ANGSTROM;
O1
X1 O1 1.0
X2 O1 1.0 X1 90.0
C1 O1 rco X1 acox X2 90.0
C2 O1 rco X1 acox X2 270.0
H1 C1 rch1 O1 ahco1 C2 180.0
H2 C2 rch1 O1 ahco1 C1 180.0
H3 C1 rch2 O1 ahco2 C2 60.0
H4 C1 rch2 O1 ahco2 C2 300.0
H5 C2 rch2 O1 ahco2 C1 300.0
H6 C2 rch2 O1 ahco2 C1 60.0
}
rco = 1.41004
rch1 = 1.09981
rch2 = 1.10938
acox = 55.285
ahco1 = 107.534
ahco2 = 111.783
hf;orbitals,2101.2
mp2; dm,8101.2; natorb,,,core=1
optg,energy=1.0e-8,grad=1.0e-5;
put,molden,methyl_ether_monomer.C2v.mp2.vdz.molden.dat;orb,8101.2
---
****** END MOLPRO INPUT ******
More information about the Molpro-user
mailing list