<div dir="ltr"><div><div><div><div><div><div><div><div><div><div>Dear <span name="Thierry Tsafack" class="">Thierry Tsafack and QE users, <br><br></span></div><span name="Thierry Tsafack" class="">I actually have similar problems as Theirry. I am recently trying to calculate phonon dispersion relation for two dimensional materials, but always got negative frequencies. <br>
<br></span></div><div><span name="Thierry Tsafack" class="">First, I am pretty sure my structure is fully relaxed. It is also stable based on the literature report (they use VASP for calculation). <br></span></div><div><span name="Thierry Tsafack" class="">I used norm-conserving type of pseudopotential. <br>
</span></div><div><span name="Thierry Tsafack" class="">The negative frequency appears near the zone-center. At zone center the frequency is exactly zero.<br></span></div><div><span name="Thierry Tsafack" class=""><br></span></div>
<span name="Thierry Tsafack" class="">I've tried to resolve the problem using the following attempts:<br></span></div><span name="Thierry Tsafack" class="">1. Change different ASR<br></span></div><span name="Thierry Tsafack" class="">2. Increase ecutrho to from 240 to 480 Ry (8 times ecut)<br>
</span></div><div><span name="Thierry Tsafack" class="">3. Increase ecutwfc from 45 to 90<br></span></div><span name="Thierry Tsafack" class="">4. Increase k-point mesh size from 16*16*1 to 32*32*1<br></span></div><span name="Thierry Tsafack" class="">5. Increase tr2_ph from 1e-18 to 1e-20<br>
</span></div><span name="Thierry Tsafack" class="">6. Increase conv_thr from 1e-8 to 1e-12<br><br></span></div><span name="Thierry Tsafack" class="">Since the calculations are pretty time consuming, it is difficult to try even smaller convergence threshold.<br>
<br></span></div><span name="Thierry Tsafack" class="">I really hope someone can help with this issue.<br><br></span></div><span name="Thierry Tsafack" class="">Thanks,<br>Hua Bao<br></span><div><div><span name="Thierry Tsafack" class=""><br>
<br></span><div><span name="Thierry Tsafack" class=""><br><br></span></div></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Nov 2, 2013 at 4:03 AM, Thierry Tsafack <span dir="ltr"><<a href="mailto:tr.thierry@rice.edu" target="_blank">tr.thierry@rice.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
Dear fellow QE users,<br>
<br>
I hope this email finds you well.<br>
<br>
Could you please help me fix the issues below?<br>
<br>
1-) The executable ph.x requires the knowledge of nq1 nq2 nq3 to<br>
determine the density of the Monkhorst-Pack grid, i.e, number of<br>
points and their coordinates.  What is the linkage between nq1 nq2 nq3<br>
and the number of points in the grid?<br>
<br>
Everything else being equal, 3x3x1 yielded 5 points, 4x4x1 yielded 10<br>
points, 5x5x1 yielded 13 points, 6x6x1 yielded 20 points, 7x7x1<br>
yielded 25 points.<br>
<br>
For benchmark purposes, how do I know how many points nq1 x nq2 x 1<br>
will produce without doing the calculations?<br>
<br>
  2-) How to go about a situation in which ph.x produces w(q_i)>0 for<br>
all i in the above grid, but q2r.x and matdyn.x produce some w(q_i)<0,<br>
for all choices of the ASR or ZASR?  Assume the geometry relaxation<br>
was done as well as the scf calculation and you know that the<br>
structure is stable.<br>
<br>
3-) Recover =.true. in QE-4.3 always fails saying:"from checkallsym :<br>
error".  How to go around it?<br>
<br>
Thanks in advance for your clues, if your answer necessitates further<br>
details, let me know.<br>
<br>
<br>
~Thierry T.<br>
MEMS Rice University, Houston, TX<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></div>