<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="DE" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Dear experts, </p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have a problem running wannier90, which only occurs sporadically. I hence assume it is some problem with inconsistent symmetry operations:
</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><i>Wannier mode is: standalone <o:p></o:p></i></p>
<p class="MsoNormal"><i> <o:p></o:p></i></p>
<p class="MsoNormal"><i> -----------------<o:p></o:p></i></p>
<p class="MsoNormal"><i> *** Reading nnkp <o:p></o:p></i></p>
<p class="MsoNormal"><i> -----------------<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p></o:p></i></p>
<p class="MsoNormal"><i> Checking info from wannier.nnkp file<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p></o:p></i></p>
<p class="MsoNormal"><i> - Real lattice is ok<o:p></o:p></i></p>
<p class="MsoNormal"><i> - Reciprocal lattice is ok<o:p></o:p></i></p>
<p class="MsoNormal"><i> Something wrong! <o:p></o:p></i></p>
<p class="MsoNormal"><i> numk= 144 iknum= 360<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<o:p></o:p></i></p>
<p class="MsoNormal"><i> Error in routine pw2wannier90 (144):<o:p></o:p></i></p>
<p class="MsoNormal"><i> Wrong number of k-points<o:p></o:p></i></p>
<p class="MsoNormal"><i>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i> stopping ...<o:p></o:p></i></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The order I run the steps in is: </p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><i>pw.x < ncpp.scf.in | tee result.scf.out<o:p></o:p></i></p>
<p class="MsoNormal"><i>pp.x < pp.input.plot.rhops.in<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i>pw.x < ncpp.nscf.in | tee ncpp.nscf.out<o:p></o:p></i></p>
<p class="MsoNormal"><i>pw2critic.x < pw2critic.in<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i>wannier90.x -pp wannier.win<o:p></o:p></i></p>
<p class="MsoNormal"><i>$MPIEXEC $FLAGS_MPI_BATCH pw2wannier90.x < pw2wan.in<o:p></o:p></i></p>
<p class="MsoNormal"><i>wannier90.x wannier.win<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal">For the k-points of the scf calculation, I usually use an automatic grid, e.g.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><i>K_POINTS automatic<o:p></o:p></i></p>
<p class="MsoNormal"><i>6 6 4 0.0 0.0 0.0<o:p></o:p></i></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">While for the nscf calculation I use an exlicit list of k-points:</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><i>K_POINTS crystal<o:p></o:p></i></p>
<p class="MsoNormal"><i>144<o:p></o:p></i></p>
<p class="MsoNormal"><i> 0.00000000 0.00000000 0.00000000 6.944444e-03 <o:p></o:p></i></p>
<p class="MsoNormal"><i> 0.00000000 0.00000000 0.25000000 6.944444e-03<o:p></o:p></i></p>
<p class="MsoNormal"><i>……..<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal">I also tried to use the explicit list in the scf calculation, but that did not help.</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">wannier.win reads: </p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><i>num_wann=30<o:p></o:p></i></p>
<p class="MsoNormal"><i>num_iter=10000 ! default: 100<o:p></o:p></i></p>
<p class="MsoNormal"><i>conv_window=3<o:p></o:p></i></p>
<p class="MsoNormal"><i>conv_tol=1d-10 <o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i>begin projections<o:p></o:p></i></p>
<p class="MsoNormal"><i>random<o:p></o:p></i></p>
<p class="MsoNormal"><i>end projections<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i>#search_shells = 24 ! default: 12<o:p></o:p></i></p>
<p class="MsoNormal"><i>kmesh_tol = 1d-4 ! default: 1d-6<o:p></o:p></i></p>
<p class="MsoNormal"><i>mp_grid : 6 6 4<o:p></o:p></i></p>
<p class="MsoNormal"><i>begin unit_cell_cart<o:p></o:p></i></p>
<p class="MsoNormal"><i>bohr<o:p></o:p></i></p>
<p class="MsoNormal"><i>8.6888 0.0 0.0<o:p></o:p></i></p>
<p class="MsoNormal"><i>-4.3444 7.5247 0.0<o:p></o:p></i></p>
<p class="MsoNormal"><i>0.0 0.0 13.7208<o:p></o:p></i></p>
<p class="MsoNormal"><i>end unit_cell_cart<o:p></o:p></i></p>
<p class="MsoNormal"><i>begin atoms_frac<o:p></o:p></i></p>
<p class="MsoNormal"><i>Mg 0.0 0.0 0.0<o:p></o:p></i></p>
<p class="MsoNormal"><i>Mg 0.3333 0.6667 0.3678<o:p></o:p></i></p>
<p class="MsoNormal"><i>Mg 0.6667 0.3333 0.6322<o:p></o:p></i></p>
<p class="MsoNormal"><i>Sb 0.3333 0.6667 0.7747<o:p></o:p></i></p>
<p class="MsoNormal"><i>Sb 0.6667 0.3333 0.2253<o:p></o:p></i></p>
<p class="MsoNormal"><i>end atoms_frac<o:p></o:p></i></p>
<p class="MsoNormal"><i>begin kpoints<o:p></o:p></i></p>
<p class="MsoNormal"><i> 0.00000000 0.00000000 0.00000000 6.944444e-03 <o:p></o:p></i></p>
<p class="MsoNormal"><i> 0.00000000 0.00000000 0.25000000 6.944444e-03<o:p></o:p></i></p>
<p class="MsoNormal"><i>…. (with the correct number of 144 kpoints listed here)<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal">While pw2wan.in just reads:</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><i>&inputpp<o:p></o:p></i></p>
<p class="MsoNormal"><i>outdir='tmp/'<o:p></o:p></i></p>
<p class="MsoNormal"><i>prefix='Mg3Sb2'<o:p></o:p></i></p>
<p class="MsoNormal"><i>seedname='wannier'<o:p></o:p></i></p>
<p class="MsoNormal"><i>write_mmn=.true.<o:p></o:p></i></p>
<p class="MsoNormal"><i>write_amn=.true.<o:p></o:p></i></p>
<p class="MsoNormal"><i>/<o:p></o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal"><i><o:p> </o:p></i></p>
<p class="MsoNormal">Has anyone encountered this error already and/or knows how to fix it? Any suggestions are welcome!</p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks a lot and best regards, </p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Carl-Friedrich Schön (PhD student), </p>
<p class="MsoNormal">RWTH Aachen University, Germany</p>
</div>
</body>
</html>