<table cellspacing='0' cellpadding='0' border='0' ><tr><td valign='top' style='font: inherit;'><P>Dear all</P>
<P>&nbsp;</P>
<P>I calculated the lattice constant of BaTiO3 using GGA. However, the calculated lattice constant ( 3.980 angstrom) is SMALLER than experimental results (4.012 angstrom). So is this result reasonable?</P>
<P>&nbsp;</P>
<P>The&nbsp;GGA results of SrTiO3 is ok, 3.928 angstrom&nbsp;vs expermental result &nbsp;3.905. As SrTiO3 and BaTiO3 both contain Ti and O,&nbsp;I calculated lattice constant of bulk metal Ba as a test, and the&nbsp;GGA results (5.014 angstrom for bcc cell) is also&nbsp;smaller than experimental one (5.019).</P>
<P>&nbsp;</P>
<P>So is there something&nbsp;wrong with my calculation? Or GGA does underestimate the lattice constant of some metals, such as reported in Computer Coupling of Phase Diagrams and Thermochemistry 28(2004)79-90 ?</P>
<P>&nbsp;</P>
<P>The input file is attached below. Thank you very much for your help </P>
<P>Ruizhi&nbsp;<BR></P>
<P>===========input file of BaTiO3===============</P>
<P>&amp;control<BR>&nbsp;&nbsp;&nbsp; calculation&nbsp; = 'vc-relax'<BR>&nbsp;&nbsp;&nbsp; restart_mode = 'from_scratch'<BR>&nbsp;&nbsp;&nbsp; prefix = 'bto'<BR>&nbsp;&nbsp;&nbsp; pseudo_dir&nbsp;&nbsp; = '~/pseudo/'<BR>&nbsp;&nbsp;&nbsp; outdir&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = './tmp/'<BR>&nbsp;/<BR>&nbsp;&amp;system<BR>&nbsp;&nbsp;&nbsp; ibrav=1<BR>&nbsp;&nbsp;&nbsp; celldm(1)=7.584<BR>&nbsp;&nbsp;&nbsp; nat=5<BR>&nbsp;&nbsp;&nbsp; ntyp=3<BR>&nbsp;&nbsp;&nbsp; ecutwfc=30.0<BR>&nbsp;&nbsp;&nbsp; occupations = 'smearing'<BR>&nbsp;&nbsp;&nbsp; degauss=0.01<BR>&nbsp;/<BR>&nbsp;&amp;electrons<BR>&nbsp;&nbsp;&nbsp; conv_thr = 1e-8<BR>&nbsp;&nbsp;&nbsp; mixing_beta=0.3<BR>&nbsp;/<BR>&amp;IONS<BR>&nbsp; ion_dynamics = 'damp',<BR>&nbsp;/<BR>&nbsp;&amp;CELL<BR>&nbsp; cell_dynamics ='damp-pr'<BR>&nbsp;<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>&nbsp; Ba&nbsp;&nbsp;&nbsp; 137.33&nbsp; Ba.pbe-nsp-van.UPF <BR>&nbsp; Ti&nbsp;&nbsp;&nbsp; 47.90&nbsp;&nbsp;
 Ti.pbe-sp-van_ak.UPF<BR>&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; 15.99&nbsp;&nbsp; O.pbe-van_ak.UPF<BR>ATOMIC_POSITIONS (crystal)<BR>&nbsp; Ba&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp; 0.000<BR>&nbsp; Ti&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp; 0.500<BR>&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp; 0.500<BR>&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp; 0.000<BR>&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp; 0.500&nbsp;&nbsp;&nbsp; 0.000&nbsp;&nbsp;&nbsp; 0.500<BR>K_POINTS {automatic}<BR>&nbsp; 8 8 8 1 1 1</P>
<P>=============end input file=============</P>
<P><BR>========================<BR>Ruizhi Zhang, PhD<BR>School of Physics,Shandong University</P>
<P>Jinan, Peopel's Republic of China</P></td></tr></table><br>