<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:10pt">scf DO affect the result rerun for scf<br><div>&nbsp;</div><font style="font-family: times new roman,new york,times,serif;" size="2">-- <br></font><font style="font-family: times new roman,new york,times,serif;" size="2"><span style="color: rgb(96, 154, 159); font-style: italic;">Rudra<br></span><span style="color: rgb(127, 63, 0);">JRF; SNBNCBS<br></span></font><font style="font-family: times new roman,new york,times,serif;" size="2"><a rel="nofollow" target="_blank" href="http://www.bose.res.in/%7Erudra%20%20%20%20%20%20%20"><span style="color: rgb(0, 0, 255); font-style: italic;">http://www.bose.res.in/~rudra</span></a><br></font><font style="font-family: times new roman,new york,times,serif;" size="2"><span style="color: rgb(0, 127, 64); font-style: italic;"><br>A bus station is where a bus

 stops. A train station is where a train</span><br style="color: rgb(0, 127, 64); font-style: italic;"><span style="color: rgb(0, 127, 64); font-style: italic;">stops. On my desk I have a work&nbsp; station.</span><span style="font-style: italic;"> </span><br style="font-style: italic;"><span style="color: rgb(255, 0, 0); font-style: italic;"><br>Please, if possible, don't  send me MS Word or PowerPoint attachments <br>Why?See:&nbsp;  </span><span style="font-style: italic;"><a rel="nofollow" target="_blank" href="http://www.gnu.org/philosophy/no-word-attachments.html%20%20%20%20%20%20%20%20%20%20%20">http://www.gnu.org/philosophy/no-word-attachments.html</a><br></span></font><div style="font-family: times new roman,new york,times,serif; font-size: 10pt;"><br><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">----- Original Message ----<br>From: herofather &lt;herofather@163.com&gt;<br>To: pw_forum@pwscf.org<br>Sent:

 Monday, 31 March, 2008 2:31:36 PM<br>Subject: [Pw_forum] I cannot get my .dyn(m-n) file<br><br><div>&nbsp;hello everybody,</div>

<div>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; i want to get the dynamics matrix file of the molecule H2O. so i first do a scf calculation and then a ph.x one. but i cannot get this file</div>

<div>&nbsp;&amp;CONTROL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; title = H2O ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; calculation = 'scf' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; wf_collect = .true. ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/home/pwscf/tmp/' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pseudo_dir = '/home/pwscf/espresso-3.2.3/pseudo/'

, &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = H2O , &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; disk_io = 'default' , &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; verbosity = 'default' , &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tstress = .false. , &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tprnfor = .true. , &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lberry = .false.

 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; gdir = 1 ,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 1,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) = 10.0,<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 = 3,<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 = 2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 25

 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 100 ,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; electron_maxstep = 100,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_mode = 'local-TF' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.8 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; diagonalization = 'cg' ,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 16.00000&nbsp; H_US.van <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp; 1.00000&nbsp; O_US.van <br>ATOMIC_POSITIONS angstrom <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; -0.364400000&nbsp;&nbsp;&nbsp; 0.906200000&nbsp;&nbsp;&nbsp; 0000000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;

 H&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.081400000&nbsp;&nbsp;&nbsp; 1.756500000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; H&nbsp;&nbsp;&nbsp;&nbsp; -0056900000&nbsp;&nbsp; -0.003200000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; <br>K_POINTS automatic <br>&nbsp; 4 4 4&nbsp;&nbsp; 0 0 0 </div>

<div>.............................................</div>

<div>the result is </div>

<div>iteration #100&nbsp;&nbsp;&nbsp;&nbsp; ecut=&nbsp;&nbsp;&nbsp; 25.00 Ry&nbsp;&nbsp;&nbsp;&nbsp; beta=0.80<br>&nbsp;&nbsp;&nbsp;&nbsp; CG style diagonalization<br>&nbsp;&nbsp;&nbsp;&nbsp; ethr =&nbsp; 1.11E-07,&nbsp; avg # of iterations =&nbsp; 3.0</div>

<div>&nbsp;&nbsp;&nbsp;&nbsp; total cpu time spent up to now is&nbsp;&nbsp;&nbsp; 382.55 secs</div>

<div>&nbsp;&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -63.94952377 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; =&nbsp;&nbsp; -63.94952369 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp; 0.00003432 Ry</div>

<div>&nbsp;&nbsp;&nbsp;&nbsp; End of self-consistent calculation</div>

<div>&nbsp;&nbsp;&nbsp;&nbsp; convergence NOT achieved, stopping</div>

<div>..............................................</div>

<div>Then i do a ph.x calculation</div>

<div>H2O&nbsp;&nbsp;&nbsp; <br>&nbsp;&amp;INPUTPH<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; outdir = '/home/pwscf/tmp/' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; prefix = H2O ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fildyn = 'H2O.dyn' ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; niter_ph=50,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; epsil = .false.,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; trans = .true.,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ldisp =

 .true.,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; elph = .false.,<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; zue = .true.,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nq1=4,nq2=4,nq3=4,<br>&nbsp;/</div>

<div>..............................................</div>

<div>but as a result i cannot get the H2O.dyn(?) file</div>

<div>I don't know whether the scf convergence will affect the result so much?and what is the reason i cann't get the file. any reply would be greatly appreciated.</div>

<div>&nbsp;</div>

<div>&nbsp;</div>

<div>&nbsp;</div>

<div>&nbsp;</div><br><br><SPAN style="width:100%;height:2px;border-bottom:1px solid rgb(212,208,200); border-top:1px solid rgb(128,128,128);background-color:black;overflow:hidden; margin:8px 0px;"></SPAN>

<a rel="nofollow" style="font-size: 14px; line-height: 15px; color: rgb(0, 0, 0); text-decoration: none;" target="_blank" href="http://event.mail.163.com/chanel/click.htm?from=NO_26&amp;domain=163"><span style="text-decoration: underline; color: blue;">中 国 最 强 网 游 --- 网 易 梦 幻 西 游 ，166 万 玩 家 同 时 在 线</span> </a></div><br></div></div><br>

      <!--10--><hr size=1></hr> Chat on a cool, new interface. No download required. <a href="http://in.rd.yahoo.com/tagline_webmessenger_10/*http://in.messenger.yahoo.com/webmessengerpromo.php">Click here.</a></body></html>