problem in ridft

[jitmukherjeechem]

Hi, I am facing some problems in the ridft program when running energy optimization of a functionalised carbon nanotube:

here is the last portiona of the job.1 file:


 data group $actual step is not empty
 due to the abend of ridft


 check reason for abend ...

 use the command  'actual -r'  to get rid of that


 MODTRACE: no modules on stack

  CONTRL dead = actual step
 rdgrad ended abnormally
 rdgrad ended abnormally
program stopped.
error in gradient step (1)


here is the statistics file:::


Wed May 30 00:17:03 IST 2012

 DSCF STATISTICS OUTPUT :

 integral neglect threshold       :  0.11E-10
 integral storage threshold THIZE :  0.10E-04
 integral storage threshold THIME :         5


 --------------------    1.ITERATION  --------------------
 norm of density matrix                :    33.3491934459
 total energy                          : -7618.0142193799
 one-electron energy                   : -101906.32232732
 two-electron energy                   :  48274.347573755
 damping factor ttr                    :           0.7000

 --------------------    2.ITERATION  --------------------

  weight factors dlin for linear combination of density matrices :
   1.0082920
 accerr=   2.017  crierr=1000.000  ilam=  0  grange= 0.0      chatol=   0.630

 norm of difference density matrix     :    5.78564850186
 norm of simple difference density     :    5.79225325207
 norm of density matrix                :    34.1198335858
 total energy                          : -7632.0901594518
 one-electron energy                   : -102075.01312581
 two-electron energy                   :  48428.962432181
 energy increment                      :           -14.08
 damping factor ttr                    :           0.6500

[jitmukherjeechem]

here is the last portion job.last file:

 Starting SCF iterations

          Overall gridpoints after grid construction =        359816

 ITERATION  ENERGY          1e-ENERGY        2e-ENERGY     NORM[dD(SAO)]  TOL
   1  -7618.0142193799    -101906.32233     48274.347574    0.000D+00 0.108D-10
                            Exc = -1050.26728708     Coul =  49324.6148608
                         Ntotal = 1207.0079601      Nspin = 1.0000144453
                             Na = 604.00398725         Nb = 603.00397280
                            current damping = 0.700

          max. resid. norm for Fia-block=  3.132D-01 for orbital    475a    alpha
          max. resid. fock norm         =  9.751D-01 for orbital   1727a    alpha

 ITERATION  ENERGY          1e-ENERGY        2e-ENERGY     NORM[dD(SAO)]  TOL
   2  -7632.0901594518    -102075.01313     48428.962432    0.579D+01 0.681D-11
                            Exc = -1057.47623741     Coul =  49486.4386696
                         Ntotal = 1207.0058170      Nspin =0.99994623233
                             Na = 604.00288164         Nb = 603.00293541
                            current damping = 0.650


 MODTRACE: no modules on stack

 opening file <errmat>
 ridft ended abnormally
 ridft ended abnormally
program stopped.
 ridft step ended abnormally
next step = ridft

[jitmukherjeechem]

Please give me suggestion. Here is my control file:::

$title
$operating system unix
$symmetry c1
$redundant    file=coord
$coord    file=coord
$user-defined bonds    file=coord
$atoms
c  1-144,168,170-171,174,176,181,184,190,192,194-195,198,200,205,208,214,216,  \
   218-219,222,224,229,232,238                                                 \
   basis =c def2-SVP                                                           \
   jbas  =c def2-SVP
h  145-166,177-180,182-183,185-186,188,201-204,206-207,209-210,212,225-228,    \
   230-231,233-234,236                                                         \
   basis =h def2-SVP                                                           \
   jbas  =h def2-SVP
n  167,169,172-173,175,187,191,193,196-197,199,211,215,217,220-221,223,235     \
   basis =n def2-SVP                                                           \
   jbas  =n def2-SVP
o  189,213,237                                                                 \
   basis =o def2-SVP                                                           \
   jbas  =o def2-SVP
$basis    file=basis
$rundimensions
   dim(fock,dens)=4748856
   natoms=238
   nshell=1281
   nbf(CAO)=3080
   nbf(AO)=2891
 dim(trafo[SAO<-->AO/CAO])=3458
   rhfshells=2
$uhfmo_alpha   file=alpha
$uhfmo_beta   file=beta
$uhf
$alpha shells
 a       1-604                                  ( 1 )
$beta shells
 a       1-603                                  ( 1 )
$scfiterlimit      200
$thize     0.10000000E-04
$thime        5
$scfdump
$scfintunit
 unit=30       size=0        file=twoint
$scfdiis
$disp
$drvopt
   cartesian  on
   basis      off
   global     off
   hessian    on
   dipole     on
   nuclear polarizability
$interconversion  off
   qconv=1.d-7
   maxiter=25
$optimize
 internal   on
   redundant  on
   cartesian  off
   global     off
   basis      off   logarithm
$coordinateupdate
   dqmax=0.3
   interpolate  on
   statistics    5
$forceupdate
   ahlrichs numgeo=0  mingeo=3 maxgeo=4 modus=<g|dq> dynamic fail=0.3
   threig=0.005  reseig=0.005  thrbig=3.0  scale=1.00  damping=0.0
$forceinit on
   diag=default
$energy    file=energy
$grad    file=gradient
$forceapprox    file=forceapprox
$lock off
$statpt
   itrvec      0
$dft
   functional b-p
   gridsize   m3
$scfconv   6
$scfdamp   start=0.700  step=0.050  min=0.050
$scforbitalshift  closedshell=.05
$ricore     2000
$rij
$jbas    file=auxbasis
$marij
$last step      ridft
$restart dscf
$end