Picture of athlete cycling

Open Access research with a real impact on health...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by Strathclyde researchers, including by researchers from the Physical Activity for Health Group based within the School of Psychological Sciences & Health. Research here seeks to better understand how and why physical activity improves health, gain a better understanding of the amount, intensity, and type of physical activity needed for health benefits, and evaluate the effect of interventions to promote physical activity.

Explore open research content by Physical Activity for Health...

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.