<div dir="ltr"><div><div><div><div><div><div><div><div><div><div>Dear WIEN2k users and developers,<br><br></div>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:</div><div><span style="font-family:monospace,monospace"><br></span></div><div><span style="font-family:monospace,monospace">1) Finish a scf calculation, then </span><span style="font-family:monospace,monospace"><code>prepare_w2wdir wannier<br></code></span></div></div><span style="font-family:monospace,monospace"><code><br></code></span></div><div><span style="font-family:monospace,monospace"><code>2) init_w2w -up</code><code><br></code></span></div><span style="font-family:monospace,monospace"><code>  At this step, I have to specify the energy windows.<br></code></span></div><span style="font-family:monospace,monospace"><code><br></code></span></div><div><span style="font-family:monospace,monospace"><code>3) write_inwf -f wannier<br></code></span></div><span style="font-family:monospace,monospace"><code>  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></code></span></div><span style="font-family:monospace,monospace"><code>Ce:d<br></code></span></div><span style="font-family:monospace,monospace"><code>Ce:d<br></code></span></div><span style="font-family:monospace,monospace"><code>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></code></span></div><span style="font-family:monospace,monospace"><code>4) Then I wrote a wannier.win file by myself, in which there <br> </code></span><div><div><div><div><div><div><div><div><div><div><span style="font-family:monospace,monospace"><code></code></span><div><span style="font-family:monospace,monospace">are 10 *spinor* Wannier-function projections, and I set spinors=.true. (is this correct?).</span></div><div><span style="font-family:monospace,monospace"><br></span></div><div><span style="font-family:monospace,monospace">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></span></div><div><span style="font-family:monospace,monospace">  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></span></div><div><span style="font-family:monospace,monospace"><br></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">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.</font></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif"><br></font></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Best,</font></span></div><div><span style="font-family:monospace,monospace"><font face="arial,helvetica,sans-serif">Jianpeng</font><br></span></div><div><span style="font-family:monospace,monospace"><br></span></div><div><br></div></div></div></div></div></div></div></div></div></div></div></div>