Picture of a sphere with binary code

Making Strathclyde research discoverable to the world...

The Strathprints institutional repository is a digital archive of University of Strathclyde research outputs. It exposes Strathclyde's world leading Open Access research to many of the world's leading resource discovery tools, and from there onto the screens of researchers around the world.

Explore Strathclyde Open Access research content

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.