TURBOMOLE Users Forum

TURBOMOLE Modules => Escf and Egrad => Topic started by: jc3 on October 19, 2008, 06:54:17 PM

Title: problems with S1 optimization
Post by: jc3 on October 19, 2008, 06:54:17 PM
Hi everyone,

I'm trying to optimize the first singlet excited state of a molecule with jobex -ex, but in the third cycle the calculation stop. The coordinates have changed to NaN, this is the end of the job.last:

          ------------------------------------------------
           cartesian gradient of the energy (hartree/bohr)
          ------------------------------------------------

  ATOM      1 c           2 c           3 c           4 c           5 s
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM      6 c           7 c           8 s           9 c          10 c
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM     11 s          12 c          13 c          14 s          15 c
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM     16 c          17 s          18 c          19 c          20 s
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM     21 c          22 c          23 s          24 s          25 c
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM     26 c          27 s          28 c          29 c          30 h
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM     31 h          32 s          33 c          34 c          35 s
 dE/dx  NaN           NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN           NaN

  ATOM     36 c          37 c          38 h          39 h
 dE/dx  NaN           NaN           NaN           NaN
 dE/dy  NaN           NaN           NaN           NaN
 dE/dz  NaN           NaN           NaN           NaN

 resulting FORCE  (fx,fy,fz) = ( NaN     , NaN     , NaN     )
resulting MOMENT (mx,my,mz) = ( NaN     , NaN     , NaN     )


 exx =  NaN           eyy =  NaN           ezz =  NaN
 eyz =  NaN           exz =  NaN           exy =  NaN


 **********************************************************************
  |maximum component of gradient| : NaN             (atom 39 h )
   gradient norm                  : NaN
 **********************************************************************
   ***  cartesian gradients written onto <gradient>  ***



 ==============================================================================
                           electrostatic moments
 ==============================================================================

              nuc           gs,e       ->  gs,tot          excit     ->  total
 ------------------------------------------------------------------------------
                          charge
 ------------------------------------------------------------------------------
         324.000000 NaN            NaN            NaN            NaN

 ------------------------------------------------------------------------------
                       dipole moment
 ------------------------------------------------------------------------------
   x     -27.030734 NaN            NaN            NaN            NaN
   y     -10.527207 NaN            NaN            NaN            NaN
   z      12.098837 NaN            NaN            NaN            NaN

   | dipole moment | =  NaN       a.u. =  NaN       debye

 ------------------------------------------------------------------------------
                     quadrupole moment
 ------------------------------------------------------------------------------
  xx    7265.998321 NaN            NaN            NaN            NaN
  yy    6525.527799 NaN            NaN            NaN            NaN
  zz    4449.608341 NaN            NaN            NaN            NaN
  xy    -264.687950 NaN            NaN            NaN            NaN
  xz    -941.560941 NaN            NaN            NaN            NaN
  yz   -1486.023271 NaN            NaN            NaN            NaN

     1/3  trace= NaN
    anisotropy= NaN

 ==============================================================================

                                                                                                                                                               2905,1        99%

I don't know that could be. I hope someone could help me.

Thank you very much.


Title: Re: problems with S1 optimization
Post by: Arnim on October 23, 2008, 01:03:43 PM
Hi,

from this part of the output it is not really possible to see, what went wrong.
You can run jobex with the -keep option in order to check, what the
underlying problem is.

Cheers,
Arnim