Skip to main content Accessibility help

Long-time molecular dynamics simulations on massively parallel platforms: A comparison of parallel replica dynamics and parallel trajectory splicing

  • Danny Perez (a1), Rao Huang (a1) and Arthur F. Voter (a1)


Molecular dynamics (MD) is one of the most widely used techniques in computational materials science. By providing fully resolved trajectories, it allows for a natural description of static, thermodynamic, and kinetic properties. A major hurdle that has hampered the use of MD is the fact that the timescales that can be directly simulated are very limited, even when using massively parallel computers. In this study, we compare two time-parallelization approaches, parallel replica dynamics (ParRep) and parallel trajectory splicing (ParSplice), that were specifically designed to address this issue for rare event systems by leveraging parallel computing resources. Using simulations of the relaxation of small disordered platinum nanoparticles, a comparative performance analysis of the two methods is presented. The results show that ParSplice can significantly outperform ParRep in the common case where the trajectory remains trapped for a long time within a region of configuration space but makes rapid structural transitions within this region.


Corresponding author

a)Address all correspondence to this author. e-mail:


Hide All

This author was an editor of this journal during the review and decision stage. For the JMR policy on review and publication of manuscripts authored by editors, please refer to


Permanent address: Department of Physics, Xiamen University, Xiamen 361005, China.

Contributing Editor: Vikram Gavini



Hide All
1.Lions, J.L., Maday, Y., and Turinici, G.: A “parareal” in time discretization of PDE’s. C. R. Acad. Sci. Paris Sér. I Math. 332, 661668 (2001).
2.Voter, A.F.: Hyperdynamics: Accelerated molecular dynamics of infrequent events. Phys. Rev. Lett. 78, 39083911 (1997).
3.Zamora, R.J., Uberuaga, B.P., Perez, D., and Voter, A.F.: The modern temperature-accelerated dynamics approach. Annu. Rev. Chem. Biomol. Eng. 7, 87110 (2016).
4.Sorensen, M.R. and Voter, A.F.: Temperature-accelerated dynamics for simulation of infrequent events. J. Chem. Phys. 112, 95999606 (2000).
5.Henkelman, G. and Jonsson, H.: Long time scale kinetic Monte Carlo simulations without lattice approximation and predefined event table. J. Chem. Phys. 115, 96579666 (2001).
6.Perez, D., Uberuaga, B.P., and Voter, A.F.: The parallel replica dynamics method—Coming of age. Comput. Mater. Sci. 100, 90103 (2015).
7.Voter, A.F.: Parallel replica method for dynamics of infrequent events. Phys. Rev. B 57, 1398513988 (1998).
8.Perez, D., Cubuk, E.D., Waterland, A., Kaxiras, E., and Voter, A.F.: Long-time dynamics through parallel trajectory splicing. J. Chem. Theory Comput. 12, 1828 (2016).
9.Le Bris, C., Lelievre, T., Luskin, M., and Perez, D.: A mathematical formalization of the parallel replica dynamics. Monte Carlo Methods Appl. 18, 119146 (2012).
10.El-Mellouhi, F., Mousseau, N., and Lewis, L.J.: Kinetic activation-relaxation technique: An off-lattice self-learning kinetic Monte Carlo algorithm. Phys. Rev. B 78, 153202 (2008).
11.Fichthorn, K.A. and Lin, Y.: A local superbasin kinetic Monte Carlo method. J. Chem. Phys. 138, 164104 (2013).
12.Chatterjee, A. and Voter, A.F.: Accurate acceleration of kinetic Monte Carlo simulations through the modification of rate constants. J. Chem. Phys. 132, 194101 (2010).
13.Novotny, M.A.: Monte Carlo algorithms with absorbing Markov chains: Fast local algorithms for slow dynamics. Phys. Rev. Lett. 74, 1 (1995).
14.Wales, D.: Energy Landscapes: Applications to Clusters, Biomolecules and Glasses (Cambridge University Press, Cambridge, United Kingdom, 2003).
15.Available at:
16.Plimpton, S.: Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 119 (1995).
17.Voter, A.F. and Chen, S.P.: Accurate interatomic potentials for Ni, Al and Ni3Al. MRS Online Proc. Libr. 82, 175180, (1986).
18.Voter, A.F.: Embedded Atom Method Potentials for Seven FCC Metals: Ni, Pd, Pt, Cu, Ag, Au, and Al; Los Alamos Unclassified Technical Report, LA-UR-93–3901 (1993).
19.Schenk, T., Holland-Moritz, D., Simonet, V., Bellissent, R., and Herlach, D.M.: Icosahedral short-range order in deeply undercooled metallic melts. Phys. Rev. Lett. 89, 075507 (2002).
20.Reichert, H., Klein, O., Dosch, H., and Donk, M.: Observation of five-fold local symmetry in liquid lead. Nature 408, 839 (2000).
21.Deuflhard, P. and Weber, M.: Robust Perron cluster analysis in conformation dynamics. Linear Algebra Appl. 398, 161184 (2005).
22.Huang, R., Lo, L-T., Wen, Y., Voter, A.F., and Perez, D.: Cluster analysis of accelerated molecular dynamics simulations: A case study of the decahedron to icosahedron transition in Pt nanoparticles. J. Chem. Phys. 147, 152717 (2017).
23.Clarke, A.S. and Jónsson, H.: Structural changes accompanying densification of random hard-sphere packings. Phys. Rev. E 47, 3975 (1993).
24.Apra, E., Baletto, F., Ferrando, R., and Fortunelli, A.: Amorphization mechanism of icosahedral metal nanoclusters. Phys. Rev. Lett. 93, 065502 (2004).
25.Doye, J.P.K., Miller, M.A., and Wales, D.J.: The double-funnel energy landscape of the 38-atom Lennard-Jones cluster. J. Chem. Phys. 110, 68966906 (1999).


Related content

Powered by UNSILO

Long-time molecular dynamics simulations on massively parallel platforms: A comparison of parallel replica dynamics and parallel trajectory splicing

  • Danny Perez (a1), Rao Huang (a1) and Arthur F. Voter (a1)


Altmetric attention score

Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed.