Picture of neon light reading 'Open'

Discover open research at Strathprints as part of International Open Access Week!

23-29 October 2017 is International Open Access Week. The Strathprints institutional repository is a digital archive of Open Access research outputs, all produced by University of Strathclyde researchers.

Explore recent world leading Open Access research content this Open Access Week from across Strathclyde's many research active faculties: Engineering, Science, Humanities, Arts & Social Sciences and Strathclyde Business School.

Explore all Strathclyde Open Access research outputs...

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 - Accepted Author Manuscript

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.