Theoretical Studies on Reaction Mechanisms of HNCS with NH （ X^3∑）

The reaction mechanisms of HNCS with NH（X^3∑ ） were theoretically investigated. The minimum energy paths （MEP） of the reaction were calculated by using the density functional theory（DFT） at the B3LYP/6-311 ＋ ＋ G^＊＊ level. The equilibrium structural parameters, the harmonic vibrational frequencies, the total energies, and the zeropoint energies（ZPE） of all the species were calculated. The single-point energies along the MEP were further refined at the QCISD（T）/6-311 ＋ ＋ G^＊ ＊ level. It was found that the mechanisms of the HNCS ＋ NH（X^3∑） reaction involve two channels producing the HNC ＋ HNS and the N2H2 ＋ CS products. Channel 1 plays a dominant role and the HNC ＋ HNS are the main preduets. The reaction is exothermie.     

Theoretical study on the reaction of VS^＋（^3∑^-, 1^Г） with CO

Two possible reaction mechanisms of VS^＋（^3∑^-, 1^Г） with CO in the gas phase have been studied by using B3LYP/TZVP and CCSD（T）/6-311＋G （3df, 3pd） methods： the O/S exchange reaction （VS^＋＋CO→VO^＋＋CS） and the S-transfer reaction （VS^＋ ＋ CO → V^＋ ＋ COS）. The two reactions proceed via two-step and one-step mechanism, respectively. The barriers of the triplet and singlet PESs are 30.6 and 50.9 kcal/mol, respectively, for O/S exchange reaction and 7.3 and 50.2 kcal/mol, respectively, for the S-transfer reaction. The results indicate that the triplet ground state reaction is more favorable, and the S-transfer reaction is more favorable than the O/S exchange reaction, which is in good agreement with the experimental observation.     

Curing kinetics and thermal stability of epoxy blends containing phosphorous-oxirane with aromatic amide-amine as curing agents

This article describes the synthesis of a series of aromatic amide-amines and their potential use as epoxy hardeners. These amines were synthesized by the reaction of L-phenylalanine（PA） with diamines of different structures i.e.1,4- phenylene diamine（PD）,1,5-diamino naphthalene（N）,4,4’-（9-fluorenyllidene）-dianiline（F）,4,4’-diaminodiphenyl sulphide （DS） and 3,4’-oxydianiline（O） in a stoichiometric ratio（1:1）.Structural characterization of synthesized amide-amines was done with the help of elemental analysis and spectroscopic techniques viz.FT-IR,¹H-NMR and ¹³C-NMR.An epoxy blend was prepared by mixing tris（glycidyloxy） phosphine oxide（TGPO） with conventional epoxy i.e.diglycidyl ether of bisphenol-A（DGEBA） in an equivalent ratio of 2:3 to incorporate phosphorous into the main chain.The curing kinetics of the epoxy blend with synthesized aromatic amide-amines was investigated by non-isothermal DSC technique using multiple heating rate method（5,10,15 and 20 K/min.）.The activation energies were determined by fitting the experimental data into Kissinger and Ozawa kinetic models.The activation energies obtained through Ozawa method were slightly higher than those of Kissinger method but were comparable.However,both the energies were found to be dependent on the structure of amines.The thermal stability and weight loss behavior of isothermally cured thermosets were also investigated using thermogravimetric analysis（TGA） in nitrogen atmosphere.All the samples showed improved thermal stability in terms of char yield than using only amines as hardeners.     

Theoretical Studies on the Reaction Mechanism of 1-Chloroethane with Hydroxyl Radical

The reaction mechanism of 1-chloroethane with hydroxyl radical has been inves- tigated by using density functional theory (DFT) B3LYP/6-31G (d, p) method. All bond dissociation enthalpies were computed at the same theoretical level. It was found that hydrogen abstraction pathway is the most favorable. There are two hydrogen abstraction pathways with activation barriers of 0.630 and 4.988 kJ/mol, respectively, while chlorine abstraction pathway was not found. It was observed that activation energies have a more reasonable correlation with the reaction enthalpy changes (△Hr) than with bond dissociation enthalpies (BDE).     

Theoretical Study on the Reaction of Thymine with Iodomehane in the Water Phase

The reaction mechanisms of carcinogenic methylating agent iodomethane (MeI) with keto and enol tautomers of thymine (K- and E-thymine) were studied by using the B3LYP/6-311+G (d, p) method in water phase. The solvent effects were examined using the polarizable continuum model (PCM). Specifically, PCM single-point calculations at the same level of theory were performed in acetone and CCl4 that represent a range in nonpolarity. The calculated results show that the reaction of K-thymine with MeI is a two-step mechanism, whereas that of E-thymine is a one-step mechanism. Our calculations reveal that K-thymine is appreciably more stable than the enol form in the water phase or in the two solvents. The K- and E-form reaction barriers are 135.6 and 222.1 kJ/mol, respectively in water phase. These findings indicate that the reactions mentioned above could not occur efficiently in biological media in the absence of catalyst. Our conclusions are in agreement with the previous studies on the reactions of guanine with methyl chloride and methyl bromide.     

A New Target for Synthesis: Theoretical Prediction for the Reaction between Boron Nitride Nanotube and Dichlorocarbene

Understanding the chemistry of BNNT is a crucial step toward their ultimate practical use. A comparative study of Reactions A （ASWCNT （5,5） and CCl2） and B （ASWBNNT （5,5） and CCl2） have been performed by using ONIOM （B3LYP/6-31G＊： AM1） method in Gaussian03 program package. The results show that （1） the two reactions are both exothermic; （2） the mechanism of Reaction B is a two-step mechanism; （3） the difference in energy barriers suggests that the reaction of CCl2 with BNNT is easier than with CNT; （4） in reaction B, CCl2 prefers to attack the boron atom of BNNT first.     

Singlet H-cCNC-C：A Potential Interstellar Molecule

A new planar isomer of HNC3 system,H-cCNC-C,is theoretically predicted by means of B3LYP and CCSD(T) methods.The suggested species can isomerize into other five kinetically more stable isomers,which have been experimentally identified,with relatively higher reaction barriers.In view of its higher kinetic stability,we can reasonably believe that the obtained species H-cCNC-C can be experimentally observed in future studies.     

A DFT Study on the Stable Structures and Dissociation Mechanism of C3O6 Cluster

Using geometrical optimization and DFT method at the B3LYP/6-31G (d) level, nineteen equilibrium geometries were identified, and three transition states of dissociation reaction of C3O6 clusters were also found. The vibrational frequencies and intrinsic reaction coordinate (IRC) verification at the same level were computed to verify the transitions. And then we calculated the dissociation energies and analyzed the dissociation channels. The computational results show that the dissociation energies of C3O6 isomers relative to three CO2 are between 1.509 × 103 and 10.61 × 103 kJ·kg-1, and the energy barriers of the reactions are 92.857, 131.138 and 185.793 kJ·mol-1. Both the high dissociation energies and high energy barriers show that C3O6 clusters studied in this paper are stable enough to be used as high-energy-density materials.     

Mechanisms and Kinetics of Radical Reaction of O（^1D,^3P） ＋ HCN System

The reaction of HCN with O（^1D, ^3p） radical has been investigated by density functional theory （DFT） and ab initio methods. The stationary points on the reaction paths （reactants, intermediates and products） were optimized at the （U）B3LYP/aug-cc-pVTZ level. Single-point calculations were performed at the （U）QCISD（T）/aug-cc-pVTZ level for the optimized structures and all the total energies were corrected by zero-point energy. It is shown that there exist three competing mechanisms of oxygen attacking nitrogen O→N, oxygen attacking carbon O→C and oxygen attacking hydrogen O→H. The rate constants were obtained via Eyring transition-state theory in the temperature range of 600～2000 K. The linear relationship between lnk and 1/T was presented. The results show that path 1 is the main reaction channel and the product of NCO ＋ H is predominant.     

二氯亚甲基锗烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet dichloromethylene germylene and ethylene has been investigated with B3LYP/6-31G＊ method, 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）//B3LYP/6-31G＊ method. On the basis of the surface energy profile obtained with CCSD（T）// B3LYP/6-31G＊ method for the cycloaddition reaction between singlet dichloromethylene germylene and ethylene, it can be predicted that the dominant reaction pathway is that an intermediate INT1 is firstly formed between the two reactants through a barrier-free exothermic reaction of 61.7 kJ/mol, and the intermediate INT1 then isomerizes to an active four-membered ring product P2.1 via a transition state TS2, an intermediate INT2 and a transition state TS2.1, in which energy barriers are 57.7 and 42.2 kJ/mol, respectively.     

含有一个非平面杂环胺配体的新型反铂抗癌药物的水解机理

用B3LYP杂化泛函和等电子聚焦极化连续模型(IEF-PCM)研究了trans-[PtCl2(NH3)(Am)](Am: 非平面哌啶或哌嗪)新型反铂抗癌药物的水解反应机理. 对经由一般的SN2机理的第一步和第二步水解反应势能面上的稳定点进行了全优化和表征. 在水解中, 最显著的结构变化发生在反应过渡态和中间体的五配位三角双锥的赤道面上. 与经典顺铂(cisplatin)比较, 反式[PtCl2(NH3)(piperazine)]的第一步和第二步水解活化能均低于顺铂, 而反式[PtCl2(NH3)(piperidine)]的第一步水解活化能稍高于顺铂, 第二步水解活化能稍低于顺铂. 计算表明, 这些含有非平面杂环胺反铂的配合物减小了赤道面上的立体效应和水解势垒.     

Hydration process as an activation of trans- and cisplatin complexes in anticancer treatment. DFT and ab initio computational study of thermodynamic and kinetic parameters

The thermodynamic and kinetic aspects of hydration reactions of cis-/transplatin were explored. The polarizable continuum model was used for estimation of solvent effects. Using the B3LYP/6-31+G(d) method, the structures were optimized and vibrational frequencies estimated. Interaction energies and activation barriers were determined at the CCSD(T)/6-31++G(d,p) level within the COSMO approach. An associative mechanism was assumed with a trigonal-bipyramidal structure of the transition state. Within the applied model, all the hydration reactions are slightly endothermic. The Gibbs energies of cisplatin hydration amount to 7.0 and 14.2 kcal/mol for the chloride and ammonium replacement, respectively. Analogous values for the transplatin reactions are 6.8 and 11.9 kcal/mol. The determined rate constants are by several (three to four) orders of magnitude larger for the dechlorination process than for deammination. The cisplatin dechlorination rate constant was established as 1.3 x 10(-4) s(-1) in excellent accord with the experiment.     

鸟嘌呤受羟基自由基损伤反应机理的量子化学研究

用密度泛函理论(DFT)研究羟基自由基与鸟嘌呤分子加成反应的过渡态, 并进行内禀反应坐标(IRC)反应路径解析, 结果表明, 羟基自由基加成到鸟嘌呤碳碳双键上. 利用B3LYP/6-31++G**对反应物、反应物络合物、过渡态以及产物络合物等反应通道上各个能量驻点的能量进行了计算, 得到反应活化能E_a=28.0867 kJ/mol. AIM计算结果显示, 过渡态结构中鸟嘌呤分子碳碳双键结构被削弱, 羟基自由基氧原子与鸟嘌呤分子碳碳双键中的C4原子具有较强的相互作用, 双键中剩余的π电子离域到了环体系中.     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Comparison of Some Computational Methods for Geometry Optimisation of Phosphorus Acid Derivatives

Several economical methods for geometry optimization, that should be applicable to larger molecules, have been evaluated for 19 phosphorus acid derivatives. MP2/cc-pVDZ geometry optimizations are used as reference points and the geometries obtained from the other methods are evaluated with respect to deviations in bond lengths and angles, from the reference geometries. The geometry optimization methods are also compared to the much used B3LYP/6-31G(d) method. Single point energies obtained by subsequent EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) calculations on the respective equilibrium geometries are also reported relative to the energies obtained from the reference geometries. The geometries from HF/MIDI! optimizations were closer to those of the references than the geometries of the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! optimizations. The EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) single point energies obtained from the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! geometries gave a mean absolute deviation (MAD) from that of the reference geometries of 1.4-3.9 kcal mol m 1 . The HF/MIDI! geometries, however, gave EDF1/6-31+G(d) and B3LYP/6-31+G(d,p) energies with a MAD of only about 0.5 and 0.55 kcal mol m 1 respectively from the energies obtained with the reference geometries. Thus, use of HF/MIDI! for geometry optimization of phosphorus acids is a method that gives geometries of near-MP2 quality, resulting in a fair accuracy of energies in subsequent single point calculations, at a much lower computational cost other methods that give similar accuracies.     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.     

过氧烷基自由基分子内氢迁移反应类速率常数的计算

过氧烷基自由基分子内氢迁移是低温燃烧反应中的一类重要基元反应. 本文用等键反应方法计算了该类反应的动力学参数. 所有反应物、过渡态、产物的几何结构均在B3LYP/6-311+G(d,p)水平下优化得到. 本文提出了用过渡态反应中心几何结构守恒作为反应类判据, 并将该分子内氢迁移反应分为四类, 包括(1,3)、(1,4)、(1,5)、(1,n) (n=6, 7, 8)氢迁移类. 分别将这4 类反应类中最小反应体系作为类反应的主反应, 并分别在B3LYP/6-311+G(d,p)低水平和CBS-QB3 高水平下得到其近似能垒和精确能垒. 其余氢迁移反应作为目标反应, 在B3LYP/6-311+G(d,p)低水下计算得到其近似能垒, 再采用等键反应方法校正得到目标反应的精确反应势垒和精确速率常数. 研究表明, 采用等键反应方法只需在低水平用从头算计算就可以得到大分子反应体系的高精度能垒和速率常数值, 且本文按等键反应本质的分类方法更能揭示反应类的本质, 并对反应类的定义给出了客观标准. 本文的研究为碳氢化合物低温燃烧模拟中重要的过氧烷基分子内氢迁移反应提供了准确的动力学参数.     

Addition reaction of adamantylideneadamantane with Br2 and 2Br2: a computational study

Ab initio calculations were carried out for the reaction of adamantylideneadamantane (Ad=Ad) with Br2 and 2Br2. Geometries of the reactants, transition states, intermediates, and products were optimized at HF and B3LYP levels of theory using the 6-31G(d) basis set. Energies were also obtained using single point calculations at the MP2/6-31G(d)//HF/6-31G(d), MP2/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31+G(d)//B3LYP/6-31G(d) levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Only one pathway was found for the reaction of Ad=Ad with one Br2 producing a bromonium/bromide ion pair. Three mechanisms for the reaction of Ad=Ad with 2Br2 were found, leading to three different structural forms of the bromonium/Br3- ion pair. Activation energies, free energies, and enthalpies of activation along with the relative stability of products for each reaction pathway were calculated. The reaction of Ad=Ad with 2Br2 was strongly favored over the reaction with only one Br2. According to B3LYP/6-31G(d) and single point calculations at MP2, the most stable bromonium/Br3- ion pair would form spontaneously. The most stable of the three bromonium/Br3- ion pairs has a structure very similar to the observed X-ray structure. Free energies of activation and relative stabilities of reactants and products in CCl4 and CH2ClCH2Cl were also calculated with PCM using the united atom (UA0) cavity model and, in general, results similar to the gas phase were obtained. An optimized structure for the trans-1,2-dibromo product was also found at all levels of theory both in gas phase and in solution, but no transition state leading to the trans-1,2-dibromo product was obtained.     

Probing furanose ring conformation by gas-phase computational methods: energy profile and structural parameters in methyl beta-D-arabinofuranoside as a function of ring conformation

The potential energy surface of methyl beta-D-arabinofuranoside (3) has been studied by ab initio molecular orbital (HF/6-31G) and density functional theory (B3LYP/6-31G) calculations via minimization of the 10 possible envelope conformers. The partial potential energy surface identified that the global minimum and lowest energy northern conformer was E(2). In the HF calculations, (2)E was the most stable southern conformer, while the density functional theory methods identified (4)E as the local minimum in this hemisphere. Additional calculations at higher levels of theory showed that the B3LYP-derived energies of many of the envelope conformers of 3 are dependent upon the basis set used. It has also been demonstrated that B3LYP/6-31+G//B3LYP/6-31G single point energies are essentially the same as those obtained from full geometry optimizations at the B3LYP/6-31+G level. The northern and southern minima of the B3LYP/6-31+G surface are, respectively, the E(2) and (2)E conformers. The B3LYP/6-31G geometries were used to study the relationship between ring conformation and various structural parameters including bond angles, dihedral angles, bond lengths, and interatomic distances.     

Mechanistic insights into triterpene synthesis from quantum mechanical calculations. Detection of systematic errors in B3LYP cyclization energies

Most quantum mechanical studies of triterpene synthesis have been done on small models. We calculated mPW1PW91/6-311+G(2d,p)//B3LYP/6-31G* energies for many C30H51O+ intermediates to establish the first comprehensive energy profiles for the cationic cyclization of oxidosqualene to lanosterol, lupeol, and hopen-3beta-ol. Differences among these 3 profiles were attributed to ring strain, steric effects, and proton affinity. Modest activation energy barriers and the ample exothermicity of most annulations indicated that the cationic intermediates rarely need enzymatic stabilization. The course of reaction is guided by hyperconjugation of the carbocationic 2p orbital with parallel C-C and C-H bonds. Hyperconjugation for cations with a horizontal 2p orbital (in the plane of the ABCD ring system) leads to annulation and ring expansion. If the 2p orbital becomes vertical, hyperconjugation fosters 1,2-methyl and hydride shifts. Transition states leading to rings D and E were bridged cyclopropane/carbonium ions, which allow ring expansion/annulation to bypass formation of undesirable anti-Markovnikov cations. Similar bridged species are also involved in many cation rearrangements. Our calculations revealed systematic errors in DFT cyclization energies. A spectacular example was the B3LYP/6-311+G(2d,p)//B3LYP/6-31G* prediction of endothermicity for the strongly exothermic cyclization of squalene to hopene. DFT cyclization energies for the 6-311+G(2d,p) basis set ranged from reasonable accuracy (mPW1PW91, TPSSh with 25% HF exchange) to underestimation (B3LYP, HCTH, TPSS, O3LYP) or overestimation (MP2, MPW1K, PBE1PBE). Despite minor inaccuracies, B3LYP/6-31G* geometries usually gave credible mPW1PW91 single-point energies. Nevertheless, DFT energies should be used cautiously until broadly reliable methods are established.