<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Dear all,<br><br>After Prof. Stefano's comments, I will like to indicate that my aim<br>was just to give a contribution, hoping that it will help.<br><span style="font-weight: bold;">If I am wrong in any of my statements, please kindly suggest.</span><br><br>Furthermore, everything in my answer has been given with reference<br>and the reader should get acquainted with in order to solve his own <br>problem.<br><br>Sorry! <br><br>B. SITAMTZE<br><br>--- En date de : <b>Dim 21.8.11, Stefano Baroni <i><baroni@sissa.it></i></b> a écrit :<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;"><br>De: Stefano Baroni <baroni@sissa.it><br>Objet: Re: [Pw_forum] Re : Thermal conductivity and Gruneisen parameters<br>À: "PWSCF Forum" <pw_forum@pwscf.org><br>Date: Dimanche 21 août 2011,
9h23<br><br><div id="yiv2142409382"><div><blockquote type="cite"></blockquote><div>On Aug 20, 2011, at 11:45 AM, Bertrand SITAMTZE wrote:</div></div><div><div><br></div><div><blockquote type="cite"><table style="" border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font: inherit;" valign="top">Thermal conductivity could be obtained through velocity autocorrelation function <br>after running molecular dynamics (CP in QE).( See the book of Frenkel and Smith<br>on page 90).<br><br></td></tr></tbody></table></blockquote><div><table style="" border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font: inherit;" valign="top">I do not have Frenkel & Smith's book handy, but I would be surprised if it stated what you say. The calculation of the thermal conductivity through MD requires a proper definition of the energy current which is well known for classical systems (though it involves something more complicated than just
velocities), but not yet available, to the best of my knowledge, in DFT. If I am wring, please correct me.<br></td></tr></tbody></table></div><blockquote type="cite"><table style="" border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font: inherit;" valign="top"><br><br>As for Gruneisen parameters, one could calculate the specific heat through QHA and the bulk modulus with PW, then use the relationship between them. (See the book of L. Anderson on equations of state, starting from page 6).<br><br></td></tr></tbody></table></blockquote><div><table style="" border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font: inherit;" valign="top">Here, I beg to differ as well. First, the calculation of the constant-volume specific heat does not involve any Gruneisen parameters (gamma). The constant-pressure specific heat (Cp) certainly does, but in the sense that you need them in the first place in order to calculate it. Second, It would
seem rather innatural to calculate Cp just to estimate some of the parameters used to obtain it! In any case, the expression for Cp in terms of the gammas is an integral from which it would be impossible to obtain back their individual values. </td></tr></tbody></table></div><blockquote type="cite"><table style="" border="0" cellpadding="0" cellspacing="0"><tbody><tr></tr></tbody></table></blockquote><div><br></div></div></div>All in all, before asking what QE (or any other code, for that matter) can or cannot do, one should stop for a while and ask oneself what are the most elementary pieces of information needed to compute what we are interested in, and only then ask if any specific code can compute those specific pieces of information. This may take some time, but it is much more instructive than hoping that computer codes can ever be a substitute for our own understanding.<div><br></div><div>In the specific case of the thermal conductivity, one
has the choice to follow the quasi-harmonic or the Kubo-Greenwood routes. Nicola Marzari has already given some hints on the former. Dr. Simatze's comments (and my present reply) apply to the latter. As for Gruneisen parameters, you just have just to inspect their definition (gamma = v/omega d (omega) / d (v)), and realize that any code able to compute vibrational frequencies can with little effort calculate them ... QE can calculate phonon frequencies, I will leave the final conclusion to you ...</div><div><br></div><div>Hope this helps</div><div>Stefano B</div><div><br><div>
<span class="yiv2142409382Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; 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; font-size: medium;"><div><div style="margin: 0px; font-size: 14px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;">---</span></font></div><div style="margin: 0px; font-size: 14px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;">Stefano Baroni - SISSA</span></font><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"> </span></font><font class="yiv2142409382Apple-style-span" size="3"><span
class="yiv2142409382Apple-style-span" style="font-size: 12px;">&</span></font><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"> </span></font><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;">DEMOCRITOS National Simulation Center - Trieste</span></font></div><div style="margin: 0px; font-size: 14px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"><a rel="nofollow" target="_blank" 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: 0px; font: 12px Helvetica; min-height: 14px;"><br style="font-size: 12px;"></div><div style="margin: 0px; font: 12px Helvetica; min-height: 14px;"><span
class="yiv2142409382Apple-style-span" style="font-size: 14px;"><div style="margin: 0px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-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: 0px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"><br></span></font></div></span></div><div style="margin: 0px; font-size: 14px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;">Please, if possible, don't</span></font><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"> </span></font><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span"
style="font-size: 12px;">send me MS Word or PowerPoint attachments</span></font></div><div style="margin: 0px; font-size: 14px;"><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;">Why? See:</span></font><font class="yiv2142409382Apple-style-span" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"> </span></font><font class="yiv2142409382Apple-style-span" color="#0023e9" size="3"><span class="yiv2142409382Apple-style-span" style="font-size: 12px;"><a rel="nofollow" target="_blank" href="http://www.gnu.org/philosophy/no-word-attachments.html">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div></div><div><font class="yiv2142409382Apple-style-span" color="#0023e9" size="3"><br></font></div></span><br class="yiv2142409382Apple-interchange-newline">
</div>
<br></div></div><br>-----La pièce jointe associée suit-----<br><br><div class="plainMail">_______________________________________________<br>Pw_forum mailing list<br><a ymailto="mailto:Pw_forum@pwscf.org" href="/mc/compose?to=Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></div></blockquote></td></tr></table>