Many chemical systems, notably those in condensed matter and in biology, are characterized by the presence of multiple low energy states which are separated by large barriers. The presence of these barriers prevents exploration of all of the configuration space during the relatively short timescales accessible in molecular dynamics (MD) simulations. Typically this means that only those configurations in a small, locally ergodic region in the vicinity of the input structure are visited.
One trick that is commonly used to resolve this problem is to add a simulation bias in order to force the interesting, rarer events to occur. This simulation bias forces the value of some particularly interesting collective variable (CV) to change more rapidly than it otherwise would. This is a very good strategy when you have a reasonable understanding of what occurs during your rare event and when you can encode what happens mathematically. For instance, to study a process such as nucleation you simply have to find a differentiable quantity that distinguishes between configurations where the ions are distributed uniformly across the simulation cell from those in which the ions are bunched together in a crystalline nucleus. To develop coordinates for complicated problems such as these we are involved with the development of a code called plumed 2.0. Already we are using this code to find new CVs that can be used to study nucleation and to study how water molecules enter the ionic channels in large protein molecules.
Researcher: Gareth Tribello
- G. A. Tribello, M Bonomi, D Branduardi, C Camilloni, G Bussi. Computer Physics Communications 185 pp 604-613 (2014)
- J Klug, C Triguero, M Del Popolo, G.A. Tribello. Journal of Physical Chemistry B (2018)
- G. A. Tribello, M. Ceriotti and M. Parrinello. Proceedings of the National Academy of Sciences, 107(41), 17509–17514 (2010).
- P. Soderhjelm, G. A. Tribello and M. Parrinello. Proceedings of the National Academy of Sciences, 109(14), 5170–5175 (2012).
- G. A. Tribello, J. Cuny, H. Eshet and M. Parrinello. Journal of Chemical Physics, 135(11), 114109 (2011).
The code we use to do these sorts of calculations can be downloaded from