[Pw_forum] Does total force include dispersion force?

2011-11-24 Thread WF
Hello everyone,

 I am doing a vc-relax calculation with PBE+D ( with .london. = true
and PBE functional). In the output, an additional part called "Dispersion
forces acting on atoms" is displayed after "Total forces acting on atoms"
after each BFGS cycle. Do the values in total force part already include
dispersion forces ? I am not sure with what I read from the codes.

Thanks very much.

 



Wu F,

College of Chemistry and Molecular Engineering, Peking University.



 

-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/2024/4c6a96f5/attachment.htm
 


[Pw_forum] Does total force include dispersion force?

2011-11-24 Thread Gabriele Sclauzero
Dear Wu,

Il giorno 24/nov/2011, alle ore 12.44, WF ha scritto:

> Hello everyone,
>  I am doing a vc-relax calculation with PBE+D ( with .london. = true 
> and PBE functional). In the output, an additional part called ?Dispersion 
> forces acting on atoms? is displayed after ?Total forces acting on atoms? 
> after each BFGS cycle. Do the values in total force part already include 
> dispersion forces ? I am not sure with what I read from the codes.

What have you read from the codes?

If you read inside PW/forces.f90 you will find:

IF ( llondon ) force(ipol,na) = force(ipol,na) + force_disp(ipol,na)

The "dispersion force" is added to the total force (before printing), therefore 
what you find in the output does actually include also the dispersion 
contribution that is printed separately. 

Indeed I think it may be a bit misleading to print by default that 
contribution, while for all other cases (e.g. Hubbard U and similar) this is 
done only with verbosity="high".


> Thanks very much.

You're welcome,


GS

>  
> 
> Wu F,
> College of Chemistry and Molecular Engineering, Peking University.
> 
>  
> ___
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum


? Gabriele Sclauzero, EPFL SB ITP CSEA
   PH H2 462, Station 3, CH-1015 Lausanne







-- next part --
An HTML attachment was scrubbed...
URL: 
http://www.democritos.it/pipermail/pw_forum/attachments/2024/4f972cfe/attachment.htm