[Pw_forum] different forces between last 2 steps

[Jaret Qi]

is this Okay or I should take an action?

 


     On Tuesday, November 3, 2015 4:26 AM, Nicola Marzari <nicola.marzari at epfl.ch> wrote:
   

 On 03/11/2015 10:13, Jaret Qi wrote:
> Dear QE users,
> Why total forces are different in the last relaxation step and the final
> step that pw.x performs once the criteria satisfies, does this mean
> relaxation not good?
> take a look :

Hi - the code tells you exactly what has happened:

>      A final scf calculation at the relaxed structure.
>      The G-vectors are recalculated for the final unit cell
>      Results may differ from those at the preceding step.


This happens because different unit cells at the same cutoff have different
seta of corresponding plane waves - for an explanation you could check
G Francis and MC Payne circa 1991 (JPCondMatt) or G-M Rignanese
et al PRB circa 1996.

                nicola

-- 
----------------------------------------------------------------------
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL
Director, National Centre for Competence in Research NCCR MARVEL, EPFL
http://theossrv1.epfl.ch/Main/Contact http://nccr-marvel.ch/en/project


  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20151103/68cc20c6/attachment.html>