Re: Colvar vs. SMD

From: Giacomo Fiorin (giacomo.fiorin_at_gmail.com)
Date: Wed Mar 16 2022 - 10:22:06 CDT

Hi Steve, I was about to send this reply before seeing Jérôme's. The
configuration is truncated, and I don't see in particular the part where
the restraint center is moved. You should have "target*" options in the
harmonic restraint config, as shown in this example:

https://urldefense.com/v3/__https://colvars.github.io/colvars-refman-namd/colvars-refman-namd.html*sec:colvars_crash_course__;Iw!!DZ3fjg!uevFQk1KZgzAwvjrItaW5uK9soc1rPU04OOKytmp1ayNO61JTN0y0Jy8tBAaVNRkCQ$

You could also enable outputCenters and outputAccumulatedWork
<https://urldefense.com/v3/__https://colvars.github.io/colvars-refman-namd/colvars-refman-namd.html*sec:colvarbias_smd__;Iw!!DZ3fjg!uevFQk1KZgzAwvjrItaW5uK9soc1rPU04OOKytmp1ayNO61JTN0y0Jy8tBDHisCbdw$ >
to get a rough "PMF" (note the quotes: it's still SMD!), to process via
your own plotting tools or this convenience script
<https://urldefense.com/v3/__https://github.com/Colvars/colvars/blob/master/colvartools/plot_colvars_traj.py?raw=true__;!!DZ3fjg!uevFQk1KZgzAwvjrItaW5uK9soc1rPU04OOKytmp1ayNO61JTN0y0Jy8tBAZoETPMw$ >

This archive was generated by hypermail 2.1.6 : Tue Dec 13 2022 - 14:32:44 CST