[molpro-user] Can't optimize C4H10 in Molpro 2009.1
Konstantin Tokarev
annulen at yandex.ru
Fri May 7 09:50:38 BST 2010
Geometry of C4H10 doesn't converge with any optimization algorithm. I've tried qsd (default), rf, diis with step interpolation (as recommended for floppy molecules) - nothing helps. I've even used pre-calculated hessian (see input below) - nothing helps. Gradient doesn't fall below 1.d-2
The same things occur with other molecules with many degrees of freedom. Level of theory doesn't matter, problem is in optimizer
basis, 6-31G**
r2 = 1.09419 ang
r3 = 1.09467 ang
a3 = 108.09176 ang
r4 = 1.09473 ang
a4 = 108.09063 ang
d4 = 117.02306 ang
r5 = 1.52031 ang
a5 = 110.25733 ang
d5 = 238.50846 ang
r6 = 1.09614 ang
a6 = 109.16630 ang
d6 = 301.55424 ang
r7 = 1.09631 ang
a7 = 109.16295 ang
d7 = 58.44923 ang
r8 = 1.52778 ang
a8 = 111.49085 ang
d8 = 180.02562 ang
r9 = 1.52064 ang
a9 = 111.48546 ang
d9 = 179.97438 ang
r10 = 1.09443 ang
a10 = 110.98794 ang
d10 = 299.73523 ang
r11 = 1.09416 ang
a11 = 110.24023 ang
d11 = 179.97438 ang
r12 = 1.09450 ang
a12 = 110.98796 ang
d12 = 60.25483 ang
r13 = 1.09610 ang
a13 = 109.86122 ang
d13 = 58.84686 ang
r14 = 1.09617 ang
a14 = 109.85801 ang
d14 = 301.13536 ang
symmetry,nosym
geometry={
H
C, 1, r2
H, 2, r3, 1, a3
H, 2, r4, 1, a4, 3, d4
C, 2, r5, 1, a5, 3, d5
H, 5, r6, 2, a6, 1, d6
H, 5, r7, 2, a7, 1, d7
C, 5, r8, 2, a8, 1, d8
C, 8, r9, 5, a9, 2, d9
H, 9, r10, 8, a10, 5, d10
H, 9, r11, 8, a11, 5, d11
H, 9, r12, 8, a12, 5, d12
H, 8, r13, 5, a13, 2, d13
H, 8, r14, 5, a14, 2, d14
}
{rks,b3lyp
wf,34,1,0}
{frequencies, forward
thermo,print,thermo}
{optg,readhess
method,diis,5,step
}
{frequencies
thermo,print,thermo}
P.S. In GAMESS C4H10 converges with default settings
--
Regards,
Konstantin
More information about the Molpro-user
mailing list