Once again about periodic colvar

ilya ily... at gmail.com
Wed May 14 13:34:40 UTC 2008


Thanks.

I see, I have chosen a good system for studying metadynamics  :-)

On May 14, 4:42 pm, Teodoro Laino <teodor... at gmail.com> wrote:
> Ciao Ilya,
>
> Thanks again.. hopefully we should soon converge to a bug-free  
> metadynamics for periodic colvar ;-)
> There were few other missing adjustments, now in CVS.
> At least now also for extended lagrangian the ss0 and ss should be  
> both (for periodic colvar) in the
> range (-pi,pi].
>
> Check it and let us know.
> Thanks,
> Teo
>
> p.s.: of course the role of the parameters for metadynamics is very  
> important. And your attitude is the
> correct one. Keep trying!
>
> On 14 May 2008, at 09:16, ilya wrote:
>
>
>
> >  Hello!
>
> > When working with periodic colvar I noticed that while the real value
> > of colvar (instantaneous colvar value) is always printed within the
> > range (-pi,pi], but the colvar value of the extended lagrangian is not
> > corrected for the period. So for example
>
> > lagrange cv      inst val
> > 5.23293          -1.08931
> > 5.23505          -1.08060
> > or
> > 13.27589         0.70540
> > 13.27910         0.70803
>
> > Should the lagrange cv be corrected for N*2*pi?
> > The same spawned hills positions are printed to restart file.
>
> > Looking at the code of the fes program I didn't found the place where
> > the lagrange vals are corrected.
>
> > The main reason I posted this is that I still can't get the correct
> > free energy profile for butane rotation. I've played a lot with
> > parameters of metadynamics simulation (according to the article "A
> > Recipe for the Computation of the Free Energy Barrier and the Lowest
> > Free Energy Path of Concerted Reactions" and some others).
>
> > Could the above thing be the source of errors or I should continue
> > playing with params.
>
> > Thanks.
>
> > Ilya.


More information about the CP2K-user mailing list