ABSTRACT
We present an open-source MLatom@XACS software ecosystem for on-the-fly surface hopping nonadiabatic dynamics based on the Landau-Zener-Belyaev-Lebedev algorithm. The dynamics can be performed via Python API with a wide range of quantum mechanical (QM) and machine learning (ML) methods, including ab initio QM (CASSCF and ADC(2)), semiempirical QM methods (e.g., AM1, PM3, OMx, and ODMx), and many types of ML potentials (e.g., KREG, ANI, and MACE). Combinations of QM and ML methods can also be used. While the user can build their own combinations, we provide AIQM1, which is based on Δ-learning and can be used out-of-the-box. We showcase how AIQM1 reproduces the isomerization quantum yield of trans-azobenzene at a low cost. We provide example scripts that, in dozens of lines, enable the user to obtain the final population plots by simply providing the initial geometry of a molecule. Thus, those scripts perform geometry optimization, normal mode calculations, initial condition sampling, parallel trajectories propagation, population analysis, and final result plotting. Given the capabilities of MLatom to be used for training different ML models, this ecosystem can be seamlessly integrated into the protocols building ML models for nonadiabatic dynamics. In the future, a deeper and more efficient integration of MLatom with Newton-X will enable a vast range of functionalities for surface hopping dynamics, such as fewest-switches surface hopping, to facilitate similar workflows via the Python API.
ABSTRACT
Diarylethene (DAE) molecular switches have continued to attract the attention of researchers for over 20 years. Their remarkable photophysical properties endow them with countless applications in photonics and molecular technologies. However, despite extensive experimental and theoretical research, the mechanism of DAE photoswitching is not yet fully rationalized. In this work, we investigate the ring closure dynamics of a popular DAE switch, 1,2-bis(3-methyl-5-phenyl-2 thienyl)perfluorocyclopentene (PT), using nonadiabatic molecular dynamics (NAMD) simulations. Employing the fewest switches surface hopping protocol, along with the semi-empirical multireference ODM2/MRCI-SD method, we investigate possible reaction pathways for this photoprocess, as well as their timescales and resulting photoproducts. Furthermore, using a dynamic configuration-space sampling procedure, we elucidate the role of triplet states in the photocyclization of PT, supporting available experimental data for the closely related DMPT molecule, which indicate an ultrafast intersystem crossing (ISC) transition competing with the singlet-driven photoswitching reaction. Our findings not only corroborate experimental studies on DAE switches, but also provide new mechanistic insights into the potential use in the rational design of DAE switches tailored for specific technological applications.
