<div dir="ltr">Hello Andrea,<div><br></div><div>I am not using 5.0.3 patch to QE 5.0.2 right now and I guess this is the issue. I will try using it today and will see how things go.</div><div><br></div><div>Thanks a lot,</div>
<div>Ankit Jain</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Apr 14, 2014 at 3:24 AM, Andrea Dal Corso <span dir="ltr"><<a href="mailto:dalcorso@sissa.it" target="_blank">dalcorso@sissa.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Did you apply the 5.0.3 patch to QE 5.0.2? Without the patch, the phonon<br>
code gives wrong frequencies for some q points.<br>
<br>
HTH,<br>
<br>
Andrea<br>
<div class=""><br>
<br>
On Mon, 2014-04-14 at 01:19 -0400, ankit jain wrote:<br>
> Hello Sanjeev,<br>
><br>
> The pseudopotentials I am using are<br>
</div>> H.pz-vbc.UPF<<a href="http://www.quantum-espresso.org/wp-content/uploads/upf_files/H.pz-vbc.UPF" target="_blank">http://www.quantum-espresso.org/wp-content/uploads/upf_files/H.pz-vbc.UPF</a>><br>
>  and C.pz-vbc.UPF<<a href="http://www.quantum-espresso.org/wp-content/uploads/upf_files/C.pz-vbc.UPF" target="_blank">http://www.quantum-espresso.org/wp-content/uploads/upf_files/C.pz-vbc.UPF</a>><br>
<div><div class="h5">><br>
><br>
> Pseudopotential type: NORMCONS<br>
> Method: Von Barth-Car (direct fit)<br>
> Functional type: Perdew-Zunger (LDA) exch-corr<br>
> scalar relativistic<br>
><br>
><br>
> I have tried kgrid of 1*1*6 and 1*1*12, but both give negative frequencies.<br>
><br>
><br>
><br>
> On Mon, Apr 14, 2014 at 12:18 AM, Sanjeev Gupta<br>
> <<a href="mailto:physics.skgupta@gmail.com">physics.skgupta@gmail.com</a>>wrote:<br>
><br>
> > what is kind of pseudopotential?<br>
> > How much k-point?<br>
> ><br>
> ><br>
> > wishes<br>
> > sanjeev<br>
> ><br>
> ><br>
> > On Sun, Apr 13, 2014 at 9:25 AM, Ankit <<a href="mailto:ankitjain.me.iitk@gmail.com">ankitjain.me.iitk@gmail.com</a>>wrote:<br>
> ><br>
> >> 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>
> >><br>
> ><br>
> ><br>
> ><br>
> > --<br>
> > With Best Regards,<br>
> ><br>
> > ------------------------------------<br>
> > Dr. Sanjeev Kumar Gupta<br>
> > Fulbright Post-Doctoral Scholar<br>
> > Dept. of Physics<br>
> > Michigan Technological University<br>
> > 1400 Townsend Drive, Houghton<br>
> > MI 49931, USA<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>
> ><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>
<br>
--<br>
</div></div>Andrea Dal Corso                    Tel. 0039-040-3787428<br>
SISSA, Via Bonomea 265              Fax. 0039-040-3787249<br>
I-34136 Trieste (Italy)             e-mail: <a href="mailto:dalcorso@sissa.it">dalcorso@sissa.it</a><br>
<div class="HOEnZb"><div class="h5"><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>
</div></div></blockquote></div><br></div>