Computational study on the reaction mechanism of atmospheric oxidation of ethanol with ozone

The reaction of C₂H₅OH and O₃ on the singlet potential energy surface is carried out using the MP2 and CCSD(T)//MP2 theoretical approaches in connection with the 6-311++G(d,p) basis set. Three pre-reactive complexes C1, C2, and C3 are formed between ethanol and ozone at atmospheric pressure and 298.15 K temperature. With variety of the complexes, seven types of product are obtained which four types of them have enough thermodynamic stability. In thermodynamic approach, the most favor product begins with the formation of pre-reactive C2 complex and produces the CH₃CH(OH)₂ + O₂ as final adduct in a process that is computed to be exothermic by ?53.759 kcal/mol and spontaneous reaction by ?51.833 kcal/mol in Gibbs free energy. In kinetic viewpoint, the formation of CH₃COH + cis-H₂O₃ as a final adducts is the most favor path.     

N2O…HOCl复合物结构和性质的理论研究

分别在DFT-B3LYP和MP2/6-311++G**水平上求得HOCl + N2O复合物势能面上的六种(S1, S2, S3, S4, S5和S6)和四种(S1, S2, S4和S5)构型. 频率分析表明,其中的S1和S3为过渡态,其它为稳定构型. 在复合物S3, S5 和S6中,HOCl 单体的σ*(5O-6H)作为质子供体,与N2O单体中作为质子受体的3O原子相互作用,形成氢键结构,而在氢键复合物S2中, 质子受体为N2O单体中的端1N原子;复合物S1中,HOCl分子的σ*(5O-4Cl)作为质子供体与N2O分子中的端1N原子(质子受体)相互作用,形成卤键结构,而复合物S4中的卤键结构的质子受体为N2O分子中的端3O原子. 经B3LYP/6-311++G**水平上的计算,考虑了基组重叠误差(BSSE)校正的单体间相互作用能在-1.56 ~ -8.73 kJ·mol-1之间. 采用自然键轨道理论(NBO)对两种单体间相互作用的本质进行了考查,并通过分子中原子理论(AIM)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.     

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     

Theoretical study of dihydrogen bonded clusters of water with tetrahydroborate

Ab initio calculations were used to analyze interactions of BH₄ ^? with 1?C4 molecules of H₂O at the MP2/6-311++G(d,p) and B3LYP/6-311++G(d,p) computational levels. The negative cooperativity for dihydrogen bond clusters containing H₂O···H₂O hydrogen bonds is more remarkable. The negative cooperativity is increased with increasing the size and also the number of hydrogen bonds in the cluster. The B?CH stretching frequencies show blue shifts with respect to cluster formation. Also greater blue shift of stretching frequencies where predicted for B?CH bonds which did not contribute in dihydrogen bonding with water molecules. The structures obtained have been analyzed with the Atoms in Molecules (AIM) methodology.     

The phonon-assisted proton transfer between hydrazinium cations in (N2H5)2HMF6·2H2O (M∈{Ga, Al, Fe}) type compounds. FT IR and quantum theoretical study

FT IR spectra of a series of compounds with a general formula (N₂H₅)₂HMF₆·2H₂O (where M∈{Ga, Al, Fe}) were recorded at variable temperatures (from ∼100 to 300 K, at 10 K intervals). The appearance of the spectral region of ν(N-N) modes due to hydrazinium cations further supports the conclusions regarding the N₂H₅⁺?H⁺?N₂H₅⁺ hydrogen bond potential well based on Raman spectroscopic data [J. Raman Spectrosc. 28 (1997) 315]. The appearance of two bands corresponding to the ν(N-N) modes in the low temperature FT IR spectra that merge into one upon heating is a clear evidence of a symmetric potential well through which a phonon-assisted proton transfer (PAPT) occurs at higher temperatures. Ab initio MP2/6-311++G(2d,p) quantum chemical study of the proton transfer potential within the N₂H₅⁺?H⁺?N₂H₅⁺ cluster confirmed its double-minimum character. The first-order saddle point found on the MP2/6-311++G(2d,p) potential energy hypersurface corresponds to a centrosymmetric structure (C_2h symmetry), with the proton placed at the inversion center. The potential energy curve along the tunnelling coordinate was calculated by the intrinsic reaction coordinate (IRC) methodology, leading to an adiabatic PT barrier height of 3.94 kcal mol⁻¹ and a tunneling rate of 1.98 s⁻¹. The corresponding MP4(SDTQ)/6-311++G(2d,p)//MP2/6-311++G(2d,p) value of the adiabatic PT barrier height is 4.26 kcal mol⁻¹.     

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.     

C3H2环丙烯基自由基与O(3P)反应机理的理论研究

本文用量子化学密度泛函方法对C3H2 (环丙烯基自由基)与O(3P)反应的机理进行了理论研究。在B3LYP/6-311++G**计算水平上优化了各驻点（过渡态,中间体,产物）的几何结构,在QCISD(T)/6-311++G**水平下计算了各物质的单点能量,在两种水平下计算了298K和600K时的能量。计算结果表明：C3H2 + O(3P) 反应可以生成P1 (C2H +HCO),P2 (C2H2 + CO) 和P3 (HC3O+H)三种产物。生成P1反应通道的能垒最低,即P1为主要产物,与实验的结果一致。产物P1可以通过路径：R→ IM1→ IM2→ P1获得。本文详细地讨论了C3H2 + O(3P) 的反应机理,并从理论上对实验结果进行了验证。研究结果有助于深入理解C3H2 + O(3P)反应机理以及C3H2在大气中的燃烧过程。     

Reaction of C2HCl2＋O2： Combined TR-FTIR Spectroscopy and Electronic Structure

The product channels and mechanisms of the C2HC12＋O2 reaction are investigated by step-scan time-resolved Fourier transform infrared emission spectroscopy and the G3MP2// B3LYP/6-311G（d,p） level of electronic structure calculations. Vibrationally excited products of HCI, CO, and CO2 are observed in the IR emission spectra and the product vibrational state distribution are determined which shows that HCI and CO are vibrationally excited with the nascent average vibrational energy estimated to be 59.8 and 51.8 kJ/mol respectively. In combination with the G3MP2//B3LYP/6-311G（d,p） calculations, the reaction mechanisms have been characterized and the energetically favorable reaction pathways have been suggested.     

Anharmonic Effect of the Unimolecular Dissociation of Ethyl Radical

The anharmonic and harmonic rate constants have been calculated for the unimolecular dissociation of ethyl radical using the method proposed by Yao and Lin (YL method) at both B3LYP/6‐311++G** and MP2/6‐311++G** levels. The different rate constants indicate that the results obtained from B3LYP and MP2 method are very close. The anharmonic and tunneling effect of the title reaction has also been examined. The comparison shows that, both in microcanonical and canonical systems, the anharmonic rate constants are higher than those for harmonic cases, especially in the case of high total energies and temperatures, which indicates that anharmonic effect of the unimolecular dissociation of ethyl into C₂H₄ and H is so significant that cannot be neglected. The tunneling effect is very small for the decomposition of C₂H₅ radical.     

Halogen switching of azacarbenes C2NH ground states at ab initio and DFT levels

Relative stabilities and singlet–triplet energy differences are calculated for 24 C₂NX azacarbenes (where X is H, F, Cl, and Br). Three skeletal arrangements are employed including azacyclopropenylidene, [(imino)methylene]carbene, and cyanocarbene. Halogens appear to alternate the electronic ground states of C₂NH azacarbenes, from triplet to singlet states, at MP3/6‐311++G**, B1LYP/6‐311++G**, B3LYP/6‐311++G**, MP2/6‐311++G**, MP4(SDTQ)/6‐311++G**, QCISD(T)/6‐311++G**, CCSD(T)/6‐311++G**, CCSD(T)/cc‐pVTZ, G1, and G2 levels of theory. The aromatic characters of singlet cyclic azacyclopropenylidenes are measured using GIAO–NICS calculations. Linear correlations are found between the B3LYP/6‐311++G** calculated LUMO–HOMO energy gaps (ΔE_{HOMO ‐ LUMO}) of the singlet carbenes versus their corresponding singlet–triplet energy separations (ΔE). Electrophilic characters are found for all singlet azacarbenes in their addition reactions to alkenes with the highest electrophilicity being exhibited for X = F. © 2008 Wiley Periodicals, Inc. Heteroatom Chem 19:377–388, 2008; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/hc.20442     

Isomerization study of C5H5NO molecules

The complex potential energy surface (PES) for the isomerization of C₅H₅NO species, including 18 isomers and 23 interconversion transition states, is probed theoretically at the B3LYP/6‐311++G(d,p) and MP2//B3LYP/6‐311++G(d,p) levels of theory. The geometries and relative energies for various stationary points were determined. The zero‐point vibrational energy (ZPVE) corrections have been made to calculate the reliable energy. We predicted a six‐membered ring structure as a global minima isomer I, which is 118.49 and 131.48 kcal · mol^?1 more stable than the least stable, four‐ and three‐membered ring isomer VIII at B3LYP and MP2//B3LYP levels of theory, respectively. The isomers and interconversion transition states have verified by frequency calculation. The intrinsic reaction coordinates (IRC) calculations have been performed to confirm that each transition state is linked by the desired reactants and products. The isomer stability has been studied using relative energies, chemical hardness, and chemical potential. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Theoretical study on the role of cooperative solvent molecules in the neutral hydrolysis of ketene

The results of a theoretical study of the reaction mechanism for the neutral hydration of ketene, H₂C=C=O + (n + 1) H₂O → CH₃COOH + nH₂O (n = 0–4), in solution are presented. All structures were optimized and characterized at the MP2(fc)/6-31 + G* level of theory, and then re-optimized by MP2(fc)/6-311 ++G**, and the effect of the bulk solvent is taken into account according to the conductor-like polarized continuum model (CPCM) using the gas MP2(fc)/6-311 ++G** geometries. Energies were refined for five-water hydration at higher level of theory, QCISD(T)(fc)/6-311 ++G**//MP2(fc)/6-311 ++G**. In the combined supermolecular/continuum model, one water molecule directly attacks the central C-atom, and the other four explicit water molecules are divided into two groups, one acting as catalyst(s) by participating in the proton transfer to reduce the tension of proton transfer ring, and the other being placed near the non-reactive oxygen or carbon atom in order to catalyze the hydration by engaging in hydrogen-bonding to the substrate (the so-called cooperative effect). Between the two possible nucleophilic addition reactions of water molecule, across the C=O bond or the C=C bond, the former one is preferred. Our calculations suggest that the favorable hydrolysis mechanism of ketene involves a sort of eight-membered ring transition structure formed by a three-water proton transfer loop, and a cooperative dimeric water near the non-reactive carbon-atom. The best-estimated in the present paper for the rate-determining barrier in solution, $ \Updelta G_{\text{sol}}^{ \ne } $ (298 K), is about 58 kJ/mol, reasonably close to the available experimental result.     

Computational study of the substitution effect on the mechanism for phospha-Wittig reaction

The phospha-Wittig reaction HP=PH₃ + O=CHX → HP=CHX + O=PH₃ (X = H, F, Cl, Me, OMe, NMe₂, CMe₃) was examined using the density functional theory calculations. All of the structures were completely optimized at the B3LYP/6-311++G** level of theory. The reactivities of various O=CHX were examined by estimating their activation energies. The main finding of this work is that the configuration mixing model can successfully predict the relative ordering of the activation energy and reaction enthalpies of the phospha-Wittig reaction. It was demonstrated that O=CHX with more electro-releasing substituents will possess a smaller singlet–triplet splitting. This will facilitate the phospha-Wittig reaction and will result in a larger exothermicity.     

Theoretical and kinetic study of the H + C2H5CN reaction

The reaction of H radical with C₂H₅CN has been studied using various quantum chemistry methods. The geometries were optimized at the B3LYP/6‐311+G(d,p) and B3LYP/6‐311++G(2d,2p) levels. The single‐point energies were calculated using G3 and BMC‐CCSD methods based on B3LYP/6‐311++G(2d,2p) geometries. Four mechanisms were investigated, namely, hydrogen abstraction, C‐addition/elimination, N‐addition/elimination and substitution. The kinetics of this reaction were studied using the transition state theory and multichannel Rice‐Ramsperger‐Kassel‐Marcus methodologies over a wide temperature range of 200–3000 K. The calculated results indicate that C‐addition/elimination channel is the most feasible over the whole temperature range. The deactivation of initial adduct C₂H₅CHN is dominant at lower temperature with bath gas H₂ of 760 Torr; whereas C₂H₅+HCN is the dominant product at higher temperature. Our calculated rate constants are in good agreement with the available experimental data. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Computational study on the mechanism and thermodynamic of atmospheric oxidation of HCN with ozone

Oxidation of hydrogen cyanide with ozone on the singlet potential energy surface is carried out using the MP2 and CCSD(T)//MP2 theoretical approaches in connection with the 6-311++G(d,p) basis set. In this reaction, energy barrier of transition states are low, so the reaction of HCN with ozone can occur easily at atmospheric condition. With variety of pre-reactive complex, five types of products are obtained at the MP2 method which four types of them have enough thermodynamic stability at the standard condition. CO₂ + HNO are final adducts in a process that is computed to be exothermic by ?132.605 kcal/mol in standard enthalpy and spontaneous by ?144.166 kcal/mol in Gibbs free energy of reaction. In kinetic viewpoint, the production of OCN + HO₂ adducts path with one low level transition state is the most favored path.     

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     

Hydrogen-Bonding Interaction in a Complex of Amino Acid with <Emphasis Type="Italic">N</Emphasis>,<Emphasis Type="Italic">N</Emphasis>-Dimethylformamide Studied by DFT Calculations

Theoretical studies on hydrogen-bonded complexes between amino acids (glycine, alanine and leucine) and N,N-dimethylformamide (DMF) in gas phase have been carried out using density functional theory (DFT) and ab initio calculations at the B3LYP/6-311++G** and MP2/6-311++G** theory levels. The structures, binding energy, stretching frequency and bond characteristics of the mentioned complexes were calculated. The NH₂ and COOH groups of amino acids form different types of hydrogen bonds with the DMF molecule, as well as alkyl side chains. High binding energy suggests multiple hydrogen bonds present in one complex. The nearly linear OH???O and NH???O contacts are stronger than a conventional hydrogen bond interaction with their H???O separation between 1.74 and 2.14 Å. The weaker CH???O H-bond is also discussed as being a crucial interaction in biological systems involving amino acids. The formation of this interaction results in a blue shift in the CH stretching frequency.     

Theoretical studies on a possible synthesis reaction pathway on N8 (CS) clusters

The potential energy surfaces of N₈ clusters were investigated by density functional theory (DFT) and a possible synthesis reaction pathway for N₈ (C_S) was suggested. The species involved were fully optimized up to the B3LYP/6‐311+G* level of theory. Relative energies were further calculated at the QCISD/6‐311+G*//B3LYP/6‐311+G* level. The reaction rate constants of these steps from the 1 (N₅⁺?N₃^?, complex, C_S) to 2 (N₈, C_S), 2 (N₈, C_S) to 3 (N₈, C_S), 3 (N₈, C_S) to 4 (N₈, D_2d), and 4 (N₈, D_2d) to 5 (N₈, C_S) reactions were predicted by the VTST theory. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1334–1339, 2001     

Synthesis, spectroscopic characterization, X-ray structure and DFT studies on 4-(2-hydroxyphenyl)-4,5,6,7-tetrahydro-1H-imidazo[4,5-c]pyridin-5-ium chloride hydrate

The title molecular salt, 4-(2-hydroxyphenyl)-4,5,6,7-tetrahydro-1H-imidazo[4,5-c]pyridin-5-ium chloride hydrate (C₁₂H₁₄N₃O⁺·Clˉ·H₂O), was synthesized and characterized by IR-NMR spectroscopy and single-crystal X-ray diffraction. In addition to the molecular geometry from X-ray experiment, the molecular geometry, vibrational frequencies and gauge-independent atomic orbital (GIAO) ¹H and ¹³C NMR chemical shift values of the title compound in the ground state have been calculated using the density functional theory (DFT/B3LYP) method with the 6-31++G(d,p) and 6-311++G(d,p) basis sets, and compared with the experimental data. Besides, molecular electrostatic potential (MEP) distribution and non-linear optical properties of the title compound were investigated by theoretical calculations at the B3LYP/6-311++G(d,p) level.     

Author index

The triplet potential energy surface for the reaction of ethylene with molecular oxygen in the ground state has been calculated at the QCISD(T, full)/6-311++G(3df,2p)//MP2(full)/6-311G(d,p) level of theory. Four intermediates and 17 transition states are located along the minimum energy reaction path. Six major radical product channels are revealed, namely, C₂H₃ + HO₂, O(³P) + C₂H₄O, H + C₂H₃O₂, OH + C₂H₃O, OH + CH₃CO, and O(³P) + CH₃CHO. In view of reaction barrier heights, the dominant channels are predicted to be C₂H₃ + HO₂, O(³P) + C₂H₄O, H + C₂H₃O₂, and C₂H₄O₂(³A″). The calculated rate constants for the abstraction reaction path are in good agreement with the available experimental data. The implication of the current results in the initiation of the combustion of ethylene is discussed.