Dear Prof. Hutter,<br><br>thanks for your response, that obviously makes sense.<br>I guess this is even more exaggerated, due to the large box size? <br>Do you (and your team) have experience what box size fits best to simulate a molecular calculation. I need those to calculate adsorption energies.<br><br><div>Thanks and best wishes,</div><div>Miguel<br></div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">jgh schrieb am Mittwoch, 3. März 2021 um 15:06:54 UTC+1:<br/></div><blockquote class="gmail_quote" style="margin: 0 0 0 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi
<br>
<br>due to the integration grid (PW grid) the QS energy is not
<br>translationally invariant. When you move a molecule rigidly
<br>through the box, the energy will oscillate. The frequency and
<br>amplitude depend on you PW cutoff.
<br>
<br>best regards
<br>
<br>Juerg Hutter
<br>--------------------------------------------------------------
<br>Juerg Hutter                         Phone : <a href="tel:+41%2044%20635%2044%2091" value="+41446354491" target="_blank" rel="nofollow">++41 44 635 4491</a>
<br>Institut für Chemie C                FAX   : <a href="tel:+41%2044%20635%2068%2038" value="+41446356838" target="_blank" rel="nofollow">++41 44 635 6838</a>
<br>Universität Zürich                   E-mail: <a href data-email-masked rel="nofollow">h...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>---------------------------------------------------------------
<br>
<br>-----<a href data-email-masked rel="nofollow">c...@googlegroups.com</a> wrote: -----
<br>To: "cp2k" <<a href data-email-masked rel="nofollow">c...@googlegroups.com</a>>
<br>From: "Miguel Steiner" 
<br>Sent by: <a href data-email-masked rel="nofollow">c...@googlegroups.com</a>
<br>Date: 03/03/2021 01:49PM
<br>Subject: [CP2K:14861] Energy differences with identical relative coordinates
<br>
<br>Dear Cp2k developers,
<br>
<br>just by coincidence I noticed that I get different electronic energies depending on the absolute coordinates of my input molecule.
<br>I tried a simple test system with O2 in a 15x15x15x angstrom cube. I set-up 2 different xyz files, one with 3.xx angstrom in the z coordinates and one with 34.xx angstrom in the z coordinate. However, both O2 molecules are identical in their bond length, they are simply translated to each other.
<br>As long as I center the coordinates, I get the identical electronic energy of
<br>-31.929620495691506
<br>But if I don't center the coordinates, I get different results:
<br>-31.929366281248029 (inside)
<br>-31.929298672694223 (outside)
<br>
<br>I mean, I can imagine that it leads to difficulties, if the molecule is placed outside of the defined cell, but especially the different energy for the O2 molecules that is placed within the cell, albeit not the exact center, confuses me.
<br>
<br>Am I doing something wrong here or is this expected behavior? If the latter is the case, why aren't all coordinates centered by default (which I guess I will do from now on)?
<br>
<br>Input and systems see attachment.
<br>
<br>Thanks in advance and best regards,
<br>Miguel
<br>  
<br>  -- 
<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 data-email-masked rel="nofollow">cp...@googlegroups.com</a>.
<br> To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/32c783a9-f7fc-49c4-b013-693df6a3ca4an%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=de&q=https://groups.google.com/d/msgid/cp2k/32c783a9-f7fc-49c4-b013-693df6a3ca4an%2540googlegroups.com&source=gmail&ust=1615293404287000&usg=AFQjCNHWsfPFBQXmcSwDId6C_L9Vs9Tfig">https://groups.google.com/d/msgid/cp2k/32c783a9-f7fc-49c4-b013-693df6a3ca4an%40googlegroups.com</a>.
<br> 
<br>
<br>[attachment "center-outside.inp" removed by Jürg Hutter/at/UZH]
<br>[attachment "<a href="http://system-outside.xyz" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=de&q=http://system-outside.xyz&source=gmail&ust=1615293404287000&usg=AFQjCNHobI9lhmJE-8Ug05_IecnKF285ww">system-outside.xyz</a>" removed by Jürg Hutter/at/UZH]
<br>[attachment "center-inside.inp" removed by Jürg Hutter/at/UZH]
<br>[attachment "<a href="http://system-inside.xyz" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=de&q=http://system-inside.xyz&source=gmail&ust=1615293404287000&usg=AFQjCNGeeG9ChkJ3uUn7L2F6YpeuX7Dyhg">system-inside.xyz</a>" removed by Jürg Hutter/at/UZH]
<br>[attachment "nocenter-inside.inp" removed by Jürg Hutter/at/UZH]
<br>[attachment "nocenter-outside.inp" removed by Jürg Hutter/at/UZH]
<br></blockquote></div>