<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div><br></div>Dear Dr. Lorenzo Paulatto,<br><br></div> Thank you for your reply. While learning about the Density Functional Perturbation Theory I am practising the phonon calculation for bcc iron with Quantum espresso 6.2. So, I need kind help from you regarding the meaning of different parameters.<br><br></div> Your reply gave me the confidence that there was no Fortran language/ or LINUX command related problem. As per your suggestion I started to change different parameters to make run the code till end. First of all I decreased the values of conv_thr to 1.0e-16 and tr2_ph to 1.0e-16.<br><br></div> I realized that fe.dyn1, fe.dyn2, fe.dyn3 ..... are calculated one after another. And if in one of them if freq(1,2,3) becomes too negative the calculation stops. Please let me know if that is the case.<br><br></div> Then I checked the discussions at pw_forum and found this topic was discussed at <a href="http://qe-forge.org/pipermail/pw_forum/2014-August/104992.html">http://qe-forge.org/pipermail/pw_forum/2014-August/104992.html</a> As per the discussion when I increased the ecutrho from 200.0 to 400.0 the code ran till end, although there were negative frequencies. So, is it true that increasing the ecutrho helps in convergence of the code?<br></div><br></div>I have another question about the parameter mixing_beta . I did calculation with a value of 0.2. Should I increase or decrease the value to have positive values of frequencies.<br><br></div>Further I changed the K_POINTS to {automatic}  8 8 8  1 1 1. In future I plan to change it to 12 12 12 1 1 1 and to 16 16 16 1 1 1.<br><br></div>Thanking you,<br></div>Yours sincerely,<br></div>Krishnendu<br><div><div><div><div><div><div><div><div><div><div><div><br>================================================================<br>Message: 12<br>
Date: Tue, 26 Dec 2017 11:57:22 +0100<br>
From: Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>><br>
Subject: Re: [Pw_forum] Problem during phonon calculation of bcc iron<br>
        ferrite<br>
To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
Message-ID: <<a href="mailto:0f6d3673-f689-e709-b6d9-4e87649d5c1a@gmail.com">0f6d3673-f689-e709-b6d9-<wbr>4e87649d5c1a@gmail.com</a>><br>
Content-Type: text/plain; charset=utf-8; format=flowed<br>
<br>
On 26/12/17 10:20, Krishnendu Mukherjee wrote:<br>
> Fortran runtime error: End of file<br>
The file fe.dyn2 is probably empty, likely because the phonon<br>
calculation did not finish correctly.<br>
<br>
I think you should start from scratch, paying attention to detail,<br>
taking the necessary time to check that every step completes correctly<br>
or understanding why it doesn't.<br>
<br>
kind regards<br>
<br>
<br>
--<br>
Lorenzo Paulatto - Paris<br>
<br>
<br>
------------------------------<br>
<br><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></div></div></div></div></div></div></div></div></div></div>