Hi,<br>According to INPUT_PH, <br><table style="border: 2px solid rgb(181, 181, 0); margin-bottom: 10px; table-layout: auto; background-color: rgb(255, 255, 255);" width="100%"><tbody><tr><th style="padding: 2px 2px 2px 10px; background: rgb(255, 255, 153) none repeat scroll 0% 0%; -moz-background-clip: -moz-initial; -moz-background-origin: -moz-initial; -moz-background-inline-policy: -moz-initial;" align="left" valign="top" width="20%">
epsil</th>
<td style="padding: 2px 2px 2px 5px; background: rgb(255, 255, 195) none repeat scroll 0% 0%; text-align: left; vertical-align: top; -moz-background-clip: -moz-initial; -moz-background-origin: -moz-initial; -moz-background-inline-policy: -moz-initial;">
LOGICAL</td>
</tr>
<tr>
<td style="padding: 2px 10px; background: rgb(255, 255, 195) none repeat scroll 0% 0%; text-align: right; vertical-align: top; -moz-background-clip: -moz-initial; -moz-background-origin: -moz-initial; -moz-background-inline-policy: -moz-initial;">
<i>Default:</i></td>
<td style="padding: 2px 2px 2px 5px; background: rgb(255, 243, 217) none repeat scroll 0% 0%; text-align: left; vertical-align: top; -moz-background-clip: -moz-initial; -moz-background-origin: -moz-initial; -moz-background-inline-policy: -moz-initial;">
 .false.
         </td>
</tr>
<tr><td colspan="2" align="left" valign="top"><blockquote><pre>If .true. in a q=0 calculation for a non metal the<br>macroscopic dielectric constant of the system is<br>computed. Do not set epsil to .true. if you have a<br>
metallic system or q/=0: the code will complain and stop.<br>         </pre></blockquote></td></tr></tbody></table><br>Therefore, your system should be a metal, and phonon code complained and stopped.<br><br>Please supply your affiliation!<br>
Regards,<br><br><div class="gmail_quote">On Sat, Apr 25, 2009 at 11:39 PM, 玮 王 <span dir="ltr">&lt;<a href="mailto:wangweiphysics@yahoo.com.cn">wangweiphysics@yahoo.com.cn</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit;" valign="top">
<div>Dear all, </div>
<div>&nbsp;&nbsp; I want to calculate phonon frequencies at Gamma for superconductor KOs2O6, according to the example02, first I run scf calculation using pw.x command without any error, but when I perform the phonon calculation, the error appears as forrows:</div>

<div>&nbsp;&nbsp; from phq_readin: error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1</div>
<div>&nbsp;&nbsp; no elec. field with metals</div>
<div>Who can tell me why? Thank you very much!</div>
<div>&nbsp;</div>
<div>The input files as follows:</div>
<div><a href="http://kos2o6.scf.in" target="_blank">kos2o6.scf.in</a></div>
<div>&nbsp;&nbsp;&amp;CONTROL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = &#39;scf&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = &#39;from_scratch&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = &#39;/home/ww/data/kos2o6/&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; wfcdir = &#39;/home/ww/data/kos2o6/&#39; ,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;/home/ww/PWscf/pseudo/&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = &#39;kos2o6&#39; ,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 19.0945,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 18,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 3,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 60 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 600
 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; occupations = &#39;smearing&#39; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; degauss = 0.03<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; smearing = &#39;gauss&#39; ,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 1.0D-8 ,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.7 ,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp; 39.10000&nbsp; K.pbe-sp.UPF<br>&nbsp;&nbsp; Os&nbsp; 190.20000&nbsp; Os.pbe-n-van.UPF<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 16.00000&nbsp; O.pbe-van_ak.UPF<br>ATOMIC_POSITIONS crystal<br>&nbsp;&nbsp;
 ....</div>
<div>K_POINTS {automatic}<br>&nbsp; 10&nbsp; 10&nbsp; 10&nbsp; 0&nbsp; 0&nbsp; 0<br></div>
<div><a href="http://kos2o6.phG.in" target="_blank">kos2o6.phG.in</a>:</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp; phonons of kos2o6 at Gamma<br>&nbsp;&amp;inputph<br>&nbsp; tr2_ph=1.0d-14,<br>&nbsp; prefix=&#39;kos2o6&#39;,<br>&nbsp; epsil=.true.,<br>&nbsp; amass(1)=39.10,<br>&nbsp; amass(2)=190.2,<br>&nbsp; amass(3)=16.00,<br>&nbsp; outdir=&#39;/home/ww/data/kos2o6/&#39;,<br>
&nbsp; fildyn=&#39;kos2o6.dynG&#39;,<br>&nbsp;/<br>0.0 0.0 0.0<br></div></td></tr></tbody></table><br>


      <hr size="1"><a href="http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/" target="_blank">  好玩贺卡等你发，邮箱贺卡全新上线！</a><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto: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>
<br></blockquote></div><br><br clear="all"><br>-- <br>Hai-Ping Lan <br>Department of Electronics ,<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a><br>