Multiple timestep integrator and nonbondedCutoff

The functionality of OpenMM will (eventually) include everything that one would need to run modern molecular simulation.
POST REPLY
User avatar
Aron Broom
Posts: 54
Joined: Tue Mar 13, 2012 11:33 am

Multiple timestep integrator and nonbondedCutoff

Post by Aron Broom » Sat Nov 25, 2017 11:31 pm

Hi,

I'm wondering if it is possible to use the Multiple timestep integrator to evaluate forces outside the nonbondedCutoff less frequently than the other forces?

Since the identity of forces outside this cutoff would change over time, it seems the forcegroups would have to be updated every X timesteps or something. Maybe this wouldn't be worth it computationally?

Thanks,

Aron

User avatar
Peter Eastman
Posts: 2541
Joined: Thu Aug 09, 2007 1:25 pm

Re: Multiple timestep integrator and nonbondedCutoff

Post by Peter Eastman » Sun Nov 26, 2017 10:38 am

You can create two different NonbondedForces, one with a short cutoff and one with a longer cutoff. The inner timesteps would use Fshort, and the outer timesteps would use Flong-Fshort. MTSIntegrator doesn't support that directly, since there's no way to tell it to use the difference between two forces. But it's just a CustomIntegrator, so you can easily create your own version that does that.

Peter

POST REPLY