Abstract:This work addresses the construction and use of low-dimensional invariant manifolds to simplify complex chemical kinetics. Typically, chemical kinetic systems have a wide range of time scales. As a consequence, reaction trajectories rapidly approach a hierarchy of attracting manifolds of decreasing dimension in the full composition space. In previous research, several different methods have been proposed to identify these low-dimensional attracting manifolds. Here we propose a new method based on an invariant constrained equilibrium edge (ICE) manifold. This manifold (of dimension nr) is generated by the reaction trajectories emanating from its (nr-1)-dimensional edge, on which the composition is in a constrained equilibrium state. A reasonable choice of the nr represented variables (e.g., nr "major" species) ensures that there exists a unique point on the ICE manifold corresponding to each realizable value of the represented variables. The process of identifying this point is referred to as species reconstruction. A second contribution of this work is a local method of species reconstruction, called ICE-PIC, which is based on the ICE manifold and uses preimage curves (PICs). The ICE-PIC method is local in the sense that species reconstruction can be performed without generating the whole of the manifold (or a significant portion thereof). The ICE-PIC method is the first approach that locally determines points on a low-dimensional invariant manifold, and its application to high-dimensional chemical systems is straightforward. The "inputs" to the method are the detailed kinetic mechanism and the chosen reduced representation (e.g., some major species). The ICE-PIC method is illustrated and demonstrated using an idealized H2O system with six chemical species. It is then tested and compared to three other dimension-reduction methods for the test case of a one-dimensional premixed laminar flame of stoichiometric hydrogen/air, which is described by a detailed mechanism containing nine species and 21 reactions. It is shown that the error incurred by the ICE-PIC method with four represented species is small across the whole flame, even in the low temperature region.

Reaction Coordinates by Nonlinear Dimensionality Reduction

Reaction Coordinates for Conformational Transitions using Linear Discriminant Analysis on Positions

Data-driven Computation of Molecular Reaction Coordinates

A global view of reactive coordinate manifolds from nonlinear dimensionality reduction

Low-dimensional projection of reactivity classes in chemical reaction dynamics using supervised dimensionality reduction

Reaction Networks Resemble Low-Dimensional Regular Lattices

From Data to Chemistry: Revealing Causality and Reaction Coordinates through Interpretable Machine Learning in Supramolecular Transition Metal Catalysis

The Invariant Constrained Equilibrium Edge Preimage Curve Method for the Dimension Reduction of Chemical Kinetics

Reaction Coordinates are Optimal Channels of Energy Flow

Reaction coordinate flows for model reduction of molecular kinetics

Low-Dimensional Projection of Reactive Islands in Chemical Reaction Dynamics Using a Supervised Dimensionality Reduction Method

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

Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide

A Mode Evolution Metric to Extract Reaction Coordinates for Biomolecular Conformational Transitions

Investigating the hyperparameter space of deep neural network models for reaction coordinates: Revisiting the solvent coordinate in alanine dipeptide isomerization

Locating Transition Zone in Phase Space

Automated handling of complex chemical structures in Z‐matrix coordinates—The chemcoord library

The Geometry of Reaction Trajectories and Attracting Manifolds in Composition Space

Accuracy of reaction coordinate based rate theories for modelling chemical reactions: insights from the thermal isomerization in retinal

Multifaceted Stereoselectivity in Polyatomic Reactions

Solvent reaction coordinate for an S$_N$2 reaction