Loading…

Determination of the solid-liquid interfacial free energy along a coexistence line by Gibbs-Cahn integration

We calculate the solid-liquid interfacial free energy γ s l for the Lennard-Jones (LJ) system at several points along the pressure-temperature coexistence curve using molecular-dynamics simulation and Gibbs-Cahn integration. This method uses the excess interfacial energy ( e ) and stress ( τ ) along...

Full description

Saved in:
Bibliographic Details
Published in:The Journal of chemical physics 2009-09, Vol.131 (11), p.114110-114110-8
Main Authors: Laird, Brian B., Davidchack, Ruslan L., Yang, Yang, Asta, Mark
Format: Article
Language:English
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We calculate the solid-liquid interfacial free energy γ s l for the Lennard-Jones (LJ) system at several points along the pressure-temperature coexistence curve using molecular-dynamics simulation and Gibbs-Cahn integration. This method uses the excess interfacial energy ( e ) and stress ( τ ) along the coexistence curve to determine a differential equation for γ s l as a function of temperature. Given the values of γ s l for the (100), (110), and (111) LJ interfaces at the triple-point temperature ( T ∗ = k T / ϵ = 0.618 ) , previously obtained using the cleaving method by Davidchack and Laird [ J. Chem. Phys. 118 , 7657 (2003) ], this differential equation can be integrated to obtain γ s l for these interfaces at higher coexistence temperatures. Our values for γ s l calculated in this way at T ∗ = 1.0 and 1.5 are in good agreement with those determined previously by cleaving, but were obtained with significantly less computational effort than required by either the cleaving method or the capillary fluctuation method of Hoyt, Asta, and Karma [ Phys. Rev. Lett. 86 , 5530 (2001) ]. In addition, the orientational anisotropy in the excess interface energy, stress and entropy, calculated using the conventional Gibbs dividing surface, are seen to be significantly larger than the relatively small anisotropies in γ s l itself.
ISSN:0021-9606
1089-7690
DOI:10.1063/1.3231693