<div dir="ltr">I would like to add that EuTiO3 seems to be a multiferroic compound with ground state being antiferromagnetic- paraelectric system (I4/mcm). Then, under some compressive or tensile epitaxial strain, several types of transition is possible (For example, antiferromagnetic- ferroelectric or ferromagnetic- ferroelectric) depending on the magnitude and type of strain. The following articles address this issue very well.<div><a href="http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.109.267602" target="_blank">http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.109.267602</a><br><div>  <br><a href="http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.97.267602" target="_blank">http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.97.267602</a><br><br><div><div><br>Hope that helps.<br><br>Sheikh<br><br>-- <br><div><div dir="ltr"><div>Sheikh Jamil Ahmed</div><div>PhD Student</div><div>Department of Material Science and Engineering</div><div><span style="font-family:Tahoma,sans-serif;line-height:15.6px;text-align:center;text-indent:9pt">McMaster University</span></div><div><span style="font-family:Tahoma,sans-serif;line-height:15.6px;text-align:center;text-indent:9pt">1280 Main Street West</span></div><div><span style="font-family:Tahoma,sans-serif;line-height:15.6px;text-align:center;text-indent:9pt">Hamilton, Ontario, Canada L8S 4L7</span></div><div><br></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 8, 2015 at 3:07 PM, Oleg Rubel <span dir="ltr"><<a href="mailto:orubel@lakeheadu.ca" target="_blank">orubel@lakeheadu.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">To address your questions...<br>
<span class=""><br>
> I4/mcm is centrosymmetric. I studied the tutorial related to GaN and it helped me very much. So I conclude to calculate the difference in polarization, two structures that at least one of them should not be centrosymmetric are needed, Is it correct?<br>
<br>
</span>That’s correct. When the structure file contains the inversion symmetry, Wien2k switches from complex to real version for wave functions. Then there is no phase information to be processed.<br>
<span class=""><br>
> Would you please guide me whether it is allowed to calculate the difference in polarization between a strained and unstrained state of a system by implementing the berryPI on each of them separately and then subtract them ?! if it going to, what does the derived polarization  describe?! Piezoelectricity(since it is the difference in polarization caused by strain) or spontaneous polarization?<br>
<br>
</span>You are heading towards piezoelectric coefficients. There are some technicalities, such as clamped vs relaxed ion approximation, proper vs improper. Here is a good reference: <a href="http://arxiv.org/pdf/cond-mat/9903137.pdf" rel="noreferrer" target="_blank">http://arxiv.org/pdf/cond-mat/9903137.pdf</a><br>
I do not see any barriers from BerryPI perspective.<br>
<br>
<br>
Best regards<br>
<span class="HOEnZb"><font color="#888888">Oleg<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">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>
</div></div></blockquote></div><br></div>