The 'x' script should take care of the selection complex/real version for you.<div><br></div><div>Oleg<br><br>On Friday, December 19, 2014, Zhu, Jianxin <<a href="mailto:jxzhu@lanl.gov">jxzhu@lanl.gov</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Elias,<br>
<br>
For the spin-orbit coupling (but non spin-poalrized) case, we usual have<br>
the following for the charge self-consistency run<br>
<br>
x lapw1<br>
x lapwso<br>
x lapw2 -so -c<br>
<br>
In the wien2wannier90 user guide, I see there is an option of "-c" for<br>
w2w. In the above situation, should we use this option<br>
<br>
x w2w -c<br>
(It is equivalent to run w2wc case.def)?<br>
I tend to do so (by thinking the wave function can be complex)  but I am<br>
not sure.<br>
<br>
<br>
In the case without inversion symmetry, since I have<br>
x lapw1 -c<br>
x lapwso -c<br>
x lapw2 -so -c<br>
for the charge self-consistency, I do use<br>
x w2w -c  (after running x lapw1 -c, and x lapwso -c)<br>
for initial projection Amn.<br>
<br>
<br>
Thank you in advance for the explanation.<br>
<br>
<br>
Jianxin<br>
<br>
<br>
--<br>
#################################<br>
Jian-Xin Zhu, Ph.D<br>
Theoretical Division, MS B262<br>
Los Alamos National Laboratory<br>
Los Alamos, New Mexico 87545<br>
Phone: (505) 667 2363 (T-4);<br>
(505) 667 6602 (CINT)<br>
Fax: (505) 665 4063<br>
Email (main): <a href="javascript:;" onclick="_e(event, 'cvml', 'jxzhu@lanl.gov')">jxzhu@lanl.gov</a><br>
Email (backup): <a href="javascript:;" onclick="_e(event, 'cvml', 'physjxzhu@gmail.com')">physjxzhu@gmail.com</a><br>
URL: <a href="http://theory.lanl.gov" target="_blank">http://theory.lanl.gov</a><br>
#################################<br>
<br>
<br>
<br>
<br>
<br>
On 12/12/14 4:02 AM, "Elias Assmann" <<a href="javascript:;" onclick="_e(event, 'cvml', 'elias.assmann@gmail.com')">elias.assmann@gmail.com</a>> wrote:<br>
<br>
>On 12/12/2014 07:34 AM, Zhu, Jianxin wrote:<br>
>> 1/ runsp_c_lapw –so  is one way for spin-orbit coupling but<br>
>> non-spin-polarized case.<br>
><br>
>runsp_c does a *spin-polarized* calculation where up and dn are<br>
>*constrained* to be equal.<br>
><br>
>> Elias has also suggested the following way after faking the spin<br>
>> polarized calculation --<br>
>> $ x w2w -so -up && x w2w -so -dn<br>
>> $ cp CASE.eig CASE.eigup<br>
>> $ cp CASE.eig CASE.eigdn<br>
>> $ x wannier90 -so<br>
><br>
>I think you are referring to this message<br>
><<a href="http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg11327.ht" target="_blank">http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg11327.ht</a><br>
>ml>.<br>
>  This is just due to a minor bug in the ‘wannier90’ wrapper script.<br>
>Anyway, it applies to non-SP SO cases, and so should not be relevant here.<br>
><br>
>> 2/ For the initial projection, the following order may be more<br>
>> consistent with the Bloch wavefunction in the spin space (spinor).<br>
><br>
>Let me note that the order of orbitals is completely up to your<br>
>preference (it will make no difference to the calculation).<br>
><br>
>> From: 李志 <<a href="javascript:;" onclick="_e(event, 'cvml', 'zhili@hfut.edu.cn')">zhili@hfut.edu.cn</a> <mailto:<a href="javascript:;" onclick="_e(event, 'cvml', 'zhili@hfut.edu.cn')">zhili@hfut.edu.cn</a>>><br>
>>     Firtsly, I do a spin polarization calculation with SOC, and zero<br>
>>     local magnetic moment, i.e. runsp_c_lapw -so<br>
>><br>
>>     Then, preparew2w_dir wannier<br>
>><br>
>>            init_w2w -up<br>
>><br>
>>     However, 26 inital projections are requied, and we just 7+5+1=13<br>
>>     orbiatls basis.<br>
><br>
>With SO, it is not possible to treat “up” and “dn” states separately.<br>
>You need to define 13*2=26 Wannier functions; Amn and Mmn will be<br>
>computed for the “up” and “dn” parts and then added, i.e.<br>
><br>
>Amn = Amn(up) + Amn(dn)<br>
>Mmn = Mmn(up) + Mmn(dn)<br>
><br>
>>     So, I should type into f,d,and s-orbitals two times,i.e.:<br>
>><br>
>>     1:f<br>
>><br>
>>     1:f<br>
>><br>
>>     1:d<br>
>><br>
>>     1:d<br>
>><br>
>>     1:s<br>
>><br>
>>     1:s<br>
>><br>
>>     correct?<br>
><br>
>That looks correct.  For brevity, I would write it as<br>
><br>
>1:s,p<br>
>1:s,p<br>
><br>
><br>
>       Elias<br>
>_______________________________________________<br>
>Wien mailing list<br>
><a href="javascript:;" onclick="_e(event, 'cvml', 'Wien@zeus.theochem.tuwien.ac.at')">Wien@zeus.theochem.tuwien.ac.at</a><br>
><a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
>SEARCH the MAILING-LIST at:<br>
><a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="javascript:;" onclick="_e(event, 'cvml', 'Wien@zeus.theochem.tuwien.ac.at')">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" 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" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br><br>-- <br><div dir="ltr">--<br>Oleg Rubel<br>Scientist, Thunder Bay Regional Research Institute<br>Adjunct Professor, Dept Physics, Lakehead University<br>290 Munro St, Thunder Bay, P7A 7T1, Ontario, Canada<br>Phone: +1-807-7663350<br>Fax: +1-807-3441948<br>E-mail: <a href="mailto:orubel@lakeheadu.ca" target="_blank">orubel@lakeheadu.ca</a><br>Homepage: <a href="http://www.tbrri.com/~orubel/" target="_blank">http://www.tbrri.com/~orubel/</a></div><br>