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

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

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超分子.     

B炸药主要组分TNT和RDX分子间相互作用的理论研究

在B3LYP/6-31G(d)水平上研究了B炸药的主要成分--2,4,6-三硝基甲苯(TNT)与环三亚甲基三硝胺(RDX)分子间的相互作用, 得到了10种TNT＋RDX的全优化构型. 讨论了稳定构型在几何参数、稳定性、红外光谱和电荷分布上的差异. 借助自然键轨道(NBO)理论揭示了TNT与RDX分子间相互作用的本质, 主要由氢键所贡献. 分子间相互作用能在－3.930～－14.652 kJ•mol^－1之间, 经基组叠加误差(BSSE)校正, 相互作用能顺序为VI＞III＞V＞IV＞X＞I＞IX＞II＞VII＞VIII. 对全优化构型进行了热力学性质的分析, 探讨了由单体分子形成混合体系的热力学性质的变化, 结果发现, 形成分子间氢键是个放热过程. 运用Kamlet-Jacobs方程基于理论密度(ρ)估算了混合体系TNT＋RDX的爆轰性质爆速(D)和爆压(p), 与文献值进行比较表明理论计算方法和结果是可靠的.     

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

TATB二聚体分子间作用力及其气相几何构型研究

采用对称性匹配微扰理论(SAPT)定量地求得TATB分子间的静电、交换排斥、诱导和色散等分子间作用能项, 从理论上揭示了TATB分子间作用本质; 在此基础上, 阐明了密度泛函在研究TATB二聚体时的适合性问题. 结果表明: (1)在有分子间氢键的TATB二聚体中, 库仑力足以与交换排斥力相抗衡, 起主导作用. (2)含分子间氢键的气相TATB二聚体的合理几何构型为平面型结构, 此结构的产生与色散力无关, 因此不管泛函是否含有近程色散作用, 均应预测到这种强极性的平面型结构. (3)在无分子间氢键的TATB二聚体中, 库仑力难以与交换排斥力相抗衡, 色散作用起到了关键作用; (4)在这种情况下, 未含有近程色散作用的密度泛函不可能给出合理构型. 恰好相反, 含有近程色散作用的密度泛函PBE0却能正确地预测到具有“平行重叠”结构且呈微弱极性的TATB二聚体, 色散力是导致这种构型产生的根本原因. “平行重叠”TATB二聚体是典型的色散体系, 其色散力占绝对主导地位并极有可能起源于两个TATB分子上π电子的相互作用. (5)对于所有TATB二聚体, 色散力或很显著或起主导作用. 由于密度泛函或未含有近程色散, 或只能部分地把近程色散表达出来, 这样使得当前所有密度泛函不可能精确求得这些二聚体的作用能.     

硝酸甲酯分子间相互作用的DFT和ab initio比较

用密度泛函理论(DFT)和从头算(ab initio)方法,分别在B3LYP/6 31G和HF/6 31G水平上求得硝酸甲酯三种二聚体的全优化几何构型和电子结构,并用6 311G和6 311＋＋G基组进行总能量计算.对HF/6 31G计算结果进行MP4SDTQ电子相关校正.在各基组下均进行基组叠加误差(BSSE)和零点能(ZPE)校正求得结合能.对6 31G优化构型作振动分析并基于统计热力学求得200～600 K温度下单体和二聚体的热力学性质.详细比较两种方法的相应计算结果,发现DFT求得的分子间距离较短,分子内键长较长,所得结合能均小于相应ab initio计算值.     

半胱氨酸阴离子与咪唑阳离子间相互作用的理论研究

采用密度泛函理论在B3LYP/6-311＋G(d,p)水平上对1-乙基-3-甲基咪唑阳离子和半胱氨酸阴离子形成的气态阴阳离子对([Emim][Cys])进行理论研究. 通过几何结构优化以及频率分析得到势能面上7个稳定的离子对构型. 计算结果表明[Emim]＋和[Cys]－之间存在较强的氢键相互作用, 其稳定化能主要来源于[Cys]－中羰基O的孤对电子lp(O) 和[Emim]＋中C—H反键轨道 s*(C—H) 之间的相互作用, lp(O)®s*(C—H). [Emim][Cys]_S1是最稳定的离子对构型, 考虑BSSE的相互作用能为－387.66 kJ/mol. 从NPA和NBO分析以及AIM (Atoms in Molecules)计算等方面阐述了半胱氨酸阴离子与咪唑阳离子之间氢键相互作用的本质, 并初步探讨了阴阳离子对相互作用对氨基酸离子液体性质的影响.     

硝酸乙酯分子间相互作用的ab initio研究

在abinitio-HF/6-31G水平上求得硝酸乙酯二聚体势能面上的四种优化构型和电子结构。经MP2电子相关校正和基组叠加误差(BSSE)以及零点能(ZPE)校正，求得二聚体的最大结合能为11.46kJ.mol^-^1，还进行HF/6-311G和HF/6-311++G水平的总能量比较计算，发现6-31G基组对计算结合能比较适合，二子体系间的电荷转移很少，对优化构型进行振动分析，并基于统计热力学求得从单体形成二聚体的热力学性质变化。     

硝仿肼离子对相互作用的密度泛函理论研究

用密度泛函理论(DFT)方法,在B3LYP/6-31+G^**水平下,求得硝仿肼离子对体系势能面上2种全优化构型.经基组叠加误差(BSSE)和零点能校正,求得离子对最大相互作用能为-420.03kJ/mol,肼和硝仿离子化所需能量可由该值得到完全补偿.离子对间键的主要贡献为库仑作用,但键鞍点上的电子密度值表明共价作用也有显著的贡献.基于统计热力学求得相关体系的热力学性质,298.2K时由自由离子形成最稳定离子对的最大焓变和最大自由能变化分别为-419.72和-376.61kJ/mol     

Theoretical Studies on Intermolecular Interactions of 4-Amino-5-nitro-1,2,3-triazole Dimers

1 INTRODUCTION Triazole, a five-membered heterocyclic compound with three nitrogen atoms, is an important interme- diate product of medicine and chemical industry as well as insecticide [1]. Due to its small volume and high nitrogenous density, triazole holds more and more attraction for the material researchers, espe- cially the researchers of high-energy insensitive explosive. It is reported that its nitro and amino deri- vatives are a sort of important high-energy mate- rials[2]. Previ…     

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.     

Density functional theory study of the interaction between 3‐nitro‐1,2,4‐triazole‐5‐one and ammonia

Two fully optimized geometries of 3‐nitro‐1,2,4‐triazol‐5‐one (NTO)–NH₃ complexes have been obtained with the density function theory (DFT) method at the B3LYP/6‐311++G** level. The intermolecular interaction energy is calculated with zero point energy (ZPE) correction and basis set superposition error (BSSE) correction. The greatest corrected intermolecular interaction of the NTO–NH₃ complexes is ?37.58 kJ/mol. Electrons in complex systems transfer from NH₃ to NTO. The strong hydrogen bonds contribute to the interaction energies dominantly. Natural bond orbital (NBO) analysis is performed to reveal the origin of the interaction. Based on vibrational analysis, the changes of thermodynamic properties from the monomer to complexes with the temperature ranging from 200 K to 800 K have been obtained using the statistical thermodynamic method. It is found that two NTO–NH₃ complexes can be produced spontaneously from NTO and NH₃ at normal temperature. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Theoretical Study on the Intermolecular Interactions of Tetrazole Dimers

1 INTRODUCTION Tetrazole and its derivatives are widely applied in the fields of agriculture, biology, chemistry, phar- macology and photographic technology, and they play significant roles in the science and technology as well as national defence[1]. In the past, the res- earches were focused on the molecular geometries, electronic structures, IR, thermodynamic properties, tautomerization, pyrogenation and sensitivity of tetrazole compounds[1～5]. However, study of tetra- zole dimers ha…     

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.     

Magnitude and orientation dependence of intermolecular interaction of perfluoropropane dimer studied by high-level ab initio calculations: comparison with propane dimer

Intermolecular interaction energies of 12 orientations of C(3)F(8) dimers were calculated with electron correlation correction by the second-order M?ller-Plesset perturbation method. The antiparallel C(2h) dimer has the largest interaction energy (-1.45 kcal/mol). Electron correlation correction increases the attraction considerably. Electrostatic energy is not large. Dispersion is mainly responsible for the attraction. Orientation dependence of the interaction energy of the C(3)F(8) dimer is substantially smaller than that of the C(3)H(8) dimer. The calculated interaction energy of the C(3)F(8) dimer at the potential minimum is 78% of that of the C(3)H(8) dimer (-1.85 kcal/mol), whereas the interaction energies of the CF(4) and C(2)F(6) dimers are larger than those of the CH(4) and C(2)H(6) dimers. The intermolecular separation in the C(3)F(8) dimer at the potential minimum is substantially larger than that in the C(3)H(8) dimer. The larger intermolecular separation due to the steric repulsion between fluorine atoms is the cause of the smaller interaction energy of the C(3)F(8) dimer at the potential minimum. The calculated intermolecular interaction energy potentials of the C(3)F(8) dimers using an all atom model OPLS-AA (OPLS all atom model) force field and a united atom model force field were compared with the ab initio calculations. Although the two force fields well reproduces the experimental vapor and liquid properties of perfluoroalkenes, the comparison shows that the united atom model underestimates the potential depth and orientation dependence of the interaction energy. The potentials obtained by the OPLS-AA force field are close to those obtained by the ab initio calculations.     

High-level ab initio computations of structures and interaction energies of naphthalene dimers: origin of attraction and its directionality

The intermolecular interaction energies of naphthalene dimers have been calculated by using an aromatic intermolecular interaction model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order M?ller-Plesset perturbation interaction energy near saturation and the CCSD(T) correction term obtained using a medium-size basis set. The estimated interaction energies of the set of geometries explored in this work show that two structures emerge as being the lowest energy, and may effectively be considered as isoenergetic on the basis of the errors inherent in out extrapolation procedure. These structures are the slipped-parallel (Ci) structure (-5.73 kcal/mol) and the cross (D2d) structure (-5.28 kcal/mol). The T-shaped (C2v) and sandwich (D2h) dimers are substantially less stable (-4.34 and -3.78 kcal/mol, respectively). The dispersion interaction is found to be the major source of attraction in the naphthalene dimer. The electrostatic interaction is substantially smaller than the dispersion interaction. The large dispersion interaction is the cause of the large binding energies of the cross and slipped-parallel dimers.