<P>
&nbsp; <BR>
Dear Sir,<BR>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; Thank you for your reply. I do understand that the sum of all the effective charges for a system should be zero but I would like to know why there is a discrepancy in calculating the ionic phase of the cubic( EXAMPLE10 OF ESPRESSO-4.0) and tetragonal PbTiO3 i.e. in cubic phase the total ionic phase should be -1.86 whereas it is 0.14(if it is sum of all the individual ionic phases of the constituent elements).<BR>
 <BR>
&nbsp; &nbsp; &nbsp; &nbsp;  In the tetra phase the total ionic phase is the sum of the individidual elements which is not so in cubic.I thus would like to know if this is the reason for me getting a negative value of effective charge for Pb for the displacement.I am attaching the ionic phase section for comparison.Please do enlighten me and thank you in advance<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>

</P>


SATHYA SHEELA.S <br>
Grad. Student <br>
Department of&nbsp;Physics  <br>
National Institute&nbsp;of Technology <br>
Tiruchirapalli -&nbsp;620015  <br>
India<br><br>
<Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><TR><td><a href='http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2205361_2197747/2200365/1?PARTNER=3&OAS_QUERY=null' target=new ><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2205361_2197747/creative_2200365.gif'  alt='Economist'  border=0></a></td></TR></Table>