Picture of two heads

Open Access research that challenges the mind...

The Strathprints institutional repository is a digital archive of University of Strathclyde research outputs. Strathprints provides access to thousands of Open Access research papers by University of Strathclyde researchers, including those from the School of Psychological Sciences & Health - but also papers by researchers based within the Faculties of Science, Engineering, Humanities & Social Sciences, and from the Strathclyde Business School.

Discover more...

Simulating fluid-solid equilibrium with the Gibbs ensemble

Sweatman, M.B. and Quirke, N. (2004) Simulating fluid-solid equilibrium with the Gibbs ensemble. Molecular Simulation, 30 (1). pp. 23-28.

[img]
Preview
PDF (strathprints013476.pdf)
strathprints013476.pdf

Download (131kB) | Preview

Abstract

The Gibbs ensemble is employed to simulate fluid-solid equilibrium for a shifted-force Lennard-Jones system. This is achieved by generating an accurate canonical Helmholtz free-energy model of the (defect-free) solid phase. This free-energy model is easily generated, with accuracy limited only by finite-size effects, by a single isothermal-isobaric simulation at a pressure not too far from coexistence for which the chemical potential is known. We choose to illustrate this method at the known triple-point because the chemical potential is easily calculated from the coexisting gas. Alternatively, our methods can be used to locate fluid-solid coexistence and the triple-point of pure systems if the chemical potential of the solid phase can be efficiently calculated at a pressure not too far from the actual coexistence pressure. Efficient calculation of the chemical potential of solids would also enable the Gibbs ensemble simulation of bulk solid-solid equilibrium and the grand-canonical ensemble simulation of bulk solids.