Trajectory Surface Hopping with Tight Binding Density Functional Theory applied to Molecular Motors

Abstract: Non-adiabatic molecular dynamics (NAMD) has become an essential computational technique for studying the photophysical relaxation of molecular systems after light absorption. These phenomena require approximations that go beyond the Born-Oppenheimer approximation, and the accuracy of the results heavily depends on the electronic structure theory employed. Sophisticated electronic methods, however, make these techniques computationally expensive, even for medium size systems. Consequently, simulations are often performed on simplified models to interpret experimental results. In this context, a variety of techniques have been developed to perform NAMD using approximate methods, particularly Density Functional Tight Binding (DFTB). Despite the use of these techniques on large systems where ab initio methods are computationally prohibitive, a comprehensive validation has been lacking. In this work, we present a new implementation of trajectory surface hopping (TSH) combined with DFTB, utilizing non-adiabatic coupling vectors (NACVs). We selected two different systems for validation, providing an exhaustive comparison with higher-level electronic structure methods. As a case study, we simulated a system from the class of molecular motors, which has been extensively studied experimentally but remains challenging to simulate with ab initio methods due to its inherent complexity. Our approach effectively captures the key photophysical mechanism of dihedral rotation after absorption of light. Additionally, we successfully reproduce the transition from the bright to dark states observed in the time dependent fluorescence experiments, providing valuable insights into this critical part of the photophysical behavior in molecular motors.