<div dir="ltr">I have not carried out this calculation without spin polarization. However, spin polarization should be supported according to the source file header. Also, I believe not doing a spin-polarized calculation would yield an unphysical result, as the NV- center in diamond has a net spin. I suspect the simulation would not converge and bounce around the many close-to-degenerate answers in the case of spin unpolarized calculations? <div><br></div><div>Thank you,</div><div>Anil Bilgin</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Dec 19, 2019 at 7:12 AM Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</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"><div dir="ltr"><div>I can confirm that the latest version of the code yields the same results. No idea why. Do you get a similar behavior without spin polarization?</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 17, 2019 at 8:05 AM Anil Bilgin <<a href="mailto:bilgin@uchicago.edu" target="_blank">bilgin@uchicago.edu</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"><div dir="ltr">Hello all,<div><br></div><div>I am trying to calculate the permanent electric dipole moment of a negatively charged NV center in diamond. All NV centers are oriented along the 111 direction of diamond, which means that (by symmetry) contributions to polarizations along all three reciprocal lattice vectors should be equal in magnitude. While this is the case for my reciprocal lattice vectors along gdir = 2, 3; gdir = 1 gives me a vastly different answer. If anyone knows why this calculation is going wrong, or if this kind of calculation can / cannot be done in this way, please let me know. Other relevant information is below. </div><div><br></div><div>Thank you,</div><div>Anil Bilgin</div><div><br></div><div>PhD candidate</div><div>Pritzker School of Molecular Engineering</div><div>University of Chicago</div><div><br></div><div><br></div><div>I'm using QE version 6.1+intelmpi-5.1+intel-16.0</div><div><br></div><div>Little information about inputs: 'scf' calculations are done on an already relaxed structure with 215 atoms. Since the supercell is sufficiently large, k-point mesh is taken to be 1x1x1 (basically the gamma point).</div><div><br></div><div><br></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>