PCM study of the solvent and substituent effects on bond dissociation energies of the C?NO bond

Quantum chemical calculations are used to estimate the equilibrium C? NO bond dissociation energies (BDEs) for eight X? NO molecule (X = CCl₃, C₆F₅, CH₃, CH₃CH₂, iC₃H₇, tC₄H₉, CH₂CHCH₂, and C₆H₅CH₂). These compounds are studied by employing the hybrid density functional theory (B3LYP, B3PW91, B3P86) methods together with 6‐31G** and 6‐311G** basis sets and the complete basis set (CBS‐QB3) method. The obtained results are compared with the available experimental results. It is demonstrated that B3P86/6‐31G** and CBS‐QB3 methods are accurate for computing the reliable BDEs for the X? NO molecule. Considering the inevitably computational cost of CBS‐QB3 method and the reliability of the B3P86 calculations, B3P86 method with 6‐31G** basis set may be more suitable to calculate the BDEs of the C? NO bond. The solvent effects on the BDEs of the C? NO bond are analyzed and it is shown that the C? NO BDEs in a vacuum computed by using B3PW91/6‐311G** method are the closest to the computed values in acetontrile and the average solvent effect is 1.48 kcal/mol. Subsequently, the substituent effects of the BDEs of the C? NO bond are further analyzed and it is found that electron denoting group stabilizes the radical and as a result BDE decreases; whereas electron withdrawing group stabilizes the group state of the molecule and thus increases the BDE from the parent molecule. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Quantum-Chemical Study of Alkyl Carbenium Ions in 100% Sulfuric Acid

A quantum-chemical study of neutral and protonated monoalkyl sulfates RHSO₄and [RH₂SO₄]⁺(where R = CH₃, C₂H₅, iso-C₃H₇, and tert-C₄H₉) is carried out. Calculations are performed using the Hartree–Fock method in the 6-31G** and 6-31++G** basis sets taking into account electron correlation according to the Müller–Plesset perturbation theory MP2/6-31+G*//6-31+G*. Protonated tert-butyl sulfate was also calculated by the DFT B3LYP/6-31++G** method. It was found that monoalkyl sulfates are covalent compounds, and the complete abstraction of alkyl carbenium ions from them has substantial energy cost: 196.4, 161.7, 150.8 and 136.0 kcal/mol, respectively. Protonated methyl and ethyl sulfates are also covalent compounds according to the calculation. They have lower but still high energies of heterolytic dissociation (65.0 and 33.5 kcal/mol, respectively). The energy of R⁺abstraction from protonated isopropyl sulfate is much lower: 23.6 kcal/mol. The main covalent state and the ion–molecular pair, which is a carbenium ion [C(CH₃)₂H]⁺solvated by the H₂SO₄molecule, have about the same energy. The ground state of protonated tert-butyl sulfate corresponds to the ion–molecular complex [C(CH₃)⁺ ₃H₂SO₄] with still lower energy of carbenium ion [C(CH₃)₃]⁺abstraction, which is equal to 10.0 kcal/mol. Calculation according to the DFT B3LYP/6-31++G** method shows the absence of a minimum for the protonated tert-butyl sulfate with a covalent structure on the potential energy surface.     

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     

Theoretical study on the reaction mechanism of (CH3)3CO+CO

The reaction mechanism of (CH₃)₃CO with CO has been theoretically investigated using density-functional theory (DFT) calculations at B3LYP/6-31G* level. In order to get more reliable energy values the single-point energy is evaluated at CCSD (T)/6-31++G** level. The results show that the reaction is multi-channel and the reaction of (CH₃)₃CO radical with CO mostly produces (CH₃)₃C + CO₂. The reaction could play a role in eliminating air pollution.     

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     

The Rearrangement of 2,2‐Diphenyl‐1‐[(E)‐2‐phenylethenyl]cyclopropane to 3,4,4‐Triphenylcyclopent‐1‐ene: a DFT Analysis

A computational study on the rearrangement of 2,2‐diphenyl‐1‐[(E)‐2‐phenylethenyl]cyclopropane ( 1 ) is presented, using density functional theory (DFT), (U)B3LYP with the 6‐31G* basis set (DFT1) and (U)M05‐2X with the 6‐311+G** basis set (DFT2). In agreement with a biradical character of the transition structure (TS) or intermediate, the potential‐energy hypersurface is lowered by the influence of three conjugated Ph groups. Surprisingly, two conformations of the geminal diphenyl group (different twist angles) induce two different minimum‐energy pathways for the rearrangement. Independent of the functional used, the first hypersurface harbors true biradical intermediates, whereas the second energy surface is a flat, slightly ascending slope from the starting material to the TS. The functional (U)M05‐2X with the basis set 6‐311+G** provides realistic energies which seem to be close to experiment. The activation energy for racemization of enantiomers of 1 is lower than that of rearrangement by 2.5 kcal mol^?1, in agreement with experiment.     

Theoretical study of stabling function of the NO to the (CH3)3CO · radical

The stabling function of the NO to the (CH₃)₃CO · radical has been theoretically investigated. Density functional theory (DFT) calculations are performed to optimize the geometries of relevant species. The single‐point energy is evaluated at CCSD(T)/6‐31++G** level. Three reaction channels of (CH₃)₃CO · + NO in the singlet state are considered. The calculations indicate that NO is a stable reagent of active radical (CH₃)₃CO. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Kinetic parameters for the reaction of hydroxyl radical with CH3OCH2F (HFE‐161) in the temperature range of 200–400 K: Transition state theory and Ab initio calculations

Rate coefficients for the reaction of the hydroxyl radical with CH₃OCH₂F (HFE‐161) were computed using transition state theory coupled with ab initio methods, viz., MP2, G3MP2, and G3B3 theories in the temperature range of 200–400 K. Structures of the reactants and transition states (TSs) were optimized at MP2(FULL) and B3LYP level of theories with 6‐31G* and 6‐311++G** basis sets. The potential energy surface was scanned at both the level of theories. Five different TSs were identified for each rotamer. Calculations of Intrinsic reaction coordinates were performed to confirm the existence of all the TSs. The kinetic parameters due to all different TSs are reported in this article. The rate coefficients for the title reaction were computed to be k = (9 ± 1.08) × 10^?13 exp [?(1,713 ± 33)/T] cm³ molecule^?1 s^?1 at MP2, k = (7.36 ± 0.42) × 10^?13 exp [?(198 ± 16)/T] cm³ molecule^?1 s^?1 at G3MP2 and k = (5.36 ± 1.57) × 10^?13 exp [?(412 ± 81)/T] cm³ molecule^?1 s^?1 at G3B3 theories. The atmospheric lifetimes of CH₃OCH₂F at MP2, G3MP2, and G3B3 level of theories were estimated to be 20, 0.1, and 0.3 years, respectively. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Methanediol decomposition mechanisms: A study considering various ab initio approaches

The alternative decomposition reactions CH₂(OH)₂ → CH₂O + H₂O and CH₂(OH)₂ + H₂O → CH₂O + 2H₂O are investigated using the semiempirical PM 3 as well as the ab initio HF /3-21G , HF /6-31G , HF /6-31G **, and MP 2/6-31G ** calculations. Reactants, products, and appropriate transition states are located on corresponding potential energy surfaces and compared with those reported in earlier studies. © 1996 John Wiley & Sons, Inc.     

Mechanism for the gas‐phase reaction between OH and 3‐methylfuran: A theoretical study

The mechanism for the OH + 3‐methylfuran reaction has been studied via ab initio calculations to investigate various reaction pathways on the doublet potential energy surface. Optimizations of the reactants, products, intermediates, and transition structures are conducted using the MP2 level of theory with the 6‐311G(d,p) basis set. The single‐point electronic energy of each optimized geometry is refined with G3MP2 and G3MP2B3 calculations. The theoretical study suggests that the OH + 3‐methylfuran reaction is dominated by the formation of HC(O)CH?C(CH₃)CHOH (P7) and CH(OH)CH?C(CH₃)C(O)H (P9), formed from two low‐lying adducts, IM1 and IM2. The direct hydrogen abstraction pathways and the S_N2 reaction may play a minor or negligible role in the overall reaction of OH with 3‐methylfuran. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Rearrangements in Model Peptide‐Type Radicals via Intramolecular Hydrogen‐Atom Transfer

Intramolecular H‐atom transfer in model peptide‐type radicals was investigated with high‐level quantum‐chemistry calculations. Examination of 1,2‐, 1,3‐, 1,5‐, and 1,6[C ? N]‐H shifts, 1,4‐ and 1,7[C ? C]‐H shifts, and 1,4[N ? N]‐H shifts (Scheme 1), was carried out with a number of theoretical methods. In the first place, the performance of UB3‐LYP (with the 6‐31G(d), 6‐31G(2df,p), and 6‐311+G(d,p) basis sets) and UMP2 (with the 6‐31G(d) basis set) was assessed for the determination of radical geometries. We found that there is only a small basis‐set dependence for the UB3‐LYP structures, and geometries optimized with UB3‐LYP/6‐31G(d) are generally sufficient for use in conjunction with high‐level composite methods in the determination of improved H‐transfer thermochemistry. Methods assessed in this regard include the high‐level composite methods, G3(MP2)‐RAD, CBS‐QB3, and G3//B3‐LYP, as well as the density‐functional methods B3‐LYP, MPWB1K, and BMK in association with the 6‐31+G(d,p) and 6‐311++G(3df,3pd) basis sets. The high‐level methods give results that are close to one another, while the recently developed functionals MPWB1K and BMK provide cost‐effective alternatives. For the systems considered, the transformation of an N‐centered radical to a C‐centered radical is always exothermic (by 25 kJ ? mol^?1 or more), and this can lead to quite modest barrier heights of less than 60 kJ ? mol^?1 (specifically for 1,5[C ? N]‐H and 1,6[C ? N]‐H shifts). H‐Migration barriers appear to decrease as the ring size in the transition structure (TS) increases, with a lowering of the barrier being found, for example when moving from a rearrangement proceeding via a four‐membered‐ring TS (e.g., the 1,3[C ? N]‐H shift, CH₃? C(O)? NH^. → ^.CH₂? C(O)? NH₂) to a rearrangement proceeding via a six‐membered‐ring TS (e.g., the 1,5[C ? N]‐H shift, ^.NH? CH₂? C(O)? NH? CH₃ → NH₂? CH₂? C(O)? NH? CH₂^.).     

1,2,3-三氮杂苯－(水)3复合物多体相互作用

The interaction between 1,2,3-triazine and three water molecules was studied using density functional theory B3LYP method at 6-31-t＋＋G^＊＊ basis set. Various structures for 1,2,3-triazine-（water）n （n= 1, 2, 3） complex were investigated and the different lower energy structures were reported. Many-body analysis was also carded out to obtain relaxation energy and many-body interaction energy （two, three, and four-body）, and the most stable conformer has the basis set superposition error corrected interaction energy of -- 102.61 kJ/mol. The relaxation energy, two- and three-body interactions have significant contribution to the total interaction energy whereas four-body interaction was very small for 1,2,3-triazine-（water）3 complex.     

Methane elimination from n‐pentane

The quantum chemical calculations at the different levels of theory were performed with the target being to determine the vibration frequencies and to estimate the barriers to internal rotations of n‐pentane molecules. In connection with the observed losses of CH₃ and CH₄ from the n‐pentane in gas phase, the calculations at the B3LYP level of theory with the 6‐31G(d) basis set were used to study the ground‐state potential energy surface of the n‐pentane. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Structural and spectroscopic studies on 2-pyranones

Experimental methods of infrared, Raman and electronic absorption spectroscopy and DFT calculations using B3LYP functionals and 6-31G** and 6-311++G** basis sets have been used to understand the structural and spectral characteristics of 2-pyranones, 6-phenyl-4-methylsulfanyl-2-oxo-2H-pyran and 6-phenyl-4-methylsulfanyl-2-oxo-2H-pyran-3-carbonitrile in the electronic ground (S₀) and first excited (S₁) states. Information about the size, shape, charge density distribution and site of chemical reactivity of the molecules has been obtained by mapping electron density isosurface with electrostatic potential surfaces (ESP). Based on TD-DFT calculations using 6-31+G**5D basis set, an assignment of absorption peaks in the UV–VIS region has been suggested. The S₁ state is found to be a ¹(π,π*) state. A complete vibrational analysis has been attempted on the basis of experimental infrared and Raman spectra and calculated frequency and intensity of the vibrational bands and potential energy distribution over the internal coordinates. Characteristic vibrational bands of the 2-pyranone ring and methylsulfanyl and carbonyl groups have been identified.     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

DFT investigation of sites of protonation of antitumor of 2‐(4‐aminophenyl)benzazoles in the gas phase and in solution

The proton affinity on each of the possible sites in the antitumor 2‐(4‐aminophenyl)benzazoles has been calculated at the B3LYP/6‐311G** level of theory in the gas phase and in solution. The N₃‐site of protonation is found to be strongly favored over the NH₂‐site for the studied compounds both in gas phase and in solution. The stability of N₃‐protonated species is explained by the resonance interaction of the NH₂‐group with the heterocyclic ring. The potential energy surface (PES) for the protonation process was studied at the density functional theory (DFT)/B3LYP/6‐311++G** level of theory. Solvent effects on the PES were also examined using two models: Onsager self‐consistent field and polarizable continuum model (PCM). © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

The 1:1 glycine zwitterion-water complex: An ab initio electronic structure study

More than a dozen stationary points on the potential energy surface for the 1:1 glycine zwitterion—water complex have been investigated at Hartree-Fock or MP2 levels of theory with basis sets ranging from split valence (4-31G) to split valence plus polarization and diffuse function (6–31 + + G**) quality. Only one true minimum (GLYZWM, C₁ symmetry) could be located on the potential energy surface. GLYZWM features a bridged water molecule acting as both a hydrogen bond acceptor and donor with the NH₃⁻ and CO₂⁻ units of the glycine zwitterion. The total hydrogen bond energy in GLYZWM is computed as 16 kcal/mol (MP2/6–31 ++ G** // 6–31 ++ G**, including corrections for basis set superpositions errors). The computed vibrational frequencies and normal mode forms of the GLYZWM complex resemble in many cases experimental assignments made for the glycine zwitterion in bulk water on the basis of Raman spectroscopy. © 1996 by John Wiley & Sons, Inc.     

Ab initio and RRKM calculations for the reaction channels of O(1D) + CH3CHF2

Ab initio and Rice–Ramsperger–Kassel–Marcus theories are carried out to study the potential energy surface and the energy‐dependent rate constants and branching ratios of the products for O(¹D) + CH₃CHF₂ reaction. Optimized geometries and vibrational frequencies have been obtained by MP2/6‐311G(d,p) method. The main products of the title reaction are CH₃CFO + HF, CH₂CFOH + HF, and CH₃ + CF₂OH at lower collision energy; and CH₃ + CF₂OH, CH₃CF₂ + OH are the main products at higher collision energy. CHF₂ + CH₂OH are the main products in the whole range of collision energy. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012