Dear All,<div><br></div><div>i`m having a problem in understanding the output of wannier90 (file hr.dat).</div><div><br></div><div>i have made a calculation for a molecular crystal.</div><div><br></div><div>the unit cell consists of molecular sites (A and B) along the stacking direction a (with a lattice vector) as illustrated in the following:</div>
<div><br></div><div>-----> a</div><div><br></div><div>|o   o | o   o | ......</div><div> A   B  A  B</div><div><br></div><div>(here i`m plotting schematically two unit cells)</div><div><br></div><div>Let label the molecules in this way:</div>
<div><br></div><div><div>-----> a</div><div><br></div><div>|o    o | o    o | ......</div><div> A1 B1 A2 B2</div><div><br></div><div>now the bonds B1-A2 are shorter while bonds A2-B2 are longer (the system is dimerized).</div>
<div><br></div><div>the HOMOs (highest occupied molecular orbitals) of the two molecules A,B contained in the unit cell hybridize, giving rise to bonding and antinbonding bands around the Fermi level. then at Fermi level i have two states.</div>
</div><div><br></div><div>the calculations performed with PW and W90 looks fine and the agreement of PW and W90 bands is good.</div><div><br></div><div>once i plot using xcrystalden the Wannier functions of the such states at Fermi one is located at A2 and the other one is located to B1 (basically one is in the unit cell and the other one is located at the NN unit cell along a).</div>
<div><br></div><div>now i go to look at the hr.dat file and i find:</div><div><br></div><div>    0    0    0    1    1    3.334228    0.000000</div><div>    0    0    0    2    1    0.222024    0.000000</div><div>    0    0    0    1    2    0.222024    0.000000</div>
<div>    0    0    0    2    2    3.334221    0.000000</div><div><br></div><div><div>    1    0    0    1    1    0.000243    0.000000</div><div>    1    0    0    2    1    0.203498    0.000000</div><div>    1    0    0    1    2    0.000122    0.000000</div>
<div>    1    0    0    2    2    0.000242    0.000000</div></div><div><br></div><div><div>   -1    0    0    1    1    0.000243    0.000000</div><div>   -1    0    0    2    1    0.000122    0.000000</div><div>   -1    0    0    1    2    0.203498    0.000000</div>
<div>   -1    0    0    2    2    0.000242    0.000000</div><div><br></div><div>i have some doubts in understanding the hopping parameters.</div><div><br></div><div>i see the kinetic terms are different 0.222024/0.203498 as one would expect in a dimerized structure. this is fine.</div>
<div><br></div><div>the values 3.334228 are the on-site energies of the tight-binding problem.</div><div><br></div><div>my problem is the following.</div></div><div>is 0.222024 referring to the the hopping between A1 and B1 (also between A2 and B2)  or between B1 and A2?</div>
<div><br></div><div>because the Wannier functions plotted in real space are sitting between adjacent different unit cells i would expect 0.222024 referring to hopping between B1 and A2 and 0.203498  referring to hopping between A1 and B1 (also between A2 and B2).</div>
<div>then at smaller distances would correspond larger hoppings and at larger distances would correspond smaller hoppings.</div><div><br></div><div>could you please confirm this way to interpret the file hr.dat? :-)</div>
<div><br></div><div>i thank you in advance for our time.</div><div><br></div><div>cheers,</div><div><br></div><div>Gianluca Giovannetti</div>