Skip to main content
Article thumbnail
Location of Repository

Direct calculation of the crystal-melt interfacial free energies for continuous potentials: Application to the Lennard-Jones system

By Ruslan L. Davidchack and Brian B. Laird


Extending to continuous potentials a cleaving wall molecular dynamics simulation method recently developed for the hard-sphere system [Phys. Rev. Lett. 85, 4751 (2000)], we calculate the crystal–melt interfacial free energies, γ, for a Lennard-Jones system as functions of both crystal orientation and temperature. At the triple point, T*=0.617, the results are consistent with an earlier cleaving potential calculation by Broughton and Gilmer [J. Chem. Phys. 84, 5759 (1986)], however, the greater precision of the current calculation allows us to accurately determine the anisotropy of γ. From our data we find that, at all temperatures studied, γ[subscript 111]<γ[subscript 110]<γ[subscript 100]. A comparison is made to the results from our previous hard-sphere calculation and to recent results for Ni by Asta, Hoyt, and Karma [Phys. Rev. B 66 100101(R) (2002)].Peer-reviewedPublisher Versio

Publisher: American Institute of Physics
Year: 2003
DOI identifier: 10.1063/1.1563248
OAI identifier:

Suggested articles


  1. (1995). The Art of Molecular Dynamics Simulation
  2. (1991). The Science of Crystallization: Microscopic Interfacial Phenomena doi
  3. (1973). The Solid–Liquid Interface doi
  4. (1986). Theory of Simple Liquids, doi
  5. This work was performed using the University of Leicester Mathematical Modelling Center’s supercomputer, which was purchased through the EPSRC strategic equipment initiative.
  6. (2002). Understanding Molecular Simulation, doi

To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.