Density functional tight binding

This paper reviews the basic principles of the density-functional tight-binding (DFTB) method, which is based on density-functional theory as formulated by Hohenberg, Kohn and Sham (KS-DFT). DFTB consists of a series of models that are derived from a Taylor series expansion of the KS-DFT total energy. In the lowest order (DFTB1), densities and potentials are written as superpositions of atomic densities and potentials. The Kohn-Sham orbitals are then expanded to a set of localized atom-centred functions, which are obtained for spherical symmetric spin-unpolarized neutral atoms self-consistently. The whole Hamilton and overlap matrices contain one-and two-centre contributions only. Therefore, they can be calculated and tabulated in advance as functions of the distance between atomic pairs. The second contributions to DFTB1, the DFT double counting terms, are summarized together with nuclear repulsion energy terms and can be rewritten as the sum of pairwise repulsive terms. The second-order (DFTB2) and third-order (DFTB3) terms in the energy expansion correspond to a self-consistent representation, where the deviation of the ground-state density from the reference density is represented by charge monopoles only. This leads to a computationally efficient representation in terms of atomic charges (Mulliken), chemical hardness (Hubbard) parameters and scaled Coulomb laws. Therefore, no additional adjustable parameters enter the DFTB2 and DFTB3 formalism. The handling of parameters, the efficiency, the performance and extensions of DFTB are briefly discussed.