Newtonian Event-Chain Monte Carlo and Collision Prediction with Polyhedral Particles
arXiv:2104.06829 · doi:10.1021/jacs.1c07281
Abstract
Polyhedral nanocrystals are building blocks for nanostructured materials that find applications in catalysis and plasmonics. Synthesis efforts and self-assembly experiments have been assisted by computer simulations that predict phase equilibra. Most current simulations employ Monte Carlo methods, which generate stochastic dynamics. Collective and correlated configuration updates are alternatives that promise higher computational efficiency and generate trajectories with realistic dynamics. One such alternative involves event-chain updates and has recently been proposed for spherical particles. In this contribution, we develop and apply event-chain Monte Carlo for hard convex polyhedra. Our simulation makes use of an improved computational geometry algorithm XenoSweep, which predicts sweep collision in a particularly simple way. We implement Newtonian event chains in the open source general-purpose particle simulation toolkit HOOMD-blue for serial and parallel simulation. The speed-up over state-of-the-art Monte Carlo is between a factor of 10 for nearly spherical polyhedra and a factor of 2 for highly aspherical polyhedra. Finally, we validate the Newtonian event-chain algorithm by applying it to a current research problem, the multi-step nucleation of two classes of hard polyhedra.
arXiv comment12 pages, 13 figures
References in corpus (10)
- Disordered, Quasicrystalline and Crystalline Phases of Densely Packed Tetrahedra
- Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods
- Event-chain Monte Carlo algorithms for hard-sphere systems
- 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
- Scalable Metropolis Monte Carlo for simulation of hard shapes
- Massively parallel Monte Carlo for many-particle simulations on GPUs
- Ligand-Mediated Interactions between Nanoscale Surfaces Depend Sensitively and Nonlinearly on Temperature, Facet Dimensions, and Ligand Coverage
- Parallelized event chain algorithm for dense hard sphere and polymer systems
- Multithreaded event-chain Monte Carlo with local times