A conservative and a hybrid early rejection schemes for accelerating Monte Carlo molecular simulation

Ahmad Salim Kadoura, Amgad Salama, Shuyu Sun

Research output: Contribution to journalArticlepeer-review

6 Scopus citations

Abstract

Molecular simulation could provide detailed description of fluid systems when compared to experimental techniques. They can also replace equations of state; however, molecular simulation usually costs considerable computational efforts. Several techniques have been developed to overcome such high computational costs. In this paper, two early rejection schemes, a conservative and a hybrid one, are introduced. In these two methods, undesired configurations generated by the Monte Carlo trials are rejected earlier than it would when using conventional algorithms. The methods are tested for structureless single-component Lennard-Jones particles in both canonical and NVT-Gibbs ensembles. The computational time reduction for both ensembles is observed at a wide range of thermodynamic conditions. Results show that computational time savings are directly proportional to the rejection rate of Monte Carlo trials. The proposed conservative scheme has shown to be successful in saving up to 40% of the computational time in the canonical ensemble and up to 30% in the NVT-Gibbs ensemble when compared to standard algorithms. In addition, it preserves the exact Markov chains produced by the Metropolis scheme. Further enhancement for NVT-Gibbs ensemble is achieved by combining this technique with the bond formation early rejection one. The hybrid method achieves more than 50% saving of the central processing unit (CPU) time.
Original languageEnglish (US)
Pages (from-to)2575-2586
Number of pages12
JournalMolecular Physics
Volume112
Issue number19
DOIs
StatePublished - Mar 17 2014

ASJC Scopus subject areas

  • Molecular Biology
  • Physical and Theoretical Chemistry
  • Biophysics
  • Condensed Matter Physics

Fingerprint Dive into the research topics of 'A conservative and a hybrid early rejection schemes for accelerating Monte Carlo molecular simulation'. Together they form a unique fingerprint.

Cite this