Based on the pull code below, I’m trying to determine whether the resulting force profile is physically reasonable. The force increases steadily in the first 800 ps, shows two more resistance peaks (around 2400 ps and 3000 ps), and then stabilizes around 500 kJ/mol/nm for the remaining part of the 5 ns simulation. There’s no sharp drop to zero, and the pulling distance increases smoothly throughout. Would this be considered a reasonable force profile for a loop-loop SMD simulation?
This is my first time encountering such profile but I am unsure because I have only ever done SMD for systems having more than 2 entities but here, I am trying to pull two different parts of the same entity (protein).
2 nm/ns is a very high pull rate. I would expect that your results are far from equilibrium. If I didn’t make a mistake, the integral of the plot above results in feeding 5000 kJ/mol into the system in those 5 ns. This is an enormous number.
The 0.002 (2nm/ns) is something I have reduced from 0.01 (10nm.ns), given in the supplied mdp file in http://www.mdtutorials.com/gmx/umbrella/Files/md_pull.mdp . This was because I wanted to reduce the pull rate. If I have to go down further, do you suggest going beyond 0.001 (i.e. 1 nm/ns)?
Do you also mean that the 5000 ps is a substantially high timescale for SMD and is unnecessary? So I should decrease that to 2000 ps or so?
Your issue, which is a very common issue, is that you often need to pull extremely slow with SMD to get reasonable results. You would need to pull one or more orders of magnitude slower than 0.002 nm/ps.
I don’t know what you are after. If you want a free-energy profile, umbrella sampling or AWH are better approaches. But even those might not converge within reasonable time.
For what it’s worth, just to provide some scale context, when we absolutely need to pull something, the rate is typically ~1 nm over 1 microsecond. For proper energy estimates etc, we usually do umbrella sampling, just like hess is suggesting.