The NRL Tight-binding method provides an efficient method for calculating properties of materials. The advantage of the NRL-TB method over classical potential simulations is that it explicitly incorporates the real electronic structure and bonding of the material, obtained by an interpolation from a database of first-principles results.
There are now nine codes used to perform these calculations:
The NRL-Tight-Binding method (TB) is a Slater-Koster-like TB method [[1]–[3]], which has been successfully applied to many single element, binary and ternary systems [[4], [5]]. The method performs very well in metals, insulators, and semiconductors [[6], [7]], and can be extended to clusters and molecules. It is built on a fitting process to Density Functional Theory results. It works in both an orthogonal and a nonorthogonal basis and computes both the energy bands and total energies. It has transferability far beyond the fitted DFT database. The NRL-TB formalism with the basic equations is outlined below. First, the total energy is determined using the fact that in DFT the total energy is given by the expression:
where is the electron density and contains the remaining parts of the DFT total energy minus the sum of the one-electron eigenvalues. In other TB approaches is a sum of pair potentials. However, based on the fact that in DFT one can uniformly shift the energy bands by a constant, in this method a shift is applied to the first-principles eigenvalues , where is the number of valence electrons in the system. This shifts the one-electron eigenvalues to new values,
and hence the total energy becomes the sum of the shifted DFT eigenvalues over the Brillouin zone, i.e.
This is a two-center TB scheme, where the on-site terms have a polynomial form as a function of the atomic density. For a single element, the density of atom i is defined as
where the sum is over all the neighboring atoms j within a range of cutoff distance of atom , is a fitting parameter, and is a smooth cutoff function. The angular-momentum-dependent on-site terms are defined by
where represents the , , and orbitals, and , , and are fitting coefficients. The two-center Slater-Koster (SK) Hamiltonian and overlap integrals consist of ten independent SK parameters, which are assumed to all have polynomial times exponential forms in terms of the neighbor distance R, given by
where indicates the type of ineractions, including , , , , , , , , , and . , , , and are our fitting coefficients. The SK overlap functions, in a non-orthogonal calculation, are defined to have the same form as the hopping parameters. Overall, there are in total 93 fitting coefficients for a single element in the on-site, hopping and overlap terms in the TB Hamiltonian with s, p, and d orbitals.
For semi-conductors, in the overlap matrix we use a polynomial of the form
where is the Kronecker delta.
In order to determine the above coefficients a least-squares procedure is used to fit to DFT total energies and energy bands as a function of volume for high symmetry structures as bcc, fcc, and simple cubic crystal structures with varying atomic volumes. In some cases, to improve the fit one uses selected volumes for the hcp structure, and for semiconductors fit the diamond lattice as well. The total energy is usually weighed at around 200-300 times over a single band energy. In general, the fitting RMS error is less than 10 mRy and 0.2 mRy for the energy bands and total energy, respectively.
Having determined the above coefficients, the method is used to predict total energies accessible to standard DFT which were not fitted, such as elastic constants, phonon spectra and surface energies. Furthermore, large scale simulations can be performed which are not practical via DFT, such as static calculations for the energetics of systems containing up to 20,000 atoms, or calculations for very large number of k-points needed for mapping Fermi surfaces and evaluating susceptibilities.
The method also has capabilities to perform molecular dynamics (MD) simulations accommodating as many as 10,000 atoms, and 10,000 MD steps, which is an impossible task for standard DFT codes. Such calculations should yield mean square displacements, thermal expansion, and vacancy formation energies [[2]–[4]].