P.S. I just tried the 2022.1 version on our HPC centre and this also gives the same errors. So it is not a bug introduced in version 2023.1.<br /><br /><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 5 July 2023 at 15:23:13 UTC+2 Léon Luntadila Lufungula wrote:<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 Jurg,<div><br></div><div>I guess it's a bug in the 2023.1 version then, because I also just tested the most recent pre-compiled binary (<a href="https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp&source=gmail&ust=1688649795611000&usg=AOvVaw0bwRIyzIUjRMuRl9nBFdOK">https://github.com/cp2k/cp2k/releases/download/v2023.1/cp2k-2023.1-Linux-gnu-x86_64.ssmp</a>) and it still gives the same errors:</div><div><br><img alt="pre-compiled binary.png" width="653px" height="783px" src="https://groups.google.com/group/cp2k/attach/b0d9f4df8e132/pre-compiled%20binary.png?part=0.1&view=1"></div><div><br></div><div>Is a new version (2023.2) coming anytime soon perhaps? I tried compiling the trunk earlier this week after you fixed the CIF parser, both through cloning the Master from Git and using the tarball, but got errors during compilation... I'm going on vacation tomorrow but after I return I'll give it another try or if a newer version has released by then I'll try that and convince the HPC centre to install it.</div><div><br></div><div>Best,</div><div>Léon</div><div><br></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Wednesday, 5 July 2023 at 15:03:44 UTC+2 Jürg Hutter wrote:<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>it works for me (CP2K Trunk version)
<br>
<br>  Precond_solver : DEFAULT
<br>  stepsize       :    0.08000000                  energy_gap     :    0.08000000
<br>  eps_taylor     :   0.10000E-15                  max_taylor     :             4
<br>  ----------------------------------- OT ---------------------------------------
<br>
<br>  Step     Update method      Time    Convergence         Total energy    Change
<br>  ------------------------------------------------------------------------------
<br>     1 OT DIIS     0.80E-01    7.1     0.10459074       -83.0126092838 -8.30E+01
<br>     2 OT DIIS     0.80E-01    7.0     0.07934045       -85.4719301173 -2.46E+00
<br>     3 OT DIIS     0.80E-01    7.3     0.04237559       -86.5088026249 -1.04E+00
<br>     4 OT DIIS     0.80E-01    7.3     0.02451275       -87.1588661747 -6.50E-01
<br>     5 OT DIIS     0.80E-01    7.3     0.01463005       -87.2604337873 -1.02E-01
<br>     6 OT DIIS     0.80E-01    7.3     0.00701148       -87.3511713376 -9.07E-02
<br>     7 OT DIIS     0.80E-01    7.3     0.00477888       -87.3692546957 -1.81E-02
<br>     8 OT DIIS     0.80E-01    7.3     0.00289219       -87.3773030540 -8.05E-03
<br>
<br>This is running your input (removed the tabs and changed the file path).
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a rel="nofollow">cp...@googlegroups.com</a> <<a rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Sent: Wednesday, July 5, 2023 2:57 PM
<br>To: cp2k
<br>Subject: Re: [CP2K:19064] Unexpected 0D wavelet errors at large vacuum
<br>
<br>Just now, I also let CP2K print the atomic coordinates and the cell information so that I could visualize it in VMD to see if perhaps something goes wrong with reading in the structure. However, as you can see below, there is nothing wrong with the structure and there is ample vacuum space for the wavelet solver to work without throwing errors. I really don't understand how these errors occur for a structure which is essentially "over-qualified" for the wavelet method  with 10 Å vacuum on all sides... Is this perhaps a bug in the code behind the wavelet solver?
<br>[3mppa-vac-10.png]
<br>
<br>Best,
<br>Léon
<br>On Wednesday, 5 July 2023 at 14:07:38 UTC+2 Léon Luntadila Lufungula wrote:
<br>Hi Jurg,
<br>
<br>I still get the same errors when I alter the input as suggested, see the files in attachment. In which version did you do the calculation? My version of CP2K is 2023.1 and it seems to still throw errors even when altering the file as suggested.
<br>
<br>Kind regards,
<br>Léon
<br>
<br>On Wednesday, 5 July 2023 at 13:57:53 UTC+2 Jürg Hutter wrote:
<br>Hi
<br>
<br>works for me. Only change was to set the CELL in the input file:
<br>
<br>&CELL
<br>ABC [angstrom] 25.000 25.000 25.000
<br>PERIODIC NONE
<br>&END CELL
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a rel="nofollow">cp...@googlegroups.com</a> <<a rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Sent: Wednesday, July 5, 2023 1:41 PM
<br>To: cp2k
<br>Subject: Re: [CP2K:19060] Unexpected 0D wavelet errors at large vacuum
<br>
<br>Hi Jurg,
<br>
<br>Thanks for the tip, but I already tried this and it didn't work. I have tried once more and the results are in attachment... Any other suggestions?
<br>
<br>Best,
<br>Léon
<br>
<br>On Wednesday, 5 July 2023 at 13:03:51 UTC+2 Jürg Hutter wrote:
<br>Hi
<br>
<br>add center coordinates to your Topolgy section
<br>
<br>&TOPOLOGY
<br>COORD_FILE_NAME mol.cif
<br>COORD_FILE_FORMAT cif
<br>&CENTER_COORDINATES
<br>&END
<br>&END TOPOLOGY
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a rel="nofollow">cp...@googlegroups.com</a> <<a rel="nofollow">cp...@googlegroups.com</a>> on behalf of Léon Luntadila Lufungula <<a rel="nofollow">Leon.luntad...@uantwerpen.be</a>>
<br>Sent: Wednesday, July 5, 2023 12:03 PM
<br>To: cp2k
<br>Subject: [CP2K:19057] Unexpected 0D wavelet errors at large vacuum
<br>
<br>Hi all,
<br>
<br>I have a quick question regarding the 0D wavelet solver. I tried using it to calculate the energy of an isolated molecule and noticed some weird behaviour regarding the vacuum padding. At 2.5 Å vacuum around the molecule, I got an error saying "density non-zero on the edges of the unit cell: wrong results in WAVELET solver", so I increased the vacuum to 5 Å and this suppressed the error message. I wanted to see the effect of the vacuum on the energy so I also increased the vacuum to compare the energies. At 7.5 Å it all worked as expected, but then at 10 Å I again get errors stating that the density is non-zero at the edges, even when I include CENTER_COORDINATES. Can someone explain why this is happening? It seems very counterintuitive as it is exactly the same molecule and the only difference is that the vacuum padding around the molecule has increased...
<br>
<br>The outputs for 5, 7.5 and 10 Å are in attachment. The one with 2.5 Å, I seem to have deleted already as it resulting in the errors stated above.
<br>
<br>Best,
<br>Léon
<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 rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/589fb184-702d-485e-b911-5a75da8dae66n%40googlegroups.com" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/589fb184-702d-485e-b911-5a75da8dae66n%2540googlegroups.com&source=gmail&ust=1688649795612000&usg=AOvVaw2lTYAyB4zQQPdzNAvA-twR">https://groups.google.com/d/msgid/cp2k/589fb184-702d-485e-b911-5a75da8dae66n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/589fb184-702d-485e-b911-5a75da8dae66n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/589fb184-702d-485e-b911-5a75da8dae66n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688649795612000&usg=AOvVaw2iag9GS1UF6c7gJLofqRm7">https://groups.google.com/d/msgid/cp2k/589fb184-702d-485e-b911-5a75da8dae66n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<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 rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/e7995875-2d28-41ca-8eaa-d69146b28775n%40googlegroups.com" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/e7995875-2d28-41ca-8eaa-d69146b28775n%2540googlegroups.com&source=gmail&ust=1688649795612000&usg=AOvVaw1HiYtiD8QaJSOChqOG_g6M">https://groups.google.com/d/msgid/cp2k/e7995875-2d28-41ca-8eaa-d69146b28775n%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/e7995875-2d28-41ca-8eaa-d69146b28775n%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/e7995875-2d28-41ca-8eaa-d69146b28775n%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688649795612000&usg=AOvVaw0zIuhbnq9uK6zMShxfAYkA">https://groups.google.com/d/msgid/cp2k/e7995875-2d28-41ca-8eaa-d69146b28775n%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<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 rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/0b51bf55-96d9-4f6a-aa3a-3f2a77e6f5afn%40googlegroups.com" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/0b51bf55-96d9-4f6a-aa3a-3f2a77e6f5afn%2540googlegroups.com&source=gmail&ust=1688649795612000&usg=AOvVaw3SkLJMte-7qV513nmn2ZRf">https://groups.google.com/d/msgid/cp2k/0b51bf55-96d9-4f6a-aa3a-3f2a77e6f5afn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/0b51bf55-96d9-4f6a-aa3a-3f2a77e6f5afn%40googlegroups.com?utm_medium=email&utm_source=footer" rel="nofollow" target="_blank" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/0b51bf55-96d9-4f6a-aa3a-3f2a77e6f5afn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1688649795612000&usg=AOvVaw1YWE4WtfjRsQUQCydddN_M">https://groups.google.com/d/msgid/cp2k/0b51bf55-96d9-4f6a-aa3a-3f2a77e6f5afn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></blockquote></div></blockquote></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 on the web visit <a href="https://groups.google.com/d/msgid/cp2k/d59d2aad-c860-40d3-8102-f588ef489904n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d59d2aad-c860-40d3-8102-f588ef489904n%40googlegroups.com</a>.<br />