[Wien] force output problem

Hua Peng penghua85 at gmail.com
Fri Jan 17 02:49:34 CET 2014


Dear Wien2k users,
I have questions about the Wien2K force output.
After I run run_lapw -fc 0.01 -p &, the calculation converged, which
means the force in the system should small than 0.01mRy/a.u..
But in end of  the scf file, the last interaction gives,

       TOTAL FORCE IN mRy/a.u. = |F|     Fx             Fy
Fz     with/without FOR in case.in2
:FOR001:   1.ATOM       0.000000       0.000000       0.000000
0.000000 total forces
:FOR002:   2.ATOM       0.000000       0.000000       0.000000
0.000000 total forces
:FOR003:   3.ATOM       6.480000       0.000000       0.000000
-6.480000 total forces

       TOTAL FORCE WITH RESPECT TO THE GLOBAL COORDINATE SYSTEM:
:FGL001:   1.ATOM                 0.000000000     0.000000000
0.000000000 total forces
:FGL002:   2.ATOM                 0.000000000     0.000000000
0.000000000 total forces
:FGL003:   3.ATOM                -6.480000000     0.000000000
0.000000000 total forces
:FSUM  : Sum of forces Fx,Fy,Fz            0.00000            0.00000
          0.00000

Besides, in the previous two  interaction, the force is 1.905.

     TOTAL FORCE IN mRy/a.u. = |F|     Fx             Fy
Fz     with/without FOR in case.in2
:FOR001:   1.ATOM          0.000          0.000          0.000
 0.000 partial forces
:FOR002:   2.ATOM          0.000          0.000          0.000
 0.000 partial forces
:FOR003:   3.ATOM          1.905          0.000          0.000
 1.905 partial forces

Is the TOTAL force in the last interaction the sum of force in the
cell? The equivalent atom for ATOM3 is 4.
Best,
Hua Peng


More information about the Wien mailing list