Bimolecular reaction rates from ring polymer molecular dynamics

Kinetic study of the CN + C2H6 hydrogen abstraction reaction based on an analytical potential energy surface

The temperature dependence of the thermal rate constants and kinetic isotope effects (KIE) of the CN + C2H6 gas-phase hydrogen abstraction reaction was theoretically determined within the 25-1000 K temperature range, i.e., from very low- to high-temperature regimes. Based on a recently developed full-dimensional analytical potential energy surface fitted to highly accurate explicitly correlated ab initio calculations, three different kinetic theories were used: canonical variational transition state theory (CVT), quasiclassical trajectory theory (QCT), and ring polymer molecular dynamics (RPMD) method for the computation of rate constants. We found that the thermal rate constants obtained with the three theories show a V-shaped temperature dependence, with a pronounced minimum near 200 K, qualitatively reproducing the experimental measurements. Among the three methods used in this work, the QCT and RPMD methods have the best agreement with the experiment at low and high temperatures, respectively, while the CVT model shows the largest discrepancies. The significant increase in the rate constant at very low temperatures in this very exothermic and practically barrierless reaction could be attributed to the large value of the impact parameter, possibly ruling out the role of the tunneling effect and the intermediate complexes in the entrance channel. The theoretical H/D KIE depicted a "normal" behaviour, i.e., values greater than unity, emulating the experimental measurements and improving previous theoretical results. Finally, the discrepancies between theory and experiments were analysed as a function of several factors, such as limitations of the kinetics theories and the potential energy surface, as well as the uncertainties in the experimental measurements.

Dynamics studies for the multi-well and multi-channel reaction of OH with C2H2 on a full-dimensional global potential energy surface

The C2H2 + OH reaction is an important acetylene oxidation pathway in the combustion process, as well as a typical multi-well and multi-channel reaction. Here, we report an accurate full-dimensional machine learning-based potential energy surface (PES) for the C2H2 + OH reaction at the UCCSD(T)-F12b/cc-pVTZ-F12 level, based on about 475 000 ab initio points. Extensive quasi-classical trajectory (QCT) calculations were performed on the newly developed PES to obtain detailed dynamic data and analyze reaction mechanisms. Below 1000 K, the C2H2 + OH reaction produces H + OCCH2 and CO + CH3. With increasing temperature, the product channels H2O + C2H and H + HCCOH are accessible and the former dominates above 1900 K. It is found that the formation of H2O + C2H is dominated by a direct reaction process, while other channels belong to the indirect mechanism involving long-lived intermediates along the reaction pathways. At low temperatures, the C2H2 + OH reaction behaves like an unimolecular reaction due to the unique PES topographic features, of which the dynamic features are similar to the decomposition of energy-rich complexes formed by C2H2 + OH collision. The classification of trajectories that undergo different reaction pathways to generate each product and their product energy distributions were also reported in this work. This dynamic information may provide a deep understanding of the C2H2 + OH reaction.

Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package

A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.

Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates

A ring polymer molecular dynamics (RPMD) method is proposed for the calculation of the dissociative chemisorption rate coefficient on surfaces. The RPMD rate theory is capable of handling quantum effects such as the zero-point energy and tunneling in dissociative chemisorption, while it relies on classical trajectories for the simulation. Applications to H2 dissociative chemisorption are demonstrated. For the highly activated process on Ag(111), strong deviations from Arrhenius behavior are found at low temperatures and attributed to tunneling. On Pt(111), where the dissociation has a barrierless pathway, the RPMD rate coefficient is found to agree with the experimentally derived thermal sticking coefficient within a factor of 2 over a large temperature range. Significant quantum effects are also found.

First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion

Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.

Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations

In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.

New Stable and Fast Ring-Polymer Molecular Dynamics for Calculating Bimolecular Rate Coefficients with an Example of OH + CH4

The accurate and efficient calculation of the rate coefficients of chemical reactions is a key issue in the research of chemical dynamics. In this work, by applying the dimension-free ultrastable Cayley propagator, the thermal rate coefficients of a prototypic high dimensional chemical reaction OH + CH₄ → H₂O + CH₃ in the temperature range of 200 to 1500 K are investigated with ring polymer molecular dynamics (RPMD) on a highly accurate full-dimensional potential energy surface. Kinetic isotope effects (KIEs) for three isotopologues of the title reaction are also studied. The results demonstrate excellent agreement with experimental data, even in the deep tunneling region. Especially, the Cayley propagator shows a high calculation efficiency with little loss of accuracy. The present results confirmed the applicability of the RPMD method, particularly the speed-up using a Cayley propagator, in theoretical calculations of bimolecular reaction rates.

Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study

We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.

Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches

The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.

Adiabatic models for the quantum dynamics of surface scattering with lattice effects

In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.

First-Principles Theory for the Reaction of Chlorine with Methane

A full-dimensional quantum dynamics simulation of the Cl + CH₄ → HCl + CH₃ reaction based on first-principles theory is reported. Accurate thermal rate constants are calculated, and perfect agreement with experiment is obtained. Despite the heavy atoms present in both reactants, the passage of the reaction barrier is found to occur within only a few tens of femtoseconds. This surprisingly short time scale results from correlated motion of the transferring hydrogen atom and the hydrogen atoms in the methyl fragment which facilitates irreversible barrier passage without relevant participation of heavy atoms. Resonance effects resulting from the heavy-light-heavy characteristics of the reaction system, which were observed in reactive scattering studies, do not affect the thermal rate constant.

Strong non-Arrhenius behavior at low temperatures in the OH + HCl → H2O + Cl reaction due to resonance induced quantum tunneling

The OH + HCl reaction possesses many Feshbach resonances trapped in the hydrogen bond well in the entrance channel, which substantially enhance the reaction rates at low temperatures.

Isotopic separation of helium through graphyne membranes: a ring polymer molecular dynamics study

Microscopic-level understanding of the separation mechanism for two-dimensional (2D) membranes is an active area of research due to potential implications of this class of membranes for various technological processes. Helium (He) purification from the natural resources is of particular interest due to the shortfall in its production. In this work, we applied the ring polymer molecular dynamics (RPMD) method to graphdiyne (Gr2) and graphtriyne (Gr3) 2D membranes having variable pore sizes for the separation of He isotopes, and compare for the first time with rigorous quantum calculations. We found that the transmission rate through Gr3 is many orders of magnitude greater than Gr2. The selectivity of either isotope at low temperatures is a consequence of a delicate balance between the zero-point energy effect and tunneling of ⁴He and ³He. In particular, a remarkable tunneling effect is reported on the Gr2 membrane at 10 K, leading to a much larger permeation of the lighter species as compared to the heavier isotope. RPMD provides an efficient approach for studying the separation of He isotopes, taking into account quantum effects of light nuclei motions at low temperatures, which classical methods fail to capture.

Dynamics and kinetics of the Si(1D) + H2/D2 reactions on a new global ab initio potential energy surface

Recent studies on the exothermic complex-forming reactions have improved our understanding of complex-forming reactions greatly, however, so far a similar level of study on endothermic ones has been rather limited. In this work, the endothermic complex-forming reaction Si(1D) + H2 → SiH + H and its deuterated isotopic variant are investigated by quantum dynamics (QD) and ring polymer molecular dynamics (RPMD) calculations on a new global ab initio potential energy surface (PES) for the ground electronic state, which is constructed based on 8996 symmetry unique points computed at the icMRCI+Q/aug-cc-pVQZ level. The PES reproduces our ab initio data very well in the dynamically important regions, on which the ro-vibrational energy levels of SiH2 are calculated and general good agreement with experiment is obtained. The integral cross sections and product angular and state distributions are computed in a wide range of collision energies, and interesting dynamics behaviors are revealed. The rate coefficients are also investigated, which display an exponential rise from 2.09 × 10-20 to 6.00 × 10-12 cm3 s-1 for the Si(1D) + H2 reaction as the temperature increases from 300 to 1500 K, in contrast to the nearly temperature-independent behavior of exothermic complex-forming reactions. In addition, the applicability of the RPMD approach is demonstrated, and the kinetic isotope effect is investigated, the ratio of which decreases from 7.89 (300 K) to 1.70 (1500 K). The effects of tunneling and initial rotational excitation are also discussed.

Theoretical H + O3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions

Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.

Kinetic study of the OH + HO2 → H2O + O2 reaction using ring polymer molecular dynamics and quantum dynamics

The reaction OH + HO2 → H2O + O2 is a prototype of radical-radical reactions. It plays an important role in interstellar/atmospheric chemistry and combustion, and considerable attention has thus been dedicated to its kinetics. In our previous work, we reported an accurate full-dimensional potential energy surface for the title reaction on the ground triplet electronic state. The quasi-classical trajectory (QCT) approach was employed to investigate its kinetics. Although the QCT rate coefficients were in good agreement with some experimental and theoretical results, QCT cannot account for the quantum mechanical effects, such as zero-point vibrational energy, recrossing, and tunneling, which may significantly affect the rate coefficients, particularly at low temperatures. In this work, the reduced-dimensional quantum dynamics and ring polymer molecular dynamics calculations were carried out to examine these effects and their impact on rate coefficients over the temperature range of 300-1300 K.

VTST and RPMD kinetics study of the nine-body X + C2H6 (X ≡ H, Cl, F) reactions based on analytical potential energy surfaces

Thermal rate constants of nine-atom hydrogen abstraction reactions, X + C2H6 → HX + C2H5 (X ≡ H, Cl, F) with qualitatively different reaction paths, have been investigated using two kinetics approaches - variational transition state theory with multidimensional tunnelling (VTST/MT) and ring polymer molecular dynamics (RPMD) - and full dimensional analytical potential energy surfaces. For the H + C2H6 reaction, which proceeds through a noticeable barrier height of 11.62 kcal mol-1, kinetics approaches showed excellent agreement between them (with differences less than 30%) and with the experiment (with differences less than 60%) in the wide temperature range of 200-2000 K. For X = Cl and F, however, the situation is very different. The barrier height is either low or very low, 2.44 and 0.23 kcal mol-1, respectively, and the presence of van der Waals complexes in the entrance channel leads to a very flat topography and, consequently, imposes theoretical challenges. For the Cl(2P) reaction, VTST/MT underestimates the experimental rate constants (with differences less than 86%), and RPMD demonstrates better agreement (with differences less than 47%), although the temperature dependence is opposite to the experiment at low temperatures. Finally, for the F(2P) reaction, available experimental information shows discrepancies, both in the absolute values of the rate constants and also in the temperature dependence. Unfortunately, kinetics theories did not resolve this discrepancy. Different possible causes of these theory/experiment discrepancies were analyzed.

An accurate potential energy surface and ring polymer molecular dynamics study of the Cl + CH4→ HCl + CH3reaction

Thermal rate coefficients for the Cl + CH₄/CD₄reactions were studied on a new full-dimensional accurate potential energy surface with the spin–orbit corrections considered in the entrance channel.

Path integral methods for reaction rates in complex systems

We shall use this introduction to the Faraday Discussion on quantum effects in complex systems to review the recent progress that has been made in using imaginary time path integral methods to calculate chemical reaction rates.

Ring Polymer Molecular Dynamics in Gas-Surface Reactions: Inclusion of Quantum Effects Made Simple

Accurately modeling gas-surface collision dynamics presents a great challenge for theory, especially in the low-energy (or temperature) regime where quantum effects are important. Here, a path integral-based nonequilibrium ring polymer molecular dynamics (NE-RPMD) approach is adapted to calculate dissociative initial sticking probabilities (S₀) of H₂ on Cu(111) and D₂O on Ni(111), revealing the distinct quantum nature in the two benchmark surface reactions. NE-RPMD successfully captures quantum tunneling in H₂ dissociation at very low energies, where the quasi-classical trajectory (QCT) method suddenly fails. Additionally, QCT substantially overestimates S₀ of D₂O because of severe zero point energy (ZPE) leakage, even at collision energies greater than the ZPE-corrected barrier. Instead, NE-RPMD predicts S₀ values of D₂O in much improved agreement with reference results obtained by the quantum wavepacket method with reasonable corrections of the thermal contribution. Our results suggest NE-RPMD as a promising approach to model quantum effects in gas-surface reactions.

Dynamics of H + HeH+(v = 0, j = 0) → H2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction

The H + HeH⁺→ He + H₂⁺ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.

New Stress Test for Ring Polymer Molecular Dynamics: Rate Coefficients of the O(3P) + HCl Reaction and Comparison with Quantum Mechanical and Quasiclassical Trajectory Results

In the past decade, ring polymer molecular dynamics (RPMD) has emerged as a very efficient method to determine thermal rate coefficients for a great variety of chemical reactions. This work presents the application of this methodology to study the O(³P) + HCl reaction, which constitutes a stringent test for any dynamical calculation due to rich resonant structure and other dynamical features. The rate coefficients, calculated on the ³A' and ³A″ potential energy surfaces (PESs) by Ramachandran and Peterson [ J. Chem. Phys. 2003 , 119 , 9590 ], using RPMD and quasiclassical trajectories (QCT) are compared with the existing experimental and the quantum mechanical (QM) results by Xie et al. [ J. Chem. Phys. 2005 122 , 014301 ]. The agreement is very good at T > 600 K, although RPMD underestimates rate coefficients by a factor between 4 and 2 in the 200-500 K interval. The origin of these discrepancies lies in the large contribution from tunneling on the ³A″ PES, which is enhanced by resonances due to quasibound states in the van der Waals wells. Although tunneling is fairly well accounted for by RPMD even below the crossover temperature, the effect of resonances, a long-time effect, is not included in the methodology. At the highest temperatures studied in this work, 2000-3300 K, the RPMD rate coefficients are somewhat larger than the QM ones, but this is shown to be due to limitations in the QM calculations and the RPMD are believed to be more reliable.

The low temperature D+ + H2→ HD + H+ reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study

The reaction between D⁺ and H₂ plays an important role in astrochemistry at low temperatures and also serves as a prototype for a simple ion-molecule reaction. Its ground X[combining tilde]¹A' state has a very small thermodynamic barrier (up to 1.8 × 10^-2 eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well with a depth of at least 0.2 eV, thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20-100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34-2.01 × 10^-9 cm³ s^-1. They are also in very good agreement with previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20-50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that the experimental rate measurements for the title reaction are highly desirable in future.

Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy

The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.

Theoretical description of quantum mechanical permeation of graphene membranes by charged hydrogen isotopes

It has been recently shown that in the presence of an applied voltage, hydrogen and deuterium nuclei can be separated from one another using graphene membranes as a nuclear sieve, resulting in a 10-fold enhancement in the concentration of the lighter isotope. While previous studies, both experimental and theoretical, have attributed this effect mostly to differences in vibrational zero point energy (ZPE) of the various isotopes near the membrane surface, we propose that multi-dimensional quantum mechanical tunneling of nuclei through the graphene membrane influences this proton permeation process in a fundamental way. We perform ring polymer molecular dynamics calculations in which we include both ZPE and tunneling effects of various hydrogen isotopes as they permeate the graphene membrane and compute rate constants across a range of temperatures near 300 K. While capturing the experimentally observed separation factor, our calculations indicate that the transverse motion of the various isotopes across the surface of the graphene membrane is an essential part of this sieving mechanism. An understanding of the multi-dimensional quantum mechanical nature of this process could serve to guide the design of other such isotopic enrichment processes for a variety of atomic and molecular species of interest.

Experimental and theoretical studies of the reactions of ground-state sulfur atoms with hydrogen and deuterium

The gas-phase kinetics of S(³P) atoms with H₂ and D₂ have been studied via the laser flash photolysis-resonance fluorescence technique. S atoms were generated by pulsed photolysis of CS₂ at 193 nm and monitored by time-resolved fluorescence at 181 nm. The rate coefficients for H₂ (k₁) and D₂ (k₂), respectively, are summarized as k₁(600-1110 K) = 3.0 × 10^-9 exp-1.317×10⁵-2.703×10⁷K/T8.314 T/K cm³ molecule^-1 s^-1 and k₂(770-1110 K) = 2.2 × 10^-14 (T/298 K)^3.55 exp(-5420 K/T) cm³ molecule^-1 s^-1. Error limits are discussed in the text. The rate coefficients for formation of SH(SD) + H(D) on a newly developed triplet potential energy surface were characterized via ring polymer molecular dynamics and canonical variational transition-state theory. There is excellent agreement above about 1000 K between theory and experiment. At lower temperatures, the experimental rate coefficient is substantially larger than the results computed for the adiabatic reaction, suggesting a significant role for intersystem crossing to the singlet potential energy surface at lower temperatures.

A combined theoretical and experimental investigation of the kinetics and dynamics of the O(1D) + D2 reaction at low temperature

The O(¹D) + H₂ reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(¹D) + D₂ reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.

Ab initio instanton rate theory made efficient using Gaussian process regression

Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into the calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but it is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but it requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

Accurate Determination of Tunneling-Affected Rate Coefficients: Theory Assessing Experiment

The thermal rate coefficients of a prototypical bimolecular reaction are determined on an accurate ab initio potential energy surface (PES) using ring polymer molecular dynamics (RPMD). It is shown that quantum effects such as tunneling and zero-point energy (ZPE) are of critical importance for the HCl + OH reaction at low temperatures, while the heavier deuterium substitution renders tunneling less facile in the DCl + OH reaction. The calculated RPMD rate coefficients are in excellent agreement with experimental data for the HCl + OH reaction in the entire temperature range of 200-1000 K, confirming the accuracy of the PES. On the other hand, the RPMD rate coefficients for the DCl + OH reaction agree with some, but not all, experimental values. The self-consistency of the theoretical results thus allows a quality assessment of the experimental data.