<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body dir="auto">
Dear <span style="-webkit-text-size-adjust: auto; caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); font-family: Arial, sans-serif; font-size: large; background-color: rgb(255, 255, 255);">Sh. Rahim,</span>
<div><font color="#000000" face="Arial, sans-serif" size="4"><span style="caret-color: rgb(0, 0, 0); -webkit-text-size-adjust: auto; background-color: rgb(255, 255, 255);"><br>
</span></font></div>
<div><font color="#000000" face="Arial, sans-serif" size="4"><span style="caret-color: rgb(0, 0, 0); -webkit-text-size-adjust: auto; background-color: rgb(255, 255, 255);">I would not worry about pi wrapping at this point and rather calculate a polarization
difference between two structures. Polarization is similar to the total energy. Only difference matters. To assure proper calculation (eliminate pi wrapping artefacts), sometimes we need to explore intermediate structures (smooth evolution between two structures).
Please check older posts in this list for more details.</span></font></div>
<div><font color="#000000" face="Arial, sans-serif" size="4"><span style="caret-color: rgb(0, 0, 0); -webkit-text-size-adjust: auto; background-color: rgb(255, 255, 255);"><br>
</span></font></div>
<div><font color="#000000" face="Arial, sans-serif" size="4"><span style="caret-color: rgb(0, 0, 0); -webkit-text-size-adjust: auto; background-color: rgb(255, 255, 255);">Best regards </span></font></div>
<div><font color="#000000" face="Arial, sans-serif" size="4"><span style="caret-color: rgb(0, 0, 0); -webkit-text-size-adjust: auto; background-color: rgb(255, 255, 255);">Oleg</span></font></div>
<div>
<div dir="ltr"><br>
<blockquote type="cite">On May 29, 2021, at 04:00, shahrbano rahimi <shahrbanorahimi60@gmail.com> wrote:<br>
<br>
</blockquote>
</div>
<blockquote type="cite">
<div dir="ltr">
<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>
<span>_______________________________________________</span><br>
<span>Wien mailing list</span><br>
<span>Wien@zeus.theochem.tuwien.ac.at</span><br>
<span>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</span><br>
<span>SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</span><br>
</div>
</blockquote>
</div>
</body>
</html>