Dear PWscf users,<br>I did some phonon calculations using pw.x and ph.x and obtained the results of phonon density of states (1.phdos). But when I calculated the free energy using fqha.x, some wrong informations always appear.<br>In the fqha.f90 file, I deleted the following line "if (abs(nu(i)-de*(i-0.5d0)).gt.1.0d-6) stop ' wrong grid'" (I think that the limitative condition is too rigorous).<br><br>the executing processing:<br>fqha<br> file containing the dos >>>1.phdos<br>1368<br>0.2048<br>1.0<br> file for the Free energy >>>1.txt<br> temperature >>>500.0<br><br>the error information:<br>"run-time error M6202: MATH<br>- log: SING error"<br><br>I also know the negative value appears in log function described as the following line,<br>"    if (T.gt.0.d0) Ftot=Ftot+dos(i)*a2*T*log(1.0d0-exp(-a3*nu(i)/T))<br>"<br>but I do not know how to resolve this problem.<br>Who can help me?<br>Thanks in advance.<br><br>Best regards,<br>Chaohao<br><br><!-- footer --><br><br><br>
<font style="font-size:14.8px">
<!--广告footer 开始-->
<!--广告footer 结束-->

<!--内部footer开始-->
<br>
<br><p style="line-height:160%;">
<a href="http://mail.163.com" target="_blank">
<font color=blue>需要一个<span style="color:red">2000兆的免费</span>邮箱吗？<br>网易免费邮箱是<span style="color:red">中国最多人使用</span>的电子邮箱。</a>
</p>
<!--内部footer结束-->
</font>