<div dir="ltr">Dear all,<div><br></div><div>When I calculate the density of the lowest conduction states of NaI (see attachment), I get different results for using k-points from the full and irreducible Brillouin zones. The difference is smaller when the k-mesh is allowed to use less symmetries, e.g. only the time reversal symmetry. There is no difference if the 'kpoint.dat' file contains the full BZ. This suggests the code I use for generating IBZ points might be in error, but I get the same results with both abinit and elk (they deliver different, but equivalent IBZ points). Could it be that the DOS calculation is very sensitive to small numerical inaccuracies and this leads to slightly different results for full and irreducible BZs? The DOS related part of my input file is the following:</div><div><br></div><div><div>dos true</div><div>dos_kmesh 20</div><div>dos_adpt_smr true<br></div><div>wanint_kpoint_file false # true for the second run<br></div></div><div><br></div><div>The difference is also present with dos_kmesh 50. Please ask for other details I'm not showing. Any suggestions are welcome.</div><div><br></div><div>Thanks,</div><div>Raul Laasner</div><div><br clear="all"><div><br></div>-- <br><div><div dir="ltr">Raul Laasner<br>Institute of Physics<br>University of Tartu<br>Ravila 14c, 50411, Estonia<br>e-mail: <a href="mailto:raullaasner@gmail.com" target="_blank">raullaasner@gmail.com</a><br>cell: <a href="tel:%28%2B372%295182268" value="+3725182268" target="_blank">(+372)5182268</a></div></div>
</div></div>