Electronic Structure, Stability, and Electrical Mobility of Cationic Silver Oxide Atomic Clusters.

Silver oxide cluster cations (AgnOm+) can readily be produced by a number of methods including atmospheric-pressure spark ablation of pure silver electrodes when trace amounts of oxygen are present in the carrier gas. Here we determine the equilibrium geometries of AgnOm+ clusters (n = 1-4; m = 1-5) using accurate coupled cluster with singles and doubles (CCSD) method, while the stabilization energies are calculated with additional perturbative triples correction (CCSD(T)). Although a number of stable states have been identified, our results show that the AgnOm+ clusters with m = 1 are more stable than those with m ≥ 2 due to the absence of the terminally attached O2 molecule, corroborating recent observations by mass spectrometry. Using the computed structures, we calculate the electrical mobilities of the AgnOm+ clusters and label the values on a respective experimentally determined spectrum in an attempt to better interpret the occurrence of the peaks and troughs in the measurements.

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 4He and 3He. 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.

Experimental and theoretical studies of the gas-phase reactions of O(1D) with H2O and D2O at low temperature.

Here we report the results of an experimental and theoretical study of the gas-phase reactions between O(1D) and H2O and O(1D) and D2O at room temperature and below. On the experimental side, the kinetics of these reactions have been investigated over the 50-127 K range using a continuous flow Laval nozzle apparatus, coupled with pulsed laser photolysis and pulsed laser induced fluorescence for the production and detection of O(1D) atoms respectively. Experiments were also performed at 296 K in the absence of a Laval nozzle. On the theoretical side, the existing full-dimensional ground X 1A potential energy surface for the H2O2 system involved in this process has been reinvestigated and enhanced to provide a better description of the barrierless H-atom abstraction pathway. Based on this enhanced potential energy surface, quasiclassical trajectory calculations and ring polymer molecular dynamics simulations have been performed to obtain low temperature rate constants. The measured and calculated rate constants display similar behaviour above 100 K, showing little or no variation as a function of temperature. Below 100 K, the experimental rate constants increase dramatically, in contrast to the essentially temperature independent theoretical values. The possible origins of the divergence between experiment and theory at low temperatures are discussed.

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.

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

The H + HeH+→ He + H2+ 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.

Experimental and Theoretical Study of the O(1D) + HD Reaction.

This work addresses the kinetics and dynamics of the gas-phase reaction between O(1D) and HD molecules down to low temperature. Here, measurements were performed by using a supersonic flow (Laval nozzle) reactor coupled with pulsed laser photolysis for O(1D) production and pulsed-laser-induced fluorescence for O(1D) detection to obtain rate constants over the 50-300 K range. Additionally, temperature-dependent branching ratios (OD + H/OH + D) were obtained experimentally by comparison of the H/D atom atom yields with those of a reference reaction. In parallel, theoretical rate constants and branching ratios were calculated by using three different techniques; mean potential phase space theory (MPPST), the statistical quantum mechanical method (SQM), and ring polymer molecular dynamics (RPMD). Although the agreement between experimental and theoretical rate constants is reasonably good, with differences not exceeding 30% over the entire temperature range, the theoretical branching ratios derived by the MPPST and SQM methods are as much as 50% larger than the experimental ones. These results are presented in the context of earlier work, while the possible origins of the discrepancies between experiment and theory are discussed.

Ring polymer molecular dynamics and active learning of moment tensor potential for gas-phase barrierless reactions: Application to S + H2.

Ring polymer molecular dynamics (RPMD) has proven to be an accurate approach for calculating thermal rate coefficients of various chemical reactions. For wider application of this methodology, efficient ways to generate the underlying full-dimensional potential energy surfaces (PESs) and the corresponding energy gradients are required. Recently, we have proposed a fully automated procedure based on combining the original RPMDrate code with active learning for PES on-the-fly using moment tensor potential and successfully applied it to two representative thermally activated chemical reactions [I. S. Novikov et al., Phys. Chem. Chem. Phys. 20, 29503-29512 (2018)]. In this work, using a prototype insertion chemical reaction S + H2, we show that this procedure works equally well for another class of chemical reactions. We find that the corresponding PES can be generated by fitting to less than 1500 automatically generated structures, while the RPMD rate coefficients show deviation from the reference values within the typical convergence error of the RPMDrate. We note that more structures are accumulated during the real-time propagation of the dynamic factor (the recrossing factor) as opposed to the previous study. We also observe that a relatively flat free energy profile along the reaction coordinate before entering the complex-formation well can cause issues with locating the maximum of the free energy surface for less converged PESs. However, the final RPMD rate coefficient is independent of the position of the dividing surface that makes it invulnerable to this problem, keeping the total number of necessary structures within a few thousand. Our work concludes that, in the future, the proposed methodology can be applied to realistic complex chemical reactions with various energy profiles.

Unimolecular Reaction Pathways of a γ-Ketohydroperoxide from Combined Application of Automated Reaction Discovery Methods.

Ketohydroperoxides are important in liquid-phase autoxidation and in gas-phase partial oxidation and pre-ignition chemistry, but because of their low concentration, instability, and various analytical chemistry limitations, it has been challenging to experimentally determine their reactivity, and only a few pathways are known. In the present work, 75 elementary-step unimolecular reactions of the simplest γ-ketohydroperoxide, 3-hydroperoxypropanal, were discovered by a combination of density functional theory with several automated transition-state search algorithms: the Berny algorithm coupled with the freezing string method, single- and double-ended growing string methods, the heuristic KinBot algorithm, and the single-component artificial force induced reaction method (SC-AFIR). The present joint approach significantly outperforms previous manual and automated transition-state searches - 68 of the reactions of γ-ketohydroperoxide discovered here were previously unknown and completely unexpected. All of the methods found the lowest-energy transition state, which corresponds to the first step of the Korcek mechanism, but each algorithm except for SC-AFIR detected several reactions not found by any of the other methods. We show that the low-barrier chemical reactions involve promising new chemistry that may be relevant in atmospheric and combustion systems. Our study highlights the complexity of chemical space exploration and the advantage of combined application of several approaches. Overall, the present work demonstrates both the power and the weaknesses of existing fully automated approaches for reaction discovery which suggest possible directions for further method development and assessment in order to enable reliable discovery of all important reactions of any specified reactant(s).

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 H2 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]1A' 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 cm3 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.

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

The O(1D) + H2 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(1D) + D2 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.

SCOPA and META-SCOPA: software for the analysis and aggregation of genome-wide association studies of multiple correlated phenotypes.

BACKGROUND: Genome-wide association studies (GWAS) of single nucleotide polymorphisms (SNPs) have been successful in identifying loci contributing genetic effects to a wide range of complex human diseases and quantitative traits. The traditional approach to GWAS analysis is to consider each phenotype separately, despite the fact that many diseases and quantitative traits are correlated with each other, and often measured in the same sample of individuals. Multivariate analyses of correlated phenotypes have been demonstrated, by simulation, to increase power to detect association with SNPs, and thus may enable improved detection of novel loci contributing to diseases and quantitative traits. RESULTS: We have developed the SCOPA software to enable GWAS analysis of multiple correlated phenotypes. The software implements "reverse regression" methodology, which treats the genotype of an individual at a SNP as the outcome and the phenotypes as predictors in a general linear model. SCOPA can be applied to quantitative traits and categorical phenotypes, and can accommodate imputed genotypes under a dosage model. The accompanying META-SCOPA software enables meta-analysis of association summary statistics from SCOPA across GWAS. Application of SCOPA to two GWAS of high-and low-density lipoprotein cholesterol, triglycerides and body mass index, and subsequent meta-analysis with META-SCOPA, highlighted stronger association signals than univariate phenotype analysis at established lipid and obesity loci. The META-SCOPA meta-analysis also revealed a novel signal of association at genome-wide significance for triglycerides mapping to GPC5 (lead SNP rs71427535, p = 1.1x10-8), which has not been reported in previous large-scale GWAS of lipid traits. CONCLUSIONS: The SCOPA and META-SCOPA software enable discovery and dissection of multiple phenotype association signals through implementation of a powerful reverse regression approach.

Kinetics study of the CN + CH4 hydrogen abstraction reaction based on a new ab initio analytical full-dimensional potential energy surface.

We have developed an analytical full-dimensional potential energy surface, named PES-2017, for the gas-phase hydrogen abstraction reaction between the cyano radical and methane. This surface is fitted using high-level ab initio information as input. Using the PES-2017 surface, a kinetics study was performed via two theoretical approaches: variational transition-state theory with multidimensional tunnelling (VTST-MT) and ring polymer molecular dynamics (RPMD). The results are compared with the experimental data. In the whole temperature range analysed, 300-1500 K, both theories agree within a factor of <2, reproducing the experimental behaviour taking into account the experimental uncertainties. At high temperatures, where the recrossing effects dominate and the RPMD theory is exact, both theories differ by a factor of about 20%; while at low temperatures this difference is larger, 45%. Note that in this temperature regime, the tunnelling effect is negligible. The CN + CH4/CD4 kinetic isotope effects are important, reproducing the scarce experimental evidence. The good agreement with the ab initio information used in the fitting process (self-consistency test) and with the kinetic behaviour in a wide temperature range gives confidence and strength to the new surface.

Low-Temperature Experimental and Theoretical Rate Constants for the O(1D) + H2 Reaction.

In the present joint experimental and theoretical study, we report thermal rate constants for the O(1D) + H2 reaction within the 50-300 K temperature range. Experimental kinetics measurements were performed using a continuous supersonic flow reactor coupled with pulsed laser photolysis for O(1D) production and pulsed laser-induced fluorescence in the vacuum ultraviolet wavelength range (VUV LIF) for O(1D) detection. Theoretical rate constants were obtained using the ring polymer molecular dynamics (RPMD) approach over the two lowest potential energy surfaces 11A' and 11A″, which possess barrierless and thermally activated energy profiles, respectively. Both the experimental and theoretical rate constants exhibit a weak temperature dependence. The theoretical results show the dominant role of the 11A' ground state and that contribution of the 11A″ excited state to the total thermal rate decreases dramatically at lower temperature. Agreement between the experimental and theoretical results is good, and the discrepancy does not exceed 25%. It is argued that these differences are likely to be due to nonadiabatic couplings between the 11A' and 21A' surfaces.

An experimental and theoretical investigation of the C(1D) + D2 reaction.

In a previous joint experimental and theoretical study of the barrierless chemical reaction C(1D) + H2 at low temperatures (300-50 K) [K. M. Hickson, J.-C. Loison, H. Guo and Y. V. Suleimanov, J. Phys. Chem. Lett., 2015, 6, 4194], excellent agreement was found between experimental thermal rate constants and theoretical estimates based on ring polymer molecular dynamics (RPMD) over the two lowest singlet potential energy surfaces (PESs). Here, we extend this work to one of its isotopologues, C(1D) + D2, over the same temperature range. Experimental and RPMD results are in good agreement when contributions from both PESs to this chemical reaction are included in the RPMD simulations. The deviation between experiment and the RPMD calculations does not exceed 25% and both results exhibit a slight negative temperature dependence. The first excited 1A'' PES plays a more important role than the ground 1A' PES as the temperature is decreased, similar to our previous studies of the C(1D) + H2 reaction but with a more pronounced effect. The small differences in temperature dependence between the earlier and present experimental studies of C(1D) + H2/D2 reactions are discussed in terms of the use of non-equilibrium populations of ortho/para-H2/D2. We argue that RPMD provides a very convenient and reliable tool to study low-temperature chemical reactions.

Thermal Rate Coefficients for the Astrochemical Process C + CH+ → C2+ + H by Ring Polymer Molecular Dynamics.

Thermal rate coefficients for the astrochemical reaction C + CH+ → C2+ + H were computed in the temperature range 20-300 K by using novel rate theory based on ring polymer molecular dynamics (RPMD) on a recently published bond-order based potential energy surface and compared with previous Langevin capture model (LCM) and quasi-classical trajectory (QCT) calculations. Results show that there is a significant discrepancy between the RPMD rate coefficients and the previous theoretical results that can lead to overestimation of the rate coefficients for the title reaction by several orders of magnitude at very low temperatures. We argue that this can be attributed to a very challenging energy profile along the reaction coordinate for the title reaction, not taken into account in extenso by either the LCM or QCT approximation. In the absence of any rigorous quantum mechanical or experimental results, the computed RPMD rate coefficients represent state-of-the-art estimates to be included in astrochemical databases and kinetic networks.

Chemical Reaction Rate Coefficients from Ring Polymer Molecular Dynamics: Theory and Practical Applications.

This Feature Article presents an overview of the current status of ring polymer molecular dynamics (RPMD) rate theory. We first analyze the RPMD approach and its connection to quantum transition-state theory. We then focus on its practical applications to prototypical chemical reactions in the gas phase, which demonstrate how accurate and reliable RPMD is for calculating thermal chemical reaction rate coefficients in multifarious cases. This review serves as an important checkpoint in RPMD rate theory development, which shows that RPMD is shifting from being just one of recent novel ideas to a well-established and validated alternative to conventional techniques for calculating thermal chemical rate coefficients. We also hope it will motivate further applications of RPMD to various chemical reactions.

Cross-Aldol Reaction of Isatin with Acetone Catalyzed by Leucinol: A Mechanistic Investigation.

Comprehensive mechanistic studies on the enantioselective aldol reaction between isatin (1 a) and acetone, catalyzed by L-leucinol (3 a), unraveled that isatin, apart from being a substrate, also plays an active catalytic role. Conversion of the intermediate oxazolidine 4 into the reactive syn-enamine 6, catalyzed by isatin, was identified as the rate-determining step by both the calculations (ΔG(≠) =26.1 kcal mol(-1) for the analogous L-alaninol, 3 b) and the kinetic isotope effect (kH /kD =2.7 observed for the reaction using [D6 ]acetone). The subsequent reaction of the syn-enamine 6 with isatin produces (S)-2 a (calculated ΔG(≠) =11.6 kcal mol(-1) ). The calculations suggest that the overall stereochemistry is controlled by two key events: 1) the isatin-catalyzed formation of the syn-enamine 6, which is thermodynamically favored over its anti-rotamer 7 by 2.3 kcal mol(-1) ; and 2) the high preference of the syn-enamine 6 to produce (S)-2 a on reaction with isatin (1 a) rather than its enantiomer (ΔΔG(≠) =2.6 kcal mol(-1) ).

Should thermostatted ring polymer molecular dynamics be used to calculate thermal reaction rates?

We apply Thermostatted Ring Polymer Molecular Dynamics (TRPMD), a recently proposed approximate quantum dynamics method, to the computation of thermal reaction rates. Its short-time transition-state theory limit is identical to rigorous quantum transition-state theory, and we find that its long-time limit is independent of the location of the dividing surface. TRPMD rate theory is then applied to one-dimensional model systems, the atom-diatom bimolecular reactions H + H2, D + MuH, and F + H2, and the prototypical polyatomic reaction H + CH4. Above the crossover temperature, the TRPMD rate is virtually invariant to the strength of the friction applied to the internal ring-polymer normal modes, and beneath the crossover temperature the TRPMD rate generally decreases with increasing friction, in agreement with the predictions of Kramers theory. We therefore find that TRPMD is approximately equal to, or less accurate than, ring polymer molecular dynamics for symmetric reactions, and for certain asymmetric systems and friction parameters closer to the quantum result, providing a basis for further assessment of the accuracy of this method.

Quantum rate coefficients and kinetic isotope effect for the reaction Cl + CH4 → HCl + CH3 from ring polymer molecular dynamics.

Thermal rate coefficients and kinetic isotope effect have been calculated for prototypical heavy-light-heavy polyatomic bimolecular reactions Cl + CH4/CD4 → HCl/DCl + CH3/CD3, using a recently proposed quantum dynamics approach: ring polymer molecular dynamics (RPMD). Agreement with experimental rate coefficients, which are quite scattered, is satisfactory. However, differences up to 50% have been found between the RPMD results and those obtained from the harmonic variational transition-state theory on one of the two full-dimensional potential energy surfaces used in the calculations. Possible reasons for such discrepancy are discussed. The present work is an important step in a series of benchmark studies aimed at assessing accuracy for RPMD for chemical reaction rates, which demonstrates that this novel method is a quite reliable alternative to previously developed techniques based on transition-state theory.

Theoretical kinetics study of the O(³P) + CH4/CD4 hydrogen abstraction reaction: the role of anharmonicity, recrossing effects, and quantum mechanical tunneling.

Using a recently developed full-dimensional accurate analytical potential energy surface [Gonzalez-Lavado, E., Corchado, J. C., and Espinosa-Garcia, J. J. Chem. Phys. 2014, 140, 064310], we investigate the thermal rate coefficients of the O((3)P) + CH4/CD4 reactions with ring polymer molecular dynamics (RPMD) and with variational transition-state theory with multidimensional tunneling corrections (VTST/MT). The results of the present calculations are compared with available experimental data for a wide temperature range 200-2500 K. In the classical high-temperature limit, the RPMD results match perfectly the experimental data, whereas VTST results are smaller by a factor of 2. We suggest that this discrepancy is due to the harmonic approximation used in the present VTST calculations, which leads to an overestimation of the variational effects. At low temperatures the tunneling plays an important role, which is captured by both methods, although they both overestimate the experimental values. The analysis of the kinetic isotope effects shows a discrepancy between both approaches, with the VTST values smaller by a factor about 2 at very low temperatures. Unfortunately, no experimental results are available to shed any light on this comparison, which keeps it as an open question.