<div dir="ltr"><div><div>Please,<br><br></div><div>I'm trying to reproduce a result obtained with WIEN2k from a publication. I'm using the same structure, which is a monoclinic cell with the C2/c space group. So, I used the program SETSTRU to transform it from the space group C12/c1 to B112/b. Regarding the WIEN2k calculation, the authors reported the MT radii, the product of the minimal MT radius and the cut-off for the PW basis set Kmax (RmtxKmax), the maximal value of l in the expansion of the PW in spherical harmonics (lmax), the highest l in the expansion of the potential in lattice harmonics (lpot), and the k-points sampling.<br><br></div><div>I tested the following parameters (in the order listed below) for the convergence of the total energy, the main component of the electric field gradient tensor, and its asymmetry parameter.<br></div><div><br>- the k-points sampling using "x kgen" and selecting different number of k-points;<br>- the value of RmtxKmax from 6.0 to 9.0 in steps of 0.5;<br>- the cutoff for the potential (and charge density) Gmax from 12 to 20 in steps of 2;<br>- lmax = 12, 13, and 14;<br></div><div>- lpot = 4, 5, and 6;<br><br></div><div>Additionally, I also tried to reduce the MT radii by 0.05% and 0.1%, and tested again the RmtxKmax. However, I'm not being able to reproduce "exactly" the reported theoretical results. I just want to know what else can be done to improve the obtained observables. If there is some document I have missed, please I'm sorry.<br></div><div><br></div>Thank you very much!<br><br></div>Ary Ferreira<br clear="all"><div><div><div><div><br>-- <br><div><div dir="ltr"><a href="http://lattes.cnpq.br/8221674673413336" target="_blank">http://lattes.cnpq.br/8221674673413336</a></div></div>
</div></div></div></div></div>