<div dir="ltr"><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">Hello All!</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">I was wondering if any of you have encountered the following issue when fitting wannier functions in systems with specific symmetry.</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">I'm dealing with a magnetic system (but in principle, this could be any system) which has PT-symmetry. Since [PT, H] = 0, but in general, PT does not commute with other observables, the spectrum of H should be of doubly-degenerate bands (much like Kramers' degeneracy). Indeed a self-consistent VASP run converges to such a spectrum remarkably well, with the splitting between bands less than 1meV.</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">The trouble begins when attempting to fit Wannier functions to this system. While overall the fit is of good quality (see attached spread information; and band-structure comparisons) in many places along the high-symmetry lines the splitting between bands amounts to >10meV which is unacceptable in terms of preserving the symmetry.</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">More to the point the Wannier90.eig file <b>DOES </b>carry the symmetry rather well (most eigenvalues differ by less than 1meV), but then throughout the band-structure this is no longer the case.</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">What would you recommend in this instance?</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif"> </p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">I was thinking as well about the possibility of imposing this symmetry directly. Since any Bloch band \phi(k) has a PT-partner PT \phi(k), and these are naturally degenerate, then it would make a lot of sense (at least, superficially) to fit the Wannier bands to <i>half</i> the spectrum, and then reconstruct the other half by applying PT on the eigenvectors composed of the Wannier basis.</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">Is such a thing at all possible?</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif"> </p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">Yours thankful,</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">Daniel Kaplan</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">Dept. of Condensed Matter Physics</p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;font-size:11pt;font-family:Calibri,sans-serif">Weizmann Institute of Science</p></div>