Crystal-melt interface stresses: atomistic simulation calculations for a Lennard-Jones binary alloy, Stillinger-Weber Si, and embedded atom method Ni

Phys Rev E Stat Nonlin Soft Matter Phys. 2007 Jun;75(6 Pt 1):061610. doi: 10.1103/PhysRevE.75.061610. Epub 2007 Jun 29.

Abstract

Molecular-dynamics and Monte Carlo simulations have been used to compute the crystal-melt interface stress (f) in a model Lennard-Jones (LJ) binary alloy system, as well as for elemental Si and Ni modeled by many-body Stillinger-Weber and embedded-atom-method (EAM) potentials, respectively. For the LJ alloys the interface stress in the (100) orientation was found to be negative and the f vs composition behavior exhibits a slight negative deviation from linearity. For Stillinger-Weber Si, a positive interface stress was found for both (100) and (111) interfaces: f{100}=(380+/-30)mJ/m{2} and f{111}=(300+/-10)mJ/m{2}. The Si (100) and (111) interface stresses are roughly 80 and 65% of the value of the interfacial free energy (gamma) , respectively. In EAM Ni we obtained f{100}=(22+/-74)mJ/m{2}, which is an order of magnitude lower than gamma. A qualitative explanation for the trends in f is discussed.