Correlated one‐electron wave functions

Using four basis sets, 6‐311G(d,p), 6‐31+G(d,p), 6‐311++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the acidic H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. By contrast with above the three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen‐bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is a noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD(T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Long‐range π‐type hydrogen bond in dimers CH2?HF,CH2O?H2O,and CH2O?NH3

Using four basis bets, (6‐311G(d,p), 6‐31+G(d,p), 6‐31++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for the dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. In contrast with the above three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD (T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

DFT studies on the mechanism of the reaction of C2H5S with NO2

The mechanisms for the reaction of C₂H₅S with NO₂ are investigated at the QCISD(T)/6‐311++G(d, p)//B3LYP/6‐311++G(d, p) level on both single and triple potential energy surfaces. The geometries, vibrational frequencies and zero‐point energy (ZPE) corrections of all stationary points involved in the title reaction are calculated at the B3LYP/6‐311++G(d, p) level. The results show that the reaction is more predominant on the single potential energy surface, while it is negligible on the triple potential energy surface. Without barrier height in the whole process, the major channel is R → C₂H₅SONO (IM1 and IM2) → P1 (C₂H₅SO+NO). With much heat released in the formation of C₂H₅SNO₂ (IM3) and the transition state involved in the subsequent step more stable than reactants, P4 (CH₃CHS + t‐HONO) is subdominant product energetically. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Density functional theory investigation of the binding interactions between phosphoryl,carbonyl, imino,and thiocarbonyl ligands and the pentaaqua nickel(II) complex: Coordination affinity and associated parameters

Density Functional Theory (UB3LYP/6‐311++G(d,p)) calculations of the affinity of the pentaaqua nickel(II) complex for a set of phosphoryl [O?P(H)(CH₃)(PhR)], imino [HN?C(CH₃)(PhR)], thiocarbonyl [S?C(CH₃)(PhR)] and carbonyl [O?C(CH₃)(PhR)] ligands were performed, where R?NH₂, OCH₃, OH, CH₃, H, Cl, CN, and NO₂ is a substituent at the para‐position of a phenyl ring.The affinity of the pentaaqua nickel(II) complex for these ligands was analized and quantified in terms of interaction enthalpy (ΔH), Gibbs free energy (ΔG²⁹⁸), geometric and electronic parameters of the resultant octahedral complexes. The ΔH and ΔG²⁹⁸ results show that the ligand coordination strength increases in the following order: carbonyl < thiocarbonyl < imino < phosphoryl. This coordination strength order is also observed in the analysis of the metal‐ligand distances and charges on the ligand atom that interacts with the Ni(II) cation. The electronic character of the substituent R is the main parameter that affects the strength of the metal‐ligand coordination. Ligands containing electron‐donating groups (NH₂, OCH₃, OH) have more exothermic ΔH and ΔG²⁹⁸ than ligands with electron‐withdrawing groups (Cl, CN, NO₂). The metal‐ligand interaction decomposed by means of the energy decomposition analysis (EDA) method shows that the electronic character of the ligand modulates all the components of the metal‐ligand interaction. The absolute softness of the free ligands is correlated with the covalent contribution to the instantaneous interaction energy calculated using the EDA method. © 2013 Wiley Periodicals, Inc.     

Isomerizations of CH3SO2: Quantum Chemistry and Topological Study

The structures, stabilities and the isomerization reactions of CH₃SO₂ isomers in a doublet electronic state have been studied at B3LYP/6‐311+ +G (d,p), MP2/6‐311++G (d,p) and CCSD(T)/6‐311++G (d,p) levels. The three different levels of calculation give the similar results: thirteen minimum isomers were located and they were connected by eleven transition states. Among the thirteen isomers, cis‐CH₃OSO, trans‐CH₃OSO and CH₃SO₂ are the most stable species, and they should be detected easily in experiment. This is well consistent with the experimental result. These isomers could isomerize to each other by chemical bond vibration, chemical bond rotation and atom migration. The non‐planar ring structure transition state (STS), which was found in this paper, extended the concept of ring STS to the non‐planar systems.     

DFT studies on the multi‐channel reaction of CH3S+NO2

The mechanisms for the reaction of CH₃S with NO₂ are investigated at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311++G(d,p) on both single and triple potential energy surfaces (PESs). The geometries, vibrational frequencies, and zero‐point energy (ZPE) correction of all stationary points involved in the title reaction are calculated at the B3LYP/6‐311++G(d,p) level. More accurate energies are obtained at the QCISD(T)/6‐311++G(d,p). The results show that 5 intermediates and 14 transition states are found. The reaction is more predominant on the single PES, while it is negligible on the triple PES. Without any barrier height for the whole process, the main channel of the reaction is to form CH₃SONO and then dissociate to CH₃SO+NO. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Complexation of Phenol and Thiophenol by Amine N‐Oxides: Isothermal Titration Calorimetry and ab Initio Calculations

To develop a new solvent‐impregnated resin (SIR) system for removal of phenols from water, the complex formation of dimethyldodecylamine N‐oxide (DMDAO), trioctylamine N‐oxide (TOAO), and tris(2‐ethylhexyl)amine N‐oxide (TEHAO) with phenol (PhOH) and thiophenol (PhSH) is studied. To this end we use isothermal titration calorimetry (ITC) and quantum chemical modeling (on B3LYP/6‐311G(d,p)‐optimized geometries: B3LYP/6‐311+G(d,p), B3LYP/6‐311++G(2d,2p), MP2/6‐311+G(d,p), and spin component scaled (SCS) MP2/6‐311+G(d,p); M06‐2X/6‐311+G(d,p)//M06‐2X/6‐311G(d,p), MP2 with an extrapolation to the complete basis set limit (MP2/CBS), as well as CBS‐Q). The complexes are analyzed in terms of structural (e.g., bond lengths) and electronic elements (e.g., charges). Furthermore, complexation and solvent effects (in benzene, toluene, and mesitylene) are investigated by ITC measurements, yielding binding constants K, enthalpies ΔH⁰, Gibbs fre energies ΔG⁰, and entropies ΔS⁰ of complex formation, and stoichiometry N. The ITC measurements revealed strong 1:1 complex formation between both DMDAO–PhOH and TOAO–PhOH. The binding constant (K=1.7–5.7×10⁴ M ^?1) drops markedly when water‐saturated toluene was used (K=5.8×10³ M ^?1), and π–π interaction with the solvent is shown to be relevant. Quantum mechanical modeling confirms formation of stable 1:1 complexes with linear hydrogen bonds that weaken on attachment of electron‐withdrawing groups to the amine N‐oxide moiety. Modeling also showed that complexes with PhSH are much weaker than those with PhOH, and in fact too weak for ITC determination. CBS‐Q incorrectly predicts equal or even higher binding enthalpies for PhSH than for PhOH, which invalidates it as a benchmark for other calculations. Data from the straightforward SCS‐MP2 method without counterpoise correction show very good agreement with the MP2/CBS values.     

Bond dissociation energies and radical heats of formation in CH3Cl,CH2Cl2, CH3Br,CH2Br2, CH2FCl,and CHFCl2

An analysis of thermochemical and kinetic data on the bromination of the halomethanes CH_4–nX_n (X = F, Cl, Br; n = 1–3), the two chlorofluoromethanes, CH₂FCl and CHFCl₂, and CH₄, shows that the recently reported heats of formation of the radicals CH₂Cl, CHCl₂, CHBr₂, and CFCl₂, and the C? H bond dissociation energies in the matching halomethanes are not compatible with the activation energies for the corresponding reverse reactions. From the observed trends in CH₄ and the other halomethanes, the following revised ΔH°_f,298 (R) values have been derived: ΔH°_f(CH₂Cl) = 29.1 ± 1.0, ΔH°_f(CHCl₂) = 23.5 ± 1.2, ΔH_f(CH₂Br) = 40.4 ± 1.0, ΔH°_f(CHBr₂) = 45.0 ± 2.2, and ΔH°_f(CFCl₂) = ?21.3 ± 2.4 kcal mol^?1. The previously unavailable radical heat of formation, ΔH°_f(CHFCl) = ?14.5 ± 2.4 kcal mol^?1 has also been deduced. These values are used with the heats of formation of the parent compounds from the literature to evaluate C? H and C? X bond dissociation energies in CH₃Cl, CH₂Cl₂, CH₃Br, CH₂Br₂, CH₂FCl, and CHFCl₂.     

Definitive heat of formation of methylenimine,CH2NH,and of methylenimmonium ion,CH2NH2+, by means of W2 theory

A long‐standing controversy concerning the heat of formation of methylenimine has been addressed by means of the W2 (Weizmann‐2) thermochemical approach. Our best calculated values, ΔH°_f,298(CH₂NH) = 21.1±0.5 kcal/mol and ΔH°_f,298(CH₂NH₂⁺) = 179.4±0.5 kcal/mol, are in good agreement with the most recent measurements but carry a much smaller uncertainty. As a byproduct, we obtain the first‐ever accurate anharmonic force field for methylenimine: upon consideration of the appropriate resonances, the experimental gas‐phase band origins are all reproduced to better than 10 cm^?1. Consideration of the difference between a fully anharmonic zero‐point vibrational energy and B3LYP/cc‐pVTZ harmonic frequencies scaled by 0.985 suggests that the calculation of anharmonic zero‐point vibrational energies can generally be dispensed with, even in benchmark work, for rigid molecules. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1297–1305, 2001     

Thermochemistry,Reaction Paths,and Kinetics on the Secondary Isooctane Radical Reaction with 3O2

2,2,4‐Trimethylpentane, also known as isooctane, is used as one of the model fuel species on spark and homogeneous charge compression ignition engines. This study presents thermochemical and kinetic properties in the oxidation of the secondary isooctane radical, which includes the peroxy radical formed from the ³O₂ association, the hydroperoxy alkyl radicals formed from the intramolecular hydrogen transfers, and the products formed from reactions of the hydroperoxy alkyl radicals. Geometries, vibration frequencies, internal rotor potentials, and thermochemical properties, Δ_fH, S°(T), and C°_p(T) (5 K ≤ T ≤ 5000 K) were calculated at the individual B3LYP/6–31G(d,p) and the composite CBS‐QB3 calculation method. The standard enthalpies of formation at 298 K were evaluated using isodesmic reaction schemes with several work reactions for each species. Entropy and heat capacities were determined using geometric parameters and frequencies from the B3LYP/6–31G(d,p) calculations for the lowest energy conformer. Internal rotor barriers were determined. Application of group additivity with comparison to calculated values is also illustrated. Transition states and kinetic parameters for intramolecular hydrogen atom transfer and molecular elimination channels were characterized to evaluate reaction paths and kinetics. Kinetic parameters were determined versus pressure and temperature for the chemical activated formation and unimolecular dissociation of the peroxide adduct. Multifrequency quantum Rice–Ramsperger–Kassel analysis was used for k(E) with master equation analysis for falloff. The kinetic analysis shows that the main reaction channels are the formation of isooctene ((CH3)₃CCH=C(CH₃)₂) + HO₂^?, and the formation of the cyclic: (CH₃)₂‐y(CCH₂CHO)‐(CH₃)₂, (CH₃)₃C‐y(CHCO)‐(CH₃)₂, and (CH₃)₃C‐y(CHCHCH₂O)‐(CH₃).     

Hetero‐π‐systems from 2 + 2 cycloreversion,part 2.1Ab initio thermochemical study of heterocyclobutanes 2 + 2 cycloreversion to form heteroethenes H2C=X (X=NH,O, SiH2, PH,S)

Ab initio and DFT thermochemical study of diradical mechanism of 2 + 2 cycloreversion of parent heterocyclobutanes and 1,3‐diheterocyclobutanes, cyclo‐(CH₂CH₂CH₂X), and cyclo‐(CH₂XCH₂X), where X = NH, O, SiH₂, PH, S, was undertaken by calculating closed‐shell singlet molecules at three levels of theory: MP4/6‐311G(d)//MP2/6‐31G(d)+ZPE, MP4/6‐311G(d,p)//MP2/6‐31G (d,p)+ZPE, and B3LYP/6‐311+G(d,p)+ZPE. The enthalpies of 2 + 2 cycloreversion decrease on going from group 14 to group 16 elements, being substantially higher for the second row elements. Normally endothermic 2 + 2 cycloreversion is predicted to be exothermic for 1,3‐diazetidine and 1,3‐dioxtane. Strain energies of the four‐membered rings were calculated via the appropriate homodesmic reactions. The enthalpies of ring opening via the every possible one‐bond homolysis that results in the formation of the corresponding 1,4‐diradical were found by subtracting the strain energies from the central bond dissociation energies of the heterobutanes CH₃CH₂—CH₂XH, CH₃CH₂—XCH₃, and HXCH₂—XCH₃. The latter energies were determined via the enthalpies of the appropriate dehydrocondensation reactions, using C—H and X—H bond energies in CH₃XH calculated at G2 level of theory. Except 1,3‐disiletane, in which ring‐opening enthalpy attains 69.7 kcal/mol, the enthalpies of the most economical ring openings do not exceed 60.7 kcal/mol. The 1,4‐diradical decomposition enthalpies found as differences between 2 + 2 cycloreversion and ring‐opening enthalpies were negative, the least exothermicity was calculated for ⋅ CH₂SiH₂CH₂CH₂. The only exception was 1,3‐disiletane, which being diradical, CH₂SiH₂CH₂SiH₂, decomposed endothermically. Since decomposition of the diradical containing two silicon atoms required extra energy, raising the enthalpy of the overall reaction to 78.9 kcal/mol, 1,3‐disiletane was predicted to be highly resisting to 2 + 2 cycloreversion. © 2007 Wiley Periodicals, Inc. Heteroatom Chem 18:704–720, 2007; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20377     

Theoretical prediction regarding structural and thermodynamical characteristics of stable CH3PO2 isomers and unimolecular decomposition mechanisms of species CH3P(O)2, CH3O?PO,and CH2P(O)OH

The detailed isomerization and dissociation reaction potential energy profile of the CH₃PO₂ system was established at the UCCSD(T)/6‐311++G(3df,2p)//UB3LYP/6‐311++G(d,p) level of theory. Seventy minimum isomers were located and connected by 93 optimized interconversion transition states. Furthermore, 32 isomers with high kinetic stability were predicted to be possible candidates for further experimental detection. The bonding nature of the suggested stable isomers was analyzed while their molecular properties including heats of formation, adiabatic ionization potentials, and adiabatic electronic affinities were calculated at the G2, G2(MP2), G3, and CBS‐Q levels. Based on the isomerization and dissociation potential energy surface, possible unimolecular decomposition mechanisms and pathways of the low‐lying molecules CH₃P(?O)₂, CH₃O? P?O, and CH₂?P(?O)OH were discussed. Furthermore, the transition state theory rate constants of the primary unimolecular dissociation channels were also calculated. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

A theoretical study on the mechanism and thermodynamics of ozone-water gas phase reaction

Ozone water reaction including a complex was studied at the MP2/6-311++G(d,p) and CCSD/6-311++G(2df,2p)//MP2/6-311++G(d,p) levels of theory. The interaction between water oxygen and central oxygen of ozone produces stable H₂O-O₃ complex with no barrier. With decomposition of this complex through H-abstraction by O₃ and O-abstraction by H₂O, three possible product channels were found. Intrinsic reaction coordinate, topological analyses of atom in molecule, and vibrational frequency calculation have been used to confirm the preferred mechanism. Thermodynamic data at T = 298.15 K and atmospheric pressure have been calculated. The results show that the production of hydrogen peroxide is the main reaction channel with ΔG = ?21.112 kJ mol^-1.     

Effect of core substituents on the intramolecular exchange interaction in N,N′‐dioxy‐2,6‐diazaadamantane biradical: DFT studies

Density‐functional theory calculations of a series of organic biradicals on the basis of the N,N′‐dioxy‐2,6‐diazaadamantane core with different substituents at carbon atoms adjacent to the nitroxyl groups have been performed by the UB3LYP/6‐311++G(2d,2p) method. Using the breaking symmetry approach, the values of the exchange interaction parameter, J, between the radical centers are calculated. It is shown that the intramolecular exchange interaction for the most part is ferromagnetic in nature, but the J parameter gradually decreases, changing its sign to antiferromagnetic interaction for the last substituent in the following sequence: CF₃(CH₃)COH > CH₂F(H)COH > CH₂OH > H > CBr₃ > CH₂F > CCl₃ > CF₃ > CH₂Br > CH₂Cl > CH₃ > C₂H₅ > C₃H₇ > i‐C₄H₉ > F > Br > OCH₃ > Cl > CH₂C₆H₅. The calculations at the UHSEH1PBE/6‐311++G(2d,2p) level with the most of substituents show nearly the same variation sequence for the J parameter. It is concluded that spin polarization effects in the diazaadamantane cage and a direct through‐space antiferromagnetic exchange interaction between the nitroxyl groups are the main mechanisms contributing to the exchange interaction parameter value in the studied series of compounds. The exchange coupling constant, J, depends on the electronic effects and geometry of the substituents, as well as on their specific interactions with the nitroxyl radical groups.     

磷叶立德CH2PH3和类磷叶立德自由基∙CHPH3优势构型和C—P键结构的研究

采用MP4/6-311++G(d,p)和B3LYP/6-311++G(d,p)对磷叶立德CH₂PH₃和类磷叶立德自由基∙CHPH₃进行构型优化，从电子密度拓扑分析的角度对C—P键的键结构进行了探讨。得到如下结论：类磷叶立德自由基和磷叶立德的C—P键性质类似，但磷叶立德中π键由两个电子形成，类磷叶立德自由基中π键由一个电子形成，所以前者的π性明显，而后者的π性不明显。类磷叶立德自由基中的这个单电子在碳原子附近，垂直于对称面的方向上运动，有p_(C→P)配键的特征，所以类磷叶立德自由基∙CHPH₃中的C—P键比相应的产物∙CH₂PH₂中的C—P键要弱一些。     

ab initio calculations and thermochemical analysis on C1 atom abstractions of chlorine from chlorocarbons and the reverse alkyl abstractions: Cl2 + R· ↔ Cl· + RCl

Thermodynamic properties (ΔH°_f(298), S°₍₂₉₈₎ and Cp(T) from 300 to 1500 K) for reactants, adducts, transition states, and products in reactions of CH₃ and C₂H₅ with Cl₂ are calculated using CBSQ//MP2/6‐311G(d,p). Molecular structures and vibration frequencies are determined at the MP2/6‐311G(d,p), with single‐point calculations for energy at QCISD(T)/6‐311 + G(d,p), MP4(SDQ)/CbsB4, and MP2/CBSB3 levels of calculation with scaled vibration frequencies. Contributions of rotational frequencies for S°₍₂₉₈₎ and Cp(T)'s are calculated based on rotational barrier heights and moments of inertia using the method of Pitzer and Gwinn [1]. Thermodynamic parameters, ΔH°_f(298), S°₍₂₉₈₎, and C_P(T), are evaluated for C₁ and C₂ chlorocarbon molecules and radicals. These thermodynamic properties are used in evaluation and comparison of Cl₂ + R· → Cl· + RCl (defined forward direction) reaction rate constants from the kinetics literature for comparison with the calculations. Data from some 20 reactions in the literature show linearity on a plot of Ea_fwd vs. ΔH_rxn,fwd, yielding a slope of (0.38 ± 0.04) and intercept of (10.12 ± 0.81) kcal/mole. A correlation of average Arrhenius preexponential factor for Cl· + RCl → Cl₂ + R· (reverse rxn) of (4.44 ± 1.58) × 10¹³ cm³/mol‐sec on a per‐chlorine basis is obtained with Ea_Rev = (0.64 ± 0.04) × ΔH_rxn,Rev + (9.72 ± 0.83) kcal/mole, where Ea_Rev is 0.0 if ΔH_rxn,Rev is more than 15.2 kcal/mole exothermic. Kinetic evaluations of literature data are also performed for classes of reactions. Ea_fwd = (0.39 ± 0.11) × ΔH_rxn,fwd + (10.49 ± 2.21) kcal/mole and average A_fwd = (5.89 ± 2.48) × 10¹² cm³/mole‐sec for hydrocarbons: Ea_fwd = (0.40 ± 0.07) × ΔH_rxn,fwd + (10.32 ± 1.31) kcal/mole and average A_fwd = (6.89 ± 2.15) × 10¹¹ cm³/mole‐sec for C₁ chlorocarbons: Ea_fwd = (0.33 ± 0.08) × ΔH_rxn,fwd + (9.46 ± 1.35) kcal/mole and average A_fwd = (4.64 ± 2.10) × 10¹¹ cm³/mole‐sec for C₂ chlorocarbons. Calculation results on the methyl and ethyl reactions with Cl₂ show agreement with the experimental data after an adjustment of +2.3 kcal/mole is made in the calculated negative Ea's. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 548–565, 2000     

The DFT study on mechanisms for NCO with C2H5 reaction

A detailed investigation has been performed at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311+G(d,p) level for the reaction of NCO with C₂H₅ by constructing singlet and triplet potential energy surfaces (PES). The results show that the title reaction is more favorable on the singlet PES than on the triplet PES. On the singlet PES, the initial addition processes are barrierless and release lots of energy. The dominant channel occurs via the fragmentations of the initial adduct C₂H₅NCO and C₂H₅OCN to form C₂H₄ + HNCO and HOCN, respectively. With higher barrier heights, other products such as CH₄ + HNC + CO, CH₃CHNH + CO, CH₃CH + HNCO, and CH₃CN + H₂ + CO are less competitive. On the triplet PES, the entrance reactions surpass significant barriers; therefore, it could be negligible at the normal atmospheric condition. However, the most feasible channel on the triplet PES is the direct hydrogen abstraction channel to form CH₂CH₂ + HNCO. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical investigations on analytical potential energy function and spectroscopic parameters for the state b3Πu of dimer 7Li2

The SAC‐CI (symmetry‐adapted‐cluster configuration‐interaction) method presented in Gaussian 03 program package is applied to investigate the adiabatic potential energy curves (PECs) of ⁷Li₂(b³Π_u). These calculations are performed at numbers of basis sets, such as 6‐311++G(3df,3pd), 6‐311++G(2df,2pd), 6‐311++G(df,pd), D95V++, D95(3df,3pd), D95(d,p), cc‐PVTZ, 6‐311++G and 6‐311++G(d,p). All the ab initio calculated points are fitted to the analytic Murrell‐Sorbie functions and then used to compute the spectroscopic parameters. The analytic potential energy function (APEF) for this b³Π_u state is reported. By comparison, the spectroscopic parameters reproduced by the APEF attained at 6‐311++G(2df,2pd) are found to be very close to the latest experimental findings. With the APEF obtained at the SAC‐CI/6‐311++G(2df,2pd) level of theory, a total of 62 vibrational states is found when J = 0. The complete vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants for these vibrational states are also reported. The reasonable dissociation limit for this state is deduced using the calculated results at present. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical study on the mechanism for the thermal rearrangement reactions of 2‐silylethyl acetate H3SiCH2CH2OCOCH3

The thermal rearrangement mechanisms of 2‐silylethylacetate H₃SiCH₂CH₂OOCCH₃ were investigated by ab initio molecular orbital theory for the first time. All structures of reactant, transition states, and products were located and fully optimized at the B3LYP/6‐311+G(d, p) levels, and harmonic vibrational frequencies for the involved stationary points on the potential energy surface were obtained. The reaction pathways were analyzed and confirmed by intrinsic reaction coordinate (IRC) calculations. Furthermore, atomic charges were determined by using the natural bond orbital (NBO) analysis. The calculational results show that H₃SiCH₂CH₂OOCCH₃ can rearrange thermally in two ways. One is [1,3] rearrangement (Reaction A), in which silyl group transfers from carbon to oxygen(in C? O? C) via a four‐membered ring transition state, forming silyl acetate and ethylene, the other way, [1,5] rearrangement (Reaction B), happens with transferring of silyl group from carbon to oxygen (in C?O) via a six‐membered ring transition state, forming the same products as in Reaction A. The energy barriers of the Reactions A and B were calculated to be 188.9 and 191.6 kJ/mol at the B3LYP/6‐311+G(d,p) levels, respectively. Changes in thermodynamic functions (ΔS, ΔH, and ΔG), equilibrium constant K(T), as well as preexponential factor A(T), and reaction rate constant k(T) in Eyring transition state theory were calculated over a temperature range of 200–1600 K, and then thermodynamic and kinetic properties of the reactions were analyzed. It can be suggested that Reactions A and B are noncompetitive, and both happen only at elevated temperature. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

A DFT study of hydride transfers to the carbonyl oxygen of DDQ

Hydride‐transfer reactions between benzylic substrates and 2,3‐dichloro‐5,6‐dicyano‐1,4‐benzoquinone (DDQ) were investigated by DFT (density functional theory) calculations. The lowest unoccupied molecular orbital of DDQ has the largest extension on two carbonyl oxygens, which comes from two‐step mixing of antisymmetric orbitals of fragment π MOs. Transition‐state (TS) geometries and activation energies of reactions of four benzylic substrates R²? CH₂‐para‐C₆H₄? R¹ (R¹, R² = H and/or OCH₃) with DDQ were calculated. M06‐2X/6‐311(+*)G* was found to be a practical computational method, giving energies and geometries similar to those of M06‐2X/6‐311++G(3df,2pd) and wB97xD/6‐311++G(3df,2pd). For toluene (R¹ = R² = H), an initiation‐propagation model was suggested, and the calculated kinetic isotope effect k(H)/k(D) = 5.0 with the tunnel correction at the propagating step is in good agreement with the experimental value 5.2. A reaction of para‐MeO? C₆H₄? CH₂(OMe) + DDQ + (H₂O)₁₄ → para‐MeO? C₆H₄? C(?O)H + HOMe + DDQH₂ + (H₂O)₁₃ was investigated by M06‐2X/6‐311(+*)G*. Four elementary processes were found and the hydride transfer (TS1) is the rate‐determining step. The hydride transfer was promoted by association with the water cluster. The size of the water cluster, (H₂O)_n, at TS1 was examined. Three models of n = 14, 20, and 26 were found to give similar activation energies. Metal‐free neutral hydride transfers from activated benzylic substrates to DDQ were proposed to be ready processes both kinetically and thermodynamically. © 2015 Wiley Periodicals, Inc.