NTO二聚体分子间相互作用的理论研究

在DFT-B3LYP/6-311＋＋G**水平上求得NTO二聚体势能面上六种优化构型和电子结构. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得分子间最大相互作用能为－53.66 kJ/mol. 二子体系间的电荷转移很少. 由自然键轨道分析揭示了相互作用的本质. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围从单体形成二聚体的热力学性质变化. 发现二聚主要由强氢键所贡献, 但结合能大小并不为氢键所完全决定. 二聚过程在较低温度或常温下能自发进行.     

3,6-二氨基-1,2,4,5-四嗪二聚体分子间相互作用的理论研究

在DFT-B3LYP/6-31G(d)水平下,求得3,6-二氨基-1,2,4,5-四嗪二聚体势能面上3种优化几何构型和电子结构。经基组叠加误差(BSSE)和零点能(ZPE)校正,求得分子间最大相互作用能为-38.88kJ/mol。电荷分布与转移分析表明,二子体系间的电荷转移很少,但接触点上氮原子和氢原子电荷变化比较大。由自然键轨道(NBO)分析揭示了分子间相互作用的本质。对优化构型进行振动分析,并基于统计热力学求得200.0~800.0 K温度范围从单体形成二聚体的热力学性质变化,发现二聚主要由强氢键所贡献,二聚过程在较低温度或常温下能自发进行。     

3-硝基-1,2,4-三唑-5-酮与NH3及H2O分子间相互作用的理论研究

在DFT—B3LYP／6—311 G^**水平上,求得3-硝基-1,2,4-三唑-5-酮(NTO)／NH3和NTO／H2O两种超分子体系势能面上5种全优化构型．经基组叠加误差(BSSE)和零点能(ZPE)校正,求得NTO与NH3和H2O的分子间最大相互作用能依次为-37．58和-30．14kJ／mol,表明NTO与NH3的分子问相互作用强于与H2O的作用．超分子体系中电子均由NH3或H2O向NTO转移,相互作用能主要由强氢键所贡献,由自然键轨道分析揭示了相瓦作用的本质．对优化构型进行振动分析,并基于统计热力学求得200．0～800．0K温度范围从单体形成超分子的热力学性质变化．发现由NTO和NH3形成超分子Ⅱ和Ⅲ在常温下可自发进行;而NTO和H2O只在低温下才能自发形成Ⅳ,Ⅴ和Ⅵ超分子．     

3-硝基-1,2,4-三唑-5-酮与NH3及H2O分子间相互作用的理论研究

在DFT-B3LYP/6-311＋＋G**水平上, 求得3-硝基-1,2,4-三唑-5-酮(NTO)/NH₃和NTO/H₂O两种超分子体系势能面上5种全优化构型. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得NTO与NH₃和H₂O的分子间最大相互作用能依次为－37.58和－30.14 kJ/mol, 表明NTO与NH₃的分子间相互作用强于与H₂O的作用. 超分子体系中电子均由NH₃或H₂O向NTO转移, 相互作用能主要由强氢键所贡献, 由自然键轨道分析揭示了相互作用的本质. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围从单体形成超分子的热力学性质变化. 发现由NTO和NH₃形成超分子II和III在常温下可自发进行; 而NTO和H₂O只在低温下才能自发形成IV, V和VI超分子.     

S-亚硝基-N-乙酰基-D,L-青霉胺二肽分子中S—NO键断裂能的测定

利用滴定量热技术并结合适当的热力学循环测定了乙腈溶液中7个S-亚硝基-N-乙酰基-D,L-青霉胺二肽化合物中S—NO键的异裂能和均裂能, 其能量范围分别为234.5—246.2 kJ/mol和101.6—122.1 kJ/mol. 结果表明, 所研究的亚硝基硫醇化合物更容易通过S—NO键的均裂释放NO自由基(NO·). 通过热力学循环对7个亚硝基硫醇化合物自由基负离子中S—NO键的异裂能和均裂能进行估算, 能量范围分别为19.2—35.5 kJ/mol和-4.2—22.6 kJ/mol, 表明这些自由基负离子在室温下不稳定, 容易通过S—NO键的异裂释放出NO-.     

聚(N-甲基-N-对乙烯基苄基乙酰胺)对水中苯酚的吸附

由大孔交联氯甲基化聚苯乙烯进行甲胺化和乙酰化反应,合成了聚(N-甲基-N-对乙烯基苄基乙酰胺)树脂,树脂的氯含量由原来的3.95mmol/g降至0.63mmol/g.测定了树脂对水中苯酚的吸附等温线,结果表明,吸附平衡数据符合Freundlich等温方程.利用热力学函数关系计算了吸附焓、吸附自由能和吸附熵,结果表明,吸附量在10～40mg/g时,吸附焓为-71.20～-24.65kJ/mol,吸附自由能为-4.64～-6.07kJ/mol,吸附熵为-62.35～-218.56J/(mol·K).讨论了吸附过程中树脂对水中苯酚吸附的氢键与疏水作用的吸附机理.     

1H-ANTA二聚体分子间相互作用的密度泛函理论研究

在DFT-B3LYP/6-311＋＋G**水平上求得1H-3-硝基-5-氨基-1,2,4-三唑(1H-ANTA)二聚体势能面上5种优化构型和电子结构. 经基组叠加误差(BSSE)和零点能(ZPE)校正, 求得分子间最大结合能为70.63 kJ/mol. 二聚体的形成使电荷向三唑环转移. 由氢键强弱推断二聚体稳定性的顺序与结合能顺序相一致, 氢键是二聚体的主要作用形式. 对优化构型进行振动分析, 并基于统计热力学求得200.0～800.0 K温度范围内单体形成二聚体的热力学性质变化. 发现在该温度范围所有二聚过程均能自发进行.     

2-巯基苯并噻唑直接键合硅胶的制备及其对Cu（Ⅱ）的吸附性能研究

采用一种简单有效的方法合成了2-巯基苯并噻唑(MBT)直接键合硅胶.将键合硅胶用于水中Cu(Ⅱ)的吸附,考查了溶液pH值的影响.结果表明当pH为6时,键合硅胶的吸附量最大.通过分批实验,研究了278K、293K和308K 3个温度条件下的吸附平衡.应用Langmuir和Freundlich两个双参数等温线模型进行拟合,MATLAB软件的误差分析表明,Freundlich模型与实验数据更加吻合.计算得到吸附过程的热力学参数△G°、△H°和△S°的值分别为-21.90kJ/mol(293K)、14.57kJ/mol和124.61J/(mol·K).结果表明,键合硅胶对Cu(Ⅱ)的吸附是一个自发进行的吸热过程.     

3,7-二硝亚胺基-2,4,6,8-四硝基-2,4,6,8-四氮杂双环[3.3.0]辛烷结构与性能的量子化学研究

设计了一种新型高能量密度化合物--3,7-二硝亚胺基-2,4,6,8-四硝基-2,4,6,8-四氮杂双环[3.3.0]辛烷, 应用密度泛函理论(DFT)的B3LYP 方法在6-31G(d,p)基组水平上对该化合物进行了结构全优化, 并计算得到其红外(IR)光谱; 通过键级分析获得热解引发键的位置为N7-N22, 同时求得校正后的键离解能为91.47 kJ/mol. 采用Monte-Carlo方法预测该化合物的理论密度为2.16 g/cm³; 基于理论密度并结合等键反应及Kamlet-Jacobs公式预测了生成焓、爆速、爆压和爆热值分别为1219.94 kJ/mol, 10.43 km/s, 53.44 GPa和7407.84 J/g. 以上性能参数显示, 该目标化合物达到了高能量密度化合物的基本要求, 是一种潜在的含能材料. 同时给出了该化合物的逆合成路线.     

烷基芳基磺酸盐分子量及其分布对微乳液结构参数及热力学函数的影响

用稀释法求得了由自制的7种烷基芳基磺酸盐复配体系(AAS)/正丁醇/正癸烷/水组成的W/O型微乳液的结构参数,还求得了醇从油相转移到界面时的标准自由能,并计算出标准焓变和标准熵变。探讨了烷基芳基磺酸盐平均分子量及其分布对结构参数及热力学函数的影响。结果表明,分散相有效半径Re,内核水半径Rw,二者之差di和表面活性剂分子在每个液滴中的平均聚集数n值均呈现为正态分布<递减分布<均匀分布<递增分布<反正态分布;分散相颗粒总数Nd和分散相界面外层总面积Ad值均呈现为正态分布>递减分布>均匀分布>递增分布>反正态分布;-ΔG0o→i值呈现为正态分布(5.36 kJ/mol)<递减分布(5.49 kJ/mol)<均匀分布(5.64 kJ/mol)<递增分布(5.78 kJ/mol)<反正态分布(6.01 kJ/mol);ΔS0o→i值呈现为正态分布(26.88 J/(mol.K))<递减分布(27.12 J/(mol.K))<均匀分布(27.60 J/(mol.K))<递增分布(28.06 J/(mol.K))<反正态分布(29.23 J/(mol.K));Rw、Re、n、di、-ΔG0o→i、ΔH0o→i和ΔS0o→i值均随磺酸盐平均分子量的增大而增大;Nd、Ad值均随磺酸盐平均分子量的增大而减小;且在实验范围内,结构参数、-ΔG0o→i、ΔS0o→i与磺酸盐平均分子量均呈线性关系;后两者分别为y=0.0586x-17.916,y=0.2203x-61.275。     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Nitroamine Dimers

Ab initio self-consistent field(SCF) and Mφller-Plesset correlation correction methods employing 6-31G^** basis set have been applied to the optimizations of nitroamine dimers.The binding energies have been corrected for the basis set superposition error (BSSE) and the zero-point energy.Theree optimized dimers have been obtained.The BSSE corrected binding energy of the most stable dimer is predicted to be -31.85kJ/mol at the MP4/6-31G^**//MP2/6-31G^** level.The energy barriers of the Walden conversion for -NH2 group are 19.7kJ/mol and 18.3kJ/mol for monomer and the most stable dimer,respectively.The molecular interaction makes the internal rotation around N1-N2 even more difficult.The thermodynamic properties of nitroamine and its dimers at different temperatures have been calculated on the basis of vibrational analyses.The change of the Gibbs free energy for the aggregation from monomer to the most stable dimer at standard pressure and 298.2 K is predicted to be 14.05kJ/mol.     

A density-functional theory investigation of 3-nitro-1,2,4-triazole-5-one dimers and crystal

Density-functional method with different basis sets was applied to the study of the highly efficient and low sensitive explosive 3-nitro-1,2,4-triazole-5-one (NTO) in both gaseous dimer and its bulk state. The binding energies have been corrected for the basis set superposition errors. Six stable dimers (II-VII) were located. The corrected binding energy of the most stable dimer VII is predicted to be -53.66 kJ/mol at the B3LYP/6-311++G(**) level. It was found that the structures of the more stable dimers (V-VII) are through the hydrogen bonding interaction between the carbonyl oxygen and the azole hydrogen of 3-nitro-1,2,4-triazole-5-one. The changes of Gibbs free energies (DeltaG) in the processes from the monomer to the dimers at 298.15 K are 8.51, 0.90, 0.35, -8.74, -10.67, and -11.06 kJ/mol for dimers from II to VII, respectively. Dimers V-VII, possessing cyclic structures, can be spontaneously produced from the isolated monomer at room temperature. The lattice energy is -156.14 kJ/mol, and this value becomes to -150.43 kJ/mol when a 50% correction of the basis set superposition error was adopted. The frontier bands are quite flat. Judged from the value of band gap of 4.0 eV, it may be predicted that 3-nitro-1,2,4-triazole-5-one is an insulator. Most atoms in NTO, with the exception of C(5) atom and the nitro atoms, make up the upper valence bands. In contrast, the lower conduction bands mainly consist of the nitro N and O atoms. The population of the C-NO(2) bond is much less than those of the other bonds and the detonation may be initiated by the breakdown of this bond.     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Dimethylnitroamine Clusters

Ab initio SCF and Mφller-Plesset correlation correction methods in combination with counterpose procedure for BSSE correction have been applied to the theroetical studying of dimethylnitroamine and its dimers and trimers.Three optimized stable dimers and two trimers have been obtained.The corrected binding energies of the most stable dimer and trimer were predicted to be -24.68kJ/mol and -47.27kJ/mol,respectively at the MP2/6-31G^*//HF/6-31G^* level.The proportion of correlated interation energies to their total interaction energies for all clusters was at least 29.3 percent,and the BSSE of ΔE(MP2) was at least 10.0kJ/mol.Dispersion and/or electrostatic force were dominant in all clusters.There exist cooperative effects in both the chain and the cyclic trimers.The vibrational frequencies associated with N-O stretches or wags exhibit slight red shifts,but the modes associated with the motion of hydrogen atoms of the methyl group show somewhat blue shifts with respect to those of monomer.Thermodynamic properties of dimethylnitroamine and its clusters at different temperatures have been calculated on the basis of vibrational analyses.The changes of the Gibbs free energies for the aggregation from monomer to the most stable dimer and trimer were predicted to be 14.37kJ/mol and 30.40kJ/mol,respectively,at 1 atm and 298.15K.     

Ab initio Studies on Intermolecular Interaction of Formamide and Hydroxyacetonitrile Dimers

The structures, the binding energies and the thermodynamic properties of formamide and hydroxyacetonitrile(HAN) dimers have been studied by means of the self-consistent ab initio Hartree-Fock and the second-order Mφller-Plesset correlation energy correction methods. The counterpoise procedure was used to check the basis set superposition error(BSSE) of the binding energies. There exist cyclic structures in a formamide dimer(Ⅰ), a HAN dimer(Ⅱ) and their heterodimer(Ⅲ). The corrected binding energies for dimers Ⅰ, Ⅱ and Ⅲ are respectively -45.53, -45.83 and -43.89 kJ/mol at the MP2/aug-cc-p VDZ//HF/aug-cc-p VDZ level. The change of the Gibbs free energies(ΔG) in the process of Ⅰ Ⅱ→2Ⅲ was predicted to be -2.74 kJ/mol at 298.15 K. Dimer Ⅲ can be spontaneously produced in the mixture of formamide and HAN, which is in agreement with the experimental fact that most cyanohydrins are capable of interacting with dipeptide cyclo-His-Phe(CHP).     

Clusters of hydrated methane sulfonic acid CH3SO3H.(H2O)n (n = 1-5): a theoretical study

Ab initio and density functional methods have been used to examine the structures and energetics of the hydrated clusters of methane sulfonic acid (MSA), CH3SO3H.(H2O)n (n = 1-5). For small clusters with one or two water molecules, the most stable clusters have strong cyclic hydrogen bonds between the proton of OH group in MSA and the water molecules. With three or more water molecules, the proton transfer from MSA to water becomes possible, forming ion-pair structures between CH3SO3- and H3O+ moieties. For MSA.(H2O)3, the energy difference between the most stable ion pair and neutral structures are less than 1 kJ/mol, thus coexistence of neutral and ion-pair isomers are expected. For larger clusters with four and five water molecules, the ion-pair isomers are more stable (>10 kJ/mol) than the neutral ones; thus, proton transfer takes place. The ion-pair clusters can have direct hydrogen bond between CH3SO3- and H3O+ or indirect one through water molecule. For MSA.(H2O)5, the energy difference between ion pairs with direct and indirect hydrogen bonds are less than 1 kJ/mol; namely, the charge separation and acid ionization is energetically possible. The calculated IR spectra of stable isomers of MSA.(H2O)n clusters clearly demonstrate the significant red shift of OH stretching of MSA and hydrogen-bonded OH stretching of water molecules as the size of cluster increases.     

尿素与鸟嘌呤相互作用的DFT研究

采用B3LYP／6—311＋G^＋方法对鸟嘌呤-尿素复合物氢键相互作朋体系进行了研究,并对该复合物的几何构型及结合能（BSSE）进行了计算．此外,采用从静电势导出原子净电荷的chelpg方法分析了体系中的电荷转移和利用分子中的原子理论（AIM）方法对相互作用的本质进行了分析．结果一共得到五个稳定的复合物构型,其中A5是最稳定的,结合能为-73．95kJ／mol．     

Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

Elimination of methane from protonated acetaldehyde: The Ab initio transition state

The transition state (TS) for loss of CH₄ from protonated acetaldehyde has been located at the second-order Moller-Plesset (MP2)/6-31G(d,p) level of theory. The activation energy is predicted to be 263.9 kJ/mol starting from the more stable form (methyl and hydrogen E) and 261.6 kJ/mol starting from the less stable form (methyl and hydrogen Z) that is required for reaction. The products (methane and the formyl ion) are predicted to lie 136.6 kJ/mol below the TS for their formation. MP2 methods underestimate the heats of formation of both the TS and the reaction products by about 40 kJ/mol when compared with experiment. Restricted Hartree-Fock (RHF) calculations give much more accurate relative energies. The MP2 TS leads directly to fragmentation and is described as a protonation of the methyl group by the acidic proton on oxygen. Under RHF theory the reaction is stepwise. An RHF TS similar to the MP2 TS leads to a nonclassical intermediate (which is stable at this level of theory) that has one of the C---H bonds protonated. This mechanism (protonation of an alkyl group) appears to be a general one for high energy 1,2 eliminations from organic cations. (J Am Soc Mass Spectrom 1994, 5, 1102-1106)     

Direct determination of the ionization energies of PtC, PtO, and PtO2 with VUV radiation

Photoionization efficiency curves were measured for gas-phase PtC, PtO, and PtO2 using tunable vacuum ultraviolet (VUV) radiation at the Advanced Light Source. The molecules were prepared by laser ablation of a platinum tube, followed by reaction with CH4 or N2O and supersonic expansion. These measurements provide the first directly measured ionization energy for PtC, IE(PtC) = 9.45 +/- 0.05 eV. The direct measurement also gives greatly improved ionization energies for the platinum oxides, IE(PtO) = 10.0 +/- 0.1 eV and IE(PtO2) = 11.35 +/- 0.05 eV. The ionization energy connects the dissociation energies of the neutral and cation, leading to greatly improved 0 K bond dissociation energies for the neutrals: D0(Pt-C) = 5.95 +/- 0.07 eV, D0(Pt-O) = 4.30 +/- 0.12 eV, and D0(OPt-O) = 4.41 +/- 0.13 eV, as well as enthalpies of formation for the gas-phase molecules DeltaH(0)(f,0)(PtC(g)) = 701 +/- 7 kJ/mol, DeltaH(0)(f,0)(PtO(g)) = 396 +/- 12 kJ/mol, and DeltaH(0)(f,0)(PtO2(g)) = 218 +/- 11 kJ/mol. Much of the error in previous Knudsen cell measurements of platinum oxide bond dissociation energies is due to the use of thermodynamic second law extrapolations. Third law values calculated using statistical mechanical thermodynamic functions are in much better agreement with values obtained from ionization energies and ion energetics. These experiments demonstrate that laser ablation production with direct VUV ionization measurements is a versatile tool to measure ionization energies and bond dissociation energies for catalytically interesting species such as metal oxides and carbides.     

On the origin of the stabilization of the zwitterionic resting state of cysteine proteases: a theoretical study

Papain-like cysteine proteases are ubiquitous proteolytic enzymes. The protonated His199/deprotonated Cys29 ion pair (cathepsin B numbering) in the active site is essential for their proper functioning. The presence of this ion pair stands in contrast to the corresponding intrinsic residue p K a values, indicating a strong influence of the enzyme environment. In the present work we show by molecular dynamics simulations on quantum mechanical/molecular mechanical (QM/MM) potentials that the ion pair is stabilized by a complex hydrogen bond network which comprises several amino acids situated in the active site of the enzyme and 2-4 water molecules. QM/MM reaction path computations for the proton transfer from His199 to the thiolate of the Cys29 moiety indicate that the ion pair is about 32-36 kJ mol (-1) more stable than the neutral form if the whole hydrogen bonding network is active. Without any hydrogen bonding network the ion pair is predicted to be significantly less stable than the neutral form. QM/MM charge deletion analysis and QM model calculations are used to quantify the stabilizing effect of the active-site residues and the L1 helix in favor of the zwitterionic form. The active-site water molecules contribute about 30 kJ mol (-1) to the overall stabilization. Disruption of the hydrogen bonding network upon substrate binding is expected to enhance the nucleophilic reactivity of the thiolate.