Hi,
It seems that there is a mistake for scaling with LDA.
In vxclm2.f (in directory SRC_lapw0) there are these
three lines:
vxcu1=vxu
vxu=(vxcu1*(1.d0+XC1)+vxd*(1.d0-XC1))/2.d0
vxc=(vxcu1*(1.d0-XC1)+vxd*(1.d0+XC1))/2.d0
in the 3rd line it should be "vxd" instead of "vxc":
vxd=(vxcu1*(1.d0-XC1)+vxd*(1.d0+XC1))/2.d0
Probably there is the same error in the WIEN2k that you
are using. So, you need to replace vxc by vxd and recompile lapw0.
F. Tran
On Wed, 30 Mar 2016, Wen Fong Goh wrote:
Hi,
I am using the spin scaling option to reduce the spin polarization in an
antiferromagnetic calculation. I tried xc1 = 0 in both LDA and GGA, but only
GGA gives
zero magnetic moment. We know that with xc1 = 0, the magnetic part of the
exchange correlation functional vanishes (equation 11 in Ortenzi et al. 2012),
so the
system should become non-magnetic. Do you know what is going on with the spin
scaling in LDA?
_______________________________________________
Wien mailing list
Wien@zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
SEARCH the MAILING-LIST at:
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
