<div dir="ltr">Hello Sanjeev,<div><br></div><div>The pseudopotentials I am using are <a class="" target="_blank" href="http://www.quantum-espresso.org/wp-content/uploads/upf_files/H.pz-vbc.UPF" style="color:rgb(175,16,29);font-size:13.333333015441895px;margin:0px;padding:0px;text-decoration:none;font-family:Arial,Helvetica,sans-serif;line-height:18px">H.pz-vbc.UPF</a> and <a class="" target="_blank" href="http://www.quantum-espresso.org/wp-content/uploads/upf_files/C.pz-vbc.UPF" style="color:rgb(175,16,29);font-size:13.333333015441895px;margin:0px;padding:0px;text-decoration:none;font-family:Arial,Helvetica,sans-serif;line-height:18px">C.pz-vbc.UPF</a></div>
<span style="font-size:13.333333015441895px;color:rgb(51,51,51);font-family:Arial,Helvetica,sans-serif;line-height:18px"></span><pre style="font-size:13.333333015441895px;margin-top:0px;margin-bottom:0px;padding:0px;color:rgb(51,51,51);line-height:18px">
 
Pseudopotential type: NORMCONS
Method: Von Barth-Car (direct fit)
Functional type: Perdew-Zunger (LDA) exch-corr
scalar relativistic</pre><pre style="font-size:13.333333015441895px;margin-top:0px;margin-bottom:0px;padding:0px;color:rgb(51,51,51);line-height:18px"><br></pre><pre style="font-size:13.333333015441895px;margin-top:0px;margin-bottom:0px;padding:0px;color:rgb(51,51,51);line-height:18px">
I have tried kgrid of 1*1*6 and 1*1*12, but both give negative frequencies.</pre></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Apr 14, 2014 at 12:18 AM, Sanjeev Gupta <span dir="ltr"><<a href="mailto:physics.skgupta@gmail.com" target="_blank">physics.skgupta@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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"><div><div class="h5">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>

</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">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></div></div>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">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><span class="HOEnZb"><font color="#888888"><br>
</font></span></blockquote></div><span class="HOEnZb"><font color="#888888"><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="line-height:19.49652862548828px;text-align:justify;color:rgb(68,68,68);font-family:'Trebuchet MS',Geneva,sans-serif">Michigan Technological University </span><br>

<span style="line-height:19.5px;text-align:justify;color:rgb(68,68,68);font-size:13px;font-family:'Trebuchet MS',Geneva,sans-serif">1400 Townsend Drive, </span><span style="line-height:19.5px;text-align:justify;color:rgb(68,68,68);font-size:13px;font-family:'Trebuchet MS',Geneva,sans-serif">Houghton </span><br style="line-height:19.5px;text-align:justify;color:rgb(68,68,68);font-size:13px;font-family:'Trebuchet MS',Geneva,sans-serif;margin:0px;padding:0px">

<span style="line-height:19.5px;text-align:justify;color:rgb(68,68,68);font-size:13px;font-family:'Trebuchet MS',Geneva,sans-serif">MI 49931, USA </span><br>------------------------------------</div>
</div></div></div>
</font></span></div>
<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></div>