C?X⋅⋅⋅O Halogen Bonding: Interactions of Trifluoromethyl Halides with Dimethyl Ether

The formation of weakly bound molecular complexes between dimethyl ether (DME) and the trifluoromethyl halides CF₃Cl, CF₃Br and CF₃I dissolved in liquid argon and in liquid krypton is investigated, using Raman and FTIR spectroscopy. For all halides evidence is found for the formation of C? X???O halogen‐bonded 1:1 complexes. At higher concentrations of CF₃Br, a weak absorption due to a 1:2 complex is also observed. Using spectra recorded at temperatures between 87 and 125 K, the complexation enthalpies for the complexes are determined to be ?6.8(3) kJ mol^?1 (DME?CF₃Cl), ?10.2(1) kJ mol^?1 (DME?CF₃Br), ?15.5(1) kJ mol^?1 (DME?CF₃I), and ?17.8(5) kJ mol^?1 [DME(?CF₃Br)₂]. Structural and spectral information on the complexes is obtained from ab initio calculations at the MP2/ 6‐311++G(d,p) and MP2/6‐311++G(d,p)+LanL2DZ* levels. By applying Monte Carlo free energy perturbation calculations to account for the solvent influences, and statistical thermodynamics to estimate the zero‐point vibrational and thermal influences, the ab initio complexation energies are converted into complexation enthalpies for the solutions in liquid argon. The resulting values are compared with the experimental data deduced from the cryosolutions.     

Reaction paths of BCl3?+?CH4?+?H2 in the chemical vapor deposition process

The reaction paths in the chemical vapor deposition preparation of boron carbides with BCl₃?CCH₄?CH₂ precursors were investigated theoretically in detail with a total number of 82 intermediates (IM) and 118 transition states (TS). The geometries of the species were optimized with B3PW91/6-311G(d,p) method and the TS as well as their linked IM were confirmed with the frequency and the intrinsic reaction coordinates analyses at the same theoretical level. The energy barriers and the reaction energies were determined with the accurate model chemistry method G3(MP2) after a diagnosis of the non-dynamic electronic correlations. The heat capacities and entropies were obtained with statistical thermodynamics. The Gibbs free energies at 298.15?K for all of the reaction steps were reported and the data at any temperature can be developed with the classical thermodynamics by using the fitted (as a function of temperature) heat capacities. All the possible elementary reactions, including both direct decomposition and the radical attacking dissociations for each reaction step were examined. It was found that there are nine reaction steps in the lowest reaction pathway to produce the final boron carbide and five steps to produce boron. The highest energy barrier in the lowest reaction pathway is 238.6?kJ?mol^?1 at 298.15?K and 346.0?kJ?mol^?1 at 1,200?K for producing BC, and is 294.7?kJ?mol^?1 at 298.15?K and 314.2?kJ?mol^?1 at 1,200?K for producing B.     

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.     

Cooperative and Diminutive Interplay Between Lithium and Dihydrogen Bonding in F3YLi…︁NCH…︁HMH and F3YLi…︁HMH…︁HCN Triads (Y=C,Si; M=Be,Mg)

The F₃YLi…NCH…HMH and F₃YLi…HMH…HCN triads (Y=C, Si; M=Be, and Mg) are connected by lithium and dihydrogen bonds. To understand the properties of the systems better, the corresponding dyads are also studied. Molecular geometries, binding energies, infrared spectra and NMR properties of monomers, dyads, and triads are investigated at the MP2/6‐311++G** computational level. Particular attention is paid to parameters, such as cooperative energies, and many‐body interaction energies. Triads with the HMH molecule located at the end of the chain, show energetic cooperativity ranging between ?3.66 to ?7.59 kJ mol^‐1. When the HMH molecule is located in the middle, the obtained cluster is diminutive with an energetic effect between 3.49 to 5.17 kJ mol^‐1. The electronic properties of the complexes are analyzed using parameters derived from the atoms in molecules (AIM) methodology.     

Cooperativity between the Halogen Bond and the Hydrogen Bond in H3N⋅⋅⋅XY⋅⋅⋅HF Complexes (X,Y=F,Cl, Br)

Ab initio calculations are used to provide information on H₃N???XY???HF triads (X, Y=F, Cl, Br) each having a halogen bond and a hydrogen bond. The investigated triads include H₃N???Br₂‐HF, H₃N???Cl₂???HF, H₃N???BrCI???HF, H₃N???BrF???HF, and H₃N???ClF???HF. To understand the properties of the systems better, the corresponding dyads are also investigated. Molecular geometries, binding energies, and infrared spectra of monomers, dyads, and triads are studied at the MP2 level of theory with the 6‐311++G(d,p) basis set. Because the primary aim of this study is to examine cooperative effects, particular attention is given to parameters such as cooperative energies, many‐body interaction energies, and cooperativity factors. The cooperative energy ranges from ?1.45 to ?4.64 kcal mol^?1, the three‐body interaction energy from ?2.17 to ?6.71 kcal mol^?1, and the cooperativity factor from 1.27 to 4.35. These results indicate significant cooperativity between the halogen and hydrogen bonds in these complexes. This cooperativity is much greater than that between hydrogen bonds. The effect of a halogen bond on a hydrogen bond is more pronounced than that of a hydrogen bond on a halogen bond.     

Ab initio Mechanism Study on the Reaction of Chlorine Atom with Formic Acid

The potential energy surface(PES) for the reaction of Cl atom with HCOOH is predicted using ab initio molecular orbital calculation methods at UQCIDS(T,full)6-311 G(3df,2p)//UMP2(full)/6-311 G(d,P) level of theory with zero-point vibrational energy (ZPVE) correction.The calculated results show that the reaction mechanism of Cl atom with formic acid is a C-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom with a 3.73kJ/mol reaction barrier height,leading to the formation of cis-HOCO radical which will reacts with Cl atom or other molecules in such a reaction system.Because the reaction barrier height of O-site hydrogen abstraction reaction from cis-HOC(H)O molecule by Cl atom which leads to the formation of HCO2 radical is 67.95kJ/mol,it is a secondary reaction channel in experiment,This is in good agreement with the prediction based on the previous experiments.     

MP2, DFT‐D,and PCM study of the HMB–TCNE complex: Thermodynamics,electric properties,and solvent effects

Geometry, thermodynamic, and electric properties of the π‐EDA complex between hexamethylbenzene (HMB) and tetracyanoethylene (TCNE) are investigated at the MP2/6‐31G* and, partly, DFT‐D/6‐31G* levels. Solvent effects on the properties are evaluated using the PCM model. Fully optimized HMB–TCNE geometry in gas phase is a stacking complex with an interplanar distance 2.87 × 10^?10 m and the corresponding BSSE corrected interaction energy is ?51.3 kJ mol^?1. As expected, the interplanar distance is much shorter in comparison with HF and DFT results. However the crystal structures of both (HMB)₂–TCNE and HMB–TCNE complexes have interplanar distances somewhat larger (3.18 and 3.28 × 10^?10 m, respectively) than our MP2 gas phase value. Our estimate of the distance in CCl₄ on the basis of PCM solvent effect study is also larger (3.06–3.16 × 10^?10 m). The calculated enthalpy, entropy, Gibbs energy, and equilibrium constant of HMB–TCNE complex formation in gas phase are: ΔH⁰ = ?61.59 kJ mol^?1, ΔS = ?143 J mol^?1 K^?1, ΔG⁰ = ?18.97 kJ mol^?1, and K = 2,100 dm³ mol^?1. Experimental data, however, measured in CCl₄ are significantly lower: ΔH⁰ = ?34 kJ mol^?1, ΔS = ?70.4 J mol^?1 K^?1, ΔG⁰ = ?13.01 kJ mol^?1, and K = 190 dm³ mol^?1. The differences are caused by solvation effects which stabilize more the isolated components than the complex. The total solvent destabilization of Gibbs energy of the complex relatively to that of components is equal to 5.9 kJ mol^?1 which is very close to our PCM value 6.5 kJ mol^?1. MP2/6‐31G* dipole moment and polarizabilities are in reasonable agreement with experiment (3.56 D versus 2.8 D for dipole moment). The difference here is due to solvent effect which enlarges interplanar distance and thus decreases dipole moment value. The MP2/6‐31G* study supplemented by DFT‐D parameterization for enthalpy calculation, and by the PCM approach to include solvent effect seems to be proper tools to elucidate the properties of π‐EDA complexes. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Ab initio study on the mechanism of reaction HNCO+NH2

Ab initio UMP2 and UQCISD(T) calculations, with 6-311G** basis sets, were performed for the titled reactions. The results show that the reactions have two product channels: NH₂+ HNCO?NH₃+NCO (1) and NH₂+HNCO?N₂H₃+CO (2), where reaction (1) is a hydrogen abstraction reaction via an H-bonded complex (HBC), lowering the energy by 32.48 kJ/mol relative to reactants. The calculated QCISD(T)//MP2(full) energy barrier is 29.04 kJ/mol, which is in excellent accordance with the experimental value of 29.09 kJ/mol. In the range of reaction temperature 2300–2700 K, transition theory rate constant for reaction (1) is 1.68×10¹¹–3.29×10¹¹ mL·mol^-1·s^-1, which is close to the experimental one of 5.0×10¹¹mL·mol^-1·s^-1or less. However, reaction (2) is a stepwise reaction proceeding via two orientation modes,cis andtrans, and the energy barriers for the rate-control step at our best calculations are 92.79 kJ/mol (forcis-mode) and 147.43 kJ/mol (fortrans-mode), respectively, which is much higher than reaction (1). So reaction (1) is the main channel for the titled reaction.     

Computational study on the mechanism for the gas‐phase reaction of dimethyl disulfide with OH

The mechanisms for the reaction of CH₃SSCH₃ with OH radical are investigated at the QCISD(T)/6‐311++G(d,p)//B3LYP/6‐311++G(d,p) level of theory. Five channels have been obtained and six transition state structures have been located for the title reaction. The initial association between CH₃SSCH₃ and OH, which forms two low‐energy adducts named as CH₃S(OH)SCH₃ (IM1 and IM2), is confirmed to be a barrierless process, The S? S bond rupture and H? S bond formation of IM1 lead to the products P1(CH₃SH + CH₃SO) with a barrier height of 40.00 kJ mol^?1. The reaction energy of Path 1 is ?74.04 kJ mol^?1. P1 is the most abundant in view of both thermodynamics and dynamics. In addition, IMs can lead to the products P2 (CH₃S + CH₃SOH), P3 (H₂O + CH₂S + CH₃S), P4 (CH₃ + CH₃SSOH), and P5 (CH₄ + CH₃SSO) by addition‐elimination or hydrogen abstraction mechanism. All products are thermodynamically favorable except for P4 (CH₃ + CH₃SSOH). The reaction energies of Path 2, Path 3, Path 4, and Path 5 are ?28.42, ?46.90, 28.03, and ?89.47 kJ mol^?1, respectively. Path 5 is the least favorable channel despite its largest exothermicity (?89.47 kJ mol^?1) because this process must undergo two barriers of TS5 (109.0 kJ mol^?1) and TS6 (25.49 kJ mol^?1). Hopefully, the results presented in this study may provide helpful information on deep insight into the reaction mechanism. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Structure of the CCl3CCl=NCH2C6H5 molecule and reorientations of its CCl3 group, according to ab initio calculations

The optimal geometry of the isolated CCl₃CCl=NCH₂C₆H₅ molecule and the intramolecular barrier to reorientations of its trichloromethyl group are calculated by the RHF/6-31G* and B3LYP/6-31G* methods. The barrier found (14.1 kJ mol^?1) is compared to that determined previously by ³⁵Cl NQR for a crystal of this compound, which allows estimation of the contribution of intermolecular interactions to braking of this motion of the CCl₃ group. The structural features of the molecule of this compound are consistently manifested in quantum-chemical calculations and NQR spectra.     

亚烷基锡烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet alkylidenestannylene and ethylene has been investigated with MP2/3-21 G^* and B3LYP/3-21 G* methods, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD(T)//MP2/3-2 IG^* and CCSD(T)//B3LYP/3-21G^* methods, respectively. The results show that the dominant reaction pathway of the cycloaddition is that an intermediate (INT) is firstly formed between the two reactants through a barrier-free exothermic reaction of 39.7 kJ/mol, and the intermediate then isomerizes to a four-membered ring product (P2.1) via a transition state TS2.1 with a barrier of 66.8 kJ/mol.     

F Cl_2→ClF Cl和Cl′F Cl→Cl′ ClF的反应机理

用密度泛函理论(DFT)B3LYP方法,取6-311G基组,计算研究了F＋Cl2→ClF＋Cl的反应机理.求得1个线形和2个三角形过渡态,反应能垒分别为1.24、46.37和105.09kJ·mol－1;同时发现F以∠FClCl为10～20°（或120～160°）进攻Cl2时,反应无能垒.此外,求得对称反应Cl′F＋Cl→Cl′＋ClF的能垒为40.57kJ·mol－1的1个过渡态.     

Fluorine Substitution Effects on Flexibility and Tunneling Pathways: The Rotational Spectrum of 2‐Fluorobenzylamine

The effect of ring fluorination on the structural and dynamical properties of the flexible model molecule 2‐fluorobenzylamine has been studied by rotational spectroscopy in free‐jet expansion and quantum chemical methods. The complete potential energy surface originating from the flexibility of the aminic side chain has been calculated at the B3LYP/6‐311++G** level of theory and the stable geometries were also characterized with MP2/6‐311++G**. The rotational spectra show the presence of two of the predicted four stable conformers: the global minimum (I), in which the side chain’s dihedral angle with the phenyl plane is almost perpendicular, is stabilized by an intramolecular hydrogen bond between the fluorine atom and one hydrogen of the aminic group; and a second conformer II (E_II?E_I≈5 kJ mol^?1) in which the dihedral angle is smaller and the amino group points towards the aromatic ortho hydrogen atom. This conformation is characterized by a tunneling motion between two equivalent positions of the amino group with respect to the phenyl plane, which splits the rotational transition. The ortho fluorination increases, with respect to benzylamine, the tunneling splitting of this motion by four orders of magnitude. The motion is analyzed with a one‐dimensional flexible model, which allows estimation of the energy barrier for the transition state as approximately 8.0 kJ mol^?1.     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Theoretical Gas Phase Study of the Gauche and Trans Conformers of 1-Bromo-2-Iodoethane and Solvent Effects

This is a gas-phase study of the gauche and trans conformers of 1-bromo-2-iodoethane. The methods used are the second-order Møller-Plesset theory (MP2) and density functional theory (DFT). The functional used for the DFT method is B3LYP and the basis sets used are 6-311++G(d,p) for all atoms except that different basis sets, namely 3-21G, LANECP, CRENBL ECP, Stuttgart RLC ECP and 6-311G(d,p), have been explored for the iodine atom. The results indicate that the trans conformer is preferred. The energy difference between the gauche and trans conformers (ΔE _g?t) and related thermodynamic parameters are reported. The ΔE _g?t values are 12.50 kJ?mol^?1 (B3LYP) and 10.00 kJ?mol^?1 (MP2) with the basis sets being 6-311++G(d,p)[C,H,Br]/6-311G(d,p)[I]. The conformers of 1-bromo-2-iodoethane have also been subjected to vibrational analysis. The results from the two theoretical levels are in good agreement but they are not much affected by the basis set of the iodine atom. The study has been extended to explore solvent effects using Self-Consistent Reaction Field methods. The structural parameters of the conformers are little affected by the polarity of the solvent but ΔE _g?t decreases and the solvation Gibbs energy increases with increasing polarity of the solvent.     

A quantum mechanical approach to the kinetics of the hydrogen abstraction reaction H2O2 + •OH → HO2 + H2O

The kinetics of the hydrogen abstraction from H₂O₂ by ^?OH has been modeled with MP2/6‐31G*//MP2/6‐31G*, MP2‐SAC//MP2/6‐31G*, MP2/6‐31+G**//MP2/6‐31+G**, MP2‐SAC// MP2/6‐31+G**, MP4(SDTQ)/6‐311G**//MP2/6‐31G*, CCSD(T)/6‐31G*//CCSD(T)/6‐31G*, CCSD(T)/6‐31G**//CCSD(T)/6‐31G**, CCSD(T)/6‐311++G**//MP2/6‐31G* in the gas phase. MD simulations have been used to generate initial geometries for the stationary points along the potential energy surface for hydrogen abstraction from H₂O₂. The effective fragment potential (EFP) has been used to optimize the relevant structures in solution. Furthermore, the IEFPCM model has been used for the supermolecules generated via MD calculations. IEFPCM/MP2/6‐31G* and IEFPCM/CCSD(T)/6‐31G* calculations have also been performed for structures without explicit water molecules. Experimentally, the rate constant for hydrogen abstraction by ^?OH drops from 1.75 × 10^?12 cm³ molecule^?1 s^?1 in the gas phase to 4.48 × 10^?14 cm³ molecule^?1 s^?1 in solution. The same trend has been reproduced best with MP4 (SDTQ)/6‐311G**//MP2/6‐31G* in the gas phase (0.415 × 10^?12 cm³ molecule^?1 s^?1) and with EFP (UHF/6‐31G*) in solution (3.23 × 10^?14 cm³ molecule^?1 s^?1). © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 502–514, 2005     

Study on Pair Correlation Energies of Isoelectronic Systems F^－, HF and H2F^＋

The intrapair and interpair correlation energies of F-, HF and H2F^ systems are calculated and analyzed using MP2-OPT2 method of MELD program with cc-PVSZ^* basis set. From the analysis of pair correlation energies of these isoelectronlc sysoterns, it is found that the 1sF^2 pair correlation energy is trans-ferable in these three isociectronic systems. According to the definition of pair correlation contribution of one electron pair to a system, the pair correlation contribution values of these three systems are calculated. The correlation contribution values of inner electron pairs and H—F bonding electron pair in HF molecule with those in H2F^ system are compared. The results indicate that the bonding effect of a molecule is one of the im-portant factors to influence electron correlation energy of the system. The comparison of correlation energy contributions in-cluding triple and quadruple excitations with those only includ-ing singles and doubles calculated with 6-311 G(d) basis set shows that the higher.excitation correlation energy contribution gives more than 2 % of the total correlation energy for these sys-tems.     

A theoretical study on the mechanism of the addition reaction between carbene and epoxyethane

The mechanism of addition reaction between carbene and epoxyethane has been investigated employing the MP2 and B3LYP/6-311+G* levels of theory. Geometry optimization, vibrational analysis, and energy property for the involved stationary points on the potential energy surface have been calculated. Based on the calculated results at the MP2/6-311+G* level of theory, it can be predicted that there are two reaction mechanisms (1) and (2). In the first reaction carbene attacks the atom O of epoxyethane to form an intermediate 1a (IM1a), which is a barrier-free exothermic reaction. Then, IM1a can isomerize to IM1b via a transition state 1a (TS1a), where the potential barrier is 48.6 kJ/mol. Subsequently, IM1b isomerizes to a product epoxypropane (Pro1) via TS1b with a potential barrier of 14.2 kJ/mol. In the second carbene attacks the atom C of epoxyethane firstly to form IM2 via TS2a. Then IM2 isomerizes to a product allyl alcohol (Pro2) via TS2b with a potential barrier of 101.6 kJ/mol. Correspondingly, the reaction energies for the reactions (1) and (2) are −448.4 and −501.6 kJ/mol, respectively. Additionally, the orbital interactions are also discussed for the leading intermediate. The results based on the B3LYP/6-311+G* level of theory are paralleled to those on the MP2/6-311+G* level of theory. Furthermore, the halogen and methyl substituent effects of H₂C: on the two reaction mechanisms have been investigated. The calculated results indicate that the introductions of halogen or methyl make the addition reaction difficult to proceed.     

Deamination of the Radical Cation of the Base Moiety of 2′‐Deoxycytidine: A Theoretical Study

Five pathways leading to the deamination of cytosine (to uracil) after formation of its deprotonated radical cation are investigated in the gas phase, at the UB3LYP/6‐311G(d,p) level of theory, and in bulk aqueous solvent. The most favorable pathway involves hydrogen‐atom transfer from a water molecule to the N3 nitrogen of the deprotonated radical cation, followed by addition of the resulting hydroxyl radical to the C4 carbon of the cytosine derivative. Following protonation of the amino group (N4), the C4? N4 bond is broken with elimination of the NH₃^?+ radical and formation of a protonated uracil. The rate‐determining step of this mechanism is hydrogen‐atom transfer from a water molecule to the cytosine derivative. The associated free energy barrier is 70.2 kJ mol^?1.     

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)分析了复合物中氢键和卤键键鞍点处的电子密度拓扑性质.