<div dir="ltr"><div>For which q-point are you getting the imaginary modes? If it is at gamma, then it may be related to the instability in the structure, if it is at other q-points, may be lattice parameter is not properly optimized. How large is the stress you are getting at the end of the vc-relax calculation (when the code recalculates the total energy, stress and forces with the optimized lattice parameters & atomic positions). Also have you checked the convergence of the phonon modes with ecutwfc because typically the phonons converge much slowly with the size of the basis set than other parameters?<br>
<br></div><div>With regards,<br><br></div>Prasenjit<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 13 April 2014 21:55, Ankit <span dir="ltr"><<a href="mailto:ankitjain.me.iitk@gmail.com" target="_blank">ankitjain.me.iitk@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello QE developers and users,<br>
<br>
I am trying to calculate phonon dispersion using ph.x and matdyn.x<br>
routines of QE and I am stuck woth negative frequencies for a while now.<br>
<br>
I have a 1-D chain of polyethylene which has 6 atoms in the unit-cell.<br>
The polymer chain is aligned in z-direction and I have added vacuum on<br>
x- and y- direction. I did the structure relaxation and I am getting<br>
total force after relaxtion as 0.000020, which I guess is sufficiently<br>
small.<br>
<br>
After relaxation, when I am trying to calculate phonon dispersion using<br>
ph.x and matdyn.x I am getting netting frequencies for some of the<br>
phonon modes. I  searched online and on QE mailing list and tried couple<br>
of different things:<br>
<br>
1. The negative frequencies are appearing for non-gamma point as well,<br>
so I guess its not a ASR issue. Anyway, I am using ASR but even after<br>
this I am getting negative frequencies.<br>
2. I reduced my tr2_ph and conv_thr to 1.0d-14 and 1.0d-15 which I guess<br>
are very small as suggested on QE FAQs.<br>
3. I checked phonon mode shapes for modes with negative frequencies and<br>
they do not look like rotational mode.<br>
4. I tried changing number of k points from 6 to 12 in z direction but<br>
this is also not the cause.<br>
<br>
I am really stuck here and not sure about the cause of negative<br>
frequencies. With vc-relax, my structure is getting relaxed which I<br>
guess means my structure is stable. But then I am not sure why I am<br>
getting negative frequencies.<br>
<br>
I really appreciate any help or suggestions.<br>
<br>
Thanks,<br>
<br>
Ankit Jain<br>
IIT Indore,<br>
Indore,<br>
India<br>
<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr">PRASENJIT GHOSH,<br>IISER Pune,<br>Dr. Homi Bhabha Road, Pashan<br>Pune, Maharashtra 411008, India<br><br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790</div>

</div>