Hello,<div><br /></div><div>I am trying to calculate the dipole moment of my non-periodic system. The system has more than one molecule, and I am interested in calculating both molecular and total dipole moment. Currently, I use two sections (section1: FORCE_EVAL/DFT/LOCALIZE, section 2: FORCE_EVAL/DFT/PRINT/MOMENTS) in my input, as elaborated at the end of this message. Section 2 outputs a total dipole moment of 18.9465 Debye. From section 1 and using the Wannier centers, I get a total dipole moment of 18.9876 Debye. Here is my first question: Is this discrepancy normal? According to CP2K documentation, the first value (from section 2) is based on "simple operator". However, I am not sure what that is, and whether it justifies the discrepancy?</div><div><br /></div><div>Also, I have another question. In section 1, I have the TOTAL_DIPOLE and MOLECULAR_DIPOLES subsections as well. The dipole values that I get from these two are completely (orders of magnitude) different from the values mentioned in the previous paragraph (~19 Debye), although I set all the reference points to zero. Moreover, for the molecular dipole, the output file shows zero for most of the atoms, and non-zero values for the rest of them. Therefore, I cannot find the "molecule" concept there. Here is my second question: Am I missing anything when using these two subsections (TOTAL_DIPOLE and MOLECULAR_DIPOLES)?</div><div><br /></div><div>Thanks in advance for your help, and sorry for the long message!</div><div><br /></div><div><b>section 1:</b></div><div>    &LOCALIZE<br />      MAX_ITER 1000<br />      METHOD JACOBI<br />      EPS_LOCALIZATION 1.0E-6<br />      &PRINT<br />        &WANNIER_CENTERS<br />           IONS+CENTERS .TRUE.<br />        &END<br />        &TOTAL_DIPOLE<br />           PERIODIC F<br />           REFERENCE ZERO<br />           FILENAME total<br />        &END<br />        &MOLECULAR_DIPOLES<br />           PERIODIC F<br />           REFERENCE ZERO<br />           FILENAME molecular<br />        &END<br />      &END<br />    &END<br /></div><div><br /></div><div><b>section 2:</b></div><div>    &PRINT<br />      &MOMENTS</div><div>            PERIODIC F<br />            REFERENCE ZERO<br />            FILENAME dipole<br />      &END<br />    &END<br /></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/11c851cf-ed72-4d76-bee8-3afcd2ed3368n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/11c851cf-ed72-4d76-bee8-3afcd2ed3368n%40googlegroups.com</a>.<br />