
<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: Times New Roman; font-size: 12pt; color: #000000'><P>Dear Suza</P>

<P>&nbsp;</P>

<P>It seems like BaTiO3 is not a simple cubic structure, so you should not choose ibrav=1.</P>

<P>&nbsp;</P>

<P>Chengyang<BR><BR></P>

<P>

<HR id=zwchr>

</P>

<BLOCKQUOTE style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: rgb(16,16,255) 2px solid"><B>From: </B>"Suza W" &lt;suza.rri@gmail.com&gt;<BR><B>To: </B>"PWSCF Forum" &lt;pw_forum@pwscf.org&gt;<BR><B>Sent: </B>Wednesday, November 30, 2011 3:32:15 PM<BR><B>Subject: </B>[Pw_forum] DC &amp; BEC<BR><BR><BR>Dear All, <BR><BR>Being an epitome of ferroelectric, BaTiO3 has extensively been <BR>studied using PWSCF. However,&nbsp; unfortunately, with several <BR>standard available USPP, I encounter&nbsp; following unphysical&nbsp; results :&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR><BR>&nbsp;&nbsp;&nbsp;&nbsp; Dielectric constant in cartesian axis <BR><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (-8142842.391767483&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000-8142842.391767481&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000-8142842.391767476 )<BR><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Effective charges (d Force / dE) in cartesian axis<BR><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; Ba <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ex&nbsp; (&nbsp;&nbsp;&nbsp; -1136.43937&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ey&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp; -1136.43937&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ez&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp; -1136.43937 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; Ti <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ex&nbsp; (&nbsp;&nbsp;&nbsp; -1874.03151&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ey&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp; -1874.03151&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ez&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp; -1874.03151 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp; O&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ex&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 931.42284&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ey&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp; 1067.39482&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ez&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp; 1067.39482 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp; O&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ex&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp; 1067.39482&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ey&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 931.42284&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ez&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp; 1067.39482 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp; O&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ex&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp; 1067.39482&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ey&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp; 1067.39482&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000 )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ez&nbsp; (&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 931.42284 )<BR><BR>Any suggestion for correcting these errors will be highly appreciated.&nbsp; <BR>Here goes the input file which renders such a blunder.&nbsp; <BR>&nbsp;&amp;system<BR>&nbsp;&nbsp;&nbsp; ibrav =&nbsp; 1,<BR>&nbsp;&nbsp;&nbsp; celldm(1) = 7.5589,<BR>&nbsp;&nbsp;&nbsp; nat = 5,<BR>&nbsp;&nbsp;&nbsp; ntyp=3,<BR>&nbsp;&nbsp;&nbsp; ecutwfc = 30.0,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 240.0,<BR>&nbsp;/<BR>&nbsp;&amp;electrons<BR>&nbsp;&nbsp;&nbsp; diagonalization = 'david'<BR>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-9<BR>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<BR>&nbsp;/<BR>&nbsp;&amp;ions<BR>&nbsp;/<BR>&nbsp;&amp;cell<BR>&nbsp;/<BR>ATOMIC_SPECIES<BR>&nbsp;Ba&nbsp; 137.327&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ba.pbe-nsp-van.UPF<BR>&nbsp;Ti&nbsp;&nbsp; 47.867&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ti.pbe-sp-van_ak.UPF<BR>&nbsp;O&nbsp;&nbsp;&nbsp; 15.9994&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; O.pbe-van_ak.UPF<BR>ATOMIC_POSITIONS (alat)<BR>Ba&nbsp; 0.0&nbsp; 0.0&nbsp; 0.0<BR>Ti&nbsp; 0.5&nbsp; 0.5&nbsp; 0.5<BR>O&nbsp; 0.0&nbsp; 0.5&nbsp; 0.5<BR>O&nbsp; 0.5&nbsp; 0.0&nbsp; 0.5<BR>O&nbsp; 0.5&nbsp; 0.5&nbsp; 0.0<BR>K_POINTS {automatic}<BR>&nbsp;8 8 8&nbsp; 1 1 1<BR><BR>Thanking you, <BR>Yours sincerely, <BR>Suza W<BR>PhD Student <BR>Department of Materials Science<BR>Bangalore<BR><BR>_______________________________________________<BR>Pw_forum mailing list<BR>Pw_forum@pwscf.org<BR>http://www.democritos.it/mailman/listinfo/pw_forum<BR></BLOCKQUOTE><BR></div></body></html>