[CP2K-user] [CP2K:22126] Re: CP2K constrained MD: O–H coordination number constraint does not prevent OH bond breaking

[Marcella Iannuzzi]

Hi

The fact that the distance between one Oh and one  hydrogen increases does 
not mean that the coordination number is not maintained.
You should monitor the coordination number not the distance. 
If I understand correctly, the definition of the coordination number in the 
input considers all the Oh and all the H. The distance between one specific 
Oh and one specific H is not that relevant.
Regards
Marcella

On Saturday, February 28, 2026 at 1:53:54 PM UTC+1 zqlea... at gmail.com wrote:

> Hello CP2K developers and users,
>
> I am running a constrained MD simulation in CP2K and trying to keep the 
> O–H coordination number (CN) of hydroxyl groups constrained during the 
> trajectory. However, the constraint does not seem to hold as expected.
>
> My setup is:
>
> &COLVAR
> &COORDINATION
> KINDS_FROM Oh
> KINDS_TO H
> R_0 [angstrom] 1.25
> NN 8
> ND 16
> &END COORDINATION
> &END COLVAR
>
> &CONSTRAINT
> &COLLECTIVE
> COLVAR 1
> INTERMOLECULAR
> TARGET 0.88
> &END COLLECTIVE
> &LAGRANGE_MULTIPLIERS
> COMMON_ITERATION_LEVELS 1
> &END
> &END CONSTRAINT
>
> My intention is to constrain the O–H CN so that each hydroxyl oxygen keeps 
> its single bonded H during the MD run.
>
> However, after about 15 ps of simulation, I found from the trajectory 
> analysis that the CN is not actually maintained. Around ~8 ps, some O atoms 
> in OH groups lose their only H, so their O–H CN effectively becomes 0.
>
> To diagnose this, I checked:
>
>    - 
>    
>    the O–H distances in the OH groups
>    - 
>    
>    the number of H atoms around each hydroxyl O as a function of time
>    
> Both analyses suggest that some hydroxyl O atoms indeed lose their bonded 
> H during the simulation, even though the collective variable is supposed to 
> be constrained.
>
> My questions are:
>
>    1. 
>    
>    Is this the correct way to constrain an intramolecular O–H 
>    coordination number in CP2K?
>    2. 
>    
>    Could the use of INTERMOLECULAR here be incorrect for this purpose?
>    3. 
>    
>    What should I do to correctly constrain the O–H coordination number to 
>    1? 
>    
> I will attach the figures and the input file showing: 
>
>    - 
>    
>    O–H bond distances vs time
>    - 
>    
>    the number of H atoms around each hydroxyl O vs time
>    
> Any suggestions would be greatly appreciated. Thank you very much for your 
> help.
>
> Best regards
>
> [image: OH_distances_vs_time.png]
>
> [image: O_H_count_vs_time.png]
>

-- 
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/df8058e1-71c6-430d-82a9-80ce320683ben%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20260228/d60fabb8/attachment.htm>