Thank you for your detailed and professional response. I understand now!<br />This computational result was obtained when I was investigating the performance of different k-point calculation settings. <div>In principle, the gamma point calculation setting is reasonable for this system.</div><div><br /></div><div>best wishes,</div><div>Dongfei<br /><br /></div><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Monday, 19 May 2025 at 22:53:17 UTC+8 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>These are the total forces for 64 waters from the ripples.
<br>(BTW no k-points here, but why do you have a 3 3 1 grid?)
<br>
<br>GPW PBE      0.011
<br>GPW LDA      0.002
<br>GAPW PBE    0.002
<br>GAPW LDA    0.0001
<br>
<br>Keep in mind that these forces are correct forces for the energy functional defined in CP2K.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of liu dongfei <<a href data-email-masked rel="nofollow">liud...@gmail.com</a>>
<br>Sent: Monday, May 19, 2025 4:28 PM
<br>To: cp2k
<br>Subject: Re: [CP2K:21478] The total force calculated by CP2K is not zero
<br>
<br>Thanks for your reply!
<br>The attached file is the result of a calculation, showing that the total force of the system is :
<br>-------------------------------------------------------
<br>....
<br>    189      2      H          -0.02270940    -0.01715719    -0.01124088
<br>    190      2      H          -0.00819690     0.01778486    -0.00220956
<br>    191      2      H          -0.00491958    -0.00982567     0.00540266
<br>    192      2      H          -0.00744376     0.00269638    -0.00803075
<br> SUM OF ATOMIC FORCES           0.00405273    -0.00782982     0.01173920     0.01468126
<br>-------------------------------------------------------
<br>The total force here is not particularly small, and it is comparable to the forces experienced by some H atoms.
<br>Is this result reasonable, even considering the ripple effect?
<br>
<br>best wishes,
<br>
<br>Dongfei
<br>
<br>On Monday, 19 May 2025 at 21:26:00 UTC+8 Jürg Hutter wrote:
<br>Hi
<br>
<br>inaccurate numerics can be a source for forces not summing up to zero. However, I assume
<br>(we don't know actually, as you don't share the forces) that the residual force is due to the
<br>'ripple effect'. CP2K uses aequidistant grids to calculate XC energies. For such grids there is no
<br>longer any translational invariance.
<br>
<br>regards
<br>JH
<br>
<br>________________________________________
<br>From: <a href data-email-masked rel="nofollow">cp...@googlegroups.com</a> <<a href data-email-masked rel="nofollow">cp...@googlegroups.com</a>> on behalf of liu dongfei <<a href data-email-masked rel="nofollow">liud...@gmail.com</a>>
<br>Sent: Monday, May 19, 2025 5:00 AM
<br>To: cp2k
<br>Subject: [CP2K:21474] The total force calculated by CP2K is not zero
<br>
<br>Hi all, I’m currently simulating a pure water system with CP2K. Despite configuring appropriate convergence criteria (e.g., CUTOFF, REL_CUTOFF, and EPS_SCF), the total forces from the calculation do not converge to zero. Attached is my input configuration.
<br>Could anyone provide insights into potential causes or recommend literature discussing this specific issue?
<br>Thanks in advance!
<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">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/ddbe59a0-1631-4bec-9854-724d785b793bn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/ddbe59a0-1631-4bec-9854-724d785b793bn%2540googlegroups.com&source=gmail&ust=1747754835094000&usg=AOvVaw3AQ2RFQd067OvD4Jduanil">https://groups.google.com/d/msgid/cp2k/ddbe59a0-1631-4bec-9854-724d785b793bn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/ddbe59a0-1631-4bec-9854-724d785b793bn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/ddbe59a0-1631-4bec-9854-724d785b793bn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1747754835094000&usg=AOvVaw18Na2R5LuSLkAIfa4WtlMo">https://groups.google.com/d/msgid/cp2k/ddbe59a0-1631-4bec-9854-724d785b793bn%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 href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a><mailto:<a href data-email-masked rel="nofollow">cp2k+uns...@googlegroups.com</a>>.
<br>To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/8c6b7571-cd57-482c-9715-0fb08598e50bn%40googlegroups.com" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/8c6b7571-cd57-482c-9715-0fb08598e50bn%2540googlegroups.com&source=gmail&ust=1747754835094000&usg=AOvVaw3l6f60G7nrzSG7eVG8UO1X">https://groups.google.com/d/msgid/cp2k/8c6b7571-cd57-482c-9715-0fb08598e50bn%40googlegroups.com</a><<a href="https://groups.google.com/d/msgid/cp2k/8c6b7571-cd57-482c-9715-0fb08598e50bn%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank" rel="nofollow" data-saferedirecturl="https://www.google.com/url?hl=en-GB&q=https://groups.google.com/d/msgid/cp2k/8c6b7571-cd57-482c-9715-0fb08598e50bn%2540googlegroups.com?utm_medium%3Demail%26utm_source%3Dfooter&source=gmail&ust=1747754835094000&usg=AOvVaw2es88cxVn_4IX3haziCnZ6">https://groups.google.com/d/msgid/cp2k/8c6b7571-cd57-482c-9715-0fb08598e50bn%40googlegroups.com?utm_medium=email&utm_source=footer</a>>.
<br></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 visit <a href="https://groups.google.com/d/msgid/cp2k/6c326058-6280-4897-a11c-fa1137398549n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/6c326058-6280-4897-a11c-fa1137398549n%40googlegroups.com</a>.<br />