YouTube video by AmsterdamDensityFunctional

Nonadiabatic molecular dynamics simulations provide a theoretical understanding of various excited-state processes in photochemistry, offering access to band widths, radiative or nonradiative relaxation and corresponding lifetimes, excited-state energies, and charge transfer. The range of method developments within the framework of time-dependent density functional theory is exceedingly large for molecular quantum chemistry. Still, it shrinks significantly when aiming to treat periodic boundary conditions. To address this gap and complement existing software packages for solid-state nonadiabatic molecular dynamics, we present an interface between the CP2K electronic structure and the NEWTON-X surface hopping codes. The interface features the generation of initial conditions, as well as adiabatic and nonadiabatic molecular dynamics, based on phenomenological or numerical time-derivative couplings. Setups are validated on gas-phase pyrazine, with electronic absorption spectra and excited-state populations for transitions between the lowest singlet states being in agreement with established molecular quantum chemistry methods. Extending the system size to crystalline pyrazine, limitations of approximate couplings are discussed, and the efficiency and applicability of the interface are demonstrated by computing broad spectra over several eV and 100 fs trajectories, considering couplings between all 80th lowest excited states, at low computational cost with a mixed semiempirical density functional theory setup.

We investigate the ability of Aufbau suppressed coupled cluster theory to act as a post-linear-response correction to widely used linear response methods for electronically excited states. We find that the theory is highly resilient to shortcomings in the underlying linear response method, with final results from less accurate starting points nearly as good as those from the best starting points. This pattern is especially stark in charge transfer states, where the approach converts starting points with multi-eV errors into post-linear-response results with errors on the order of 0.1 eV. These findings highlight the ability of Aufbau suppressed coupled cluster to perform its own orbital relaxations and raise the question of whether initializing it with an orbital relaxed reference is worth the trouble.

Modelling photochemical reactions remains a significant challenge due to the need for accurate descriptions of multiple excited states and their couplings with nuclear vibrations. Despite considerable...

We introduce a state-interaction approach for computing g-matrices within time-dependent density functional theory (TDDFT) and the Tamm–Dancoff approximation (TDA), applied here for the first time. Th...