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

9 views
Skip to first unread message

Ibrahim Mohamed

unread,
Sep 3, 2023, 4:06:54 AM9/3/23
to westpa-users
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

unread,
Sep 4, 2023, 10:54:08 AM9/4/23
to westpa...@googlegroups.com
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
Reply all
Reply to author
Forward
0 new messages