<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; = 'scf'<BR>
&nbsp; &nbsp; restart_mode = 'from_scratch'<BR>
&nbsp; &nbsp; pseudo_dir&nbsp;  = '/home/nbaski/espresso-4.0.1/pseudo/'<BR>
&nbsp; &nbsp; outdir&nbsp; &nbsp; &nbsp;  = '/home/nbaski/tmp/'<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 = 'fixed'<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 border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><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&OAS_QUERY=null' target=new ><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2206641_2199021/creative_2201650.gif'  alt='Rediff Shopping '  border=0></a></td></TR></Table>