<DIV>Dear all,</DIV>
<DIV>  I have a question about zero point energy calculation. According to a article called "Tutorial on calculating phonons: comparing the linear response and the small displacement methods", author states that we can calculate zero point energy use the results obtained with linear response method, i.e phonon code in pwscf. Calculated steps as follows(for example Al):</DIV>
<DIV>(1) > pw.x < Al.scf (This generates the self-consistent wave-function)</DIV>
<DIV>(2) > ph.x < Al.ph.in (This generate the dynamic matrices)</DIV>
<DIV>(3) > q2r.x < q2r.in (Obtaining force constant matrix)</DIV>
<DIV>(4) > matdyn.x < matdyn.in (Calculating the phonon dispersions)</DIV>
<DIV>(5) > matdyn.x < phdos.in (Calculating phonon density of states)<BR>and this will produce a file called "Al.phdos" which contains the phonon density of states. Now by typing:</DIV>
<DIV>    awk '{e=e+2*1/2}END{print e/33.357*4.1357/1000}' Al.phdos<BR>and you should see that the zero point energy. The question is, I do not understand the meaning of "{e=e+2*1/2}END{print e/33.357*4.1357/1000}", is this appropriate for the calculation of zero point energy? </DIV>
<DIV>    In addition, In mail list, by calculating phonons all over the Brillouin zone (in practice, on a regular grid), and zero point energy by calculating {\hbar\over 2} \sum_{{\bf q},\nu} \omega_{{\bf q},\nu}. I do not understand the meaning of {\hbar\over 2} \sum_{{\bf q},\nu} \omega_{{\bf q},\nu} well. Anyone who help me will be appreciated.</DIV>
<DIV> </DIV>
<DIV>Many thanks </DIV>
<DIV>Best regards</DIV>
<DIV>Lihui</DIV>
<DIV><BR>--<BR>
<DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV>
<DIV><FONT face="Arial" size="3">Lihui Ou</FONT></DIV>
<DIV><FONT face="Arial" size="3">PH.D Candidate in Electrochemistry                      </FONT></DIV>
<DIV><FONT face="Arial" size="3">College of Chemistry and Molecular Science            </FONT></DIV>
<DIV><FONT face="Arial" size="3">Wuhan University,<EM>430072</EM>,Hubei Province,China </FONT></DIV>
<DIV><EM><FONT face="Arial" size="3">E-mail:</FONT><A href="mailto:oulihui666@126.com" target="_blank"><FONT face="Arial" size="3">oulihui666@126.com</FONT></A></EM></DIV>
<DIV><FONT face="Arial" size="3">======================================</FONT></DIV></DIV></DIV><br><br><span title="neteasefooter"/><hr/>
<a href="http://www.yeah.net/?from=footer">没有广告的终身免费邮箱,www.yeah.net</a>
</span>