Ab initio molecular dynamics (AIMD) simulation, analyzed in terms of vibrational normal modes, is a widely used technique that facilitates understanding of complex structural motions and coupling between electronic and nuclear degrees of freedom. Usually, only a subset of vibrations is directly involved in the process of interest. The impact of these vibrations can be evaluated by performing AIMD simulations by selectively freezing certain motions. Herein, we present frozen normal mode (FrozeNM), a new algorithm to apply normal-mode constraints in AIMD simulations, as implemented in the nonadiabatic excited state molecular dynamics code. We further illustrate its capacity by analyzing the impact of normal-mode constraints on the photoinduced energy transfer between polyphenylene ethynylene dendrimer building blocks. Our results show that the electronic relaxation can be significantly slowed down by freezing a well-selected small subset of active normal modes characterized by their contributions in the direction of energy transfer. The application of these constraints reduces the nonadiabatic coupling between electronic excited states during the entire dynamical simulations. Furthermore, we validate reduced dimensionality models by freezing all the vibrations, except a few active modes. Altogether, we consider FrozeNM as a useful tool that can be broadly used to underpin the role of vibrational motion in a studied process and to formulate reduced models that describe essential physical phenomena.

We utilize first-principles theory to investigate photo-induced excited-state dynamics of functionalized perylene diimide. This class of materials is highly suitable for solar energy conversion because of the strong optical absorbance, efficient energy transfer, and chemical tunability. We couple time-dependent density functional theory to a recently developed time-resolved non-adiabatic dynamics approach based on a semi-empirical description. By studying the monomer and dimer, we focus on the role stacking plays on the time-scales associated with excited-state non-radiative relaxation from a high excitonic state to the lowest energy exciton. We predict that the time-scale for energy conversion in the dimer is significantly faster than that in the monomer when equivalent excited states are accounted for. Additionally, for the dimer, the decay from the second to the nearly degenerate lowest energy excited-state involves two time-scales: a rapid decay on the order of ∼10 fs followed by a slower decay of ∼100 fs. Analysis of the spatial localization of the electronic transition density during the internal conversion process points out the existence of localized states on individual monomers, indicating that the strength of thermal fluctuations exceeds electronic couplings between the states such that the exciton hops between localized states throughout the simulation.

The role of spin-orbit coupling in the optical spectroscopy of atomic sodium isolated in solid xenon Low Temperature Physics 45, 715 (2019);