Dynamics of the C(1D)+H2reaction: A comparison of crossed molecular beam experiments with quantum mechanical and quasiclassical trajectory calculations on the first two singlet (11A′ and 11A″) potential energy surfaces

Gas-phase formation of silicon monoxide via non-adiabatic reaction dynamics and its role as a building block of interstellar silicates

Silicon monoxide (SiO) is classified as a key precursor and fundamental molecular building block to interstellar silicate nanoparticles, which play an essential role in the synthesis of molecular building blocks connected to the Origins of Life. In the cold interstellar medium, silicon monoxide is of critical importance in initiating a series of elementary chemical reactions leading to larger silicon oxides and eventually to silicates. To date, the fundamental formation mechanisms and chemical dynamics leading to gas phase silicon monoxide have remained largely elusive. Here, through a concerted effort between crossed molecular beam experiments and electronic structure calculations, it is revealed that instead of forming highly-stable silicon dioxide (SiO₂), silicon monoxide can be formed via a barrierless, exoergic, single-collision event between ground state molecular oxygen and atomic silicon involving non-adiabatic reaction dynamics through various intersystem crossings. Our research affords persuasive evidence for a likely source of highly rovibrationally excited silicon monoxide in cold molecular clouds thus initiating the complex chain of exoergic reactions leading ultimately to a population of silicates at low temperatures in our Galaxy.

Gas-Phase Preparation of Subvalent Germanium Monoxide (GeO, X1+) via Non-Adiabatic Reaction Dynamics in the Exit Channel

The subvalent germanium monoxide (GeO, X¹Σ⁺) molecule has been prepared via the elementary reaction of atomic germanium (Ge, ³P_j) and molecular oxygen (O₂, X³Σ_g^-) with each reactant in its electronic ground state by means of single-collision conditions. The merging of electronic structure calculations with crossed beam experiments suggests that the formation of germanium monoxide (GeO, X¹Σ⁺) commences on the singlet surface through unimolecular decomposition of a linear singlet collision complex (GeOO, i1, C_∞v, ¹Σ⁺) via intersystem crossing (ISC) yielding nearly exclusively germanium monoxide (GeO, X¹Σ⁺) along with atomic oxygen in its electronic ground state [p1, O(³P)]. These results provide a sophisticated reaction mechanism of the germanium-oxygen system and demonstrate the efficient "heavy atom effect" of germanium in ISC yielding (nearly) exclusive singlet germanium monoxide and triplet atomic oxygen compared to similar systems (carbon dioxide and dinitrogen monoxide), in which non-adiabatic reaction dynamics represent only minor channels.

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.

Low-Temperature Rate Constants and Product-Branching Ratios for the C(1D) + H2O Reaction

The gas-phase reaction between atomic carbon in its first electronically excited ¹D state and water has been studied over the 50-296 K temperature range using a supersonic flow apparatus. C(¹D) atoms were produced by pulsed ultraviolet multiphoton dissociation of carbon tetrabromide; a process that also generates ground-state atomic carbon C(³P). The reaction was followed by detecting product H-atoms by pulsed vacuum ultraviolet laser-induced fluorescence. Two types of experiment were performed. First, temperature-dependent rate constants were derived by recording H-atom formation curves at various gas-phase water concentrations at each temperature. Secondly, temperature-dependent H-atom yields were extracted by comparing the H-atom fluorescence intensities generated by the target C(¹D) + H₂O reaction with those of a reference reaction. The second-order rate constants are large and increase to low temperature, whereas the measured H-atom yields are close to the theoretical maximum value of 2 above 100 K. At 50 K, neither rate constants nor H-atom yields could be derived because of H-atom formation by quantum tunneling in the activated C(³P) + H₂O reaction. The present results are discussed in the context of earlier work on the C(¹D)/C(³P) + H₂O reactions.

Conical intersection-regulated intermediates in bimolecular reactions: Insights from C(1D) + HD dynamics

The importance of conical intersections (CIs) in electronically nonadiabatic processes is well known, but their influence on adiabatic dynamics has been underestimated. Here, through combined experimental and theoretical studies, we show that CIs induce a barrier and regulate conversion from a precursor metastable intermediate (CI-R) to a deep well one. This results in bond-selective activation, influencing the adiabatic dynamics markedly in the C(1D) + HD reaction. Theory is validated by experiment; quantum dynamics calculations on highly accurate ab initio potential energy surfaces yield rate coefficients and product branching ratios in excellent agreement with the experiment. Quasi-classical trajectory calculations reveal that the CI-R intermediate leads to unusual reaction mechanisms (designated as C─H activation complex conversion and cyclic complex), which are responsible for large branching ratios. We also reveal that CI-R intermediates exist in other reactive systems, and the dynamical effects uncovered here may have general significance.

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.

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

Rate constants derived from ring polymer molecular dynamics calculations confirm the validity of this method for studying low-temperature complex-forming reactions

A Combined Experimental and Theoretical Study on the Formation of the 2-Methyl-1-silacycloprop-2-enylidene Molecule via the Crossed Beam Reactions of the Silylidyne Radical (SiH; X(2)Π) with Methylacetylene (CH3CCH; X(1)A1) and D4-Methylacetylene (CD3CCD; X(1)A1)

The bimolecular gas-phase reactions of the ground-state silylidyne radical (SiH; X(2)Π) with methylacetylene (CH3CCH; X(1)A1) and D4-methylacetylene (CD3CCD; X(1)A1) were explored at collision energies of 30 kJ mol(-1) under single-collision conditions exploiting the crossed molecular beam technique and complemented by electronic structure calculations. These studies reveal that the reactions follow indirect scattering dynamics, have no entrance barriers, and are initiated by the addition of the silylidyne radical to the carbon-carbon triple bond of the methylacetylene molecule either to one carbon atom (C1; [i1]/[i2]) or to both carbon atoms concurrently (C1-C2; [i3]). The collision complexes [i1]/[i2] eventually isomerize via ring-closure to the c-SiC3H5 doublet radical intermediate [i3], which is identified as the decomposing reaction intermediate. The hydrogen atom is emitted almost perpendicularly to the rotational plane of the fragmenting complex resulting in a sideways scattering dynamics with the reaction being overall exoergic by -12 ± 11 kJ mol(-1) (experimental) and -1 ± 3 kJ mol(-1) (computational) to form the cyclic 2-methyl-1-silacycloprop-2-enylidene molecule (c-SiC3H4; p1). In line with computational data, experiments of silylidyne with D4-methylacetylene (CD3CCD; X(1)A1) depict that the hydrogen is emitted solely from the silylidyne moiety but not from methylacetylene. The dynamics are compared to those of the related D1-silylidyne (SiD; X(2)Π)-acetylene (HCCH; X(1)Σg(+)) reaction studied previously in our group, and from there, we discovered that the methyl group acts primarily as a spectator in the title reaction. The formation of 2-methyl-1-silacycloprop-2-enylidene under single-collision conditions via a bimolecular gas-phase reaction augments our knowledge of the hitherto poorly understood silylidyne (SiH; X(2)Π) radical reactions with small hydrocarbon molecules leading to the synthesis of organosilicon molecules in cold molecular clouds and in carbon-rich circumstellar envelopes.

Global analytical ab initio ground-state potential energy surface for the C((1)D)+H2 reactive system

A new global ab initio potential energy surface (called ZMB-a) for the 1(1)A' state of the C((1)D)+H2 reactive system has been constructed. This is based upon ab initio calculations using the internally contracted multireference configuration interaction approach with the aug-cc-pVQZ basis set, performed at about 6300 symmetry unique geometries. Accurate analytical fits are generated using many-body expansions with the permutationally invariant polynomials, except that the fit of the deep well region is taken from our previous fit. The ZMB-a surface is unique in the accurate description of the regions around conical intersections (CIs) and of van der Waals (vdW) interactions. The CIs between the 1(1)A' and 2(1)A' states cause two kinds of barriers on the ZMB-a surface: one is in the linear H-CH dissociation direction with a barrier height of 9.07 kcal/mol, which is much higher than those on the surfaces reported before; the other is in the C((1)D) collinearly attacking H2 direction with a barrier height of 12.39 kcal/mol. The ZMB-a surface basically reproduces our ab initio calculations in the vdW interaction regions, and supports a linear C-HH vdW complex in the entrance channel, and two vdW complexes in the exit channel, at linear CH-H and HC-H geometries, respectively.