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.

Signature of shape resonances on the differential cross sections of the S(1D)+H2 reaction

Shape resonances appear when the system is trapped in an internuclear potential well after tunneling through a barrier. They manifest as peaks in the collision energy dependence of the cross section (excitation function), and in many cases, their presence can be observed experimentally. High-resolution crossed-beam experiments on the S(¹D) + H₂(j = 0) reaction in the 0.81-8.5 meV collision energy range reaction revealed non-monotonic behavior and the presence of oscillations in the reaction cross section as a function of the collision energy, as predicted by quantum mechanical (QM) calculations. In this work, we have analyzed the effect of shape resonances on the differential cross sections for this insertion reaction by performing additional QM calculations. We have found that, in some cases, the resonance gives rise to a large enhancement of extreme backward scattering for specific final states. Our results also show that, in order to yield a significant change in the state-resolved differential cross section, the resonance has to be associated with constructive interference between groups of partial waves, which requires not getting blurred by the participation of many product helicity states.

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.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Dynamics study on the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction: insertion-abstraction mechanism

Time-dependent wave packet calculations are carried out for two reaction channels of the non-adiabatic Na(3p) + HD → NaH/NaD + D/H reaction. The potential well on the excited state potential energy surface makes the reaction preferable to proceed through the insertion reaction path. The dominance of the NaD + H reaction channel and product rotational state distributions are found to be in agreement with the characteristics of typical adiabatic insertion reactions. However, significant forward scattering peaks in the differential cross sections (DCS) are found to be inconsistent with the forward-backward symmetric scattering characteristic of typical adiabatic insertion reactions, which indicate that the Na(3p) + HD reaction is dominated by a direct reaction mechanism. The comparison between adiabatic and non-adiabatic calculated DCSs reveals that the non-adiabatic couplings in the reaction could reduce the lifetime of the intermediate complex. Finally, the insertion-abstraction mechanism is put forward for the non-adiabatic Na(3p) + HD reaction.

An Empirical Dynamical Barrier for Statistical Theory of Low-Energy Reactive S(1D) + HD(j = 0), H2(j = 0) Collisions

A simple model potential is proposed to describe the dynamical barrier in the mean interaction potential at small distances between the reactants in S(¹D) + HD(¹Σ, v = 0, j = 0) reaction. The statistical theory of collision complex formation and complex decay is applied to calculate the total reaction cross sections and the cross sections for SH and SD productions in the range of low collision energies E_c = (0.4-60) meV. The results are compared with measured cross sections and results of hyperspherical close coupling calculations. As a check of consistency the same comparisons are presented for the case of S(¹D) + H₂(¹Σ, v = 0, j = 0) reaction.

A new potential energy surface for the H2S system and dynamics study on the S((1)D) + H2(X(1)Σg(+)) reaction

We constructed a new global potential energy surface (PES) for the electronic ground state ((1)A') of H2S based on 21,300 accurate ab initio energy points over a large configuration space. The ab initio energies are obtained from multireference configuration interaction calculations with a Davidson correction using basis sets of quadruple zeta quality. The neural network method is applied to fit the PES, and the root mean square error of fitting is small (1.68 meV). Time-dependent wave packet studies for the S((1)D) + H2(X(1)Σg(+)) → H((2)S) + SH(X(2)Π) reaction on the new PES are conducted to study the reaction dynamics. The calculated integral cross sections decrease with increasing collision energy and remain fairly constant within the high collision energy range. Both forward and backward scatterings can be observed as expected for a barrierless reaction with a deep well on the PES. The calculated integral cross sections and differential cross sections are in good agreement with the experimental results.

Wave packet and statistical quantum calculations for the He + NeH⁺ → HeH⁺ + Ne reaction on the ground electronic state

A real wave packet based time-dependent method and a statistical quantum method have been used to study the He + NeH(+) (v, j) reaction with the reactant in various ro-vibrational states, on a recently calculated ab initio ground state potential energy surface. Both the wave packet and statistical quantum calculations were carried out within the centrifugal sudden approximation as well as using the exact Hamiltonian. Quantum reaction probabilities exhibit dense oscillatory pattern for smaller total angular momentum values, which is a signature of resonances in a complex forming mechanism for the title reaction. Significant differences, found between exact and approximate quantum reaction cross sections, highlight the importance of inclusion of Coriolis coupling in the calculations. Statistical results are in fairly good agreement with the exact quantum results, for ground ro-vibrational states of the reactant. Vibrational excitation greatly enhances the reaction cross sections, whereas rotational excitation has relatively small effect on the reaction. The nature of the reaction cross section curves is dependent on the initial vibrational state of the reactant and is typical of a late barrier type potential energy profile.

State-resolved time-dependent wave packet and quasiclassical trajectory studies of the adiabatic reaction S(3P) + HD on the (1(3)A″) state

Time-dependent wave packet (TDWP) and quasiclassical trajectory (QCT) calculations have been carried out for the reaction S(3P) + HD(X1Σg+) at the lowest 13A″ state with both rotational and vibrational excitations of reactant HD. The calculated integral cross sections from QCT agree fairly well with the TDWP calculations. The reaction probability results from TDWP show that the reaction displays a strong tendency to the SD channel. When the reactant HD is vibrationally excited, both channels are promoted apparently. The vibration of the HD bond tends to reduce the difference of reactivity between the two channels. The detailed state-to-state differential cross sections (DCSs) are calculated. These distributions show some significant characters of the barrier-type reactions. At the same time, the scattering width of product SD has a certain relationship with its rotation excitation. For the vector properties, P(θr), P(r), and P(θr,r) distributions are calculated by QCT, and the increased collision energy weakens the rotational polarization of the SD molecule.

INVESTIGATION OF THE EXCHANGE REACTION H + H′S → HS + H′ ON THE 1A′ STATE POTENTIAL ENERGY SURFACE

The quantum time dependent wave packet (TDWP) and quasiclassical trajectory (QCT) calculations were carried out to study the exchange reaction H(2S) + H′S(2Π) → HS(2Π) + H′(2S) on the 1A′ potential energy surface (PES). The integral cross sections of the H + H′S (v = j = 0) → HS + H′ reaction calculated by the two methods were presented. The results reveal that the integral cross sections (ICS) decrease with the collision energy increasing. The result of the QCT calculations is reasonably consistent with the time-dependent wave packet. Moreover, the differential cross sections (DCS) were calculated by the QCT method at the four different collision energies, which display a forward–backward symmetry. A long-lifetime H2S intermediate complex of the exchange reaction was found according to the trajectories. In the stereodynamics investigation, the polar and dihedral angle distribution functions were calculated, which have the distinct oscillations. The oscillations could be attributed to the deep well on the 1A′ PES. However, based on the polar-angle and dihedral angle distribution functions, it could be predicted that the main product rotational angular momentum preferentially point to the positive or negative direction of y-axes.

Rate coefficients from quantum and quasi-classical cumulative reaction probabilities for the S(1D) + H2 reaction

Cumulative reaction probabilities (CRPs) at various total angular momenta have been calculated for the barrierless reaction S((1)D) + H(2) → SH + H at total energies up to 1.2 eV using three different theoretical approaches: time-independent quantum mechanics (QM), quasiclassical trajectories (QCT), and statistical quasiclassical trajectories (SQCT). The calculations have been carried out on the widely used potential energy surface (PES) by Ho et al. [J. Chem. Phys. 116, 4124 (2002)] as well as on the recent PES developed by Song et al. [J. Phys. Chem. A 113, 9213 (2009)]. The results show that the differences between these two PES are relatively minor and mostly related to the different topologies of the well. In addition, the agreement between the three theoretical methodologies is good, even for the highest total angular momenta and energies. In particular, the good accordance between the CRPs obtained with dynamical methods (QM and QCT) and the statistical model (SQCT) indicates that the reaction can be considered statistical in the whole range of energies in contrast with the findings for other prototypical barrierless reactions. In addition, total CRPs and rate coefficients in the range of 20-1000 K have been calculated using the QCT and SQCT methods and have been found somewhat smaller than the experimental total removal rates of S((1)D).

Dynamics of the S(1D2)+HD(j=0) reaction at collision energies approaching the cold regime: a stringent test for theory

We report integral cross sections for the S(1D2)+HD(j=0)→DS+H and HS+D reaction channels obtained through crossed-beam experiments reaching collision energies as low as 0.46 meV and from adiabatic time-independent quantum-mechanical calculations. While good overall agreement with experiment at energies above 10 meV is observed, neither the product channel branching ratio nor the low-energy resonancelike features in the HS+D channel can be theoretically reproduced. A nonadiabatic treatment employing highly accurate singlet and triplet potential energy surfaces is clearly needed to resolve the complex nature of the reaction dynamics.