<div dir="ltr">Hi Pablo,<div><br></div><div>Upon a quick glance your inputs look reasonable to me (although your pw.x convergence thresholds are really really unreasonably too tight) so this is just a wild guess but it could happen that the relaxation of the cluster converged to a saddle point of the potential energy surface instead of a local minimum (maybe due to the high initial symmetry of your structure).</div><div><br></div><div>If this is indeed the reason you could try to displace the cluster along the mode with imaginary frequency to gently push it down the hill. Another approach would be to do a short MD or an even simpler approach would be to apply random and small displacements to the initial atomic positions just to remove the high symmetry. I'm not sure if it is related but there are hydrogens on the surface of your cluster that repel each other: it was suggested in some early papers that surface reconstruction (=removal of those hydrogens and formation of long Si-Si bonds) could be favorable on these tiny clusters, too (similar to the bulk case).</div><div><br></div><div><div>Again, I'm not sure if this is indeed the reason but maybe it is worth trying!</div></div><div><br></div><div>Cheers,</div><div>Marton Voros</div><div><br></div><div>--</div><div>Materials Science Division</div><div>Argonne National Laboratory</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Apr 14, 2017 at 8:47 AM, Pablo García Risueño <span dir="ltr"><<a href="mailto:garcia.risueno@gmail.com" target="_blank">garcia.risueno@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"><div><div><div>Dear Espresso community<br><br></div>I come to you because I find a problem that I cannot solve. I am running phonon calculations in different Si-H clusters. Some of them (the smallest ones) work well, but for the rest I always find negative frequencies beyond the 6th frequency (at q=0, I deal with isolated systems). <br><br></div>I obtain very nice relaxations (e.g. << Total force = 0.000093 >>) but ph.x gives some negative frequencies beyond the 6th one:<br><br> freq ( 6) = -2.042467 [THz] = -68.129354 [cm-1]<br> freq ( 7) = -2.040393 [THz] = -68.060182 [cm-1]<br> freq ( 8) = -2.040381 [THz] = -68.059778 [cm-1]<br> freq ( 9) = 0.396144 [THz] = 13.213941 [cm-1]<br><br></div>When I run dynmat.x, the results get worse instead of improving:<br><br> 10 -134.72 -4.0388 0.0499<br> 11 -134.72 -4.0388 0.0499<br> 12 -77.44 -2.3217 0.0000<br> 13 -77.44 -2.3217 0.0000<br> 14 -76.18 -2.2838 0.0000<br> 15 -76.18 -2.2838 0.0000<br> 16 -76.18 -2.2838 0.0000<br> 17 -71.69 -2.1492 0.0201<br> 18 -71.69 -2.1492 0.0201<br> 19 -71.69 -2.1492 0.0201<br> 20 -58.72 -1.7604 0.0000<br> 21 -52.87 -1.5849 0.0000<br> 22 -45.50 -1.3642 0.0000<br> 23 -45.50 -1.3642 0.0000<br> 24 -45.50 -1.3642 0.0000<br> 25 -38.77 -1.1622 0.0000<br> 26 -38.77 -1.1622 0.0000<br> 27 2.10 0.0628 0.0103<br> 28 2.10 0.0628 0.0103<br> 29 2.10 0.0628 0.0103<br> 30 31.92 0.9568 0.0000<br><br><div><div><br>I have tried many recommendations that I found at the forum (increase the box size, reduce the conv_thr, try other pseudopotentials...) but everything failed. Could you give me a clue on how this problem could be solved? I am attaching a document with the input files and part of the output informations.<br><br></div><div>Thank you very much. Best regards.<span class="HOEnZb"><font color="#888888"><br clear="all"></font></span></div><span class="HOEnZb"><font color="#888888"><div><div><div><br>-- <br><div class="m_-1922306817625235043gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">--<br><br>Dr. Pablo García Risueño<br><br>Institut für Physikalische Chemie, Universität Hamburg, Grindelallee 117, 20146 Hamburg<br><br>Tel. +49 040 42 83 84 82 7</div></div></div></div></div></div>
</div></div></div></font></span></div></div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div>