• Citation: R. Gröger, and J. Fikar (2026), "Variable-Charge Model of GaN Within the Streitz-Mintmire Formalism". DOI: 10.2139/ssrn.6118127.
    Abstract: An ionic-covalent interatomic potential for GaN is developed by combining the Tersoff potential with electrostatic energy, where charges are determined self-consistently with atomic positions using the charge equilibration method of Streitz and Mintmire. The electrostatic parameters of both elements are derived from their experimental ionization energies and electron affinities, whereas the rates of decay of electron densities are represented by Slater's 1s orbitals and optimized to reproduce net atomic charges obtained from DFT data using the DDEC6 charge partitioning scheme. The incorporation of variable charges resolves a long-standing error in degenerate ground state between the wurtzite and zincblende structures and correctly predicts wurtzite as the ground state. As a consequence, the model predicts finite stacking fault energy, which is shown to be in excellent agreement with DFT calculations. The potential demonstrates transferability across wurtzite, zincblende and rocksalt structures, and yields the formation energies of native point defects consistent with DFT results. Under high pressure, the model predicts a phase transition from wurtzite to rocksalt at 48 GPa, consistent with measurements based on X-ray diffraction.

    Notes: This is ionic-covalent potential with self-consistent calculation of ionic charges. The covalent part of bonds is described by the Tersoff potential [J. Nord et al., J. Phys. Cond. Mat. 15, 5649 (2003)] and the electrostatic part is treated within the Streitz-Mintmire method [F. H. Streitz and J. W. Mintmire, Phys. Rev. B 50, 11996 (1994)].

  • LAMMPS pair_style hybrid/overlay tersoff/mod coul/streitz (2026--Groger-R--Ga-N--LAMMPS--ipr1)
    See Computed Properties
    Notes: These files were provided by Roman Gröger on June 1, 2026, and match with similarly named files in the LAMMPS/potentials directory starting with LAMMPS Feature release 11 February 2026. Using this potential requires LAMMPS Feature release 11 February 2026 or newer. It cannot be used with previous releases, because its implementation required a modification of source code to ensure smooth termination of Coulomb integrals at the cut-off. An example on how to use it is provided in the example/streitz subdirectory of LAMMPS.
    "GaN.streitz" contains the electrostatic parameters. "GaN.streitz+tersoff.mod" contains the Tersoff parameters in the tersoff/mod format. "scmin.mod" provides a routine for energy minimization with self-consistent charge equilibration that should be used instead of calling minimize. "pot.in" provides an example of the LAMMPS pair_style and pair_coeff lines to use. "in.gan" provides a full example of using the potential.
    File(s):
  • Citation: A. Béré, and A. Serra (2006), "On the atomic structures, mobility and interactions of extended defects in GaN: dislocations, tilt and twin boundaries", Philosophical Magazine 86(15), 2159-2192. DOI: 10.1080/14786430600640486.
    Abstract: Results obtained by atomic computer simulation based on an adapted Stillinger–Weber (SW) potential concerning the structure and relative stability of lattice dislocations, tilt and twin boundaries in GaN are discussed. The method used for the search and description of all possible atomic configurations depends on the crystallographic structure; consequently it is of general application and the results are transferable to the wurtzite binary compounds. On the contrary, the relaxed structures and their relative energetic stability are potential dependent. The results presented here correspond to a GaN model described by a pair potential. Whenever it has been possible our results have been compared with experiments or with ab initio calculations. We present the core shape and energy of a and c crystal dislocations of both edge and screw character; [0001] tilt boundaries of misorientation angles from 9.3° (corresponding to Σ37) to θ = 44.8° (corresponding to Σ43) and (10-1n) twin boundaries (n = 1, 2, 3) [1, 2, 3, 4]. The atomic structures of the tilt boundaries can be described in terms of the three stable structures of the prism a-edge dislocation core. The (10-13) twin boundary is entirely described by 6-coordinated channels whereas the other twin boundaries present more complex structural units.

    Related Models:
  • See Computed Properties
    Notes: This file was taken from the August 22, 2018 LAMMPS distribution.
    File(s):
  • See Computed Properties
    Notes: Listing found at https://openkim.org. This KIM potential corresponds to the GaN.sw distributed with the LAMMPS package, but the parameter file format is different.
    Link(s):
  • Citation: J. Nord, K. Albe, P. Erhart, and K. Nordlund (2003), "Modelling of compound semiconductors: analytical bond-order potential for gallium, nitrogen and gallium nitride", Journal of Physics: Condensed Matter 15(32), 5649-5662. DOI: 10.1088/0953-8984/15/32/324.
    Abstract: An analytical bond-order potential for GaN is presented that describes a wide range of structural properties of GaN as well as bonding and structure of the pure constituents. For the systematic fit of the potential parameters reference data are taken from total-energy calculations within the density functional theory if not available from experiments. Although long-range interactions are not explicitly included in the potential, the present model provides a good fit to different structural geometries including defects and high-pressure phases of GaN.

    Related Models:
  • LAMMPS pair_style tersoff (2003--Nord-J--Ga-N--LAMMPS--ipr1)
    See Computed Properties
    Notes: This file was created and verified by Lucas Hale. The parameter values are comparable to those in the GaN.tersoff file in the August 22, 2018 LAMMPS distribution with this file using higher precision for the derived parameters. The parameter values are identical to the ones in the parameter file used by openKIM model MO_612061685362_001.
    File(s):
  • See Computed Properties
    Notes: Listing found at https://openkim.org. This KIM potential is based on a parameter file with identical parameter values as 2003--Nord-J--Ga-N--LAMMPS--ipr1.
    Link(s):
Date Created: October 5, 2010 | Last updated: June 04, 2026