<div dir="ltr"><div><div><div>Hi,<br><br></div>The experimental values for the quadrupolar frequencies and the asymmetry parameter could be obtained with high accuracy for the material under study. The point is that in this work, these properties are being reported for the experimental structure and also for the equilibrium structure with forces smaller than 0.1 mRy. Actually the quadrupolar frequencies at experimental positions are in a very good agreement with their WIEN2k results, but the asymmetry parameters differ by 0.05321 and 0.11776, for the same structure. And more, when I try to reproduce the results at the supposed relaxed structure, both properties are quite different, the quadrupolar frequencies differ by more than 20 kHz. Additionally, with the relaxed structure, the forces in my case.scf are higher than 40 mRy after converge everything. Since I'm using the same structures, I was expecting to have the "same" numbers.<br><br></div>Thank you very much for your attention,<br><br></div>Ary Ferreira<br><div><div><div><br><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Apr 5, 2015 at 10:36 PM,  <span dir="ltr"><<a href="mailto:tran@theochem.tuwien.ac.at" target="_blank">tran@theochem.tuwien.ac.at</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
What do you mean by "exactly"? What are the differences between the<br>
published results and yours?<br>
<br>
F. Tran<div><div class="h5"><br>
<br>
On Sun, 5 Apr 2015, Ary Ferreira wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Please,<br>
<br>
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<br>
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<br>
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<br>
(lpot), and the k-points sampling.<br>
<br>
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<br>
parameter.<br>
<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>
- lpot = 4, 5, and 6;<br>
<br>
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<br>
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>
<br>
Thank you very much!<br>
<br>
Ary Ferreira<br>
<br>
--<br>
<a href="http://lattes.cnpq.br/8221674673413336" target="_blank">http://lattes.cnpq.br/<u></u>8221674673413336</a><br>
<br>
<br>
</blockquote></div></div>
______________________________<u></u>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.<u></u>at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.<u></u>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" target="_blank">http://www.mail-archive.com/<u></u>wien@zeus.theochem.tuwien.ac.<u></u>at/index.html</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><a href="http://lattes.cnpq.br/8221674673413336" target="_blank">http://lattes.cnpq.br/8221674673413336</a></div></div>
</div>