<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.<br clear="all"></div><div><div><div><br>-- <br><div class="gmail_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></div></div>