[Pw_forum] dE0s is positive which should never happen

Guntram Schmidt guntram.schmidt at chemie.uni-halle.de
Fri Jul 8 00:38:20 CEST 2011


Dear All,

I tried to optimize a crystal structure of a molecular compound and get 
this strange error. Is there something I've forgotten in the input file?
I've chosen a small ecutwfc to make some "preoptimization" - is it okay 
to do so (I'm still not able to converge the cutoffs - see my other post)?

&control
    calculation = 'relax'
    restart_mode='from_scratch',
    prefix='i2t0480',
    tstress = .true.
    tprnfor = .true.
    pseudo_dir = '/gpfs/home/aasfu/espresso-4.3.1/pseudo/',
    outdir='/gpfs/home/aasfu/tmp/'
    verbosity = 'high'
/
&system
  ibrav = 14,
  A = 8.2555,
  B = 17.0986,
  C = 16.3082,
  cosAB = 0,
  cosAC = -.4968073528,
  cosBC = 0,
  nat = 196,
  ntyp = 5,
  occupations = 'fixed',
  ecutwfc = 20,
  ecutrho = 80,
  input_dft = 'vdW-DF'
/
&electrons
     conv_thr        = 1.0d-5
     electron_maxstep= 1000
/
&ions
/

ATOMIC_SPECIES
  P 39.974 P.rpb-nc.UPF
  O 15.999 O.rpb-nc.UPF
  N 14.007 N.rpb-nc.UPF
  C 12.011 C.rpb-nc.UPF
  H 1.0079 H.rpb-nc.UPF


ATOMIC_POSITIONS angstrom
  H      -4.163902    3.036540    7.286078 1 1 1
  C      -3.730208    5.027159    7.959787 1 1 1
  H      -2.669409    4.951703    7.725352 1 1 1
  C      -4.231111    6.183196    8.511764 1 1 1
  H      -3.568842    7.024960    8.708593 1 1 1
and 210 more atoms...



The output says after 41 cycles of geometry optimizations in the 
subroutine stress:



    0:     number of scf cycles    =  30
    0:     number of bfgs steps    =  28
    0:
    0:     energy   old            =   -1592.2501932245 Ry
    0:     energy   new            =   -1547.5620256423 Ry
    0:
    0:     CASE: energy  _new > energy  _old
    0:
    0:
    0: 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    0:     from bfgs : error #         1
    0:     dE0s is positive which should never happen
    0: 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    0:
    0:     stopping ...

The strange thing is, that the energy was not -1592... before at all:
grep '!' i2t0480.xyz.relax.20.rpb-nc.UPF.out
    0:!    total energy              =   -1542.41391161 Ry
    0:!    total energy              =   -1544.87249615 Ry
    0:!    total energy              =   -1546.30689086 Ry
    0:!    total energy              =   -1546.89859418 Ry
    0:!    total energy              =   -1547.07908001 Ry
    0:!    total energy              =   -1547.12438175 Ry
    0:!    total energy              =   -1547.16304879 Ry
    0:!    total energy              =   -1547.21850898 Ry
    0:!    total energy              =   -1547.26799568 Ry
    0:!    total energy              =   -1547.31147437 Ry
    0:!    total energy              =   -1547.33202922 Ry
    0:!    total energy              =   -1547.35963149 Ry
    0:!    total energy              =   -1547.36650754 Ry
    0:!    total energy              =   -1547.38045447 Ry
    0:!    total energy              =   -1547.38597654 Ry
    0:!    total energy              =   -1547.39351928 Ry
    0:!    total energy              =   -1547.40184806 Ry
    0:!    total energy              =   -1547.41383843 Ry
    0:!    total energy              =   -1547.43030750 Ry
    0:!    total energy              =   -1547.45110560 Ry
    0:!    total energy              =   -1547.47294386 Ry
    0:!    total energy              =   -1547.48424959 Ry
    0:!    total energy              =   -1547.48774725 Ry
    0:!    total energy              =   -1547.49458349 Ry
    0:!    total energy              =   -1547.50098365 Ry
    0:!    total energy              =   -1547.50487139 Ry
    0:!    total energy              =   -1547.50802242 Ry
    0:!    total energy              =   -1547.51133404 Ry
    0:!    total energy              =   -1547.51564694 Ry
    0:!    total energy              =   -1547.52054790 Ry
    0:!    total energy              =   -1547.52614229 Ry
    0:!    total energy              =   -1547.53597920 Ry
    0:!    total energy              =   -1547.54614655 Ry
    0:!    total energy              =   -1547.55021920 Ry
    0:!    total energy              =   -1547.55191818 Ry
    0:!    total energy              =   -1547.55360256 Ry
    0:!    total energy              =   -1547.55527582 Ry
    0:!    total energy              =   -1547.55697509 Ry
    0:!    total energy              =   -1547.55871207 Ry
    0:!    total energy              =   -1547.56042619 Ry
    0:!    total energy              =   -1547.56202564 Ry


Thanks a lot for help,
Guntram



More information about the users mailing list