Re: Rep: [CP2K:3667] reaction coordinate constrained MD

Sandeep Kumar Reddy kuma... at
Sun Dec 11 02:39:08 UTC 2011

Dear Ross,
                   To use thermodynamic integration, you use Lagrange
multipliers in put your input file.

input file should look like

      COLVAR 1
      TARGET [angstrom] 1.0
      TARGET_GROWTH [angstrom*fs^-1] 0.1
      TARGET_LIMIT [angstrom] 1.2

The average of lagrange multiplier gives the mean force acting on reaction
coordinate.   Read "Computation of the pK of liquid water using coordination
constraints" by Sprik  (Chemical Physics 258 (2000) 139,150)  for details
about this method.  This will also tell you how to remove the error
introduced during the thermodynamic integration.

Carry out the simulation at many no. of  distances.  Integration of these
forces give you the free energy.


On Sun, Dec 11, 2011 at 7:02 AM, shoutian sun
<sunshou... at>wrote:

> Dear Sandeep,
> Thanks for your reply. Actually, I want to use the
> Thermodynamic integration, and I have no experience
> on this method. Could give me some advice?
> Look forward to your reply.
> Best regards,
> Ross,sun
>  --
> You received this message because you are subscribed to the Google Groups
> "cp2k" group.
> To post to this group, send email to cp... at
> To unsubscribe from this group, send email to
> cp2k+uns... at
> For more options, visit this group at
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <>

More information about the CP2K-user mailing list