Quasiclassical trajectory study of H+SiH4 reactions in full-dimensionality reveals atomic-level mechanisms

Superalkali nature of the Si9M5 (M = Li, Na, and K) Zintl clusters: a theoretical study on electronic structure and dynamic nonlinear optical properties

Zintl clusters have attracted widespread attention because of their intriguing bonding and unusual physical properties. We explore the Si9 and Si9M5 (where M = Li, Na, and K) Zintl clusters using the density functional theory combined with other methods. The exothermic nature of the Si9M5 cluster formation is disclosed, and the interactions of alkali metals with pristine Si9 are shown to be noncovalent. The reduced density gradient analysis is performed, in which increased van der Waals interactions are observed with the enlargement of the size of alkali metals. The influence of the implicit solvent model is considered, where the hyperpolarizability (βo) in the solvent is found to be about 83 times larger than that in the gas phase for Si9K5. The frequency-dependent nonlinear optical (NLO) response for the dc-Kerr effect is observed up to 1.3 × 1011 au, indicating an excellent change in refractive index by an externally applied electric field. In addition, natural bonding orbitals obtained from the second-order perturbation analysis show the charge transfer with the donor-acceptor orbitals. Electron localization function and localized orbital locator analyses are also performed to better understand the bonding electrons in designed clusters. The studied Zintl clusters demonstrate the superalkali character in addition to their remarkable optical and nonlinear optical properties.

Deciphering Dynamics of the Cl + SiH4 → H + SiH3Cl Reaction on a Machine Learning Made Globally Accurate Full-Dimensional Potential Energy Surface

Chemical reaction dynamics needs the joint effort from both experiment and theory, and theory is useful to rationalize the experimental results by offering intimate details of chemical reaction dynamics and to explore new reaction pathways. With the aid of machine learning, we develop here an accurate full-dimensional potential energy surface (PES) for the reaction between Cl + SiH₄. This PES can describe well the hydrogen abstraction channel to HCl + SiH₃. It can also give a good description for the hydrogen substitution channel to H + SiH₃Cl, which is the focus of the current study and has never been reported by theory. The dynamics of this substitution channel is revealed in detail by calculating ample quasi-classical trajectories (QCTs) on the new PES. The computed product angular distributions are in good agreement with the only crossed molecular beam experiment. Both theory and experiment suggest that this channel takes place mainly via the typical S_N2 inversion mechanism. Theory reveals that there also exists a novel torsion mechanism for the substitution channel. Two dynamic mechanisms are analyzed in detail. The present detailed theoretical dynamics study sheds insightful and novel understanding for this fundamentally important chemical reaction.

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.

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.

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.

Ab initio conical intersections for the Si(1D) + H2 reaction system: a lowest five singlet states study

Conical intersections and geometric phase effects of the Si(¹D) + H₂ system were clarified intuitively, and important features of them are revealed.

A global full-dimensional potential energy surface and quasiclassical trajectory study of the O(1D) + CH4 multichannel reaction

The QCT calculations based on an accurate global full-dimensional PES are capable of reproducing the experimental dynamic features for O(¹D) + CH₄.

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.

Imaging the O((1)D) + CD4 → OD + CD3 Reaction Dynamics: The Threshold of Abstraction Pathway

The O((1)D) + CD4 → OD + CD3 reaction was investigated using the crossed molecular beam technique with sliced velocity map imaging at four different collision energies: 1.6, 2.8, 4.6, and 6.8 kcal/mol. The vibrational ground state product CD3 was detected using a (2 + 1) resonance-enhanced multiphoton ionization (REMPI). Remarkably different features were found in the forward and backward scatterings, and gradually changed with the collision energy. These features were attributed to two distinctive reaction mechanisms-insertion and abstraction-that occur on the ground and excited state surfaces, respectively. Contributions from the two mechanisms were extracted from the experiment results, and a positive correlation was found between the abstraction proportion and the collision energy. The threshold for the abstraction pathway was determined and compared with results from calculations.