<div dir="ltr"><p class="MsoNormal" style="margin:0in 0in 8pt;line-height:107%;font-family:Calibri,sans-serif"><font size="4"><span style="line-height:107%;font-family:Arial,sans-serif;background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial">Dear Prof. Oleg Rubel,</span><span style="line-height:107%"></span></font></p>
<p class="MsoNormal" style="margin:0in 0in 8pt;line-height:107%;font-family:Calibri,sans-serif"><font size="4"><span style="line-height:107%;font-family:Arial,sans-serif;background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial">Thank you very much for your
valuable and complete reply. I found that the ionic and electronic
polarizations are summed but, in my opinion, if the phases are near to pi, it is
possible that the total polarization is bigger than half of the quantum of the
polarization. Thus, maybe the total polarization also needs to be wrapped.</span><span style="line-height:107%;font-family:Arial,sans-serif"><br>
<span style="background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial">Sincerely,</span><br>
<span style="background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial">Sh. Rahimi</span></span></font><span style="font-size:24pt;line-height:107%;font-size:24pt;line-height:107%"></span></p></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, May 29, 2021 at 12:15 AM <<a href="mailto:sjalali@sci.ui.ac.ir">sjalali@sci.ui.ac.ir</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><u></u>
<div style="font-family:Arial;font-size:14px">
</div>
<br><br><br>---------- Forwarded message ----------<br>From: "Rubel, Oleg" <<a href="mailto:rubelo@mcmaster.ca" target="_blank">rubelo@mcmaster.ca</a>><br>To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>Cc: <br>Bcc: <br>Date: Fri, 28 May 2021 18:46:41 +0000<br>Subject: Re: [Wien] wrap<br>You refer to an older version of BerryPI, which I do not have in place. I am inclined to think that the total phase is _not_ wrapped, but it would help to see your output tables (as below) before giving the exact answer. (please post to the list)<br>
<br>
As of WIEN2k v21, ionic and electronic phases are not added up, instead the corresponding polarizations are after converting into Cartesian coordinate. It is needed for generic (non-orthogonal) lattice vectors. This is how the new output (BaTiO3 tutorial) looks like:<br>
....<br>
CALCULATION OF ELECTRONIC POLARIZATION (primitive lattice coordinates)<br>
=======================================================================================<br>
Value | spin | dir(1) | dir(2) | dir(3)<br>
---------------------------------------------------------------------------------------<br>
Berry phase wrapped (rad) sp(1) [-2.660094e-13, 3.450573e-13, 1.258913e+00]<br>
Berry phase (rad) up+dn [-5.320189e-13, 6.901146e-13, 2.517826e+00]<br>
Berry phase wrapped (rad) up+dn [-5.320189e-13, 6.901146e-13, 2.517826e+00]<br>
Electronic polarization (C/m2) sp(1) [-8.415856e-14, 1.091673e-13, 4.024760e-01]<br>
=======================================================================================<br>
<br>
The electronic polarization vector is presented in this coord. system:<br>
dir(1) = [ 7.547566e+00, 4.621551e-16, 4.621551e-16] bohr<br>
dir(2) = [ 0.000000e+00, 7.547566e+00, 4.621551e-16] bohr<br>
dir(3) = [ 0.000000e+00, 0.000000e+00, 7.626934e+00] bohr<br>
and will be transformed into Cartesian coordinates.<br>
<br>
<br>
CALCULATION OF IONIC POLARIZATION (conventional lattice coordinates)<br>
=======================================================================================<br>
Elem.| Fractional coord. | spin | Zion | dir(1) | dir(2) | dir(3)<br>
---------------------------------------------------------------------------------------<br>
+------------ Ionic phase (rad) ------------+<br>
Ba (0.0000, 0.0000, 0.0000) sp(1) 10.00 [ 0.000000e+00, 0.000000e+00, 0.000000e+00]<br>
Ti (0.5000, 0.5000, 0.5152) sp(1) 12.00 [ 3.769911e+01, 3.769911e+01, 3.884323e+01]<br>
O (0.5000, 0.5000, 0.9736) sp(1) 6.00 [ 1.884956e+01, 1.884956e+01, 3.670240e+01]<br>
O (0.5000, 0.0000, 0.4834) sp(1) 6.00 [ 1.884956e+01, 0.000000e+00, 1.822516e+01]<br>
O (0.0000, 0.5000, 0.4834) sp(1) 6.00 [ 0.000000e+00, 1.884956e+01, 1.822516e+01]<br>
---------------------------------------------------------------------------------------<br>
Total ionic phase (rad) sp(1) [ 7.539822e+01, 7.539822e+01, 1.119960e+02]<br>
Total ionic phase wrap. (rad) sp(1) [ 3.552714e-15, 3.552714e-15, -1.101384e+00]<br>
Ionic polarization (C/m2) sp(1) [ 5.619937e-16, 5.619937e-16, -1.760570e-01]<br>
=======================================================================================<br>
The ionic positions and associated polarization vector is presented in this coord. system:<br>
dir(1) = [ 7.547566e+00, 4.621551e-16, 4.621551e-16] bohr<br>
dir(2) = [ 0.000000e+00, 7.547566e+00, 4.621551e-16] bohr<br>
dir(3) = [ 0.000000e+00, 0.000000e+00, 7.626934e+00] bohr<br>
<br>
<br>
SUMMARY OF POLARIZATION CALCULATION IN CARTESIAN COORDINATES<br>
=======================================================================================<br>
Value | spin | X | Y | Z<br>
---------------------------------------------------------------------------------------<br>
Electronic polarization (C/m2) sp(1) [-8.415856e-14, 1.091673e-13, 4.024760e-01]<br>
Ionic polarization (C/m2) sp(1) [ 5.619937e-16, 5.619937e-16, -1.760570e-01]<br>
Tot. spin polariz.=Pion+Pel (C/m2) sp(1) [-8.359656e-14, 1.097293e-13, 2.264190e-01]<br>
---------------------------------------------------------------------------------------<br>
TOTAL POLARIZATION (C/m2) both [-8.359656e-14, 1.097293e-13, 2.264190e-01]<br>
=======================================================================================<br>
<br>
As one can see, there is no "Total phase" only "Tot. spin polariz.=Pion+Pel".<br>
<br>
I hope it will help<br>
Oleg<br>
<br>
________________________________________<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of shahrbano rahimi <<a href="mailto:shahrbanorahimi60@gmail.com" target="_blank">shahrbanorahimi60@gmail.com</a>><br>
Sent: Friday, May 28, 2021 12:31<br>
To: <a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a><br>
Subject: [Wien] wrap<br>
<br>
Dear Prof. Rubel Oleg and WIEN2k users,<br>
I have a general question concerning Berry phase calculations.<br>
Would you, please, let me know that in the BerryPI package whether only the ionic and electronic parts of the phase are wrapped into the [-pi, pi] or the total phase is also wrapped into this interval. By this I mean that if the total phase after summing the wrapped ionic and electronic parts became larger than pi or smaller than -pi, whether the BerryPI package also wraps the total phase into the interval [-pi, pi] or leaves it to be located outside of the interval?<br>
The BerryPI version: 1.4.0 (May 25, 2018)<br>
The Wien2k version:18.2<br>
Sincerely,<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>