<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Dear all, I have done phonon calculations at the Gamma point to find the vibrational frequencies of a system I am working on and I got three negative frequencies; the results are:<div><br></div><div><div>&nbsp;q = ( &nbsp; &nbsp;0.000000000 &nbsp; 0.000000000 &nbsp; 0.000000000 )</div><div><br></div><div>&nbsp;**************************************************************************</div><div>&nbsp; &nbsp; &nbsp;omega( 1) = &nbsp; &nbsp; -11.303890 [THz] = &nbsp; &nbsp;-377.057178 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 2) = &nbsp; &nbsp; -11.228798 [THz] = &nbsp; &nbsp;-374.552397 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 3) = &nbsp; &nbsp; &nbsp;-1.493780 [THz] = &nbsp; &nbsp; -49.827127 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 4) = &nbsp; &nbsp; &nbsp; 1.499866 [THz] = &nbsp; &nbsp; &nbsp;50.030148 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 5) = &nbsp; &nbsp; &nbsp; 2.192955 [THz] = &nbsp; &nbsp; &nbsp;73.149113 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 6) = &nbsp; &nbsp; &nbsp;11.690342 [THz] = &nbsp; &nbsp; 389.947822 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 7) = &nbsp; &nbsp; &nbsp;15.929343 [THz] = &nbsp; &nbsp; 531.345701 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 8) = &nbsp; &nbsp; &nbsp;17.762801 [THz] = &nbsp; &nbsp; 592.503255 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega( 9) = &nbsp; &nbsp; &nbsp;17.814756 [THz] = &nbsp; &nbsp; 594.236307 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(10) = &nbsp; &nbsp; &nbsp;22.128875 [THz] = &nbsp; &nbsp; 738.139807 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(11) = &nbsp; &nbsp; &nbsp;24.754227 [THz] = &nbsp; &nbsp; 825.712121 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(12) = &nbsp; &nbsp; &nbsp;25.174421 [THz] = &nbsp; &nbsp; 839.728307 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(13) = &nbsp; &nbsp; &nbsp;25.229402 [THz] = &nbsp; &nbsp; 841.562251 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(14) = &nbsp; &nbsp; &nbsp;31.677488 [THz] = &nbsp; &nbsp;1056.647257 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(15) = &nbsp; &nbsp; &nbsp;32.931458 [THz] = &nbsp; &nbsp;1098.475192 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(16) = &nbsp; &nbsp; &nbsp;32.974208 [THz] = &nbsp; &nbsp;1099.901170 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(17) = &nbsp; &nbsp; &nbsp;37.529033 [THz] = &nbsp; &nbsp;1251.833794 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(18) = &nbsp; &nbsp; &nbsp;37.585396 [THz] = &nbsp; &nbsp;1253.713860 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(19) = &nbsp; &nbsp; &nbsp;38.689108 [THz] = &nbsp; &nbsp;1290.529726 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(20) = &nbsp; &nbsp; &nbsp;44.468725 [THz] = &nbsp; &nbsp;1483.317012 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(21) = &nbsp; &nbsp; &nbsp;44.490793 [THz] = &nbsp; &nbsp;1484.053106 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(22) = &nbsp; &nbsp; 100.618488 [THz] = &nbsp; &nbsp;3356.271501 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(23) = &nbsp; &nbsp; 100.705119 [THz] = &nbsp; &nbsp;3359.161186 [cm-1]</div><div>&nbsp; &nbsp; &nbsp;omega(24) = &nbsp; &nbsp; 103.337467 [THz] = &nbsp; &nbsp;3446.966862 [cm-1]</div><div><br></div><div>To check whether the first three freqnecies are the accoustic ones and not instabilities &nbsp;i applied dynmat.x with asr='crystal' and got:</div><div><br></div><div><div>&nbsp; mode &nbsp; [cm-1] &nbsp; &nbsp; [THz] &nbsp; &nbsp; &nbsp; IR</div><div>&nbsp; &nbsp; 1 &nbsp; -377.06 &nbsp;-11.3039 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 2 &nbsp; -374.55 &nbsp;-11.2287 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 3 &nbsp; &nbsp; &nbsp;0.00 &nbsp; &nbsp;0.0000 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 4 &nbsp; &nbsp; &nbsp;0.00 &nbsp; &nbsp;0.0000 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 5 &nbsp; &nbsp; &nbsp;0.00 &nbsp; &nbsp;0.0000 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 6 &nbsp; &nbsp;406.11 &nbsp; 12.1749 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 7 &nbsp; &nbsp;531.35 &nbsp; 15.9293 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 8 &nbsp; &nbsp;592.50 &nbsp; 17.7628 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp; 9 &nbsp; &nbsp;594.24 &nbsp; 17.8148 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;10 &nbsp; &nbsp;738.13 &nbsp; 22.1286 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;11 &nbsp; &nbsp;828.85 &nbsp; 24.8482 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;12 &nbsp; &nbsp;839.68 &nbsp; 25.1730 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;13 &nbsp; &nbsp;841.62 &nbsp; 25.2311 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;14 &nbsp; 1056.65 &nbsp; 31.6775 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;15 &nbsp; 1099.09 &nbsp; 32.9498 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;16 &nbsp; 1099.25 &nbsp; 32.9547 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;17 &nbsp; 1251.32 &nbsp; 37.5135 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;18 &nbsp; 1253.47 &nbsp; 37.5781 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;19 &nbsp; 1290.53 &nbsp; 38.6891 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;20 &nbsp; 1483.05 &nbsp; 44.4608 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;21 &nbsp; 1485.30 &nbsp; 44.5283 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;22 &nbsp; 3356.09 &nbsp;100.6131 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;23 &nbsp; 3359.58 &nbsp;100.7178 &nbsp; &nbsp;0.0000</div><div>&nbsp; &nbsp;24 &nbsp; 3446.96 &nbsp;103.3373 &nbsp; &nbsp;0.0000</div></div><div><br></div><div>As it is seen, the freqencies are very close but the thing and the system is stable! Howevere,I could not understand is about the first 5 frequencies: i still got the first two negative but in addition i have 3 more zero frequencies; does it mean we have five accoustic modes due to the symmetry of this particular system?</div><div><br></div><div><br></div><div>Thanks in advance</div><div><br></div><div>Elie Koujaes</div><div>University of Notts</div><div>NG7 2RD</div><div>UK</div><div><br></div></div>                                               </div></body>
</html>