RE: Maximum Timestep for NAMD

From: Vermaas, Joshua (Joshua.Vermaas_at_nrel.gov)
Date: Thu Feb 07 2019 - 15:02:20 CST

Yup. The 1-2 fs for atomic scale simulation comes from the fastest vibrational modes in the system. Bonds to hydrogen have oscillatory periods on the order of 10fs, with heavy atom bonds having periods that are only somewhat larger. In order to have your system not routinely blow itself apart, you need to keep your timestep somewhere between 8-10 times smaller than the vibrational periods in your system. This limits classical atomic simulations to 1fs. Unless you fix bond-lengths to hydrogen with SHAKE/RATTLE/LINCS, which will get you to 2fs. Or you can also restrain heavy-atom bond lengths and some angles, and get up to 4fs with virtual sites in gromacs. Coarse grained systems allow for larger timesteps because the particle masses are larger, and 20 fs timesteps aren't unheard of. Basically, if we could take 100fs timesteps without a loss in accuracy, we would already be doing it routinely. :)

-Josh

On 2019-02-07 13:49:27-07:00 owner-namd-l_at_ks.uiuc.edu wrote:
Hello,

The first thing to consider is why the 1-2 fs timestep is most often used.
You might also consider various accelerated MD methods and coarse grained MD.

Brian

From: owner-namd-l_at_ks.uiuc.edu <owner-namd-l_at_ks.uiuc.edu> On Behalf Of Alex Saad-Falcon
Sent: Thursday, February 7, 2019 12:25 PM
To: namd-l_at_ks.uiuc.edu
Subject: namd-l: Maximum Timestep for NAMD

Hi all,

What is the maximum acceptable timestep for NAMD? I would like to run longer (1ns to 1ms) simulations, and I am wondering how much I can increase the timestep (to 10fs or 100fs) without detracting terribly from accuracy. What do I need to consider?

Thanks,

Alex

This archive was generated by hypermail 2.1.6 : Tue Dec 31 2019 - 23:20:29 CST