Warning! Note that elemental potentials taken from alloy descriptions may not work well for the pure species. This is particularly true if the elements were fit for compounds instead of being optimized separately. As with all interatomic potentials, please check to make sure that the performance is adequate for your problem.
Citation: X.W. Zhou, R.E. Jones, and K. Chu (2017), "Polymorphic improvement of Stillinger-Weber potential for InGaN", Journal of Applied Physics, 122(23), 235703. DOI: 10.1063/1.5001339.
Abstract: A Stillinger-Weber potential is computationally very efficient for molecular dynamics simulations. Despite its simple mathematical form, the Stillinger-Weber potential can be easily parameterized to ensure that crystal structures with tetrahedral bond angles (e.g., diamond-cubic, zinc-blende, and wurtzite) are stable and have the lowest energy. As a result, the Stillinger-Weber potential has been widely used to study a variety of semiconductor elements and alloys. When studying an A-B binary system, however, the Stillinger-Weber potential is associated with two major drawbacks. First, it significantly overestimates the elastic constants of elements A and B, limiting its use for systems involving both compounds and elements (e.g., an A/AB multilayer). Second, it prescribes equal energy for zinc-blende and wurtzite crystals, limiting its use for compounds with large stacking fault energies. Here, we utilize the polymorphic potential style recently implemented in LAMMPS to develop a modified Stillinger-Weber potential for InGaN that overcomes these two problems.