烃类化合物碳氢键解离焓的密度泛函理论研究

通过比较10种密度泛函方法对烃类化合物碳氢键解离焓的计算精度, 发现新型密度泛函BMK方法具有最高的计算精度. 利用该方法计算了包含饱和链烃,、不饱和链烃、脂环烃和芳香烃在内的172个烃类化合物的碳氢键解离焓,计算均方根误差仅为7.95 kJ•mol^－1, 线性拟合常数为0.985. 通过自然键轨道法分析发现, 烃类物质的碳氢键解离焓与母体的碳氢键杂化轨道成分p%, 自由基奇电子轨道杂化成分p%及自由基的自旋密度三个参数之间存在较好的定量关系. 此外, 饱和链烷烃及不饱和链烃的碳氢键解离焓与碳氢键键长之间也存在较好的线性关系.     

密度泛函法计算C1-C14正构烃生成焓及C-C键裂解能

采用密度泛函理论(DFT)方法研究了费托石脑油裂解反应中涉及到C₁-C₁₄正构烃和自由基中间体的生成焓及其C-C键解离能(BDE)。 结果表明,在所有评价的密度泛函理论方法（B97-1、BB1K、B1B95、MPWB1K和MPW1B95）中,MPW1B95/6-311G(d,p)方法计算最精确。 以此方法为基准,进一步对高碳烃及其裂解产物的标准生成焓和C-C键解离能进行了预测。 与可得到的实验数据相比,MPW1B95/6-311G(d,p)方法预测的烃和自由基的平均生成焓分别为0.8和2.7 kJ/mol,C-C键解离能的平均绝对误差只有3.1 kJ/mol,表明此方法不仅可准确计算正构烃标准生成焓和C-C键解离能,而且还能正确预测C-C键解离能变化趋势。     

取代基效应对褐煤模型化合物离解焓影响的理论研究

为了探究褐煤热解过程中氧桥键C-O均裂这一重要反应, 选取α-O-4和β-O-4类结构单元作为褐煤模型化合物, 运用不同密度泛函计算了部分模型化合物中C-O的离解焓, 并以CBS-QB3作为理论基准值进行比较, 最后选取M05-2X进行离解焓计算. 结果显示, 对于选定的α-O-4和β-O-4类模型化合物, 其平均离解焓分别为51.0 kcal/mol和66.1 kcal/mol. 周围取代环境能显著影响C-O离解焓, 芳环上存在给电子基团(OH, OCH₃和CH₃)能降低C-O离解焓, 而吸电子基团COOH则能增加其离解焓. 然后深层次分析了取代基效应对C-O离解焓的影响. 此外, 分子内氢键的形成对离解焓也有很大的影响. C-O的离解焓与其键长没有特定的相关性, 不能简单的通过C-O键长来预测其离解焓.     

N—O键解离焓的密度泛函理论研究

在高精度计算方法G3和G3B3的基础上,比较了密度泛函理论(DFT)十几种方法对N—O键解离焓(BDE)相对于实验值的计算精度,发现用B3P86方法计算15种化合物N—O键的BDE,均方根误差最小,仅为6.36kJ·mol-1,计算值与实验值的线性相关系数为0.991.在此基础上,用该方法分别计算了非芳香化合物及芳香化合物的N—O键BDE.通过自然键轨道分析,发现部分N—O键的BDE与N—O键的键长、原子电荷密度及键级之间存在定量关系.此外,在B3P86方法的基础上预测了几种典型的杂环芳香化合物N—O键BDE值.     

槲皮素抗氧化活性的密度泛函理论研究

采用杂化密度泛函理论(DFT)方法, 预测了黄酮类化合物槲皮素分子的几何结构、电子结构和脱氢解离焓, 分析了这些性质与分子活性位的关系, 探讨了槲皮素分子的抗氧化活性, 即与活性氧自由基•OH, •OOH和 的反应机理. 在B3LYP/6-31＋G(d)水平下, 计算得到的槲皮素分子脱氢自由基的相对稳定性、脱氢解离焓和氢提取过程的活化能都表明, 槲皮素中的4 -羟基活性最高, 最有可能参与自由基的清除. 4 -羟基位的这种反应活性主要来源于相邻羟基之间的弱氢键相互作用. 深入研究槲皮素分子的抗氧化机理, 有助于更合理地设计和合成新的抗氧化剂.     

应用理论化学方法预测有机分子的标准摩尔生成焓

介绍一个针对高年级本科生开设的计算化学创新性实验,采用高斯理论计算一系列有机小分子的构型和热化学性质,并通过原子化反应法计算其标准摩尔生成焓.该实验可让学生初步掌握相关计算软件的使用,了解不同计算方法的适用范围,认识理论与实验相结合的实践途径.     

金红石型TiO2表面嵌入簇模型的密度泛函研究

应用金属氧化物簇模型的选取原则和密度泛函方法,对一系列嵌入族模型金红石型(TiO₂)n(n=2～15)的表面体系进行了研究,计算结果表明,合理选取的嵌入簇模型可以得到与实验一致的固体表面电子结构,通过描述表面的性质能够给出表面金属吸附的电子行为的定性解释,进一步证实了该模型的密度泛函研究能够用于固体金属氧化物的簇-表面类比.揭示了选择合理的簇模型是保证族表面类比方法成功的主要因素.     

均四嗪衍生物结构、生成焓及爆轰性能的理论研究

运用DFT-w B97/6-31+G**方法对23种1,2,4,5-四嗪衍生物的几何结构、自然键轨道(NBO)和生成焓(EOF)进行研究,并在此基础上运用Kamlet-Jacobs方程估算衍生物的爆轰性能,得到其爆速在6.69~9.37 km/s之间;基于统计热力学,求得部分标题化合物在200~800 K温度范围内的热力学性质,随温度T升高,热容Cp、熵Sm及焓Hm逐渐增大。根据最小键级理论,C-R(取代基)键和N-R键可能是1,2,4,5-四嗪衍生物高温裂解的热引发键。综合分析,基团-NO2、-N3和-N=N-有助于提高四嗪衍生物的生成焓和爆轰性能,3,6-二硝基-1,2,4,5-四嗪和3,6-二偶氮基-二硝基-1,2,4,5-四嗪从能量、爆轰性能上可以作为高能量密度材料候选物。     

取代基对木质素三聚体模型化合物醚键均裂反应影响的密度泛函理论研究

利用密度泛函理论M062X/6-31++G(d,p)方法,对27种具有不同取代基(甲基、羟甲基和甲氧基)的木质素三聚体模型化合物的Cα-O和Cβ-O键均裂解离能进行了理论计算,探究了不同位置取代基对醚键解离能的影响规律。结果表明,当R2或R3位氢原子仅有一个被甲氧基取代时,Cβ-O键解离能变化很小;当R2、R3位氢原子均被甲氧基取代时,Cβ-O键解离能明显降低;且R4、R5位甲氧基能强化R2、R3位甲氧基对Cβ-O键解离能的降低程度,而不受R1位取代基的影响。当R4、R5位氢原子相继被甲氧基取代时,Cα-O键解离能逐渐降低,且R2、R3位甲氧基也能强化R4、R5位甲氧基对Cα-O键解离能的降低程度。当R1位氢原子相继被甲基、羟甲基取代时,Cα-O键解离能逐渐升高,然而R2、R3位甲氧基会弱化R1位甲基、羟甲基对Cα-O键解离能的升高程度;R1位甲基不会影响Cβ-O键解离能,羟甲基却能明显提高Cβ-O键解离能。     

N2O在焦炭表面异相生成和分解机理的密度泛函理论研究

选用合理简化的焦炭模型,对煤焦燃烧过程中N₂O的异相生成和分解机理进行了分子水平上的研究。采用UB3LYP/6-31G(d)密度泛函理论方法优化得到了反应路径上反应物、产物、中间体和过渡态的几何构型和各中间反应的活化能和反应焓变。NO与其预先吸附在焦炭表面解离生成的表面氮组分反应生成N₂O的路径有两个,需要克服的势垒分别为69.3kJ/mol和200.0kJ/mol;NO亦可直接与焦炭中的吡啶氮结合释放出N₂O,该反应路径所需克服的最大势垒为418.0kJ/mol。N₂O可在焦炭表面分解释放出N₂,异相分解反应为一步反应,计算所得活化能为100.8kJ/mol。N₂O的异相生成和异相分解反应均为放热反应。采用经典过渡态理论计算得到了各路径中速率控制步骤的反应速率常数。低温条件下,N₂O的异相分解反应速率略低于其异相生成速率,随着温度的升高,两者逐渐接近,说明高温条件有利于N₂O的异相分解。     

碳氟键均裂解离能的密度泛函理论研究

运用6种密度泛函方法(B3LYP, B3P86, B3PW91, PBE1PBE, MPW1B95, MPW1K)对15个含氟有机化合物的碳氟键均裂解离能进行理论计算, 得到的理论值与实验值比较, 发现B3P86方法用于碳氟键均裂解离能的计算相对可靠. 使用验证后的理论方法对含氟杂环有机化合物和卤氟烃中的碳氟键均裂解离能进行了预测和分析, 并进一步讨论了α-取代基效应以及Hammett型取代基效应对碳氟键均裂解离能的影响.     

Studying the Chemistry of Cationized Triacylglycerols Using Electrospray Ionization Mass Spectrometry and Density Functional Theory Computations

Analysis of triacylglycerols (TAGs), found as complex mixtures in living organisms, is typically accomplished using liquid chromatography, often coupled to mass spectrometry. TAGs, weak bases not protonated using electrospray ionization, are usually ionized by adduct formation with a cation, including those present in the solvent (e.g., Na⁺). There are relatively few reports on the binding of TAGs with cations or on the mechanisms by which cationized TAGs fragment. This work examines binding efficiencies, determined by mass spectrometry and computations, for the complexation of TAGs to a range of cations (Na⁺, Li⁺, K⁺, Ag⁺, NH₄ ⁺). While most cations bind to oxygen, Ag⁺ binding to unsaturation in the acid side chains is significant. The importance of dimer formation, [2TAG + M]⁺ was demonstrated using several different types of mass spectrometers. From breakdown curves, it became apparent that two or three acid side chains must be attached to glycerol for strong cationization. Possible mechanisms for fragmentation of lithiated TAGs were modeled by computations on tripropionylglycerol. Viable pathways were found for losses of neutral acids and lithium salts of acids from different positions on the glycerol moiety. Novel lactone structures were proposed for the loss of a neutral acid from one position of the glycerol moiety. These were studied further using triple-stage mass spectrometry (MS³). These lactones can account for all the major product ions in the MS³ spectra in both this work and the literature, which should allow for new insights into the challenging analytical methods needed for naturally occurring TAGs.

Accurate Prediction of IrH Bond Dissociation Enthalpies by Density Functional Theory Methods

The iridium hydride complexes have been extensively used in organic reactions, such as oxidation and hydrogenation reactions. In many of these reactions, the dissociation or formation of Ir? H bond plays an important role in determining the overall reaction rates and yields. In the present study, the accuracy of different theoretical methods for prediction of Ir? H bond strengths has been examined on the basis of the previously reported Ir? H BDEs of 17 different complexes. Comparing the performance of different DFT functionals (e.g. B3LYP, TPSS, M06), different basis sets (including the different effective core potentials (ECP) on Ir and I atoms, and the total electron basis sets on the other atoms), and different solvation models (SMD, CPCM, and IEFPCM) in solution phase single point calculations, we found that the gas‐phase calculation with TPSS/(LanL2DZ: 6‐31G(d)) method is relatively more accurate than the other gas‐phase calculation methods, and can well simulate the Ir? H BDEs in low‐polarity solvents (such as chlorobenzene and dichloroethane). Finally, efforts were put in analyzing the structure‐activity relationships between the ligand structure (around Ir center) and the Ir? H BDEs. We wish the present study could benefit future studies on the Ir‐H complexes involved organic reactions.     

Polynomial Dependence for the Calculation of Vaporization Enthalpies of Organophosphorus Compounds

Enthalpies of vaporization of organic phosphorus compounds are related to the normal-pressure boiling points through a polynomial dependence. Each class of organophosphorus compounds is characterized by its own specific parameters of the general dependence. The polynomial dependences may be used for the calculation of vaporization enthalpies of organophosphorus compounds with different modes of phosphorus coordination.     

硝基烃光异构化反应的密度泛函理论计算

采用DFT(B3LYP)计算方法, 在6-31G*水平上获得了反式-β-硝基苯乙烯、硝基乙烯和硝基甲烷基态异构化反应时的过渡态分子结构, 并计算了异构化能垒及激发态电子跃迁能. 结果显示, 反式-β-硝基苯乙烯和硝基苯与硝基甲烷相比具有较短的过渡态C—N键长, 较低的异构化能垒, 并且随着不饱和度的增加, 硝基苯和反式-β-硝基苯乙烯电子垂直跃迁能与基态异构化反应过渡态之间能量的差值ΔE迅速减小. 从能量的角度分析, 取代基的不饱和度越大, 越有利于激发态势能面与异构化反应势能面发生锥型或漏斗交叉, 因而越有利于光化学反应沿光异构化通道进行. 激发态分子的初始电子运动的定域或离域特征的差别可能是导致硝基苯等硝基芳烃与硝基甲烷等硝基烷烃光解通道不同的一个重要原因.     

Dynamic Combinatorial Chemistry: A New Methodology Comes of Age

Dynamic combinatorial chemistry (DCC) has repeatedly proven to be an effective approach to generate directed ligand libraries for macromolecular targets. In the absence of an external stimulus, a dynamic library forms from reversibly reacting building blocks and reaches a stable thermodynamic equilibrium. However, upon addition of a macromolecular host which can bind and stabilize certain components of the library, the equilibrium composition changes and induces an evolution-like selection and enrichment of high-affinity ligands. A valuable application of this so-called target-directed DCC (tdDCC) is the identification of potent ligands for pharmacologically relevant targets. Over time, the term tdDCC has been applied to describe a number of different experimental setups, leading to some ambiguity concerning its definition. This article systematically classifies known procedures for tdDCC and related approaches, with a special focus on the methods used for analysis and evaluation of experiments.     

Constructing Homodesmic Reactions for Calculating the Enthalpies of Formation of Organic Compounds

The mathematics and software of the homodesmic approach to calculating enthalpies of formation Δ_fH° of organic compounds are described. The program generates a set of independent homodesmic reactions for a studied compound, increasing the reliability of theoretical determinations of the standard enthalpy of formation. The basic steps of program construction, using the algorithm to determine the basis of homodesmic reactions for organic compounds, and calculations for the enthalpies of formation for a random test set comprising 53 CHNO-containing organic compounds of different classes are described. It is shown that the homodesmic approach provides highly reliable theoretical values of Δ_fH° virtually regardless of the quantumchemical approach that is used. The average absolute error of MAD calculations for Δ_fH° of a test set is 2.2 kJ/mol for the simple B3LYP/6-31G(d) approach and 2.1 kJ/mol for the moderately complicated M06-2X/cc-pVTZ approach, while using the complicated method of G3 calculations reduces the MAD value to 1.7 kJ/mol.     

有机二硫化物润滑添加剂抗磨极压性能的密度泛函理论研究

用量子化学的密度泛函理论计算了12种有机二硫化物和铁原子簇的分子轨道指数及其与铁原子簇的化学吸附作用能, 探讨了这种作用能与抗磨性能的关系; 运用轨道能量近似原则讨论了有机二硫化物与铁原子的作用方式; 以前线电子密度、超离域性指数和原子净电荷作为判据分析了12种有机二硫化物与铁原子间键合的强弱、反应性的大小等表征有机二硫化物与金属作用强弱的参数。结果表明: 有机二硫化物与铁接触时, 在较缓和条件下, SS键优先断裂与金属发生化学吸附形成配位键, 起到抗磨作用; 在高负荷下, 与金属发生常规条件下不能发生的化学反应, 即CS键断裂生成无机膜, 起到极压作用; 且随着碳链的增长, 有机二硫化物的抗磨性能愈来愈好, 但极压性能愈来愈差; 运用量子化学计算得到的预测结果与摩擦学试验结果具有良好的一致性, 可为同类极压添加剂化合物的分子设计提供较为可靠的参考依据和理论方法。     

过渡金属类立方烷原子簇化合物电子结构和二阶非线性光学性质的密度泛函理论研究

采用基于第一原理的含时密度泛函理论(TDDFT)对一系列具有类立方烷簇芯结构的过渡金属簇合物二阶非线性光学性质进行了研究。结果证明, 由于簇芯结构的对称性的影响, 这一类簇分子的二阶非线性光学系数的数值要小于三核欠完整类立方烷体系。通过对电子结构的分析, 发现二阶非线性光学性质主要是由簇芯内电荷的迁移造成的。轨道分析显示, -S原子对于电荷的迁移起主要的传递作用。定域化轨道分析证明簇分子中存在的多中心键有利于簇芯内电荷的迁移。分子模拟的研究表明:虽然类立方烷结构簇分子的值比较小, 但是通过合理的配体设计, 获得具有较大值非线性光学晶体是可能的。