Dear All,<br><br>I have noticed that  calculation of polarization using Berry Phase method and pwscf-4.0 doesn't work in parallel for gdir=1 and gdir=2. For example, I took tetragonal phase of BaTiO3 with atomic displacements along 'x' direction and (c,a,a) lattice vectors. For nprocs =1, it works fine:<br>
<br>VALUES OF POLARIZATION<br>                             ~~~~~~~~~~~~~~~~~~~~~~<br><br>        The calculation of phases done along the direction of vector 1<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.2164982  (mod  15.2155768)  (e/Omega).bohr<br><br>           P =   0.0051349  (mod   0.0352497)  e/bohr^2<br>
<br>           P =   0.2935730  (mod   2.0152882)  C/m^2<br><br>        The polarization direction is:  ( 1.00000 , 0.00000 , 0.00000 )<br><br clear="all"><br>but for nprocs > 1 (for example, 16) it doesn't:<br><br>
==================================================<br>                            POLARIZATION CALCULATION<br>                         !!! NOT THOROUGHLY TESTED !!!<br>               --------------------------------------------------<br>
<br>  error: translated G=  -0.2483941319370271         0.000000000000000      <br>    0.000000000000000       with crystal coordinates           -1            0 <br>            0  corresponds to ng=            0  but G(ng)= <br>
   2.8071327799762095E-319    0.000000000000000        1.4927122354773060E-316<br>  probably because G_par is NOT a reciprocal lattice vector <br>  Possible choices as smallest  G_par:<br>  i=            1    G=    0.000000000000000         0.000000000000000      <br>
    0.000000000000000     <br>  i=            2    G=    0.000000000000000         0.000000000000000      <br>  -0.2508780732563974     <br>  i=            3    G=    0.000000000000000         0.000000000000000      <br>   0.2508780732563974     <br>
  i=            4    G=    0.000000000000000         0.000000000000000      <br>  -0.5017561465127947     <br>  i=            5    G=    0.000000000000000         0.000000000000000      <br>   0.5017561465127947     <br>  i=            6    G=   0.4967882638740543       -0.2508780732563974      <br>
    0.000000000000000     <br>  i=            7    G=   0.4967882638740543       -0.2508780732563974      <br>  -0.2508780732563974     <br>.......<br><br>is there a way to solve this problem?<br><br> Best regards,<br>   <br>
-- <br>===================================== <br>Jess Kondor, PhD student<br><br>PICYT, Av. Venustiano Carranza 2425-A, <br>San Luis Potosí 78210, México<br><br>=====================================