van der waals interactions in the Cl + HD reaction

Significant Isotope Effects from the Nonadiabatic Couplings in the Cl(2P) + HD(v = 0, j = 0) Reaction

The impact of non-Born-Oppenheimer couplings on the isotopic effects in the reaction of the Cl(2P) atom with the HD (v = 0, j = 0) molecule is investigated with our recently developed nonadiabatic time-independent quantum scattering methods, where the full open-shell characteristics are included in the six-state model, and also with the recently developed two-state model solving by time-independent methods, where part of the open-shell characteristic is included. The same reaction is also calculated with the simple adiabatic model using the lowest adiabatic potential energy surface. Compared with the results from different models, it is found that the reactivity of the Cl + HD → HCl + D channel is significantly overestimated in the adiabatic model. In contrast, the reactivity of the other channel agrees well with the nonadiabatic models. This is due to the van der Waals well in the reactant channel being changed a lot by including the nonadiabatic couplings. These quantum dynamics calculations suggest that sometimes the adiabatic model should be used with caution; otherwise, it may result in significant deviations for some reactions.

The Decay of Dispersion Interaction and Its Remarkable Effects on the Kinetics of Activation Reactions Involving Alkyl Chains

The importance of dispersion interactions in many chemical processes is well recognized. It is known that the dispersion strength would decay with the increasing separation between the interacting groups; however, its effects on chemical reactivity have not been well understood. Here we reveal the decay law of dispersion interactions along the n-alkyl chain, its effective interaction ranges for common functional groups, and their remarkable effects on the kinetics of activation reactions involving alkyl chains. This is achieved by DLPNO-CCSD(T) calculations and the local energy decomposition analysis and is supported by experimental findings. In particular, our calculations indicate that the lifetime of alkyl-substituted cis-azobenzenes increases with the alkyl chain length but reaches a steady value when alkyl chains are longer than butyl groups, which is in satisfactory accordance with experimental measurements. We also propose a concise expression to describe the dispersion decay, which shows excellent agreement with our computed results.

Accurate determination of reaction energetics and kinetics of the HO2˙ + O3 → OH˙ + 2O2 reaction

In the present work, we have studied the HO₂˙ + O₃ → HO˙ + 2O₂ reaction using chemical kinetics and quantum chemical calculations. We have employed the post-CCSD(T) method to estimate the barrier height and reaction energy for the title reaction. In the post-CCSD(T) method, we have included zero point energy corrections, contributions from full triple excitations and partial quadratic excitations at the coupled-cluster level, and core corrections. We have also computed the reaction rate in the temperature range of 197-450 K and found good agreement with all the available experimental results. In addition, we have also fitted the computed rate constants with the Arrhenius expression and obtained an activation energy of 1.0 ± 0.1 kcal mol^-1, almost identical to the value recommended by IUPAC and JPL.

Production of ultracold polyatomic molecules with strong polarity by laser cooling: A detailed theoretical study on CaNC and SrNC

Laser cooling molecules to the ultracold regime is the prerequisite for many novel science and technologies. It is desirable to take advantage of theoretical approaches to explore polyatomic molecular candidates, which are capable of being cooled to the ultracold regime. In this work, we explore two polyatomic candidates, CaNC and SrNC, which are suitable for laser cooling. These molecules possess impressively large permanent dipole moments (∼6 Debye), which is preferred for applications using an external electric field. High-level ab initio calculations are carried out to reveal electronic structures of these molecules, and the calculated spectroscopic constants agree very well with the available experimental data. For each molecule, the Franck-Condon factor matrix is calculated and shows a diagonal distribution. The radiative lifetimes for CaNC and SrNC are estimated to be 15.5 and 15.8 ns, respectively. Based upon the features of various electronic states and by choosing suitable spin-orbit states, we construct two feasible laser cooling schemes for the two molecules, each of which allows scattering nearly 10000 photons for direct laser cooling. These indicate that CaNC and SrNC are excellent ultracold polyatomic candidates with strong polarity.

Accurate Calculation of Rate Constant and Isotope Effect for the F + H2 Reaction by the Coupled 3D Time-Dependent Wave Packet Method on the Newly Constructed Ab Initio Ground Potential Energy Surface

We employ coupled three-dimensional (3D) time dependent wave packet formalism in hyperspherical coordinates for reactive scattering problem on the newly constructed ab initio calculated ground adiabatic potential energy surface for the F + H₂/D₂ reaction. The convergence profiles for various reactive channels are depicted at low collision energy regimes with respect to the total angular momentum (J) quantum numbers. For two different reactant diatomic molecules (H₂ and D₂) initially at their respective ground roto-vibrational state (v = 0, j = 0), calculated state-to-state as well as total integral cross sections as a function of collision energy, temperature dependent rate constants, and the kinetic isotope effect for various reactivity profiles of F + H₂ and F + D₂ reactions are presented along with previous theoretical and experimental results.

Mechanism and Kinetics of Methane Combustion. Part II: Potential Energy Surface for Hydrogen-Abstraction Reaction of CH4 + O(3P)

Methane combustion plays an important role in various fields such as combustion chemistry and atmospheric chemistry of the stratosphere. Highly accurate study of its initial reaction remains a key challenge. Here, through extensive studies with a state-of-the-art ab initio and neural network method, we present a potential energy surface of the O(³P) + CH₄ → OH + CH₃ reaction on the ground state 1³A and the first excited state 2³A. In this work, the energies of 10 167 points covering all important regions are obtained with state-averaged complete active space self-consistent field calculations and then fitted using the Levenberg-Marquardt algorithm with a root-mean-square error of 0.391 and 0.442 kcal/mol for the 1³A and 2³A states, respectively. This study explores the characteristics of the radical van der Waals (VdW) complex and reveals a detailed mechanism of the methane combustion initial reaction. Within the scope of this mechanism, this surface gives a fairly accurate description of the regions around the saddle point, conical intersection, and vdW wells in the entrance for efficient computational simulations. As a theoretical study on a prototypical polyatomic reaction, it is hopeful that this work will modify our understanding of the primary process in hydrocarbon combustion.

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.

Charge Transfer Processes for H + H2+ Reaction Employing Coupled 3D Wavepacket Approach on Beyond Born-Oppenheimer Based Ab Initio Constructed Diabatic Potential Energy Surfaces

The dynamics of the H + H₂⁺ reaction has been analyzed from the electronically first excited state of diabatic potential energy surfaces constructed by employing the Beyond Born-Oppenheimer theory [J. Chem. Phys. 2014, 141, 204306]. We have employed the coupled 3D time-dependent wavepacket formalism in hyperspherical coordinates for multisurface reactive scattering problems. To be specific, the charge transfer processes have been investigated extensively by calculating state-to-state as well as total reaction probabilities and integral cross sections, when the reaction process is initiated from the first excited electronic state (2¹A'). We have depicted the convergence profiles of reaction probabilities for the competing charge transfer processes, namely, reactive charge transfer (RCT) and nonreactive charge transfer (NRCT) processes for different total energies with respect to total angular momentum, J. Total and state-to-state integral cross sections are calculated as a function of total energy for the initial rovibrational state, namely, v = 0, j = 0 level of H₂⁺ (²Σ_g⁺) molecule and are compared with previous theoretical calculations. Finally, we have calculated temperature-dependent rate constants using our presently evaluated cross sections and compared their average with the experimentally measured one.

On the formation of van der Waals complexes through three-body recombination

In this work, we show that van der Waals molecules X-RG (where RG is the rare gas atom) may be created through direct three-body recombination collisions, i.e., X + RG + RG → X-RG + RG. In particular, the three-body recombination rate at temperatures relevant for buffer gas cell experiments is calculated via a classical trajectory method in hyperspherical coordinates [Pérez-Ríos et al., J. Chem. Phys. 140, 044307 (2014)]. As a result, it is found that the formation of van der Waals molecules in buffer gas cells (1 K ≲ T ≲ 10 K) is dominated by the long-range tail (distances larger than the LeRoy radius) of the X-RG interaction. For higher temperatures, the short-range region of the potential becomes more significant. Moreover, we notice that the rate of formation of van der Walls molecules is of the same order of the magnitude independent of the chemical properties of X. As a consequence, almost any X-RG molecule may be created and observed in a buffer gas cell under proper conditions.

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.

Probing the location of the unpaired electron in spin-orbit changing collisions of NO with Ar

Understanding the molecular forces that drive a reaction or scattering process lies at the heart of molecular dynamics. Here, we present a combined experimental and theoretical study of the spin-orbit changing scattering dynamics of oriented NO molecules with Ar atoms. Using our crossed molecular beam apparatus, we have recorded velocity-map ion images and extracted differential and integral cross sections of the scattering process in the side-on geometry. We observe an overall preference for collisions close to the N atom in the spin-orbit changing manifold, which is a direct consequence of the location of the unpaired electron on the potential energy surface. In addition, a prominent forward scattered feature is observed for intermediate, even rotational transitions when the atom approaches the molecule from the O-end. The appearance of this peak originates from an attractive well on the A' potential energy surface, which efficiently directs high impact parameter trajectories towards the region of high unpaired electron density near the N-end of the molecule. The ability to orient molecules prior to collision, both experimentally and theoretically, allows us to sample different regions of the potential energy surface(s) and unveil the associated collision pathways.

Roaming Reactions and Dynamics in the van der Waals Region

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

Quantum Dynamics Study of the C(1D) + HD Reaction on the ã1A' and b̃1A″ Potential Energy Surfaces

We present an in-depth theoretical study of the C(¹D) + HD (v = 0, j = 0) → CD (CH) (v', j') + H (D) reaction using a time-dependent wave packet method with full Coriolis coupling on the Zhang-Ma-Bian potential energy surfaces (PESs) recently constructed by our group. The integral cross sections (ICS), differential cross sections, CD/CH branching ratios, and product state distributions are calculated over a wide range of collision energies. We find that the vibrational branching ratio defined as ICS(v'=1)/ICS(v'=0) obtained from the b̃¹A″ PES is much smaller than that from the ã¹A' PES for both product channels, which may be attributed to the dynamical effects of the conical intersection regulated (CI-R) intermediate on the b̃¹A″ PES. The collision energy dependence of CD/CH branching ratios displays oscillatory structures, which may be caused by the resonance states supported by the wells on the PESs. The high-temperature rate coefficients are also obtained and compared with previous results. The role of the excited-state PESs is also discussed.

Universal crossed beam imaging studies of polyatomic reaction dynamics

Crossed-beam imaging studies of polyatomic reactions show surprising dynamics not anticipated by extrapolation from smaller model systems.

Stereodynamical control of product branching in multi-channel barrierless hydrogen abstraction of CH3OH by F

Comprehensive dynamical simulations of a prototypical multi-channel reaction on a globally accurate potential energy surface show that the non-statistical product branching is dictated by unique stereodynamics in the entrance channels.

Hydrogen abstraction from methanol (CH₃OH) by F atoms presents an ideal proving ground to investigate dynamics of multi-channel reactions, because two types of hydrogen can be abstracted from the methanol molecule leading to the HF + CH₃O and HF + CH₂OH products. Using the quasi-classical trajectory approach on a globally accurate potential energy surface based on high-level ab initio calculations, this work reports a comprehensive dynamical investigation of this multi-channel reaction, yielding measurable attributes including integral and differential cross sections, as well as branching ratios. It is shown that while complex-forming and direct mechanisms coexist at low collision energies, these barrierless reaction channels are dominated at high energies by the direct mechanism, in which the reaction is only possible for trajectories entering into the respective dynamical cones of acceptance. Perhaps more importantly, the non-statistical product branching is found to be dictated by unique stereodynamics in the entrance channels.

Isotope-selective chemistry in the Be+(2S1/2) + HOD → BeOD+/BeOH+ + H/D reaction

Low temperature reactions between laser-cooled Be+(2S1/2) ions and partially deuterated water (HOD) molecules have been investigated using an ion trap and interpreted with zero-point corrected quasi-classical trajectory calculations on a highly accurate global potential energy surface for the ground electronic state. Both product channels have been observed for the first time, and the branching to BeOD+ + H is found to be 0.58 ± 0.14. The experimental observation is reproduced by both quasi-classical trajectory and statistical calculations. Theoretical analyses reveal that the branching to the two product channels is largely due to the availability of open states in each channel.

Enhancing radical molecular beams by skimmer cooling

A high-intensity supersonic beam source has been a key component in studies of molecular collisions, molecule-surface interaction, chemical reactions, and precision spectroscopy. However, the molecular density available for experiments in a downstream science chamber is limited by skimmer clogging, which constrains the separation between a valve and a skimmer to at least several hundred nozzle diameters. A recent experiment (Sci. Adv., 2017, 3, e1602258) has introduced a new strategy to address this challenge: when a skimmer is cooled to a temperature below the freezing point of the carrier gas, skimmer clogging can be effectively suppressed. We go beyond this proof-of-principle work in several key ways. Firstly, we apply the skimmer cooling approach to discharge-produced radical and metastable beams entrained in a carrier gas. We also identify two different processes for skimmer clogging mitigation-shockwave suppression at temperatures around the carrier gas freezing point and diffusive clogging at even lower temperatures. With the carrier clogging removed, we now fully optimize the production of entrained species such as hydroxyl radicals, resulting in a gain of 30 in density over the best commercial devices. The gain arises from both clogging mitigation and favorable geometry with a much shorter valve-skimmer distance.

First-Principles Computed Rate Constant for the O + O2 Isotopic Exchange Reaction Now Matches Experiment

We show, by performing exact time-independent quantum molecular scattering calculations, that the quality of the ground electronic state global potential energy surface appears to be of utmost importance in accurately obtaining even as strongly averaged quantities as kinetic rate constants. The oxygen isotope exchange reaction, ¹⁸O + ³²O₂, motivated by the understanding of a complex long-standing problem of isotopic ozone anomalies in the stratosphere and laboratory experiments, is explored in this context. The thermal rate constant for this key reaction is now in quantitative agreement with all experimental data available to date. A significant recent progress at the frontier of three research domains, advanced electronic structure calculations, ultrasensitive spectroscopy, and quantum scattering calculations, has therefore permitted a breakthrough in the theoretical modeling of this crucial collision process from first principles.

Dynamical resonances in chemical reactions

The transition state is a key concept in the field of chemistry and is important in the study of chemical kinetics and reaction dynamics.

Universality and chaoticity in ultracold K+KRb chemical reactions

A fundamental question in the study of chemical reactions is how reactions proceed at a collision energy close to absolute zero. This question is no longer hypothetical: quantum degenerate gases of atoms and molecules can now be created at temperatures lower than a few tens of nanokelvin. Here we consider the benchmark ultracold reaction between, the most-celebrated ultracold molecule, KRb and K. We map out an accurate ab initio ground-state potential energy surface of the K₂Rb complex in full dimensionality and report numerically-exact quantum-mechanical reaction dynamics. The distribution of rotationally resolved rates is shown to be Poissonian. An analysis of the hyperspherical adiabatic potential curves explains this statistical character revealing a chaotic distribution for the short-range collision complex that plays a key role in governing the reaction outcome.

Studying chemical reactions near zero temperature in detail is challenging both in theory and practice. Here the authors report an explicit quantum mechanical study of the benchmark ultracold reaction between a K atom and a KRb molecule, important for future controlled chemistry experiments.

Competition between the H- and D-atom transfer channels in the H2O+ + HD reaction: reduced-dimensional quantum and quasi-classical studies

The ion-molecule reaction between a water cation and a hydrogen molecule has recently attracted considerable interest due to its importance in astrochemistry. In this work, the intramolecular isotope effect of the H₂O⁺ + HD reaction is investigated using a seven-dimensional initial state-selected time-dependent wave packet approach as well as a full-dimensional quasi-classical trajectory method on a full-dimensional ab initio global potential energy surface. The calculated branching ratios for the formation of H₃O⁺ and H₂DO⁺via H- and D-transfer agree reasonably well with the experimental values. The preference to the formation of the H₃O⁺ product observed using the experiment at low collision energies is reproduced by theoretical calculations and explained by a one-dimensional effective potential model.

A review of dynamical resonances in A + BC chemical reactions

The concept of the transition state has played an important role in the field of chemical kinetics and reaction dynamics. Reactive resonances in the transition-state region can dramatically enhance the reaction probability; thus investigation of the reactive resonances has attracted great attention from chemical physicists for many decades. In this review, we mainly focus on the recent progress made in probing the elusive resonance phenomenon in the simple A  +  BC reaction and understanding its nature, especially in the benchmark F/Cl  +  H₂ and their isotopic variants. The signatures of reactive resonances in the integral cross section, differential cross section (DCS), forward- and backward-scattered DCS, and anion photodetachment spectroscopy are comprehensively presented in individual prototype reactions. The dynamical origins of reactive resonances are also discussed in this review, based on information on the wave function in the transition-state region obtained by time-dependent quantum wave-packet calculations.

Dynamical importance of van der Waals saddle and excited potential surface in C(1D)+D2 complex-forming reaction

Encouraged by recent advances in revealing significant effects of van der Waals wells on reaction dynamics, many people assume that van der Waals wells are inevitable in chemical reactions. Here we find that the weak long-range forces cause van der Waals saddles in the prototypical C(¹D)+D₂ complex-forming reaction that have very different dynamical effects from van der Waals wells at low collision energies. Accurate quantum dynamics calculations on our highly accurate ab initio potential energy surfaces with van der Waals saddles yield cross-sections in close agreement with crossed-beam experiments, whereas the same calculations on an earlier surface with van der Waals wells produce much smaller cross-sections at low energies. Further trajectory calculations reveal that the van der Waals saddle leads to a torsion then sideways insertion reaction mechanism, whereas the well suppresses reactivity. Quantum diffraction oscillations and sharp resonances are also predicted based on our ground- and excited-state potential energy surfaces.

It is commonly held that van der Waals wells are inevitable in chemical reactions. Here, the authors show that weak van der Waals forces in the entrance channel of a prototypical complex-forming reaction cause a van der Waals saddle instead, with different dynamical effects from a well at low collision energies.

Quasiclassical trajectory study of the C(1D) + HD reaction

Isotopic branching ratios are investigated by detailed quasiclassical trajectory calculations on our recent singlet ground and excited potential energy surfaces.

Recent advances in quantum scattering calculations on polyatomic bimolecular reactions

This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.

Crossed beam polyatomic reaction dynamics: recent advances and new insights

This review summarizes the developments in polyatomic reaction dynamics, focusing on reactions of unsaturated hydrocarbons with O-atoms and methane with atoms/radicals.