CH3OCF2CF2OCH3+Cl的反应机理及动力学性质

利用密度泛函理论直接动力学方法研究了反应CH₃OCF₂CF₂OCH₃+Cl的微观机理和动力学性质. 在BB1K/6-31+G(d,p)水平上获得了反应的势能面信息, 计算中考虑了反应物CH₃OCF₂CF₂OCH₃两个稳定构象(SC1和SC2)的氢提取通道和取代反应通道. 利用改进的正则变分过渡态理论结合小曲率隧道效应(ICVT/SCT)计算了各氢提取通道的速率常数, 进而根据Boltzmann配分函数得到总包反应速率常数(k_T)以及每个构象对总反应的贡献. 结果表明296 K温度下计算的k_T(ICVT/SCT)值与已有实验值符合得很好. 由于缺乏其他温度速率常数的实验数据, 我们预测了该反应在200-2000 K温度区间内反应速率常数的三参数表达式: k_T=0.40×10^-14T^1.05exp(-206.16/T).     

N-异丙基丙烯酰胺温度敏感性水凝胶相转变动力学研究及其应用

 研究指出表观二级动力学方程可以很好地描述N-异丙基丙烯酰胺水凝胶的溶胀和消溶胀动力学.即溶胀动力学方程为dR/dt=k₁(R_e-R)²,消溶胀动力学方程为-dR/dt=k_c(R-R_e)².把这种水凝胶用于分离高分子水溶液时可引入“单位溶张比分离循环的合理时耗”这样一个参量.它根据溶胀和消溶胀过程中的起始溶胀比、平衡溶胀比、表观溶胀动力学常数和表观消溶胀动力学常数求出.具体公式为△t₁(T_s,T_c)=2/[R_c(T_s)-R₀(T_s)]²k_s(T_s)+15/[R₀(T_c)- R_s(T_c)]²k_c(T_c)^1/2在理想情况下,分离过程的“总合理时耗”与△t_1成正比,比例系数为分离过程中的除水总量与干凝胶用量的比值,即△t_r=W_W/WG·△t₁.当根据二个动力学方程求得的总时耗计算值处于(0.9△t_r,1.1△t_r)范围内时,表明所选干凝胶用量和循环溶胀比区段均合适.     

层状钙钛矿化合物 (n-CnH2n+1NH3)2ZnCl4非等温固-固相变动力学的确定

确定了具有层状钙钛矿结构,脂肪链C原子数从12到18四氯合锌酸二烷基铵(C_nH_2n+1NH₃)₂ZnCl₄系列化合物的非等温固-固相变动力学。采用Kissinger和Ozawa两种动力学模型,对不同温度下测定的每个样品DSC热分析曲线进行数据处理,计算固-固相变过程的活化能E_a。实验结果表明,随着C原子数的增加,脂肪链的有序度和刚度降低,导致固-固相变活化能E_a随C原子数的增大而降低。两种模型方法的计算活化能E_a的结果相一致。并且各个化合物的固-固相变反应级数不随升温速率和烷基链长变化而变化,总保持为1。     

1,3,3-三硝基氮杂环丁烷的热安全性

借助不同加热速率(β)的非等温DSC曲线离开基线的初始温度(T₀)、onset温度(T_e)和峰顶温度(T_p), Kissinger法和Ozawa法求得的热分解反应的表观活化能(E_k和E_O)和指前因子(A_k), Hu-Zhao-Gao方程ln β_i＝ln[A₀/(b_{e0 or p0}G(α))]＋   b_{e0 or p0}T_{ei or pi}求得的b_{e0 or p0}, Zhao-Hu-Gao方程ln β_i＝ln[A₀/((a_{e0 or p0}＋1)G(α))]＋(a_{e0 or p0}＋1) ln T_{ei or pi}求得的a_{e0 or p0}, 微热量法确定的比热容(C_p), 以及密度(ρ)、热导率(λ)和分解热(Q_d, 取爆热之半)数据, Zhang-Hu-Xie-Li公式、Hu-Yang-Liang-Xie公式、Hu-Zhao-Gao公式、Zhao-Hu-Gao公式、Smith方程、Friedman公式和Bruckman-Guillet公式, 计算了TNAZ在β→0时的T₀, T_e和T_p值(T₀₀, T_eo和T_p0)、热爆炸临界温度(T_be和T_bp)、绝热至爆时间(t_Tlad)、撞击感度50%落高(H₅₀)和热点起爆临界温度(T_cr), 得到了评价TNAZ热安全性的结果: T_SADT＝T_e0＝485.81 K, T_p0＝497.38 K, T_beo＝499.50 K, T_bp0＝513.45 K, t_Tlad＝8.90 s (n＝0), t_Tlad＝8.96 s (n＝1), t_Tlad＝9.01 s (n＝2), H₅₀＝28.88 cm, T_cr＝641.46 K (T_room＝293.15 K), T_cr＝658.89 K (T_room＝300 K), 表明: (1) TNAZ对热是稳定的; (2)撞击感度好于环三亚甲基三硝胺(RDX); (3)热点起爆临界温度高于RDX, 而界于1,3,5-三氨基-2,4,6-三硝基苯(TATB)和六硝基茋(HNS)之间.     

蝎型铜配合物：Cu2[ μ-pz]2[HB(pz)3]2和Cu[B(pz)4]2的合成、结构及热分解动力学性质

本文设计合成了两种以聚吡唑硼酸盐、吡唑为配体的铜配合物Cu₂[ μ-pz]₂[HB(pz)₃]₂(1)和Cu[B(pz)₄]₂(2)(pz：吡唑(C₃H₄N₂))。运用元素分析、红外光谱对配合物进行了表征，并用X-ray衍射测定了它们的晶体结构。非等温热分解动力学研究表明：配合物1的热分解反应分两步，配合物2的热分解反应一步进行。通过计算，配合物1热分解的第一步反应的可能机理为成核与生长，n=1/4；第二步反应的可能机理为化学反应。其非等温动力学方程分别为：dα/dT=A/β e^-E/RT·1/4(1-α)[-ln(1-α)]^-3和dα/dT=A/β e^-E/RT·(1-α)²。分解反应的表观活化能分别是520.37 kJ·mol^-1和149.65 kJ·mol^-1；指前因子lnA分别是118.06 s^-1和28.10 s^-1。配合物2热分解的可能机理为化学反应。其非等温动力学方程为：dα/dT=A/β e^-E/RT·(1-α)²。分解反应的表观活化能是111.41 kJ·mol^-1；指前因子lnA是21.20 s^-1。     

超分子化合物 [Eu(C10H9N2O4)(C10H8N2O4)(H2O)3]2·phen·4H2O 的热分解非等温动力学

用非等温热重法研究了Eu(III)的化合物[Eu(C₁₀H₉N₂O₄)(C₁₀H₈N₂O₄)(H2O)₃]₂·phen·4H₂O 的热分解及其动力学，并用Kissinger和Ozawa 方法计算了第一分解阶段的活化能.     

高聚物玻璃化转变温度和分子参数的关系 Ⅰ.高聚物链柔性与Tg的关系

 本文采用晶格模型,以动力学链段长度作为统计单元大小,推导了高聚物玻璃化温度T_8和链静态刚性因子σ²(T₈),链动态刚性因子β(T₈)以及聚合度DP等分子参数之间的关系。具体讨论了链柔性对T₈的影响。理论预测和几十种聚合物的实验数据能较好吻合,分析结果表明T₈值基本上取决于高聚物链σ(T₈)大小。     

BTATz-Pb复合物对双基和RDX-改性双基推进剂的热行为、非等温动力学及燃烧性能的影响

制备了含3,6-双(1-氢-1,2,3,4-四唑-5-氨基)-1,2,4,5-四嗪(BTATz)铅复合物(LCBTATZ)的双基推进剂和改性双基推进剂. 采用热重-微商热重法(TG-DTG)及差示扫描量热法(DSC)研究了其热分解行为和非等温分解动力学并在此基础上评价了其热安全性. 结果表明, LCBTATz-DB复合物中在350-540 K之间只存在一个放热分解峰, LCBTATz-CMDB复合物中存在两个连续的放热分解峰在390-540 K温度范围内, 其机理方程分别为: f(α)=α^-1/2和f(α)=2(1-α)^3/2. 计算了热加速分解温度(T_SADT)、热爆炸临界温度(T_b)、热点火温度(T_TIT)和绝热至爆时间(t_Tlad),其值分别为: DB001复合物T_SADT=444.50 K, T_TITT=453.96 K, T_b=471.84 K; t_Tlad=39.36 s; CMDB100复合物, T_SADT=442.38 K, T_TITT=452.89 K,T_b=464.13 K,t_Tlad=21.3 s,并以此来评价化合物的热安全性. 考察了LCBTATz-DB以及LCBTATz-CMDB的燃烧性能, 结果表明LCBTATZ 是一种高效的双基燃烧催化剂, 在较大的压力范围内可以显著的提高燃速并且大幅度的降低压力指数. 对于双基推进剂在2-8 MPa压力范围内出现了明显的超燃速现象, 8-12 MPa出现了“麦撒”效应, 对于改性双基推进剂的压力指数降到0.18.     

纳米结晶体Ni0.5Zn0.5Fe2O4对高氯酸铵热行为及分解反应动力学的影响

采用差示扫描量热法(DSC)、热重和微分热重(TG-DTG)及固相原位反应池/快速扫描傅立叶变换红外联用技术(hyphenated in situ thermolysis/RSFTIR)研究了纳米结晶体Ni_0.5Zn_0.5Fe₂O₄与高氯酸铵(AP)组成的混合物的热行为和分解反应动力学。结果表明：Ni_0.5Zn_0.5Fe₂O₄使得AP的低、高温分解放热峰温分别提前17.44 K和27.74 K,并使得对应的分解热分别增加3.7 J·g^-1和193.7 J·g^-1。Ni_0.5Zn_0.5Fe₂O₄并不影响AP的晶转温度和晶转热。Ni_0.5Zn_0.5Fe₂O₄使得AP的TG曲线出现3个阶段,并使得后2个失重阶段的初始和终止温度都有所提前。凝聚相分解产物分析表明Ni_0.5Zn_0.5Fe₂O₄加速了凝聚相AP的分解及氨气的释放。含Ni_0.5Zn_0.5Fe₂O₄的AP的高温分解反应的动力学参数E_a=238.88 kJ·mol^-1,A=1018.59 s^-1,动力学方程可表示为dα/dt=10^18.99(1-α)[-ln(1-α)]^3/5e^-2.87×104T。始点温度(T_e)和峰顶温度(T_p)计算得出AP的热爆炸临界温度值分别为:574.83 K和595.41 K。分解反应的活化熵(ΔS^≠)、活化焓(ΔH^≠)和活化能(ΔG^≠)分别为:109.61 J·mol^-1·K^-1、236.49 kJ·mol^-1及172.58 kJ·mol^-1。     

从DSC曲线数据计算/确定含能材料自催化分解反应动力学参数和热爆炸临界温升速率的方法

为应用热爆炸临界温升速率(dT/dt)_Tb评价含能材料(EMs)的热安全性, 得到计算(dT/dt)_Tb值的基本数据, 用合理的假设, 由Semenov的热爆炸理论和9 个自催化反应速率方程[dα/dt=Aexp(-E/RT)α(1-α) (I), dα/dt=Aexp(-E/RT)(1-α)ⁿ(1+K_catα) (II), dα/dt=Aexp(-E/RT)[α^a-(1-α)ⁿ)] (III), dα/dt=A₁exp(-E_a1/RT)(1-α)+A₂exp(-E_a2/RT)α(1-α) (IV), dα/dt=A₁exp(-E_a1/RT)(1-α)^m+A₂exp(-E_a2/RT)αⁿ(1-α)^p (V), dα/dt=Aexp(-E/RT)(1-α) (VI), dα/dt=Aexp(-E/RT)(1-α)ⁿ (VII), dα/dt=A₁exp(-E_a1/RT)+A₂exp(-E_a2/RT)(1-α) (VII), dα/dt=A₁exp(-E_a1/RT)+A₂exp(-E_a2/RT)α(1-α) (IX)]导出了计算(dT/dt)_Tb值的9 个表达式. 提出了从不同恒速升温速率(β)条件下的差示扫描量热(DSC)曲线数据计算/确定EMs自催化分解反应的动力学参数和自催化分解转向热爆炸时的(dT/dt)_Tb的方法. 由DSC曲线数据的分析得到了用于计算(dT/dt)_Tb值的β→0 时的onset 温度(T_e0),热爆炸临界温度(T_b)和相应于T_b时的转化率(α_b). 分别用线性最小二乘法和信赖域方法得到方程(I)和(VI)及方程(II)-(V)和方程(VII)-(IX)中的自催化分解反应动力学参数. 用上述基础数据得到了EMs的(dT/dt)_Tb值. 结果表明: (1) 在非等温DSC条件下硝化棉(NC, 13.54% N)分解反应可用表观经验级数自催化反应速率方程dα/dt=10^15.82exp(-170020/RT)(1-α)^1.11+10^15.82exp(-157140/RT)α^1.51(1-α)^2.51描述; (2) NC (13.54% N)自催化分解转向热爆炸时的(dT/dt)_Tb值为0.103 K·s^-1.     

Deviation of activation energy caused by neglecting a temperature term in Ozawa Equation

In the non-isothermal kinetic equation, the temperature integral from 0 to T ₀ never equals to zero, so a deviation of activation energy is introduced when the term is neglected. We propose a new evaluation method on Ozawa Equation. The results show that neglecting this term in Ozawa Equation leads to large activation energy errors, when small absolute values of Δx = x − x ₀(x = E/RT). The deviation of the activation energy is less than 0.55% when |Δx| > 5, but may exceed 10% if |Δx| < 2. The application of this method on the dehydration of CaC₂O₄· H₂O shows that the theoretical error was close to the actual error.     

Non-isothermal kinetics in solids

The integral methods, which are obtained from the various approximations for the temperature integral, have been extensively used in the non-isothermal kinetic analysis. In order to obtain the precision of the integral methods for the determination of the activation energy, several authors have calculated the relative errors of the activation energy obtained from the integral methods. However, in their calculations, the temperature integral at the starting temperature was neglected. In this work, we have performed a systematic analysis of the precision of the activation energy calculated by the integral methods without doing any simplifications. The results have shown that the relative error involved in the activation energy determined from the integral methods depends on two dimensionless quantities: the normalized temperature θ=T/T ₀, and the dimensionless activation energy x ₀=E/RT ₀ (where E is the activation energy, T is the temperature, T ₀ is the starting temperature, R is the gas constant).     

The correctional kinetic equation for the peak temperature in the differential thermal analysis

The peak temperature (T _p) and different temperature (ΔT) are the basic information in the differential thermal analysis (DTA). Considering the kinetic relation and the heat equilibrium in DTA, a correctional differential kinetic equation (containing T _p and ΔT parameter) is proposed. In the dehydration reaction of CaC₂O₄·H₂O, the activation energy calculated from the new equation showed some smaller than that from Kissinger equation, but some bigger than that from Piloyan equation.     

Thermogravimetric investigation of the pyrolysis of pitch materials. A compensation effect and variation in kinetic parameters with heating rate

Thermogravimetric measurements of weight loss accompanying the pyrolysis of four pitches have been made over a range of linear heating rates. For three of the samples, the data at each heating rate could be described by an integral and a differential method of analysis, assuming a simple order function for f(α), with the result that the apparent activation energy increased with heating rate. The data for all four samples could also be satisfactorily described by the Ozawa or Friedman multiple heating rate methods, and these resulted in apparent activation energies (E_a) which increased with the value of β at which they were determined. It is suggested that this tendency for the apparent activation energy to increase, as the temperature is raised, is due to a change in the relative importance of the different reactions which lead to weight loss in this system. The apparent kinetic parameters all fall on a common compensation plot which is used to explain the relative magnitude of E_a values from Ozawa and Doyle methods of analysis. The higher values of E_a from Friedman than from Ozawa analyses are also explained.     

Kinetic analysis of solid‐state reactions: Precision of the activation energy calculated by integral methods

The integral methods are extensively used for the kinetic analysis of solid‐state reactions. As the Arrhenius integral function [p(x)] does not have an exact analytical solution, different approximated equations have been proposed in the literature for performing the kinetic analysis of experimental integral data. Since the first approximation of Van Krevelen, a large number of equations have been proposed with the objective of increasing the precision in the determination of the Arrhenius integral, as checked from the standard deviation of the approximated function with regard to the real exact value of the integral. However, the main application of these equations is the determination of the kinetic parameters, in particular activation energies, and not the computation of the Arrhenius integral. A systematic analysis of the errors involved in the determination of the activation energy from these integral methods is still missing. A comparative study of the precision of the activation energy as a function of x and T computed from the different integral methods has been carried out. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 658–666, 2005     

Study on kinetics of combustion of brick-shaped carbonaceous materials

Combustion of brick-shaped carbonaceous materials (carbon deposits from coke oven, coke and electrographite) was carried out in thermobalance in static air. Analysis of kinetics of the process was carried out using both classical (Arrhenius law) and newer (three-parametric equation) methods. In classical approach two types of kinetic equations were used in calculations: differential and integral. The results obtained show that, independently on kinetic variables (α – conversion degree or m – mass of sample) used in differential equations, kinetics of combustion of brick-shaped carbonaceous materials is characterized by only one pair of Arrhenius coefficients: activation energy (E) and pre-exponential constant (A). At the same time the integral equation demonstrates distinction in relation to methods based on differential equations, generating higher activation energies and separate isokinetic effect (IE). Parallel IE shows that kinetic analysis has to encompass activation energy in connection to second coefficient, pre-exponential constant A, depending on assumptions made for kinetic equations. On the other hand three-parametric equation allows describing kinetic of combustion in alternative way using only one experimental value – initial temperature in form of point of initial oxidation (PIO) – and also offers new methods of interpretation of the process.     

Thermal analysis and gelation property of multifunctional epoxy resins used for pultruded composite ropes

Non-isothermal curing reactions of three different multifunctional epoxy resin systems were investigated by differential scanning calorimetry. The Kissinger equation was applied to calculate the apparent activation energy, and the Levenberg–Marquardt algorithm was used to fit the curing kinetic data. It was observed that the two-parameter model was in good match with the curing kinetics. In addition, dynamic mechanical thermal analysis was used to obtain the glass transition temperature (T _g). Furthermore, the thermal stabilities of the systems were studied by thermogravimetric (TG) analysis, the integral procedure decomposition temperature and temperature index T _s were used to characterize the thermal stability. Finally, the gelation time was measured by plate–knife method of a home-made device, and the relationship between gelation time and temperature was established, according to which the pultrusion process parameters were predicted.     

Crystallization kinetics of Ti20Zr20Cu60 metallic glass by isoconversional methods using modulated differential scanning calorimetry

The study of crystallization kinetics of amorphous alloys has been a matter of great interest for material researchers for past few decades, since it provides information about the kinetic parameters i.e., activation energy of crystallization and the frequency factor. These kinetic parameters can be calculated by model-free isoconversional methods. Isoconversional methods allow calculating the activation energy as a function of degree of conversion, α. Hence, these methods provide accurate results for multistep processes like crystallization. Model-free methods are categorized as linear and non-linear isoconversional methods. Linear methods are further classified as linear differential and linear integral isoconversional methods. In present work, we have used these isoconversional methods to study the effect of non-linear heating rate, employed by modulated differential scanning calorimetry (MDSC), on the non-isothermal crystallization kinetics of Ti₂₀Zr₂₀Cu₆₀ metallic glass. For Ti₂₀Zr₂₀Cu₆₀, MDSC curves clearly indicate a two-step crystallization process. Both crystallization peaks were studied based on the modified expressions for isoconversional methods by non-linear heating rate. The term corresponding to non-linearity comes out to be (A _T ω/2β)². The effect of non-linear heating rate on measurement of kinetic parameters by isoconversional methods is studied. The activation energy of crystallization is calculated for Ti₂₀Zr₂₀Cu₆₀ metallic glass for various degrees of conversion by linear integral isoconversional methods i.e., Ozawa–Flynn–Wall, Kissinger–Akahira–Sunose, and also with Friedman method which is a linear differential isoconversional method.     

Thermal behavior and kinetic modeling of (NH4)4UO2(CO3)3 decomposition under non-isothermal conditions

The thermal behavior and kinetic analysis of ammonium uranyl carbonate decomposition has been studied in inert gas, O₂, and 90%Ar–10%H₂ atmospheres under non-isothermal conditions. The results showed a dependence on specific surface area with the decomposition temperature of ammonium uranyl tri-carbonate (AUC). Specific surface area increases and reaches a maximum between 300 and 400 °C and decreases at T > 400 °C. The reaction paths of AUC decomposition under the three atmospheres were proposed. The integral methods Flynn–Wall–Ozawa (FWO) and Kissinger–Akahira–Sunose (KAS) were used for the kinetic analysis. The activation energy averages are 58.01 and 56.19 kJ/mol by KAS and FWO methods, respectively.

     

Kinetic parameters determination in non-isothermal conditions for the crystallisation of a silica-soda-lead glass

Two integral isoconversional methods (Flynn–Wall–Ozawa and Kissinger–Akahira–Sunose) and the invariant kinetic parameters method (IKP) were used in order to examine the kinetics of the non-isothermal crystallisation of a silica-soda-lead glass. The objective of the paper is to show the usefulness of the IKP method to determine both the activation parameters and the kinetic model of the investigated process. Thismethod associated with the criterion of coincidence of kinetic parameters for all heating rates and some procedures of the evaluation of the parameter from Johnson–Mehl–Avrami–Erofeev–Kolmogorov (JMAEK) equation led us to the following kinetic triplet: activation energy, E=170.5±2.5 kJ mol^–1 , pre-exponential factor, A=1.178±0.350·10 10 min^–1 and JMAEK model (A _m) m=1.5.