含能材料中键离解能的密度泛函理论计算

通过计算一个包括臭氧,硝基甲烷以及1,3,5 -三硝基-1,3,5-三氮杂环已烷(RDX)在内的典型系统的键离解能,对由四种交换/相关函数(BLYP,B3LYP,B3PW91和B3P86)加上不同的基函数组合而成的多种密度泛函方法的准确性进行了比较研究.结果表明:B3P86/6-31G**是计算该系统C-NO2,O-O和N-NO2键离解能的最可靠的方法.     

凝聚态硝基甲烷分解机理的密度泛函研究

用密度泛函理论在B3LYP/6-311++G(2d,2P)计算水平上对最低单态和最低三态的硝基甲烷分子进行了分子动力学计算分析,发现:基态硝基甲烷分子沿C-N键分解生成硝基和甲基反应通道上不存在过渡态,只能是在能量足够高的时候造成C-N键的断裂,键离解能为53.4kcal/mol;硝基甲烷分子在最低三态沿C-N键分解生成硝基和甲基的反应通道上,有一个活化能为87.8kcal/mol的能垒.计算得到硝基甲烷分子从基态到最低三态分解反应发生所需要的总能量为144.58kcal/mol.这个数值与硝基甲烷材料的电子碰撞实验在193nm处有强吸收峰的结论相符合.依据多声子迁移理论,结合硝基甲烷分子在最低三态动力学分解的可能性,可以认为在相同条件下,硝基甲烷材料在撞击条件下,分子沿CN键分裂生成硝基和甲基的反应在最低三态分子分解的可能性较大.文章用量化计算从分子构型、频率分析和势能面扫描方面对分析结论进行了加强和确定,并且,依据多声子迁移理论对硝基甲烷分子基态键离解过程、基态到三态激发过程和最低三态活化过程中的声子迁移进行了初步分析.     

凝聚态硝基甲烷分解机理的密度泛函研究

用密度泛函理论在B3LYP/6-311++G（2d,2P）计算水平上对最低单态和最低三态的硝基甲烷分子进行了分子动力学计算分析,发现,硝基甲烷分子基态沿C-N键分解生成硝基和甲基反应通道上不存在过渡态,只能是在能量足够高的时候造成C-N键的断裂,键离解能为53.4kcal/mol;硝基甲烷分子在最低三态沿C-N键分解生成硝基和甲基的反应通道上,有一个活化能为87.8kcal/mol的能垒。计算得到硝基甲烷分子从基态到最低三态份分解反应发生所需要的总能量为144.58kcal/mol.这个数值与硝基甲烷材料的电子碰撞实验在193nm处的有强吸收峰的结论相符合。依据多声子迁移理论,结合硝基甲烷分子在最低三态动力学分解的可能性,可以认为在相同条件下,硝基甲烷材料在撞击条件下,分子沿C-N键分裂生成硝基和甲基的反应在最低三态分子分解的可能性较大。文章用量化计算从分子构型、频率分析和势能面扫描方面对分析结论进行了加强和确定。并且,依据多声子迁移理论对硝基甲烷分子基态键离解过程、基态到三态激发过程和最低三态活化过程中的声子迁移进行了初步分析。     

氘代甲烷几何构型及物性的量子化学研究

用HF/6-31G^**、密度泛函方法B3LYP/31G^**、二级微扰MP2/6-31G^**、四级微扰MP4/6-31G^**方法对甲烷和氘代甲烷进行几何构型全优化,并将优化的结果与实验值进行比较.用上述4种方法对甲烷和氘代甲烷分子进行分子的振动基频计算.密度泛函、二级微扰、四级微扰优于HF/6-31G^**,尤其是密度泛函、四级微扰方法.密度泛函方法所用的机时远小于微扰方法.不同方法计算所得的氘代甲烷振动频率值与实验值的最大误差为10.4%,最小误差为2.0%.     

多硝基苯酸酯炸药分子键离解能与撞击感度关系研究

用密度泛函方法,在B3LYP-6-31G*、B3P86/6-31G*和B3LYP-6-311G*三种理论水平对四个含硝基烷基的苯酸酯炸药分子进行了几何结构全优化、能量和频率计算。并对这些炸药分子苯环上的C-NO2和烷基上的C-NO2键离解能分别进行了三种理论水平的计算。结果表明,这类分子中的最弱键是烷基上的C-NO2键。进一步分析实验撞击感度与分子中最弱键离解能量的关系,结果表明,分子最弱键离解能与分子总能量的比值BDE/E和实验撞击感度h50%之间存在一个几乎线性的关联关系。BDE/E是表征炸药撞击感度的一个实用的合理指针。     

含硝基烷基的一些CHNO多硝基化合物的C-NO2键离解能计算——应用于撞击感度

用密度泛函方法,在B3LYP/6-31G*水平对五个含硝基烷基的硝基芳香族炸药分子和四个含硝基烷基的苯酸酯炸药分子进行了几何结构全优化、能量和频率计算.并分别对这两类炸药分子苯环上的C-NO2和烷基上的C-NO2键离解能进行了同等水平的计算.结果表明,这两类分子中的最弱键均是烷基上的C-NO2键.进一步分析实验撞击感度与分子中最弱键离解能量的关系.结果表明,最弱键离解能BDE是表征炸药撞击感度的重要指针,但不是唯一指针.除了BDE,炸药的撞击感度可能还受其它因素或反应途径的影响.     

普适的基于能量的分块方法计算甲烷水合物簇的精确结合能和拉曼光谱

甲烷水合物在化学、能源和环境科学等领域中都具有重要作用. 本文采用普适的基于能量的分块(GEBF)方法计算了多种甲烷水合簇的结合能和拉曼光谱. 首先使用这些甲烷水合簇的在显相关耦合簇[CCSD(T)(F12^*)]水平下得到的GEBF结合能,评估了一系列密度泛函计算的结果. 计算结果表明B3PW91-D3和B97D泛函表现最佳,与GEBF-CCSD(T)(F12^*)基准相比的平均绝对误差分别仅为0.27和0.47 kcal/mol. 然后用GEBF-B3PW91-D3方法计算得到了单、双笼甲烷水合簇的结构和拉曼光谱,得到的甲烷C-H键伸缩拉曼振动峰与实验值的偏差小于3 cm^-1,说明B3PW91-D3泛函可以很好地重现实验结果. 随着水笼尺寸的增加. 甲烷C-H键伸缩拉曼峰发生红移,该现象与实验中提出的“松笼-紧笼”模型吻合. 此外,甲烷分子邻近环境(水笼)的改变对拉曼光谱的影响很小,环境从单笼变为双笼导致C-H键伸缩拉曼振动峰的蓝移不超过3 cm^-1. 甲烷水合簇的理论拉曼光谱与实验拉曼光谱结合可以用来研究海底或星际冰体内部的甲烷水合物的结构. 结合B3PW91-D3或B97D泛函和机器学习模型,可以进一步应用分子动力学模拟研究甲烷水合物的成核/生长机制和相变过程.     

二聚物X2 和 XK (X= Mn, Fe, Co, Ni, Cu)吸附NO的密度泛函理论研究

为有效降低NO对环境以及人体造成的危害,着重从催化剂入手,研究过渡金属Mn₂, Fe₂, Co₂, Ni₂, Cu₂及K掺杂的MnK, FeK, CoK, NiK, CuK的二聚物对NO的吸附性能.采用密度泛函理论的B3LYP, B3PW91, B3P86, B1B95以及PBE1PBE方法结合LANL2DZ, SDD, CEP-121G基组,采用B3LYP, B3PW91, B3P86以及PBE1PBE方法结合6-31++G(3df, 3pd), 6-311++G(3df, 3pd), 6-31G(d, p), LANL2DZ, SDD基组分别系统研究了X₂和XK (X=Mn, Fe, Co, Ni, Cu)团簇和NO的几何结构,计算出其键长,频率,离解能,再与对应的实验值进行对比.进一步采用PBE1PBE方法结合LANL2DZ和6-31G(d, p)基组研究了X₂NO和KXNO团簇的几何结构,计算出各原子间键长和吸附能.结果表明,...     

甘氨酸和丙氨酸二肽、三肽中N-H…O=C分子内氢键键能的理论研究

提出一种计算多肽中N-H…O=C分子内氢键键能的新方法.并将新方法应用于计算甘氨酸和丙氨酸二肽、三肽中N-H…O=C分子内氢键键能.利用密度泛函理论B3LYP/6-31G(d)方法优化几何构型和计算频率.对全部结构计算MP2/6-311 G(3df,2p)水平上的单点能量.结果表明:在甘氨酸二肽中氢键键能为-6.38 kcal/mol,在丙氨酸二肽中氢键键能为-7.09-、6.25 kcal/mol;在甘氨酸三肽中氢键键能为-5.62 kcal/mol,在丙氨酸三肽中氢键键能为-5.37,-5.74 kcal/mol.     

乙烯基自由基与IB金属小团簇的相互作用研究

采用密度泛函理论的B3LYP和PBE1PBE方法,对IB金属小团簇与C2H3自由基间的相互作用进行了研究.结果发现:C2H3M0/-n(M=Cu、Ag、Au;n=1-3)最稳定及其异构体的结构都相同,只是键长、键角有所不同.在C2H3M0/-n(M=Cu、Ag、Au;n=1-3)的最稳定结构中,C2H3自由基和金属团簇都作为整体单元存在.由于金属团簇与C2H3自由基的相互作用,在红外吸收光谱中,C2H3自由基中C=C、C-H键的伸缩振动发生了红移,而随着团簇中金属原子数目的增多,C=C键及C-H键的伸缩振动,红移程度减小.基于含时密度泛函(TDDFT)计算,模拟了C2H3M-n(M=Cu、Ag、Au;n=1-3)最稳定结构的光电子能谱(PES).     

Electrocyclic [1,5] hydrogen shift in the thermal elimination kinetics of phenyl acetate and p-tolyl acetate in the gas phase: a density functional theory study

The kinetics and mechanisms of thermal decomposition of phenyl acetate and p-tolyl acetate in the gas phase were studied by means of electronic structure calculations using density functional theory methods: B3LYP/6-31G(d,p), B3LYP/6-31++G(d,p), B3PW91/6-31G(d,p), B3PW91/6-31++G(d,p), MPW1PW91/6-31G(d,p), MPW1PW91/6-31++G(d,p), PBE/6-31G(d,p) and PBE/6-31++G(d,p). Two possible mechanisms have been considered: mechanism A is a stepwise process involving electrocyclic [1,5] hydrogen shift to eliminate ketene through concerted six-membered cyclic transition-state structure, followed by tautomerisation of cyclohexadienone or by 4-methyl cyclohexadienone intermediate to give the corresponding phenol. Mechanism B is a one-step concerted [1,3] hydrogen shift through a four-membered cyclic transition-state geometry, to produce ketene and phenol or p-cresol. Theoretical calculations showed reasonable agreement with experimental activation parameters when using the Perdew, Burke and Ernserhof (PBE)functional, through the stepwise [1,5] hydrogen-shift mechanism. For mechanism B, large deviation for the entropy of activation was observed. No experimental data were available for p-tolyl acetate; however, theoretical calculations showed similar results to phenyl acetate, thus supporting the stepwise mechanism for both phenyl acetate and p-tolyl acetate.     

硝酸丙酯结构、振动频率和热力学性质的密度泛函理论研究

用密度泛函理论(DFT)研究硝酸丙酯化合物的分子结构、振动光谱和热力学等基本性质.取BLYP、B3LYP方法和6-31G*、6-31G**、6-311G*、6-311G**基组,对硝酸丙酯分子的几何构型进行全优化计算并分析其电子结构性质.和考虑了二级相关能校正的MP2/6-311G*计算结果比较表明,B3LYP/6-31G*是研究许多较大体系化合物卓有成效和颇有前途的方法.在B3LYP/6-31G*的水平上对优化后的结构进行了正则振动频率分析,用因子0.95校正后的振动光谱和实验结果比较,符合较好.进一步     

Quantum chemical calculations of bond dissociation energies for COOH scission and electronic structure in some acids

Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies(BDEs) for 15 acids.These compounds are studied by utilizing the hybrid density functional theory(DFT)(B3LYP,B3PW91,B3P86,PBE1PBE) and the complete basis set(CBS-Q) method in conjunction with the 6311G** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital(HOMO) and the lowest unoccupied molecular orbital(LUMO) of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.     

聚对苯二甲酸丁二醇酯二聚体键离能的理论研究

采用密度泛函理论方法对聚对苯二甲酸丁二醇酯(PBT)二聚体的键离能进行了计算.为了选取较为精确的方法来计算PBT各个键的键离能,以与PBT具有相同的酯基官能团的乙酸乙酯为模型参照物.采用M062X, B3P86, M06, PBE0, wB97xD方法分别在基组6-31G(d), 6-311G(d), 6-311+G(d, p), 6-311++G(d, p), cc-pVDZ, cc-pVTZ水平下对乙酸乙酯的键离能进行计算.通过对比计算结果与iBonD数据库的乙酸乙酯实验测定值可知,M062X在基组6-311G(d)水平下计算结果与实验值最为接近.因此,本研究采用M062X方法在基组6-311G(d)水平下对聚对苯二甲酸丁二醇酯(PBT)二聚体的键离能进行计算.计算结果表明：在PBT的各键中C-C_arcmatic键的键离能最大,主链上的C-C键离能最小,为370.9 kJ/mol.其次就是C-O键,为404.6 kJ/mol.基于PBT键离能的计算结果,设计了3条PBT二聚体热降解过程可能形成的反应路径,分析了热解产物的形成机理.结果表明PBT二聚体热解过程可...     

A study of the role played by the Hartree-Fock orbital exchange in the formation of the energy of the first singlet charge-transfer excited state by the example of JK-62 and JK-201 sensitizing dye molecules

Based on the time-dependent density functional theory with the use of the functionals B3LYP, B97-2, BHandHLYP, BMK, MPWB1K, PBE1W, PBE1PBE, τ-HCTHh, and TPSS in the 6–31G(d) basis set of atomic orbitals, we have performed a quantum-chemical investigation of electronic and spectral properties of JK-62 and JK-201 bis-dimethylfluorenyl dye sensitizers for Grätzel photoelectric converters. In terms of the Bader theory, we have done a complete analysis of the electron-density distribution function in the dye molecules under study, which proves the occurrence of intramolecular nonvalent interactions, which, in turn, stabilize a planar mutual arrangement of structural fragments of dye molecules. The role that the Hartree-Fock orbital exchange plays in the energy formation of the first and most intense electronic transition, which is responsible for the primary current generation in a solar cell, has been elucidated.     

Quantum chemical study of a new class of sensitisers: influence of the substitution of aromatic rings on the properties of copper complexes

We present a computational study of new copper complexes with potential applications as sensitisers for solar cells. The applied methodology for this study is based on the density functional theory (DFT) and time-dependent DFT, using the B3LYP, PBE0 and M06 functionals with the LANL2DZ (D95V on first row), 6-31G(d,p), 6-311G(d,p) and DZVP basis sets. Optimised molecular structure, the absorption spectra, the molecular orbitals energies and the chemical reactivity parameters that arise from conceptual DFT were calculated. Solvent effects have been taken into account by an implicit approach, namely, the polarisable continuum model (PCM), using the non-equilibrium version of the integral equation formalism of the PCM model. Interesting work for experimentalists in the dye sensitised solar cells’ field.     

Theoretical Study of the C-Cl Bond Dissociation Enthalpy and Electronic Structure of Substituted Chlorobenzene Compounds

利用密度泛函(DFT)三种交换/相关函数(B3LYP, B3PW91,B3P86)结合6-31G**和6-311G**基组,计算了13个取代氯苯化合物的键离解能. 结果表明B3P86/6-311G**方法是计算取代氯苯化合物键离解能的可信方法,研究发现C-Cl键的键离解能与所使用的基组和计算方法密切相关,取代基对C-Cl键的键离解能的影响不明显. 研究了目标化合物的前线轨道能级差,并对取代氯苯化合物的热稳定性做了评估.     

The mechanism of the homogeneous,unimolecular gas‐phase elimination kinetic of 1,1‐dimethoxycyclohexane: experimental and theoretical studies

The gas‐phase elimination of 1,1‐dimethoxycyclohexane yielded 1‐methoxy‐1‐cyclohexene and methanol. The kinetics were determined in a static system, with the vessels deactivated with allyl bromide, and in the presence of the free radical inhibitor cyclohexene. The working temperature was 310–360 °C and the pressure was 25–85 Torr. The reaction was found to be homogeneous, unimolecular, and follows a first‐order rate law. The temperature dependence of the rate coefficients is given by the following Arrhenius equation: log k(s^?1) = [(13.82 ± 0.07) – (193.9 ± 1.0)(kJ mol^?1)](2.303RT)^?1; r = 0.9995. Theoretical calculations were carried out using density functional theory (DFT) functionals B3LYP, MPW1PW91, and PBE with the basis set 6‐31G(d,p) and 6‐31G++(d,p). The calculated values for the energy of activation and enthalpy of activation are in reasonably good agreement with the experimental values using the PBE/6‐31G (d,p) level of theory. Both experimental results and theoretical calculations suggest a molecular mechanism involving a concerted polar four‐membered cyclic transition state. The transition state structure of methanol elimination from 1,1‐dimethoxycyclohexane is characterized by a significantly elongated C? O bond, while the C_β? H bond is stretched to a smaller extent, as compared to the reactant. The process can be described as moderately asynchronic with some charge separation in the TS. Copyright © 2010 John Wiley & Sons, Ltd.     

A coupled cluster treatment of intramonomer electron correlation within symmetry-adapted perturbation theory: benchmark calculations and a comparison with a density-functional theory description

Symmetry-adapted perturbation theory (SAPT) with intramonomer electron correlation described by coupled cluster theory limited to single and double excitations was applied to 21 noncovalent complexes in their minimum geometries. The resulting benchmark contributions to the interaction energy were utilized to examine the accuracy of a more approximate variant of the SAPT method, where interacting molecules are described by density functional theory (DFT) with different functionals, like LDA, PBE, B3LYP, PBE0, M05, M05-2X, M06, and M06-2X (in all cases the asymptotic correction for the exchange-correlation functional has been utilized). Average errors for individual energy components of SAPT(DFT) are not larger than 10% for best functionals under study. Among the tested functionals PBE0, M05, and B3LYP should be especially recommended for the SAPT(DFT) approach. The M06 functional gives the largest errors with respect to SAPT(CCSD) and should not be used for describing intramonomer correlation in SAPT.