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.     

带电组氨酸侧链与DNA碱基间非键作用强度的理论研究

采用MP2方法和6-31+G(d,p)基组优化得到了带有一个正电荷的组氨酸侧链与4个DNA碱基间形成的18个氢键复合物的气相稳定结构, 从文献中获取了组氨酸侧链与DNA碱基间形成的12个堆积和T型复合物的气相稳定结构, 使用包含基组重叠误差(BSSE)校正的MP2方法和aug-cc-pVTZ基组及密度泛函理论M06-2X-D3方法和aug-cc-pVDZ基组计算了这些复合物的结合能. 研究结果表明, 包含BSSE校正的M06-2X-D3方法和aug-cc-pVDZ基组能够给出较准确的结合能; 气相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用明显强于堆积作用和T型作用, 组氨酸侧链最易通过离子氢键与胞嘧啶C和鸟嘌呤G作用形成氢键复合物, 组氨酸与胞嘧啶C和鸟嘌呤G间的T型作用强于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用; 水相条件下, 组氨酸侧链与同种DNA碱基间的离子氢键作用仍明显强于堆积作用和T型作用, 组氨酸侧链更易与胞嘧啶C和鸟嘌呤G相互作用形成氢键复合物, 但是最强的组氨酸侧链与胞嘧啶C间的T型作用明显弱于与腺嘌呤A和胸腺嘧啶T间的离子氢键作用, 说明水相条件下组氨酸侧链与DNA碱基间主要通过离子氢键作用形成氢键复合物.     

Interaction of anions with perfluoro aromatic compounds

The complexes formed by a variety of anions with perfluoro derivatives of benzene, naphthalene, pyridine, thiophene, and furan have been calculated using DFT (B3LYP/6-31++G**) and MP2 (MP2/6-31++G** and MP2/6-311++G**) ab initio methods. The minimum structures show the anion interacting with the pi-cloud of the aromatic compounds. The interaction energies obtained range between -8 and -19 kcal mol(-1). The results obtained at the MP2/6-31++G** and MP2/6-311++G** levels are similar. However, the B3LYP/6-31++G** results provide longer interaction distances and smaller interaction energies than do the MP2 results. The interaction energies have been partitioned using an electrostatic, polarization, and van der Waals scheme. The AIM analysis of the electron density shows a variety of topologies depending on the aromatic system considered.     

Calculated structures and relative stabilities of furoxan,some 1,2-dinitrosoethylenes and other isomers

The structures and relative stabilities of furoxan and some of its isomers, e.g., the 1,2-dinitrosoethylenes, have been determined by means of ab initio Hartee–Fock and Møller–Plesset calculations. Geometries were optimized at the HF/3-21G, HF/6-31G* and MP2/6-31G* levels, and subsequently used for computing MP2/6-31G*, MP3/6-31G*, and MP4/6-31G* energies. The results are markedly affected by the inclusion of electronic correlation, which renders three of the isomers unstable. It also emphasizes the importance of a zwitterionic contribution to the structure of furoxan, which promotes ring-opening through a cis 1,2-dinitrosoethylene intermediate/transition state that has an MP4/6-31G*//MP2/6-31G* energy that is 31.6 kcal/mol above furoxan.     

Binding energies of the proton-bound amino Acid dimers gly.gly, ala.ala, gly.ala, and lys.lys measured by blackbody infrared radiative dissociation

Arrhenius activation energies in the zero-pressure limit for dissociation of gas-phase proton-bound homodimers of N,N-dimethylacetamide (N,N-DMA), glycine, alanine, and lysine and the heterodimer alanine.glycine were measured using blackbody infrared radiative dissociation (BIRD). In combination with master equation modeling of the kinetic data, binding energies of these dimers were determined. A value of 1.25 +/- 0.05 eV is obtained for N,N-DMA and is in excellent agreement with that reported in the literature. The value obtained from the truncated Boltzmann model is significantly higher, indicating that the assumptions of this model do not apply to these ions. This is due to the competitive rates of photon emission and dissociation for these relatively large ions. The binding energies of the amino acid dimers are ~1.15 +/- 0.05 eV and are indistinguishable despite the difference in their gas-phase basicity and structure. The threshold dissociation energies can be accurately modeled using a range of dissociation parameters and absorption/emission rates. However, the absolute values of the dissociation rates depend more strongly on the absorption/emission rates. For N,N-DMA and glycine, an accurate fit was obtained using frequencies and transition dipole moments calculated at the ab initio RHF/2-31G* and MP2/2-31G* level, respectively. In order to obtain a similar accuracy using values obtained from AM1 semiempirical calculations, it was necessary to multiply the transition dipole moments by a factor of 3. These results demonstrate that in combination with master equation modeling, BIRD can be used to obtain accurate threshold dissociation energies of relatively small ions of biological interest.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。     

Effects of extending the computational model on DNA-protein T-shaped interactions: the case of adenine-histidine dimers

The MP2/6-31G*(0.25) π-π or π(+)-π T-shaped (edge-to-face) interactions between neutral or protonated histidine and adenine were considered using computational models of varying size to determine the effects of the protein and DNA backbones on the preferred dimer structure and binding strength. The overall consequences of the backbones are reasonably subtle for the neutral adenine-histidine T-shaped dimers. Furthermore, the minor changes in the binding strengths of these dimers upon model extension arise from additional (attractive) backbone-π (bb-π) contacts and changes in the preferred π-π orientations, which is verified by the quantum theory of atoms in molecules (QTAIM). Since the binding strength of the extended dimer equals the sum of the individual backbone-π and π-π contributions, the π-π component is not appreciably affected by polarization of the ring upon inclusion of the biological backbone. In contrast, the larger effect of the backbone on the protonated histidine dimers cannot simply be predicted as the sum of changes in the π-π and bb-π components regardless of the dimer type or model. This suggests, and QTAIM qualitatively supports, that the magnitude of the π(+)-π contribution changes, which is likely due to alterations in the electrostatic landscape of the monomer rings upon inclusion of the biological backbone that largely affect T-shaped dimers. These findings differ from those previously reported for (neutral) π-π stacked and (metallic) cation-π interactions, which highlights the distinct properties of each (π-π, π(+)-π, and cation-π) classification of noncovalent interaction. Furthermore, these results emphasize the importance of considering backbone-π interactions when analyzing contacts that appear in experimental crystal structures and cautions the use of truncated models when evaluating the magnitude of the π(+)-π contribution present in large biological complexes.     

A new method for empirical force field calculations on localized and delocalized carbocations

A new empirical force field method for localized and delocalized carbocations is described. Additional geometry parameters for carbocations were added to Allinger's MMP2 molecular mechanics program, which can treat delocalized π-systems. The effect of hyperconjugation in carbocations is introduced via a quantum chemical term into force field calculations for the first time. The calculated heats of formation are in excellent agreement with a wide range of experimental data; the largest deviations are about 3.5 kcal/mol. The calculated structures agree very well with those computed at correlated ab initio levels (MP2(full)/6-31G*). The relative energies and geometries of different conformations of representative carbocations also were in good agreement with MP4/6-31G*//MP2(full)/6-31G* results. © 1996 by John Wiley & Sons, Inc.     

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.     

Coordination of methanol clusters to benzene: a computational study

Benzene-methanol cluster structures were investigated with theoretical chemistry methods to describe the microsolvation of benzene and the benzene-methanol azeotrope. Benzene-methanol (MeOH) clusters containing up to six methanol molecules have been calculated by ab initio [MP2/6-311++G(d,p)//MP2/6-31+G(d,p) + BSSE correction] method. The BSSE was found quite large with this basis set, hence, different extrapolation schemes in combination with the aug-cc-pVxZ basis sets have been used to estimate the complete basis set limit of the MP2 interaction energy [ΔE(MP2/CBS)]. For smaller clusters, n ≤ 3, DFT procedures (DFTB+, MPWB1K, M06-2X) have also been applied. Geometries obtained for these clusters by M06-2X and MP2 calculations are quite similar. Based on the MP2/CBS results, the most stable C(6)H(6)(MeOH)(3) cluster is characterized by a hydrogen bonded MeOH trimer chain interacting with benzene via π···H-O and O···H-C(benzene) hydrogen bonds. Larger benzene-MeOH clusters with n ≥ 4 consist of cyclic (MeOH)(n) subclusters interacting with benzene by dispersive forces, to be denoted by C(6)H(6) + (MeOH)(n). Interaction energies and cooperativity effects are discussed in comparison with methanol clusters. Besides MP2/CBS calculations, for selected larger clusters the M06-2X/6-311++G(d,p)//M06-2X/6-31+G(d,p) procedure including the BSSE correction was also used. Interaction energies obtained thereby are usually close to the MP2/CBS limit. To model the benzene-MeOH azeotrope, several structures for (C(6)H(6))(2)(MeOH)(3) clusters have been calculated. The most stable structures contain a tilted T-shaped benzene dimer interacting by π···H-O and O···H-C (benzene) hydrogen bonds with a (MeOH)(3) chain. A slightly less negative interaction energy results for a parallel displaced benzene sandwich dimer with a (MeOH)(3) chain atop of one of the benzene molecules.     

Mixed Aggregates of 1-Methoxyallenyllithium with Lithium Chloride

A combined computational and ¹³C NMR study was used to investigate the formation of mixed aggregates of 1-methoxyallenyllithium and lithium chloride in tetrahydrofuran (THF) solution. The observed and calculated chemical shifts, as well as the calculated free energies of mixed aggregate formation (MP2/6-31+G(d)), are consistent with the formation of a mixed dimer as the major species in solution. Free energies of mixed dimer, trimer, and tetramer formation were calculated by using the B3LYP and MP2 methods and the 6-31+G(d) basis set. The two methods generated different predictions of which mixed aggregates will be formed, with B3LYP/6-31+G(d) favoring mixed trimers and tetramers in THF solution, and MP2/6-31+G(d) favoring mixed dimers. Formation of the sterically unhindered mixed dimers is also consistent with the enhanced reactivity of these compounds in the presence of lithium chloride. The spectra are also consistent with some residual 1-methoxyallenyllithium tetramer, as well as small amounts of higher mixed aggregates. Although neither computational method is perfect, for this particular system, the calculated free energies derived using the MP2 method are in better agreement with experimental data than those derived using the B3LYP method.     

Estimation on the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides

Computation of accurate intramolecular hydrogen-bonding energies for peptides is of great importance in understanding the conformational stabilities of peptides and developing a more accurate force field for proteins. We have proposed a method to determine the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies in glycine and alanine peptides. In this article, the method is further applied to evaluate the intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies in peptides. The optimal structures of the intramolecular 10-membered ring N-H...O=C hydrogen bonds in glycine and alanine tripetide molecules are obtained at the MP2 level with 6-31G(d), 6-31G(d,p), and 6-31+G(d,p) basis sets. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are then evaluated based on our method at the MP2/6-311++G(3df,2p) level with basis set superposition error correction. The intramolecular 10-membered ring N-H...O=C hydrogen-bonding energies are calculated to be in the range of -6.84 to -7.66, -4.44 to -4.98, and -6.95 to -7.88 kcal/mol. The method is also applied to estimate the individual intermolecular hydrogen-bonding energies in the dimers of amino-acetaldehyde, 2-amino-acetamide, formamide, and oxalamide, each dimer having two identical intermolecular hydrogen bonds. According to our method, the individual intermolecular hydrogen-bonding energies in the four dimers are calculated to be -1.77, -1.67, -6.35, and -4.82 kcal/mol at the MP2/6-311++G(d,p) level, which are in good agreement with the values of -1.84, -1.72, -6.23, and -4.93 kcal/mol predicted by the supermolecular method.     

Investigation of the use of B3LYP zero-point energies and geometries in the calculation of enthalpies of formation

The use of B3LYP/6–31G^* zero-point energies and geometries in the calculation of enthalpies of formation has been investigated for the enlarged G2 test set of 148 molecules [J. Chem. Phys. 106 (1997) 1063]. A scale factor of 0.96 for the B3LYP zero-point energies gives an average absolute deviation nearly the same as scaled HF/6–31G^* zero-point energies for G2, G2(MP2), and B3LYP/6–311 + G(3df,2p) enthalpies. A scale factor of 0.98, which has been recommended in some studies, increases the average absolute deviation by about 0.2 kcal/mol. Geometries from B3LYP/6–31G^* are found to do as well as MP2/6–31G^* geometries in the calculation of the enthalpies of formation.     

Theoretical study of hydrogen bonding in homodimers and heterodimers of amide, boronic acid, and carboxylic acid, free and in encapsulation complexes

The homodimers and the heterodimers of two amides, two boronic acids, and two carboxylic acids have been calculated in the gas phase and in N,N-dimethylformamide (DMF) and CCl(4) solvents using the DFT (M06-2X and M06-L) and the MP2 methods in conjunction with the 6-31G(d,p) and 6-311+G(d,p) basis sets. Furthermore, their pairwise coencapsulation was studied to examine its effect on the calculated properties of the hydrogen bonds at the ONIOM[M06-2X/6-31G(d,p);PM6], ONIOM[MP2/6-31G(d,p); PM6], and M06-2X/6-31G(d,p) levels of theory. The present work is directed toward the theoretical rationalization and interpretation of recent experimental results on hydrogen bonding in encaptulation complexes [D. Ajami et al. J. Am. Chem. Soc. 2011, 133, 9689-9691]. The calculated dimerization energy (ΔE) values range from 0.74 to 0.35 eV for the different dimers in the gas phase, with the ordering carboxylic homodimers > amide-carboxylic dimers > amide homodimers > boronic-carboxylic dimers > amide-boronic dimers > boronic homodimers. In solvents, generally smaller ΔE values are calculated with only small variations in the ordering. In the capsule, the ΔE values range between 0.67 and 0.33 eV with practically the same ordering as in the gas phase. The calculated % distributions of the encapsulated dimers, taking into account statistical factors, are in agreement with the experimental distribution, where the occurrence of boronic homodimer dominates, even though it is calculated to have the smallest ΔE.     

Ab initio study of hydrogen-bond formation between aliphatic and phenolic hydroxy groups and selected amino acid side chains

Hydrogen bonding was studied in 24 pairs of isopropyl alcohol and phenol as one partner, and water and amino-acid mimics (methanol, acetamide, neutral and protonated imidazole, protonated methylalamine, methyl-guanidium cation, and acetate anion) as the other partner. MP2/6-31+G* and MP2/aug-cc-pvtz calculations were conducted in the gas phase and in a model continuum dielectric environment with dielectric constant of 15.0. Structures were optimized in the gas phase with both basis sets, and zero-point energies were calculated at the MP2/6-31+G* level. At the MP2/aug-cc-pvtz level, the BSSE values from the Boys-Bernardi counterpoise calculations amount to 10-20 and 5-10% of the uncorrected binding energies of the neutral and ionic complexes, respectively. The geometry distortion energy upon hydrogen-bond formation is up to 2 kcal/mol, with the exception of the most strongly bound complexes. The BSSE-corrected MP2/aug-cc-pvtz binding energy of -27.56 kcal/mol for the gas-phase acetate...phenol system has been classified as a short and strong hydrogen bond (SSHB). The CH3NH3+...isopropyl alcohol complex with binding energy of -22.54 kcal/mol approaches this classification. The complete basis set limit (CBS) for the binding energy was calculated for twelve and six complexes on the basis of standard and counterpoise-corrected geometry optimizations, respectively. The X...Y distances of the X-H...Y bridges differ by up to 0.03 A as calculated by the two methods, whereas the corresponding CBS energy values differ by up to 0.03 kcal/mol. Uncorrected MP2/aug-cc-pvtz hydrogen-bonding energies are more negative by up to 0.35 kcal/mol than the MP2/CBS values, and overestimate the CCSD(T)/CBS binding energies generally by up to 5% for the eight studied complexes in the gas phase. The uncorrected MP2/aug-cc-pvtz binding energies decreased (in absolute value) by 11-18 kcal/mol for the ionic species and by up to 5 kcal/mol for the neutral complexes when the electrostatic effect of a polarizable model environment was considered. The DeltaECCSD(T) - DeltaEMP2 corrections still remained close to their gas-phase values for four complexes with 0, +/-1 net charges. Good correlations (R2 = 0.918-0.958) for the in-environment MP2/aug-cc-pvtz and MP2/6-31+G* hydrogen-bonding energies facilitate the high-level prediction of these energies on the basis of relatively simple MP2/6-31+G* calculations.     

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.     

Comparison of levels of electronic structure theory in direct dynamics simulations of C2H5F --> HF + C2H4 product energy partitioning

Direct dynamics simulations at the MP2/6-311++G** level of theory were performed to study C(2)H(5)F --> HF + C(2)H(4) product energy partitioning. The simulation results are compared with experiment and a previous MP2/6-31G* simulation. The current simulation with the larger basis set releases more energy to HF vibration and less to HF + C(2)H(4) relative translation as compared to the previous simulation with the 6-31G* basis set. The HF rotation and vibration energy distributions determined from the current simulation are in overall very good agreement with previous experimental studies of C(2)H(5)F dissociation by chemical activation and IRMPA. A comparison of the simulations with experiments suggests there may be important mass effects for energy partitioning in HX elimination from haloalkanes. The transition state (TS) structures and energies calculated with MP2 and the 6-31G* and 6-311++G** basis sets are compared with those calculated using CCD, CCSD, CCSD(T), and the 6-311++G** basis set.     

乙醛二聚体内蓝移型氢键的理论研究

利用理论方法研究了乙醛二聚体内的氢键. 在MP2/6-31＋G(d), B3LYP/6-31＋G(d), B3LYP/6-311＋＋G(d,p)和B3LYP/6-311＋＋G(3df,2p)水平上, 利用常规方法和均衡校正方法对3种稳定的乙醛二聚体进行了几何优化和振动频率计算. 计算结果表明: 在二聚体A和C中乙醛中C—H键强烈收缩, 存在显著的C—H…O蓝移型氢键. 自然键轨道(NBO)分析表明, 电子供体轨道和电子受体轨道之间相互作用的稳定化能、分子内电子密度重排、轨道再杂化和结构重组是决定氢键红移和蓝移的主要因素. 其中, 轨道间稳定化能属于键伸长效应, 分子内电子密度重排、轨道再杂化和电子受体内部结构重组属于键收缩效应. 在二聚体A和C中, 由于键收缩效应处于优势地位导致C—H…O蓝移氢键存在.     

乙烯、乙炔与双卤分子间π型卤键的电子密度拓扑研究

运用DFT和MP2(full)在6-311++G(d, p)和aug-cc-pvdz基组水平上, 对一系列简单的分子间π型卤键体系C2H4(C2H2)-XY(XY= F2、Cl2、Br2、ClF、BrF、BrCl) 进行构型全优化, 得到了T型卤键复合物. 结果表明MP2(full)/ 6-311++G(d, p)计算结果与实验结果较吻合. 并在MP2水平上计算了分子间的相互作用能, 用标准Counterpoise procedure (CP)方法对基函数迭加误差(BSSE)进行了校正. 利用电子密度拓扑分析方法对卤键复合物的拓扑性质进行了分析研究.