Efficient Equilibration of Hard Spheres with Newtonian Event Chains
arXiv:1903.11392 · doi:10.1063/1.5090882
Abstract
An important task in the simulation of hard spheres and other hard particles is structure prediction via equilibration. Event-driven molecular dynamics is efficient because its Newtonian dynamics equilibrates fluctuations with the speed of sound. Monte Carlo simulation is efficient if performed with correlated position updates in event chains. Here, we combine the core concepts of molecular dynamics and event chains into a new algorithm involving Newtonian event chains. Measurements of the diffusion coefficient, nucleation rate, and melting speed demonstrate that Newtonian event chains outperform other algorithms. Newtonian event chains scale well to large systems and can be extended to anisotropic hard particles without approximations.
arXiv comment8 pages, 6 figures, 3 tables
References in corpus (17)
- Strong scaling of general-purpose molecular dynamics simulations on GPUs
- Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods
- Crystal nucleation of hard spheres using molecular dynamics, umbrella sampling and forward flux sampling: A comparison of simulation techniques
- Strict Detailed Balance is Unnecessary in Monte Carlo Simulation
- Event-chain Monte Carlo algorithms for hard-sphere systems
- Avoiding unphysical kinetic traps in Monte Carlo simulations of strongly attractive particles
- Phase diagram of colloidal hard superballs: from cubes via spheres to octahedra
- Generalized event-chain Monte Carlo: Constructing rejection-free global-balance algorithms from infinitesimal steps
- Rejection-free Geometric Cluster Algorithm for Complex Fluids
- Rejection-free Monte-Carlo sampling for general potentials
- A parameter-free, solid-angle based, nearest-neighbor algorithm
- Scalable Metropolis Monte Carlo for simulation of hard shapes
- Simple and efficent algorithm for large scale molecular dynamics simulation in hard disk systems
- Complex Crystals from Size-disperse Spheres
- Hard-sphere melting and crystallization with event-chain Monte Carlo
- Event-chain Monte Carlo algorithms for three- and many-particle interactions
- Parallelized event chain algorithm for dense hard sphere and polymer systems
Cited by in corpus (19)
- Colloidal Hard Spheres: Triumphs, Challenges and Mysteries
- Event-chain Monte Carlo: foundations, applications, and prospects
- Spontaneous Crystallization in Systems of Binary Hard Sphere Colloids
- Hard-disk dipoles and non-reversible Markov chains
- PDMP characterisation of event-chain Monte Carlo algorithms for particle systems
- Molecular simulation from modern statistics: Continuous-time, continuous-space, exact
- Large-scale dynamics of event-chain Monte Carlo
- Sparse hard-disk packings and local Markov chains
- Lifted TASEP: a Bethe ansatz integrable paradigm for non-reversible Markov chains
- Direction-sweep Markov chains
- Sampling algorithms in statistical physics: a guide for statistics and machine learning
- Fast event-driven simulations for soft spheres: from dynamics to Laves phase nucleation
- Velocity trapping in the lifted TASEP and the true self-avoiding random walk
- Microscopic Mechanisms of Diffusion Dynamics: A Comparative Efficiency Study of Event-Chain Monte Carlo Variants in Dense Hard Disk Systems
- Markov-chain sampling for long-range systems without evaluating the energy
- Generalization of the Wall Theorem to Out-of-equilibrium Conditions
- Kinetic Event-Chain Algorithm for Active Matter
- Necessary and sufficient symmetries in Event-Chain Monte Carlo with generalized flows and Application to hard dimers
- Lifted TASEP: long-time dynamics,generalizations, and continuum limit