Theoretical Study on the Intermolecular Interactions between 1,1‐Diamino‐2,2‐Dinitroethylene and H2O

The density functional method was applied to the study of 1,1‐diamino‐2,2‐dinitroethylene (Fox‐7)/H₂O dimer. All the possible dimers ( 1, 2 and 3 ), as well as the monomers, were fully optimized with the DFT method at the B3LYP/6‐311++G^** level. The basis set superposition errors (BSSE) are 4.62, 4.07 and 3.45 kJ/mol, and the zero point energy (ZPE) corrections for the interaction energies are 7.94, 5.66 and 6.40 kJ/mol for 1, 2 and 3 , respectively. Dimer 1 is the most stable, judged by binding energy. After BSSE and ZPE corrections, the greatest corrected intermolecular interaction energy of dimer 1 was predicted to be ?29.36 kJ/mol. The charge redistribution mainly occurs on the adjacent N–H··· O atoms and N–O··· H atoms between submolecules. The oxygen in the nitro group acts as a moderate hydrogen acceptor as compared to water oxygen. Based on the statistical thermodynamic method, the standard thermodynamic functions, heat capacities (C⁰_P), entropies (S⁰_T) and thermal corrections to enthalpy (H⁰_T), and the changes of thermodynamic properties on going from monomer to dimer over the temperature range 200.00‐700.00 K were predicted. It is energetically or thermodynamically favorable for Fox‐7 to bind with H₂O and to form dimer 1 at room temperature.     

有机π共轭配体溶剂化效应与分子间相互作用的理论研究

采用可极化的连续介质模型(PCM), 运用密度泛函理论(DFT), 在B3LYP/6-31+G^**水平下研究了溶剂极性对有机π共轭配体N,N'-Bis-(3-pyridyl)ethylene-bis-urea(BPEBU)中syn-anti构象的分子几何和电子结构的影响, 并借助分子动力学模拟的方法, 采用明确溶剂模型研究了溶质-溶剂分子间的相互作用. 密度泛函理论计算结果表明, 随着溶剂极性的增强, BPEBU中尿素基上的CO键和N-H键以及吡啶环上的C-N键被明显极化, 使羰基氧原子和吡啶氮原子的电负性明显增强, 尿素基的N-H键上氢原子的正电荷也显著增加. 分子动力学模拟统计的结果表明, 在极性较强的乙醇溶液中, 有明确的O…H-O, N…H-O和N-H…O等3种氢键作用存在, 而在丙酮溶液中, 只有N…H-O一种氢键作用存在, 而且与乙醇溶液中的N…H-O作用相比要弱些. 另外, 采用密度泛函理论方法结合连续/明确的混合溶剂模型, 优化得到了溶质-溶剂三聚体的超分子簇结构, 与分子动力学模拟的第一溶剂层中的超分子结构相比, 两者定性一致.     

4-羟甲基吡啶－水红移氢键的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-311++G**基函数对4-羟甲基吡啶与水形成的1:1和1:2（摩尔比）氢键复合物进行了理论计算研究,分别得到稳定的4-羟甲基吡啶－H₂O和4-羟甲基吡啶－(H₂O)₂氢键复合物3个和8个。经基组重叠误差和零点振动能校正后,最稳定的1:1和1:2氢键复合物的相互作用能分别为-20.536和-44.246 kJ/mol。振动分析显示O-H···N(O)氢键的形成使复合物中O-H键对称伸缩振动频率红移(减小)。自然键轨道分析表明,4-羟甲基吡啶与水形成最稳定的1:1和1:2氢键复合物时,分子间电荷转移分别为0.02642 e 和0.03813 e 。含时密度泛函理论TD-B3LYP/ 6-311++G**计算显示,相对于4-羟甲基吡啶单体分子,氢键H-OH···N和H-OH···OH的形成分别使最大吸收光谱波长兰移8~16纳米和红移4~11纳米。     

Density Functional Study on the Reactivity of Carbenes Toward 1,2‐H Shifts

The 1,2‐H shift reactions of simple carbenes (CH₂Y‐C‐X) have been studied using density functional theory (DFT). The influence of the substituent X and Y groups on the activation energy (Ea) of 1,2‐H shifts were examined. The ‘by stander’ Y substituents lower E_a in the order of Me, F > Cl, Br > H. Our analysis shows that the X effect is more significant than the by stander Y effect. X substitutions increase E_a of carbenes in the order of F > Cl > Br > Me > H. The influence of X on E_a is governed by the singlet‐triplet energy separation (ΔE_S‐T) of the carbene, i.e., E_a of a carbene is larger as its ΔE_S‐T in creases due to an X substitution. The X effect was also found to be related to the magnitude of the exothermicity: Ea of reaction is smaller when the reaction is more exothermic. Origin of the Y effect is attributed to the inter play between two factors: ‘lateness’ of transition state on the potential energy surface, and the exothermicity of the reaction.     

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温度范围内单体形成二聚体的热力学性质变化. 发现在该温度范围所有二聚过程均能自发进行.     

NH3在硼纳米管表面吸附的密度泛函理论研究

利用密度泛函理论研究了NH₃在完整和含有缺陷的硼纳米管上的吸附行为以及相关电子性质. 计算结果表明, 对于α硼纳米管, 在不同的直径和手性条件下, NH₃均倾向于吸附在配位数为6的顶位上. 电子结构计算结果表明, NH₃能够吸附在纳米管表面主要是由于N和B原子产生了较强的相互作用. 表明硼纳米管是一种潜在的NH₃气气敏材料.     

HCN(HNC)与NH3, H2O和HF分子间相互作用的理论研究

在MP2/aug-cc-pVTZ水平上, 对HCN(HNC)与NH3, H2O和HF分子间可能存在的氢键型复合物进行了全自由度能量梯度优化, 通过在相同水平上的频率验证分析发现了稳定的分子间相互作用形式是HCN(HNC)作为质子供体或作为质子受体形成的复合物. 基组重叠误差对总相互作用能的影响均小于3.34 kJ/mol. 通过自然键轨道(NBO)分析, 研究了单体和复合物中的原子电荷和电荷转移对分子间相互作用的影响. 对称性匹配微扰理论(SAPT, Symmetry Adapted Perturbation Theory)能量分解结果表明, 在分子间相互作用中, 静电作用与诱导作用占主导地位, 而诱导作用与复合物的电荷转移之间具有良好的正相关性.     

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只在低温下才能自发形成Ⅳ,Ⅴ和Ⅵ超分子．     

中性团簇(CeO2)m(m=1~3)活化甲烷C—H的密度泛函理论计算

采用密度泛函理论(DFT)计算了CH₄在电中性(CeO₂)_m(m=1~3)团簇上的活化情况, 并对其机理进行了探讨. 计算结果表明, 甲烷C—H键在团簇上的活化为亲核加成模式, 电子由团簇流向甲烷C—H反键轨道, 使甲烷C—H键削弱而得以活化, 反应的过渡态为四中心结构. 团簇的桥氧位活化甲烷C—H键的活性大于端氧位, 而三重桥氧位的活性高于二重桥氧位. 团簇中作用位点Ce和O原子的电荷布居与其活化甲烷C—H的能力密切相关. 溶剂的存在不仅降低了甲烷C—H活化自由能垒, 而且使与甲烷作用的团簇各位点的活性差异缩小.     

Synthesis,Crystal Structure and Density Functional Studies on 1N‐Phenyl‐3‐(2,4‐Dichlorophenyl)‐5‐(4‐Chlorophenyl)‐2‐Pyrazoline

1N‐Phenyl‐3‐(2,4‐dichlorophenyl)‐5‐(4‐chlorophenyl)‐2‐pyrazoline has been synthesized and characterized by elemental analysis, IR, UV‐Vis and X‐ray single crystal diffraction. Density functional calculations have been carried out for the title compound by using the B3LYP method with a 6‐311G** basis set. The calculated results show that the predicted geometry can reproduce well the structural parameters. The electronic absorption spectra calculated in the gas phase are better than those calculated in EtOH solvent to model the experimental electronic spectra. Natural Bond Orbital (NBO) analyses suggest that the above electronic transitions are mainly assigned to π → π* transitions. On the basis of vibrational analyses, the thermodynamic properties of the compound at different temperatures have been calculated, revealing the correlations between C⁰_{p, m}, S⁰_m, H⁰_m and temperature.     

Density Functional Theory Studies on the Intermolecular Interactions of Five Aza-calix [6] arene Host with HMX

Five fully optimized structures of complexes between aza-calix[6]arene host monomers(Ma～Me) and complexes(a～e) have been obtained at the B3LYP/6-31G(d) level.Natural bond orbital(NBO) analysis was performed to reveal the origin of the interaction.The intermolecular interaction energy was evaluated with basis set superposition error correction(BSSE) and zero point energy correction(ZPEC).The B3LYP/6-31G(d) calculations on the five complexes have shown that the greatest interaction(-13.98 kJ/mol) is found in the complex between HMX and hexa-aza-calix[3]-p-tri-arene[3]-2-amido-1,3,5-tri-azine.The results have indicated that intermolecular interaction energies of aza-calix[6]arenes with substituted group are stronger than those without substituted group,and those with amido are greater than with nitryl.Thus,hexa-azacalix[3]-p-tri-arene[3]-2-amido-1,3,5-tri-azine is rather equal to eliminate HMX from explosive waste water.     

(CH2)2N和(CH3)2NH+的密度泛函理论计算

用密度泛函理论UB3LYP／6-31G（d,p）方法研究了二甲胺自由基（CH3）2N及其质子化离子（CH3）2NH 的构象和超精细结构．在由构象研究得到的两种自由基的最稳定结构上,用密度泛函的UB3LYP和UB3PW91方法及从头计算的UHF、UMP2（FULL）方法计算了α-质子、β-质子和N核上的超精细偶合常数A（Hα）、A（Hβ）和A（N）结果表明：两种自由基中甲基内旋转的位垒均很低,分别为0．46kJ·mol-1（（CH3）2NH ）和2．6kJ·mol-1（（CH3）2N）．UB3LYP／6-31G（d,p）和UB3PW91／6-31G（d,p）计算的A（Hα）、A（Hβ）和A（N）与ESR实验结果符合得很好,UMP2／6-31G（d,p）方法的计算值与实验值符合得也较好．     

弱氢键作用对硫脲水溶液拉曼光谱影响的密度泛函理论计算

将拉曼光谱和密度泛函理论(DFT)计算结合能够把拉曼谱峰和分子结构及分子间氢键作用之间的变化联系起来, 反映分子周围的结构信息. 本文通过理论计算水分子和硫脲形成复合物的拉曼光谱, 并结合实验报道的拉曼光谱, 探究硫脲基频对分子局域结构的依赖性. 通过轨道分析, 发现水分子的氢键作用可以引起硫脲的前线轨道对易, 这将影响到硫脲的拉曼光谱性质. 最后计算也表明在中性水溶液中由于大的正Gibbs 自由能变, 硫脲不易发生异构转变.     

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温度范围从单体形成二聚体的热力学性质变化,发现二聚主要由强氢键所贡献,二聚过程在较低温度或常温下能自发进行。     

弱氢键作用对硫脲水溶液拉曼光谱影响的密度泛函理论计算

将拉曼光谱和密度泛函理论(DFT)计算结合能够把拉曼谱峰和分子结构及分子间氢键作用之间的变化联系起来,反映分子周围的结构信息.本文通过理论计算水分子和硫脲形成复合物的拉曼光谱,并结合实验报道的拉曼光谱,探究硫脲基频对分子局域结构的依赖性.通过轨道分析,发现水分子的氢键作用可以引起硫脲的前线轨道对易,这将影响到硫脲的拉曼光谱性质.最后计算也表明在中性水溶液中由于大的正Gibbs自由能变,硫脲不易发生异构转变.     

Study of the [Cr(H2O)5NO]2+ Complex via Density Functional Theory

On the basis of density functional theory, the spin ground state of chromium‐nitrosyl complex [Cr(H₂O)₅NO]²⁺ (S = 1/2) is studied via B3LYP hybrid method. Its vibrational frequencies, atomic charges, and spin densities are analyzed. The excitation energies are evaluated using the CIS method. Our calculated N‐O stretching frequency and excitation energies are in good agreement with the IR and UV‐vis data. The related Cr^I(H₂O)₆⁺, Cr^II(H₂O)₆²⁺, and Cr^III(H₂O)₆³⁺ complexes are employed as the reference compounds to determine the characteristics of the central Cr. Results indicate that the effective Cr oxidation state is close to Cr(I).     

DFT Study of the Group Interactions in 1,3,5‐Triamino‐2,4,6‐Trinitrobenzene and 1,3,5‐Triamino‐2,4,6‐Tridifluoroaminobenzene

Density functional calculations on isodesmic disproportionation reactions of 1,3,5‐triamino‐2,4,6‐trinitrobenzene (TATB) and 1,3,5‐triamino‐2,4,6‐tridifluoroaminobenzene (TATDB) indicate that the interaction between nitro groups on meta carbons of TATB, which brings about unstability to the molecule, is surprisingly larger than that between difluroamino groups in TATDB. The electron‐withdrawing and electron‐donating groups generate large positive and very small negative values of E_{disproportion}, respectively. When both electron‐withdrawing and electron‐donating groups are attached to the benzene skeleton at the same time, large negative disproportionation energy is produced, which stabilizes the derivatives. The values of E_{disproportion} for TATB and TATDB are predicted to be ‐48.03 kJ/mol and ‐63.54 kJ/mol, respectively, indicating that the total interaction among groups with stabilization effects in TATDB is larger than that in TATB. The large difference of the E_{disproportion} values between TATB and TATDB is derived from the large difference between the interactions of the meta‐nitro group and those of meta‐difluoroamino groups. The energy barriers for the C‐N internal rotation of NO₂ group and NF₂ groups are 74.7 kJ/mol and 185.8 kJ/mol for TATB and TATDB, respectively. The large energy barrier for the rotation of the NF₂ group is caused by its stabilization interaction with neighbor amino groups, instead of steric effects. When the number of pairs of amino‐nitro or amino‐difluoroamino groups increases, there are more negative charges on the NO₂/NF₂ groups and on the O/F atoms.     

Study on Chromium‐Containing Aluminophosphate CrAPO‐5 via Density Functional Theory

Various structures of CrAPO‐5 clusters are studied via density functional B3LYP exchange‐correlation method. The optimized structures are compared with data from X‐ray absorption. Their total energies and atomic net charges are also analyzed. Results indicate that the substitution of the aluminum site of an AFI framework by chromium is in general not feasible. The chromium ion is more likely docked in between two neighboring 12‐membered rings of the framework of AFI. To further verify our claim, the excitation energies of the representative chromium structures of CrAPO‐5 clusters are calculated via the TDDFT method. The results for excitation energies further support that Cr³⁺ is not incorporated into the framework.     

A density functional theory study on diazotization and nitration of 3,5‐diamino‐1,2,4‐triazole

The reaction mechanism, thermodynamic and kinetic properties for diazotization and nitration of 3,5‐diamino‐1,2,4‐triazole were studied by a density functional theory. The geometries of the reactants, transition states, and intermediates were optimized at the B3LYP/6‐31G (d, p) level. Vibrational analysis was carried out to confirm the transition state structures, and the intrinsic reaction coordinate (IRC) method was used to explore the minimum energy path. The single‐point energies of all stagnation points were further calculated at the B3LYP (MP2)/6‐311+G (2d, p) level. The statistical thermodynamic method and Eyring transition state theory with Wigner correction were used to study the thermodynamic and kinetic characters of all reactions within 0–25°C. Two reaction channels are computed, including the diazotization and nitration of 3‐NH₂ or 5‐NH₂, and there are six steps in each channel. The reaction rate in each step is increased with temperature. The last step in each channel is the slowest step. The first, second, and fifth steps are exothermic reactions, and are favored at lower temperature in the thermodynamics. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011