<div dir="ltr"><div>Dear Users,<br><br></div>I am trying to perform berry phase calculations on pristine as well as doped BaTiO3 (<span class=""><span></span></span>tetragonal). I get reasonable results for the value of polarization for pristine BaTiO3. <br><br>------------------------------------------------------------------------------------------------------------------------<br>                           <span style="font-family:arial,helvetica,sans-serif">   SUMMARY OF PHASES<br>                               ~~~~~~~~~~~~~~~~~<br><br>                          Ionic Phase: -0.12600 (mod 2)<br>                     Electronic Phase:  0.40671 (mod 2)<br>                          TOTAL PHASE:  0.28071 (mod 2)<br><br><br>                             VALUES OF POLARIZATION<br>                             ~~~~~~~~~~~~~~~~~~~~~~<br><br>        The calculation of phases done along the direction of vector 3<br>        of the reciprocal lattice gives the following contribution to<br>        the polarization vector (in different units, and being Omega<br>        the volume of the unit cell):<br><br>           P =   2.1352035  (mod  15.2130988)  (e/Omega).bohr<br><br>           P =   0.0049373  (mod   0.0351775)  e/bohr^2<br><br>           P =   0.2822719  (mod   2.0111573)  C/m^2<br><br>        The polarization direction is:  ( 0.00000 , 0.00000 , 1.00000 )</span><br><br>------------------------------------------------------------------------------------------------------------------------<br><br>But in case of doped BaTiO3, I get negative value of polarization. I used a 2x2x3 supercell for these calculations.<br><br>-------------------------------------------------------------------------------------------------------------------------<br><br>                              SUMMARY OF PHASES<br>                               ~~~~~~~~~~~~~~~~~<br><br>                          Ionic Phase: -0.50398 (mod 2)<br>                      Electronic Phase: -0.37438 (mod 2)<br>                          TOTAL PHASE: -0.87836 (mod 2)<br><br><br>                             VALUES OF POLARIZATION<br>                             ~~~~~~~~~~~~~~~~~~~~~~<br><br>        The calculation of phases done along the direction of vector 3<br>        of the reciprocal lattice gives the following contribution to<br>        the polarization vector (in different units, and being Omega<br>        the volume of the unit cell):<br><br>           P = -19.9950561  (mod  45.5280337)  (e/Omega).bohr<br><br>           P =  -0.0038804  (mod   0.0088355)  e/bohr^2<br><br>           P =  -0.2218486  (mod   0.5051413)  C/m^2<br><br>        The polarization direction is:  ( 0.00000 , 0.00000 , 1.00000 )<br>--------------------------------------------------------------------------------------------------------------------------<br><br><div>I also tried using a finer k-mesh, but landed up with the same value of polarization. Does the negative polarization have any physical meaning or is it a numerical error? <br><br></div><div>Any help appreciated. Thanks in advance.<br><br></div><div>Best,<br></div><div>DS.<div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div class="gmail_signature">Deepashri Saraf<br></div><div class="gmail_signature">PhD Student<br></div><div class="gmail_signature">Department of Physics, SPPU<br></div><div class="gmail_signature">Pune (India)<br></div></div></div>
</div></div></div>