Ideal chemical potential contribution in molecular dynamics simulations of the grand canonical ensemble

Samantha Weerasinghe, B. Montgomery Pettitt

Research output: Contribution to journalArticle

38 Scopus citations


An extended system Hamiltonian for the grand canonical ensemble that includes the number dependence of the ideal chemical potential is investigated. Use of the ideal contribution explicitly in the equations of motion provides simpler and more stable equations of motion than previous grand molecular dynamics methods. We find the equations of motion remain quite stable even in gaseous conditions where mean field treatments of the ideal contribution provide a trivial result. The equations of motion are solved in real variable space as opposed to using virtual variables. Application of this simulation method with a Lennard-Jones fluid in the gas, fluid and solid phases is demonstrated.

Original languageEnglish (US)
Pages (from-to)897-912
Number of pages16
JournalMolecular Physics
Issue number5
StatePublished - Aug 10 1994
Externally publishedYes


ASJC Scopus subject areas

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

Cite this