<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><blockquote type="cite"><div>Dear all,</div></blockquote><div><br></div>Hi Lihui:</div><div><br><blockquote type="cite"> <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></div></blockquote><div><br></div>I am close to null in AWK, nor did I know anything about ".phdos" files and it took me a while to understand what's going on here.</div><div><br></div><div>1) I see that ".phdos" files contain two columns of data, both starting from zero, both positive, data in the first column equally spaced. Let's assume the first column is a table of frequencies, the second the corresponding dos.</div><div><br></div><div>2) You did not say what article you are referring to. You could have, but in any case GOOGLE helps. Here is, I think, the source of your troubles ...</div><div><a href="http://www.mcc.uiuc.edu/summerschool/2007/qmc/tutorials/Tutorial_phonons_alfe.pdf">http://www.mcc.uiuc.edu/summerschool/2007/qmc/tutorials/Tutorial_phonons_alfe.pdf</a></div><div><br></div><div>3) AWK parses text files line by line: "$1" and "$2" is the name that AWK givs to the first and second text fields it encounters (I did not remember this: I re-learn AWK once every 5 years, the last time was 5 mins ago here: </div><div><a href="http://www.gnu.org/software/gawk/manual/gawk.html">http://www.gnu.org/software/gawk/manual/gawk.html</a></div><div><br><blockquote type="cite"><div>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></blockquote><div><br></div>the awk line reported in the lecture notes you are reading is supposed to make an integral (sum) over frequency of the product of the frequency (first column, "$1"), and dos (second column, "$2"). (see Eq. 7 in the lecture notes and take the T->0 limit). can you see now what happened? (HINT the lecture notes contain a typo: dollars "$" are not displayed ...; the instruction in curly brackets following "print" must be a change of variables / normalization ...</div><div> <br><blockquote type="cite"> <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></blockquote><div><br></div>I do not understand if what you do not understand are basic concepts (how the zero-point energy is related to the dos, how BZ integrals are related to frequency integrals ...) or details, such as the syntax of AWK/LaTeX ... Please, try to clarify this to yourself in the frst place. If it is physics concepts that you do not understand, make an effort to formulate a specific question for each specific point you do not understand, and then revert to us with a specific question. If it is AWK/LaTeX, do like me: browse the net with google and/or take a walk to the nearest library and learn the tools you need (occasionally, we would also be pleased to help for such matters, but the help is in this case even less guaranteed the with physics).</div><div><br></div><div>SB</div><div><br></div><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">---</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Stefano Baroni - SISSA</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">&</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">DEMOCRITOS National Simulation Center - Trieste</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://stefano.baroni.me">http://stefano.baroni.me</a> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><br style="font-size: 12px; "></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><span class="Apple-style-span" style="font-size: 14px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br></span></font></div></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Please, if possible, don't</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">send me MS Word or PowerPoint attachments</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Why? See:</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://www.gnu.org/philosophy/no-word-attachments.html">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div><div><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br class="webkit-block-placeholder"></span></font></div></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"> </div><br></body></html>