[Pw_forum] relax question

lan haiping lanhaiping at gmail.com
Fri Jul 13 15:51:24 CEST 2007


Dear yuehua

you should do a look at $QE/Doc/INPUT_PW, you can find all that you want to
know.

Regards,

On 7/13/07, xu yuehua <njuxuyuehua at gmail.com> wrote:
>
> 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<https://mail.google.com/mail?view=cm&tf=0&to=dingxunlei@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<https://mail.google.com/mail?view=cm&tf=0&to=Pw_forum@pwscf.org>
> > > http://www.democritos.it/mailman/listinfo/pw_forum
> > >
> >
> > = = = = = = = = = = = = = = = = = = = =
> >
> > _______________________________________________
> > Pw_forum mailing list
> > Pw_forum at pwscf.org<https://mail.google.com/mail?view=cm&tf=0&to=Pw_forum@pwscf.org>
> > http://www.democritos.it/mailman/listinfo/pw_forum
> >
>
>
>
> --
> Xu Yuehua
> physics Department of Nanjing university
> China
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org<https://mail.google.com/mail?view=cm&tf=0&to=Pw_forum@pwscf.org>
> http://www.democritos.it/mailman/listinfo/pw_forum
>
>


-- 
Hai-Ping Lan
Department of Electronics ,
Peking University , Bejing, 100871
lanhaiping at gmail.com, hplan at pku.edu.cn
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20070713/6a9c85c5/attachment.html>


More information about the users mailing list