The Anisotropic Free Energy of the Lennard-Jones Crystal-Melt Interface
We have calculated the free energy of the crystal-melt interface for the Lennard-Jones system as a function of crystal orientation, near zero pressure, by examining the roughness of the interface using molecular dynamic simulations. The anisotropy is weak, but can be accurately resolved using this approach due to the sensitivity of the fluctuations on the anisotropy. We find that the anisotropy can be described well using two parameters, based upon a low-order expansion satisfying cubic symmetry. The results are in good agreement with previous calculations of the free energies, based upon simulations used to calculate the reversible work required to create the interfaces. The weak anisotropy is also in reasonable agreement: The work here and the work of Davidchack and Laird [R. L. Davidchack and B. B. Laird, J. Chem. Phys. 118, 7651 (2003)] both predict γ100>γ110>γ111. The only discrepancy is that we find a smaller value for the difference γ100−γ111 by an amount larger than the combined error bars.
The following article appeared in Journal of Chemical Physics 116 (2003): 3920, doi:10.1063/1.1591725.