<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>