Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. I. Algorithmic details

Aleksandar Donev, Salvatore Torquato, Frank H. Stillinger

Research output: Contribution to journalArticle

Abstract

In this first part of a series of two papers, we present in considerable detail a collision-driven molecular dynamics algorithm for a system of non-spherical particles, within a parallelepiped simulation domain, under both periodic or hard-wall boundary conditions. The algorithm extends previous event-driven molecular dynamics algorithms for spheres, and is most efficient when applied to systems of particles with relatively small aspect ratios and with small variations in size. We present a novel partial-update near-neighbor list (NNL) algorithm that is superior to previous algorithms at high densities, without compromising the correctness of the algorithm. This efficiency of the algorithm is further increased for systems of very aspherical particles by using bounding sphere complexes (BSC). These techniques will be useful in any particle-based simulation, including Monte Carlo and time-driven molecular dynamics. Additionally, we allow for a non-vanishing rate of deformation of the boundary, which can be used to model macroscopic strain and also alleviate boundary effects for small systems. In the second part of this series of papers we specialize the algorithm to systems of ellipses and ellipsoids and present performance results for our implementation, demonstrating the practical utility of the algorithm.

Original languageEnglish (US)
Pages (from-to)737-764
Number of pages28
JournalJournal of Computational Physics
Volume202
Issue number2
DOIs
StatePublished - Jan 20 2005

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. I. Algorithmic details'. Together they form a unique fingerprint.

Cite this