Exploring the Global Reaction Coordinate for Retinal Photoisomerization: A Graph Theory-Based Machine Learning Approach

Unraveling the reaction pathway of photoinduced reactions poses a great challenge owing to its complexity. Recently, graph theory-based machine learning combined with nonadiabatic molecular dynamics (NAMD) has been applied to obtain the global reaction coordinate of the photoisomerization of azobenzene. However, NAMD simulations are computationally expensive as they require calculating the nonadiabatic coupling vectors at each time step. Here, we showed that ab initio molecular dynamics (AIMD) can be used as an alternative to NAMD by choosing an appropriate initial condition for the simulation. We applied our methodology to determine a plausible global reaction coordinate of retinal photoisomerization, which is essential for human vision. On rank-ordering the internal coordinates, based on the mutual information (MI) between the internal coordinates and the HOMO energy, NAMD and AIMD give a similar trend. Our results demonstrate that our AIMD-based machine learning protocol for retinal is 1.5 times faster than that of NAMD to study reaction coordinates.