Single-Valued DMBE Potential Energy Surface for HSO: A Distributed n-Body Polynomial Approach

Theoretical Investigation on Electronic Excited States of the HSO Radical

HSO radicals play an important role in the photochemical processes in combustion, the atmosphere, and the interstellar medium. In this work, we perform a high-level ab initio study on the electronic excited states of HSO using the internally contracted multireference configuration interaction methods including Davidson correction (icMRCI + Q) in combination with the correlation-consistent basis sets. The molecular geometries, vertical transition energies, oscillator strengths, and electronic configurations of 19 electronic states of HSO are computed. The electronic potential energy curves of HSO along the bond lengths and bond angles are presented. Based on our calculations, the interactions between the electronic states involved in the ultraviolet region and the mechanism of photodissociation are discussed, which will lay a foundation for revealing the dissociation dynamics of gas-phase HSO molecules in outer space and the earth's atmosphere.

Accurate High-Level Ab Initio-Based Global Potential Energy Surface and Quantum Dynamics Calculation for the First Excited State of CH2. J Phys Chem A 2021;125:5490-5498. [PMID: 34137628 DOI: 10.1021/acs.jpca.1c02413]

A full three-dimensional global potential energy surface (PES), covering the whole configuration space, is reported first for the title system by fitting high-level ab initio energies at the multireference configuration interaction level with the aug-cc-pV6Z basis set. In this work, the many-body expansion method is invoked to fit the innate character of the CH₂⁺(1²A″) PES. The topographical features are examined in detail based on the new global PES and in accordance with the other calculations from the ab initio energies, which show the correct behavior at the C⁺(²P) + H₂(X¹Σ_g⁺) and CH⁺(a³Π) + H(²S) dissociation limits. Using a time-dependent wave packet method, we provide insights into the dynamics behavior for reaction of C⁺(²P) + H₂(X¹Σ_g⁺) → CH⁺(a³Π) + H(²S). The integral cross sections and reaction probabilities increase monotonically in terms of the collision energy.

Accurate global potential energy surface by extrapolation to the complete basis set limit and dynamics studies for ground state of H2S. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2019;219:267-273. [PMID: 31048256 DOI: 10.1016/j.saa.2019.04.051]

The atomic ions with hydrogen reactions play important role in biological chemistry, atmospheric and combustion chemistry. In this work, a full three-dimensional global many-body expansion potential energy surface (PES) for the ground state of H₂S^-(²A') is reported for the first time, since there is almost no investigations can accurately and completely describe the PES. It is obtained using many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit, for improving the accuracy of the potential energy surface. The topographical features of the new three-dimensional global potential energy surface are examined in detail, and found to be in good agreement with theoretical and experimental results. To better test the reaction, such a surface was applied to quasi classical trajectory calculations for S^-(²P)+H₂(¹Σ_g⁺)→SH^-(¹Σ)+H(²S) reaction. The integral cross sections, differential cross sections and the rate coefficients have been computed. Results indicate that such a work can be recommended for corresponding atmospheric dynamics studies and as building blocks for constructing the many body expansion potential energy surface of larger S^-/H containing systems.

Cn (n=2-4): current status

The major aspects of the C₂, C₃ and C₄ elemental carbon clusters are surveyed. For C₂, a brief analysis of its current status is presented. Regarding C₃, the most recent results obtained in our group are reviewed with emphasis on modelling its potential energy surface which is particularly complicated due to the presence of multiple conical intersections. As for C₄, the most stable isomeric forms of both triplet and singlet spin states and their possible interconversion pathways are examined afresh by means of accurate ab initio calculations. The main strategies for modelling the ground triplet C₄ potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C₃, an approximate four-body term is calibrated from the ab initio energies. The final six-dimensional global DMBE form so obtained reproduces all known topographical aspects while providing an accurate description of the C₄ linear-rhombic isomerization pathway. It is therefore commended for both spectroscopic and reaction dynamics studies.This article is part of the theme issue 'Modern theoretical chemistry'.

Accurate global potential energy surface for the ground state of CH2+ by extrapolation to the complete basis set limit

A full three-dimensional global potential energy surface is reported for the ground state of CH₂⁺ by fitting accurate multireference configuration interaction energies calculated using aug-cc-pVQZ and aug-cc-pV5Z basis sets with extrapolation of the electron correlation energy to the complete basis set limit.

Accurate high level ab initio-based global potential energy surface and dynamics calculations for ground state of CH2(+)

A global many-body expansion potential energy surface is reported for the electronic ground state of CH2 (+) by fitting high level ab initio energies calculated at the multireference configuration interaction level with the aug-cc-pV6Z basis set. The topographical features of the new global potential energy surface are examined in detail and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. In turn, in order to validate the potential energy surface, a test theoretical study of the reaction CH(+)(X(1)Σ(+))+H((2)S)→C(+)((2)P)+H2(X(1)Σg (+)) has been carried out with the method of time dependent wavepacket on the title potential energy surface. The total integral cross sections and the rate coefficients have been calculated; the results determined that the new potential energy surface can both be recommended for dynamics studies of any type and as building blocks for constructing the potential energy surfaces of larger C(+)/H containing systems.

Accurate ab initio-based adiabatic global potential energy surface for the 2(2)A" state of NH2 by extrapolation to the complete basis set limit

A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N((2)D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N((2)D)+H2(X(1)Σg (+))(ν=0,j=0)→NH(a(1)Δ)+H((2)S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction.

A quasi-classical trajectory study of the OH + SO reaction: the role of ro-vibrational energy

The reaction cross section and the reaction probability after forming the complex as a function of vibrational energy of reactants.

A global ab initio potential energy surface for the X2A' ground state of the Si + OH → SiO + H reaction

We report the first global potential energy surface (PES) for the X(2)A' ground electronic state of the Si((3)P) + OH(X(2)Π) → SiO(X(1)Σg(+)) + H((2)S) reaction. The PES is based on a large number of ab initio energies obtained from multireference configuration interaction calculations plus Davidson correction (MRCI+Q) using basis sets of quadruple zeta quality. Corrections were applied to the ab initio energies in the reactant channel allowing a proper description of long-range interactions between Si((3)P) and OH(X(2)Π). An analytical representation of the global PES has been developed by means of the reproducing kernel Hilbert space method. The reaction is found barrierless. Two minima, corresponding to the SiOH and HSiO isomers, and six saddle points, among which the isomerization transition state, have been characterized on the PES. The vibrational spectra of the SiOH/HSiO radicals have been computed from second-order perturbation theory and quantum dynamics methods. The structural, energetic, and spectroscopic properties of the two isomers are in good agreement with experimental data and previous high quality calculations.

Accurate double many-body expansion potential energy surface by extrapolation to the complete basis set limit and dynamics calculations for ground state of NH2

An accurate single-sheeted double many-body expansion potential energy surface is reported for the title system. A switching function formalism has been used to warrant the correct behavior at the H2(X1Σg+)+N(2D) and NH (X3Σ-)+H(2S) dissociation channels involving nitrogen in the ground N(4S) and first excited N(2D) states. The topographical features of the novel global potential energy surface are examined in detail, and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. The novel surface can be using to treat well the Renner-Teller degeneracy of the 12A″ and 12A' states of NH 2. Such a work can both be recommended for dynamics studies of the N(2D)+H2 reaction and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen-containing systems. In turn, a test theoretical study of the reaction N(2D)+H2(X1Σg+)(ν=0,j=0)→NH (X3Σ-)+H(2S) has been carried out with the method of quantum wave packet on the new potential energy surface. Reaction probabilities, integral cross sections, and differential cross sections have been calculated. Threshold exists because of the energy barrier (68.5 meV) along the minimum energy path. On the curve of reaction probability for total angular momentum J = 0, there are two sharp peaks just above threshold. The value of integral cross section increases quickly from zero to maximum with the increase of collision energy, and then stays stable with small oscillations. The differential cross section result shows that the reaction is a typical forward and backward scatter in agreement with experimental measurement result.

Accurate ab initio based DMBE potential energy surface for the ground electronic state of N2H2

A global single-sheeted double many-body expansion potential energy surface is reported for the ground electronic state of N(2)H(2). Starting from an approximate cluster expansion of the molecular potential that utilizes previously reported functions of the same family for the triatomic fragments, four-body energy terms have been calibrated from extensive accurate ab initio data so as to reproduce the main features of the title system. The switching function formalism previously suggested for three-body systems [A. J. C. Varandas and L. Poveda, Theor. Chem. Acc. 116, 404 (2006)] has been generalized to approximate the true multisheeted nature of N(2)H(2) potential energy surface, thus allowing the correct behavior at the N((2)D) + NH(2)((2)A(")) and N((4)S) + NH(2)((4)A(")) dissociation limits. The resulting fully six-dimensional potential energy function reproduces the correct symmetry under permutation of identical atoms and predicts the main stationary points of the molecule in the valence and long-range regions in good agreement with available experimental and theoretical data on the diazene molecule.

On the vibrational spectra of HSO and SOH

The harmonic and fundamental vibrational frequencies, rotational constants, vibration-rotation corrections and Zero point energies of HSO and SOH and their deuterated isomers (DSO and SOD) have been determined employing the CCSD(T) methodology in conjunction with the aug-cc-pV(X+d)Z and cc-pwCVQZ basis sets. The calculated fundamental frequencies of SOH are 830, 1150 and 3577cm(-1), for the upsilon(SO), upsilon(bend) and upsilon(OH), respectively. In the case of HSO the computed fundamentals are 1002, 1077 and 2335cm(-1), for the upsilon(SO), upsilon(bend) and upsilon(SH), respectively. The values are discussed in terms of the experimental determinations available. The rotational constants reported for HSO are in reasonable agreement with experiment; the computed values are 301,271, 20,557 and 19,192MHz. In the case of SOH, for the first time, we report the rotational constants including vibration-rotation corrections, they are: 654,236, 16,621, 16,178MHz. The force fields calculated allowed as to estimate accurate ZPEs, the suggested values are 6.48 and 8.19kcal/mol for HSO and SOH, respectively. Finally, we recommend the following structural parameters for HSO r(SO)=1.4924A, r(SH)=1.3649A and angleHSO=104.76 degrees ; whereas for SOH we recommend r(SO)=1.6302A, r(OH)=0.9629A and angleSOH=107.97 degrees .

Accurate ab initio based multisheeted double many-body expansion potential energy surface for the three lowest electronic singlet states of H3+

The authors present diabatic and adiabatic potential energy surfaces for the three lowest electronic singlet states of H3+. The modeling of the surfaces is based on the multi-sheeted double many-body expansion method which consists of dressing the various matrix elements of the diatomics-in-molecules potential matrix with three-body terms. The avoided crossing between the two lowest states and the conical intersection between the second and the third state are accurately represented by construction.

Recalibrated Double Many-Body Expansion Potential Energy Surface and Dynamics Calculations for HN2

A single-sheeted double many-body expansion potential energy surface is reported for the lowest doublet state of HN2 by fitting additional multireference configuration interaction energies in the N...NH channel. A stratified analysis of the root-mean-squared error indicates an accuracy superior to that achieved for the previously reported form. Detailed dynamical tests are also performed for the N + NH reaction using both the quasi-classical trajectory method and the capture theory, and the results are compared with available empirical data. The vibrational resonances of the HN2 metastable radical are also calculated and compared with previous theoretical predictions.

Theoretical Studies on the Low-Lying Electronic States of the HSO Neutral Radical and Its Cation

Using the complete active space self-consistent field (CASSCF) method with large atomic natural orbital (ANO-L) basis set, four electronic states of the HSO neutral radical are optimized. The vertical transitions of the HSO neutral radical are investigated by using the same method under the basis set of ANO-L functions augmented with a series of adapted 1s1p1d Rydberg functions, through which eight valence states and eight Rydberg states are probed. Ionic states of the HSO neutral radical are extensively studied in both cases of the adiabatic and vertical ionization, from which the relatively complete understanding of ionization energies is given. To include further correlation effects, the second-order perturbation method (CASPT2) is implemented, and the comparison between CASSCF and CASPT2 methods is performed.

New Double Many-Body Expansion Potential Energy Surface for Ground-State HCN from a Multiproperty Fit to Accurate ab Initio Energies and Rovibrational Calculations

An accurate single-sheeted double many-body expansion potential energy surface has been obtained for the ground electronic state of the hydrogen cyanide molecule via a multiproperty fit to ab initio energies and rovibrational data. This includes 106 rovibrational levels and 2313 discrete points, which are fit with a rmsd of 4 cm(-1) and 2.42 kcal mol(-1), respectively, and seven zero first-derivatives that are reproduced at three stationary points. Since the potential also describes accurately the appropriate asymptotic limits at the various dissociation channels, it is commended both for the simulation of rovibrational spectra and reaction dynamics.

Repulsive double many-body expansion potential energy surface for the reactions N(4S)+ H2⇌ NH(X3Σ–)+ H from accurate ab initio calculations

A single-sheeted DMBE potential energy surface is reported for the reactions N(4S)+H2<-->NH(X3Sigma-)+H based on a fit to accurate multireference configuration interaction energies. These have been calculated using the aug-cc-pVQZ basis set of Dunning and the full valence complete active space wave function as reference, being semi-empirically corrected by scaling the two-body and three-body dynamical correlation energies. The topographical features of the novel global potential energy surface are examined in detail, including a conical intersection involving the two first 4A'' potential energy surfaces which has been transformed into an avoided crossing in the present single-sheeted representation.