Theoretical studies on the structures and isomerization of the LiSiF3 system

Various possible isomers of LiSiF₃ system and isomerization between them have been studied at G2(MP2) level usingab initio calculations. The relative energies of four minimum points on the potential energy surface are-128.6,-194.3,-12.7 and-122.8 kJ/mol (taking the sum of the energies of LiF and SiF₂ as zero). The structural energy of the four-membered ring that contains three F-Si-F-Li four-membered rings with C_3v symmetry is the lowest. The highest potential barrier for the isomerization of the remaining three- or four-membered structure is 12.5 kJ/ mol. Project supported by the National Natural Science Foundation of China (Grant No. 29673026).     

Theoretical study on the mechanism of cycloaddition reaction between dimethyl germylidene and formaldehyde

The cycloaddition mechanism of the reaction between singlet dimethyl germylidene and formaldehyde has been investigated with MP2/6-31G* method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. The energies of the different conformations are calculated with CCSD (T)//MP2/6-31G* method. From the potential energy profile, we predict that the cycloaddition reaction between singlet dimethyl germylidene and formaldehyde has two dominant reaction pathways. First dominant reaction pathway consists of three steps: (1) the two reactants (R1, R2) firstly form an intermediate INT1a through a barrier-free exothermic reaction of 43.0 kJ/mol; (2) INT1a then isomerizes to a four-membered ring compound P1 via a transition state TS1a with an energy barrier of 24.5 kJ/mol; (3) P1 further reacts with formaldehyde(R2) to form a germanic heterocyclic compound INT3, which is also a barrier-free exothermic reaction of 52.7 kJ/mol; Second dominant reaction pathway is as following: (1) the two reactants (R1, R2) firstly form a planar four-membered ring intermediate INT1b through a barrier-free exothermic reaction of 50.8 kJ/mol; (2) INT1b then isomerizes to a twist four-membered ring intermediate INT1.1b via a transition state TS1b with an energy barrier of 4.3 kJ/mol; (3) INT1.1b further reacts with formaldehyde(R2) to form an intermediate INT4, which is also a barrier-free exothermic reaction of 46.9 kJ/mol; (4) INT4 isomerizes to a germanic bis-heterocyclic product P4 via a transition state TS4 with an energy barrier of 54.1 kJ/mol.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

Theoretical study of the mechanism of cycloaddition reaction between dichloro-germylidene and acetaldehyde

The mechanism of the cycloadditional reaction between singlet dichloro-germylidene(R1) and (acetaldehyde(R2) has been investigated with MP2/6-31G* method, including geometry optimization, vibrational analysis and energies for the involved stationary points on the potential energy surface. From the potential energy profile, we predict that the cycloaddition reaction between singlet dichloro-germylidene and acetaldehyde has two competitive dominant reaction pathways. Going with the formation of two side products (INT3 and INT4), simultaneously. The two competitive reactions both consist of two steps: (1) two reactants firstly form a three-membered ring intermediate (INT1) and a twisted four-membered ring intermediate (INT2), respectively, both of which are barrier-free exothermic reactions of 44.5 and 63.0 kJ/mol; (2) then INT1 and INT2 further isomerize to a four-membered ring product (P1) and a chlorine-transfer product (P2) via transitions (TS1 and TS2), respectively, with the barriers of 9.3 and 1.0 kJ/mol; simultaneously, P1 and INT2 react further with acetaldehyde(R2) to give two side products (INT3 and INT4), respectively, which are also barrier-free exothermic reaction of 65.4 and 102.7 kJ/mol.     

Geometries, thermodynamic properties and reactions of methylzinc alkoxide clusters studied by density functional theory calculations

Methylzinc alkoxide complexes are precursors for the preparation of nanosized zinc oxide particles, which in turn are catalysts or reagents in important industrial processes such as methanol synthesis and rubber vulcanization. We report for the first time the structures, energies, atomic charges, dipole moments, and vibrational spectra of more than 20 species of the type [(MeZnOR')n] with R' = H, Me, tBu and n = 1-6, calculated by density functional theory methods, mostly at the B3LYP/6-31+G* level of theory. For R' = Me, the global minimum structure of the tetramer (n = 4) is a highly symmetrical heterocubane but a ladder-type isomer is by only 70.9 kJ mol(-1) less stable. The corresponding trimer is most stable as a rooflike structure; a planar six-membered ring of relative energy 13.5 kJ mol(-1) corresponds to a saddle point connecting two equivalent rooflike trimer structures. All dimers form planar four-membered Zn2O2 rings whereas the monomer has a planar CZnOC backbone. A hexameric drumlike structure represents the global minimum on the potential energy hypersurface of [(MeZnOMe)6]. The enthalpies and Gibbs energies of the related dissociation reactions hexamer --> tetramer --> trimer --> dimer --> monomer as well as of a number of isomerization reactions have been calculated. The complexes [(MeZnOMe)n] (n = 1-3) form adducts with Lewis bases such as tetrahydrofuran (thf) and pyridine (py). The binding energy of py to the zinc atoms is about 65% larger than that of thf but is not large enough to break up the larger clusters. The bimolecular disproportionation of [(MeZnOMe)4] with formation of the dicubane [Zn{(MeZn)3(OMe)4}2] and Me2Zn is less endothermic than any isomerization or dissociation reaction of the heterocubane, but for steric reasons this reaction is not possible if R' = tBu. A novel reaction mechanism for the reported interconversion, disproportionation and ligand exchange reactions of zinc alkoxide complexes is proposed.     

The lowest-lying electronic singlet and triplet potential energy surfaces for the HNO-NOH system: energetics, unimolecular rate constants, tunneling and kinetic isotope effects for the isomerization and dissociation reactions

The lowest-lying electronic singlet and triplet potential energy surfaces (PES) for the HNO-NOH system have been investigated employing high level ab initio quantum chemical methods. The reaction energies and barriers have been predicted for two isomerization and four dissociation reactions. Total energies are extrapolated to the complete basis set limit applying focal point analyses. Anharmonic zero-point vibrational energies, diagonal Born-Oppenheimer corrections, relativistic effects, and core correlation corrections are also taken into account. On the singlet PES, the (1)HNO → (1)NOH endothermicity including all corrections is predicted to be 42.23 ± 0.2 kcal mol(-1). For the barrierless decomposition of (1)HNO to H + NO, the dissociation energy is estimated to be 47.48 ± 0.2 kcal mol(-1). For (1)NOH → H + NO, the reaction endothermicity and barrier are 5.25 ± 0.2 and 7.88 ± 0.2 kcal mol(-1). On the triplet PES the reaction energy and barrier including all corrections are predicted to be 7.73 ± 0.2 and 39.31 ± 0.2 kcal mol(-1) for the isomerization reaction (3)HNO → (3)NOH. For the triplet dissociation reaction (to H + NO) the corresponding results are 29.03 ± 0.2 and 32.41 ± 0.2 kcal mol(-1). Analogous results are 21.30 ± 0.2 and 33.67 ± 0.2 kcal mol(-1) for the dissociation reaction of (3)NOH (to H + NO). Unimolecular rate constants for the isomerization and dissociation reactions were obtained utilizing kinetic modeling methods. The tunneling and kinetic isotope effects are also investigated for these reactions. The adiabatic singlet-triplet energy splittings are predicted to be 18.45 ± 0.2 and 16.05 ± 0.2 kcal mol(-1) for HNO and NOH, respectively. Kinetic analyses based on solution of simultaneous first-order ordinary-differential rate equations demonstrate that the singlet NOH molecule will be difficult to prepare at room temperature, while the triplet NOH molecule is viable with respect to isomerization and dissociation reactions up to 400 K. Hence, our theoretical findings clearly explain why (1)NOH has not yet been observed experimentally.     

Ab initio studies on mechanisms of cycloaddition reaction between dichloro-vinylidene and ethylene

Mechanisms of cycloaddition reaction between singlet dichloro-vinylidene (R1) and ethylene (R2) have been investigated with MP2 and B3LYP /6-31G* methods, including geometry optimization, vibrational analysis, and energy for the involved stationary points on the potential energy surface. CCSD(T)/6-31G* single point calculations are also applied to the geometries from both methods. CCSD(T) relative energies for the stationary points predicted by MP2 and B3LYP are in a very good agreement. Three reaction pathways are located. Along the first one, one intermediate (INT1) is firstly generated, which then rearranges into a three-membered ring compound (P1) via a small barrier of 5.4 kJ/mol; the other two paths share the other intermediate (INT2), which isomerizes into two four-membered ring compounds (P2 and P3) via a chlorine and a hydrogen transfer, respectively. The barriers of the latter two paths are significantly higher by approximately 25 kJ/mol than that of the former (27.2 and 28.8 vs 5.4 kJ/mol), the major reaction is therefore the formation of P1.     

Theoretical study on the structures, isomerization and stability of SiC3H isomers

The calculations of the geometry optimizations, energies, dipole moments, vibrational spectra, rotational constants, and isomerization of doublet SiC₃H species were performed using density functional theory and ab initio methods. Four types of isomers, a total of 18 minima, connected by 16 interconversion transition states, were located on the potential energy surface (PES) at the B3LYP/6-311G (d, p) level. More accurate energies were obtained at the CCSD(T)/6-311G(2df, 2p), and G3(MP2) levels. With the highest isomerization barrier, the lowest lying structure, linear A1 possesses the largest kinetic stability. Besides, the isomerization barriers of A2, A4, C2, F1, F4 and F5 are over 10 kcal/mol, and these isomers are also considered to be higher kinetically stable. Other isomers cannot be kinetically stabilized with considerably low isomerization barriers. Investigation on the bonding properties and the computations of vibrational spectra, dipole moments, and rotational constants for SiC₃H isomers are helpful for understanding their structures and also valuable for their detections in the interstellar space and laboratory.     

Ab initio study on the mechanism of forming a germanic hetero-polycyclic compound between alkylidenegermylene and ethylene

The mechanism of the cycloaddition reaction of forming a germanic hetero-polycyclic compound between singlet alkylidenegermylene and ethylene has been investigated with MP2/6-31G* method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. The energies of the different conformations are calculated by CCSD(T)//MP2/6-31G* method. From the surface energy profile, it can be predicted that the dominant reaction pathway for this reaction consists of three steps: the two reactants first form a three-membered ring intermediate INT1 through a barrier-free exothermic reaction of 35.4 kJ/mol; this intermediate then isomerizes to an active four-membered ring product P2.1 via a transition-state TS2.1 with a barrier of 57.6 kJ/mol; finally, P2.1 further reacts with ethylene to form the germanic hetero-polycyclic compound P3, for which the barrier is only 0.8 kJ/mol. The rate of this reaction path considerably differs from other competitive reaction paths, indicating that the cycloaddition reaction has an excellent selectivity.     

Theoretical Study on Mechanism of Cycloadditional Reaction Between Dichloro-Germylidene and Formaldehyde

Mechanism of the cycloadditional reaction between singlet dichloro-germylidene and formaldehyde has been investigated with MP2/6-31G^* method, including geometry opti-mization, vibrational analysis and energies for the involved stationary points on the poten-tial energy surface. From the potential energy profile, we predict that the cycloaddition reaction between singlet dichloro-germylidene and formaldehyde has two competitive dom-inant reaction pathways, going with the formation of two side products (INT3 and INT4), simultaneously. Both of the two competitive reactions consist of two steps, two reactants firstly form a three-membered ring intermediate INT1 and a twisted four-membered ring intermediate INT2, respectively, both of which are barrier-free exothermic reactions of 41.5 and 72.3 kJ/mol; then INT1 isomerizes to a four-membered ring product P1 via transition state TS1, and INT2 isomerizes to a chlorine-transfer product P2 via transition state TS2,with the barriers of 2.9 and 0.3 kJ/mol, respectively. Simultaneously, P1 and INT2 further react with formaldehyde to form INT3 and INT4, respectively, which are also barrier-free exothermic reaction of 74.9 and 88.1 kJ/mol.     

A Theoretical Study on the Equilibrium Structures and Relative Stabilities of H_2SO

1INTRODUCTIONTheconformationandvariousphysicalprope-rtiesofdichalcogencompoundsR-X-Y-R?(X,Y=O,S,Se,andTe)havebeenatopicofseveralexperi-mentalandtheoreticalinvestigations[1].CompoundscontainingS=SandS=Obondshavelongbeenpro-posedasintermediatesinorganicsynthesisand,onoccasion,asstableentities[2～18].Despitethecom-monlyusedrepresentationofS–Obondinsulfoxi-desandothersulfinylderivativesasS=O,sulfoxidesareinmanywaysbestdescribedasylideswithhighlypolarizedS–Oσ-bondbecauseofelectro-sta…     

Ab Initio Study on the Mechanism of Cycloaddition Reaction between Silylene Carbene （H_2Si=C：） and Acetone

The mechanism of cycloaddition reaction between singlet silylene carbene and acetone has been investigated with CCSD(T)//MP2/6-31G method. From the potential energy profile, it can be predicted that the reaction has two competitive dominant reaction pathways. One consists of two steps: (1) the two reactants (R1, R2) firstly form a four-membered ring intermediate (INT4) through a barrier-free exothermic reaction of 585.9 kJ/mol; (2) Then intermediate (INT4) isomerizes to CH3-transfer product (P4.1) via a transition state (TS4.1) with energy barrier of 5.3 kJ/mol. The other is as follows: on the basis of intermediate (INT4) created between R1 and R2, intermediate (INT4) further reacts with acetone (R2) to form the intermediate (INT5) through a barrier-free exothermic reaction of 166.3 kJ/mol; Then, intermediate (INT5) isomerizes to a silicic bis-heterocyclic product (P5) via a transition state (TS5), for which the barrier is 54.9 kJ/mol. The presented rule of this reaction: the [2+2] cycloaddition effect between the π orbital of silylene carbene and the π orbital of π-bonded compounds leads to the formation of a four-membered ring intermediate (INT4); The unsaturated property of C atom from carbene in the four-membered ring intermediate (INT4) results in the generation of CH3-transfer product (P4.1) and silicic bis-heterocyclic compound (P5).     

Theoretical study of mechanism of cycloaddition reaction between singlet dichlorosilylene carbene (Cl<Subscript>2</Subscript>Si=C:) and formaldehyde

The mechanism of the cycloaddition reaction between singlet dichlorosilylene carbene (Cl₂Si=C:) and formaldehyde has been investigated with MP2/6-31G* method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. The energies of the different conformations are calculated by Zero-point energy and CCSD (T)//MP2/6-31G* method. From the potential energy profile, it can be predicted that the reaction has two competitive dominant reaction pathways. The first dominant reaction pathway consists of two steps: (1) the two reactants (R1, R2) firstly form a four-membered ring intermediate (INT4) through a barrier-free exothermic reaction of 387.9 kJ/mol; (2) intermediate (INT4) then isomerizes to H-transfer product (P4.2) via a transition state (TS4.2) with energy barrier of 4.7 kJ/mol. The second dominant reaction pathway as follows: on the basis of intermediate (INT4) created between R1 and R2, intermediate (INT4) further reacts with formaldehyde (R2) to form the intermediate (INT5) through a barrier-free exothermic reaction of 158.3 kJ/mol. Then, intermediate (INT5) isomerizes to a silicic bis-heterocyclic product (P5) via a transition state (TS5), for which the barrier is 40.1 kJ/mol.     

1-Germavinylidene (Ge═CH2), germyne (HGeCH), and 2-germavinylidene (H2Ge═C) molecules and isomerization reactions among them: anharmonic rovibrational analyses

Theoretical investigations of three equilibrium structures and two associated isomerization reactions of the GeCH(2) - HGeCH - H(2)GeC system have been systematically carried out. This research employed ab initio self-consistent-field (SCF), coupled cluster (CC) with single and double excitations (CCSD), and CCSD with perturbative triple excitations [CCSD(T)] wave functions and a wide variety of correlation-consistent polarized valence cc-pVXZ and cc-pVXZ-DK (where X = D, T, Q) basis sets. For each structure, the total energy, geometry, dipole moment, harmonic vibrational frequencies, and infrared intensities are predicted. Complete active space SCF (CASSCF) wave functions are used to analyze the effects of correlation on physical properties and energetics. For each of the equilibrium structures, vibrational second-order perturbation theory (VPT2) has been utilized to obtain the zero-point vibration corrected rotational constants, centrifugal distortion constants, and fundamental vibrational frequencies. The predicted rotational constants and anharmonic vibrational frequencies for 1-germavinylidene are in good agreement with available experimental observations. Extensive focal point analyses, including CCSDT and CCSDT(Q) energies and basis sets up to quintuple zeta, are used to obtain complete basis set (CBS) limit energies. At all levels of theory employed in this study, the global minimum of the GeCH(2) potential energy surface (PES) is confirmed to be 1-germavinylidene (GeCH(2), 1). The second isomer, germyne (HGeCH, 2) is predicted to lie 40.4(41.1) ± 0.3 kcal mol(-1) above the global minimum, while the third isomer, 2-germavinylidene (H(2)GeC, 3) is located 92.3(92.7) ± 0.3 kcal mol(-1) above the global minimum; the values in parentheses indicate core-valence and zero-point vibration energy (ZPVE) corrected energy differences. The barriers for the forward (1→2) and reverse (2→1) isomerization reactions between isomers 1 and 2 are 48.3(47.7) ± 0.3 kcal mol(-1) and 7.9(6.6) ± 0.3 kcal mol(-1), respectively. On the other hand, the barriers of the forward (2→3) and reverse (3→2) isomerization reactions between isomers 2 and 3 are predicted to be 55.2(53.2) ± 0.3 kcal mol(-1) and 3.3(1.6) ± 0.3 kcal mol(-1), respectively.     

Quantum mechanical study of the potential energy surface of the ClO + NO2 reaction

In the study of the reaction pathways of the ClO + NO2 reaction including reliable structures of the reactants, products, intermediates, and transition states as well as energies the MP2/6-311G(d), B3LYP/6-311G(d), and G2(MP2) methods have been employed. Chlorine nitrate, ClONO2, is formed by N-O association without an entrance barrier and is stabilized by 29.8 kcal mol(-1). It can undergo either a direct 1,3 migration of Cl or OCl rotation to yield an indistinguishable isomer. The corresponding barriers are 45.8 and 7.1 kcal mol(-1), respectively. ClONO2 can further decompose into NO3 + Cl with an endothermicity of 46.4 kcal mol(-1). The overall endothermicity of the NO2 + ClO --> NO3 + Cl reaction is calculated to be 16.6 kcal mol(-1). The formation of cis-perp and trans-perp conformer of chlorine preoxynitrite, ClOONO(cp) and ClOONO(tp), are exothermic by 5.4 and 3.8 kcal mol(-1), respectively. Calculations on the possible reaction pathways for the isomerization of ClOONO to ClONO2 showed that the activation barriers are too high to account for appreciable nitrate formation from peroxynitrite isomerization. All quoted relative energies are related to G2(MP2) calculations.     

Thermal behavior of [2.1.1]propellane: a DFT/ab initio study

Density functional and ab initio molecular orbital calculations have been used to search for the low energy path of the thermal isomerization of [2.1.1]propellane 1. Three reaction modes were considered: ring opening of the bicyclo[1.1.0]butane unit in 1 to give 1,2-dimethylenecyclobutane 21, opening of the four-membered ring of 1 to afford 1,3-dimethylenecyclobutane 20, and breaking of the [2.1.1]propellane central bond and one of the bicyclo[1.1.0]butane side bonds to form carbene 17. At the CAS(12,12)PT2N/6-31G(d) level of theory, the activation barrier of the latter route was lowest in energy. Further investigation of this process at the QCISD(T)/6-311G(d,p)//QCISD/6-31G(d) and B3PW91/6-311G(d,p)// B3PW91/6-311G(d,p) level of theory indicated that the barrier of isomerization of 1 --> 17 amounts to 29 kcal/mol and that 17 is stabilized by hydrogen migration to give dienes 18 and 19.     

Density functional theory (DFT) and ab initio investigations of the electronic and molecular structures of the monomer and the dimer of trimethylaluminium

The electronic and molecular structures of the monomer and dimer of trimethylalu-minium have been studied using density functional theory and ab initio MP2 method. The optimized geometry of the monomer Al(CH3)3 is of C3h symmetry, whereas that of the dimer [A1(CH3)3]2 contains a carbon-bridged four-membered ring structure with C2h symmetry. The hydrogen-bridged six-membered ring structure is found to be unstable. The calculated dimerization energy for the four-membered ring structure is 78 kJ/mol, in close proximity to the experimental value of 85.27 kJ/mol.     

Theoretical study on the structures, isomerization and stability of SiC4 isomers

The structures, energetics, dipole moments, vibrational spectra, rotational constants, and isomerization of singlet SiC₄ isomers were explored using ab initio methods. Five types of isomers, a total of 11 minima, connected by 11 interconversion transition states, were located on the potential energy surface at the MP2/6-311G(d, p) level. More accurate energies were obtained at the G3(MP2) level. With the highest isomerization barrier, a C2v tetra-angular cone possesses the largest kinetic stability. The lowest-lying structure, linear SiCCCC is also highly kinetically stabilized. Besides, D2d bicyclic c-Si(CC)₂, C2v five-membered ring c-SiCCCC, another C2v tetra-angular cone isomer and C3v trigonal bipyramid isomer are also considered to be kinetically stable, because their isomerization barriers are all over 10 kcal/mol. Other isomers cannot be kinetically stabilized with considerably low isomerization barriers. Investigation on the vibrational spectra, dipole moments, and rotational constants for SiC₄ isomers are valuable for their detections in the interstellar space and laboratory.     

Theoretical study on the ring-opening isomerization reaction mechanism of the ring isomers of N8H8

Ring-opening isomerization from ring-shaped isomers to chain-shaped isomers of N(8)H(8) has been studied by a density function B3LYP method at 6-311+ +G** level. 20 ring-shaped isomers have been found to be able to transform into chain-shaped isomers, with 20 possible transition states got by ring-opening structure optimization. Furthermore, the ring-openings have been found in the longer N-N single bond by analyzing the length change of N-N bond of ring-shaped isomers in ring-opening processes. In addition, with the activation energies in ring-opening processes, the differences of the activation energies in isomerization between the isomers have been found according to the classification of rings. The activation energies in ring-opening isomerization of six-membered ring-shaped isomers are higher than that of the four-membered ring-shaped isomers. It indicates that six-membered ring-shaped isomers difficult in ring-opening in the isomerization are the steadiest ring-shaped isomers of N(8)H(8) while four-membered ring-shaped isomers easy in ring-opening are the most unstable.     

Theoretical study on reactions of nitroethylene with the Si(100)-2 x 1 surface

Possible reaction pathways of nitroethylene with the Si(100)-2 x 1 surface have been investigated by unrestricted density functional theory. The facile occurrence of the studied reactions was demonstrated by the low activation energies of the rate-determining steps (1.07-5.23 kcal/mol). It was found that the [4 + 2] cycloaddition reaction of nitroethylene is most kinetically favorable. The isomerization reactions of the addition products were also investigated. The [3 + 2] cycloaddition product may further undergo a rearrangement by overcoming a 12.37 kcal/mol activation energy barrier into an isomer, with an oxygen atom of the nitryl group inserted between two silicon atoms of the Si(100) surface.