7-Ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A): a DFT study of the antioxidant mechanism 2.* The structure of monosodium salts of echinochrome A and their reactions with the hydroperoxyl radical

The relative gas-phase acidities of all OH groups of 7-ethyl-2,3,5,6,8-pentahydroxy-1,4-naphthoquinone (echinochrome A, 1) were evaluated by the B3LYP/6-311G(d) and B3LYP/6-311G(d,p) methods. Calculations predict that gB-OH groups at the C(2) and C(6) atoms are the most acidic in molecule 1and their acidity is higher than that of o-nitrophenol. Conformational analysis of undissociated monosodium salts of 1and their radicals was performed. It was shown that gas-phase quenching reactions of the hydroperoxyl radical by mono-sodium salts of 1are exothermic.     

2-溴丙酸气相热消除反应的机理

用密度泛函（DFT）方法,在B3LYP/6-31G**水平上对2-溴丙酸气相消除反应机理进行了研究.计算表明,反应主要是通过半极化五元环结构过渡态进行的,羧基上的氢原子协助溴原子离去,羧基氧原子帮助稳定过渡态.在B3LYP/6-311＋＋G(3df,3pd)水平上对B3LYP/6-31G**优化的几何构型进行了单点能计算,计算所得反应的速度控制步骤的活化能为189.461 kJ•mol－1,偏离实验值((180.3±3.4) kJ•mol－1)5.08％.     

椅式(8,8)单壁碳纳米管内偶氮苯的顺反异构化

采用ONIOM(B3LYP/6-31+G*:UFF)方法对偶氮苯受限于椅式单壁碳纳米管CNT(8,8)内的结构、电子光谱和热致顺反异构化势能面进行了计算.结果表明,反式偶氮苯进入CNT(8,8)碳纳米管内是一个放热过程.由于碳纳米管的限制作用,偶氮苯分子的苯环绕CN键有一定的旋转,反式偶氮苯的平面结构发生扭曲,但其它结构参数变化不明显.受限于碳纳米管内的偶氮苯的顺反异构体能量差比非受限状态下增加了8.1 kJ.mol-1,表明碳纳米管的限制作用对偶氮苯顺反异构体的相对热稳定性有一定影响.光谱计算表明,受限于碳纳米管内的偶氮苯最低的三个单线态吸收波长仅蓝移1-5 nm.异构化势能面计算发现,偶氮苯在CNT(8,8)碳纳米管内发生反式到顺式异构化的能垒增加,由顺式回复到反式异构体的能垒无明显变化,表明CNT(8,8)碳纳米管的限制作用可以抑制偶氮苯从反式到顺式异构体的热致异构化过程.此外,受限状态下的偶氮苯主要通过CNN键角反转发生热致异构化.     

Cyclic vinyl p-tolyl sulfilimines as chiral dienophiles: theoretical study on the stereoselectivity, Lewis acid catalysis, and solvent effects in their Diels-Alder reactions

The theoretical study reported in the present work deals with chiral cyclic vinyl sulfilimines and their reactivity as dienophiles in [4 + 2] cycloaddition reactions, using B3LYP/6-31G(d)//AM1 and B3LYP/6-31G(d)//B3LYP/6-31G(d) model chemistries. Consideration of Lewis acid catalysis, illustrated by BF(3), decreases the activation energies of the cycloaddition process while the charge transfer from the diene to the sulfilimine is augmented. The [4 + 2] cycloaddition reactions of sulfilimines with both furan and cyclopentadiene occur in the gas phase with endo stereoselectivity, which is more pronounced with the latter diene. Endo-exo energy differences in the gas phase with the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31G(d)//AM1 model chemistries are almost the same. Solvent effects are responsible for the inversion of the stereoselectivity in the reactions of sulfilimines with furan because of the great difference in the dipole moments in endo and exo approaches.     

Predicting Rate Constants for Nucleophilic Reactions of Amines with Diarylcarbenium Ions Using an ONIOM Method

The rate constants of the nucleophilic reactions between amines and benzhydrylium ions were calculated using first-principles theoretical methods. Solvation models including PCM, CPCM, and COSMORS, as well as different types of atomic radii including UA0, UAKS, UAHF, Bondi, and UFF, and several single-point energy calculation methods (B3LYP, B3P86, B3PW91, BHANDH, PBEPBE, BMK, M06, MP2, and ONIOM method) were exam-ined. By comparing the correlation between experimental rate constants and the calculated values, the ONIOM(CCSD(T)/6-311++G(2df,2p):B3LYP/6-311++G(2df,2p))// B3LYP/6-31G(d)/PCM/UFF) method was found to perform the best. This method was then employed to calculate the rate constants of the reactions between diverse amines and diarylcarbenium ions. The calculated rate constants for 65 reactions of amines with diarylcarbenium ions are in agreement with the experimental values, indicating that it is feasible to predict the rate constant of a reaction between an amine and a diarylcarbenium ion through ab initio calculation.     

适用于液晶大分子的量子化学计算方法以及六烷氧基苯并菲取代的1,3,5-三酰胺苯液晶的研究

盘状液晶分子之间的相互作用决定液晶的性质. 为了选择合适的计算方法, 以便用量子化学研究液晶大分子, 设计了对位取代苯和三酰胺苯作为模型分子, 用高水平的ONIOM [MP2/6-31G*(0.25):HF/6-31G(d,p)]计算提供了与晶体结构资料相符合的较准确的二体相互作用结果. 然后与各种较低级别的计算进行比较, 说明ONIOM (HF/STO-3G: AM1:UFF)水平比较合适. 盘重叠部分的相互作用主要是色散作用, 用UFF力场处理, 氢键主要是静电作用, 用HF/STO-3G处理, 其余部分用AM1处理. 通过对苯取代的1,3,5-三酰胺苯的双分子构型优化, 并与晶体结构数据进行比较, 进一步说明ONIOM (HF/STO-3G:AM1:UFF)水平计算的适用性. 在此基础上, 对六烷氧基苯并菲取代的1,3,5-三酰胺苯的双分子构型进行了优化, 为解释它所形成的液晶具有较高的电荷传输能力提供了有用的信息.     

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.     

Thermochemistry of 2-amino-3-quinoxalinecarbonitrile-1,4-dioxide. Evaluation of the mean dissociation enthalpy of the (N-O) bond

The standard enthalpy of formation of the 2-amino-3-quinoxalinecarbonitrile-1,4-dioxide compound in the gas-phase was derived from the enthalpies of combustion of the crystalline solid measured by static bomb combustion calorimetry and its enthalpy of sublimation determined by Knudsen mass-loss effusion at T= 298.15 K. This value is (383.8 +/- 5.4) kJ mol(-1) and was subsequently combined with the experimental gas-phase enthalpy of formation of atomic oxygen and with the computed gas-phase enthalpy of formation of 2-amino-3-quinoxalinecarbonitrile, (382.0 +/- 6.3) kJ mol(-1), in order to estimate the mean (N-O) bond dissociation enthalpy in the gas-phase of 2-amino-3-quinoxalinecarbonitrile-1,4-dioxide. The result obtained is (248.3 +/- 8.3) kJ mol(-1), which is in excellent agreement with the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d) computed value.     

Ab initio and density functional theory studies of the structure,gas-phase acidity and aromaticity of tetraselenosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure,acidity (free energy of deprotonation,G0) and aro-maticity of tetraselenosquaric acid (3,4-diselenyl-3-cy-dobutene-1,2-diselenone,H2C4Se4) are reported.The global minimum found on the potential energy surface of tetraselenosquaric acid presents a planar conformation.The ZZ iso-mer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very dose in energy.The optimized geometric parameters exhibit a bond length equalization relative to reference compounds,cyclobu-tanediselenone,and cydobutenediselenol.The computed aromatic stabilization energy (ASE) by homodesmotic reaction is -77.4 (MP2(fu)/6 - 311 G //RHF/6 - 311 G) and - 54.8 kJ/mol (B3LYP/6 - 311 G //B3LYP/6 -311 G).The aromaticity of tetraselenosquaric add is indicated by the calculated diamagnetic susceptibility exaltation (A) - 19.13 (CSGT(IGAEM) - RHF/6 - 311 G// RHF/6-     

Intermolecular Interaction of HMX： an Application of ONIOM Methodology

IntroductionExtensive studies on intermolecular interac-tions have been made in the past decades due totheir importance in a wide range of physical,chem-ical and biological fields.Researches on the weakintermolecular interactions began with hydrogenbonds.Scheiner summarized the ab initio investiga-tions on hydrogen bonding in detail[1] .With the abinitio method supermolecular structures and bind-ing energies can be predicted notonly for H- bondedsystems but also for other systems[2 _ 5] when …     

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.     

Theoretical predictions of the structure, gas-phase acidity and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid

Results of ab initio self-consistent-field (SCF) and density functional theory (DFT) calculations of the gas-phase structure, acidity (free energy of deprotonation, ΔG^o), and aromaticity of 1,2-diseleno-3,4-dithiosquaric acid (3,4-dithiohydroxy-3-cyclobutene-1,2-diselenone, H₂C₄Se₂S₂) are reported. The global minimum found on the potential energy surface of 1,2-diseleno-3,4-dithiosquaric acid presents a planar conformation. The ZZ isomer was found to have the lowest energy among the three planar conformers and the ZZ and ZE isomers are very close in energy. The optimized geometric parameters exhibit a bond length equalization relative to reference compounds, cyclobutanediselenone, and cyclobutenedithiol. The computed aromatic stabilization energy (ASE) by homodesmotic reaction (Eq 1) is −20.1 kcal/mol (MP2(fu)/6-311+G** //RHF/6-311+G**) and −14.9 kcal/mol (B3LYP//6-311+G**//B3LYP/6-311+G**). The aromaticity of 1,2-diseleno-3,4-dithiosquaric acid is indicated by the calculated diamagnetic susceptibility exaltation (Λ) −17.91 (CSGT(IGAIM)-RHF/6-311+G**//RHF/6-311+G**) and −31.01 (CSGT(IGAIM)-B3LYP/6-311+G**//B3LYP/6-311+G**). Thus, 1,2-diseleno-3,4-dithiosquaric acid fulfils the geometric, energetic and magnetic criteria of aromaticity. The calculated theoretical gas-phase acidity is ΔG^o _1(298K)=302.7 kcal/mol and ΔG^o _2(298K)=388.4 kcal/mol. Hence, 1,2-diseleno-3,4-dithiosquaric acid is a stronger acid than squaric acid(3,4-dihydroxy-3-cyclobutene-1,2-dione, H₂C₄O₄). Received: 11 April 2000 / Accepted: 7 July 2000 / Published online: 27 September 2000     

Phosphinofenchol or metastable phosphorane? Phosphorus derivatives of fenchol

Not the expected phosphinofenchol 1 but phosphorane 2 is obtained after reaction of 2-lithio(diphenylphosphino)benzene with (-)-fenchone. Surprisingly, ONIOM(B3LYP/6-31G*:UFF) computations of 1 and 2 as well as B3LYP analyses of smaller model systems point to a lower thermodynamic stability of phosphoranes relative to their isomeric alkoxyphosphines. An analogue inherent instability is computed for the methylphosphorane 10, which is also synthesized and characterized by X-ray analysis. Decreasing ring size in cyclic phosphoranes, that is, from five- to four-membered ring systems, destabilizes cyclic phosphoranes even more. This computational prediction is verified experimentally by reaction of lithiomethyl(diphenylphosphine) with (-)-fenchone and subsequent isolation of the corresponding phosphinofenchol. Protonation or alkylation of phosphoranide intermediates can account for the formation of metastable phosphoranes.     

Intra- and intermolecular 1,3-dipolar cycloaddition of sugar ketonitrones with mono-, di-, and trisubstituted dipolarophiles

The 1,3-dipolar cycloaddition of sugar ketonitrones is a useful synthetic procedure to build up nitrogenated quaternary centers in terms of scope (substrate, dipolarophile, inter- and intramolecular versions), yield, and regio- and stereoselectivity. The hybrid ONIOM (B3LYP/6-31G(d):AM1) theoretical method followed by single-point energy calculations at the B3LYP/6-31G(d) level adequately perform to model this cycloaddition for the relatively large ketosugar precursors commonly used.     

Proton and sodium cation affinities of harpagide: a computational study

The aim of this work was to estimate the proton and sodium cation affinities of harpagide (Har), an iridoid glycoside responsible for the antiinflammatory properties of the medicinal plant Harpagophytum. Monte Carlo conformational searches were performed at the semiempirical AM1 level to determine the most stable conformers for harpagide and its protonated and Na+-cationized forms. The 10 oxygen atoms of the molecule were considered as possible protonation and cationization sites. Geometry optimizations were then refined at the DFT B3LYP/6-31G level from the geometries of the most stable conformers found. Final energetics were obtained at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G level. The proton and sodium ion affinities of harpagide have been estimated at 223.5 and 66.0 kcal/mol, respectively. Since harpagide mainly provides HarNa+ ions in electrospray experiments, the DeltarG298 associated with the reaction of proton/sodium exchange between Har and methanol, MeOHNa+ + HarH+ --> MeOH2+ + HarNa+ (1), has been calculated; it has been estimated to be 1.9 kcal/mol. Complexing a methanol molecule to each reagent and product of reaction 1 makes the reaction become exothermic by 1.7 kcal/mol. These values are in the limit of the accuracy of the method and do not allow us to conclude definitely whether the reaction is endo- or exothermic, but, according to these very small values, the cation exchange reaction is expected to proceed easily in the final stages of the ion desolvation process.     

Oxygen atom transfer energetics: assessment of the effect of method and solvent

Several density functional methods, the semiempirical methods AM1 and PM3, Hartree-Fock, and Gaussian3 theories were applied to compute the oxygen atom transfer enthalpies for 14 X/XO couples (inorganic and organic systems, charged and neutral species, light and heavy main group element containing molecules). The calculated reaction enthalpies were compared to available experimental data. The G3 method alone was found to perform within the experimental error, while the popular B3LYP and BLYP functionals provided inadequate results. Solvent effects were estimated for 19 neutral and anionic X/XO couples by using the conductor-like polarizable continuum model and several cavity models coupled with the B3LYP/6-31++G(2d,2p) level of theory. Surprisingly, the magnitude of the aqueous solvent correction was found to vary significantly for different solute cavity models, occasionally giving larger errors than the gas-phase calculation.     

Ga/HZSM-5分子筛上乙烯二聚反应的理论研究

应用理论计算方法研究了Ga/HZSM-5及Al/HZSM-5 分子筛上乙烯二聚生成1-丁烯的反应历程, 比较了分子筛酸性对反应能量的影响. 计算采用分为两层的76T簇模型, 应用量子力学和分子力学联合的 ONIOM2 (B3LYP/6-31G(d, p):UFF) 方法. 乙烯二聚过程可按照分步机理和协同机理进行, 均得到表面丁基烷氧中间产物. 计算结果表明, 与在Al/HZSM-5分子筛上的反应过程相比, 乙烯在 Ga/HZSM-5分子筛上的吸附能低 20.62 kJ·mol^-1, 但质子化反应的活化能只高出1.26 kJ·mol^-1; 而乙基烷氧中间体与乙烯分子结合过程的活化能高出 62.55 kJ·mol^-1, 原因是Ga 原子半径大, 降低了六元环过渡态的稳定性. 若按协同机理, 质子转移和C―C键聚合同时进行, 在 Ga/HZSM-5分子筛上的活化能较Al/HZSM-5的高16.44 kJ·mol^-1. 因此乙烯二聚按照协同机理有利. 研究还表明, 表面丁基烷氧中间体脱质子, 生成1-丁烯并吸附在复原的分子筛酸性位上. 该反应在两种酸中心上的活化能几乎相同, 但明显高于其他各步的活化能, 因此成为整个反应的速度控制步骤.     

Polar substituent effects in the bicyclo[1.1.1]pentane ring system: acidities of 3-substituted bicyclo[1.1.1]pentane-1-carboxylic acids

Experimental gas-phase acidities are reported for a series of 3-substituted (X) bicyclo [1.1.1]pent-1-yl carboxylic acids (1, Y = COOH). A comparison with available calculated data (MP2/6-311++G**// B3LYP/6-311+G**) reveals good agreement. The relative substituent effects are shown to be adequately described by a much lower level of theory (B3LYP/6-31+G*). Various correlations are presented which clearly point to polar field effects as being the origin of the relative acidities.     

An experimental and theoretical study of the basicity of tetra-tert-butyltetrahedrane

The gas-phase basicity (GB) of tetra-tert-butyltetrahedrane (tBu4THD) was determined by FT-ICR mass spectrometry and comparison with reference compounds of known basicity. Its GB, 1035+/-10 kJ x mol(-1), makes tetra-tert-butyltetrahedrane one of the strongest bases reported so far. Ab initio calculations [B3LYP/6-31G(d) and B3LYP/6-311 + G(d,p)//6-31G(d)] have been carried out in order to compare the high experimental basicity of tBu4THD with that estimated theoretically. Both B3LYP/6-31G(d) and QCISD(T) calculations were used to determine the reaction path which connects the initial tetrahedrane-ammonium complex with the final products, protonated cyclobutadiene (CBDH+) and ammonia.     

Structures and reaction mechanisms of propene oxide isomerization on H-ZSM-5: an ONIOM study

The isomerization mechanisms of propene oxide over H-ZSM-5 zeolite have been investigated via the utilization of 5T and 46T cluster models calculated by the B3LYP/6-31G(d,p) and the ONIOM(B3LYP/6-31G(d,p):UFF) methods, respectively. The reactions are considered to proceed through a stepwise mechanism: (1) the epoxide ring protonation, and concurrently the ring-opening, and (2) the 1,2-hydride shift forming the adsorbed carbonyl compound. Because of the asymmetric structure of propene oxide, two different C-O bonds (more or less substituted carbon atom sides) can be broken leading to two different types of products, propanal and propanone. The ring-opening step of these mechanisms is found to be the rate-determining step with an activation barrier of 38.5 kcal/mol for the propanal and of 42.4 kcal/mol for the propanone. Therefore, the propanal is predicted to be the main product for this reaction.