Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe role of depth and flatness of a potential energy surface in chemical reaction dynamics
173
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this study, we analyze how changes in the geometry of a potential energy surface in terms of depth and flatness can affect the reaction dynamics. We formulate depth and flatness in the context of one and two degree-of-freedom (DOF) Hamiltonian normal form for the saddle-node bifurcation and quantify their influence on chemical reaction dynamics. In a recent work, Garcia-Garrido, Naik, and Wiggins illustrated how changing the well-depth of a potential energy surface (PES) can lead to a saddle-node bifurcation. They have shown how the geometry of cylindrical manifolds associated with the rank-1 saddle changes en route to the saddle-node bifurcation. Using the formulation presented here, we show how changes in the parameters of the potential energy control the depth and flatness and show their role in the quantitative measures of a chemical reaction. We quantify this role of the depth and flatness by calculating the ratio of the bottleneck-width and well-width, reaction probability (also known as transition fraction or population fraction), gap time (or first passage time) distribution, and directional flux through the dividing surface (DS) for small to high values of total energy. The results obtained for these quantitative measures are in agreement with the qualitative understanding of the reaction dynamics.
rate research
Read More
We examine the phase space structures that govern reaction dynamics in the absence of critical points on the potential energy surface. We show that in the vicinity of hyperbolic invariant tori it is possible to define phase space dividing surfaces that are analogous to the dividing surfaces governing transition from reactants to products near a critical point of the potential energy surface. We investigate the problem of capture of an atom by a diatomic molecule and show that a normally hyperbolic invariant manifold exists at large atom-diatom distances, away from any critical points on the potential. This normally hyperbolic invariant manifold is the anchor for the construction of a dividing surface in phase space, which defines the outer or loose transition state governing capture dynamics. We present an algorithm for sampling an approximate capture dividing surface, and apply our methods to the recombination of the ozone molecule. We treat both 2 and 3 degree of freedom models with zero total angular momentum. We have located the normally hyperbolic invariant manifold from which the orbiting (outer) transition state is constructed. This forms the basis for our analysis of trajectories for ozone in general, but with particular emphasis on the roaming trajectories.
A model Hamiltonian for the reaction CH$_4^+ rightarrow$ CH$_3^+$ + H, parametrized to exhibit either early or late inner transition states, is employed to investigate the dynamical characteristics of the roaming mechanism. Tight/loose transition states and conventional/roaming reaction pathways are identified in terms of time-invariant objects in phase space. These are dividing surfaces associated with normally hyperbolic invariant manifolds (NHIMs). For systems with two degrees of freedom NHIMS are unstable periodic orbits which, in conjunction with their stable and unstable manifolds, unambiguously define the (locally) non-recrossing dividing surfaces assumed in statistical theories of reaction rates. By constructing periodic orbit continuation/bifurcation diagrams for two values of the potential function parameter corresponding to late and early transition states, respectively, and using the total energy as another parameter, we dynamically assign different regions of phase space to reactants and products as well as to conventional and roaming reaction pathways. The classical dynamics of the system are investigated by uniformly sampling trajectory initial conditions on the dividing surfaces. Trajectories are classified into four different categories: direct reactive and non reactive trajectories,which lead to the formation of molecular and radical products respectively, and roaming reactive and non reactive orbiting trajectories, which represent alternative pathways to form molecular and radical products. By analysing gap time distributions at several energies we demonstrate that the phase space structure of the roaming region, which is strongly influenced by non-linear resonances between the two degrees of freedom, results in nonexponential (nonstatistical) decay.
We study reaction dynamics on a model potential energy surface exhibiting post-transition state bifurcation in the vicinity of a valley ridge inflection point. We compute fractional yields of products reached after the VRI region is traversed, both with and without dissipation. It is found that apparently minor variations in the potential lead to significant changes in the reaction dynamics. Moreover, when dissipative effects are incorporated, the product ratio depends in a complicated and highly non-monotonic fashion on the dissipation parameter. Dynamics in the vicinity of the VRI point itself play essentially no role in determining the product ratio, except in the highly dissipative regime.
A variety of natural phenomena comprises a huge number of competing reactions and short-lived intermediates. Any study of such processes requires the discovery and accurate modeling of their underlying reaction network. However, this task is challenging due to the complexity in exploring all the possible pathways and the high computational cost in accurately modeling a large number of reactions. Fortunately, very often these processes are dominated by only a limited subset of the networks reaction pathways. In this work we propose a novel computationally inexpensive method to identify and select the key pathways of complex reaction networks, so that high-level ab-initio calculations can be more efficiently targeted at these critical reactions. The method estimates the relative importance of the reaction pathways for given reactants by analyzing the accelerated evolution of hundreds of replicas of the system and detecting products formation. This acceleration-detection method is able to tremendously speed up the reactivity of uni- and bimolecular reactions, without requiring any previous knowledge of products or transition states. Importantly, the method is efficiently iterative, as it can be straightforwardly applied for the most frequently observed products, therefore providing an efficient algorithm to identify the key reactions of extended chemical networks. We verified the validity of our approach on three different systems, including the reactivity of t-decalin with a methyl radical, and in all cases the expected behavior was recovered within statistical error.
A new potential energy surface (PES) and dynamical study are presented of the reactive process between H2CO + OH towards the formation of HCO + H2O and HCOOH + H. In this work a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, what prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many body term has been used fitted with a novel artificial neural networks scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing a good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated to excitations of the normal modes of the ring polymer.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
