<div dir="auto"><div>How different are the total energies?</div><div><br></div><div data-smartmail="gmail_signature">___<br>Emeritus Professor Laurence Marks (Laurie)<br>Department of Materials Science and Engineering, Northwestern University<br><a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a><br>"Research is to see what everybody else has seen, and to think what nobody else has thought" Albert Szent-Györgyi</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Oct 19, 2024, 04:21 pluto via Wien <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.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">Dear Prof. Blaha, dear All,<br>
<br>
My sequence is:<br>
<br>
instgen_lapw -ask<br>
init_lapw -sp -b (-prec 1 -ecut 0.999)<br>
<br>
Converge without SOC runsp -I (-ec 0.0001 -cc 0.001), save_lapw<br>
<br>
init_so_lapw with M=001<br>
<br>
Converge calculation with SOC, runsp -p -NI -so (-ec 0.0001 -cc 0.001), <br>
save_lapw<br>
<br>
At this point one can do another init_so_lapw (symmetso it included in <br>
there) and set M=00-1, and rerun runsp -p -NI -so. However, this <br>
produces something that is not right, eigenvalues at Gamma are different <br>
compared to the M=001.<br>
<br>
What works is to re-start from the saved non-SOC calculation. Restore <br>
non-SOC calculation (maybe in a new folder), do init_so_lapw with <br>
M=00-1, and runsp -orb -p -NI -so. Then the eigenvalues at Gamma are the <br>
same as in the M001 calculation and things look reasonable.<br>
<br>
All this seems to imply that one should not change the direction of M on <br>
top of previous SOC calculation.<br>
<br>
Best,<br>
Lukasz<br>
<br>
<br>
<br>
<br>
<br>
On 2024-10-18 16:12, Peter Blaha wrote:<br>
> It is not clear what you actually did. For sure you should not just<br>
> modify case.inso.<br>
> <br>
> You should rerun symmetso (or init_so) when changing the directions<br>
> and accept the changes.<br>
> <br>
> Most likely different symmetry operations will survive for the 2 cases.<br>
> <br>
> <br>
> Am 18.10.2024 um 16:05 schrieb pluto via Wien:<br>
>> Dear All,<br>
>> <br>
>> I am calculating a slab with 44 atoms (22 inequivalent atoms). I have <br>
>> two equivelent manganese atoms in the slab.<br>
>> <br>
>> Space group discovered by sgroup is 164, it has inversion.<br>
>> <br>
>> I run "instgen_lapw -ask" and I set Mn spin to u, and all other atoms <br>
>> to "n" (other atoms are non-magnetic). I monitor this during the SCF <br>
>> cycle by "grep :MM case.scf", and I am getting around 4.33 for the Mn <br>
>> moment, a reasonable number.<br>
>> <br>
>> The calculation is spin-polarized with SOC. During init_so_lapw I <br>
>> accept changes made by symmetso.<br>
>> <br>
>> I run two separate calculations with magnetization 001 and with 00-1. <br>
>> This is set in case.inso. Energy (0.0001) and charge (0.001) is <br>
>> converged.<br>
>> <br>
>> The two calculations give different eigenvalues, also at the Gamma <br>
>> point. Differences are not tiny, they are of the order of 50 meV and <br>
>> bands look slightly different. Correct me if I am wrong, but I think <br>
>> this is not physical in this setup.<br>
>> <br>
>> Could you please suggest what can be wrong and what can be checked.<br>
>> <br>
>> Am I allowed to converge with M = 001, save_lapw, then re-run <br>
>> init_so_lapw with M=00-1, and then re-run the SCF cycle with -NI? Or <br>
>> should I always run M=00-1 from the scratch, that is on top of <br>
>> converged non-SOC calculation?<br>
>> <br>
>> Best,<br>
>> Lukasz<br>
>> <br>
>> <br>
>> <br>
>> <br>
>> <br>
>> _______________________________________________<br>
>> Wien mailing list<br>
>> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">Wien@zeus.theochem.tuwien.ac.at</a><br>
>> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer 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/" rel="noreferrer noreferrer" target="_blank">http://www.mail-archive.com/</a> <br>
>> <a href="http://wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer noreferrer" target="_blank">wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank" rel="noreferrer">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer 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 noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>