用绝热测试的压力数据计算分解反应的表观活化能

绝热量热技术的广泛应用促进了绝热动力学研究的发展.至今为止，绝热动力学的研究仍是以Townsend 和 Tou在1980年提出的绝热动力学方程为基础.他们利用反应物浓度与温度的关系建立动力学模型，从而求得了反应的活化能.但这种方法在反应体系的温度效应不明显时的应用，就受到很大的限制.因此，本文通过建立用压力表示的动力学方程，充分利用绝热量热中的压力数据，提出了一种计算反应活化能的新方法.同时用这种方法计算了过氧化二叔丁基和2-甲基-5-硝基苯磺酸的绝热分解活化能值，并与文献值和其他计算方法所得结果进行了比较.     

砷原子团簇结构的量子化学密度泛函理论研究

采用密度泛函理论的三种方法:局域自旋密度近似SVWN、梯度修正BLYP、杂化密度泛函B3LYP,优化了中性Asn、负离子Asn-(n=2～5)的结构,在优化结构基础上计算了它们的振动光谱,获得它们稳定的最低能量态的结构.其中中性Asn(n=2～5)的稳定结构的计算结果,与已有的理论结果以及实验数据进行了比较.而对负离子Asn-(n=2～5))的稳定结构作了预言.同时计算了Asn(n=2～5)的绝热电子亲和能(EAa),与有关光电子谱学的实验值符合较好.     

二甲基二硫醚对Ni/Al_2O_3催化剂加氢性能的影响

考察了二甲基二硫醚(CH3SSCH3)对Ni/Al2O3催化剂上苯、环己烯和苯乙烯加氢活性的影响,并采用BET、XRD、H2-TPR、XPS、SEM和EA等手段对催化剂进行表征。实验结果表明,在CH3SSCH3存在下,Ni/Al2O3催化剂对苯和环己烯加氢迅速失活,且环己烯加氢对CH3SSCH3的耐硫性要略强于苯加氢,而苯乙烯中共轭烯烃的加氢转化率则维持100%长时间不变。CH3SSCH3的影响顺序为芳环单烯烃共轭烯烃。此外,通过设计实验研究了CH3SSCH3对催化剂的毒化机理,发现CH3SSCH3分子首先吸附在催化剂的表面,并发生氢解生成甲烷随尾气逸出,故CH3SSCH3分子中碳对催化剂的失活影响较小,而留下的硫原子则与镍活性组分发生相互作用,毒化催化剂。     

氟代苯阳离子的理论研究

用B3LYP方法及6-311G(d，p)和6-311 G(d，p)基组，对十二种氟代苯阳离子做了理论研究，优化了它们的电子基态的结构，计算了对应分子的垂直电离势(VIP)和绝热电离势(AIP)．依据Jahn．Teller理论，计算确定了1，3，5-C6H3F^ 3和C6F^ 6离子分别具有C2v(2↑B)和D2h(2↑B2g)结构(对应分子分别为D3h和D6h结构)．其余十个离子的构型的对称点群与对应分子相同，但构型参数值有明显差别．自然布居分析计算表明这些离子的正电荷主要分布在与F原子相连的C原子和各H原子上．B3LYP／6-311 G(d，p)级别上计算的各氟代苯分子的VIP和AIP值和实验值符合得很好．     

Li+与线性碳链HC2n-体系相互作用的理论研究

应用密度泛函理论(DFT)和多体微扰理论(MP2),对Li 与线性碳链HC2n-形成的体系HC2nLi(n=1~8,C∞v)的平衡几何构型、谐振动频率、第一绝热电离能和结合能进行了研究并计算了HC2n-与Li 之间相互作用的势能曲线。计算结果表明,Li 可与线性碳链HC2n-形成稳定的化合物HC2nLi;化合物的电离能和结合能与体系大小n之间存在非线性关系。基于计算结果,得到了体系的绝热电离能和结合能与体系大小n的解析表达式。     

右旋布洛芬的低温热容

在80～370 K温度范围内, 用精密自动绝热量热计准确测量了右旋布洛芬的摩尔热容.其固态右旋布洛芬测量值对折和温度X[X=f(T)]的拟和方程为:Cp,m(S)=－39.483X4－66.649X3＋95.196X2＋210.84X＋172.98;相应的液态的拟和方程为 :Cp,m(L)=7.191X3＋4.2774X2＋56.365X＋498.5.并计算得到右旋布洛芬相对于室温(298.15 K)的摩尔焓和摩尔熵.右旋布洛芬的熔点为(324.15±0.02) K.基于摩尔热容的测量,还可获得右旋布洛芬的纯度为99.44％.并对右旋布洛芬和消旋布洛芬的热容进行了对比研究.     

非绝热聚合波的数学模拟

对非绝热波聚合体系进行了数学模拟, 将反应器尺寸与聚合波关联, 建立了非绝热条件下蔓延聚合波波速和最大波温的数学模型; 利用所建模型计算不同引发剂浓度、 不同环境温度、 不同试管直径下淀粉接枝聚合波的波速和波温, 并与实验测定值进行了比较. 结果表明, 非绝热聚合波的数学模型对无相变体系的拟合结果令人满意, 对聚合波的发生条件及反应器尺寸的影响有很好的预见性; 另一方面, 对于存在相变及与环境有质量交换的波聚合体系, 所建模型的预测性受到了限制.     

三聚氰胺的生成焓、热容和熵

通过多种热化学方法研究了三聚氰胺的热力学性质. 首先用氧弹式燃烧热量计测定了三聚氰胺在298.15 K 时的燃烧热, 根据燃烧热结果, 计算出三聚氰胺的标准摩尔燃烧焓和标准摩尔生成焓, 分别为: △cHΘm=(-2455.17±4.65) kJ·mol-1; △fHΘm =(-763.38±5.16) kJ·mol-1. 然后根据键焓与燃烧焓之间的关系, 估算出三聚氰胺中的C≈N(此键介于单键与双键之间)键能为458.30 kJ·mol-1, 此值介于碳氮单键键能和双键键能之间. 通过绝热热量计测定了三聚氰胺从80到385 K的低温热容. 根据热容值, 计算了此温度区间的标准摩尔生成焓, 其与温度呈线性关系. 另外, 三聚氰胺的热稳定性也用热重-差示扫描量热(TG-DSC)法进行了分析, 确定其分解的DSC曲线的峰顶温度为603.37 K.     

两种量热模式下物质热分解的动力学补偿效应

热分析量热仪主要包括动态、等温、恒温及绝热四种操作模式。很多学者基于动态及等温模式的测试结果,采用Arrhenius速率常数进行动力学计算,进而发现了所谓的“动力学补偿效应”。为了解绝热模式下是否也存在动力学补偿效应,分别采用绝热加速量热法(ARC)及动态差示扫描量热法(DSC)研究了过氧化二异丙苯(DCP)、40%(质量分数,下同)DCP溶液、葡萄糖、45%葡萄糖溶液的热分解特性,在此基础上基于Arrhenius公式计算了对应的表观活化能E和指前因子A,并对计算结果进行了分析。结果表明:绝热模式下,不同质量的同种样品及其溶液的最佳动力学参数,或者同一组数据采用不同的反应级数获得的lnA和E之间均存在明显的线性关系。此外,尽管由动态DSC数据计算获得的E和lnA普遍小于绝热模式的结果,但两种模式下获得的lnA和E之间仍然存在动力学补偿效应。由此可以推断,具有相同或类似反应机理的反应,虽然实验模式不同,但其E和lnA之间存在明显的动力学补偿效应。     

用绝热量热法研究自由基聚合动力学

该文用Dewar聚合釜研究均相自由基聚合动力学 ,并测定动力学和热力学有关参数 .先用动力学机理已知的甲基丙烯酸甲酯 (MMA)或苯乙烯 (St) / 2 ,2′ 偶氮二 ( 2 ,4 二甲基戊腈 ) (DAVN) /甲苯系统验证绝热量热法的可靠性 ,再用该法研究了动力学机理未知的系统 ,氯化二烯丙基二甲基铵 (DADMAC) / 2 ,2′ 偶氮双(氯化 2 丙基铵 ) /水体系 .用所测动力学和热力学参数计算而得的结果与实验值吻合良好 ,且所测MMA和St的有关参数与文献值也能吻合 .     

SiH_4与AB型卤素互化物之间的反转氢键

采用从头算方法对SiH4与AB型卤素互化物(ClF、BrF、IF、ICl、IBr、BrCl)形成的复合物的结构特征及本质进行了探讨.在MP2/6-311++G(3d,3p)水平上优化复合物的分子结构,并进行频率验证.通过分子间距离、自然键轨道(NBO)净电荷迁移数及分子图,确认SiH4与卤素互化物形成反转氢键型复合物.在MP2/6-311++G(3d,3p)水平上进行基组重叠误差(BSSE)校正后的作用能为-5.113--9.468kJ·mol-1.用对称匹配微扰理论(SAPT)对作用能进行分解,结果显示,诱导能对总吸引能的贡献在55.0%到72.2%之间,是最主要的贡献部分,静电能和色散能对总吸引能的贡献都小于25.0%.     

Multireference composite approaches for the accurate study of ground and excited electronic states: C2, N2, and O2

A multireference analog of the correlation consistent composite approach (MR-ccCA) based on complete active space with second-order perturbation theory (CASPT2) has been utilized in an investigation of the ground and valence excited states of C(2), N(2), and O(2). The performance of different second-order multireference perturbation theory methods including second-order n-electron valence state perturbation theory, second-order multireference M?ller-Plesset, and second-order generalized van Vleck perturbation theory has been analyzed as potential alternatives to CASPT2 within MR-ccCA. The MR-ccCA-P predicts spectroscopic constants with overall mean absolute deviations from experimental values of 0.0006 A?, 7.0 cm(-1), and 143 cm(-1) for equilibrium bond length (r(e)), harmonic frequency (ω(e)), and term values (T(e)), respectively, which are comparable to the predictions by more computationally costly multireference configuration interaction-based methods.     

Multireference perturbation theory with optimized partitioning. II. Applications to molecular systems

The second‐order multireference perturbation theory using an optimized partitioning, denoted as MROPT(2), is applied to calculations of various molecular properties—excitation energies, spectroscopic parameters, and potential energy curves—for five molecules: ethylene, butadiene, benzene, N₂, and O₂. The calculated results are compared with those obtained with second‐ and third‐order multireference perturbation theory using the traditional partitioning techniques. We also give results from computations using the multireference configuration interaction (MRCI) method. The presented results show very close resemblance between the new method and MRCI with renormalized Davidson correction. The accuracy of the new method is good and is comparable to that of second‐order multireference perturbation theory using Møller‐Plesset partitioning. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1390–1400, 2003     

羟烷基胺功能化离子液体吸收SO_2的量子化学计算(英文)

采用量子化学中的密度泛函理论(DFT)对羟烷基胺离子液体(HyAAILs)与二氧化硫(SO2)的相互作用进行了研究.通过几何结构优化,电荷分布和热力学参数计算等来确定离子液体中能够有效吸收SO2的官能团.HyAAILs与SO2反应形成平均距离为0.240nm的S—N键,导致电荷从ILs转移到SO2以及S—O键长和O—S—O键角的改变.气态和液态模型的计算结果表明,标准吉布斯函数变(△G苓)主要取决于阳离子的结构和分子质量.阳离子结构影响了吸收反应能垒,对于三种阳离子体系的反应活化能顺序为:Ea(secondary)Ea(tertiary)Ea(primary).理论计算结果得到了实验数据的验证,羟乙基伯胺离子液体吸收的SO2理论摩尔分数与文献中的实验数据非常接近.本研究提供了一种预测和验证功能离子液体性质的有效方法.     

He-N2O的从头算势能面及振转能级

采用超分子MP4方法和较大的基组计算得到了He-N₂O体系的分子间势能面,发现该势能面有3个极小值点,分别对应T形构型及两个线性He-ONN和He-NNO构型.同时采用离散变量表象方法预测了体系的振转能级,计算结果表明,MP4势能面支持5个振动束缚态.     

Spin component scaling in multiconfiguration perturbation theory

We investigate a term-by-term scaling of the second-order energy correction obtained by perturbation theory (PT) starting from a multiconfiguration wave function. The total second-order correction is decomposed into several terms, based on the level and the spin pattern of the excitations. To define individual terms, we extend the same spin/different spin categorization of spin component scaling in various ways. When needed, identification of the excitation level is facilitated by the pivot determinant underlying the multiconfiguration PT framework. Scaling factors are determined from the stationary condition of the total energy calculated up to order 3. The decomposition schemes are tested numerically on the example of bond dissociation profiles and energy differences. We conclude that Grimme's parameters determined for single-reference M?ller-Plesset theory may give a modest error reduction along the entire potential surface, if adopting a multireference based PT formulation. Scaling factors obtained from the stationary condition show relatively large variation with molecular geometry, at the same time they are more efficient in reducing the error when following a bond dissociation process.     

Simplified reference wave functions for multireference perturbation theory

A number of simplifications in defining the reference wave functions used in multireference second-order M?ller-Plesset perturbation theory (MRMP2) calculations are studied. The usual multiconfigurational orbital optimization is avoided by using Hartree-Fock or Kohn-Sham orbitals; the complete configuration expansion in the active-space orbitals is replaced by a severely truncated expansion, and the spin-component-scaling idea is applied to the multireference perturbation expansion. We assess these approximations to the full procedure by calculating the barrier heights for 15 processes taken from the Zhao-Gonzalez-Garcia-Truhlar database. Our results suggest that reliable and relatively cheap reference wave functions for MRMP2 calculations can be obtained from the simplifications introduced here. We hope that this will enable the application of the MRMP2 method to a larger range of chemical systems.     

Theoretical study on electronic and spin structures of [Fe2S2](2+,+) cluster: reference interaction site model self-consistent field (RISM-SCF) and multireference second-order Møller-Plesset perturbation theory (MRMP) approach

Electronic structures of [Fe(2)S(2)(SCH(3))(4)](2-,3-) in DMSO solution are calculated using reference interaction site model complete active space self-consistent field (RISM-CASSCF)/multireference second-order M?ller-Plesset perturbation theory (MRMP) method. For the reduced state, we obtain both the low-spin Fe(3+)Fe(2+) localized and high-spin Fe(2.5+)Fe(2.5+) delocalized forms, which are very close in energy. The spin interaction constants obtained from the energies of states with various spin multiplicities are in good agreement with the available experimental estimates both for the oxidized and for the reduced states. The dynamic electron correlation effect is found to be important in estimating the spin interaction between the Fe ions. The redox potentials are calculated to be 2.87 and 2.78 eV for the localized and delocalized reduced states, respectively, which are close to the experimental values. We devise a simple model for calculating the free energy curves of the reduction process based on the RISM-SCF theory. The activation barrier height is calculated to be 7.4 kcal/mol at the equilibrium geometry of oxidized state, indicating that the reduction reaction will occur efficiently in DMSO solvent. The effect of solvent fluctuation on the free energy profiles is discussed on the basis of the present calculations.     

Theoretical study of thermal rearrangements of 1-hexen-5-yne, 1,2,5-hexatriene, and 2-methylenebicyclo[2.1.0]pentane

In this research, a comprehensive theoretical investigation of the thermal rearrangements of 1-hexen-5-yne, 1,2,5-hexatriene, and 2-methylenebicyclo[2.1.0]pentane is carried out employing density functional theory (DFT) and high level ab initio methods, such as the complete active space self-consistent field (CASSCF), multireference second-order M?ller-Plesset perturbation theory (MRMP2), and coupled-cluster singles and doubles with perturbative triples [CCSD(T)]. The potential energy surface (PES) for the relevant system is explored to provide a theoretical account of pyrolysis experiments by Huntsman, Baldwin, and Roth on the target system. The rate constants and product distributions are calculated using theoretical kinetic modelings. The rate constant for each isomerization reaction is computed using the transition state theory (TST). The simultaneous first-order ordinary-differential equations are solved numerically for the relevant system to obtain time-dependent concentrations, hence the product distributions at a given temperature. Our computed energy values (reaction energies and activation parameters) are in agreement with Roth's experiments and the product distributions of Huntsman's experiments at 340 and 385 °C with various reaction times, while simulated product fractions are in qualitative accordance with Baldwin's experiment.     

Computational studies of nonstoichiometric sodium auride clusters

The molecular structures of low-lying isomers of anionic and neutral sodium auride clusters have been studied computationally at the second-order M?ller-Plesset perturbation theory level using quadruple-ζ basis sets augmented with a double set of polarization functions. The first vertical detachment energies were calculated at the M?ller-Plesset level as the energy difference between the cluster anion and the corresponding neutral cluster. The photodetachment energies of higher-lying ionization channels were calculated by adding electronic excitation energies of the neutral clusters to the first vertical detachment energy. The excitation energies were calculated at the linear response approximate coupled-cluster singles and doubles level using the anionic cluster structures. The obtained ionization energies for NaAu(-), NaAu(2)(-), NaAu(3)(-), NaAu(4)(-), Na(2)Au(2)(-), Na(2)Au(3)(-), Na(3)Au(3)(-), and Na(2)Au(4)(-) were compared to values deduced from experimental photoelectron spectra. Comparison of the calculated photoelectron spectra for a few energetically low-lying isomers shows that the energetically lowest cluster structures obtained in the calculations do not always correspond to the clusters produced experimentally. Spin-component-scaled second-order M?ller-Plesset perturbation theory calculations shift the order of the isomers such that the observed clusters more often correspond to the energetically lowest structure, whereas the spin-component-scaled approach does not improve the photodetachment energies of the sodium aurides. The potential energy surface of the sodium aurides is very soft, with several low-lying isomers requiring an accurate electron correlation treatment. The calculations show that merely the energetic criterion is not a reliable means to identify the structures of the observed sodium auride clusters; other experimental information is needed to ensure a correct assignment of the cluster structures. The cluster structures of nonstoichiometric anionic sodium aurides have been determined by comparing calculated ionization energies for low-lying structures of the anionic clusters with experimental data.