• Citation: T. Kumagai, S. Izumi, S. Hara, and S. Sakai (2007), "Development of bond-order potentials that can reproduce the elastic constants and melting point of silicon for classical molecular dynamics simulation", Computational Materials Science, 39(2), 457-464. DOI: 10.1016/j.commatsci.2006.07.013.
    Abstract: The Tersoff potential is one of the most widely used interatomic potentials for silicon. However, its poor description of the elastic constants and melting point of diamond silicon is well known. In this research, three bond-order type interatomic potentials have been developed: the first one is fitted to the elastic constants by employing the Tersoff potential function form, the second one is fitted to both the elastic constants and melting point by employing the Tersoff potential function form and the third one is fitted to both the elastic constants and melting point by employing the modified Tersoff potential function form in which the angular-dependent term is improved. All of developed potentials well reproduce the elastic constants of diamond silicon as well as the cohesive energies and equilibrium bond lengths of silicon polytypes. The third potential can reproduce the melting point, while the second one cannot reproduce that. The elastic constants and melting point calculated using the third potential turned out to be C11 = 166.4 GPa, C12 = 65.3 GPa, C44 = 77.1 GPa and Tm = 1681 K. It was also found that only elastic constants can be reproduced using the original Tersoff potential function, and that our proposed angular-dependent term is a key to reproducing the melting point.

    Related Models:
  • LAMMPS pair_style tersoff/mod (2007--Kumagai-T--Si--LAMMPS--ipr1)
    See Computed Properties
    Notes: This file was taken from the August 22, 2018 LAMMPS distribution.
    File(s):
Date Created: October 5, 2010 | Last updated: November 20, 2024