<div dir="ltr">Dear Wannier90 developers,<div><br></div><div>I have a strange problem I cannot figure out concerning the execution of wannier90 with gamma_only = true.  I am forming localized molecular orbitals for perylene dimer pairs, with one orbital (either the monomeric HOMO or LUMO) on each monomer.  I am using only the gamma point for these calculations.  For the electronic structure, I use VASP.6.3.0, and this is linked with wannier90 v3.1.0.</div><div><br></div><div>If I use vasp_std followed by a normal wannier calculation without gamma_only = true, the calculation completes in full: disentanglement and wannierisation both complete without issue, and I can plot the resultant molecular orbitals in the form of *.xsf files.  However, since I only sample with the gamma point, I figured I could compare this calculation with vasp_gam and wannier90 computed with gamma_only = true.  If I do this, the disentanglement and wannierisation phases still complete just fine, and a similar final spread is reached for the wannierisation phase; however, I cannot plot the orbital that comes out of the calculation within an *.xsf file.  The code crashes and leaves an error message.</div><div><br></div><div>I have attached some salient example files in ZIP folders to this email.  Folder 'std' corresponds to the initial start using vasp_std and with gamma_only = false in wannier90.  Folder 'gam' corresponds to the initial start using vasp_gam with gamma_only = true in wannier90.  Within these folders are e.* files; for the folder 'gam' you will find the error message from above in this e.* file.  Note that I had to remove the *.mmn files in order to meet the file size requirement.</div><div><br></div><div>I could not include the UNK* files for these calculations due to their immense sizes, but I will say that the UNK* file for the gamma_only = true calculation is half that of the gamma_only = false calculation, which makes sense if only the former is purely real.  For the std calculation, since I could plot the resultant *.xsf file, I have provided an image of the target orbital in lieu of the actual *.xsf file.  The target orbital is correct.  You might notice that the spreads of these calculations are quite large, but that's because I'm finding localized molecular orbitals, not atomic-like orbitals.  So this is not an error, if that was something you might think.</div><div><br></div><div>Any help would be appreciated.  It seems strange to me that I cannot plot the orbitals if I use the gamma_only version of the code, so I'm wondering if I'm missing an input flag, or if there is some compilation error on my end.  Or maybe I'm mistaken that plotting is possible for gamma_only = true.</div><div><br></div><div>Many thanks,</div><div>Peyton Cline<br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><font size="1"><br></font></div><div><font size="1">-------------------------------------------------------</font><span style="font-size:x-small">-------------</span></div><div><font size="1">R. Peyton Cline</font></div><div><font size="1">Pronouns: he/him</font></div><div><div style="color:rgb(34,34,34)"><font size="1">Postdoctoral Associate</font></div><div style="color:rgb(34,34,34)"><font size="1">Prof. Joel Eaves' Group</font></div><div style="color:rgb(34,34,34)"><font size="1">University of Colorado Boulder</font></div><div style="color:rgb(34,34,34)"><font size="1">Department of Chemistry</font></div></div><div dir="auto" style="color:rgb(34,34,34)"><font color="#8B8B8B" face="Helvetica" size="1"></font></div><div style="margin:2px 0px 0px;color:rgb(34,34,34)"></div></div></div></div></div>