<div dir="ltr">what is kind of pseudopotential?<div>How much k-point?</div><div><br></div><div><br></div><div>wishes</div><div>sanjeev</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Apr 13, 2014 at 9:25 AM, 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"><div><br></div>-- <br><div><span style="font-size:13px;font-family:Tahoma">With Best Regards,</span></div><div><span style="font-size:13px;font-family:Tahoma"><br></span></div><div>
<span style="font-size:13px;font-family:Tahoma">------------------------------</span><span style="font-size:13px;font-family:Tahoma">------</span></div><div><div style="font-family:Tahoma;font-size:13px"><div style="font-size:13px">
<div style="font-size:13px">Dr. Sanjeev Kumar Gupta<br>Fulbright Post-Doctoral Scholar<br>Dept. of Physics<br><span style="color:rgb(68,68,68);font-family:'Trebuchet MS',Geneva,sans-serif;line-height:19.49652862548828px;text-align:justify;background-color:rgb(255,255,255)">Michigan Technological University </span><br>
<span style="font-size:13px;color:rgb(68,68,68);font-family:'Trebuchet MS',Geneva,sans-serif;line-height:19.5px;text-align:justify;background-color:rgb(255,255,255)">1400 Townsend Drive, </span><span style="font-size:13px;color:rgb(68,68,68);font-family:'Trebuchet MS',Geneva,sans-serif;line-height:19.5px;text-align:justify;background-color:rgb(255,255,255)">Houghton </span><br style="font-size:13px;margin:0px;padding:0px;color:rgb(68,68,68);font-family:'Trebuchet MS',Geneva,sans-serif;line-height:19.5px;text-align:justify;background-color:rgb(255,255,255)">
<span style="font-size:13px;color:rgb(68,68,68);font-family:'Trebuchet MS',Geneva,sans-serif;line-height:19.5px;text-align:justify;background-color:rgb(255,255,255)">MI 49931, USA </span><br>------------------------------------</div>
</div></div></div>
</div>