<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dear Miguel, <div class=""><br class=""></div><div class="">the „egg-box“ you are observing is mainly a matter of the PW cutoff, which by default is only 280 Ry, </div><div class="">but also due to the G-space differentiation within the XC-potential. Instead of the latter, this can be </div><div class="">ameliorated by the use of nearest neighbor smoothing and smoothed finite difference operators, e.g. </div><div class=""><br class=""></div><div class=""><pre class="code" style="margin-top: 0px; margin-bottom: 1.4em; padding: 0.7em 1em; font-family: Consolas, "Andale Mono WT", "Andale Mono", "Bitstream Vera Sans Mono", "Nimbus Mono L", Monaco, "Courier New", monospace; font-size: 14px; direction: ltr; background-color: rgb(251, 250, 249); color: rgb(51, 51, 51); box-shadow: rgb(204, 204, 204) -4px -4px 0.5em -0.3em inset; border-top-left-radius: 2px; border-top-right-radius: 2px; border-bottom-right-radius: 2px; border-bottom-left-radius: 2px; overflow: auto; word-wrap: normal; border: 1px solid rgb(204, 204, 204); caret-color: rgb(51, 51, 51);"> &XC_GRID
XC_SMOOTH_RHO NN50
XC_DERIV NN50_SMOOTH
&END XC_GRID</pre><div class="">Cheers, </div><div class="">Thomas</div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">Am 08.03.2021 um 13:41 schrieb Miguel Steiner <<a href="mailto:stein...@gmail.com" class="">stein...@gmail.com</a>>:</div><br class="Apple-interchange-newline"><div class="">Dear Prof. Hutter,<br class=""><br class="">thanks for your response, that obviously makes sense.<br class="">I guess this is even more exaggerated, due to the large box size? <br class="">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 class=""><br class=""><div class="">Thanks and best wishes,</div><div class="">Miguel<br class=""></div><br class=""><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 class=""></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 class="">
<br class="">due to the integration grid (PW grid) the QS energy is not
<br class="">translationally invariant. When you move a molecule rigidly
<br class="">through the box, the energy will oscillate. The frequency and
<br class="">amplitude depend on you PW cutoff.
<br class="">
<br class="">best regards
<br class="">
<br class="">Juerg Hutter
<br class="">--------------------------------------------------------------
<br class="">Juerg Hutter Phone : <a href="tel:+41%2044%20635%2044%2091" value="+41446354491" target="_blank" rel="nofollow" class="">++41 44 635 4491</a>
<br class="">Institut für Chemie C FAX : <a href="tel:+41%2044%20635%2068%2038" value="+41446356838" target="_blank" rel="nofollow" class="">++41 44 635 6838</a>
<br class="">Universität Zürich E-mail: <a href="" data-email-masked="" rel="nofollow" class="">h...@chem.uzh.ch</a>
<br class="">Winterthurerstrasse 190
<br class="">CH-8057 Zürich, Switzerland
<br class="">---------------------------------------------------------------
<br class="">
<br class="">-----<a href="" data-email-masked="" rel="nofollow" class="">c...@googlegroups.com</a> wrote: -----
<br class="">To: "cp2k" <<a href="" data-email-masked="" rel="nofollow" class="">c...@googlegroups.com</a>>
<br class="">From: "Miguel Steiner"
<br class="">Sent by: <a href="" data-email-masked="" rel="nofollow" class="">c...@googlegroups.com</a>
<br class="">Date: 03/03/2021 01:49PM
<br class="">Subject: [CP2K:14861] Energy differences with identical relative coordinates
<br class="">
<br class="">Dear Cp2k developers,
<br class="">
<br class="">just by coincidence I noticed that I get different electronic energies depending on the absolute coordinates of my input molecule.
<br class="">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 class="">As long as I center the coordinates, I get the identical electronic energy of
<br class="">-31.929620495691506
<br class="">But if I don't center the coordinates, I get different results:
<br class="">-31.929366281248029 (inside)
<br class="">-31.929298672694223 (outside)
<br class="">
<br class="">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 class="">
<br class="">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 class="">
<br class="">Input and systems see attachment.
<br class="">
<br class="">Thanks in advance and best regards,
<br class="">Miguel
<br class="">
<br class=""> --
<br class=""> You received this message because you are subscribed to the Google Groups "cp2k" group.
<br class=""> To unsubscribe from this group and stop receiving emails from it, send an email to <a href="" data-email-masked="" rel="nofollow" class="">cp...@googlegroups.com</a>.
<br class=""> 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" class="">https://groups.google.com/d/msgid/cp2k/32c783a9-f7fc-49c4-b013-693df6a3ca4an%40googlegroups.com</a>.
<br class="">
<br class="">
<br class="">[attachment "center-outside.inp" removed by Jürg Hutter/at/UZH]
<br class="">[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" class="">system-outside.xyz</a>" removed by Jürg Hutter/at/UZH]
<br class="">[attachment "center-inside.inp" removed by Jürg Hutter/at/UZH]
<br class="">[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" class="">system-inside.xyz</a>" removed by Jürg Hutter/at/UZH]
<br class="">[attachment "nocenter-inside.inp" removed by Jürg Hutter/at/UZH]
<br class="">[attachment "nocenter-outside.inp" removed by Jürg Hutter/at/UZH]
<br class=""></blockquote></div><div class=""><br class="webkit-block-placeholder"></div>
-- <br class="">
You received this message because you are subscribed to the Google Groups "cp2k" group.<br class="">
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp...@googlegroups.com" class="">cp...@googlegroups.com</a>.<br class="">
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d4fe0744-7591-4ee5-b8c1-12062cd7fd8bn%40googlegroups.com?utm_medium=email&utm_source=footer" class="">https://groups.google.com/d/msgid/cp2k/d4fe0744-7591-4ee5-b8c1-12062cd7fd8bn%40googlegroups.com</a>.<br class="">
</div></blockquote></div><br class=""></div></div><br class=""><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">==============================</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Thomas D. Kühne</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Dynamics of Condensed Matter</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Chair of Theoretical Chemistry</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">University of Paderborn</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Warburger Str. 100</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">D-33098 Paderborn</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">Germany</div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class=""><a href="mailto:tdku...@mail.upb.de" class="">tdku...@mail.upb.de</a></div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; orphans: 2; widows: 2;" class="">+49/(0)5251/60-5726</div></div></div></div>
</div>
<br class=""></body></html>