is it correct to use the leap frog integrator in GROMACS with WESTPA?

[Ibrahim Mohamed]

Ii am using GROMACS to perform the dynamics and usually (without using WE) i set integrator = md, tcoupl = Nose-Hoover, and pcoupl = Parrinello-Rahman. However, in the additional tutorials of WESTPA, i noticed that the mdp file contained integrator = sd.
So is it correct to use md as an integrator?

Thanks

[Lillian Chong]

Ibrahim,

WE simulation requires the use of stochastic dynamics. For MD simulations, this is typically achieved by using a stochastic thermostat, i.e. a Langevin or Bussi's stochastic velocity rescaling thermostat rather than a deterministic thermostat like Nose-Hoover. So you should use the Gromacs settings for either a Langevin or stochastic velocity rescaling thermostat, including the recommended barostat for that thermostat.

For the choice of coupling constant to use with the thermostat, I suggest that you consult this paper by Basconi & Shirts.

All the best,

Lillian

--
You received this message because you are subscribed to the Google Groups "westpa-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to westpa-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/westpa-users/4595ca35-b363-4578-a202-102b95844616n%40googlegroups.com.

--

Lillian T. Chong, Ph.D.

Professor of Chemistry

University of Pittsburgh

Pronouns: she/her/hers