Dear developers,<br><br>According Changelog for CVS version, it is possible now to perform Berry Phase calculation for gdir=1 and gdir=2 . I&#39;ve tried it for rhombohedral unit cell, but got the following message:<br><br>
....<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&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;&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; -1.00993858493275&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.583088313876004&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; 0.406263994517323&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; with crystal coordinates&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp; corresponds to ng=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; but G(ng)=&nbsp; 0.000000000000000E+000
<br>&nbsp; 0.000000000000000E+000&nbsp; 4.362056180567782E-319<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; 1&nbsp;&nbsp;&nbsp; G=&nbsp;&nbsp; 1.00993858493275&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.583088313876004
&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; 0.406263994517323&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; G=&nbsp; 0.000000000000000E+000&nbsp; -1.16617662775201&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; 0.812527989034646&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp; G=&nbsp; -1.00993858493275&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.74926494162801&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; 5.467848396278896E-014
<br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp; G=&nbsp;&nbsp; 2.01987716986549&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000000000E+000<br>&nbsp; 0.000000000000000E+000<br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp; G=&nbsp; -1.00993858493275&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.74926494162801&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 1.21879198355197&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp; G=&nbsp; 
0.000000000000000E+000&nbsp;&nbsp; 2.33235325550402&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;-0.406263994517268&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 7&nbsp;&nbsp;&nbsp; G=&nbsp; -2.01987716986549&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.16617662775201&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; 0.406263994517378&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp; G=&nbsp;&nbsp; 3.02981575479824
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.583088313876004&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;-0.406263994517323&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 9&nbsp;&nbsp;&nbsp; G=&nbsp;&nbsp; 1.00993858493275&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.91544156938002&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;-0.812527989034591&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp; G=&nbsp; -3.02981575479824&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.583088313876004
&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; 0.812527989034701&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp; G=&nbsp; -2.01987716986549&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.33235325550402&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 1.62505597806929&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 12&nbsp;&nbsp;&nbsp; G=&nbsp;&nbsp; 1.00993858493275&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.91544156938002&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; -2.84384796162148
&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 13&nbsp;&nbsp;&nbsp; G=&nbsp;&nbsp; 2.01987716986549&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.49852988325603&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; -1.21879198355191&nbsp;&nbsp;&nbsp;&nbsp; <br>&nbsp; i=&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 14&nbsp;&nbsp;&nbsp; G=&nbsp; -4.03975433973098&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000000000E+000<br>.....<br clear="all"><br>Lattice/reciprocal vectors are:
<br>&nbsp;&nbsp;&nbsp;&nbsp; crystal axes: (cart. coord. in units of a_0)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(1) = (&nbsp; 0.495080 -0.285834&nbsp; 0.820485 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(2) = (&nbsp; 0.000000&nbsp; 0.571669&nbsp; 0.820485 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; a(3) = ( -0.495080 -0.285834
&nbsp; 0.820485 )&nbsp; <br><br>&nbsp;&nbsp;&nbsp;&nbsp; reciprocal axes: (cart. coord. in units 2 pi/a_0)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(1) = (&nbsp; 1.009939 -0.583088&nbsp; 0.406264 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(2) = (&nbsp; 0.000000&nbsp; 1.166177&nbsp; 0.406264 )&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; b(3) = ( -
1.009939 -0.583088&nbsp; 0.406264 )&nbsp; <br><br>There is no error in case of gdir=3 .&nbsp; Does anybody have such problem before?<br><br>Best wishes,<br>&nbsp;<br>-- <br>===================================== <br>Jess Kondor<br><br>PICYT, Av. Venustiano Carranza 2425-A, 
<br>San Luis Potosí 78210, México<br><br>=====================================