<div dir="ltr"><div>Dear Nt,</div><div><br></div><div>Looking back at your previous questions, I think I might have misunderstood what you meant. If your question is about _printing_ the dipole moment using a non-periodic position operator, then you can simply choose:<br></div><div><br></div><div>PERIODIC .FALSE.</div><div><br></div><div>in this section:</div><div>https://manual.cp2k.org/cp2k-7_1-branch/CP2K_INPUT/FORCE_EVAL/PROPERTIES/LINRES/LOCALIZE/PRINT/TOTAL_DIPOLE.html</div><div><br></div><div>Of course you can always do this on snapshots extracted from the trajectory. Whether this is eventually consistent or not with your choice of boundary conditions is upon you to decide.</div><div><br></div><div>Best regards,</div><div>Hossam<br></div><br>On Tuesday, July 14, 2020 at 8:39:54 PM UTC+2, Nt wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Hossam,<div><br></div><div>Thank you very much indeed.  I'm relieved there's a solution.  Would you mind if I prepared an input file and ran it by you next week?  </div><div><br></div><div>Regards, Nt.</div></div></blockquote></div>