Nonadiabatic Molecular Dynamics for Thousand Atom Systems: A Tight-Binding Approach toward PYXAID

Excited state dynamics at the nanoscale requires treatment of systems involving hundreds and thousands of atoms. In the majority of cases, depending on the process under investigation, the electronic structure component of the calculation constitutes the computation bottleneck. We developed an efficient approach for simulating nonadiabatic molecular dynamics (NA-MD) of large systems in the framework of the self-consistent charge density functional tight binding (SCC-DFTB) method. SCC-DFTB is combined with the fewest switches surface hopping (FSSH) and decoherence induced surface hopping (DISH) techniques for NA-MD. The approach is implemented within the Python extension for the ab initio dynamics (PYXAID) simulation package, which is an open source NA-MD program designed to handle nanoscale materials. The accuracy of the developed approach is tested with ab initio DFT and experimental data, by considering intraband electron and hole relaxation, and nonradiative electron-hole recombination in a CdSe quantum dot and the (10,5) semiconducting carbon nanotube. The technique is capable of treating accurately and efficiently excitation dynamics in large, realistic nanoscale materials, employing modest computational resources.