Bonding and correlation analysis of various Si2CO isomers on the potential energy surface

At various levels of theory, singlet and triplet potential energy surfaces (PESs) of Si₂CO, which has been studied using matrix isolation infrared spectroscopy, are investigated in detail. A total of 30 isomers and 38 interconversion transition states are obtained at the B3LYP/6‐311G(d) level. At the higher CCSD(T)/6‐311+G(2d)//QCISD/6‐311G(2d)+ZPVE level, the global minimum ¹1 (0.0 kcal/mol) corresponds to a three‐membered ring singlet O‐cCSiSi (¹A′). On the singlet PES, the species ¹2 (0.2 kcal/mol) is a bent SiCSiO structure with a ¹A′ electronic state, followed by a three‐membered ring isomer Si‐cCSiO (¹A′) ¹3 (23.1 kcal/mol) and a linear SiCOSi ¹4 (¹Σ⁺) (38.6 kcal/mol). The isomers ¹1, ¹2, ¹3 , and ¹4 possess not only high thermodynamic stabilities, but also high kinetic stabilities. On the triplet PES, two isomers ³1 (³B₂) (18.8 kcal/mol) and ³7 (³A″) (23.3 kcal/mol) also have high thermodynamic and kinetic stabilities. The bonding natures of the relevant species are analyzed. The similarities and differences between C₃O, C₃S, SiC₂O, and SiC₂S are discussed. The present results are also expected to be useful for understanding the initial growing step of the CO‐doped Si vaporization processes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Theoretical investigation of the cyclic GaO2 and GaS2 molecules at DFT and correlated wave function levels

The geometries and the bonding properties have been predicted for cyclic GaO₂ and GaS₂ species at density functional theory (DFT), MPn (n=2,3,4 with different substitutions), QCISD(T), and CCSD(T) all‐electron correlation levels with 6‐311+G* basis set. The geometrical optimizations and the harmonic vibrational frequency analysis are performed using DFT and second‐order Møller–Plesset (MP2) methods. The relevant energy quantities are also calibrated at the high‐order electron correlation levels [MP3, MP4, quadratic configuration interaction (QCI), and coupled cluster (CC)]. Each species possesses a ²A₂ ground state with a higher energy level ²A₁ state. The corresponding state–state separations are about 32 kcal/mol for GaO₂ species and about 20 kcal/mol for GaS₂ species at the QCISD(T)/6‐311+G* level. The QCISD(T) and CCSD(T) calculations yield dissociation energies of 42.0 and 59.0 kcal/mol for two species, respectively, and other methods yield dissociation energies within ∼5 kcal/mol. Result analysis has indicated that the cyclic GaO₂ should be classified as superoxide and the GaS₂ species should be classified as supersulfide in their ground state, and those in the excited state (²A₁) should not be. However, the cyclic GaS₂ (²A₂) is less ionic than the GaO₂ (²A₂) and they are far less ionic than NaO₂. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 222–231, 2001     

The nitric oxide adsorption on gold neutral,cation, and anion atoms: A comparative ab initio MRCI—MRPT2 studies

The Au^ν + NO(²Π) → AuNO^ν reaction, for ν = ?1, 0, +1, anion, neutral, and cation are calculated and predicted at multireference configuration interaction (MRCI) and multireference second order perturbation (MRPT2) levels of theory, in this way the main parameters: reaction path surfaces, total and adsorption energies, optimized geometries, and Mulliken charges distribution are presented and compared. The AuNO (X ¹A′) complex is created spontaneously with ?11.32 and ?13.14 kcal/mol adsorption energies with the MRCI and MRPT2 approaches, respectively. The AuNO bonding in the neutral gold nitrosyl complex has a covalent character and the nitric oxide (NO) molecule is not dissociated. The others excited states (a ³A″, b ³A′, and A ¹A″) do not present bonding. The gold nitrosyl cationic (X ²A′, A ²A″, and a ⁴A″) and anionic (X ²A″ and a ⁴A″) are bonding and present a dative covalent bond. The Mulliken analysis done for ionic species show that the binding is done through soft electrostatic interactions, due to that there is some charge transfer, delocalized onto the NO molecule for the AuNO^± ionic species whereas the AuNO (X ¹A′) neutral complex presents a little charge transfer. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

High-level ab initio studies of unimolecular dissociation of the ground-state N3 radical

A comprehensive study of the unimolecular dissociation of the N(3) radical on the ground doublet and excited quartet potential energy surfaces has been carried out with multireference single and double excitation configuration interaction and second-order multireference perturbation methods. Two forms of the N(3) radical have been located in the linear and cyclic region of the lowest doublet potential energy surface with an isomerization barrier of 62.2 kcal/mol above the linear N(3). Three equivalent C(2v) minima of cyclic N(3) are connected by low barrier, meaning the molecule is free to undergo pseudorotation. The cyclic N(3) is metastable with respect to ground state products, N((4)S)+N(2), and dissociation must occur via intersystem crossing to a quartet potential energy surface. Minima on the seams of crossing between the doublet and quartet potential surfaces are found to lie substantially higher in energy than the cyclic N(3) minima. This strongly suggests that cyclic N(3) possesses a long collision-free lifetime even if formed with substantial internal excitation.     

A theoretical comparison of the lowest-lying singlet and triplet states of H2CBe and of HCBeH

The low-lying singlet and triplet states of H₂CBe and HCBeH are examined using ab inito molecular orbital theory. In agreement with earlier results, the lowest-lying structure of H₂CBe has C_2v symmetry and is a triplet with one π electron (³ B₁). The results presented here suggest that the lowest-energy singlet structure is the (¹B₁) open-shell singlet, also with C_2v symmetry, at least 2.5 kcal/mol higher in energy. The singlet C_2v structure with two π electrons (¹A₁) is 15.9 kcal/mol higher than ³B₁. All of these structures are bound with respect to the ground state of methylene and the beryllium atom. In HCBeH, linear equilibrium geometries are found for the triplet (³Σ) and singlet (¹Δ) states. The triplet is more stable than the singlet (¹Δ) by 35.4 kcal/mol, and is only 2.9 kcal/mol higher in energy than triplet H₂ CBe. Since the transition structure connecting these two triplet molecules is found to be 50.2 kcal/mol higher in energy than H₂ CBe, both triplet equilibrium species might exist independently. The harmonic vibrational frequencies of all structures are also reported.     

A Valence Bond Model for Electron‐Rich Hypervalent Species: Application to SFn (n=1, 2, 4), PF5, and ClF3

Some typical hypervalent molecules, SF₄, PF₅, and ClF₃, as well as precursors SF (⁴Σ^? state) and SF₂ (³B₁ state), are studied by means of the breathing‐orbital valence bond (BOVB) method, chosen for its capability of combining compactness with accuracy of energetics. A unique feature of this study is that for the first time, the method used to gain insight into the bonding modes is the same as that used to calculate the bonding energies, so as to guarantee that the qualitative picture obtained captures the essential physics of the bonding system. The ⁴Σ^? state of SF is shown to be bonded by a three‐electron σ bond assisted by strong π back‐donation of dynamic nature. The linear ³B₁ state of SF₂, as well as the ground states of SF₄, PF₅ and ClF₃, are described in terms of four VB structures that all have significant weights in the range 0.17–0.31, with exceptionally large resonance energies arising from their mixing. It is concluded that the bonding mode of these hypervalent species and isoelectronic ones complies with Coulson’s version of the Rundle–Pimentel model, but assisted by charge‐shift bonding. The conditions for hypervalence to occur are stated.     

Ab initio spectroscopy and thermochemistry of the BN molecule

The lowest¹Σ⁺ and³Π states of the BN molecule have been studied using the quadratic configuration interaction method and (spdf) basis sets. The lowest¹Σ⁺ and³Π states lie extremely closely (T _e≈100 cm^?1) together; it is not clear which is the ground state. The very small separation should form a useful benchmark for basis sets and electron correlation methods. The dissociation energyD ₀ is computed to be 103.9±2 kcal/mol. A self-consistent set of spectroscopic constants is derived from a combination of ab initio and experimental data. JANAF-style thermodynamic functions in the range 100–6000 K, including anharmonic, rovibrational coupling, centrifugal stretching, and spin-orbit coupling effects are computed using direct numerical summation over the 25 lowest electronic states. A modified procedure for the latter is outlined that reduces computer time by one or two orders of magnitude without compromise in accuracy.     

Density functional study of guanine and uracil quartets and of guanine quartet/metal ion complexes

The structures and interaction energies of guanine and uracil quartets have been determined by B3LYP hybrid density‐functional calculations. The total interaction energy ΔE^T of the C_4h‐symmetric guanine quartet consisting of Hoogsteen‐type base pairs with two hydrogen bonds between two neighbor bases is −66.07 kcal/mol at the highest level. The uracil quartet with C6 H6O4 interactions between the individual bases has only a small interaction energy of −20.92 kcal mol⁻¹, and the interaction energy of −24.63 kcal/mol for the alternative structure with N3 H3O4 hydrogen bonds is only slightly more negative. Cooperative effects contribute between 10 and 25% to all interaction energies. Complexes of metal ions with G‐quartets can be classified into different structure types. The one with Ca²⁺ in the central cavity adopts a C_4h‐symmetric structure with coplanar bases, whereas the energies of the planar and nonplanar Na⁺ complexes are almost identical. The small ions Li⁺, Be²⁺, Cu⁺, and Zn²⁺ prefer a nonplanar S₄‐symmetric structure. The lack of coplanarity prevents probably a stacking of these base quartets. The central cavity is too small for K⁺ ions and, therefore, this ion favors in contrast to all other investigated ions a C₄‐symmetric complex, which is 4.73 kcal/mol more stable than the C_4h‐symmetric one. The distance 1.665 Å between K⁺ and the root‐mean‐square plane of the guanine bases is approximately half of the distance between two stacked G‐quartets. The total interaction energy of alkaline earth ion complexes exceeds those with alkali ions. Within both groups of ions the interaction energy decreases with an increasing row position in the periodic table. The B3LYP and BLYP methods lead to similar structures and energies. Both methods are suitable for hydrogen‐bonded biological systems. Compared with the before‐mentioned methods, the HCTH functional leads to longer hydrogen bonds and different relative energies for two U‐quartets. Finally, we calculated also structures and relative energies with the MMFF94 forcefield. Contrary to all DFT methods, MMFF94 predicts bifurcated C HO contacts in the uracil quartet. In the G‐quartet, the MMFF94 hydrogen bond distances N2 H22N7 are shorter than the DFT distances, whereas the N1 H1O6 distances are longer. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 109–124, 2001     

Electronic transition moment of the AΠr-XΣ system of TiN

The A²Π_r-X²Σ⁺ transition of TiN was observed by the dispersed laser induced fluorescence (DLIF) spectroscopy. The relative intensities of the DLIF spectra were analyzed to determine the dependence of the electronic transition moment, R_e(r), on the internuclear distance, r, as R_e(r)∝{1−0.281(26)r} (1.380 Å≤r≤1.823 Å). This r-dependence was analyzed simultaneously with the reported values of the spin-orbit constants for A²Π_r and the hyperfine-coupling constants for X²Σ⁺ to evaluate the ionic character of the TiN bond, the 4s atomic character in the 9σ orbital of X²Σ⁺, and the 4p atomic character in the 4π orbital of A²Π_r. These characters were confirmed to be in accordance with the reported theoretical prediction. A strong r-dependence was indicated for the 3d-4p mixing in the A²Π_r state due to the configuration mixing of the Ti(3d⁴) and Ti(3d³4p) states at a large internuclear distance.     

Singlet-triplet separation in borylene. comparison with the methylidyne cation,methylene and the amino cation

Borylene (BH), which is known to have a ¹Σ⁺ ground state, is predicted to have a low-lying ³Π triplet state with an energy separation of 31.9 kcal/mol. The ab initio molecular-orbital method used is shown to give results which are consistent with theoretical and experimental data on CH⁺, CH₂ and NH⁺₂.     

An analysis of the hydrazine—fluorine flame. The A2A′-X2A]″ emission spectrum of HNF and its relation to other HAB co

The chemiluminescent reaction N₂H₄ + F₂ is shown to be dominated by visible emission from NH* (A³Π-X³Σ^?), HNF^* (A²A′-X²A″) and vibrationally excited HF2(X¹ Σ⁺). Possible mechanisms for formation of these species are discussed with reference to previous results from similar flame reactions. The A–X emission spectrum of HNF contains the heretofore unreported progression (0, 0, 0) → (0, υ₂^″, 0) to vibrational levels υ₂^″ = 0–3 of the X²A″ state. Each vibrational transition shows well resolved K-type subbands characteristic of a near symmetric rotor. A partial rotational analysis gives band origins and effective rotational constants for the first four vibrational levels of the ground state. The measured ground state vibrational constants are ω₀ = 1441.1 ± 1.0 cm^?1 with ω₀χ₀ = 10.9 ± 0.5 cm^?1. The substantial differences in the A–X emission spectra of HNF and its isoelectronic analog HOO are discussed as they relate to bonding trends in HAB molecules. The concepts presented are extended to other HAB systems.     

Importance of Pd and Pt excited states in N2O capture and activation: A comparative study with Rh and Au atoms

Nitrous oxide (N₂O) is an intermediate compound formed during catalysis occurring in automobile exhaust pipes. In this work, the N₂O capture and activation by Pt and Pd atoms in the ground and excited states of many multiplicities are studied. Pt and Pd + N₂O reactions are studied at multireference second‐order perturbation level of theory using C_s symmetry. The PtN₂O (¹A′, ⁵A′, and ⁵A″) species are spontaneously created from excited states. Only the ⁵A′ and ⁵A″ states exhibit N₂O activation reaction paths when N₂O approaches Pt end‐on by the N or O atoms side or side‐on yielding NO or N₂ as products, respectively. Pt⁺ cations ground and excited states, capture N₂O, although only Pt⁺ (⁶A′ and ⁶A″) states show N₂O activation yielding O and N₂ as products. In the Pd atom case, PdN₂O (¹A′ and ⁵A″) species are also spontaneously created from excited states. The ⁵A″ state exhibits N₂O activation yielding N₂ + O as products. Pd⁺ cations in both ground and excited states capture N₂O; however, only the [PdN₂O]⁺ (⁴A′, ⁴A″, ⁶A′, and ⁶A″) states in side‐on approaches and (⁶A′) in end‐on approach activate the N₂O and yield the N₂ bounded to the metal and O as product. The results obtained in this work are discussed and compared with previous calculations of Rh and Au atoms. The reaction paths show a metal–gas dative covalent bonding character. Löwdin charge population analyses for Pt and Pd active states show a binding done through charge donation and retrodonation between the metals and N₂O. © 2013 Wiley Periodicals, Inc.     

A Valence Bond Study of the Low‐Lying States of the NF Molecule

The electronic structures of the three lowest‐lying states of NF are investigated by means of modern valence bond (VB) methods such as the VB self‐consistent field (VBSCF), breathing orbital VB (BOVB), and VB configuration interaction (VBCI) methods. The wave functions for the three states are expressed in terms of 9–12 VB structures, which can be further condensed into three or four classical Lewis structures, whose weights are quantitatively estimated. Despite the compactness of the wave functions, the BOVB and VBCI methods reproduce the spectroscopic properties and dipole moments of the three states well, in good agreement with previous computational studies and experimental values. By analogy to the isoelectronic O₂ molecule, the ground state ³Σ^? possesses both a σ bond and 3‐electron π bonds. However, here the polar σ bond contributes the most to the overall bonding. It is augmented by a fractional (19 %) contribution of three‐electron π bonding that arises from π charge transfer from fluorine to nitrogen. In the singlet ¹Δ and ¹Σ⁺ excited states the π‐bonding component is classically covalent, and it contributes 28 % and 37 % to the overall bonding picture for the two states, respectively. The resonance energies are calculated and reveal that π bonding contributes at least 24, 35 and 42 kcal mol^?1 to the total bonding energies of the ³Σ^?, ¹Δ and ¹Σ⁺ states, respectively. Some unusual properties of the NF molecule, like the equilibrium distance shortening and bonding energy increasing upon excitation, the counterintuitive values of the dipole moments and the reversal of the dipole moments as the bond is stretched, are interpreted in the light of the simple valence bond picture. The overall polarity of the molecule is very small in the ground state, and is opposite to the relative electronegativity of N vs F in the singlet excited states. The values of the dipole moments in the three states are quantitatively accounted for by the calculated weights of the VB structures.     

Insight into spin–orbital interaction using MCSCF method: A special analysis of the 1Σg+ electronic state in C2 and the linear polyacetylenic C4 and C6

The symmetry-broken wave function can transform the ¹Σ_g⁺ state of C₂ from the classic double bonding to the quadruple bonding, where the transformed wave functions of ϕ _L and ϕ _R are singly occupied by two opposite-spinning electrons. In this article, the effective bond order (EBO) contribution of the fourth bond in C₂ is assessed through the overlap integral between ϕ _L and ϕ _R , namely the value (0.60) is the EBO contribution of the fourth bond in the transformed scheme. Hence, the new EBO is 3.36, which is more equitable than the original EBO (2.15) in the traditional scheme. In addition, the singlet diradical character of the linear polyacetylenic C₄ and C₆ in the ¹Σ_g⁺ state is addressed for the first time. No spin-polarized bonding exists in other linear C_2n clusters, because the ionic interaction in the polyacetylenic ¹Σ_g⁺ state of C₄ is negligible. Moreover, the coupling energy between α and β single electrons in C₄ is only 4.0 kcal mol⁻¹ based on the electron spin-flip energy. © 2019 Wiley Periodicals, Inc.     

Dehydrophenylnitrenes: quartet versus doublet states

The geometries and energies of 4-, 3-, and 2-dehydrophenylnitrenes (3, 4, and 5) are investigated using complete active space self-consistent field (CASSCF), multiconfiguration quasi-degenerate second-order perturbation (MCQDPT), and internally contracted multiconfiguration-reference configuration interaction (MRCI) theories in conjunction with a correlation consistent triple-zeta basis set. 4-Dehydrophenylnitrene 3 has a quartet ground state ((4)A(2)). The adiabatic excitation energies to the (2)A(2), (2)B(2), (2)A(1), and (2)B(1) states are 5, 21, 34, and 62 kcal mol(-1), respectively. The (2)B(2) state has pronounced closed-shell carbene/iminyl radical character, while the lowest-energy (2)B(1) state is a combination of a planar allene and a 2-iminylpropa-1,3-diyl. The MCQDPT treatment overestimates the excitation energy to (2)B(2) significantly as compared to CASSCF and MRCI+Q. Among quartet states, (4)A(2)-3 is the most stable one, while those of 4 and 5 (both (4)A') are 3 and 1 kcal mol(-1) higher in energy. 5 also has a quartet ground state and a (2)A' ' state 7 kcal mol(-1) higher in energy. On the other hand, the doublet-quartet energy splitting is -6 kcal mol(-1) for 4 in favor of the doublet state ((2)A'). Hence, (2)A'-4 is the most stable dehydrophenylnitrene, 3.5 kcal mol(-1) below (4)A(2) of 3. The geometry of (2)A'-4 shows the characteristic features of through-bond interaction between the in-plane molecular orbitals at N and at C3. The (2)A' state of 4 resembles the (2)A(1) state of 3 and lies 32 kcal mol(-1) above (4)A'-4. The lowest-energy (2)A' state of 5, on the other hand, resembles the (2)B(2) state of 3 and lies 22 kcal mol(-1) above (4)A'-5.     

Ab Initio Study on the Potential Energy Surfaces of B4 Cluster

The singlet and triplet potential energy surfaces (PES) for the isomerization and dissociation reactions of B₄ isomers have been investigated using ab initio methods. Ten B4 isomers have been identified and of these 10 species, 4 have not been reported previously. The singlet rhombic structure ¹1 is found to be the most stable on the B₄ surface, in agreement with the results of previous reports. Several isomerization and dissociation pathways have been found. On the singlet PES, the linear ¹3b can rearrange to rhombus ¹1 directly, while ¹3c rearranges to ¹1 through two‐step reactions involving a cyclic intermediate. On the triplet PES, the capped triangle structure ³2 undergoes ring opening to the linear isomer ³3b with a barrier of 34.8 kcal/mol and 44.9 kcal/mol, and the latter undergoes ring closure to the square structure ³1 with a barrier of 30.4 kcal/mol and 33.0 kcal/mol at the MP4/6–311+G(3df)//MP2/6–311G(d) and CCSD/aug‐cc‐pVTZ//MP2/6–311G(d) levels of theory, respectively. The direct decomposition of singlet B₄ yielding to B₃+B is shown to have a large endothermicity of 87.3 kcal/mol (CCSD), and that producing 2B₂ to have activation energy of 133.4 kcal/mol (CCSD).     

Configuration-Interaction study of lower excited states of O2: Valence and rydberg characters of the two lowest 3Σu − states

Configuration-interaction calculations, with an extended basis, are carried out on the ground and lower excited states of O₂ and O₂⁺ at and near the equilibrium internuclear distance (R = 2.3 a.u.) of the ground state of O₂. Particular attention has been paid to the two lowest ³Σ_u^? states, and the mixing of the valence and Rydberg characters in these states are studied. The lowest ³Σ_u^? state is a Rydberg-type state for R < 2.3 a.u., but becomes valence-type for R ? 2.3 a.u. The second ³Σ_u^? state, which is 1.6 eV above the lowest ³Σ_u^? at R = 2.3 a.u., changes its character from Rydberg to valence, valence to Rydberg, and then to valence again when R increases from 1.9 to 3.1 a.u. Satisfactory agreement between the calculated and experimental vertical excitation energies is obtained.     

Molecular beam study of the chemiluminescent reaction C + SO2 → SO(A 3Π) + CO

The reaction C(³P) + SO₂ → SO(A ³Π) + CO has been identified and studied in a beam-cell configuration. Chemiluminescence (240–360 nm) indicates that the SO(A ³Π) molecules are formed in the lowest vibrational levels (v = 0–6). The predominance of SO(A ³Π) production over SO(B ³Σ^?) formation observed is explained by adiabatic correlation arguments.     

Theoretical study of the mechanism for C-H bond activation in spin-forbidden reaction between Ti<Superscript>+</Superscript> and C<Subscript>2</Subscript>H<Subscript>4</Subscript>

The mechanism of the spin-forbidden reaction Ti⁺(⁴F, 3d²4s¹) + C₂H₄→TiC₂H₂ ⁺ (²A₂) + H₂ on both doublet and quartet potential energy surfaces has been investigated at the B3LYP level of theory. Crossing points between the potential energy surfaces and the possible spin inversion process are discussed by means of spin-orbit coupling (SOC) calculations. The strength of the SOC between the low-lying quartet state and the doublet state is 59.3 cm⁻¹ in the intermediate complex IM1-⁴B₂. Thus, the changes of its spin multiplicity may occur from the quartet to the doublet surface to form IM1-²A₁, leading to a sig-nificant decrease in the barrier height on the quartet PES. After the insertion intermediate IM2, two distinct reaction paths on the doublet PES have been found, i.e., a stepwise path and a concerted path. The latter is found to be the lowest energy path on the doublet PES to exothermic TiC₂H₂ ⁺(²A₂) + H₂ products, with the active barrier of 4.52 kcal/mol. In other words, this reaction proceeds in the following way: Ti⁺+C₂H₄→⁴IC→IM1-⁴B₂→^4,2ISC→IM1-²A₁→[²TS_ins]→IM2→[²TS_MCTS]→IM5→TiC₂H₂ ⁺(²A₂)+H₂. Supported by ‘Qinglan’ Talent Engineering Funds by Tianshui Normal University.