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.

Activation of Methane by Zr+: A Deep-Dive into the Potential Surface via Pressure- and Temperature-Dependent Kinetics with Statistical Modeling.

The kinetics of Zr+ + CH4 are measured using a selected-ion flow tube apparatus over the temperature range 300-600 K and the pressure range 0.25-0.60 Torr. Measured rate constants are small, never exceeding 5% of the Langevin capture value. Both collisionally stabilized ZrCH4+ and bimolecular ZrCH2+ products are observed. A stochastic statistical modeling of the calculated reaction coordinate is used to fit the experimental results. The modeling indicates that an intersystem crossing from the entrance well, necessary for the bimolecular product to be formed, occurs faster than competing isomerization and dissociation processes. That sets an upper limit on the lifetime of the entrance complex to crossing of 10-11 s. The endothermicity of the bimolecular reaction is derived to be 0.09 ± 0.05 eV, in agreement with a literature value. The observed ZrCH4+ association product is determined to be primarily HZrCH3+ not Zr+(CH4), indicating that bond activation has occurred at thermal energies. The energy of HZrCH3+ relative to separated reactants is determined to be -0.80 ± 0.25 eV. Inspection of the statistical modeling results under best-fit conditions reveals reaction dependences on impact parameter, translation energy, internal energy, and angular momentum. Reaction outcomes are heavily affected by angular momentum conservation. Additionally, product energy distributions are predicted.

Stereodynamical Control of Cold Collisions between Two Aligned D_{2} Molecules.

Resonant scattering of optically state-prepared and aligned molecules in the cold regime allows the most detailed interrogation and control of bimolecular collisions. This technique has recently been applied to collisions of two aligned ortho-D_{2} molecules prepared in the j=2 rotational level of the v=2 vibrational manifold using the Stark-induced adiabatic Raman passage technique. Here, we develop the theoretical formalism for describing four-vector correlations in collisions of two aligned molecules and apply our approach to state-prepared D_{2}(v=2,j=2)+D_{2}(v=2,j=2)→D_{2}(v=2,j=2)+D_{2}(v=2,j=0) collisions, making possible the simulations of the experimental results from first principles. Key features of the experimental angular distributions are reproduced and attributed primarily to a partial wave resonance with orbital angular momentum ℓ=4.

A Universal and Modular Scaffold for Heparanase Activatable Probes and Drugs.

Heparanase (HPSE) is an endo-ß-glucuronidase involved in extracellular matrix remodeling in rapidly healing tissues, most cancers and inflammation, and viral infection. Its importance as a therapeutic target warrants further study, but such is hampered by a lack of research tools. To expand the toolkits for probing HPSE enzymatic activity, we report the design of a substrate scaffold for HPSE comprised of a disaccharide substrate appended with a linker, capable of carrying cargo until being cleaved by HPSE. Here exemplified as a fluorogenic, coumarin-based imaging probe, this scaffold can potentially expand the availability of HPSE-responsive imaging or drug delivery tools using a variety of imaging moieties or other cargo. We show that electronic tuning of the scaffold provides a robust response to HPSE while simplifying the structural requirements of the attached cargo. Molecular docking and modeling suggest a productive probe/HPSE binding mode. These results further support the hypothesis that the reactivity of these HPSE-responsive probes is predominantly influenced by the electron density of the aglycone. This universal HPSE-activatable scaffold will greatly facilitate future development of HPSE-responsive probes and drugs.

Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory: The Lowest Doublet States Responsible for the N(4Su) + C2(a 3Πu) → CN(X 2Σ+) + C(3Pg) Reaction.

Global potential energy surfaces (PESs) for the 1 2A' and 1 2A″ states of the C2N system responsible for the N(4Su) + C2(a 3Πu) → CN(X 2Σ+) + C(3Pg) reaction are mapped using compressed-state multistate pair-density functional theory (CMS-PDFT), which is a multi-state version of multiconfiguration pair-density functional theory (MC-PDFT). Calculations are also performed at selected geometries by explicitly correlated multireference configuration interaction with quadruple corrections, MRCI-F12+Q, and the comparison of the two sets of calculations shows that CMS-PDFT describes the globally reactive PESs well, including the bond-breaking asymptotes. We conclude that CMS-PDFT is an efficient method for constructing PESs for strongly correlated reactive systems. The PESs for producing CN + C are found to be barrierless and proceed through intermediate complexes. The CMS-PDFT PESs were fitted with a neural network method, and quasiclassical trajectories were computed on the resulting analytic PESs. These trajectories predict that the reaction produces vibrationally excited CN.

Quantum effects in thermal reaction rates at metal surfaces.

There is wide interest in developing accurate theories for predicting rates of chemical reactions that occur at metal surfaces, especially for applications in industrial catalysis. Conventional methods contain many approximations that lack experimental validation. In practice, there are few reactions where sufficiently accurate experimental data exist to even allow meaningful comparisons to theory. Here, we present experimentally derived thermal rate constants for hydrogen atom recombination on platinum single-crystal surfaces, which are accurate enough to test established theoretical approximations. A quantum rate model is also presented, making possible a direct evaluation of the accuracy of commonly used approximations to adsorbate entropy. We find that neglecting the wave nature of adsorbed hydrogen atoms and their electronic spin degeneracy leads to a 10× to 1000× overestimation of the rate constant for temperatures relevant to heterogeneous catalysis. These quantum effects are also found to be important for nanoparticle catalysts.

Differential Cross Sections for Cold, State-to-State Spin-Orbit Changing Collisions of NO(v = 10) with Neon.

Inelastic scattering processes have proven a powerful means of investigating molecular interactions, and much current effort is focused on the cold and ultracold regime where quantum phenomena are clearly manifested. Studies of collisions of the open shell nitric oxide (NO) molecule have been central in this effort since the pioneering work of Houston and co-workers in the early 1990s. State-to-state scattering of vibrationally excited molecules in the cold regime introduces challenges that test the suitability of current theoretical methods for ab initio determination of intermolecular potentials, and concomitant electronically nonadiabatic processes raise the bar further. Here we report measurements of differential cross sections for state-to-state spin-orbit changing collisions of NO (v = 10, Ω″ = 1.5, and j″ = 1.5) with neon from 2.3 to 3.5 cm-1 collision energy using our recently developed near-copropagating beam technique. The experimental results are compared with those obtained from quantum scattering calculations on a high-level set of coupled cluster potential energy surfaces and are shown to be in good agreement. The theoretical results suggest that distinct backscattering in the 2.3 cm-1 case arises from overlapping resonances.

Dissection of the Multichannel Reaction O(3P) + C2H2: Differential Cross-Sections and Product Energy Distributions.

The O(3P) + C2H2 reaction plays an important role in hydrocarbon combustion. It has two primary competing channels: H + HCCO (ketenyl) and CO + CH2 (triplet methylene). To further understand the microscopic dynamic mechanism of this reaction, we report here a detailed quasi-classical trajectory study of the O(3P) + C2H2 reaction on the recently developed full-dimensional potential energy surface (PES). The entrance barrier TS1 is the rate-limiting barrier in the reaction. The translation of reactants can greatly promote reactivity, due to strong coupling with the reaction coordinate at TS1. The O(3P) + C2H2 reaction progress through a complex-forming mechanism, in which the intermediate HCCHO lives at least through the duration of a rotational period. The energy redistribution takes place during the creation of the long-lived high vibrationally (and rotationally) excited HCCHO in the reaction. The product energy partitioning of the two channels and CO vibrational distributions agree with experimental data, and the vibrational state distributions of all modes of products present a Boltzmann-like distribution.

Full-Dimensional Potential Energy Surface for Ro-vibrationally Inelastic Scattering between H2 Molecules.

We report a new full-dimensional potential energy surface (PES) for the inelastic scattering between ro-vibrationally excited H2 molecules. The new PES is based on 39,462 multi-reference configuration interaction points in dynamically relevant regions. The analytic form of the PES consists of a short-range term fit with the permutational invariant polynomial-neural network method and a long-range term with a physically correct asymptotic functional form accounting for both electrostatic and dispersion terms, which are connected smoothly with a switching function. The PES compares favorably with existing accurate PESs near the H2 equilibrium geometries but covers a much larger configuration space for H2 with up to 10 vibrational quanta. Full-dimensional quantum scattering calculations on the new PES reproduce the recent Stark-induced adiabatic Raman passage results for the HD(v = 1) + H2 scattering near 1 K, validating its accuracy. These calculations also revealed significant differences with existing PESs in describing scattering of vibrationally excited molecules, underscoring the ability of the new PES in handling such dynamics.

Full-Dimensional Global Potential Energy Surface for the KRb + KRb → K2Rb2* → K2 + Rb2 Reaction with Accurate Long-Range Interactions and Quantum Statistical Calculation of the Product State Distribution under Ultracold Conditions.

A full-dimensional global potential energy surface (PES) for the KRb + KRb → K2Rb2* → K2 + Rb2 reaction is reported based on high-level ab initio calculations. The short-range part of the PES is fit with the permutationally invariant polynomial-neural network method, while the long-range parts of the PES in both the reactant and product asymptotes are represented by an asymptotically correct form. The long- and short-range parts are connected with intermediate-range parts to make them smooth. Within a statistical quantum model, this PES reproduces both the measured loss rates of ultracold KRb molecules and the K2 and Rb2 product state distributions, underscoring the important role of tunneling in ultracold chemistry. The PES also correctly predicts the lifetime of the K2Rb2* intermediate complex within the Rice-Ramsperger-Kassel-Marcus limit. It thus provides a reliable platform for future dynamical studies of the prototypical reaction.

Rainbow scattering in rotationally inelastic collisions of HCl and H2.

We examine rotational transitions of HCl in collisions with H2 by carrying out quantum mechanical close-coupling and quasi-classical trajectory (QCT) calculations on a recently developed globally accurate full-dimensional ab initio potential energy surface for the H3Cl system. Signatures of rainbow scattering in rotationally inelastic collisions are found in the state resolved integral and differential cross sections as functions of the impact parameter (initial orbital angular momentum) and final rotational quantum number. We show the coexistence of distinct dynamical regimes for the HCl rotational transition driven by the short-range repulsive and long-range attractive forces whose relative importance depends on the collision energy and final rotational state, suggesting that the classification of rainbow scattering into rotational and l-type rainbows is effective for H2 + HCl collisions. While the QCT method satisfactorily predicts the overall behavior of the rotationally inelastic cross sections, its capability to accurately describe signatures of rainbow scattering appears to be limited for the present system.

Time-independent quantum theory on vibrational inelastic scattering between atoms and open-shell diatomic molecules: Applications to NO + Ar and NO + H scattering.

A full-dimensional rigorous quantum mechanical treatment of non-reactive inelastic scattering of an open-shell diatom [e.g., NO(2Π)] with a structureless and spinless atom is presented within the time-independent close-coupling framework. The inclusion of the diatomic vibrational degree of freedom allows the investigation of transitions between different vibrational manifolds, in addition to those between different rotational, spin-orbit, and Λ-doublet states. This method is applied to the scattering of vibrationally excited NO(2Π) with Ar and H (with its spin ignored). The former has negligible vibrational inelasticity, thanks to the weak interaction between the two collisional partners. This conclusion justifies the commonly used two-dimensional approximation in treating NO scattering with rare gas atoms. The latter, on the other hand, is shown to undergo significant vibrational relaxation, even in the ultra-cold regime, owing to a chemically bonded (HNO) complex on the lowest-lying singlet potential energy surfaces.

Theoretical Investigations of Rate Coefficients for H + O3 and HO2 + O Reactions on a Full-Dimensional Potential Energy Surface.

In this work, a machine learning method is used to construct a high-fidelity multichannel global reactive potential energy surface (PES) for the HO3 system from 21452 high-level ab initio calculations at the explicitly correlated multireference configuration interaction (MRCI-F12) level of theory. The permutation invariance of the PES with respect to the three identical oxygen atoms is enforced using permutation invariant polynomials (PIPs) in the input layer of a neural network (NN). This PIP-NN representation is highly faithful to the ab initio points, with a root-mean-square error of 0.20 kcal/mol. Using this PES, the kinetics of H + O3 → OH + O2 (R1) and HO2 + O → OH + O2 (R2) reactions were investigated using a quasi-classical trajectory method over a wide temperature range (200-2000 K). It was found that the calculated thermal rate coefficients of R1 and R2, exhibiting positive and negative temperature dependences, respectively, are in reasonably good agreement with most experimental measured values. These temperature dependences can be attributed to the presence and absence of an entrance channel potential barrier.

State-to-state scattering of highly vibrationally excited NO at broadly tunable energies.

Experimental developments continue to challenge the theoretical description of molecular interactions. One key arena in which these advances have taken place is in rotationally inelastic scattering. Electric fields have been used with great success to select the initial quantum state and slow molecules for scattering studies, revealing novel stereodynamics, diffraction oscillations and scattering resonances. These have enjoyed excellent agreement with quantum scattering calculations performed on state-of-the-art coupled-cluster potential energy surfaces. To date these studies have largely employed reactants in the ground vibrational state (v = 0) and the lowest low-field-seeking quantum state. Here we describe the use of stimulated emission pumping to prepare NO molecules in arbitrary single rotational and parity states of v = 10 for inelastic scattering studies. These are employed in a near-copropagating molecular beam geometry that permits the collision energy to be tuned from above room temperature to 1 K or below, with product differential cross-sections obtained by velocity map imaging. This extremely nonequilibrium condition, not found in nature, tests current theoretical methods in a new regime.

A Global Full-Dimensional Potential Energy Surface for the K2Rb2 Complex and Its Lifetime.

A full-dimensional global potential energy surface for the KRb + KRb → K2 + Rb2 reaction is developed from 20 759 ab initio points calculated using a coupled cluster singles, doubles, and perturbative triples (CCSD(T)) method with effective core potentials, extrapolated to the complete basis set limit. The ab initio points are represented with high fidelity (root-mean-square error of 1.86 cm-1) using the permutation-invariant polynomial-neural network method, which enforces the permutation invariance of the potential with respect to exchange of identical nuclei. The potential energy surface features two D2h minima and one Cs minimum connected by the isomerization saddle points. The Rice-Ramsperger-Kassel-Marcus lifetime of the K2Rb2 reaction intermediate estimated using the potential energy surface is 227 ns, in reasonable agreement with the latest experimental measurement.

An ab initio based full-dimensional potential energy surface for OH + O2⇄ HO3 and low-lying vibrational levels of HO3.

To provide an in-depth understanding of the HO3 radical and its dissociation to OH + O2, a six-dimensional potential energy surface (PES) has been constructed by fitting 2087 energy points for the electronic ground state of HO3 (X2A'') using the permutation invariant polynomial-neural network (PIP-NN) approach. The energy points were calculated using an explicitly-correlated and Davidson-corrected multi-reference configuration interaction method with the correlation-consistent polarized valence double zeta basis (MRCI(Q)-F12/VDZ-F12). On the PES, the trans-HO3 isomer is found to be the global minimum, 33.0 cm-1 below the cis-HO3 conformer, which is consistent with previous high-level theoretical investigations. The dissociation to the OH + O2 asymptote from both conformers is shown to be barrierless. As a benchmark from a recently developed high-accuracy thermochemistry protocol, D0 for trans-HO3 is calculated to be 2.29 ± 0.36 kcal mol-1, only slightly deeper than the value of 2.08 kcal mol-1 obtained using the PES, and in reasonable agreement with the experimentally estimated value of 2.93 ± 0.07 kcal mol-1. Using this PES, low-lying vibrational energy levels of HO3 are determined using an exact quantum Hamiltonian and compared with available experimental results.

Dissection of the multichannel reaction of acetylene with atomic oxygen: from the global potential energy surface to rate coefficients and branching dynamics.

The O(3P) + C2H2 reaction is the first step in acetylene oxidation. The accurate kinetic data and the understanding of the reaction dynamics is of great importance. To this end, a full-dimensional global potential energy surface (PES) for the ground triplet state of the O(3P) + C2H2 reaction is constructed based on approximately 85 000 ab initio points calculated at the level of explicitly correlated unrestricted coupled cluster single, double, and perturbative triple excitations with the explicitly correlated polarized valence triple zeta basis set (UCCSD(T)-F12b/VTZ-F12). The PES is fit using the permutation invariant polynomial-neural network (PIP-NN) approach with a total root mean square error of 0.21 kcal mol-1. The key topographic features of the PES, including multiple potential wells and saddle points along different reaction pathways, are well represented by this fit PES. The kinetics and dynamics of the O(3P) + C2H2 reaction are investigated using the quasi-classical trajectory (QCT) method. The calculated rate coefficients are in good agreement with experimental data over a wide temperature range, especially when the temperature is lower than 1500 K. The product branch ratio has also been determined, which indicates the H + HCCO channel as the dominant reaction pathway at 298-3000 K, accounting for 80-90% of the overall rate coefficient, in agreement with experimental observations. The dynamics of the reaction is analyzed in detail.

A full-dimensional potential energy surface and quantum dynamics of inelastic collision process for H2-HF.

A full-dimensional ab initio potential energy surface for the H2-HF van der Waals complex was constructed by employing the coupled-cluster singles and doubles with noniterative inclusion of connected triples with augmented correlation-consistent polarised valence quadruple-zeta basis set plus bond functions. Using the improved coupled-states approximation including the nearest neighbor Coriolis couplings, we calculated the state-to-state scattering dynamics for pure rotational and ro-vibrational energy transfer processes. For pure rotational energy transfer, our results showed a different dynamical behavior for para-H2 and ortho-H2 in collision with hydrogen fluoride (HF), which is consistent with the previous study. Interestingly, some strong resonant peaks were presented in the cross sections for ro-vibrational energy transfer. In addition, the calculated vibrational-resolved rate constant is in agreement with the experimental results reported by Bott et al. These dynamics data can be further applied to the numerical simulation of HF chemical lasers.

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.