<div dir="ltr"><div><div><div><div>Hi Oleg,<br><br></div>Thanks for testing the procedures on GaAs. As to difference between the Wannier centers with and without SOC: did you carry out maximal localization procedure in both cases (with and without SOC)? What if you simply do a single-step projection, i.e., set num_iter=0 ? Then the projected localized functions will not hybridize with each other, and should have the symmetry of atomic orbitals (as long as these initial projections are well localized after some disentanglement procedure). <br><br></div>As to the meaning of the tight-binding model: I think it is meaningful as long as the low-energy bandstructure is the same as the DFT bandstructure and the Wannier functions are well localized. But I agree it may be more convenient in some cases if the Wannier functions are atomic-orbital like.<br><br></div>Thanks!<br><br></div>Jianpeng <br><div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 25, 2017 at 12:49 PM, Rubel, Oleg <span dir="ltr"><<a href="mailto:rubelo@mcmaster.ca" target="_blank">rubelo@mcmaster.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Jianpeng:<br>
<br>
I applied your procedure to WF in zinc-blende GaAs. Your procedure works. I had to do two extra steps<br>
$ cp wannier.inwf wannier.inwfup<br>
$ cp wannier.inwf wannier.inwfdn<br>
before running<br>
$ x w2w -so -up<br>
$ x w2w -so -dn<br>
<br>
The projections are:<br>
begin projections<br>
1:sp3(u,d)<br>
2:sp3(u,d)<br>
end projections<br>
<br>
At the end I got the following output with -so and -sp options:<br>
<br>
Final State<br>
WF centre and spread 1 ( -0.457506, 0.457588, 0.457476 ) 4.36038947<br>
WF centre and spread 2 ( 0.457535, 0.457378, -0.457522 ) 4.36341214<br>
WF centre and spread 3 ( -0.457477, -0.457464, -0.457434 ) 4.36424964<br>
WF centre and spread 4 ( 0.457409, -0.457465, 0.457515 ) 4.36405300<br>
WF centre and spread 5 ( 1.800179, -1.026445, -1.026433 ) 2.89020535<br>
WF centre and spread 6 ( 1.800146, -1.800199, -1.800158 ) 2.89044140<br>
WF centre and spread 7 ( 1.026459, -1.026417, -1.800189 ) 2.88989157<br>
WF centre and spread 8 ( 1.026408, -1.800139, -1.026430 ) 2.89026491<br>
WF centre and spread 9 ( -0.457482, 0.457621, 0.457463 ) 4.36089512<br>
WF centre and spread 10 ( 0.457601, 0.457520, -0.457595 ) 4.35805760<br>
WF centre and spread 11 ( -0.457628, -0.457579, -0.457578 ) 4.35685014<br>
WF centre and spread 12 ( 0.457550, -0.457579, 0.457661 ) 4.35678662<br>
WF centre and spread 13 ( 1.800180, -1.026454, -1.026421 ) 2.89005317<br>
WF centre and spread 14 ( 1.800147, -1.800197, -1.800166 ) 2.88979441<br>
WF centre and spread 15 ( 1.026464, -1.026412, -1.800185 ) 2.89031329<br>
WF centre and spread 16 ( 1.026412, -1.800135, -1.026420 ) 2.88997190<br>
Sum of centres and spreads ( 11.306397,-11.306377,-11.<wbr>306414 ) 58.00562975<br>
<br>
It can be compared to the calculation without -so and -sp<br>
<br>
Final State<br>
WF centre and spread 1 ( -0.000000, 0.000000, -0.000000 ) 1.91775879<br>
WF centre and spread 2 ( 0.000000, -0.000000, 0.000000 ) 5.85946098<br>
WF centre and spread 3 ( -0.000000, -0.000000, -0.000000 ) 5.85946098<br>
WF centre and spread 4 ( 0.000000, 0.000000, 0.000000 ) 5.85946098<br>
WF centre and spread 5 ( 1.413300, -1.413300, -1.413300 ) 1.61203280<br>
WF centre and spread 6 ( 1.413300, -1.413300, -1.413300 ) 3.81968345<br>
WF centre and spread 7 ( 1.413300, -1.413300, -1.413300 ) 3.81968345<br>
WF centre and spread 8 ( 1.413300, -1.413300, -1.413300 ) 3.81968345<br>
Sum of centres and spreads ( 5.653200, -5.653200, -5.653200 ) 32.56722487<br>
<br>
It looks like Wannier centres with SOC are not centred at atomic positions, which implies that their interpretation as atomic orbitals will be a stretch. I am not sure how useful the Hamiltonian will be then.<br>
<br>
<br>
Thank you<br>
Oleg<br>
--<br>
Oleg Rubel (PhD, PEng)<br>
Department of Materials Science and Engineering<br>
McMaster University<br>
JHE 359, 1280 Main Street West, Hamilton, Ontario L8S 4L8, Canada<br>
Email: <a href="mailto:rubelo@mcmaster.ca">rubelo@mcmaster.ca</a><br>
Tel: <a href="tel:%2B1-905-525-9140%2C%20ext.%2024094" value="+19055259140">+1-905-525-9140, ext. 24094</a><br>
Web: <a href="http://olegrubel.mcmaster.ca" rel="noreferrer" target="_blank">http://olegrubel.mcmaster.ca</a><br>
<div><div class="h5"><br>
> On Sep 25, 2017, at 10:52, Jianpeng Liu <<a href="mailto:ljp198795@gmail.com">ljp198795@gmail.com</a>> wrote:<br>
><br>
> Dear WIEN2k users and developers,<br>
><br>
> I am using wien2wannier to generate a tight-binding model for some rare-earth compound including Ce. Both spin-orbit coupling and spin-polarization have been included in my calculation. The following is my workflow:<br>
><br>
> 1) Finish a scf calculation, then prepare_w2wdir wannier<br>
><br>
> 2) init_w2w -up<br>
> At this step, I have to specify the energy windows.<br>
><br>
> 3) write_inwf -f wannier<br>
> At this step, I have to input the minimal and maximal band indices exported from the previous step, then specify the projections of the Wannier functions. I want to project onto 5 Ce d orbitals, but I made two identical copies of projections, i.e.:<br>
> Ce:d<br>
> Ce:d<br>
> One for spin-up, and the other for spin-dn, so in total there are 10 projections. I copied the generated file wannier.inwf to wannier.inwfup and wannier.inwfdn. Now both wannier.inwfup and .inwfdn have 10 projections, but actually there are only 5 orbitals without the spin degrees of freedom. Did I do it correctly?<br>
><br>
> 4) Then I wrote a wannier.win file by myself, in which there<br>
> are 10 *spinor* Wannier-function projections, and I set spinors=.true. (is this correct?).<br>
><br>
> 5)With all the inputs generated from the previous steps, I re-run the lapw1 calculation on the full BZ mesh, run wannier90 -pp to generate .nnkp file, run w2w, and then finally wannier90:<br>
> x lapw1 -up -p<br>
> x lapw1 -dn -p<br>
> x lapwso -up -p<br>
> x wannier90 -pp<br>
> x w2w -so -up -p<br>
> x w2w -so -dn -p<br>
> x wannier90 -so<br>
><br>
> The above is my workflow of wien2wannier. The code seems to be running well, but I am not sure if I did everything correctly. In particular I am worried that there are 2 identical copies of the orbital projections in both wannier.inwfup and .inwfdn. I would like to make sure that it is correct. I would greatly appreciate your help.<br>
><br>
> Best,<br>
> Jianpeng<br>
><br>
><br>
</div></div>> ______________________________<wbr>_________________<br>
> Wien mailing list<br>
> <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>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/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
<br>
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>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/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
</blockquote></div><br></div>