Acceleration of Reaction Space Projector Analysis Using Combinatorial Optimization: Application to Organic Chemical Reactions

In recent years, automated reaction path search methods have established the concept of a reaction route network. The Reaction Space Projector (ReSPer) visualizes the potential energy hypersurface into a lower-dimensional subspace using principal coordinates. The main time-consuming process in ReSPer is calculating the structural distance matrix, making it impractical for complex organic reaction route networks. We implemented the Alternate Optimization (AO) algorithm, one of the combinatorial optimizations, in ReSPer to reduce computational costs. Evaluations using gold clusters and the Au5 several reaction route networks showed that ReSPer-AO accurately computes distances with lower computational costs. Applying ReSPer-AO to the C5H8O reaction route network clarified dynamic conformation changes in its potential energy landscape. The ReSPer-AO method enables analysis of chemical reactions and dynamic conformations in a low-dimensional reaction space that accurately represents hydrocarbon reaction route networks.

Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation

The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.

Non-IRC Mechanism of Bimolecular Reactions with Submerged Barriers: A Case Study of Si+ + H2O Reaction

Chemical reactions with submerged barriers may feature interesting dynamic behaviors that are distinct from those with substantial barriers or those entirely dominated by capture. The Si+ + H2O reaction is a prototypical example, involving even two submerged saddle points along the reaction path: one for the direct dissociation of H (H-dissociation SP) and another for H migration from the O-side to the Si-side (H-migration SP). We investigated the intricacies of this process by employing quasi-classical trajectory calculations on an accurate, full-dimensional ab initio potential energy surface. Through careful trajectory analysis, an interesting nonintrinsic reaction coordinate mechanism was found to play an important role in producing SiOH+ and H. This new pathway is featured as that the H atoms do not form HSiOH+ complexes along the minimum-energy path but directly dissociate into the products after passing through the H-migration SP. Furthermore, based on artificially scaled potential energy surfaces (PESs), the impact of barrier height on the reaction is also explored. This work provides new insights into the dynamics of the Si+ + H2O reaction and enriches our understanding of reactions with submerged barriers.

Challenges for Kinetics Predictions via Neural Network Potentials: A Wilkinson's Catalyst Case

Ab initio kinetic studies are important to understand and design novel chemical reactions. While the Artificial Force Induced Reaction (AFIR) method provides a convenient and efficient framework for kinetic studies, accurate explorations of reaction path networks incur high computational costs. In this article, we are investigating the applicability of Neural Network Potentials (NNP) to accelerate such studies. For this purpose, we are reporting a novel theoretical study of ethylene hydrogenation with a transition metal complex inspired by Wilkinson's catalyst, using the AFIR method. The resulting reaction path network was analyzed by the Generative Topographic Mapping method. The network's geometries were then used to train a state-of-the-art NNP model, to replace expensive ab initio calculations with fast NNP predictions during the search. This procedure was applied to run the first NNP-powered reaction path network exploration using the AFIR method. We discovered that such explorations are particularly challenging for general purpose NNP models, and we identified the underlying limitations. In addition, we are proposing to overcome these challenges by complementing NNP models with fast semiempirical predictions. The proposed solution offers a generally applicable framework, laying the foundations to further accelerate ab initio kinetic studies with Machine Learning Force Fields, and ultimately explore larger systems that are currently inaccessible.

Multi-state Energy Landscape for Photoreaction of Stilbene and Dimethyl-stilbene

We have recently developed the reaction space projector (ReSPer) method, which constructs a reduced-dimensionality reaction space uniquely determined from reference reaction paths for a polyatomic molecular system and projects classical trajectories into the same reaction space. In this paper, we extend ReSPer to the analysis of photoreaction dynamics and relaxation processes of stilbene and present the concept of a "multi-state energy landscape," incorporating the ground- and excited-state reaction subspaces. The multi-state energy landscape successfully explains the previously established photoreaction processes of cis-stilbene, such as the cis-trans photoisomerization and photocyclization. In addition, we discuss the difference in the excited-state reaction dynamics between stilbene and 1,1'-dimethyl stilbene based on a common reaction subspace determined from the framework part of reference structures with different number of atoms. This approach allows us to target any molecule with a common framework, greatly expanding the applicability of the ReSPer analysis. The multi-state energy landscape provides fruitful insight into photochemical reactions, exploring the excited- and ground-state potential energy surfaces, as well as comprehensive reaction processes with nonradiative transitions between adiabatic states, within the stage of a reduced-dimensionality reaction space.

The principal component analysis of the ring deformation in the nonadiabatic surface hopping dynamics

The analysis of the leading active molecular motions in the on-the-fly trajectory surface hopping simulation provides the essential information to understand the geometric evolution in nonadiabatic dynamics. When the ring deformation is involved, the identification of the key active coordinates becomes challenging. A "hierarchical" protocol based on the dimensionality reduction and clustering approaches is proposed for the automatic analysis of the ring deformation in the nonadiabatic molecular dynamics. The representative system keto isocytosine is taken as the prototype to illustrate this protocol. The results indicate that the current hierarchical analysis protocol is a powerful way to clearly clarify both the major and minor active molecular motions of the ring distortion in nonadiabatic dynamics.

Reaction Space Projector (ReSPer) for Visualizing Dynamic Reaction Routes Based on Reduced-Dimension Space

To analyze chemical reaction dynamics based on a reaction path network, we have developed the "Reaction Space Projector" (ReSPer) method with the aid of the dimensionality reduction method. This program has two functions: the construction of a reduced-dimensionality reaction space from a molecular structure dataset, and the projection of dynamic trajectories into the low-dimensional reaction space. In this paper, we apply ReSPer to isomerization and bifurcation reactions of the Au₅ cluster and succeed in analyzing dynamic reaction routes involved in multiple elementary reaction processes, constructing complicated networks (called "closed islands") of nuclear permutation-inversion (NPI) isomerization reactions, and elucidating dynamic behaviors in bifurcation reactions with reference to bundles of trajectories. Interestingly, in the second application, we find a correspondence between the contribution ratios in the ability to visualize and the symmetry of the morphology of closed islands. In addition, the third application suggests the existence of boundaries that determine the selectivity in bifurcation reactions, which was discussed in the phase space. The ReSPer program is a versatile and robust tool to clarify dynamic reaction mechanisms based on the reduced-dimensionality reaction space without prior knowledge of target reactions.

Visualization of reaction route map and dynamical trajectory in reduced dimension

In the quantum chemical approach, chemical reaction mechanisms are investigated based on a potential energy surface (PES). Automated reaction path search methods enable us to construct a global reaction route map containing multiple reaction paths corresponding to a series of elementary reaction processes. The on-the-fly molecular dynamics (MD) method provides a classical trajectory exploring the full-dimensional PES based on electronic structure calculations. We have developed two reaction analysis methods, the on-the-fly trajectory mapping method and the reaction space projector (ReSPer) method, by introducing a structural similarity to a pair of geometric structures and revealed dynamic aspects affecting chemical reaction mechanisms. In this review, we will present the details of these analysis methods and discuss the dynamics effects of reaction path curvature and reaction path bifurcation with applications to the CH₃OH + OH^- collision reaction and the Au₅ cluster branching and isomerization reactions.

Mode-specific dynamics in multichannel reaction NH+ + H2

The vibrational- and rotational-mode specificity of the multichannel NH⁺ + H₂ reaction is studied on a recently constructed ab initio-based global potential energy surface using an initial state selected quasi-classical trajectory method, and the trajectories are analyzed using an isometric feature mapping and k-means approach. All excitation modes promote two reactions (R1: NH'⁺ + H₂ → NH⁺ + HH' and R4: NH'⁺ + H₂ → NH₂⁺ + H') where both NH and HH bonds are broken, but reduce the reactivity of the proton-transfer reaction R2 (NH'⁺ + H₂ → N + H'H₂⁺) at low collision energies. For the hydrogen-transfer reaction R3 (NH'⁺ + H₂ → HNH'⁺ + H), the rotational excitation of NH⁺ enhances the reactivity remarkably, while its vibrational excitation has an inhibiting effect on the reaction. The trajectory analyses show that the vibrational and rotational excitations of NH⁺ make R3 tend to go over a submerged saddle point instead of extracting hydrogen atoms directly. On the other hand, the motions of the H₂ reactant facilitate the enhancement of the reactivity but they do not affect the mechanism of R3. In addition, the results suggest that the coupling of the isometric feature mapping and the k-means approach in the trajectory analysis is an appropriate tool for reaction-dynamics studies.

General Many-Body Framework for Data-Driven Potentials with Arbitrary Quantum Mechanical Accuracy: Water as a Case Study

We present a general framework for the development of data-driven many-body (MB) potential energy functions (MB-QM PEFs) that represent the interactions between small molecules at an arbitrary quantum-mechanical (QM) level of theory. As a demonstration, a family of MB-QM PEFs for water is rigorously derived from density functionals belonging to different rungs across Jacob's ladder of approximations within density functional theory (MB-DFT) and from Møller-Plesset perturbation theory (MB-MP2). Through a systematic analysis of individual MB contributions to the interaction energies of water clusters, we demonstrate that all MB-QM PEFs preserve the same accuracy as the corresponding ab initio calculations, with the exception of those derived from density functionals within the generalized gradient approximation (GGA). The differences between the DFT and MB-DFT results are traced back to density-driven errors that prevent GGA functionals from accurately representing the underlying molecular interactions for different cluster sizes and hydrogen-bonding arrangements. We show that this shortcoming may be overcome, within the MB formalism, by using density-corrected functionals (DC-DFT) that provide a more consistent representation of each individual MB contribution. This is demonstrated through the development of a MB-DFT PEF derived from DC-PBE-D3 data, which more accurately reproduce the corresponding ab initio results.

Analysis of bath motion in MM-SQC dynamics via dimensionality reduction approach: Principal component analysis

The system-plus-bath model is an important tool to understand the nonadiabatic dynamics of large molecular systems. Understanding the collective motion of a large number of bath modes is essential for revealing their key roles in the overall dynamics. Here, we applied principal component analysis (PCA) to investigate the bath motion in the basis of a large dataset generated from the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian nonadiabatic dynamics for the excited-state energy transfer in the Frenkel-exciton model. The PCA method clearly elucidated that two types of bath modes, which either display strong vibronic coupling or have frequencies close to that of the electronic transition, are important to the nonadiabatic dynamics. These observations were fully consistent with the physical insights. The conclusions were based on the PCA of the trajectory data and did not involve significant pre-defined physical knowledge. The results show that the PCA approach, which is one of the simplest unsupervised machine learning dimensionality reduction methods, is a powerful one for analyzing complicated nonadiabatic dynamics in the condensed phase with many degrees of freedom.

Neural Network Potential Energy Surfaces for Small Molecules and Reactions

We review progress in neural network (NN)-based methods for the construction of interatomic potentials from discrete samples (such as ab initio energies) for applications in classical and quantum dynamics including reaction dynamics and computational spectroscopy. The main focus is on methods for building molecular potential energy surfaces (PES) in internal coordinates that explicitly include all many-body contributions, even though some of the methods we review limit the degree of coupling, due either to a desire to limit computational cost or to limited data. Explicit and direct treatment of all many-body contributions is only practical for sufficiently small molecules, which are therefore our primary focus. This includes small molecules on surfaces. We consider direct, single NN PES fitting as well as more complex methods that impose structure (such as a multibody representation) on the PES function, either through the architecture of one NN or by using multiple NNs. We show how NNs are effective in building representations with low-dimensional functions including dimensionality reduction. We consider NN-based approaches to build PESs in the sums-of-product form important for quantum dynamics, ways to treat symmetry, and issues related to sampling data distributions and the relation between PES errors and errors in observables. We highlight combinations of NNs with other ideas such as permutationally invariant polynomials or sums of environment-dependent atomic contributions, which have recently emerged as powerful tools for building highly accurate PESs for relatively large molecular and reactive systems.