On the calculation of the chemical potential using the particle deletion scheme

Georgios C. Boulougouris, Ioannis Economou, Doros N. Theodor Ou

Research output: Contribution to journalArticle

52 Citations (Scopus)

Abstract

A new formulation is presented for the calculation of the chemical potential from molecular simulation based on a test particle removal (inverse Widom) scheme. The new formulation introduces an intermediate stage in the calculation where the test particle to be removed is substituted by a hard particle. Chemical potential calculations at high densities from NVT and NPT Monte Carlo simulations, using a fast analytical algorithm for the computation of volume accessible to the hard core particle, are presented for the pure Lennard-Jones fluid and found to be in excellent agreement with predictions from an accurate equation of state and from simulations using the test particle insertion scheme (Widom insertion scheme). Binary mixture simulations are presented also and the new scheme is applied in the Gibbs ensemble. In all cases where the system is highly dense, the new inverse Widom scheme affords a significant reduction in CPU time compared with the widely used test particle insertion method of Widom. In addition, the new scheme is compared with a recently proposed scheme based on staged insertion.

Original languageEnglish
Pages (from-to)905-913
Number of pages9
JournalMolecular Physics
Volume96
Issue number6
Publication statusPublished - 1999
Externally publishedYes

Fingerprint

deletion
Chemical potential
insertion
simulation
Binary mixtures
Equations of state
Program processors
formulations
binary mixtures
Fluids
equations of state
fluids
predictions

ASJC Scopus subject areas

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

Cite this

On the calculation of the chemical potential using the particle deletion scheme. / Boulougouris, Georgios C.; Economou, Ioannis; Theodor Ou, Doros N.

In: Molecular Physics, Vol. 96, No. 6, 1999, p. 905-913.

Research output: Contribution to journalArticle

Boulougouris, GC, Economou, I & Theodor Ou, DN 1999, 'On the calculation of the chemical potential using the particle deletion scheme', Molecular Physics, vol. 96, no. 6, pp. 905-913.
Boulougouris, Georgios C. ; Economou, Ioannis ; Theodor Ou, Doros N. / On the calculation of the chemical potential using the particle deletion scheme. In: Molecular Physics. 1999 ; Vol. 96, No. 6. pp. 905-913.
@article{7f8aaff647aa4aca9ff2049d12b52d54,
title = "On the calculation of the chemical potential using the particle deletion scheme",
abstract = "A new formulation is presented for the calculation of the chemical potential from molecular simulation based on a test particle removal (inverse Widom) scheme. The new formulation introduces an intermediate stage in the calculation where the test particle to be removed is substituted by a hard particle. Chemical potential calculations at high densities from NVT and NPT Monte Carlo simulations, using a fast analytical algorithm for the computation of volume accessible to the hard core particle, are presented for the pure Lennard-Jones fluid and found to be in excellent agreement with predictions from an accurate equation of state and from simulations using the test particle insertion scheme (Widom insertion scheme). Binary mixture simulations are presented also and the new scheme is applied in the Gibbs ensemble. In all cases where the system is highly dense, the new inverse Widom scheme affords a significant reduction in CPU time compared with the widely used test particle insertion method of Widom. In addition, the new scheme is compared with a recently proposed scheme based on staged insertion.",
author = "Boulougouris, {Georgios C.} and Ioannis Economou and {Theodor Ou}, {Doros N.}",
year = "1999",
language = "English",
volume = "96",
pages = "905--913",
journal = "Molecular Physics",
issn = "0026-8976",
publisher = "Taylor and Francis Ltd.",
number = "6",

}

TY - JOUR

T1 - On the calculation of the chemical potential using the particle deletion scheme

AU - Boulougouris, Georgios C.

AU - Economou, Ioannis

AU - Theodor Ou, Doros N.

PY - 1999

Y1 - 1999

N2 - A new formulation is presented for the calculation of the chemical potential from molecular simulation based on a test particle removal (inverse Widom) scheme. The new formulation introduces an intermediate stage in the calculation where the test particle to be removed is substituted by a hard particle. Chemical potential calculations at high densities from NVT and NPT Monte Carlo simulations, using a fast analytical algorithm for the computation of volume accessible to the hard core particle, are presented for the pure Lennard-Jones fluid and found to be in excellent agreement with predictions from an accurate equation of state and from simulations using the test particle insertion scheme (Widom insertion scheme). Binary mixture simulations are presented also and the new scheme is applied in the Gibbs ensemble. In all cases where the system is highly dense, the new inverse Widom scheme affords a significant reduction in CPU time compared with the widely used test particle insertion method of Widom. In addition, the new scheme is compared with a recently proposed scheme based on staged insertion.

AB - A new formulation is presented for the calculation of the chemical potential from molecular simulation based on a test particle removal (inverse Widom) scheme. The new formulation introduces an intermediate stage in the calculation where the test particle to be removed is substituted by a hard particle. Chemical potential calculations at high densities from NVT and NPT Monte Carlo simulations, using a fast analytical algorithm for the computation of volume accessible to the hard core particle, are presented for the pure Lennard-Jones fluid and found to be in excellent agreement with predictions from an accurate equation of state and from simulations using the test particle insertion scheme (Widom insertion scheme). Binary mixture simulations are presented also and the new scheme is applied in the Gibbs ensemble. In all cases where the system is highly dense, the new inverse Widom scheme affords a significant reduction in CPU time compared with the widely used test particle insertion method of Widom. In addition, the new scheme is compared with a recently proposed scheme based on staged insertion.

UR - http://www.scopus.com/inward/record.url?scp=0000769793&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0000769793&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0000769793

VL - 96

SP - 905

EP - 913

JO - Molecular Physics

JF - Molecular Physics

SN - 0026-8976

IS - 6

ER -