<div dir="ltr"><div><br></div>Dear Experts,<br><div><br> I am interested in phonon calculation for Zr. I found the following useful post in the QE-forum (please see below).<br><br> Can you please let me know where to put the values for nr1, nr2 and nr3. I mean is it to be put in <a href="http://scf.in">scf.in</a>, <a href="http://ph.in">ph.in</a> or <a href="http://q2r.in">q2r.in</a>?<br><br> And also if possible please let me know what does nr represents?<br><br>Thanking you,<br>Yours sincerely,<br>Krishnendu<br><br>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>-------<br><br><br><h1>[Pw_forum] Bug in ph.x sometimes occur in hexagonal phases</h1>
    <b>Uri Argaman</b> 
    <a href="mailto:pw_forum%40pwscf.org?Subject=Re:%20%5BPw_forum%5D%20Bug%20in%20ph.x%20sometimes%20occur%20in%20hexagonal%20phases&In-Reply-To=%3CCA%2B6Rk3xVk4TsRhe9RSNxUFRS2huAP_w6fbUDWFh6cmkm7bqt6w%40mail.gmail.com%3E" title="[Pw_forum] Bug in ph.x sometimes occur in hexagonal phases" target="_blank">argamanu at post.bgu.ac.il
       </a><br>
    <i>Thu Aug 27 12:56:57 CEST 2015</i>
    <ul><li>Previous message: <a href="http://qe-forge.org/pipermail/pw_forum/2015-August/107726.html" target="_blank">[Pw_forum] Bug in ph.x sometimes occur in hexagonal phases
</a></li><li>Next message: <a href="http://qe-forge.org/pipermail/pw_forum/2015-August/107732.html" target="_blank">[Pw_forum]  XSpectra and US or PAW pseudos?
</a></li><li> <b>Messages sorted by:</b> 
              <a href="http://qe-forge.org/pipermail/pw_forum/2015-August/date.html#107729" target="_blank">[ date ]</a>
              <a href="http://qe-forge.org/pipermail/pw_forum/2015-August/thread.html#107729" target="_blank">[ thread ]</a>
              <a href="http://qe-forge.org/pipermail/pw_forum/2015-August/subject.html#107729" target="_blank">[ subject ]</a>
              <a href="http://qe-forge.org/pipermail/pw_forum/2015-August/author.html#107729" target="_blank">[ author ]</a>
         </li></ul>
    <hr>  

<pre>Thank you Andrea
I think the phonon calculated with 13 dynamical matrices is wrong because
when I calculate the anisotropic coefficients of thermal expansion in
zirconium I do PW calculations as well as phonon calculations in several
lattice parameters. In some of them (not in all of them), QE calculate 13
dynamical matrices and in the other 12. When I use this data to calculate
the anisotropic coefficient of thermal expansion I got very strange results
that do not agree with some publications which also calculate it from DFT.
I do not know what is the problem. Maybe there is some problem in q2r.x.
Using nr1=30, nr2=30, nr3=48 indeed solve the problem. In addition,
different ecut and ecutrho also solve the problem.
I now calculate the coefficient of anisotropic thermal expansion in
zirconium again with ecut=50 and ecutrho=250. These calculations always
produce 12 dynamical matrices. These calculations takes a few days on my
machines. When I will get the results I will report if it makes any
difference.
Uri Argaman
Ben-Gurion University
Israel
</pre><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div>Dr. Krishnendu Mukherjee,</div></div><div><br></div><div>Principal Scientist,</div><div>CSIR-NML,</div><div>Jamshedpur.</div></div></div></div></div></div></div></div></div>
</div></div>