# [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}

method,diis,5,step
}

{frequencies
thermo,print,thermo}

P.S. In GAMESS C4H10 converges with default settings

--
Regards,
Konstantin

```