<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><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 cellspacing="0" cellpadding="0" border="0" style="position: relative; z-index: 0; "><tbody><tr><td valign="top" style="font: inherit; ">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 cellspacing="0" cellpadding="0" border="0" style="position: relative; z-index: 0; "><tbody><tr><td valign="top" style="font: inherit; ">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 cellspacing="0" cellpadding="0" border="0" style="position: relative; z-index: 0; "><tbody><tr><td valign="top" style="font: inherit; "><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 cellspacing="0" cellpadding="0" border="0" style="position: static; z-index: auto; "><tbody><tr><td valign="top" style="font: inherit; ">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 cellspacing="0" cellpadding="0" border="0" style="position: static; z-index: auto; "><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="Apple-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-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: 0px; font-size: medium; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-size: 14px; "><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-size: 14px; "><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-size: 14px; "><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-size: 14px; "><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-size: 14px; "><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><div><font class="Apple-style-span" color="#0023E9" size="3"><br></font></div></span><br class="Apple-interchange-newline">
</div>
<br></div></body></html>