[Pw_forum] relax question

Stefano de Gironcoli degironc at sissa.it
Fri Jul 13 17:24:54 CEST 2007


a tipical bond stiffness, if I recall correctly, is in the range of 0.1-1.0
Ry/bohr^2 thus a threshold of 1.d-6 Ry/bohr in the force means that you are
asking a numerical accuracy in your interatomic distances of the order of
1.d-5-1.d-6 bohr... are you sure you really need it ?
best regards,
stefano deGironcoli, SISSA & DEMOCRITOS

Quoting xu yuehua <njuxuyuehua at gmail.com>:

> hello all
> i am doing relax.
> in my inputfile,
> i use  etot_conv_thr=1.D-6
>   forc_conv_thr=1.D-6
> and the output is
>    atom   1 type  1   force =     0.00000290    0.00000095   -0.00000364
>      atom   2 type  1   force =     0.00000081   -0.00000108    0.00000033
>      atom   3 type  1   force =    -0.00000435   -0.00000425   -0.00000160
>      atom   4 type  1   force =    -0.00000312   -0.00000321   -0.00000080
>      atom   5 type  1   force =    -0.00000446   -0.00000036   -0.00000273
>      atom   6 type  1   force =    -0.00000048    0.00000230    0.00000112
>      atom   7 type  1   force =     0.00000164    0.00000118    0.00000060
>      atom   8 type  1   force =     0.00000131   -0.00000040    0.00000145
>      atom   9 type  2   force =    -0.00000362    0.00000154   -0.00000201
>      atom  10 type  2   force =    -0.00000029    0.00000191    0.00000302
>      atom  11 type  2   force =     0.00000576   -0.00000134    0.00000194
>      atom  12 type  2   force =     0.00000390    0.00000276    0.00000233
> 
> Total force =     0.000015     Total SCF correction =     0.000054
>      SCF correction compared to forces is too large, reduce conv_thr
> 
>      The maximum number of steps has been reached.
> my qustion is ,force at every atom is ok ,but why  the program asks to
> reduce conv_thr
> my conv_thr is not enough small ?
> or my total relax step (50) is not adequate?
> is there any principle about these parameter?
> need your help
> thanks a lot
> -- 
> Xu Yuehua
> physics Department of Nanjing university
> China
> 




----------------------------------------------------------------
  SISSA Webmail https://webmail.sissa.it/
  Powered by Horde http://www.horde.org/




More information about the users mailing list