Exact quantum calculations of the kinetic isotope effect: cross sections and rate constants for the F+HD reaction and role of tunneling

Quantum Dynamics Studies of the Significant Intramolecular Isotope Effects on the Nonadiabatic Be+(2P) + HD → BeH+/BeD+ + D/H Reaction

Quantum time-dependent wave packet dynamics studies on the nonadiabatic Be⁺(²P) + HD → BeH⁺/BeD⁺ + D/H reaction are performed for the first time employing recently constructed diabatic potential energy surfaces. Strong intramolecular isotope effects and unusual results are presented, which are attributed to the dynamic effects of shallow wells induced by avoided crossing on the diagonal V₂₂^d surface. The BeH⁺ + D and BeD⁺ + H channels are dominated by high-J and low-J partial waves, respectively. The BeD⁺/BeH⁺ branching ratio is larger than 10 at low energy and gradually decreases with increasing collision energy. The BeH⁺ product is primarily distributed at low vibrational states, whereas there exists an obvious population inversion of vibrational states on the BeD⁺ product. The results of differential cross sections suggest that the formation of the BeH⁺ + D channel favors a direct reaction process, while the BeD⁺ + H channel is mainly generated by the complex-forming mechanism.

Dynamics study on the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction: insertion-abstraction mechanism

Time-dependent wave packet calculations are carried out for two reaction channels of the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction. The potential well on the excited state potential energy surface makes the reaction preferable to proceed through the insertion reaction path. The dominance of the NaD + H reaction channel and product rotational state distributions are found to be in agreement with the characteristics of typical adiabatic insertion reactions. However, significant forward scattering peaks in the differential cross sections (DCS) are found to be inconsistent with the forward-backward symmetric scattering characteristic of typical adiabatic insertion reactions, which indicate that the Na(3p) + HD reaction is dominated by a direct reaction mechanism. The comparison between adiabatic and non-adiabatic calculated DCSs reveals that the non-adiabatic couplings in the reaction could reduce the lifetime of the intermediate complex. Finally, the insertion-abstraction mechanism is put forward for the non-adiabatic Na(3p) + HD reaction.

Benchmark Quantum Kinetics at Low Temperatures toward Absolute Zero and Role of Entrance Channel Wells on Tunneling, Virtual States, and Resonances: The F + HD Reaction

This paper reports a study of the quantum reaction dynamics and kinetics of the F + HD reaction at low and ultralow temperatures, focusing on the range from the Wigner limit up to 50 K. Close coupling time-independent quantum reactive scattering calculations for the production of HF and DF molecules have been carried out on two potential energy surfaces differing in the description of the reaction entrance channel. This case is computationally more demanding than the cases of F with H₂ and D₂ ( De Fazio et al. Frontiers in Chemistry 2019 , 7 , 328 ) but offers a wider phenomenology regarding the roles of quantum mechanical effects of tunneling, of virtual states, and of resonances. The results show that at the temperatures in the cold and ultracold regimes small changes in the entrance channel long-range interaction induce surprising near threshold features. The presence of a virtual state close to the reactive threshold gives rise to a marked anti-Arrhenius behavior of the rate constants below 100 mK. This effect enhances reaction rates by about 2 orders of magnitude, making them of the same order as those at room temperature and confining the onset of the Wigner regime in the microkelvin region.

Temperature Dependence of Rate Processes Beyond Arrhenius and Eyring: Activation and Transitivity

Advances in the understanding of the dependence of reaction rates from temperature, as motivated from progress in experiments and theoretical tools (e. g., molecular dynamics), are needed for the modeling of extreme environmental conditions (e.g., in astrochemistry and in the chemistry of plasmas). While investigating statistical mechanics perspectives (Aquilanti et al., 2017b, 2018), the concept of transitivity was introduced as a measure for the propensity for a reaction to occur. The Transitivity plot is here defined as the reciprocal of the apparent activation energy vs. reciprocal absolute temperature. Since the transitivity function regulates transit in physicochemical transformations, not necessarily involving reference to transition-state hypothesis of Eyring, an extended version is here proposed to cope with general types of transformations. The transitivity plot permits a representation where deviations from Arrhenius behavior are given a geometrical meaning and make explicit a positive or negative linear dependence of transitivity for sub- and super-Arrhenius cases, respectively. To first-order in reciprocal temperature, the transitivity function models deviations from linearity in Arrhenius plots as originally proposed by Aquilanti and Mundim: when deviations are increasingly larger, other phenomenological formulas, such as Vogel-Fulcher-Tammann, Nakamura-Takayanagi-Sato, and Aquilanti-Sanches-Coutinho-Carvalho are here rediscussed from the transitivity concept perspective and with in a general context. Emphasized is the interest of introducing into this context modifications to a very successful tool of theoretical kinetics, Eyring's Transition-State Theory: considering the behavior of the transitivity function at low temperatures, in order to describe deviation from Arrhenius behavior under the quantum tunneling regime, a "d-TST" formulation was previously introduced (Carvalho-Silva et al., 2017). In this paper, a special attention is dedicated to a derivation of the temperature dependence of viscosity, making explicit reference to feature of the transitivity function, which in this case generally exhibits a super-Arrhenius behavior. This is of relevance also for advantages of using the transitivity function for diffusion-controlled phenomena.

Kinetics of low-temperature transitions and a reaction rate theory from non-equilibrium distributions

This article surveys the empirical information which originated both by laboratory experiments and by computational simulations, and expands previous understanding of the rates of chemical processes in the low-temperature range, where deviations from linearity of Arrhenius plots were revealed. The phenomenological two-parameter Arrhenius equation requires improvement for applications where interpolation or extrapolations are demanded in various areas of modern science. Based on Tolman's theorem, the dependence of the reciprocal of the apparent activation energy as a function of reciprocal absolute temperature permits the introduction of a deviation parameter d covering uniformly a variety of rate processes, from those where quantum mechanical tunnelling is significant and d < 0, to those where d > 0, corresponding to the Pareto-Tsallis statistical weights: these generalize the Boltzmann-Gibbs weight, which is recovered for d = 0. It is shown here how the weights arise, relaxing the thermodynamic equilibrium limit, either for a binomial distribution if d > 0 or for a negative binomial distribution if d < 0, formally corresponding to Fermion-like or Boson-like statistics, respectively. The current status of the phenomenology is illustrated emphasizing case studies; specifically (i) the super-Arrhenius kinetics, where transport phenomena accelerate processes as the temperature increases; (ii) the sub-Arrhenius kinetics, where quantum mechanical tunnelling propitiates low-temperature reactivity; (iii) the anti-Arrhenius kinetics, where processes with no energetic obstacles are rate-limited by molecular reorientation requirements. Particular attention is given for case (i) to the treatment of diffusion and viscosity, for case (ii) to formulation of a transition rate theory for chemical kinetics including quantum mechanical tunnelling, and for case (iii) to the stereodirectional specificity of the dynamics of reactions strongly hindered by the increase of temperature.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

Deformed transition-state theory: Deviation from Arrhenius behavior and application to bimolecular hydrogen transfer reaction rates in the tunneling regime

A formulation is presented for the application of tools from quantum chemistry and transition-state theory to phenomenologically cover cases where reaction rates deviate from Arrhenius law at low temperatures. A parameter d is introduced to describe the deviation for the systems from reaching the thermodynamic limit and is identified as the linearizing coefficient in the dependence of the inverse activation energy with inverse temperature. Its physical meaning is given and when deviation can be ascribed to quantum mechanical tunneling its value is calculated explicitly. Here, a new derivation is given of the previously established relationship of the parameter d with features of the barrier in the potential energy surface. The proposed variant of transition state theory permits comparison with experiments and tests against alternative formulations. Prescriptions are provided and implemented to three hydrogen transfer reactions: CH₄  + OH → CH₃  + H₂ O, CH₃ Cl + OH → CH₂ Cl + H₂ O and H₂  + CN → H + HCN, widely investigated both experimentally and theoretically. © 2016 Wiley Periodicals, Inc.

Stereodirectional Origin of anti-Arrhenius Kinetics for a Tetraatomic Hydrogen Exchange Reaction: Born-Oppenheimer Molecular Dynamics for OH + HBr

Among four-atom processes, the reaction OH + HBr → H2O + Br is one of the most studied experimentally: its kinetics has manifested an unusual anti-Arrhenius behavior, namely, a marked decrease of the rate constant as the temperature increases, which has intrigued theoreticians for a long time. Recently, salient features of the potential energy surface have been characterized and most kinetic aspects can be considered as satisfactorily reproduced by classical trajectory simulations. Motivation of the work reported in this paper is the investigation of the stereodirectional dynamics of this reaction as the prominent reason for the peculiar kinetics: we started in a previous Letter ( J. Phys. Chem. Lett. 2015 , 6 , 1553 - 1558 ) a first-principles Born-Oppenheimer "canonical" molecular dynamics approach. Trajectories are step-by-step generated on a potential energy surface quantum mechanically calculated on-the-fly and are thermostatically equilibrated to correspond to a specific temperature. Here, refinements of the method permitted a major increase of the number of trajectories and the consideration of four temperatures -50, +200, +350, and +500 K, for which the sampling of initial conditions allowed us to characterize the stereodynamical effect. The role is documented of the adjustment of the reactants' mutual orientation to encounter the entrance into the "cone of acceptance" for reactivity. The aperture angle of this cone is dictated by a range of directions of approach compatible with the formation of the specific HOH angle of the product water molecule; and consistently the adjustment is progressively less effective the higher the kinetic energy. Qualitatively, this emerging picture corroborates experiments on this reaction, involving collisions of aligned and oriented molecular beams, and covering a range of energies higher than the thermal ones. The extraction of thermal rate constants from this molecular dynamics approach is discussed and the systematic sampling of the canonical ensemble is indicated as needed for quantitative comparison with the kinetic experiments.

Complementarity between Quantum and Classical Mechanics in Chemical Modeling. The H + HeH+ → H2 + + He Reaction: A Rigourous Test for Reaction Dynamics Methods

In this work we present a dynamical study of the H + HeH+ → H2+ + He reaction in a collision energy range from 0.1 meV to 10 eV, suitable to be used in applicative models. The paper extends and complements a recent work [ Phys. Chem. Chem. Phys. 2014, 16, 11662] devoted to the characterization of the reactivity from the ultracold regime up to the three-body dissociation breakup. In particular, the accuracy of the quasi-classical trajectory method below the three-body dissociation threshold has been assessed by a detailed comparison with previous calculations performed with different reaction dynamics methods, whereas the reliability of the results in the high energy range has been checked by a direct comparison with the available experimental data. Integral cross sections for several HeH+ roto-vibrational states have been analyzed and used to understand the extent of quantum effects in the reaction dynamics. By using the quasi-classical trajectory method and quantum mechanical close coupling data, respectively, in the high and low collision energy ranges, we obtain highly accurate thermal rate costants until 15 000 K including all (178) the roto-vibrational bound and quasi-bound states of HeH+. The role of the collision-induced dissociation is also discussed and explicitly calculated for the ground roto-vibrational state of HeH+.

Complex angular momentum theory of state-to-state integral cross sections: resonance effects in the F + HD → HF(v′ = 3) + D reaction

State-to-state reactive integral cross sections (ICSs) are often affected by quantum mechanical resonances, especially in the neighborhood of a reactive threshold.

Theoretical reaction kinetics astride the transition between moderate and deep tunneling regimes: the F + HD case

For the reaction between F and HD, giving HF + D and DF + H, the rate constants, obtained from rigorous quantum scattering calculations at temperatures ranging from 350 K down to 100 K, show deviations from the Arrhenius behavior that have been interpreted in terms of tunneling of either H or D atoms through a potential energy barrier. The interval of temperature investigated extends from above to below a crossover value Tc, a transition temperature separating the moderate and deep quantum tunneling regimes. Below Tc, the rate of the H or D exchange reaction is controlled by the prevalence of tunneling over the thermal activation mechanism. In this temperature range, Bell's early treatment of quantum tunneling, based on a semiclassical approximation for the barrier permeability, provides a reliable tool to quantitatively account for the contribution of the tunneling effect. This treatment is here applied for extracting from rate constants properties of the effective tunneling path, such as the activation barrier height and width. We show that this is a way of parametrizing the dependence of the apparent activation energy on temperature useful for both calculated and experimental rate constants in an ample interval of temperature, from above to below Tc, relevant for modelization of astrophysical and in general very low-temperature environments.

The He + H2+ → HeH+ + H reaction: ab initio studies of the potential energy surface, benchmark time-independent quantum dynamics in an extended energy range and comparison with experiments

In this work we critically revise several aspects of previous ab initio quantum chemistry studies [P. Palmieri et al., Mol. Phys. 98, 1835 (2000); C. N. Ramachandran et al., Chem. Phys. Lett. 469, 26 (2009)] of the HeH(2)(+) system. New diatomic curves for the H(2)(+) and HeH(+) molecular ions, which provide vibrational frequencies at a near spectroscopic level of accuracy, have been generated to test the quality of the diatomic terms employed in the previous analytical fittings. The reliability of the global potential energy surfaces has also been tested performing benchmark quantum scattering calculations within the time-independent approach in an extended interval of energies. In particular, the total integral cross sections have been calculated in the total collision energy range 0.955-2.400 eV for the scattering of the He atom by the ortho- and para-hydrogen molecular ion. The energy profiles of the total integral cross sections for selected vibro-rotational states of H(2)(+) (v = 0,...,5 and j = 1,...,7) show a strong rotational enhancement for the lower vibrational states which becomes weaker as the vibrational quantum number increases. Comparison with several available experimental data is presented and discussed.

THEORETICAL STUDY ON PARTIAL POTENTIAL SURFACE AND SCATTERING RESONANCE STATE OF THE ASYMMETRICAL H EXCHANGING X + H2O → XH + OH (X = Cl, F, H) REACTIVE SYSTEM

In the course of an extensive investigation aimed at understanding the detailed mechanism of a prototypical polyatomic reaction, several remarkable observations were uncovered. To interpret these findings, we surmise the existence of a reactive resonance in this polyatomic reaction. The system of concern is X + H 2 O → XH + OH ( X = Cl, F, H ), and it belongs to a type of asymmetrical H exchanging reaction. Because the concerned electronic number is quite small, it is usually regarded as a template constructed the model of theory or experiment. In this paper, compared with the symmetrical systems, we constructed the partial potential surfaces of asymmetrical systems in order to research their resonance states. All results are in good agreement with previous theoretical and experimental works.

Product pair correlation in bimolecular reactions

The vast majority of chemical reactions involve polyatomic species as reactants and/or products. The added degree of complexity offers opportunities to address dynamical questions other than those already encountered in a typical atom + diatom reaction. Product pair correlation is one of them. This article introduces the basic concept, outlines the experimental approach we developed and then highlights some of the applications to bimolecular reaction dynamics. Particular emphasis is placed on the information contents and unique insights gained from this type of measurements, which otherwise would have been lost by the conventional approach.