<div><div dir="auto">Dear Dr Fabien Tran,</div></div><div dir="auto"><br></div><div dir="auto">Thanks for your reply. Very clear explanation as it is always.</div><div dir="auto"><br></div><div dir="auto">I got another question about the difference charge density calculation.</div><div dir="auto">For magnetic materials (including ferromagnetic, anti ferromagnetic or quadrupolar-antiferromagntic materials etc.), what I can think of to calculate the difference charge density is </div><div dir="auto"><br></div><div dir="auto">1) calculate the valence charge density (option “rho”) for spin-up electrons (x lapw5 -up)</div><div dir="auto"><br></div><div dir="auto">2) calculate the valence charge density (Option “rho” for spin-down electrons (x lapw5 -dn)</div><div dir="auto"><br></div><div dir="auto">3) sum 1 and 2 and subtract them by the superposed atomic valence densities (option “Over”.</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">Can I ask if this is correct way to do the calculation of the difference charge density for magnetic materials? and also if there is an automatic option to directly calculate it in Wien2k “, just like option “Diff” for non-magnetic materials ?</div><div dir="auto"><br></div><div dir="auto">Many thanks,</div><div dir="auto">Ding</div><div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 3 May 2020 at 12:16 am, Tran, Fabien <<a href="mailto:fabien.tran@tuwien.ac.at">fabien.tran@tuwien.ac.at</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">As I already said in my previous email, the last column is the occupancy.<br>
An occupancy of zero means that the corresponding band is not included<br>
in the valence electron density (stored case.clmval, that is used for plotting).<br>
"Energy to separate semicore and valence states" seems to be a prediction<br>
of what it means.<br>
<br>
<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of Ding Peng <<a href="mailto:ding.peng@monash.edu" target="_blank">ding.peng@monash.edu</a>><br>
Sent: Saturday, May 2, 2020 5:23 PM<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] Questions about difference electron density map<br>
<br>
Dear Dr Fabien Tran,<br>
<br>
<br>
Thanks for your reply. Actually, I set the Emin as -0.2 Ry, which means I discarded the states that are below -0.2 Ry. Could I ask you if the last number of the "band lines" is the occupancy? If it is 0.0000000D+00 (double precision of 0), does it mean that the corresponding band is included as the valence state? Also I noticed there is a line "Energy to separate semicore and valence states: -3.971792594507138E-002". Could I ask if the value "-3.971792594507138E-002" in this line the predicted Emin by Wien2k?<br>
<br>
<br>
Sorry for asking so many questions in this email.<br>
<br>
<br>
Best regards,<br>
Ding<br>
<br>
<br>
<br>
On 2 May 2020, at 10:12 pm, Tran, Fabien <<a href="mailto:fabien.tran@tuwien.ac.at" target="_blank">fabien.tran@tuwien.ac.at</a>> wrote:<br>
<br>
<br>
<br>
>From your output2, we can see that you discarded the states that are below -0.3 Ry<br>
<br>
<br>
<br>
(the occupancy in the last column is zero for them) to calculate the electron density.<br>
You have to check that the number of remaining electrons used to calculate the<br>
valence electron density is the same as the number of electrons obtained from lstart to<br>
create the valence electron density from superimposed atomic densities.<br>
<br>
<br>
<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of Ding Peng <<a href="mailto:ding.peng@monash.edu" target="_blank">ding.peng@monash.edu</a>><br>
Sent: Saturday, May 2, 2020 1:11 PM<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] Questions about difference electron density map<br>
<br>
<br>
<br>
<br>
<br>
Dear Prof Peter Blaha and other WIEN2k experts,<br>
<br>
<br>
<br>
Many thanks for your advice. Your explanation is very clear and I can understand most of them. But I still get some questions to ask:<br>
<br>
<br>
<br>
<br>
<br>
By following your advices, I've checked the band ranges in output2 and they are : <br>
<br>
<br>
<br>
Compared to the information in case.scf1, I was thinking the Emin should be set slightly above the energy band 10, which is from -0.3524 to -0.2916 (please correct me if I'm wrong). Does that mean the Emin of the whole valence electrons should be set slightly above Emax of energy band 10, which is -0.2916? <br>
<br>
<br>
Many thanks,<br>
Ding _______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div></div>