[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