According to my experience, phonon calculation is very sensitive to the groud state. Therefore, a strict relax (by pw.x) is required before using ph.x, for example, higher cut-off, more k-points and smaller threshold.<br><div>--<br>GAO Zhe<br>CMC Lab, MSE, SNU, Seoul, S.Korea<br>        
</div><br>At 2011-05-26 17:15:53£¬yhho&nbsp;&lt;yhho@pub.iams.sinica.edu.tw&gt; wrote:<br> <blockquote id="isReplyContent" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">

<div><font size="2">Deae all,</font></div>
<div><font size="2">I'm trying to do the <font color="#ff0000">phonon 
dispersion</font> calaulation&nbsp;of Graphene.</font></div>
<div><font size="2">Although the jobs&nbsp;have done sucessfully, the result are 
not great.</font></div>
<div><font size="2">I think it may be&nbsp;coused by the wrong setting of my <font color="#ff0000">q2r.in</font> or <font color="#ff0000">matdyn.in</font>.</font></div>
<div><font size="2">But some of terms in thoese input file are unfamiliar to 
me.</font></div>
<div><font size="2">Please give me some suggestions&nbsp;about setting them. 
Thanks.</font></div>
<div><font size="2"></font>&nbsp;</div>
<div>
<div><font color="#ff0000" size="2">&lt;ph.in&gt;</font></div>
<div><font size="2"><font color="#0000ff">phonons of 
AuC<br>&nbsp;&amp;inputph<br>&nbsp; tr2_ph=1.0d-14,<br>&nbsp; 
ldisp=.true.,<br>&nbsp; nq1=4, nq2=4, nq3=4<br>&nbsp; amass(1)= 
12.0107<br>&nbsp; prefix='AuC',<br>&nbsp; outdir='./'<br>&nbsp; 
fildyn='AuC.dyn',<br>&nbsp;/<br></font></font></div></div>
<div><font color="#ff0000" size="2">&lt;q2r.in&gt;</font></div>
<div><font color="#0000ff" size="2">&amp;input fildyn='AuC.dyn', zasr='crystal', 
flfrc='AuC444.fc' /</font></div>
<div><font size="2"></font>&nbsp;</div>
<div><font color="#ff0000" size="2">&lt;matdyn.in&gt;</font></div>
<div><font color="#0000ff" size="2">&nbsp;&amp;input<br>&nbsp;&nbsp;&nbsp; 
asr='simple',&nbsp; amass(1)=12.0107,<br>&nbsp;&nbsp;&nbsp; flfrc='AuC444.fc', 
flfrq='AuC.freq'<br>&nbsp;/<br>&nbsp;&nbsp; 13<br>&nbsp; 0.00000000&nbsp; 
0.00000000&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 1<br>&nbsp; 0.00000000&nbsp; 
0.12500000&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 2<br>&nbsp; 0.00000000&nbsp; 
0.25000000&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 3<br>&nbsp; 0.00000000&nbsp; 
0.37500000&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 4<br>&nbsp; 0.00000000&nbsp; 
0.50000000&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 5<br>&nbsp; 0.08333250&nbsp; 
0.45833250&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 6<br>&nbsp; 0.16666500&nbsp; 
0.41666500&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 7<br>&nbsp; 0.24999750&nbsp; 
0.37499750&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 8<br>&nbsp; 0.33333000&nbsp; 
0.33333000&nbsp; 0.00000000&nbsp;&nbsp;&nbsp; 9<br>&nbsp; 0.24999750&nbsp; 
0.24999750&nbsp; 0.00000000&nbsp;&nbsp; 10<br>&nbsp; 0.16666500&nbsp; 
0.16666500&nbsp; 0.00000000&nbsp;&nbsp; 11<br>&nbsp; 0.08333250&nbsp; 
0.08333250&nbsp; 0.00000000&nbsp;&nbsp; 12<br>&nbsp; 0.00000000&nbsp; 
0.00000000&nbsp; 0.00000000&nbsp;&nbsp; 13<br></font></div>
<div><font size="2">Otherwise, I have another question about the last part of 
matdyn.in.</font></div>
<div><font size="2">Does that mean the points such as gamma, M, K 
,gamma?</font></div>
<div><font size="2">I'm not sure about that...</font></div>
<div><font size="2">Thanks a lot.</font></div>
<div><font size="2"></font>&nbsp;</div>
<div><font size="2">Best regarts,</font></div>
<div><font size="2"></font>&nbsp;</div>
<div><font size="2">Yu Han Ho</font></div>
<div><font size="2"><font color="#0000ff"></font>&nbsp;</font></div>
</blockquote><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>