[CP2K-user] [CP2K:21974] Why doesn't the force eventually tend to zero when using CP2K for Slow-growth simulation

[ice ZQing]

Hello everyone,

I am currently learning how to perform slow-growth simulations in CP2K to 
obtain free-energy profiles. I followed the exercise *“NaCl Free Energy 
Dissociation”* from the CP2K website (
exercises:2015_ethz_mmm:nacl_free_energy 
<https://www.cp2k.org/exercises:2015_ethz_mmm:nacl_free_energy>), which 
uses PMF to compute the free-energy surface.

Based on that, I tried to reproduce a similar NaCl dissociation process 
using the following constraint setup:

&CONSTRAINT

    &COLLECTIVE

      COLVAR 1

      INTERMOLECULAR

      TARGET [angstrom] 2.4

      TARGET_GROWTH [angstrom*fs^-1] 1E-4

      TARGET_LIMIT [angstrom] 6.0

    &END COLLECTIVE

    &LAGRANGE_MULTIPLIERS

      COMMON_ITERATION_LEVELS 1

    &END

 &END CONSTRAINT

This means that the Na–Cl distance is slowly increased from 2.4 Å to 6.0 Å 
at a rate of 1E–4 Å/fs.
The simulation itself runs smoothly, but the results are somewhat puzzling.

As shown in the figure below, as the distance increases, the recorded 
constraint force gradually approaches a small negative value close to zero, 
instead of approaching zero from either side. Because of this, when I 
integrate the force to obtain the free-energy profile, the resulting free 
energy continues to increase in the long-distance region rather than 
converging to a plateau — which does not seem physically reasonable for a 
dissociating ion pair.

[image: shake_smooth1.png]

[image: shake_integral.png]

For reference, my system contains 1 Na⁺, 1 Cl⁻, and 63 H₂O molecules, 
exactly the same as in the CP2K exercise mentioned above. My full input 
files are attached below.

I would greatly appreciate any suggestions or insights into why the 
constraint force at large separations does not decay to zero, and how to 
correctly obtain a converged free-energy profile.

Thank you very much!

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/27df279d-fcd6-427b-8405-8c1bb9ee4a77n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20251118/519bd8ae/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NaCl_MD.inp
Type: chemical/x-gamess-input
Size: 20307 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20251118/519bd8ae/attachment-0001.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: shake_smooth1.png
Type: image/png
Size: 234702 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20251118/519bd8ae/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: shake_integral.png
Type: image/png
Size: 202165 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20251118/519bd8ae/attachment-0003.png>