**From:** sruthi c k (*cksruthikvr_at_gmail.com*)

**Date:** Fri Jan 08 2016 - 06:14:08 CST

**Next message:**Downs, Andrew S CTR USARMY ARL (US): "Compilation error across Intel 14,15 and 16 (UNCLASSIFIED)"**Previous message:**sruthi c k: "Fwd: Regarding system forces in NAMD"**In reply to:**Giacomo Fiorin: "Re: Regarding system forces in NAMD"**Next in thread:**sruthi c k: "Re: Regarding system forces in NAMD"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Hi Jerome and Giacomo,

Thanks for your clarifications.

I think my question is not clear enough, so let me clarify it.

1) I did 200ns long UNBIASED NPT simulation.

2) We do not do ABF. We use COLVAR module and system forces to postprocess

the forces for another purpose. But before doing so, we ran into confusions

discussed below.

3) Since we do not use ABF, HideJacobian option was NOT activated and from

the manual I understand that its default value is "no"

4) Hence I assume the system forces output includes Jacobian forces also.

5) I then calculated the histogram. The range of Rg was divided into 30

small bins and the number of Rg appeared in each bin was counted to get the

histogram.

6) The forces corresponding to Rg values falling in each bin was averaged

to get the <F_system(Rg)>

7) Regarding convergence of the simulation: The average forces look

*converged* when I calculate it using data obtained from the 150ns

simulation and compare it with average forces obtained in 200ns simulation.

*After following this procedure, this is the problem I am facing:A1)*. PMF1

defined as -kbT log(histogram) and that the PMF2 defined by -(integral of

<F_system(Rg)>) are NOT MATCHING.

I am attaching the plot of histogram corresponding to PMF1 and the *<F_system>

*corresponding to PMF2 with this mail.

I have another question related to system force which I had NOT mentioned

in the previous mail.

*B1)* I tried to calculate the system forces by taking the numerical

derivative of colvar velocities appearing in the traj file.

*B2)* Later I found that I am missing a term "m_ksi" as defined in Eqn (7)

in the paper "Adaptive biasing force method for scalar and vector free

energy calculations" by Eric Darve et al.

m_ksi is defined as

m_ksi = SUM((1/m_k)* (derivative of ksi wrt cartesian

coordinates)^2)

*B3)* When including m_ksi in my calculations, I was not sure how to handle

the derivatives relative to all possible cartesian coordinates, as some of

them are related by the covalent bond lengths in the system (rigid bonds

option "ALL" was used in the namd input file)

*B4)* I have read in the paper titled "The Adaptive Biasing Force Method:

Everything You Always Wanted To Know but Were Afraid To Ask"

by Jeffrey Comer et al. that "If only the heavy atom of the bonded pair is

involved in the transition coordinate, then the adaptive biasing force

algorithm will include the constraint force emanating from the hydrogen in

addition to the thermodynamic force, thus contaminating its estimate. A

straightforward solution is to include both the hydrogen and its parent

atom in the collective variable(s) defining the transition coordinate,

causing the constraint forces to cancel each other"

*B5)* My question is that since I calculate Radius of gyration even

including hydrogen atoms can I be not worried about the constraint forces?

This question also came up because the PMF in Question *A1* above was not

coming out as expected

Thank you very much,

Sruthi

On Thu, Jan 7, 2016 at 8:48 PM, Giacomo Fiorin <giacomo.fiorin_at_gmail.com>

wrote:

*> If you're using a radius of gyration, the Jacobian force is (3N-4)/r,
*

*> where r is the radius of gyration, so it's easy to add or subtract while
*

*> plotting.
*

*>
*

*> On Thu, Jan 7, 2016 at 10:13 AM, sruthi c k <cksruthikvr_at_gmail.com> wrote:
*

*>
*

*>> Hi Jerome,
*

*>> Thanks for your quick reply.
*

*>> I am looking at radius of gyration of protein. Hide Jacobian is not
*

*>> activated. Can we have that option even if we are not doing ABF? I thought
*

*>> it is an option given for ABF. My simulation is UNBIASED.
*

*>> Thanks
*

*>> Sruthi
*

*>>
*

*>> On Thu, Jan 7, 2016 at 8:07 PM, Jérôme Hénin <jerome.henin_at_ibpc.fr>
*

*>> wrote:
*

*>>
*

*>> > Hi Sruthi,
*

*>> >
*

*>> > If you are looking at a radius of gyration, keep in mind that that
*

*>> > coordinate experiences a strong Jacobian force especially at small
*

*>> values
*

*>>
*

*>
*

*>
*

*>
*

*> --
*

*> Giacomo Fiorin
*

*> Assistant Professor of Research
*

*> Institute for Computational Molecular Science (ICMS)
*

*> College of Science and Technology, Temple University
*

*> 1925 North 12th Street (035-07), Room 704D
*

*> Philadelphia, PA 19122-1801
*

*> Phone: +1-215-204-4213
*

*>
*

*> Scholar: http://goo.gl/Q3TBQU
*

*> Personal: http://giacomofiorin.github.io/
*

*> Lab page: https://icms.cst.temple.edu/members.html
*

*>
*

*>
*

*>
*

- application/postscript attachment: Histogram.eps

- application/postscript attachment: Av_system_force.eps

**Next message:**Downs, Andrew S CTR USARMY ARL (US): "Compilation error across Intel 14,15 and 16 (UNCLASSIFIED)"**Previous message:**sruthi c k: "Fwd: Regarding system forces in NAMD"**In reply to:**Giacomo Fiorin: "Re: Regarding system forces in NAMD"**Next in thread:**sruthi c k: "Re: Regarding system forces in NAMD"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Tue Dec 27 2016 - 23:21:43 CST
*