
<html>

<head>

<style><!--

.hmmessage P

{

margin:0px;

padding:0px

}

body.hmmessage

{

font-size: 10pt;

font-family:微软雅黑

}

--></style>

</head>

<body class='hmmessage'>

<br><br><br>Dear Paolo, <br><br>The re-generated phonon from the force constant matrix are negative. <br>&nbsp; 0.177215&nbsp; 0.000000&nbsp; 0.000000<br>&nbsp; -45.4648&nbsp; -45.4648&nbsp; 119.2232&nbsp; 626.1196&nbsp; 626.1196&nbsp; 645.7030<br><br>comparing to the frequency got by ph.x <br><br>&nbsp;&nbsp; q = (&nbsp;&nbsp;&nbsp; 0.177215000&nbsp;&nbsp; 0.000000000&nbsp;&nbsp; 0.000000000 ) <br><br>&nbsp;**************************************************************************<br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 1) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.431524 [THz] =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14.394177 [cm-1]<br>&nbsp;(&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.129619 -0.022092&nbsp; 0.597312 -0.371775 ) <br>&nbsp;(&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.118762&nbsp; 0.048530&nbsp; 0.686479&nbsp; 0.000000 ) <br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 2) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.431524 [THz] =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14.394177 [cm-1]<br>&nbsp;( -0.000000&nbsp; 0.000000 -0.693560&nbsp; 0.118207&nbsp; 0.111631 -0.069481 ) <br>&nbsp;( -0.000000 -0.000000 -0.635471 -0.259673&nbsp; 0.128295 -0.000000 ) <br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 3) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.535384 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 117.928499 [cm-1]<br>&nbsp;( -0.661130 -0.445106&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 -0.000000 ) <br>&nbsp;( -0.247111 -0.551111&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 ) <br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 4) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 18.715499 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 624.285984 [cm-1]<br>&nbsp;( -0.000000 -0.000000&nbsp; 0.000000 -0.000000 -0.106238&nbsp; 0.066124 ) <br>&nbsp;(&nbsp; 0.000000&nbsp; 0.000000 -0.000000 -0.000000&nbsp; 0.992140&nbsp; 0.000000 ) <br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 5) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 18.715499 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 624.285984 [cm-1]<br>&nbsp;( -0.000000 -0.000000 -0.110727&nbsp; 0.058294&nbsp; 0.000000 -0.000000 ) <br>&nbsp;(&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.989559&nbsp; 0.071511&nbsp; 0.000000 -0.000000 ) <br>&nbsp;&nbsp;&nbsp;&nbsp; omega( 6) =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19.507679 [THz] =&nbsp;&nbsp;&nbsp;&nbsp; 650.710441 [cm-1]<br>&nbsp;(&nbsp; 0.080574&nbsp; 0.054886 -0.000000 -0.000000 -0.000000 -0.000000 ) <br>&nbsp;( -0.402244 -0.910327&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000&nbsp; 0.000000 ) <br>&nbsp;**************************************************************************<br><br>Maybe something wrong in the generation of interatomic force constant, but how to check it? Thanks a lot.<br><br>Sincerely,<br>Pan Hao <br>Tulane University <br>


<meta http-equiv="Content-Type" content="text/html; charset=unicode">

<meta name="Generator" content="Microsoft SafeHTML">

<style>

.ExternalClass .ecxhmmessage P

{padding:0px;}

.ExternalClass body.ecxhmmessage

{font-size:10pt;font-family:微软雅黑;}


</style>


<pre>  <br>Quoting HaoPan &lt;<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">adehaopan at hotmail.com</a>&gt;:<br><br>&gt;<i> On each q point I got, the phone frequencies are positive. But after <br></i>&gt;<i> the Fourier interpolation process, There exist negative frequencies.<br></i><br>hard to say what is happening. When you re-generate phonon from<br>the force constant matrix, you should find exactly the same phonon<br>frequencies for the wave-vector grid used to generate force constants<br>in real space. If not, something has gone wrong in the generation<br>of interatomic force constants<br><br>Paolo<br><i><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank"></a></i></pre>                                               </body>

</html>