DFTB+ calculation#