Photocyclization play a key role in the deactivation mechanism of tetraphenylethylene.

Trajectory surface hopping (TSH) methods have been widely used to study photoinduced nonadiabatic processes. In the present study, nonadiabatic dynamics simulations with the widely used Tully’s fewest switches surface hopping (FSSH) algorithm and a Landau-Zener-type TSH (LZSH) algorithm have been performed for the internal conversion dynamics of pyrazine. The accuracy of the two TSH algorithms has been critically evaluated by a direct comparison with exact quantum dynamics calculations for a model of pyrazine. The model comprises the three lowest excited electronic states (B_3u(nπ^*), A_1u(nπ^*), and B_2u(ππ^*)) and the nine most relevant vibrational degrees of freedom. Considering photoexcitation to the diabatic B_2u(ππ^*) state, we examined the time-dependent diabatic and adiabatic electronic population dynamics. It is found that the diabatic populations obtained with both TSH methods are in good agreement with the exact quantum results. Fast population oscillations between the B_3u(nπ^*) and A_1u(nπ^*) states, which reflect nonadiabatic electronic transitions driven by coherent dynamics in the normal mode Q_8a, are qualitatively reproduced by both TSH methods. In addition to the model study, the TSH methods have been interfaced with the second-order algebraic diagrammatic construction ab initio electronic-structure method to perform full-dimensional on-the-fly nonadiabatic dynamics simulations for pyrazine. It is found that the electronic population dynamics obtained with the LZSH method is in excellent agreement with that obtained by the FSSH method using a local diabatization algorithm. Moreover, the electronic populations of the full-dimensional on-the-fly calculations are in excellent agreement with the populations of the three-state nine-mode model, which confirms that the internal conversion dynamics of pyrazine is accurately represented by this reduced-dimensional model on the time scale under consideration (200 fs). The original FSSH method, in which the electronic wave function is propagated in the adiabatic representation, yields less accurate results. The oscillations in the populations of the diabatic B_3u(nπ^*) and A_1u(nπ^*) states driven by the mode Q_8a are also observed in the full-dimensional dynamics simulations.

Simulating the coupled electronic and nuclear response of a molecule to light excitation requires the application of nonadiabatic molecular dynamics. However, when faced with a specific photophysical or photochemical problem, selecting the most suitable theoretical approach from the wide array of available techniques is not a trivial task. The challenge is further complicated by the lack of systematic method comparisons and rigorous testing on realistic molecular systems. This absence of comprehensive molecular benchmarks remains a major obstacle to advances within the field of nonadiabatic molecular dynamics. A CECAM workshop, Standardizing Nonadiabatic Dynamics: Towards Common Benchmarks, was held in May 2024 to address this issue. This Perspective highlights the key challenges identified during the workshop in defining molecular benchmarks for nonadiabatic dynamics. Specifically, this work outlines some preliminary observations on essential components needed for simulations and proposes a roadmap aiming to establish, as an ultimate goal, a community-driven, standardized molecular benchmark set.

The excitation wavelength dependent photodynamics of pyrrole are investigated by surface-hopping nonadiabatic dynamics simulations. The results are explained in terms of correct Rydberg–valence interaction in the lowest B₂ state.