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