<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title></title>
</head>
<style type="text/css">
<!--
#send_mail {
font-family: "Courier New";
background:#FFFFFF;
color:#000000;
}
#send_mail a {
color:#999999
}
#send_mail a:hover {
color:#0000FF
}
#send_mail div {
color:#000000;
font-family: "Courier New";
font-size:14px;
line-height:150%;
}
#send_mail .usersign {
line-height:100%;
}
-->
</style>
<body id="send_mail">
<div>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<title></title>
<style type="text/css">
<!--
#send_mail {
font-family: "Courier New";
background:#FFFFFF;
color:#000000;
}
#send_mail a {
color:#999999
}
#send_mail a:hover {
color:#0000FF
}
#send_mail div {
color:#000000;
font-family: "Courier New";
font-size:14px;
line-height:150%;
}
#send_mail .usersign {
line-height:100%;
}
-->
</style>
<div><p>Dear QEers:</p><p> When I was calculating vibrational frequencies for a system with nat >= 100, I noticed a small bug of QE. The calculation steps were:</p><p>(1) pw.x for scf electronic calculation;</p><p>(2) ph.x for linear-response calculation;</p><p>(3) q2r.x for the reciprocal-to-real transformation;</p><p>Then, QE claimed that:</p><p>######################</p><p>from read_f: error # 1</p><p>wrong na read <br></p><p>######################</p><p>the error message may be different at some computers, like:</p><p>#################</p><p>end of file during read ...... <br></p><p>################# </p><p>Although this error information has been reported many times on the PW_FORUM, this time, it is caused by "nat >=100".</p><p>(1) The cause <br></p><p> When nat >=100, some lines in *.dyn[123...] file (which includes the reciprocal dynamical matrix) will be like that:</p><p>**************************************</p><p> 34 99<br> 0.00000113 0.00000000 -0.00000895 0.00000000 0.00000474 0.00000000<br> -0.00002272 0.00000000 -0.00002000 0.00000000 -0.00001948 0.00000000<br> 0.00002366 0.00000000 -0.00002707 0.00000000 0.00006012 0.00000000<br> 34100<br> 0.00006484 0.00000000 -0.00000418 0.00000000 -0.00003671 0.00000000<br> -0.00000097 0.00000000 -0.00000993 0.00000000 0.00001364 0.00000000<br> 0.00003728 0.00000000 -0.00000805 0.00000000 0.00002762 0.00000000 <br></p><p>************************************** </p><p>the sencond pair of atomic indices above get close to each other. When execute q2r.x, these two numbers are read together as one number, which is the cause of the error.</p><p> </p><p>(2) The solution</p><p>vi espresso/PH/q2r.f90</p><p>go to the Line 586, which is:</p><p>****************************</p><p> READ(1,*) i,j <br>**************************** </p><p>change it to be</p><p>****************************</p><p> READ(1,'(2i3)') i,j <br>
</p><p>**************************** </p><p>then, compile QE again. </p><p> </p><p> May this could be of a little help to you.</p><p> </p><p>Best Wishes!</p><p>Yours Sincerely</p><p>L. F. Huang <br></p><p class="usersign" align="left">------<br><br>======================================================================<br>
L.F.Huang(黄良锋) DFT and phonon physics<br>
======================================================================<br>
Add: Research Laboratory for Computational Materials Sciences,<br>
Instutue of Solid State Physics,the Chinese Academy of Sciences,<br>
P.O.Box 1129, Hefei 230031, P.R.China<br>
Tel: 86-551-5591464-326(office)<br>
Fax: 86-551-5591434<br>
Our group: http://theory.issp.ac.cn<br>
======================================================================<br>
</p></div>
</div>
</body>
</html>