<br>Dear All, <br><br>Being an epitome of ferroelectric, BaTiO3 has extensively been <br>studied using PWSCF. However, unfortunately, with several <br>standard available USPP, I encounter following unphysical results : <br>
<br> Dielectric constant in cartesian axis <br><br> (-8142842.391767483 0.000000000 0.000000000 )<br> ( 0.000000000-8142842.391767481 0.000000000 )<br> ( 0.000000000 0.000000000-8142842.391767476 )<br>
<br> Effective charges (d Force / dE) in cartesian axis<br><br> atom 1 Ba <br> Ex ( -1136.43937 0.00000 0.00000 )<br> Ey ( 0.00000 -1136.43937 0.00000 )<br>
Ez ( 0.00000 0.00000 -1136.43937 )<br> atom 2 Ti <br> Ex ( -1874.03151 0.00000 0.00000 )<br> Ey ( 0.00000 -1874.03151 0.00000 )<br> Ez ( 0.00000 0.00000 -1874.03151 )<br>
atom 3 O <br> Ex ( 931.42284 0.00000 0.00000 )<br> Ey ( 0.00000 1067.39482 0.00000 )<br> Ez ( 0.00000 0.00000 1067.39482 )<br> atom 4 O <br>
Ex ( 1067.39482 0.00000 0.00000 )<br> Ey ( 0.00000 931.42284 0.00000 )<br> Ez ( 0.00000 0.00000 1067.39482 )<br> atom 5 O <br> Ex ( 1067.39482 0.00000 0.00000 )<br>
Ey ( 0.00000 1067.39482 0.00000 )<br> Ez ( 0.00000 0.00000 931.42284 )<br><br>Any suggestion for correcting these errors will be highly appreciated. <br>Here goes the input file which renders such a blunder. <br>
&system<br> ibrav = 1,<br> celldm(1) = 7.5589,<br> nat = 5,<br> ntyp=3,<br> ecutwfc = 30.0,<br> ecutrho = 240.0,<br> /<br> &electrons<br> diagonalization = 'david'<br> conv_thr = 1.0d-9<br>
mixing_beta = 0.7<br> /<br> &ions<br> /<br> &cell<br> /<br>ATOMIC_SPECIES<br> Ba 137.327 Ba.pbe-nsp-van.UPF<br> Ti 47.867 Ti.pbe-sp-van_ak.UPF<br> O 15.9994 O.pbe-van_ak.UPF<br>ATOMIC_POSITIONS (alat)<br>
Ba 0.0 0.0 0.0<br>Ti 0.5 0.5 0.5<br>O 0.0 0.5 0.5<br>O 0.5 0.0 0.5<br>O 0.5 0.5 0.0<br>K_POINTS {automatic}<br> 8 8 8 1 1 1<br><br>Thanking you, <br>Yours sincerely, <br>Suza W<br>PhD Student <br>Department of Materials Science<br>
Bangalore<br>