<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Thank you for all response. I have fixed the proplem. epsil=.false.is the right option for metal sysytem.</DIV>

<DIV>And I am very sorry I forget giving my affiliation last time:</DIV>

<DIV>Wang wei,</DIV>

<DIV>Southeast University,</DIV>

<DIV>Nanjing,</DIV>

<DIV>P. R. China<BR>2009.4.29</DIV>

<DIV>--- <B>09年4月26日，周日, pw_forum-request@pwscf.org <I><pw_forum-request@pwscf.org></I></B> 写道：<BR></DIV>

<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><BR>发件人: pw_forum-request@pwscf.org <pw_forum-request@pwscf.org><BR>主题: Pw_forum Digest, Vol 22, Issue 61<BR>收件人: pw_forum@pwscf.org<BR>日期: 2009年4月26日,周日,上午5:54<BR><BR>

<DIV class=plainMail>Send Pw_forum mailing list submissions to<BR>    <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum@pwscf.org" ymailto="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A><BR><BR>To subscribe or unsubscribe via the World Wide Web, visit<BR>    <A href="http://www.democritos.it/mailman/listinfo/pw_forum" target=_blank>http://www.democritos.it/mailman/listinfo/pw_forum</A><BR>or, via email, send a message with subject or body 'help' to<BR>    <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum-request@pwscf.org" ymailto="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</A><BR><BR>You can reach the person managing the list at<BR>    <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum-owner@pwscf.org" ymailto="mailto:pw_forum-owner@pwscf.org">pw_forum-owner@pwscf.org</A><BR><BR>When replying, please edit your Subject line so it is

 more specific<BR>than "Re: Contents of Pw_forum digest..."<BR><BR><BR>Today's Topics:<BR><BR>   1. from phq_readin : error # 1 , no elec. field with metals (? ?)<BR>   2. Re: from phq_readin : error # 1 ,    no elec. field with     metals<BR>      (lan haiping)<BR>   3. Re: from phq_readin : error # 1 ,    no elec. field with metals<BR>      (Eyvaz Isaev)<BR>   4. -nimage -npool -ntg -ndiag (Eduardo Ariel Menendez Proupin)<BR><BR><BR>----------------------------------------------------------------------<BR><BR>Message: 1<BR>Date: Sat, 25 Apr 2009 23:39:28 +0800 (CST)<BR>From: ? ? <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=wangweiphysics@yahoo.com.cn" ymailto="mailto:wangweiphysics@yahoo.com.cn">wangweiphysics@yahoo.com.cn</A>><BR>Subject: [Pw_forum] from phq_readin : error # 1 , no elec. field

with     metals To: <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum@pwscf.org" ymailto="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A> Message-ID: <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=60951.11396.qm@web15208.mail.cnb.yahoo.com" ymailto="mailto:60951.11396.qm@web15208.mail.cnb.yahoo.com">60951.11396.qm@web15208.mail.cnb.yahoo.com</A>> Content-Type: text/plain; charset="utf-8" Dear all,  ?? 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: ?? from phq_readin: error #?????? 1 ?? no elec. field with metals Who can tell me why? Thank you very much! ? The input files as follows: kos2o6.scf.in ??&CONTROL ???????????????? calculation = 'scf' , ???????????????

 restart_mode = 'from_scratch' ,<BR>????????????????????? outdir = '/home/ww/data/kos2o6/' ,<BR>????????????????????? wfcdir = '/home/ww/data/kos2o6/' ,<BR>????????????????? pseudo_dir = '/home/ww/PWscf/pseudo/' ,<BR>????????????????????? prefix = 'kos2o6' ,<BR>?/<BR>?&SYSTEM<BR>?????????????????????? ibrav = 2,<BR>?????????????????? celldm(1) = 19.0945,<BR>???????????????????????? nat = 18,<BR>??????????????????????? ntyp = 3,<BR>???????????????????? ecutwfc = 60 ,<BR>???????????????????? ecutrho = 600 ,<BR>???????????????? occupations = 'smearing' ,<BR>???????????????????? degauss = 0.03<BR>??????????????????? smearing = 'gauss' ,<BR>?/<BR>?&ELECTRONS<BR>??????????????????? conv_thr = 1.0D-8 ,<BR>???????????????? mixing_beta = 0.7 ,<BR>?/<BR>ATOMIC_SPECIES<BR>??? K?? 39.10000? K.pbe-sp.UPF<BR>?? Os? 190.20000? Os.pbe-n-van.UPF<BR>??? O?? 16.00000? O.pbe-van_ak.UPF<BR>ATOMIC_POSITIONS crystal<BR>?? ....<BR>K_POINTS {automatic}<BR>? 10? 10? 10?

 0? 0? 0<BR><BR>kos2o6.phG.in:<BR>???? phonons of kos2o6 at Gamma<BR>?&inputph<BR>? tr2_ph=1.0d-14,<BR>? prefix='kos2o6',<BR>? epsil=.true.,<BR>? amass(1)=39.10,<BR>? amass(2)=190.2,<BR>? amass(3)=16.00,<BR>? outdir='/home/ww/data/kos2o6/',<BR>? fildyn='kos2o6.dynG',<BR>?/<BR>0.0 0.0 0.0<BR><BR><BR><BR>      ___________________________________________________________ <BR>  ????????????????? <BR><A href="http://card.mail.cn.yahoo.com/" target=_blank>http://card.mail.cn.yahoo.com/</A><BR>-------------- next part --------------<BR>An HTML attachment was scrubbed...<BR>URL: <A href="http://www.democritos.it/pipermail/pw_forum/attachments/20090425/0cb9cf45/attachment-0001.htm" target=_blank>http://www.democritos.it/pipermail/pw_forum/attachments/20090425/0cb9cf45/attachment-0001.htm</A> <BR><BR>------------------------------<BR><BR>Message: 2<BR>Date: Sun, 26 Apr 2009 00:56:47 +0800<BR>From: lan haiping <<A

 href="http://cn.mc152.mail.yahoo.com/mc/compose?to=lanhaiping@gmail.com" ymailto="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</A>><BR>Subject: Re: [Pw_forum] from phq_readin : error # 1 ,    no elec. field<BR>    with     metals<BR>To: PWSCF Forum <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum@pwscf.org" ymailto="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A>><BR>Message-ID:<BR>    <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=c92002fa0904250956u2f436576mce31e465a7bc9158@mail.gmail.com" ymailto="mailto:c92002fa0904250956u2f436576mce31e465a7bc9158@mail.gmail.com">c92002fa0904250956u2f436576mce31e465a7bc9158@mail.gmail.com</A>><BR>Content-Type: text/plain; charset="gb2312"<BR><BR>Hi,<BR>According to INPUT_PH,<BR>epsil LOGICAL  *Default:* .false.<BR><BR>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><BR><BR>Therefore, your system should be a metal, and phonon code complained and<BR>stopped.<BR><BR>Please supply your affiliation!<BR>Regards,<BR><BR>On Sat, Apr 25, 2009 at 11:39 PM, ? ? <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=wangweiphysics@yahoo.com.cn" ymailto="mailto:wangweiphysics@yahoo.com.cn">wangweiphysics@yahoo.com.cn</A>> wrote:<BR><BR>> Dear all,<BR>>    I want to calculate phonon frequencies at Gamma for superconductor<BR>> KOs2O6, according to the example02, first I run scf calculation using pw.x<BR>> command without any error, but when I perform the phonon calculation, the<BR>> error appears as forrows:<BR>>    from phq_readin: error #       1<BR>>    no elec. field with metals<BR>> Who can tell me why? Thank

 you very much!<BR>><BR>> The input files as follows:<BR>> kos2o6.scf.in<BR>>   &CONTROL<BR>>                  calculation = 'scf' ,<BR>>                 restart_mode = 'from_scratch' ,<BR>>                       outdir = '/home/ww/data/kos2o6/' ,<BR>>                       wfcdir = '/home/ww/data/kos2o6/' ,<BR>>                   pseudo_dir = '/home/ww/PWscf/pseudo/' ,<BR>>                       prefix = 'kos2o6' ,<BR>>  /<BR>>  &SYSTEM<BR>>                        ibrav =

 2,<BR>>                    celldm(1) = 19.0945,<BR>>                          nat = 18,<BR>>                         ntyp = 3,<BR>>                      ecutwfc = 60 ,<BR>>                      ecutrho = 600 ,<BR>>                  occupations = 'smearing' ,<BR>>                      degauss = 0.03<BR>>                     smearing = 'gauss' ,<BR>>  /<BR>>  &ELECTRONS<BR>>                     conv_thr

 = 1.0D-8 ,<BR>>                  mixing_beta = 0.7 ,<BR>>  /<BR>> ATOMIC_SPECIES<BR>>     K   39.10000  K.pbe-sp.UPF<BR>>    Os  190.20000  Os.pbe-n-van.UPF<BR>>     O   16.00000  O.pbe-van_ak.UPF<BR>> ATOMIC_POSITIONS crystal<BR>>    ....<BR>> K_POINTS {automatic}<BR>>   10  10  10  0  0  0<BR>> kos2o6.phG.in:<BR>>      phonons of kos2o6 at Gamma<BR>> 

 &inputph<BR>>   tr2_ph=1.0d-14,<BR>>   prefix='kos2o6',<BR>>   epsil=.true.,<BR>>   amass(1)=39.10,<BR>>   amass(2)=190.2,<BR>>   amass(3)=16.00,<BR>>   outdir='/home/ww/data/kos2o6/',<BR>>   fildyn='kos2o6.dynG',<BR>>  /<BR>> 0.0 0.0 0.0<BR>><BR>> ------------------------------<BR>> ?????????????????<<A href="http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/" target=_blank>http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/</A>><BR>> _______________________________________________<BR>> Pw_forum mailing list<BR>> <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=Pw_forum@pwscf.org" ymailto="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>><BR><BR><BR>-- <BR>Hai-Ping Lan<BR>Department of Electronics ,<BR>Peking University , Bejing, 100871<BR><A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=lanhaiping@gmail.com" ymailto="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</A>, <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=hplan@pku.edu.cn" ymailto="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</A><BR>-------------- next part --------------<BR>An HTML attachment was scrubbed...<BR>URL: <A href="http://www.democritos.it/pipermail/pw_forum/attachments/20090426/edc33cb1/attachment-0001.htm" target=_blank>http://www.democritos.it/pipermail/pw_forum/attachments/20090426/edc33cb1/attachment-0001.htm</A> <BR><BR>------------------------------<BR><BR>Message: 3<BR>Date: Sat, 25 Apr 2009 09:59:58 -0700 (PDT)<BR>From: Eyvaz Isaev <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=eyvaz_isaev@yahoo.com"

 ymailto="mailto:eyvaz_isaev@yahoo.com">eyvaz_isaev@yahoo.com</A>><BR>Subject: Re: [Pw_forum] from phq_readin : error # 1 ,    no elec. field<BR>    with metals<BR>To: PWSCF Forum <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum@pwscf.org" ymailto="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A>><BR>Message-ID: <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=950057.79544.qm@web65704.mail.ac4.yahoo.com" ymailto="mailto:950057.79544.qm@web65704.mail.ac4.yahoo.com">950057.79544.qm@web65704.mail.ac4.yahoo.com</A>><BR>Content-Type: text/plain; charset=utf-8<BR><BR><BR>Dear unknown user,<BR><BR>First of all, all questions put in this forum should be signed properly, i.e. you should provide your name and affiliation.<BR><BR>As concerns your question it is quite simple, and very basic. Can you please tell us what is macroscopic dielectric constant for metals? If you know the answer, you can fix

 the problem quite easily.<BR><BR>By the way, the message from phq_readin.f90 almost contains the answer what you should do. Read carefully INPUT_PH.txt (in /Doc) to learn more about parameters you included in pn.in file. <BR><BR>Bests,<BR>Eyvaz.<BR><BR>-------------------------------------------------------------------<BR>Prof. Eyvaz Isaev, <BR>Theoretical Physics Department, Moscow State Institute of Steel & Alloys, Russia, <BR>Department of Physics, Chemistry, and Biology (IFM), Linkoping University, Sweden <BR>Condensed Matter Theory Group, Uppsala University, Sweden <BR><A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=Eyvaz.Isaev@fysik.uu.se" ymailto="mailto:Eyvaz.Isaev@fysik.uu.se">Eyvaz.Isaev@fysik.uu.se</A>, <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=isaev@ifm.liu.se" ymailto="mailto:isaev@ifm.liu.se">isaev@ifm.liu.se</A>, <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=eyvaz_isaev@yahoo.com"

 ymailto="mailto:eyvaz_isaev@yahoo.com">eyvaz_isaev@yahoo.com</A><BR><BR><BR>--- On Sat, 4/25/09, ? ? <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=wangweiphysics@yahoo.com.cn" ymailto="mailto:wangweiphysics@yahoo.com.cn">wangweiphysics@yahoo.com.cn</A>> wrote:<BR><BR>> From: ? ? <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=wangweiphysics@yahoo.com.cn" ymailto="mailto:wangweiphysics@yahoo.com.cn">wangweiphysics@yahoo.com.cn</A>><BR>> Subject: [Pw_forum] from phq_readin : error # 1 , no elec. field with metals<BR>> To: <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum@pwscf.org" ymailto="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A><BR>> Date: Saturday, April 25, 2009, 7:39 PM<BR>> Dear all, <BR>> ?? I want to calculate phonon frequencies at Gamma for<BR>> superconductor KOs2O6, according to the example02, first I<BR>> run scf calculation using pw.x command without any error,<BR>>

 but when I perform the phonon calculation, the error appears<BR>> as forrows:<BR>> ?? from phq_readin: error #?????? 1<BR>> ?? no elec. field with metals<BR>> Who can tell me why? Thank you very much!<BR>> ?<BR>> The input files as follows:<BR>> kos2o6.scf.in<BR>> ??&CONTROL<BR>> ???????????????? calculation =<BR>> 'scf' ,<BR>> ??????????????? restart_mode =<BR>> 'from_scratch' ,<BR>> ????????????????????? outdir =<BR>> '/home/ww/data/kos2o6/' ,<BR>> ????????????????????? wfcdir =<BR>> '/home/ww/data/kos2o6/' ,<BR>> ????????????????? pseudo_dir =<BR>> '/home/ww/PWscf/pseudo/' ,<BR>> ????????????????????? prefix =<BR>> 'kos2o6' ,<BR>> ?/<BR>> ?&SYSTEM<BR>> ?????????????????????? ibrav = 2,<BR>> ?????????????????? celldm(1) = 19.0945,<BR>> ???????????????????????? nat = 18,<BR>> ??????????????????????? ntyp = 3,<BR>> ???????????????????? ecutwfc = 60 ,<BR>>

 ???????????????????? ecutrho = 600 ,<BR>> ???????????????? occupations =<BR>> 'smearing' ,<BR>> ???????????????????? degauss = 0.03<BR>> ??????????????????? smearing =<BR>> 'gauss' ,<BR>> ?/<BR>> ?&ELECTRONS<BR>> ??????????????????? conv_thr = 1.0D-8 ,<BR>> ???????????????? mixing_beta = 0.7 ,<BR>> ?/<BR>> ATOMIC_SPECIES<BR>> ??? K?? 39.10000? K.pbe-sp.UPF<BR>> ?? Os? 190.20000? Os.pbe-n-van.UPF<BR>> ??? O?? 16.00000? O.pbe-van_ak.UPF<BR>> ATOMIC_POSITIONS crystal<BR>> ?? ....<BR>> K_POINTS {automatic}<BR>> ? 10? 10? 10? 0? 0? 0<BR>> <BR>> kos2o6.phG.in:<BR>> ???? phonons of kos2o6 at Gamma<BR>> ?&inputph<BR>> ? tr2_ph=1.0d-14,<BR>> ? prefix='kos2o6',<BR>> ? epsil=.true.,<BR>> ? amass(1)=39.10,<BR>> ? amass(2)=190.2,<BR>> ? amass(3)=16.00,<BR>> ? outdir='/home/ww/data/kos2o6/',<BR>> ? fildyn='kos2o6.dynG',<BR>> ?/<BR>> 0.0 0.0 0.0<BR>> <BR>>

 <BR>> <BR>>      <BR>> ___________________________________________________________ <BR>>   ????????????????? <BR>> <A href="http://card.mail.cn.yahoo.com/_______________________________________________" target=_blank>http://card.mail.cn.yahoo.com/_______________________________________________</A><BR>> Pw_forum mailing list<BR>> <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=Pw_forum@pwscf.org" ymailto="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><BR>      <BR><BR><BR>------------------------------<BR><BR>Message: 4<BR>Date: Sat, 25 Apr 2009 17:48:27 -0400<BR>From: Eduardo Ariel Menendez Proupin <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=eariel99@gmail.com"

 ymailto="mailto:eariel99@gmail.com">eariel99@gmail.com</A>><BR>Subject: [Pw_forum] -nimage -npool -ntg -ndiag<BR>To: <A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=pw_forum@pwscf.org" ymailto="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A><BR>Message-ID:<BR>    <<A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=be564dcb0904251448t122bb093g6b9b304dab055263@mail.gmail.com" ymailto="mailto:be564dcb0904251448t122bb093g6b9b304dab055263@mail.gmail.com">be564dcb0904251448t122bb093g6b9b304dab055263@mail.gmail.com</A>><BR>Content-Type: text/plain; charset="iso-8859-1"<BR><BR>Hi, I found in the new, or maybe not so new, users guide, that 1000 atoms or<BR>so can be calculated, and new ways to paralelize.<BR><BR>The example in the manual is<BR><BR>mpirun -np 4096 ./pw.x -nimage 8 -npool 2 -ntg 8 -ndiag 144 -input myinput.in<BR><BR>I have played a bit, but not with a massive computer, and I have found that<BR>the default options

are always better than my unexpert  choices. So, I would like to see some hints, in addition to what  is reproduced below (from the users guide), about the good choices of -ntg and -ndiag. Maybe some examples is enough to understand it. >From the users guide: This execute the PWscf code on 4096 processors, to simulate a system with 8 images, each of which is distributed across 512 processors. K-points are distributed across 2 pools of 256 processors each, 3D FFT is performed using 8 task groups (64 processors each, so the 3D real-space grid is cut into 64 slices), and the diagonalization of the subspace Hamiltonian is distributed to a square grid of 144 processors (12x12). Default values are: -nimage 1 -npool 1 -ntg 1 ; ndiag is chosen by the code as the fastest n^2 (n integer) that fits into the size of each pool. *Massively parallel calculations*: For very large jobs (i.e.

 O(1000) atoms<BR>or so) or for very long jobs to be run on massively parallel machines (e.g.<BR>IBM BlueGene) it is crucial to use in an effective way both the "task group"<BR>and the "ortho group" parallelization. Without a judicious choice of<BR>parameters, large jobs will find a stumbling block in either memory or CPU<BR>requirements. In particular, the "ortho group" parallelization is used in<BR>the diagonalization of matrices in the subspace of Kohn-Sham states (whose<BR>dimension is as a strict minumum equal to the number of occupied states).<BR>These are stored as block-distributed matrixes (distributed across<BR>processors) and diagonalized using custom-taylored diagonalization<BR>algorithms that work on block-distributed matrixes.<BR>Thanks<BR><BR>-- <BR>Eduardo Menendez<BR>Departamento de Fisica<BR>Facultad de Ciencias<BR>Universidad de Chile<BR>Phone: (56)(2)9787439<BR>URL: <A href="http://fisica.ciencias.uchile.cl/~emenendez"

 target=_blank>http://fisica.ciencias.uchile.cl/~emenendez</A><BR>-------------- next part --------------<BR>An HTML attachment was scrubbed...<BR>URL: <A href="http://www.democritos.it/pipermail/pw_forum/attachments/20090425/f69ac51d/attachment.htm" target=_blank>http://www.democritos.it/pipermail/pw_forum/attachments/20090425/f69ac51d/attachment.htm</A> <BR><BR>------------------------------<BR><BR>_______________________________________________<BR>Pw_forum mailing list<BR><A href="http://cn.mc152.mail.yahoo.com/mc/compose?to=Pw_forum@pwscf.org" ymailto="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><BR>End of Pw_forum Digest, Vol 22, Issue 61<BR>****************************************<BR></DIV></BLOCKQUOTE></td></tr></table><br>

      <hr size=1><a href="http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/">  好玩贺卡等你发，邮箱贺卡全新上线！</a>