Thermal unimolecular decomposition mechanism of 2,4,6-trinitrotoluene: a first-principles DFT study

Simple C–NO₂ homolysis, 4,6-dinitroanthranil (DNAt) production by dehydration, and the nitro-nitrite rearrangement–homolysis for gas-phase TNT decomposition were recently studied by Cohen et al. (J Phys Chem A 111:11074, 2007), based on DFT calculations. Apart from those three pathways, other possible initiation processes were suggested in this study, i.e., CH₃ removal, O elimination, H escape, OH removal, HONO elimination, and nitro oxidizing adjacent backbone carbon atom. The intermediate, 3,5-dinitro-2(or 4)-methyl phenoxy, is more favor to decompose into CO and 3,5-dinitro-2(or 4)-methyl-cyclopentadienyl than to loss NO following nitro-nitrite rearrangement. Below ~1,335 K, TNT condensing to DNAt by dehydration is kinetically the most favor process, and the formations of substituted phenoxy and following cyclopentadienyl include minor contribution. Above ~1,335 K, simple C–NO₂ homolysis kinetically dominates TNT decomposition; while the secondary process changes from DNAt production to CH₃ removal above ~2,112 K; DNAt condensed from TNT by dehydration yields to that by sequential losses of OH and H above ~1,481 K and to nitro-nitrite rearrangement–fragmentation above ~1,778 K; O elimination replaces DNAt production above ~2,491 K, playing the third role in TNT decomposition; H escaping directly from TNT thrives in higher temperature (above ~2,812 K), as the fourth largest process. The kinetic analysis indicates that CH₃ removal, O elimination, and H escape paths are accessible at the suggested TNT detonation time (~100–200 fs), besides C–NO₂ homolysis. HONO elimination and nitro oxidizing adjacent backbone carbon atom paths are negligible at all temperatures. The calculations also demonstrated that some important species observed by Rogers and Dacons et al. are thermodynamically the most favor products at all temperatures, possibly stemmed from the intermediates including 4,6-dinitro-2-nitroso-benzyl alcohol, 3,5-dinitroanline, 2,6-dinitroso-4-nitro-phenylaldehyde, 3,5-dinitro-1-nitrosobenzene, 3,5-dinitroso-1-nitrobenzene, and nitrobenzene. All transition states, intermediates, and products have been indentified, the structures, vibrational frequencies, and energies of them were verified at the uB3LYP/6-311++G(d,p) level. Our calculated energies have mean unsigned errors in barrier heights of 3.4–4.2 kcal/mol (Lynch and Truhlar in J Phys Chem A 105:2936, 2001), and frequencies have the recommended scaling factors for the B3-LYP/6-311+G(d,p) method (Andersson and Uvdal in J Phys Chem A 109:2937, 2005; Merrick et al. in J Phys Chem A 111:11683, 2007). All calculations corroborate highly with the previous experimental and theoretical results, clarifying some pertinent questions.     

TNT高温热解及含碳团簇形成的反应分子动力学模拟

ReaxFF-MD模拟三硝基甲苯(TNT)高温热解显示增加了伦敦耗散力项(E_lg)的ReaxFF/lg 势函数在含能材料平衡密度计算方面具有优越性. 产物识别分析得出TNT热解的主要产物为NO₂、NO、H₂O、N₂、CO₂、CO、OH以及HONO,且最终产物为H₂O、N₂和CO₂. 使用ReaxFF势函数模拟同样过程进行比较性分析显示,在主要产物和最终产物方面与ReaxFF/lg 作用结果具有一致性,但在化学反应动力学方面表现出一些差异. ortho-NO₂键断裂和C―NO₂→C―ONO重排布-断裂形成NO₂和NO是TNT热解的主要初级反应,且前者产生速率大于后者,NO₂和NO形成后很快参与次级反应并最终形成N₂. 高温热解中形成OH等小分子会促进H₂O的形成. 环上基团相互反应或直接脱落后,主环间C―C键才发生断裂,但温度升高会加快主环断裂,并进一步分解形成CO₂,这也是高温条件下CO₂分布产生波动的一个重要原因. 并且当晶胞中的TNT分子几乎完全分解时,系统的势能开始明显衰减. 与温度相比,密度对热解中最大含碳团簇形成的影响更明显. 并且,模拟结果显示,在TNT完全分解前已经出现含碳中间体的聚合现象. 此项工作表明使用ReaxFF/lg 反应力场研究TNT高温热解可以提供具体的动力学和化学方面的信息,并有助于理解含能材料的爆轰问题并可进行安全评估.     

Application of ab initio molecular dynamics for a priori elucidation of the mechanism in unimolecular decomposition: the case of 5-nitro-2,4-dihydro-3H-1,2,4-triazol-3-one (NTO).

We have tested a new and general approach for the theoretical study of unimolecular decomposition. By combining the power of the ab initio molecular dynamics (MD) and ab initio molecular orbital (MO) methods, our approach requires no prior experimental knowledge or intuitive assumptions about the decomposition. Instead, the reaction channels are first sampled theoretically by simulating a molecule at high temperature in a number of trajectories, using the density functional theory (DFT) based ab initio MD method with a planewave basis set and pseudopotentials. Each type of these channels is then further examined by well-established ab initio MO method to locate the energy barrier and transition structure and to verify the ab initio MD results. The power of such an approach is demonstrated in a case study for the complicated unimolecular thermal decomposition of NTO (5-nitro-2,4-dihydro-3H-1,2,4-triazol-3-one), with several interesting new features uncovered. The C-NO2 homolysis is indeed the dominant channel at high temperature, while the departing NO2 could capture a H atom from the NTO ring to form HONO, by either a concerted bond breaking mechanism or by a bimolecular reaction between the NO2 group and the triazol ring. At lower temperature, the dissociation channels initiated by hydrogen migrations should be activated first. The channel with hydrogen migration followed by ring opening and then by HONO loss has an energy barrier of 38.0 kcal/mol at the rate-determining step, being the lowest among all the investigated dissociation paths and much lower than previously thought. The energy barrier for nitro-nitrite rearrangement is lower than that for the C-NO2 homolysis but makes only a minor contribution due to the entropy factor. And the NTO ring could rupture in the two C-N bonds connected to the carbonyl carbon, and the energy barriers for such processes are only 2-4 kcal/mol higher than that for the C-NO2 homolysis.     

Thermal decomposition of 2-butanol as a potential nonfossil fuel: a computational study

The thermochemistry and kinetics of the pyrolysis of 2-butanol have been conducted using ab initio methods (CBS-QB3 and CCSD(T)) and density functional theory (DFT). The enthalpies of formation and bond dissociation energies of some alcohols including 2-butanol and its derived radicals have been calculated. A variety of simple and complex dissociations have been examined. The results indicated that dehydration to 1- and 2-butene through four-center transition states is the most dominant channel at low to moderate temperatures (T ≤ 700 K), where formation of butenes is kinetically and thermodynamically more favorable than other complex and simple bond scission reactions. Although the C-C bond fission channels require more energy than needed for some complex decomposition reactions, the former pathways predominate at higher temperatures (T ≥ 800 K) due to the higher values of the pre-exponential factors. The progress of the complex decomposition reactions has been followed through intrinsic reaction coordinate (IRC) calculations to understand the mechanism of transformation of 2-butanol to different products.     

水分子对α相CL-20热分解机理影响的分子动力学研究

研究六硝基六氮杂异伍兹烷(CL-20)晶体不同晶型在不同温度下的反应机理, 对于深入认识含能材料在极端条件下的冲击起爆、冲击点火和爆轰过程等具有重要意义. 基于反应力场, 研究水分子在纯α相CL-20及其水合物的晶体结构中数量随时间的变换, 分析水分子对两种体系的初始分解和第二阶段的分解路径的影响. 计算结果表明: CL-20 分子的初始分解路径与水分子无关, 第二阶段的分解反应与水分子有关. 在低温(T<1500 K)下, 水分子对两种体系没有影响, 二者的初始分解路径均为N-NO₂键生成NO₂自由基; 在1500 K≤T≤2500 K时, 水分子作为反应物或与NO₂、、OH自由基等组成催化体系, 生成O₂、H₂O₂等产物, 加速水合物体系在高温下的第二阶段反应, 使得高温下水合物体系的化学反应速率和反应生成的NO₂自由基的数量比纯CL-20体系的化学反应速率和反应生成的NO₂自由基的数量大; 在T>2500 K时, 水分子的催化反应抑制CL-20初始分解反应, 使得在3000 K时纯CL-20体系的反应速率大于水合物体系中CL-20的反应速率.     

Thermal decomposition of pentacene oxyradicals

The energetics and kinetics of the thermal decomposition of pentacene oxyradicals were studied using a combination of ab initio electronic structure theory and energy-transfer master equation modeling. The rate coefficients of pentacene oxyradical decomposition were computed for the range of 1500-2500 K and 0.01-10 atm and found to be both temperature and pressure dependent. The computational results reveal that oxyradicals with oxygen attached to the inner rings are kinetically more stable than those with oxygen attached to the outer rings. The latter decompose to produce CO at rates comparable to those of phenoxy radical, while CO is unlikely to be produced from oxyradicals with oxygen bonded to the inner rings.     

NMR and calculational studies on the regioselective lithiation of 1-methoxynaphthalene

1-Methoxynaphthalene (1) undergoes regioselective lithiation in position 2 (n-BuLi/TMEDA) or in position 8 (t-BuLi), respectively. The detected formation of a n-BuLi/1 complex (1:1 n-BuLi/1 mixture) appears to have only minor influence on the regioselectivity (both products are obtained). The exchange of hydrogen atom H2 by deuterium results in a remarkably reduced reaction rate for the lithiation with n-BuLi in THF-d(8). This isotope effect and the formation of the thermodynamically less favorable 2-lithio compound suggest a kinetically controlled mechanism. The lack of an isotope effect for the reaction of 8-deuterio-1-methoxynaphthalene with t-BuLi and the formation of the thermodynamically preferred 8-lithiated product indicate a thermodynamically controlled mechanism. Slow conversion of the 2- into the 8-lithiated species (at higher temperatures) gives further evidence that n-BuLi and t-BuLi afford the kinetically and thermodynamically preferred products, respectively.     

单重态CCl~2与O~3反应的热力学及动力学性质研究

在量化计算的基础上运用统计热力学和Wigner校正的Eyring过渡态理论研究了不同温度下单重态CCl~2和臭氧O~3反应的热力学及动力学性质。计算结果表明该反应在低温下具有热力学优势,而在高温下具有动力学优势。     

New insights into the polymerization of methyl methacrylate initiated by rare-earth borohydride complexes: a combined experimental and computational approach

Polymerization of methyl methacrylate (MMA) initiated by the rare-earth borohydride complexes [Ln(BH(4))(3)(thf)(3)] (Ln=Nd, Sm) or [Sm(BH(4))(Cp*)(2)(thf)] (Cp*=eta-C(5)Me(5)) proceeds at ambient temperature to give rather syndiotactic poly(methyl methacrylate) (PMMA) with molar masses M(n) higher than expected and quite broad molar mass distributions, which is consistent with a poor initiation efficiency. The polymerization of MMA was investigated by performing density functional theory (DFT) calculations on an eta-C(5)H(5) model metallocene and showed that in the reaction of [Eu(BH(4))(Cp)(2)] with MMA the borate [Eu(Cp)(2){(OBH(3))(OMe)C=C(Me)(2)}] (e-2) complex, which forms via the enolate [Eu(Cp)(2){O(OMe)C=C(Me)(2)}] (e), is calculated to be exergonic and is the most likely of all of the possible products. This product is favored because the reaction that leads to the formation of carboxylate [Eu(Cp)(2){OOC-C(Me)(=CH(2))}] (f) is thermodynamically favorable, but kinetically disfavored, and both of the potential products from a Markovnikov [Eu(Cp)(2){O(OMe)C-CH(Me)(CH(2)BH(3))}] (g) or anti-Markovnikov [Eu(Cp)(2){O(OMe)C-C(Me(2))(BH(3))}] (h) hydroboration reaction are also kinetically inaccessible. Similar computational results were obtained for the reaction of [Eu(BH(4))(3)] and MMA with all of the products showing extra stabilization. The DFT calculations performed by using [Eu(Cp)(2)(H)] to model the mechanism previously reported for the polymerization of MMA initiated by [Sm(Cp*)(2)(H)](2) confirmed the favorable exergonic formation of the intermediate [Eu(Cp)(2){O(OMe)C=C(Me)(2)}] (e') as the kinetic product, this enolate species ultimately leads to the formation of PMMA as experimentally observed. Replacing H by BH(4) thus prevents the 1,4-addition of the [Eu(BH(4))(Cp)(2)] borohydride ligand to the first incoming MMA molecule and instead favors the formation of the borate complex e-2. This intermediate is the somewhat active species in the polymerization of MMA initiated by the borohydride precursors [Ln(BH(4))(3)(thf)(3)] or [Sm(BH(4))(Cp*)(2)(thf)].     

Studies on 3-Amino-4-(1H-tetrazol-5-yl)-furazan:Crystal Structure,Thermal Behavior and Energetic Performance

The structure of 4-amino-3-(5-tetrazolate)-fiirazan(HAFT) was characterized by single crystal X-ray diffraction. The thermal decomposition process of HAFT was investigated by MS-FTIR-DSC-TG coupling technique. The result shows that the exothermic process occurs from 278.7-350℃, with a peak temperature of 324.7℃. The thennal decomposition gaseous products of HAFT are NO2, CO2, HCN, CO, NH3 and H2O. The detonation velocity and detonation pressure of HAFT were calculated by the nitrogen equivalent equation. The detonation velocity of HAFT is 7727.46 m/s, which is higher than that of TNT(7178 m/s). The detonation pressure of HAFT(25.27 GPa) is satisfactory. The sensitivity tests reveal HAFT possesses excellent insensitivities to impact and friction.     

Pyrolysis of tert-butyl tert-butanethiosulfinate, t-BuS(O)St-Bu: a computational perspective of the decomposition pathways

A systematic theoretical study has been performed on the low pressure thermal decomposition pathways of t-BuS(O)St-Bu using the CCSD(T)/cc-pV(D+d)Z//B3LYP/6-311++G(2d,2p), CCSD(T)/cc-pV(D+d)Z//PBEPBE/6-311++G(2d,2p), and G3B3 level of theories. Rate constants for the unimolecular decomposition pathways are calculated using Rice?Ramsperger?Kassel?Marcus (RRKM) theory. On the basis of the experimental observation and theoretical predictions, the pyrolysis channels are considered as primary and secondary pyrolysis reactions. The primary decomposition via a five-membered transition state leads to the formation of tert-butanethiosulfoxylic acid (t-BuSSOH) and 2-methylpropene (C4H8) almost exclusively having low-pressure limit rate constant k(1)(0) = 4.67 × 10(?6)T(?4.67) exp(?11.64 kcal mol(?1)/RT) cm3 mol(?1) s(?1) (T = 500?800 K). The primary decomposition via a six-membered transition state is also identified, and that leads to the tert-butanethiosulfinic acid t-BuS(OH)S, which is the branched chain isomer of t-BuSSOH. The formation of t-BuSSOH is thermodynamically as well as kinetically favorable over t-BuS(OH)S formation, and therefore the second product could not be found experimentally. Furthermore, calculation on secondary pyrolysis pathways involving the decomposition of t-BuSSOH leads to the formation of 1-oxatrisulfane (trans-HSSOH and cis-HSSOH) and their branched isomer S(SH)OH. These three secondary product formation rates are competitive, but thermodynamics do not favor the formation of the branched isomer. Among the secondary pyrolysis products, trans-HSSOH is the most stable one, and its formation rate constant at low pressure is calculated to be k(3)(0) = 5.49 × 10(28)T(?10.70) exp(?36.22 kcal mol(?1)/RT) cm3 mol(?1) s(?1) (T = 800?1500 K). Finally, the secondary pyrolysis pathway from less stable product t-BuS(OH)S is also predicted, and that leads to trans-HSSOH and cis-HSSOH products with almost equal rates. A bond-order analysis using Wiberg bond indexes obtained by natural bond orbital (NBO) calculation predicts that the primary and secondary pyrolysis of t-BuS(O)St-Bu occur via E1-like mechanism.     

Ab initio kinetics of gas phase decomposition reactions

The thermal and kinetic aspects of gas phase decomposition reactions can be extremely complex due to a large number of parameters, a variety of possible intermediates, and an overlap in thermal decomposition traces. The experimental determination of the activation energies is particularly difficult when several possible reaction pathways coexist in the thermal decomposition. Ab initio calculations intended to provide an interpretation of the experiment are often of little help if they produce only the activation barriers and ignore the kinetics of the decomposition process. To overcome this ambiguity, a theoretical study of a complete picture of gas phase thermo-decomposition, including reaction energies, activation barriers, and reaction rates, is illustrated with the example of the β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) molecule by means of quantum-chemical calculations. We study three types of major decomposition reactions characteristic of nitramines: the HONO elimination, the NONO rearrangement, and the N-NO(2) homolysis. The reaction rates were determined using the conventional transition state theory for the HONO and NONO decompositions and the variational transition state theory for the N-NO(2) homolysis. Our calculations show that the HMX decomposition process is more complex than it was previously believed to be and is defined by a combination of reactions at any given temperature. At all temperatures, the direct N-NO(2) homolysis prevails with the activation barrier at 38.1 kcal/mol. The nitro-nitrite isomerization and the HONO elimination, with the activation barriers at 46.3 and 39.4 kcal/mol, respectively, are slow reactions at all temperatures. The obtained conclusions provide a consistent interpretation for the reported experimental data.     

Selective Insertion in Copolymerization of Ethylene and Styrene Catalyzed by Half‐Titanocene System Bearing Ketimide Ligand: A Theoretical Study

The copolymerization of ethylene and styrene can be efficiently carried out by using Cp*TiCl₂ (N = C^t Bu₂)/MAO (Cp*=η ⁵‐C₅Me₅ ) system, yielding the poly(ethylene‐co ‐styrene)s with isolated styrene units. In order to investigate the reasons for formation of the structure, the mechanism of copolymerization, especially the selective insertion of ethylene and styrene, is studied in detail by density functional theory (DFT ) method. At the initiation stage, insertion of ethylene is kinetically more favorable than insertion of styrene, and insertion of styrene kinetically and thermodynamically prefers 2,1‐insertion. That is different from the conventional half‐titanocene system, in which the 1,2‐insertion is favorable. At chain propagation stage, the computational results suggest that the continuous insertion of styrene is hard to occur at room temperature due to the high free energy barriers (28.90 and 35.04 kcal/mol for 1,2‐insertion, and 29.15 and 34.00 kcal/mol for 2,1‐insertion) and thermodynamically unfavorable factors in two different conditions. That is mainly attributed to the steric hindrance between the coming styrene and chain‐end styrene or ketimide ligand. The computational results are in good agreement with the experimental data.     

Peroxynitrate and peroxynitrite: a complete basis set investigation of similarities and differences between these NOx species

Peroxynitric acid/peroxynitrate (PNA) rivals peroxynitrous acid/peroxynitrite (PNI) in importance as a reactive oxygen species. These species possess similar two-electron oxidative behavior. On the other hand, stark differences exist in the stability of these molecules as a function of pH and in the presence of CO(2), and also in the types of bond homolysis reactions that PNA and PNI may undergo. Using CBS-QB3 theory, we examine these similarities and differences. The activation barriers for two-electron oxidation of NH(3), H(2)S, and H(2)C=CH(2) by PNA and PNI are found to be generally similar. The O-O BDE of O(2)NOOCO(2-) is predicted to be 26 kcal/mol greater than that of ONOOCO(2-). This accounts for the insensitivity of PNA to the presence of CO(2). Likewise, the O-O BDE of O(2)NOOH is predicted to be 19 kcal/mol greater than that of ONOOH, in excellent agreement with experiment. The fundamental principle underlying the large differences in O-O BDEs between PNA and PNI species is that the NO(2) that is formed from PNI can relax from the (2)B(2) excited state to the (2)A(1) ground state, whereas no such comparable state change occurs with NO(3) from PNA. Decomposition of the anions O(x)NOO(-) is more complex, with the energetics influenced by solvation. ONOO(-) can homolyze to yield NO/O(2-); however, this pathway represents a thermodynamic "dead end" since the radical pairs generated by homolysis should mostly revert to starting material. However, decomposition of O(2)NOO(-) yields the stable products NO(2-)/(3)O(2), a couple that is more stable than starting material. This may occur either by initial formation of NO(2)/O(2-) or NO(2-)/(1)O(2), with the latter intermediates thermodynamically favored both in the gas phase and in solution. Given the extremely fast back-reaction of the homolysis products, heterolysis probably dominates the observed O(2)NOO(-) decomposition kinetics. This is in agreement with the first of two "kinetically indistinguishable" mechanistic possibilities proposed for O(2)NOO(-) decomposition (Goldstein, S.; Czapski, G.; Lind, J.; Merényi, G. Inorg. Chem. 1998, 37, 3943-3947).     

Assessment of Thermal Stability and Detonation Performance of 4-Amino-1,2,4-triazolium Nitrate as compared to 2,4,6-Trinitrotoluene for Melt-cast Explosives

4-Amino-1,2,4-triazolium nitrate (4-ATN) is an energetic and non-sensitive ionic liquid, which was introduced as a good candidate in previous works for the replacement of 2,4,6-trinitrotoluene (TNT) in melt-cast explosives. Since previous studies used pure nitric acid for nitration of 4-ATN, the effect of the use of low price industrial nitric acids (50 %, 70 % and 98 %) is investigated on the percent yields of 4-ATN. The thermogravimetric and differential scanning calorimetry (TGA/DSC) are done on the synthesized 4-ATN with impure nitric acid at a heating rate of 10 K · min^–1 by the vacuum system. The obtained TGA/DSC curves confirm decomposition of 4-ATN involving melting and dissociation. Derivative thermogravimetric (DTG) curves of 4-ATN at various heating rates are applied to obtain activation energy of thermolysis by several model-free techniques. The calculated activation energies are in the range 78.7–87.7 kJ · mol^–1, which are about 10 kJ · mol^–1 more than the reported activation energy of industrial TNT (purity 98.2 %), i.e. 66–70 kJ · mol^–1. Assessments of detonation performance of 4-ATN are also compared with TNT, which show higher detonation performance of 4-ATN. Thus, 4-ATN can be used with nitramine compounds as melt-cast explosives with higher thermal stability and detonation performance than corresponding nitramine compound/TNT explosives.     

Theoretical study of the competitive decomposition and isomerization of 1-hexyl radical

The ground-state potential energy surface of the 1-hexyl system, including the main decomposition and isomerization processes, has been calculated with the MPW1K, BB1K, MPWB1K, MPW1B95, BMK, M05-2X and CBS-QB3 methods. On the basis of these data, thermal rate coefficients of different reaction channels and branching ratios were then calculated using the master equation formulation at 250–2,500 K. The results clearly point out that the 1,5 H atom transfer reaction of 1-hexyl radical with exothermicity proceeds through the lowest reaction barrier, whereas the decomposition processes are thermodynamically unfavorable with large endothermicity. The temperature effect is important on the relative importance of different reactions in the 1-hexyl system. In the low-temperature range of 250–900 K, isomerization reactions, especially 1,5 H atom transfer reaction of 1-hexyl radical, are dominating and responsible for over 82.17% of all the reactions, due to their smaller reaction barriers than those of the decomposition reactions. Furthermore, an equilibrium process involving the isomeric forms of the hexyl radicals appearing at relative low temperature was validated theoretically. However, isomerization and decomposition processes are kinetically competitive and simultaneously important under normal pyrolysis conditions.     

Bond character of thiophene on Ge(100): Effects of coverage and temperature

We have studied the adsorption and decomposition of thiophene (C4H4S) on Ge(100) using scanning tunneling microscopy (STM), high-resolution core-level photoemission spectroscopy (HRPES), and density functional theory (DFT) calculation. Analysis of S 2p core-level spectra reveals three adsorption geometries, which we assign to a Ge-S dative bonding state, a [4 + 2] cycloaddition bonding state, and a decomposed bonding state (desulfurization reaction product). Furthermore, we found that the number ratio of the three adsorption geometries depended on the molecular coverage and the annealing temperature. At low coverages, the kinetically favorable dative bonding state is initially formed at room temperature. As the molecular coverage increases, thermodynamically stable [4 + 2] cycloaddition reaction products are additionally produced. In addition, we found that as the surface temperature increased, the [4 + 2] cycloaddition reaction product either possibly desorbed as molecular thiophene or decomposed to form a metallocycle-like species (C4H4Ge2) and a sulfide (Ge2S). We systematically elucidate the changes in the bonding states of adsorbed thiophene on Ge(100) according to the thiophene coverage and annealing temperature.     

Thermodynamic and Kinetic Studies on the SiH + XH_3 (X=N, P) Reactions

Silicon chemistry has been attracted more attention because of its applications to the production of thin silicon films and the etching of silicon wafers in micro-electronics1. Silylidyne (SiH), which plays an important role in plasma chemical vapor deposition (CVD) processes, has been investigated in experimental research2-5. The reaction mechanism of SiH insertion reaction with small molecules such as XH3 (X=N, P) was recently studied by means of M豯ler-Plesset perturbation theory6.…     

活性炭纤维催化氧化NO的量子化学研究

利用密度泛函理论的B3LYP方法,6-31G(d)基组,在zigzag型的四并苯模型上对NO、O2分子在活性炭纤维(ACFs)表面的吸附行为进行研究,并探讨了ACFs催化氧化NO的主要机理路径。研究结果表明,环境气氛中的O2分子可以先吸附于ACFs表面形成两个半醌基(C-O),之后C-O和吸附态的NO(C-NO)发生氧化反应生成-NO2;游离态的O2也可以经过ACFs表面的催化作用形成活性氧原子(O*)从而直接和吸附态的NO反应生成-NO2。与NO相比,O2分子的吸附能大,在同NO的竞争吸附中占据优势,结合统计热力学分析,吸附态的NO和游离态的O2所产生的活性氧原子发生氧化反应是NO转化为NO2的主要途径。     

Decomposition reaction rate of BCl3-C3H6(propene)-H2 in the gas phase

The decomposition reaction rate in the BCl(3)-C(3)H(6)-H(2) gas phase reaction system in preparing boron carbides was investigated based on the most favorable reaction pathways proposed by Jiang et al. [Theor. Chem. Accs. 2010, 127, 519] and Yang et al. [J. Theor. Comput. Chem. 2012, 11, 53]. The rate constants of all the elementary reactions were evaluated with the variational transition state theory. The vibrational frequencies for the stationary points as well as the selected points along the minimum energy paths (MEPs) were calculated with density functional theory at the B3PW91/6-311G(d,p) level and the energies were refined with the accurate model chemistry method G3(MP2). For the elementary reaction associated with a transition state, the MEP was obtained with the intrinsic reaction coordinates, while for the elementary reaction without transition state, the relaxed potential energy surface scan was employed to obtain the MEP. The rate constants were calculated for temperatures within 200-2000 K and fitted into three-parameter Arrhenius expressions. The reaction rates were investigated by using the COMSOL software to solve numerically the coupled differential rate equations. The results show that the reactions are, consistent with the experiments, appropriate at 1100-1500 K with the reaction time of 30 s for 1100 K, 1.5 s for 1200 K, 0.12 s for 1300 K, 0.011 s for 1400 K, or 0.001 s for 1500 K, for propene being almost completely consumed. The completely dissociated species, boron carbides C(3)B, C(2)B, and CB, have very low concentrations, and C(3)B is the main product at higher temperatures, while C(2)B is the main product at lower temperatures. For the reaction time 1 s, all these concentrations approach into a nearly constant. The maximum value (in mol/m(3)) is for the highest temperature 1500 K with the orders of -13, -17, and -23 for C(3)B, C(2)B, and CB, respectively. It was also found that the logarithm of the overall reaction rate and reciprocal temperature have an excellent linear relationship within 700-2000 K with a correlation coefficient of 0.99996. This corresponds to an apparent activation energy 337.0 kJ/mol, which is comparable with the energy barrier 362.6 kJ/mol of the rate control reaction at 0 K but is higher than either of the experiments 208.7 kJ/mol or the Gibbs free energy barrier 226.2 kJ/mol at 1200 K.