FC(O)O自由基与NO2反应的量子化学研究

用量子化学DFT, MP2, G3和G3MP2方法对FC(O)O自由基与NO₂的反应机理进行了理论研究. 优化了反应势能面上各驻点的几何结构, 通过内禀反应坐标(IRC)计算和振动分析, 确认了反应中的过渡态, 并用过渡态理论(TST)计算了相关反应的速率常数.     

N(4S)+CH3X(X=Cl、Br)反应机理的理论研究

The reaction of N(4S)+CH3X(X=Cl、Br) was studied by the ab initio method. The geometries of the reactants, transition states and products were optimized at the MP2/6-311+G(d,p) level. The corresponding vibration frequencies were calculated at the same level. The single-point calculations for all the stationary points were carried out at the MP2/6-311++G(3df,2p) and the QCISD(T)/6-311+G(d,p) levels using the MP2/6-311+G(d,p) optimized geometries. The energies of all the stationary points were calculated by the G2MP2 method. The results of this theoretical study indicate that the reaction has three reaction channels: H abstraction reaction channel a, Cl or Br abstraction reaction channel b and substitution reaction channel c. For the N(4S)+CH3Cl reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have a slight contribution in the reaction. For the N(4S)+CH3Br reaction, reaction channel a is the main reaction channel. Reaction channels b and c may have some contribution in the reaction.     

Theoretical Studies on the Aminolysis of Ester: Effects of the Catalysis and Substituent to the Reaction of Methyl Indole-3-acetate with Ammonia

A comprehensive exploration of the aminolysis mechanism for methyl indole-3-acetate with ammonia is carried out by employing the B3 LYP/6-311++G(d,p), M06-2 X/6-311++G(d,p) and MP2/6-311++G(d,p)//M06-2 X/6-311++G(d,p) levels. Two alterative reaction channels of the concerted and addition/elimination stepwise processes including the uncatalyzed, base-catalyzed reactions are taken into consideration. Subsequently, the substituent effects and solvent effects in methanol are also evaluated at the M06-2 X/6-311++G(d,p) level. The calculated results indicate that the calculated values of M06-2 X level are quite close to those of MP2, the stepwise pathway has more advantages to the concerted one for all of the reaction processes and the catalyst facilitates the proton migration and decreases the energy barriers as well. It is shown that the most preferred mechanism is the based-catalyzed stepwise process, the substituent of NH2 group slightly accelerates all the aminolysis reaction processes, and the solvent effect does not remarkably change the mechanism of the reaction.     

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.     

CH2F自由基与HNCO反应机理的理论研究

The reaction mechanism of CH2F radical with HNCO was investigated by density functional theory (DFT)at the B3LYP/6-311++G(d,p) level. The geometries of the reactants, the intermediates, the transition states and the products were optimized. The transition states were verified through the vibration analysis.The relative energies were calculated at the QCISD(T)/6-311++G**//B3LYP/6-311++G(d,p) level. Seven feasible reaction pathways of the reaction were studied. The results indicate that the pathway (5) is the most favorable to occur, so it is the main pathway of the reaction.     

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

CH3与NO2的反应

The elementary reaction of the CH3 radical with NO2 was investigated by time-resolved FTIR spectroscopy and quantum chemical calculations. The CH3 radical was produced by laser photolysis of CH3Br or CH3I at 248 nm. Vibrationally excited products OH, HNO and CO2 were observed by the time-resolved spectroscopy for the first time. The formation of another product NO was also verified. According to these observations, the product channels leading to CH3O+NO, CH2NO+OH and HNO+H2CO were identified. The channel of CH3O+NO was the major one. The reaction mechanisms of the above channels were studied by quantum chemical calculations at CCSD(T)/6-311++G(df,p)//MP2/6-311G(d,p) level. The calculated results fit with the experimental observations well.     

CH3SS与XO (X＝F,Cl, Br)反应过程的电子密度拓扑研究

利用MP2/6-311＋＋G(d,p)//B3LYP/6-311＋＋G(d,p)对CH3SS与XO (X＝F, Cl, Br)的反应机理进行了研究. 着重从电子密度拓扑分析角度讨论了化学键的生成和断裂. 计算结果表明单线态反应为主要反应通道, 且由于该通道的反应能垒低、放热明显, 说明CH3SS与XO在大气中比较容易进行. 电子密度拓扑分析表明, 在单线态抽氢反应通道中存在着四元环状过渡结构, 随着反应进行, 此四元环状过渡结构通过一个T-型结构变为三元环状过渡结构, 最后环状结构消失得到产物.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Ab initio and experimental thermodynamic and kinetic study of proton-assisted bond activation in gaseous hydrocarbons: deconvolution of reaction efficiencies in the case of adamantane

1) Protonation at all possible sites of adamantane (C(10)H(16)) was studied at the MP2/6-311++G(3df,2p)//MP2/6-311++G(d,p) level. This provided values of the changes in the thermodynamic state functions for these processes. Whenever direct comparison was possible, the agreement with experimental data was very good. 2) By the same means, the reaction paths linking the various species obtained in these reactions were analyzed. 3) Fourier transform ion cyclotron resonance (FT-ICR) spectroscopy was used to determine the rate constants for proton transfer from 16 protonated reference bases to adamantane in the gas phase. Also, the rate constants for the formation of ionic products in these reactions were determined. 4) The experimental reaction rates were successfully predicted and refined on the basis of a simple mechanistic model based on the reaction profiles indicated above. 5) Our results hint at the potential usefulness of this approach for mechanistic studies.     

A theoretical study of the thermodynamic and hydrogen-bond basicity of TEMPO radical and related nitroxides

The use of B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) calculations on TEMPO and four related nitroxide radicals having another oxygen functionality (ketone, hydroxyl and ether) has allowed to determine the relative basicities (thermodynamic and hydrogen-bonded) of the oxygen lone pair relative to the nitroxide radical. The differences are small, especially the B3LYP ones, but in all cases they favor the radical. This is consistent with experimental results in the case of the hydroxyl group but not in the case of the keto group.     

H2CO与双卤分子间卤键的电子密度分析

运用B3LYP和MP2方法在6-311++G(d,p)基组水平上, 对H2CO-XY(XY=F2、Cl2、Br2、ClF、BrF、BrCl)卤键体系进行构型全优化, 得到了O…X—Y型卤键复合物. 结果表明, MP2/6-311++G(d,p)计算结果与实验值较吻合. 并在MP2水平下计算了分子间的相互作用能, 用完全均衡校正CP(counterpoise procedure)方法对基函数重叠误差(BSSE)进行了校正. 利用电子密度拓扑分析方法对卤键复合物的电子密度拓扑性质进行了分析研究.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

Quantum mechanical investigation of the atmospheric reaction CH3O2 + NO

The important stationary points on the potential energy surface of the reaction CH(3)O(2) + NO have been investigated using ab initio and density functional theory techniques. The optimizations were carried out at the B3LYP/6-311++G(d,p) and MP2/6-311++G(d,p) levels of theory while the energetics have been refined using the G2MP2, G3//B3LYP, and CCSD(T) methodologies. The calculations allow the proper characterization of the transition state barriers that determine the fate of the nascent association conformeric minima of methyl peroxynitrite. The main products, CH(3)O + NO(2), are formed through either rearrangement of the trans-conformer to methyl nitrate and its subsequent dissociation or via the breaking of the peroxy bond of the cis-conformer to CH(3)O + NO(2) radical pair. The important consequences of the proposed mechanism are (a) the allowance on energetic grounds for nitrate formation parallel to radical propagation under favorable external conditions and (b) the confirmation of the conformational preference of the homolytic cleavage of the peroxy bond, discussed in previous literature.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

Theoretical study and rate constant computation on the reaction HFCO + OH --> CFO + H2O

The potential energy surface, including the geometries and frequencies of the stationary points, of the reaction HFCO + OH is calculated using the MP2 method with 6-31+G(d,p) basis set, which shows that the direct hydrogen abstraction route is the most dominating channel with respect to addition and substitution channels. For the hydrogen abstraction reaction, the single-point energies are refined at the QCISD(T) method with 6-311++G(2df,2pd) basis set. The calculated standard reaction enthalpy and barrier height are -17.1 and 4.9 kcal mol(-1), respectively, at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The reaction rate constants within 250-2500 K are calculated by the improved canonical variational transition state theory (ICVT) with small-curvature tunneling (SCT) correction at the QCISD(T)/6-311++G(2df,2pd)//MP2/6-31+G(d,p) level of theory. The fitted three-parameter formula is k = 2.875 x 10(-13) (T/1000)1.85 exp(-325.0/T) cm(3) molecule(-1) s(-1). The results indicate that the calculated ICVT/SCT rate constant is in agreement with the experimental data, and the tunneling effect in the lower temperature range plays an important role in computing the reaction rate constants.     

On the flexibility of beta-peptides

The full conformational space was explored for an achiral and two chiral beta-peptide models: namely For-beta-Ala-NH2, For-beta-Abu-NH2, and For-beta-Aib-NH2. Stability and conformational properties of all three model systems were computed at different levels of theory: RHF/3-21G, B3LYP/6-311++G(d,p)//RHF/3-21G, B3LYP/6-311++G(d,p), MP2//B3LYP/6-311++G(d,p), CCSD//B3LYP/6-311++G(d,p), and CCSD(T)//B3LYP/6-311++G(d,p). In addition, ab initio E = E(phi, micro, psi) potential energy hypersurfaces of all three models were determined, and their topologies were analyzed to determine the inherent flexibility properties of these beta-peptide models. Fewer points were found and assigned than expected on the basis of Multidimensional Conformational Analysis (MDCA). Furthermore, it has been demonstrated, that the four-dimensional surface, E = E(phi, mu, psi), can be reduced into a three-dimensional one: E = E[phi, f(phi), psi]. This reduction of dimensionality of freedom of motion suggests that beta-peptides are less flexible than one would have thought. This agrees with experimental data published on the conformational properties of peptides composed of beta-amino acid residues.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

Theoretical studies on dynamics and thermochemistry of the reactions CF(3)CHCl(2)+Cl-->CF(3)CCl(2)+HCl and CF(3)CHFCl+Cl-->CF(3)CFCl+HCl

A dual-level direct dynamics study has been carried out for the two hydrogen abstraction reactions CF(3)CHCl(2)+Cl and CF(3)CHFCl+Cl. The geometries and frequencies of the stationary points are optimized at the BHLYP/6-311G(d,p), B3LYP/6-311G(d,p), and MP2/6-31G(d) levels, respectively, with single-point calculations for energy at the BHLYP/6-311++G(3df,2p), G3(MP2), and QCISD(T)/6-311G(d,p) levels. The enthalpies of formation for the species CF(3)CHCl(2), CF(3)CHFCl, CF(3)CCl(2), and CF(3)CFCl are evaluated at higher levels. With the information of the potential energy surface at BHLYP/6-311++G(3df,2p)//6-311G(d,p) level, we employ canonical variational transition-state theory with small-curvature tunneling correction to calculate the rate constants. The agreement between theoretical and experimental rate constants is good in the measured temperature range 276-382 K. The effect of fluorine substitution on reactivity of the C-H bond is discussed.