Fast event-driven simulations for soft spheres: from dynamics to Laves phase nucleation
arXiv:2403.12755 · doi:10.1063/5.0209178
Abstract
Conventional molecular dynamics (MD) simulations struggle when simulating particles with steeply varying interaction potentials, due to the need to use a very short time step. Here, we demonstrate that an event-driven Monte Carlo (EDMC) approach first introduced by Peters and de With [Phys. Rev. E 85, 026703 (2012)] represents an excellent substitute for MD in the canonical ensemble. In addition to correctly reproducing the static thermodynamic properties of the system, the EDMC method closely mimics the dynamics of systems of particles interacting via the steeply repulsive Weeks-Chandler-Andersen (WCA) potential. In comparison to time-driven MD simulations, EDMC runs faster by over an order of magnitude at sufficiently low temperatures. Moreover, the lack of a finite time step in EDMC circumvents the need to trade accuracy against simulation speed associated with the choice of time step in MD. We showcase the usefulness of this model to explore the phase behavior of the WCA model at extremely low temperatures, and to demonstrate that spontaneous nucleation and growth of the Laves phases is possible at temperatures significantly lower than previously reported.
SI is included as a separate file
References in corpus (18)
- Accurate determination of crystal structures based on averaged local bond order parameters
- Event-chain Monte Carlo algorithms for hard-sphere systems
- DynamO: A free O(N) general event-driven molecular-dynamics simulator
- Generalized event-chain Monte Carlo: Constructing rejection-free global-balance algorithms from infinitesimal steps
- Rejection-free Monte-Carlo sampling for general potentials
- Event-Driven Molecular Dynamics in Parallel
- Tetrahedrality dictates dynamics in hard spheres
- Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids
- Cell-veto Monte Carlo algorithm for long-range systems
- Efficient Equilibration of Hard Spheres with Newtonian Event Chains
- Quasi-crystalline order in vibrated granular matter
- Discontinuous Molecular Dynamics for Rigid Bodies: Applications
- Structure and stability of helices in square-well homopolymers
- Comparing four hard-sphere approximations for the low-temperature WCA melting line
- A simple and accurate method to determine fluid-crystal phase boundaries from direct coexistence simulations
- Concepts in Monte Carlo sampling
- Mapping continuous potentials to discrete forms
- Optimizing Event-Driven Simulations