The density functional tight-binding (DFTB) formalism applied in this work has been extensively described in the literature, for a complete formulation we address to references 21-24. DFTB is based on a second-order expansion of the density functional theory (DFT) energy with respect to a reference electron density fluctuation. The resulting DFTB energy can be written as