Dear Sathya,<br><br>First of all I would like you to note that the method used in the PRL is different from that used in QE, while the former is based on LAPW method, Pwscf is a plane wave code. So the structural parameters given in the PRL is a good guess to start your calculation. However, you should first obtain the optimized values of cell parameters and atomic positions using pwscf and then proceed to do other calculations. These values should be comparable to those mentioned in the PRL, provided you use same level of theory (eg. same exchange correlation, same pseudopotential, same k-points, etc.) and you do the things correctly.<br>
<br>As mentioned by others, you should test your pseudopotentials and also perform convergence tests for k-point mesh used, ecutwfc, etc. Moreover, since you are using ultrasoft pseudo potential, I would suggest you to use ecutrho=8-10*ecutwfc, presently u r using ecutrho=4.0*ecutwfc, which is often not sufficient for US psps. At least&nbsp; for structural relaxations,&nbsp; &quot;conv_thr=1e-08&quot; is sufficient, although you may need to use a more strict condition when you calculate the effective charges.<br>
<br>Hope this helps.<br><br>Prasenjit.<br><br>PRASENJIT GHOSH,<br>POST-DOC,<br>ROOM NO: 265, MAIN BUILDING,<br>CM SECTION, ICTP,<br>STRADA COSTERIA 11,<br>TRIESTE, 34104,<br>ITALY<br>PHONE: +39 040 2240 369 (O)<br> +39 3807528672 (M)
<br><br><div class="gmail_quote">2008/7/11 sathya subramanyan &lt;<a href="mailto:sathya_sheela1985001@rediffmail.com">sathya_sheela1985001@rediffmail.com</a>&gt;:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<p>
&nbsp; <br>
&nbsp; <br>
&nbsp; <br>
Dear Prof. Eyvaz,<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  Thank you very much for your reply In response to your advice regarding the inaccuracy in structure input for tetragonal PbTiO3, I checked the literature and thus have taken&nbsp; the structural input parameteres(lattice constant,c/a and atomic positions)from PRL 80,1998,4321 and done the convergence test. conv_threshold used now is 1e-12<br>

&nbsp; The result is as follows.Could you please let me know whether the convergence is good or not.<br>
<br>
THE INPUT FILE FOR TETRAGONAL PbTiO3<br>
<br>
&amp;control<br>
&nbsp; &nbsp; calculation&nbsp; = &#39;scf&#39;<br>
&nbsp; &nbsp; restart_mode = &#39;from_scratch&#39;<br>
&nbsp; &nbsp; pseudo_dir&nbsp;  = &#39;/home/nbaski/espresso-4.0.1/pseudo/&#39;<br>
&nbsp; &nbsp; outdir&nbsp; &nbsp; &nbsp;  = &#39;/home/nbaski/tmp/&#39;<br>
 /<br>
 &amp;system<br>
&nbsp; &nbsp; ibrav=6<br>
&nbsp; &nbsp; celldm(1)=7.373,<br>
&nbsp; &nbsp; celldm(3)=1.065<br>
&nbsp; &nbsp; nat=5<br>
&nbsp; &nbsp; ntyp=3<br>
&nbsp; &nbsp; nbnd=25<br>
&nbsp; &nbsp; ecutwfc=30.0<br>
&nbsp; &nbsp; occupations = &#39;fixed&#39;<br>
&nbsp; &nbsp; degauss=0.00<br>
 /<br>
 &amp;electrons<br>
&nbsp; &nbsp; conv_thr = 1e-12,<br>
&nbsp; &nbsp; mixing_beta=0.3,<br>
 /<br>
ATOMIC_SPECIES<br>
&nbsp; Pb&nbsp;  207.2&nbsp; &nbsp; &nbsp; Pb.vdb.UPF<br>
&nbsp; Ti&nbsp; &nbsp; 47.867&nbsp; &nbsp; Ti.vdb.UPF<br>
&nbsp; O&nbsp; &nbsp;  15.9994&nbsp;  O.vdb.UPF<br>
ATOMIC_POSITIONS<br>
&nbsp; Pb&nbsp; &nbsp; 0.000&nbsp; &nbsp; 0.000&nbsp; &nbsp; 0.000<br>
&nbsp; Ti&nbsp; &nbsp; 0.500&nbsp; &nbsp; 0.500&nbsp; &nbsp; 0.538<br>
&nbsp; O&nbsp; &nbsp;  0.500&nbsp; &nbsp; 0.000&nbsp; &nbsp; 0.612<br>
&nbsp; O&nbsp; &nbsp;  0.000&nbsp; &nbsp; 0.500&nbsp; &nbsp; 0.612<br>
&nbsp; O&nbsp; &nbsp;  0.500&nbsp; &nbsp; 0.500&nbsp; &nbsp; 0.117<br>
K_POINTS {automatic}<br>
&nbsp; 4 4 4 1 1 1<br>
<br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; **************<br>
THE OUTPUT FILE(a part only as it is a heavy file) <br>
<br>
 highest occupied, lowest unoccupied level (ev):&nbsp; &nbsp;  9.6746&nbsp;  12.5554<br>
<br>
!&nbsp; &nbsp; total energy&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; -333.71409381 Ry<br>
&nbsp; &nbsp;  Harris-Foulkes estimate&nbsp;  =&nbsp; -333.71409381 Ry<br>
&nbsp; &nbsp;  estimated scf accuracy&nbsp; &nbsp; &lt;&nbsp; &nbsp; &nbsp; &nbsp; 5.8E-13 Ry<br>
<br>
&nbsp; &nbsp;  The total energy is the sum of the following terms:<br>
<br>
&nbsp; &nbsp;  one-electron contribution =&nbsp;  -94.69940998 Ry<br>
&nbsp; &nbsp;  hartree contribution&nbsp; &nbsp; &nbsp; =&nbsp; &nbsp; 73.25826560 Ry<br>
&nbsp; &nbsp;  xc contribution&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  =&nbsp;  -49.63383021 Ry<br>
&nbsp; &nbsp;  ewald contribution&nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; -262.63911922 Ry<br>
<br>
&nbsp; &nbsp;  convergence has been achieved in&nbsp; 19 iterations<br>
<br>
&nbsp; &nbsp; 
</p>


SATHYA SHEELA.S
<br>
Grad. Student
<br>
Department of&nbsp;Physics 
<br>
National Institute&nbsp;of Technology
<br>
Tiruchirapalli -&nbsp;620015 
<br>
India<br><br>
<table style="font-family: Verdana; font-size: 11px; line-height: 15px;" border="0" cellpadding="0" cellspacing="0" height="57" width="644"><tbody><tr><td><a href="http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2206641_2199021/2201650/1?PARTNER=3&amp;OAS_QUERY=null" target="_blank"><img alt="Rediff Shopping " border="0"></a></td>
</tr></tbody></table><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><br>