Neural network approach to time-dependent dividing surfaces in classical reaction dynamics

Phase space geometry of isolated to condensed chemical reactions

The complexity of gas and condensed phase chemical reactions has generally been uncovered either approximately through transition state theories or exactly through (analytic or computational) integration of trajectories. These approaches can be improved by recognizing that the dynamics and associated geometric structures exist in phase space, ensuring that the propagator is symplectic as in velocity-Verlet integrators and by extending the space of dividing surfaces to optimize the rate variationally, respectively. The dividing surface can be analytically or variationally optimized in phase space, not just over configuration space, to obtain more accurate rates. Thus, a phase space perspective is of primary importance in creating a deeper understanding of the geometric structure of chemical reactions. A key contribution from dynamical systems theory is the generalization of the transition state (TS) in terms of the normally hyperbolic invariant manifold (NHIM) whose geometric phase-space structure persists under perturbation. The NHIM can be regarded as an anchor of a dividing surface in phase space and it gives rise to an exact non-recrossing TS theory rate in reactions that are dominated by a single bottleneck. Here, we review recent advances of phase space geometrical structures of particular relevance to chemical reactions in the condensed phase. We also provide conjectures on the promise of these techniques toward the design and control of chemical reactions.

The Middle Science: Traversing Scale In Complex Many-Body Systems

A roadmap is developed that integrates simulation methodology and data science methods to target new theories that traverse the multiple length- and time-scale features of many-body phenomena.

Thermal decay rates of an activated complex in a driven model chemical reaction

Recent work has shown that in a nonthermal, multidimensional system, the trajectories in the activated complex possess different instantaneous and time-averaged reactant decay rates. Under dissipative dynamics, it is known that these trajectories, which are bound on the normally hyperbolic invariant manifold (NHIM), converge to a single trajectory over time. By subjecting these dissipative systems to thermal noise, we find fluctuations in the saddle-bound trajectories and their instantaneous decay rates. Averaging over these instantaneous rates results in the decay rate of the activated complex in a thermal system. We find that the temperature dependence of the activated complex decay in a thermal system can be linked to the distribution of the phase space resolved decay rates on the NHIM in the nondissipative case. By adjusting the external driving of the reaction, we show that it is possible to influence how the decay rate of the activated complex changes with rising temperature.

Neural network approach for the dynamics on the normally hyperbolic invariant manifold of periodically driven systems

Chemical reactions in multidimensional systems are often described by a rank-1 saddle, whose stable and unstable manifolds intersect in the normally hyperbolic invariant manifold (NHIM). Trajectories started on the NHIM in principle never leave this manifold when propagated forward or backward in time. However, the numerical investigation of the dynamics on the NHIM is difficult because of the instability of the motion. We apply a neural network to describe time-dependent NHIMs and use this network to stabilize the motion on the NHIM for a periodically driven model system with two degrees of freedom. The method allows us to analyze the dynamics on the NHIM via Poincaré surfaces of section (PSOS) and to determine the transition-state (TS) trajectory as a periodic orbit with the same periodicity as the driving saddle, viz. a fixed point of the PSOS surrounded by near-integrable tori. Based on transition state theory and a Floquet analysis of a periodic TS trajectory we compute the rate constant of the reaction with significantly reduced numerical effort compared to the propagation of a large trajectory ensemble.