Dear CP2K developers and users,<br /><div><p>I'm running into an issue with <span dir="ltr">MOLECULAR_MOMENTS</span> after Wannier localization and was wondering if anyone has seen something similar.</p><p>I'm using UKS (also reproduced with DFT+U enabled) together with:</p><span dir="ltr"><span dir="ltr">&LOCALIZE
<br /> METHOD JACOBI
<br /> USE_HISTORY F
<br /> &PRINT
    &MOLECULAR_MOMENTS LOW
<br />    ORDER 1
<br /> &END MOLECULAR_MOMENTS
<br />   &WANNIER_CENTERS
<br />     IONS+CENTERS
<br />   &END WANNIER_CENTERS
<br /> &END PRINT
&END LOCALIZE</span></span><p>The SCF converges normally, the Jacobi localization converges, and the printed Wannier centers for both spin channels are finite and appear reasonable. However, the <span dir="ltr">MOLECULAR_MOMENTS</span> output consists entirely of <span dir="ltr">NaN</span> values for every atom/molecule.</p><p>As a sanity check, I wrote a small post-processing script that assigns the printed Wannier centers to the nearest atom (using minimum-image PBC) and reconstructs per-atom electronic dipoles which are all reasonable.</p><p>Has anyone encountered this before, or is there a known issue with <span dir="ltr">MOLECULAR_MOMENTS</span> for localized UKS wavefunctions? </p></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/92857f51-b8d0-4551-8467-e45d87cd39can%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/92857f51-b8d0-4551-8467-e45d87cd39can%40googlegroups.com</a>.<br />