An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

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

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

甲醛与乙醛,甲醚,硝基甲烷相互作用的从头算研究

用6-31G、全构型优化,研究了甲醛与乙醛、甲醚、硝基甲烷的相互作用。结果表明所有超分子稳定构型都包含2个C—H—O氢键的平面环状结构。H—O距离为0.228～0.264 nm,作用能为—19～—24 kJ/mol,与二聚水的作用能接近。稳定性主要取决于甲基上取代基Y吸电子能力以及环状结构中氢键张力。STO-3G不很适用于研究这类分子的相互作用。     

Pi-pi interaction in pyridine

pi-pi Interaction in pyridine dimer and trimer has been investigated in different geometries and orientations at the ab initio (HF, MP2) and DFT (B3LYP) levels of theory using various basis sets (6-31G, 6-31G, 6-311++G) and corrected for basis set superposition error (BSSE). While the HF and DFT calculations show the pyridine dimer and the trimer to be unstable with respect to the monomer, the MP2 calculations show them to be clearly stable, thus emphasizing the need to include electron correlation while determining stacking interaction in such systems. The calculated MP2/6-311++G binding energy (100% BSSE corrected) of the parallel-sandwich, antiparallel-sandwich, parallel-displaced, antiparallel-displaced, T-up and T-down geometries for pyridine dimer are 1.53, 3.05, 2.39, 3.97, 1.91, 1.47 kcal/mol, respectively. The results show the antiparallel-displaced geometry to be the most stable. The binding energies for the trimer in parallel-sandwich, antiparallel-sandwich, and antiparallel-displaced geometry are found to be 3.18, 6.14, and 8.04 kcal/mol, respectively.     

硝酸甲酯分子间相互作用的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计算值.     

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.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

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 study of new Pt(II) and Pd(II) complexes with functionalized nitrogen-containing tertiary phosphine oxides. Ab initio study of ortho-, meta- and para-dimethylphosphinylmethyleneoxyaniline

Vibrational study of new Pt(II) and Pd(II) complexes of functionalized nitrogen-containing tertiary phosphine oxides, namely ortho-, meta- and para-dimethylphosphinylmethyleneoxyaniline (o-, m- and p-dpmoa), (CH3)2P(O)CH2OC6H4NH2, have been presented. Geometry optimization of the ligands was performed at HF/6-31G* and B3LYP/6-31G* levels of the theory. Harmonic frequencies were calculated at HF/6-31G* optimized geometries. Relative gas-phase and solution-phase (H2O and CH3CN) basicities of o-, m- and p-dpmoa ligands have been determined by ab initio calculations at STO-3G level with the Onsager reaction field model. On the basis of the vibrational study, physical and analytical data it was suggested that the ligands in the complexes studied coordinate through the amino group and form square-planar platinum and palladium complexes of the general formula ML2Cl2 (M = Pt, Pd, L = o-, m- and p-dpmoa).     

Theoretical studies of host-guest interaction in gas hydrates

Ab initio calculations and atoms-in-molecules (AIM) analysis have been used to investigate the host-guest interaction in dodecahedral water cages using a variety of guest species that include monatomic (He, Ne, Ar, Kr, and Xe), diatomic (CO, H(2), N(2), O(2), and NO), triatomic (CO(2), NO(2), and O(3)), and polyatomic (CH(4) and NH(3)) molecules. Geometry optimization for the guest species, host cage, and their complexes was carried out using the second order M?ller-Plesset perturbation method with the 6-31G** basis set. Single point energy calculations using the same method but different basis sets (6-31++G**, 6-311++G**, aug-cc-pVDZ, and aug-cc-pVTZ) were carried out for the MP2/6-31G** optimized geometries. The interaction energy between the guest species and the host cage has been obtained in the complete basis set limit by basis set extrapolation.     

NH3…H2O体系氢键团簇结构的从头计算方案比较

对H2O，NH3单分子以及反式线性NH3…H2O双分子团簇应用HF，DFT，MR2，CCD ，QCISD等从头计算方法在不同的基组水平上进行几何结构优化和氨—水分子间相 互作用能计算。比较实验值和计算结果表明：对氨—水体系，MP2／6-31+G^**方法 能够满足氢键N-0键长误差为0．004nm，键角∠NOH误差为1°，NH3…H2O相互作用 能△E与CCD和QCISD偏差2．1kJ／mol计算精度的要求，并有较合理的计算代价。     

高能体系分子间相互作用研究: 含NNO~2和NH~2混合物

以abinitioHF/6-31G^*计算求得NH~3+NH~2NO~2的两种优化构型,经MP4电子相关能校正和Boys-Bernardi方案校正基组叠加误差求得精确的分子间相互作用能。还用PM3方法计算研究TATB(均三氨基三硝基苯)分别与HMX(奥克托金)和RDX(黑索金)的混合体系,经色散能校正电子相关近似地求得分子间相互作用能。结果表明,NH~3与NH~2NO~2之间的最大结合能为-38.32kJ/mol;分子间相互作用增强了N-NO~2键强度;TATB与HMX,RDX的结合能远大于石墨与HMX或RDX的结合能,表明TATB对HMX和RDX的润湿和钝感作用较石墨更强。     

Ab initio molecular orbital calculations of the infrared spectra of hydrogen-bonded complexes of water,ammonia and hydroxylamine

Summary The geometries of three hydrogen-bonded dimers of hydroxylamine have been optimized, at the MP2 level of theory, using the 6-31G** basis set. These calculations yielded three separate local minima on the dimer potential energy surface. The interaction energies of these three species have been calculated, and corrected for basis set superposition error. The infrared band wavenumbers and intensities have been computed, and the monomer-dimer wavenumber shifts and intensity enhancements rationalized in terms of the types and strengths of hydrogen bonds present. The predicted wavenumbers have been correlated with those measured in a recent matrix isolation spectroscopic study, and an argument for the structure of the preferred dimer has been presented.     

Ab initio Study on the Intermolecular Interaction and Thermodynamic Properties of Methyl Nitrate Dimer

Three stable dimers of methyl nitrate have been obtained and their geometries have been fully optimized at the HF/6‐31G,. level. Binding energies have been calculated with correction for the basis set superposition error (BSSE) and zero point energy (ZPE). The cyclic overlap‐type structure, the binding energy of which is 11.97 kJ/mol at the MP4SDTQ/6‐31G. / HF/6‐31G. level, is the most stable. No intermolecular hydrogen bond was found, and the charge transfer between two subsystems is minute. The thermodynamic properties of methyl nitrate and its dimers have been calculated based on the vibrational analysis and statistical thermodynamics.     

Conformational study of the structure of free 18-crown-6

A conformational search was performed for 18-crown-6 using the CONLEX method at the MM3 level. To have a more accurate energy order of the predicted conformations, the predicted conformations were geometry optimized at the HF/STO-3G level and the 198 lowest energy conformations, according to the HF/STO-3G energy order, were geometry optimized at the HF/6-31+G level. In addition, the 47 nonredundant lowest energy conformations, according to the MP2/6-31+G energy order at the HF/6-31+G optimized geometry, hereafter the MP2/6-31+G//HF/6-31+G energy order, were geometry optimized at the B3LYP/6-31+G level. According to the MP2/6-31+G//B3LYP/6-31+G energy order, three conformations had energies lower than the experimentally known Ci conformation of 18c6. At the MP2/6-31+G//B3LYP/6-31+G level, the S6 lowest energy conformation is more stable by 1.96 kcal/mol than this Ci conformation. This was confirmed by results at the MP2/6-31+G level with an energy difference of 1.84 kcal/mol. Comparison between the structure of the S6 conformation of 18c6 and the S4 lowest energy conformation of 12-crown-4, as well as other important conformations of both molecules, is made. It is concluded that the correlation energy is necessary to have an accurate energy order of the predicted conformations. A rationalization of the conformational energy order in terms of the hydrogen bonding and conformational dihedral angles is given. It is also suggested that to have a better energy order of the predicted conformations at the MM3 level, better empirical force fields corresponding to the hydrogen bond interactions are needed.     

叠氮化氢二聚体的分子间相互作用

通过ａｂ ｉｎｉｔｉｏＨＦ／６－３１Ｇ＊计算求得叠氮化氢二聚体（ＨＮ３）２势能面上四种优化结构型,经ＭＰ４ＳＤＴＱ电子相关校正和基组叠加误差２校正求得这些构型下的分子间相互作用能。结果表明,在四种优化构型中Ｎ３Ｈ．．．ＮＨＮ２最稳定,其分子间相互作用能为－１６．０７ＫＪ．ｍｏｌ＾－１。     

CH3F分子间相互作用的ab initio研究

在MP2/6-311+ +G(3d,3p)电子相关校正水平上,对CH₃F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证,发现3种势能面上有极小点的构型.进一步在高级电子相关校正的MP4S-DTQ、CCSD(T)/6-311+ +G(3df,3pd)方法水平上,对其中总能量最小的构型进行精确计算,得到二聚体的结合能为-9.707kJ/mol.研究结果支持了由光谱实验结果推测的构型,解释了CH₃F二聚体的谐振频率的多样性.     

Cl+CH3OCl反应机理的理论研究

采用从头算方法,对多通道反应体系Cl+CH3OCl的反应机理进行了理论研究.在MP2/6-31+G(d,p)水平下优化了反应物、络合物、产物和过渡态的几何构型,并对得到的平衡几何构型进行了简谐振动频率分析.在相同水平下以过渡态为出发点,通过内禀反应坐标(IRC)理论计算了反应的最小能量路径.并且在MC-QCISD(T)/6-31G(d)高水平下进行了单点能量校正.研究结果表明,该反应存在4条可行的反应通道,其中生成HCl和CH2OCl的通道为主反应通道,其他反应通道为次反应通道.     

Quantum chemical studies on structure, conformation and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene

The molecular structure, conformational stability and isomerization of nitroso, nitro substituted benzene and 1,3-cyclopentadiene in gas phase have been investigated using ab initio and density functional theory methods. The molecular geometries and energetics of possible conformers were obtained by employing MP2, B3LYP and B3PW91 levels of theory implementing 6-31G* basis set. The relative stabilities of the conformations were evaluated from the energy differences of the structure. Chemical hardness (η) and chemical potential (μ) were calculated at HF/6-31G* level of theory for all the positional and geometrical isomers to study the maximum hardness principle. Each optimized structure has been tested against the imaginary frequencies at MP2/6-31G* level of theory in order to be sure they are located at energy minimum.     

甲醛与乙腈相互作用的理论研究

在HF/6-311G(d,p)、 MP2/6-311G(d,p)和B3LYP/6-311G(d,p)水平上,对H2CO和CH3CN以及设计的4种结构H2CO…CH3CN复合物等进行几何全优化和振动频率计算,排除振动频率为负值的非局域极小点结构,并对稳定的环状构型复合物结合能进行基组重叠误差校正和零点振动能校正.分子间相互作用的能量分解分析显示,静电能在H2CO...CH3CN相互作用能量中占主导地位,电荷转移能居第二位.