[CP2K-user] [CP2K:18543] Problems with metadynamics simulations

[Nacho Sanz Martinez]

I am trying to explore the potential energy surface in the reaction 
mechanism of an enzyme.

I am able to scan based on different collective variables, but  I am having 
some problems.

For the reaction to happen, an O of a serine has to attack a C-O-P... 
almost simultaneously the H of the Serine (OH) is transferred to an O(-1) 
of a nearby glutamic acid.

When I look at the geometries, the metadynamics seem to sweep well until 
the H from the serine passes to glutamate. When this happens, the hydrogens 
of the CH2 groups begin to pass from one group to another, compensating for 
the negative charge generated, but never giving rise to the reaction.
This has to be much more energetic than the reaction itself.

I know it's a very particular hard probability question, but can someone 
help me?

in this post you can see one of my input files

Thanks.

Nacho

-- 
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 on the web visit https://groups.google.com/d/msgid/cp2k/4a263801-ea9a-4e5f-a1c2-1cbc0d1338een%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230315/11c6f43e/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: metadyn.inp
Type: chemical/x-gamess-input
Size: 3674 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230315/11c6f43e/attachment.inp>