[Wien] Question regarding spin scaling (xc1) implementation in wien2k

[Wen Fong Goh]

Dear wien2k developers,


I would like to know does the wien2k calculation with spin scaling (xc1) give a consistent total energy value? I looked up the subroutines in SRC_lapw0, particularly in vxclm2.f and only found the scaled exchange potential:

            vhelp=vxu
            vxu=(vxu*(1.d0+xc1)+vxd*(1-d0-xc1))/2.d0
            vxd=(vhelp*(1.d0-xc1)+vxd*(1-d0+xc1))/2.d0
I wonder if it also gives the total energy. Also, I have problem finding where the total energy is calculated, can anyone point it out to me? Thanks.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20160707/d659850d/attachment.html>