[Pw_forum] relax question

xu yuehua njuxuyuehua at gmail.com
Fri Jul 13 15:44:59 CEST 2007


thanks for above  Mr Ding and Mr  Lan answers.
and another related question is  whether conv_thr  is in &ELECTRONS? is it
about scf ?
i don not understand how a flow chart when relaxing ions and its details?
 hope your kind help
2007/7/13, Ding Xunlei <dingxunlei at gmail.com>:
>
> Dear Yuehua,
>
> The default number is :
> etot_conv_thr=1.D-4
> forc_conv_thr=1.D-3
>
> So, 1. Your setting for forc_conv_thr= 1.D-6 is too small.
>        2. If you do need to set force like that, you need to set
> etot_conv_thr even smaller, saying 1.D-7.  And also, you need to set TRUE
> smaller than etot_conv_thr, saying 1.D-8.
>
> Hope it help.
>
> Yours sincerely,
> Ding
> 2007-07-13
> ______________________________________________
>
> ======= 2007-07-13 16:27:54 You have written: =======
>
> >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
> >
> >_______________________________________________
> >Pw_forum mailing list
> >Pw_forum at pwscf.org
> >http://www.democritos.it/mailman/listinfo/pw_forum
> >
>
> = = = = = = = = = = = = = = = = = = = =
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>



-- 
Xu Yuehua
physics Department of Nanjing university
China
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20070713/72a067be/attachment.html>


More information about the users mailing list