<DIV>Date: Mon,&nbsp;16 Mar 2009 15:36:14 +0800<BR>From:&nbsp;cdm_mail&lt;<A href="mailto:cdm_mail@163.com">cdm_mail@163.com</A>&gt;<BR>Subject: [Pw_forum] help to look over<BR>To: PWSCF Forum &lt;<A href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</A>&gt;<BR>Message-ID:<BR>&lt;<A href="mailto:chendongming19@gmail.com">chendongming19@gmail.com</A>&gt;<BR>Content-Type: text/plain</DIV>
<DIV>Dear Giovanni,<BR>&nbsp; I'm testing our system with more or less the same configuration.</DIV>
<DIV>&nbsp;I&nbsp;find that the x (y,z) coordinates of most&nbsp;ions&nbsp;are more&nbsp;larger than&nbsp; the value of&nbsp; crystall lattice( orthorhombic P)&nbsp;. This process&nbsp;is after my car-parrinello dynamics.I think it is something wrong with setting the parameter of electronic&nbsp;&nbsp;minimazation&nbsp; process&nbsp;.So&nbsp;the following part is&nbsp;the setting of the begining&nbsp;test that i have reset again.</DIV>
<DIV>&nbsp;&nbsp; It is kind of you to help me for correcting the unsuitable&nbsp; parameters &nbsp;and tell me&nbsp;some key points&nbsp;to&nbsp;complete&nbsp;CPMD work.&nbsp;</DIV>
<DIV>&nbsp;&nbsp; Best wishes.<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &amp;CONTROL<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation='cp'<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix='silica'<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode='from_scratch'<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nstep=100<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iprint=20</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;dt=5.0 </DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ndr=50, ndr=51<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir='./tmp/'<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir='/home/111/Quantum-ESPRESSO/pseudo/'</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;etot_conv_thr=1d-5</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ekin_conv_thr=1d-9&nbsp;&nbsp;&nbsp;&nbsp; <BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &amp;SYSTEM<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav=8,</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;celldm(1) = 18.788, celldm(2) =&nbsp;0.866,&nbsp;&nbsp;celldm(3) = 1.0972<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat=72<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp=2</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbnd=192</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; spin=1</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nbr1=20, nbr2=20, nbr3=20</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &amp;ELECTRONS</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; emass=700.,emass_cutoff=5<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecectron_dynamics='damp'</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecectron_dynamics=damping=0.2<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; startingwfc='random'</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /</DIV>
<DIV>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &amp;IONS<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ion_dynamics='none'<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp; ATOMIC_SPECIES<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp; 16.00&nbsp;&nbsp;&nbsp; O.pz-rrkjus.UPF<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Si&nbsp;&nbsp;&nbsp;28.00&nbsp;&nbsp;&nbsp; Si.vbc.UPF<BR>&gt;&nbsp;&nbsp;&nbsp;&nbsp; ATOMIC_POSITIONS {&nbsp;bohr }<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; 3.5406783&nbsp; 15.0448465&nbsp;&nbsp; 1.4126108<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; 3.8071642&nbsp; 15.4043722&nbsp; 11.9216204<BR>&nbsp;O&nbsp;&nbsp;&nbsp; 13.0894146&nbsp; 15.3948822&nbsp;&nbsp; 1.5781510<BR>O&nbsp;&nbsp;&nbsp; 13.2522097&nbsp; 15.1658554&nbsp; 11.9361830<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 8.2785587&nbsp;&nbsp; 7.2548118&nbsp;&nbsp; 1.5474833<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 8.5399103&nbsp;&nbsp; 7.0021067&nbsp; 11.6012125<BR>O&nbsp;&nbsp;&nbsp; 17.9968777&nbsp;&nbsp; 6.9402800&nbsp;&nbsp; 1.2878237<BR>O&nbsp;&nbsp;&nbsp; 17.6457939&nbsp;&nbsp; 7.1132650&nbsp; 11.9548235<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 2.4405613&nbsp;&nbsp; 6.3799176&nbsp;&nbsp; 4.7482071<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 2.6654053&nbsp;&nbsp; 6.3422036&nbsp; 15.0544825<BR>O&nbsp;&nbsp;&nbsp; 12.2024078&nbsp;&nbsp; 6.3211975&nbsp;&nbsp; 4.7142868<BR>O&nbsp;&nbsp;&nbsp; 11.8810415&nbsp;&nbsp; 6.4080858&nbsp; 15.3200722<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 7.1940126&nbsp; 14.4133511&nbsp;&nbsp; 4.9239106<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 7.3077059&nbsp; 14.4793634&nbsp; 15.0265102<BR>O&nbsp;&nbsp;&nbsp; 16.9006786&nbsp; 14.3864670&nbsp;&nbsp; 5.0179539<BR>O&nbsp;&nbsp;&nbsp; 16.6793461&nbsp; 14.4224520&nbsp; 15.3299093<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 4.7042437&nbsp; 11.5882969&nbsp;&nbsp; 8.4145012<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 4.8058362&nbsp; 11.7974415&nbsp; 18.5207825<BR>O&nbsp;&nbsp;&nbsp; 14.0459738&nbsp; 11.6190434&nbsp;&nbsp; 8.4245062<BR>O&nbsp;&nbsp;&nbsp; 14.1501627&nbsp; 11.7520666&nbsp; 18.6828041<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 9.4120684&nbsp;&nbsp; 3.5592277&nbsp;&nbsp; 8.1602020<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 9.4435863&nbsp;&nbsp; 3.4508531&nbsp; 18.8181076<BR>O&nbsp;&nbsp;&nbsp; 18.8548489&nbsp;&nbsp; 3.3799570&nbsp;&nbsp; 8.2510424<BR>O&nbsp;&nbsp;&nbsp; 18.8457680&nbsp;&nbsp; 3.6430509&nbsp; 18.4409237<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 4.6323986&nbsp;&nbsp; 4.9528089&nbsp;&nbsp; 9.3928719<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 4.6253457&nbsp;&nbsp; 5.0198221&nbsp; 19.3735943<BR>O&nbsp;&nbsp;&nbsp; 14.1270161&nbsp;&nbsp; 5.0775709&nbsp;&nbsp; 9.1806593<BR>O&nbsp;&nbsp;&nbsp; 14.2851896&nbsp;&nbsp; 5.0468230&nbsp; 19.6615734<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 9.3298197&nbsp; 13.1670465&nbsp;&nbsp; 9.3121605<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 9.6129112&nbsp; 13.2867460&nbsp; 19.5585632<BR>O&nbsp;&nbsp;&nbsp; 18.7938595&nbsp; 13.0344448&nbsp;&nbsp; 9.2933855<BR>O&nbsp;&nbsp;&nbsp; 18.6957073&nbsp; 13.2711039&nbsp; 19.5695648<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 3.7519393&nbsp;&nbsp; 1.4702886&nbsp;&nbsp; 5.7478361<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 3.7406268&nbsp;&nbsp; 1.3143321&nbsp; 16.1506748<BR>O&nbsp;&nbsp;&nbsp; 13.1751394&nbsp;&nbsp; 1.4241230&nbsp;&nbsp; 5.8916039<BR>O&nbsp;&nbsp;&nbsp; 13.2039194&nbsp;&nbsp; 1.6296371&nbsp; 16.2666016<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 8.5475645&nbsp;&nbsp; 9.7668839&nbsp;&nbsp; 5.8549895<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 8.3606243&nbsp;&nbsp; 9.4664469&nbsp; 16.0193081<BR>O&nbsp;&nbsp;&nbsp; 18.0083542&nbsp;&nbsp; 9.4706783&nbsp;&nbsp; 5.8703651<BR>O&nbsp;&nbsp;&nbsp; 17.8562260&nbsp;&nbsp; 9.4727955&nbsp; 15.9772587<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 2.7527096&nbsp; 10.6110735&nbsp;&nbsp; 2.3466952<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 2.6490324&nbsp; 10.4995565&nbsp; 12.5345230<BR>O&nbsp;&nbsp;&nbsp; 11.8460188&nbsp; 10.5480690&nbsp;&nbsp; 2.2722039<BR>O&nbsp;&nbsp;&nbsp; 11.9239492&nbsp; 10.4055338&nbsp; 12.5277109<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 7.4533672&nbsp;&nbsp; 2.5166929&nbsp;&nbsp; 2.2388670<BR>O&nbsp;&nbsp;&nbsp;&nbsp; 7.3775740&nbsp;&nbsp; 2.1962755&nbsp; 12.8475285<BR>O&nbsp;&nbsp;&nbsp; 16.7454376&nbsp;&nbsp; 2.4293642&nbsp;&nbsp; 2.5411322<BR>O&nbsp;&nbsp;&nbsp; 16.6413174&nbsp;&nbsp; 2.5064237&nbsp; 12.7736511<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 0.4496534&nbsp;&nbsp; 8.2138128&nbsp;&nbsp; 3.6896179<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 0.5399753&nbsp;&nbsp; 8.1889200&nbsp; 13.9777460<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 9.8792067&nbsp;&nbsp; 8.2886143&nbsp;&nbsp; 3.7623391<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 9.8359041&nbsp;&nbsp; 8.5231533&nbsp; 14.1156855<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 5.2412066&nbsp;&nbsp; 0.0271945&nbsp;&nbsp; 3.5170169<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 5.1725807&nbsp;&nbsp; 0.0363324&nbsp; 13.8097649<BR>Si&nbsp;&nbsp;&nbsp; 14.5211735&nbsp;&nbsp; 0.0875002&nbsp;&nbsp; 3.5420916<BR>Si&nbsp;&nbsp;&nbsp; 14.3368998&nbsp;&nbsp; 0.0028480&nbsp; 14.0652027<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 2.4036527&nbsp; 12.8000431&nbsp;&nbsp; 0.2863795<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 2.3965807&nbsp; 12.7267065&nbsp; 10.3207769<BR>Si&nbsp;&nbsp;&nbsp; 11.6870117&nbsp; 12.6637812&nbsp;&nbsp; 0.1414142<BR>Si&nbsp;&nbsp;&nbsp; 11.6615248&nbsp; 12.7831507&nbsp; 10.3878822<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 7.2631865&nbsp;&nbsp; 4.6357923&nbsp;&nbsp; 0.2741336<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 6.9719872&nbsp;&nbsp; 4.4316664&nbsp; 10.6666307<BR>Si&nbsp;&nbsp;&nbsp; 16.4332924&nbsp;&nbsp; 4.5068045&nbsp;&nbsp; 0.2921807<BR>Si&nbsp;&nbsp;&nbsp; 16.6266994&nbsp;&nbsp; 4.4123287&nbsp; 10.3654356<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 2.5121100&nbsp;&nbsp; 4.1013050&nbsp;&nbsp; 7.1919127<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 2.2047310&nbsp;&nbsp; 4.1441984&nbsp; 17.4181461<BR>Si&nbsp;&nbsp;&nbsp; 11.9457951&nbsp;&nbsp; 3.9283929&nbsp;&nbsp; 7.2228098<BR>Si&nbsp;&nbsp;&nbsp; 11.7723646&nbsp;&nbsp; 3.9849200&nbsp; 17.5345020<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 7.1349611&nbsp; 12.1880169&nbsp;&nbsp; 7.1159267<BR>Si&nbsp;&nbsp;&nbsp;&nbsp; 7.1111145&nbsp; 12.1099644&nbsp; 17.5283604<BR>Si&nbsp;&nbsp;&nbsp; 16.3560925&nbsp; 12.0432882&nbsp;&nbsp; 6.9771399<BR>Si&nbsp;&nbsp;&nbsp; 16.5531750&nbsp; 12.0766220&nbsp; 17.5176411<BR><BR><BR><BR><BR></DIV><br><!-- footer --><br><span title="neteasefooter"/><hr/>
<a href="http://www.yeah.net">网易邮箱，中国第一大电子邮件服务商</a>
</span>