<P>
&nbsp; <BR>
Dear Developers,<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  I am a new user of Espresso-4.0 and when I ran the example10 (polarizayion of cubic PbTiO3 via Berry phase calculation)I found that the total ionic phase is not the sum of the ionic phases of the individiual constituents i.e Pb ,Ti ,O but in the tetragonal PbTiO3 the total ionic phase is -0.22375 which is the individual sum of Pb ,Ti, O.<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; I am attaching the O/P files (only a section as the complete file is bulky)of cubic and tetra phase of PbTiO3 ( polarization calculation ) for comparison. <BR>
<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; I would be really thankful for any explanations or reasoning in this context.<BR>
<BR>
<BR>
FOR CUBIC PHASE<BR>
<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; IONIC POLARIZATION<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ~~~~~~~~~~~~~~~~~~<BR>
<BR>
&nbsp; &nbsp; &nbsp; &nbsp; Note: (mod 1) means that the phases (angles ranging from<BR>
&nbsp; &nbsp; &nbsp; &nbsp; -pi to pi) have been mapped to the interval [-1/2,+1/2) by<BR>
&nbsp; &nbsp; &nbsp; &nbsp; dividing by 2*pi; (mod 2) refers to the interval [-1,+1)<BR>
<BR>
&nbsp; ============================================================================<BR>
&nbsp; &nbsp; Ion&nbsp; &nbsp; Species&nbsp; &nbsp; Charge&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Position&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Phase<BR>
&nbsp; ----------------------------------------------------------------------------<BR>
&nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp; &nbsp; Pb&nbsp; &nbsp; &nbsp; 14.000&nbsp; &nbsp; &nbsp; 0.0000&nbsp; 0.0000&nbsp; 0.0100&nbsp; &nbsp; 0.14000 (mod 2)<BR>
&nbsp; &nbsp; 2&nbsp; &nbsp; &nbsp; &nbsp; Ti&nbsp; &nbsp; &nbsp; 12.000&nbsp; &nbsp; &nbsp; 0.5000&nbsp; 0.5000&nbsp; 0.5000&nbsp; &nbsp; 0.00000 (mod 2)<BR>
&nbsp; &nbsp; 3&nbsp; &nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; &nbsp; 6.000&nbsp; &nbsp; &nbsp; 0.0000&nbsp; 0.5000&nbsp; 0.5000&nbsp; &nbsp; -1.00000 (mod 2)<BR>
&nbsp; &nbsp; 4&nbsp; &nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; &nbsp; 6.000&nbsp; &nbsp; &nbsp; 0.5000&nbsp; 0.5000&nbsp; 0.0000&nbsp; &nbsp; 0.00000 (mod 2)<BR>
&nbsp; &nbsp; 5&nbsp; &nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; &nbsp; 6.000&nbsp; &nbsp; &nbsp; 0.5000&nbsp; 0.0000&nbsp; 0.5000&nbsp; &nbsp; -1.00000 (mod 2)<BR>
&nbsp; ----------------------------------------------------------------------------<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; IONIC PHASE:&nbsp; 0.14000 (mod 2)<BR>
&nbsp; ============================================================================<BR>
<BR>
<BR>
<BR>
FOR TETRAGONAL PbTiO3<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; IONIC POLARIZATION<BR>
<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ~~~~~~~~~~~~~~~~~~<BR>
<BR>
&nbsp; &nbsp; &nbsp; &nbsp; Note: (mod 1) means that the phases (angles ranging from<BR>
&nbsp; &nbsp; &nbsp; &nbsp; -pi to pi) have been mapped to the interval [-1/2,+1/2) by<BR>
&nbsp; &nbsp; &nbsp; &nbsp; dividing by 2*pi; (mod 2) refers to the interval [-1,+1)<BR>
<BR>
&nbsp; ============================================================================<BR>
&nbsp; &nbsp; Ion&nbsp; &nbsp; Species&nbsp; &nbsp; Charge&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Position&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Phase<BR>
&nbsp; ----------------------------------------------------------------------------<BR>
&nbsp; &nbsp; 1&nbsp; &nbsp; &nbsp; &nbsp; Pb&nbsp; &nbsp; &nbsp; 14.000&nbsp; &nbsp; &nbsp; 0.0000&nbsp; 0.0000&nbsp; 0.0100&nbsp; &nbsp; 0.13142 (mod 2)<BR>
&nbsp; &nbsp; 2&nbsp; &nbsp; &nbsp; &nbsp; Ti&nbsp; &nbsp; &nbsp; 12.000&nbsp; &nbsp; &nbsp; 0.5000&nbsp; 0.5000&nbsp; 0.5377&nbsp; &nbsp; 0.05689 (mod 2)<BR>
&nbsp; &nbsp; 3&nbsp; &nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; &nbsp; 6.000&nbsp; &nbsp; &nbsp; 0.5000&nbsp; 0.5000&nbsp; 0.1118&nbsp; &nbsp; 0.62968 (mod 2)<BR>
&nbsp; &nbsp; 4&nbsp; &nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; &nbsp; 6.000&nbsp; &nbsp; &nbsp; 0.0000&nbsp; 0.5000&nbsp; 0.6174&nbsp; &nbsp; -0.52267 (mod 2)<BR>
&nbsp; &nbsp; 5&nbsp; &nbsp; &nbsp; &nbsp; O&nbsp; &nbsp; &nbsp; &nbsp; 6.000&nbsp; &nbsp; &nbsp; 0.5000&nbsp; 0.0000&nbsp; 0.6174&nbsp; &nbsp; -0.52267 (mod 2)<BR>
&nbsp; ----------------------------------------------------------------------------<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; IONIC PHASE:&nbsp; -0.22735 (mod 2)<BR>
&nbsp; ============================================================================<BR>
<BR>
<BR>
Thanking You<BR>
<BR>
<BR>
SATHYA SHEELA.S <BR>
Grad. Student <BR>
Department of Physics <BR>
National Institute of Technology <BR>
Tiruchirapalli - 620015 <BR>
India<BR>
<BR>
 <BR>
&nbsp; <BR>
 <BR>
<BR>
&nbsp;  <BR>
 <BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <BR>
&nbsp; <BR>
 <BR>
 <BR>
 <BR>

</P>
<br><br>