硝酸酯的热分解性能

分别在B3LYP/6-31G*和MP2/6-31G*的理论水平下,计算得到硝酸正丙酯、硝酸异丙酯、硝酸异辛酯和二缩三乙二醇二硝酸酯4种硝酸酯的O-NO2键离解能(BDE)。 采用常压DSC和高压DSC实验方法,研究了4种硝酸酯的热分解过程,获得热分解反应动力学参数。 常压下4种物质的热分解反应发生在气相区,当压力增大至2 MPa时,4种物质直接发生液相分解反应。 4种硝酸酯的O-NO2键离解能在很大程度上符合由实验分析得到的活化能,表明4种硝酸酯的热分解反应只是单分子O-NO2键的均裂反应。     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

Al6(OH)18(H2O)6的结构及成键方位的从头算及密度泛函分析Ⅱ

用RHF法6-31G、6-31G*、6-31G**水平的量子化学从头算及RB3LYP/6-31G、B3LYP/6-31G*,B3LYP/6-31G**水平的密度泛函方法,采用SCRF=Dipole溶剂模型,计算了三水铝石有利生长基元最优结构Al6(OH)18(H2O)6的总能量、布居数、原子静电荷等.计算结果表明,桥联OH基团更易成键,Al6(OH)18(H2O)6较为有利的成键方位是桥联OH基团方位.     

特殊氢方法预测丙氨酸-α-四肽构象稳定性

使用B3LYP/6-31G*方法优化得到了丙氨酸-α-四肽分子的48个稳定构象.定义了丙氨酸-α-多肽中的特殊氢原子,对构象中与特殊氢原子有关的主要非键相互作用进行分析,提出了使用与特殊氢原子有关的非键相互作用来预测多肽构象的相对稳定性,并将之称为特殊氢方法.基于丙氨酸-α-二肽分子的5个构象和三肽分子的7个构象,确定了特殊氢方法中的与特殊氢原子有关的非键相互作用参数.利用特殊氢方法预测丙氨酸-α-四肽分子的48个构象的相对稳定性,与B3LYP/6-31G*方法的相对能量比较,得到了满意的结果.特殊氢方法得到的相对能量(Y)和B3LYP/6-31G*方法得到的相对能量(X)的线性相关方程为Y=0.9296X 2.2041,相关系数R=0.9532,标准偏差为3.0kJ/mol,偏差绝对值≤4.18kJ/mol的构象占87.5%.     

硝酸酯增塑剂力学性能和界面相互作用的分子动力学模拟

运用分子动力学(MD)方法, 模拟研究了硝化甘油(NG)及其与硝化三乙二醇(TEGDN)组成的硝酸酯增塑剂的低温力学性能. 结果表明, NG/TEGDN混合体系较NG单组分体系的刚性减弱, 延展性和各向同性增强. 结合能计算和径向分布函数分析揭示了混合型硝酸酯增塑剂组分之间的相互作用及其本质.     

1, 2-二硒方酸气相酸性和芳香性的理论研究

在RHF/6-311G^*^*水平优化得到1,2-二硒方酸(3,4-二羟基-3-环丁烯-1,2-二硒酮)三种平面构象异构体的平衡几何构型。进一步用MP2(full)/6-311G^*//RHF/6-311G^*^*方法计算三种异构体的单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近。用优化的最稳定构象ZZ型异构体在RHF/6-311G^*^*//RHF/6-311G^*^*,RHF/6-311+G^*^*//RHF/6-311+G^*^*,MP2(full)/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其气相酸性[ΔGⅲ~(~2~9~8~K~)]和同键反应芳香性稳定化能(HASE)。用基团加和法(groupincrementapproach)在RHF/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其磁化率增量(Λ)。计算结果指出标题化合物的键长发生了平均化,同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

六硝基六氮杂异伍兹烷结构和性质的理论研究

用abinitio和DFT方法,分别在HF/6-31G^*和B3LYP/6-31G^*水平下全优化计算了高能量密度材料六硝基六氮杂异伍兹烷(HNIW)的α(γ),β和ε型构象的分子几何构型、电子结构、IR谱和298～1000K温度下的热力学性质,细致分析比较了两种方法和相关的实验结果。理论计算几何参数与实验值相一致。分子中N—N键较长,N—N键Mulliken集居数较小,预示该键为热解和起爆的引发键。所得的IR谱形符合实验、指纹区频率与实验的平均绝对差值小于45cm^-1。由前线MO能级及其差值预示的热力学稳定性次序[ε＞α(γ)＞β]与实验排序相吻合。     

电场作用下分子导线的理论研究

摘要利用从头计算法分别在HF/6-31G, HF/6-31G*, HF/6-31G**, HF/6-31+G, HF/6-31++G, HF/6-31+G*, HF/6-31+G**, HF/D95+*, B3LYP/6-31G*和B3LYP/6-31+G*水平上计算了5个单体的聚乙炔分子导线, 从几何构型、 SCF能量和分子轨道能级三个方面讨论了外电场对分子导线的影响, 给出了聚乙炔分子导线性质与外电场变化的定量关系.     

Si2Br6的分子振动光谱的理论研究

用量化从头算方法(HF/6-31G^*)和密度泛函方法(B3LYP/6-31G^*)以6-31G标准基组加一个极化函数,对Si₂Br₆分子的平衡几何构型和振动频率分别进行优化和计算,优化的结果与实验结果吻合得较好.按照Pulay的建议对HF/6-31G^*水平上所计算的谐性力场进行标度(标度因子取0.9).用HF/6-31G^*SQM力场所计算的基频预测值和实验值的平均误差为9.4cm^-1,最大误差为23.6cm^-1;用B₃LYP/6-31G^*未标度力场所计算的基频预测值和实验值的平均误差为8.6cm^-1,最大误差为16.6cm^-1;用该密度泛函方法所计算的基频预测值比用HF/6-31G^*的标度后的SQM力场所计算的基频预测值和实验值(除Si-Si键扭转振动基频之外的11条振动基频)吻合得更好.HF/6-31G^*和B3LYP/6-31G^*计算给出Si-Si键扭转振动基频的预测值分别为14cm^-1和9cm^-1.     

A theoretical study of explosive sensitizers

Quantum chemical calculations at the HF/6-31G^* and B3LYP/6-31G^* levels have been performed on five explosive sensitizers, ethyl nitrate (EN), n-propyl nitrate (NPN), isopropyl nitrate (IPN), 2-ethylhexyl nitrate (EHN) and tetraethylene glycol dinitrate (TEGDN). Theoretical study has made a detailed molecular-level investigation of the title compounds. Based on the Mulliken populations and bond lengths, the fission of the O₂–N₃ can be acceptable reasonably. Charge distribution analysis indicates that the five nitrates produce NO₂ gas during the dissociation of the O₂–N₃ weak bond. We also order the relative thermal stability of five nitrates on the basis of frontier orbital energy (E _HOMO, E _LUMO) and energy gap (ΔE = E _HOMO − E _LUMO).     

A theoretical study of five nitrates: Electronic structure and bond dissociation energies

Three density-functional methods (B3P86, B3PW91, and B3LYP) are employed to investigate the O–NO₂ bond lengths, frontier orbital energies, and O–NO₂ bond dissociation energies (BDEs) of n-propyl nitrate (NPN), isopropyl nitrate (IPN), 2-ethylhexyl nitrate (EHN), triethylene glycol dinitrate (Tri-EGDN), and tetraethylene glycol dinitrate (Tetra-EGDN). It is found that the O–NO₂ bond lengthens (destabilizes) in the order of IPN, NPN, EHN, Tetra-EGDN, and Tri-EGDN. From the data of frontier orbital energies (E_HOMO, E_LUMO), and energy gaps (ΔE), we estimate the relative thermal stability ordering of five nitrates and their corresponding radicals. The predicted BDEs of O–NO₂ bond in NPN, IPN, EHN, Tri-EGDN, and Tetra-EGDN, are 176.6, 174.5, 168.1, 156.1, and 159.3 kJ mol⁻¹, respectively. Based on the finding that the present results of BDEs are well coincident with the experimental results of apparent activation energies from the literature, we can draw a conclusion that the experimental thermolysis of five nitrates is only unimolecular homolytical cleavage of the O–NO₂ bonds.     

Experimental and theoretical studying the thermolysis of nitrate

This paper presents calorimeter measurement for the thermal decomposition of n-propyl nitrate (NPN), isopropyl nitrate (IPN) and 2-ethylhexyl nitrate (EHN). Similar experimental results of triethylene glycol dinitrate (tri-EGDN) and tetraethylene glycol dinitrate (tetra-EGDN) are included for comparison. The potential energy surfaces (PESs) along O-NO₂ bond stretch are investigated using the DFT (B3P86, B3PW91 and B3LYP), ab initio Hartree-Fock and PM3 methods. The good coincidence of experimental with theoretical results indicates that initial stage in the thermal decomposition of five nitrates is only unimolecular homolytical dissociation of the O-NO₂ bonds and the activation energies of thermolysis by DSC correspond to the energies of O-NO₂ bond scission of nitrates.     

Gas-phase structures of acetyl peroxynitrate and trifluoroacetyl peroxynitrate

The molecular structures and conformational properties of acetyl peroxynitrate (PAN, CH3C(O)OONO2) and trifluoroacetyl peroxynitrate (FPAN, CF3C(O)OONO2) were investigated in the gas phase by electron diffraction (GED), microwave spectroscopy (MW), and quantum chemical methods (HF/3-21G, HF/6-31G*, MP2/6-31G*, B3PW91/6-31G*, and B3PW91/6-311+G*). All experimental and theoretical methods show the syn conformer (C=O bond of acetyl group syn to O-O bond) to be strongly predominant relative to the anti conformer. The O-NO2 bonds are extremely long, 1.492(7) A in PAN and 1.526(10) A in FPAN, which correlates with their low bond energy and the easy formation of CX3C(O)OO* and *NO2 radicals in the atmosphere. The O-O bonds (1.418(12) A in PAN and 1.408(8) A in FPAN) are shorter than that in hydrogen peroxide (1.464 A). In both compounds the C-O-O-N dihedral angle is close to 85 degrees.     

Kinetics of the multichannel reaction of methanethiyl radical (CH3S*) with 3O2

The CH3S* + O2 reaction system is considered an important process in atmospheric chemistry and in combustion as a pathway for the exothermic conversion of methane-thiyl radical, CH3S*. Several density functional and ab initio computational methods are used in this study to determine thermochemical parameters, reaction paths, and kinetic barriers in the CH3S* + O2 reaction system. The data are also used to evaluate feasibility of the DFT methods for higher molecular weight oxy-sulfur hydrocarbons, where sulfur presents added complexity from its many valence states. The methods include: B3LYP/6-311++G(d,p), B3LYP/6-311++G(3df,2p), CCSD(T)/6-311G(d,p)//MP2/6-31G(d,p), B3P86/6-311G(2d,2p)//B3P86/6-31G(d), B3PW91/6-311++G(3df,2p), G3MP2, and CBS-QB3. The well depth for the CH3S* + 3O2 reaction to the syn-CH3SOO* adduct is found to be 9.7 kcal/mol. Low barrier exit channels from the syn-CH3SOO* adduct include: CH2S + HO2, (TS6, E(a) is 12.5 kcal/mol), CH3 + SO2 via CH3SO2 (TS2', E(a) is 17.8) and CH3SO + O (TS17, E(a) is 24.7) where the activation energy is relative to the syn-CH3SOO* stabilized adduct. The transition state (TS5) for formation of the CH3SOO adduct from CH3S* + O2 and the reverse dissociation of CH3SOO to CH3S* + O2 is relatively tight compared to typical association and simple bond dissociation reactions; this is a result of the very weak interaction. Reverse reaction is the dominant dissociation path due to enthalpy and entropy considerations. The rate constants from the chemical activation reaction and from the stabilized adduct to these products are estimated as functions of temperature and pressure. Our forward rate constant and CH3S loss profile are in agreement with the experiments under similar conditions. Of the methods above, the G3MP2 and CBS-QB3 composite methods are recommended for thermochemical determinations on these carbon-sulfur-oxygen systems, when they are feasible.     

((Fluoroformyl)imido)sulfuryl difluoride, FC(O)N=S(O)F2: structural, conformational, and configurational properties in the gaseous and condensed phases

Structural, conformational, and configurational properties of the gaseous molecule ((fluoroformyl)imido)sulfuryl difluoride, FC(O)N=S(O)F(2), have been studied by vibrational spectroscopy (IR (gas) and Raman (liquid)) and quantum chemical calculations (HF, MP2, and B3LYP with 6-31+G* and 6-311+G* basis sets); in addition, the solid-state structure has been determined by X-ray crystallography. FC(O)N=S(O)F(2) exists in the gas phase as a mixture of a favored antiperiplanar-synperiplanar form (the S=O double bond antiperiplanar with respect to the C-N single bond, and the C=O group synperiplanar with respect to the S=N double bond) in equilibrium with less abundant antiperiplanar-antiperiplanar, synclinal-synperiplanar, and synclinal-antiperiplanar structures. The crystalline solid at 163 K (monoclinic, P2(1)/c, a = 5.1323(7) A, b = 15.942(2) A, c = 16.798(2) A, beta = 95.974(3) degrees , Z = 12) consists of three similar antiperiplanar-synperiplanar forms.     

Theoretical study on the hydrogen bonding of five-membered heteroaromatics with water

The hydrogen-bonding ability of five-membered heteroaromatic molecules containing one chalcogen and two heteroatoms with nitrogen in addition to chalcogen, respectively, have been analyzed using density functional and molecular orbital methods through adduct formation with water. The stabilization energies for all the adducts are established at B3LYP/6-31+G* and MP2/6-31+G* levels after correcting for the basis set superposition error by using the counterpoise method and also corrected for zero-point vibrational energies. A natural bond orbital analysis at B3LYP/6-31+G* level and natural energy decomposition analysis at HF/6-31+G* using MP2/6-31+G* geometries have been carried out to understand the nature of hydrogen-bonding interaction in monohydrated heterocyclic adducts. Nucleus-independent chemical shift have been evaluated to understand the correlation between hydrogen bond formation and aromaticity.     

Vibrational analysis of isopropyl nitrate and isobutyl nitrate

Raman and infrared spectra of isopropyl nitrate and isobutyl nitrate are reported. These spectra are used in combination with computational studies employing density functional theory at the B3-LYP/6-31G* level to assign the vibrational transitions to their corresponding normal coordinates. Similar to other alkyl nitrates, the frequency of the NO2 symmetric stretch remains relatively unchanged while the asymmetric stretch shifts to lower frequency with increasing alpha-carbon substitution. The mode assignments involving the photochemically relevant -ONO2 chromophore agree well with those from previous infrared work. Raman depolarization ratios are also presented, and provide evidence that the condensed phase, ground-state molecular structure of isobutyl nitrate is of Cs symmetry. In contrast, the minimum energy structure of isopropyl nitrate is predicted to contain a pronounced twist around the C-O bond relative to the Cs-symmetry structure that lies 2.6 kcal/mol higher in energy. Infrared intensities of isopropyl nitrate are consistent with the twisted geometry, demonstrating that this conformer is favored in solution.