Re: smd force is too high

From: Yi Wang, PhD (dexterwy_at_gmail.com)
Date: Sat Mar 16 2013 - 09:35:44 CDT

Jackie,

There was one time I got a sudden increase of force: I restarted the simulation from a snapshot where the ligand is entering the box from the other side because of PBC. Have you checked the structures at those timesteps where you got ultrahigh force?
___________________________

Yi (Yves) Wang, PhD
Chinese Academy of Science

On 2013-3-16, at 上午1:35, Jacqueline Schmidt wrote:

> Dear Aron,
> I just learnt the units of force printed by NAMD are PN not kcal/mol/Ang. Converting the units, I actually get something close to -0.6...I appreciate your comments in any case
>
> This amount is just the amount of work done for pulling the peptide for ~0.1 or less angstrom away from the surface. I expect the overall work or free energy ends up to less than 10-20 kcal/mol...
>
> jackie
>
> From: Aron Broom <broomsday_at_gmail.com>
> To: Jacqueline Schmidt <jacqueline_schmidt_1983_at_yahoo.com>
> Cc: namd <namd-l_at_ks.uiuc.edu>
> Sent: Friday, March 15, 2013 1:10 PM
> Subject: Re: namd-l: smd force is too high
>
> Hi Jackie,
>
> I'm not really very familiar with SMD, but have you tried a smaller pulling velocity? I can imagine that if the pulling velocity is too fast to allow the system to be in a constant pseudo-equilibrium, you might appear to have a steeper energy barrier to pulling the protein because you have slowly relaxing orthogonal degrees of freedom that aren't being given time to relax. That is, you're crossing the energy barriers associated with those slowly relaxing degrees of freedom using a path that isn't the lowest energy one. Again, I'm not terribly familiar with SMD, but this kind of problem is fairly general to any biased MD and it will generally present itself as the energy barrier appearing to be much greater than you expected.
>
> The above being said, is the value you are getting really that extreme? I suppose -0.6 kcal/mol/Ang is effectively the gradient for the binding energy? If so, what would you be expecting for the actual binding energy over the distance you want to pull? It seems like a gradient of this magnitude would only yield a binding energy of at most a few kcal/mol or am I confused?
>
> ~Aron
>
>
> On Fri, Mar 15, 2013 at 12:25 PM, Jacqueline Schmidt <jacqueline_schmidt_1983_at_yahoo.com> wrote:
> Dear NAMD users,
>
> I have a question regarding the smd simulation using NAMD. I am pulling off a peptide from a surface (surface is fixed) using smd with force constant of 6 kcal/mol/Ang^2 and pulling velocity of 5x10^-6 Ang/fs. My problem is the forces that are printed in the output are unreasonably highe!
>
> e.g. at step 100, the center of mass of peptide has been moved ~0.1 Ang in z-direction. The smd force printed in the output is fz~-43 and I assume the units are kcal/mol/Ang
>
> but this is so high! Doing a rough calculation for force Fz=k(vt-z+z0)= 6(5*10^-6 -0.1) ~-0.6 kcal/mol/Ang
>
> can anyone please help me with this?
>
> I will appreciate it
> jackie
>
>
>
> --
> Aron Broom M.Sc
> PhD Student
> Department of Chemistry
> University of Waterloo
>
>

This archive was generated by hypermail 2.1.6 : Wed Dec 31 2014 - 23:21:02 CST