Hard-sphere melting and crystallization with event-chain Monte Carlo

Abstract

We simulate crystallization and melting with local Monte Carlo (LMC), event-chain Monte Carlo (ECMC), and with event-driven molecular dynamics (EDMD) in systems with up to one million three-dimensional hard spheres. We illustrate that our implementations of the three algorithms rigorously coincide in their equilibrium properties. We then study nucleation in the NVE ensemble from the fcc crystal into the homogeneous liquid phase and from the liquid into the homogeneous crystal. ECMC and EDMD both approach equilibrium orders of magnitude faster than LMC. ECMC is also notably faster than EDMD, especially for the equilibration into a crystal from a disordered initial condition at high density. ECMC can be trivially implemented for hard-sphere and for soft-sphere potentials, and we suggest possible applications of this algorithm for studying jamming and the physics of glasses, as well as disordered systems.

7 pages, 7 figures

References in corpus (4)

• 2D Melting: From Liquid-Hexatic Coexistence to Continuous Transitions
• Parallelized event chain algorithm for dense hard sphere and polymer systems
• Generalized Bond Order Parameters to Characterize Transient Crystals
• Crystallization: Colloidal suspense

Cited by in corpus (27)

• Models and algorithms for the next generation of glass transition studies
• Phase-field modeling of crystal nucleation in undercooled liquids -- A review
• Colloidal Hard Spheres: Triumphs, Challenges and Mysteries
• Piecewise-Deterministic Markov Chain Monte Carlo
• Event-chain Monte Carlo: foundations, applications, and prospects
• Event-chain Monte Carlo for classical continuous spin models
• Event-chain Monte Carlo algorithms for three- and many-particle interactions
• Applicability of dynamic facilitation theory to binary hard disk systems
• Crystallization of hard spheres revisited. II. Thermodynamic modeling, nucleation work, and the surface of tension
• Event-chain algorithm for the Heisenberg model: Evidence for $z\simeq 1$ dynamic scaling
• Efficient Equilibration of Hard Spheres with Newtonian Event Chains
• Sampling efficiency of transverse forces in dense liquids
• Irreversible Monte Carlo algorithms for hard disk glasses: from event-chain to collective swaps
• Phase diagram of a system of hard cubes on the cubic lattice
• Population Annealing Simulations of a Binary Hard Sphere Mixture
• Hard core lattice gas with third next-nearest neighbor exclusion on triangular lattice: one or two phase transitions?
• High-fugacity expansion, Lee-Yang zeros and order-disorder transitions in hard-core lattice systems
• Cores for Piecewise-Deterministic Markov Processes used in Markov Chain Monte Carlo
• First Principle Computation of Random Pinning Glass Transition, Glass Cooperative Length-Scales and Numerical Comparisons
• Configuration spaces of hard spheres
• High-fugacity expansion and crystallization in non-sliding hard-core lattice particle models without a tiling constraint
• Asymptotic properties of short-range interaction functionals
• Microscopic Mechanisms of Diffusion Dynamics: A Comparative Efficiency Study of Event-Chain Monte Carlo Variants in Dense Hard Disk Systems
• Event-chain Monte Carlo simulations of the liquid to solid transition of two-dimensional decagonal colloidal quasicrystals
• Bosonized one-dimensional quantum systems through enhanced event-chain Monte Carlo
• Crystalline ordering and large fugacity expansion for hard core lattice particles
• A geometric conjecture about phase transitions