Algorithm optimization in molecular dynamics simulation

Di Bao Wang, Fei Bin Hsiao, Cheng Hsin Chuang, Yung-Chun Lee

Research output: Contribution to journalArticlepeer-review

9 Citations (Scopus)

Abstract

Establishing the neighbor list to efficiently calculate the inter-atomic forces consumes the majority of computation time in molecular dynamics (MD) simulation. Several algorithms have been proposed to improve the computation efficiency for short-range interaction in recent years, although an optimized numerical algorithm has not been provided. Based on a rigorous definition of Verlet radius with respect to temperature and list-updating interval in MD simulation, this paper has successfully developed an estimation formula of the computation time for each MD algorithm calculation so as to find an optimized performance for each algorithm. With the formula proposed here, the best algorithm can be chosen based on different total number of atoms, system average density and system average temperature for the MD simulation. It has been shown that the Verlet Cell-linked List (VCL) algorithm is better than other algorithms for a system with a large number of atoms. Furthermore, a generalized VCL algorithm optimized with a list-updating interval and cell-dividing number is analyzed and has been verified to reduce the computation time by 30 ∼ 60 % in a MD simulation for a two-dimensional lattice system. Due to similarity, the analysis in this study can be extended to other many-particle systems.

Original languageEnglish
Pages (from-to)551-559
Number of pages9
JournalComputer Physics Communications
Volume177
Issue number7
DOIs
Publication statusPublished - 2007 Oct 1

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Physics and Astronomy(all)

Fingerprint Dive into the research topics of 'Algorithm optimization in molecular dynamics simulation'. Together they form a unique fingerprint.

Cite this