Analysis of Trajectory Similarity and Configuration Similarity in On-the-Fly Surface-Hopping Simulation on Multi-Channel Nonadiabatic Photoisomerization Dynamics
arXiv:1806.07655 · doi:10.1063/1.5048049
Abstract
We propose an 'automatic' approach to analyze the results of the on-the-fly trajectory surface hopping simulation on the multi-channel nonadiabatic photoisomerization dynamics by considering the trajectory similarity and the configuration similarity. We choose a representative system phytochromobilin chromophore model to illustrate the analysis protocol. After a large number of trajectories are obtained, it is possible to define the similarity of different trajectories by the Fréchet distance and to employ the trajectory clustering analysis to divide all trajectories into several clusters. Each cluster in principle represents a photoinduced isomerization reaction channel. This idea provides an effective approach to understand the branching ratio of the multi-channel photoisomerization dynamics. For each cluster the dimensionality reduction is employed to understand the configuration similarity in the trajectory propagation, which provides the understanding of the major geometry evolution features in each reaction channel. The results show that this analysis protocol not only assigns all trajectories into different photoisomerization reaction channels, but also extracts the major molecular motion without the requirement of the pre-known knowledge of the active photoisomerization site. As a side product of this analysis tool, we can also easily to find the so-called 'typical' or 'representative' trajectory for each reaction channel.
References in corpus (9)
- Scikit-learn: Machine Learning in Python
- On representing chemical environments
- Deep Potential Molecular Dynamics: a scalable model with the accuracy of quantum mechanics
- API design for machine learning software: experiences from the scikit-learn project
- On-the-fly CASPT2 surface hopping dynamics
- Landau-Zener type surface hopping algorithms
- Path Similarity Analysis: a Method for Quantifying Macromolecular Pathways
- Nonadiabatic nuclear dynamics of the ammonia cation studied by surface hopping classical trajectory calculations
- Constructing grids for molecular quantum dynamics using an autoencoder