Ab initio base-pairing energies of an oxidized thymine product, 5-formyluracil, with standard DNA bases at the BSSE-free DFT and MP2 theory levels

Oxidation of the thymine methyl group produces two stable products, non-mutagenic 5-hydroxymethyluracil and highly mutagenic 5-formyluracil. We have calculated the interaction energy of base-pair formation involving 5-formyluracil bound to the natural DNA bases adenine (A), cytosine (C), guanine (G), and thymine (T), and discuss the effects of the 5-formyl group with respect to similar base-pairs containing uracil, 5-hydroxyuracil, thymine (5-methyluracil), and 5-hydroxycytosine. The interaction geometries and energies were calculated four ways: (a) using density functional theory (DFT) without basis set super-position error (BSSE) corrections, (b) using DFT with BSSE correction of geometries and energies, (c) using M?ller-Plesset second order perturbation theory (MP2) without BSSE correction, and (d) using MP2 with BSSE geometry and energy correction. All calculations used the 6-311G(d,p) basis set. Notably, we find that the A:5-formyluracil base-pair is more stable than the precursor A:T base-pair. The relative order of base-pair stabilities is A:5-Fo-U > G:5-Fo-U > C:5-Fo-U > T:5-Fo-U.     

DNA甲基化-非甲基化碱基间堆积作用的理论研究

运用二级Mфller-Plesset(MP2)理论方法和cc-pVDZ基组优化了6-甲基鸟嘌呤(O6-MethylG),4-甲基胸腺嘧啶(O4-MethylT)以及5-甲基胞嘧啶(C5-MethylC)与DNA碱基鸟嘌呤(G),腺嘌呤(A),胞嘧啶(C),胸腺嘧啶(T)之间的堆积构型.在MP2/aug-cc-pVXZ//MP2/cc-pVDZ(X=D,T)水平上,采用完全基组外推方法校正了堆积碱基对间的相互作用能,并用完全均衡校正法(CP)校正了基组重叠误差(BSSE).MP2计算结果表明,DNA碱基甲基化使得嘧啶-嘧啶、嘧啶-嘌呤堆积碱基间的平行旋转角发生明显改变,并使堆积碱基间的相互作用能增大.在MP2/cc-pVDZ计算级别上得到了各堆积碱基对的全电子波函数,并用分子中的原子理论(AIM)分析了堆积碱基对间的弱相互作用.AIM分析结果显示,甲基化增强了堆积碱基间的π-π作用,且甲基氢与相邻碱基间形成H2CH…X(X=O,N,CH3,NH2)等类型的氢键.甲基化损伤使碱基间重叠程度增大、π-π作用增强以及堆积碱基间形成多个氢键,是堆积作用能增加的主要原因.     

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.     

MP2 study on the hydrogen bonding interaction between 5-hydroxy-5-methylhydantoin and DNA bases: A,C, G,T

The 5-hydroxy-5-methylhydantoin (5-OH-5-Me-dHyd) is a nucleobase lesion induced by the action of ionizing radiation on thymine residue in DNA. In this study, we present the hydrogen bonding base pairs involving 5-OH-5-Me-dHyd bound to the four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys–Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs were characterized using NBO analysis and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: 5-OH-5-Me-dHyd:G>5-OH-5-Me-dHyd:A>5-OH-5-Me -dHyd:C~5-OH-5-Me-dHyd:T.     

羰基镍簇Ni(CO)n(n=1-4)的结构和Ni-CO键解离性质的密度泛函理论研究

在密度泛函理论框架下, 应用不同泛函计算了配合物Ni(CO)n(n=1~4)的平衡几何构型和振动频率. 考察了泛函和基组重叠误差对预测Ni—CO键解离能的影响. 计算结果表明, 用杂化泛函能得到与实验一致的优化几何构型和较合理的振动频率. 对Ni(CO)n(n=2~4)体系, 用“纯”泛函, 如BP86和BPW91, 可得到与CCSD(T)更符合、 并与实验值接近的解离能. 当解离产物出现单个金属原子或离子(如金属羰基配合物的完全解离)时, BSSE校正项的计算中应保持金属部分的电子结构一致. 只有考虑配体基组和不考虑配体基组两种情况下金属的电子构型与配合物中金属的构型一致时, 才能得到合理的BSSE校正, 从而预测合理的解离能.     

The effect of oxidation on the stability of G:C base pair: a MP2 study

The base-pairing energies of eleven oxidized G:C base pairs were characterized by Møller–Plesset perturbation theory. The analysis was focused on the base pairs consisting of guanine and one of the following cytosine derivatives: G:C1 (5-hydroxycytosine); G:C2 (5-hydroxyuracil); G:C3 (5,6-dihydroxy-cytosine); G:C4 (5,6-dihydroxy-uracil); G:C5 (cytosine glycol);G:C6 (isodialuric acid); G:C7 (uracil glycol), and the base pairs between cytosine and one of the following guanine derivatives: G8:C (8-oxo-guanine); G9:C (6-enol-8-keto-guanine); G10:C (xantho-sine); and G11:C (8-hydroxy-guanine). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys-Bernardi counterpoise correction scheme. The results obtained show that the interaction energies of the base pairs decrease in the following order: G8:C > G:C5 ~ G:C1 > G:C3 ~ G:C ~ G11:C > G:C2 > G:C4 > G:C7 > G:C6 > G10:C > G9:C.     

Rapid evaluation of the binding energies between peptide amide and DNA base

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 20 hydrogen‐bonded amide–base dimers are evaluated from the analytic potential energy function established in our laboratory recently. The analytic potential energy function is used to calculate the N? H···N, N? H···O?C, C? H···N, and C? H···O?C dipole–dipole attractive interaction energies and C?O···O?C, N? H···H? N, and N? H···H? C dipole–dipole repulsive interaction energies in the 20 dimers composed of DNA bases adenine, guanine, cytosine, or thymine and peptide amide. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the basis set superposition error (BSSE) correction. For all the 20 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 19 dimers, only one difference is larger than 0.50 kcal/mol and the difference is only 0.61 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.030 Å for all the 20 dimers. The analytic potential energy function is further applied to four more complicated DNA base‐peptide amide systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded peptide amide–DNA base dimers quickly and accurately. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Accurate ab initio binding energies of alkaline earth metal clusters

The effects of basis set superposition error (BSSE) and core-correlation on the electronic binding energies of alkaline earth metal clusters Y(n) (Y = Be, Mg, Ca; n = 2-4) at the Moller-Plesset second-order perturbation theory (MP2) and the single and double coupled cluster method with perturbative triples correction (CCSD(T)) levels are examined using the correlation consistent basis sets cc-pVXZ and cc-pCVXZ (X = D, T, Q, 5). It is found that, while BSSE has a negligible effect for valence-electron-only-correlated calculations for most basis sets, its magnitude becomes more pronounced for all-electron-correlated calculations, including core electrons. By utilizing the negligible effect of BSSE on the binding energies for valence-electron-only-correlated calculations, in combination with the negligible core-correlation effect at the CCSD(T) level, accurate binding energies of these clusters up to pentamers (octamers in the case of the Be clusters) are estimated via the basis set extrapolation of ab initio CCSD(T) correlation energies of the monomer and cluster with only the cc-pVDZ and cc-pVTZ sets, using the basis set and correlation-dependent extrapolation formula recently devised. A comparison between the CCSD(T) and density functional theory (DFT) binding energies is made to identify the most appropriate DFT method for the study of these clusters.     

A theoretical study of nonlinearity in heterocyclic hydrogen-bonded complexes

Ab initio calculations are performed at the MP2/6-311++G(d,p) and DFT/B3LYP/6-311++G(d,p) theoretical levels to obtain geometries, H-bond energies and harmonic infrared vibrational properties for the Cs symmetry structures of heterocyclic hydrogen-bonded complexes, CnHmY-HX. The H-bond lengths in DFT/B3LYP calculation level are in better agreement with the experimental values than the MP2 results. The geometry optimization are interpreted in terms of hydrogen bond nonlinearity represented by theta; and phi angles, once the hydrogen bond is formed among n-electrons pairs of the heteroatom in heterocyclic and the hydrogen atom in HX. The hydrogen bond energy after of the zero-point vibrational energy (ZPE) and basis set superposition error (BSSE) corrections are overestimated at DFT/B3LYP, whereas the MP2 BSSE corrections are very large than corresponding DFT/B3LYP. For example, the BSSE corrections for the C2H4S-HNC complex are 7.60 and 0.09 kJ mol(-1) in MP2 and DFT/B3LYP calculations levels, respectively. The new vibrational modes in infrared harmonic spectrum arising from complexation show several interesting features, especially the intermolecular stretching mode.     

MP2 study on the hydrogen-bonding interaction between 5-hydroxymethyl-uracil and DNA bases: A,C, G,T

The 5-hydroxymethyl-uracil (HmU) is a product of oxidative attack on the methyl group of thymine in DNA. In this work, we present the hydrogen bonding complexes formation involving HmU bound to the four bases in DNA: adenine (A), cytosine (C), guanine (G), and thymine (T). Full geometry optimizations have been performed for the studied complexes by MP2 method. The interaction energies were corrected for the basis-set superposition error (BSSE), using the full Boys-Bernardi counterpoise correction scheme. Hydrogen bonding patterns of these base pairs were characterized using NBO analysis and AIM analysis. According to the calculated binding energies and structural parameters, the stability of the base pairs decrease in the following order: HmU:A > HmU:G > HmU:C > HmU:T.     

Rapid evaluation of the binding energies in hydrogen-bonded amide-thymine and amide-uracil dimers in gas phase

The binding energies and the equilibrium hydrogen bond distances as well as the potential energy curves of 48 hydrogen‐bonded amide–thymine and amide–uracil dimers are evaluated from the analytic potential energy function established in our lab recently. The calculation results show that the potential energy curves obtained from the analytic potential energy function are in good agreement with those obtained from MP2/6‐311+G** calculations by including the BSSE correction. For all the 48 dimers, the analytic potential energy function yields the binding energies of the MP2/6‐311+G** with BSSE correction within the error limits of 0.50 kcal/mol for 46 dimers, only two differences are larger than 0.50 kcal/mol and the largest one is only 0.60 kcal/mol. The analytic potential energy function produces the equilibrium hydrogen bond distances of the MP2/6‐311+G** with BSSE correction within the error limits of 0.050 Å for all the 48 dimers. The analytic potential energy function is further applied to four more complicated hydrogen‐bonded amide–base systems involving amino acid side chain and β‐sheet. The values of the binding energies and equilibrium hydrogen bond distances obtained from the analytic potential energy function are also in good agreement with those obtained from MP2 calculations with the BSSE correction. These results demonstrate that the analytic potential energy function can be used to evaluate the binding energies in hydrogen‐bonded amide–base dimers quickly and accurately. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Base stacking in cytosine dimer. A comparison of correlated ab initio calculations with three empirical potential models and density functional theory calculations

Ab initio MP2/6-31G* interaction energies were calculated for more than 80 geometries of stacked cytosine dimer. Diffuse polarization functions were used to properly cover the dispersion energy. The results of ab initio calculations were compared with those obtained from three electrostatic empirical potential models, constructed as the sum of a Lennard-Jones potential (covering dispersion and repulsion contributions) and the electrostatic term. Point charges and point multipoles of the electrostatic term were also obtained at the MP2/6-31G* level of theory. The point charge MEP model (atomic charges derived from molecular electrostatic potential) satisfactorily reproduced the ab initio data. Addition of π-charges localized below and above the cytosine plane did not affect the calculated energies. The model employing the distributed multipole analysis gave worse agreement with the ab initio data than the MEP approach. The MP2 MEP charges were also derived using larger sets of atomic orbitals: cc-pVDZ, 6-311 + G(2d, p), and aug-cc-pVDZ. Differences between interaction energies calculated using these three sets of point charges and the MP2/6-31G* charges were smaller than 0.8 kcal/mol. The correlated ab initio calculations were also compared with the density functional theory (DFT) method. DFT calculations well reproduced the electrostatic part of interaction energy. They also covered some nonelectrostatic short-range effects which were not reproduced by the empirical potentials. The DFT method does not include the dispersion energy. This energy, approximated by an empirical term, was therefore added to the DFT interaction energy. The resulting interaction energy exhibited an artifact secondary minimum for a 3.9-4.0 vertical separation of bases. This defect is inherent in the DFT functionals, because it is not observed for the Hartree-Fock + dispersion interaction energy.© 1996 John Wiley & Sons, Inc.     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

A theoretical study on three conformational structures of 2,6‐distyrylpyridine

Ab initio methods at the levels HF/cc‐pVDZ, HF/6‐31G(d,p), MP2/cc‐pVDZ, and MP2/6‐31G(d,p), as well as methods based on density functional theory (DFT) employing the hybrid functional B3LYP with the basis sets cc‐pVDZ and 6‐31G(d,p), have been applied to study the conformers of 2,6‐distyrylpyridine. Bond distances, bond angles, and dihedral angles have been calculated at the B3LYP level. The calculated values were in good agreement with those measured by X‐ray diffraction analysis of 2,6‐distyrylpyridine. The values calculated using the Hartree‐Fock method and second‐order perturbation theory (MP2) were inconsistent. The optimized lowest‐energy geometries were calculated from the reported X‐ray structural data by the B3LYP/cc‐pVDZ method. Three conformations, A, B, and C, were proposed for 2,6‐distyrylpyridine. Calculations at the three levels of theory indicated that conformation A was the most stable structure, with conformations C and B being higher in energy by 1.10 and 2.57 kcal/mol, respectively, using the same method and basis function. The same trend in the relative energies of the three possible conformations was observed at the two levels of theory and with the different basis sets employed. The reported X‐ray data were utilized to optimize total molecular energy of conformation A at the different calculation levels. The bond lengths, bond angles, and dihedral angles were then obtained from the optimized geometries by ab initio methods and by applying DFT using the two basis functions cc‐pVDZ and 6‐31G(d,p). The values were analyzed and compared. The calculated total energies, the relative energies of the molecular orbitals, the gap between them, and the dipole moment for each conformational structure proposed for 2,6‐distyrylpyridine are also reported. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Ab initio and analytic intermolecular potentials for Ar-CF4

Ab initio calculations at the CCSD(T) level of theory were performed to characterize the Ar + CF4 intermolecular potential. Potential energy curves were calculated with the aug-cc-pVTZ basis set, and with and without a correction for basis set superposition error (BSSE). Additional calculations were performed with other correlation consistent basis sets to extrapolate the Ar-CF4 potential energy minimum to the complete basis set (CBS) limit. Both the size of the basis set and BSSE have substantial effects on the Ar + CF4 potential. Calculations with the aug-cc-pVTZ basis set, and with a BSSE correction, appear to give a good representation of the BSSE corrected potential at the CBS limit. In addition, MP2 theory is found to give potential energies in very good agreement with those determined by the much higher level CCSD(T) theory. Two model analytic potential energy functions were determined for Ar + CF4. One is a fit to the aug-cc-pVTZ calculations with a BSSE correction. The second was derived by fitting an average BSSE corrected potential, which is an average of the CCSD(T)/aug-cc-pVTZ potentials with and without a BSSE correction. These analytic functions are written as a sum of two-body potentials and excellent fits to the ab initio potentials are obtained by representing each two-body interaction as a Buckingham potential.     

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.     

Explicit correlation and basis set superposition error: the structure and energy of carbon dioxide dimer

We have investigated the slipped parallel and t-shaped structures of carbon dioxide dimer [(CO(2))(2)] using both conventional and explicitly correlated coupled cluster methods, inclusive and exclusive of counterpoise (CP) correction. We have determined the geometry of both structures with conventional coupled cluster singles doubles and perturbative triples theory [CCSD(T)] and explicitly correlated cluster singles doubles and perturbative triples theory [CCSD(T)-F12b] at the complete basis set (CBS) limits using custom optimization routines. Consistent with previous investigations, we find that the slipped parallel structure corresponds to the global minimum and is 1.09 kJ mol(-1) lower in energy. For a given cardinal number, the optimized geometries and interaction energies of (CO(2))(2) obtained with the explicitly correlated CCSD(T)-F12b method are closer to the CBS limit than the corresponding conventional CCSD(T) results. Furthermore, the magnitude of basis set superposition error (BSSE) in the CCSD(T)-F12b optimized geometries and interaction energies is appreciably smaller than the magnitude of BSSE in the conventional CCSD(T) results. We decompose the CCSD(T) and CCSD(T)-F12b interaction energies into the constituent HF or HF CABS, CCSD or CCSD-F12b, and (T) contributions. We find that the complementary auxiliary basis set (CABS) singles correction and the F12b approximation significantly reduce the magnitude of BSSE at the HF and CCSD levels of theory, respectively. For a given cardinal number, we find that non-CP corrected, unscaled triples CCSD(T)-F12b/VXZ-F12 interaction energies are in overall best agreement with the CBS limit.     

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

High-level ab initio versus DFT calculations on (H2O2)2 and H2O2–H2O complexes as prototypes of multiple hydrogen bond systems

The performance of B-LYP, B-P86, B3-LYP, B3-P86, and B3-PW91 density functionals to describe multiple hydrogen bond systems was studied. For this purpose we have chosen the dimers of hydrogen peroxide and the hydrogen peroxide–water complexes. The geometries and vibrational frequencies obtained with a 6-311+G(d,p) basis set were compared with those obtained at the MP2 level using the same basis set expansion. The corresponding dimerization energies were obtained using a 6-311+G(3df,2p) basis set and compared with those obtained using the G2(MP2) theory. Red shiftings of the OH donor stretching frequencies were predicted by all approaches investigated; however, in all cases, the DFT values were sizably larger than the MP2 ones. Similarly, the blue shifting of the torsion of the hydrogen peroxide subunit was larger when evaluated at the DFT level. All functionals reproduced the G2(MP2) relative stabilities of the different local minima quite well. With the exception of the B-LYP and B3-PW91 approaches, all functionals yielded binding energies which deviated from the G2(MP2) values by less than 0.5 kcal/mol, provided that G2-type basis sets were used and that the corresponding BSSE corrections were included. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1124–1135     

Intermolecular potentials of the silane dimer calculated with Hartree-Fock theory, Møller-Plesset perturbation theory, and density functional theory

We have calculated the intermolecular interaction potentials of the silane dimer at the D3d conformation using the Hartree-Fock (HF) self-consistent theory, the correlation-corrected second-order M?ller-Plesset (MP2) perturbation theory, and the density functional theory (DFT) with 108 functionals chosen from the combinations of 9 exchange and 12 correlation functionals. Single-point coupled cluster [CCSD(T)] calculations have also been carried out to calibrate the correlation effect. The HF calculations yield unbound potentials largely because of the exchange-repulsion interaction. In the MP2 calculations, the basis set effects on the repulsion exponent, the equilibrium bond length, the binding energy, and the asymptotic behavior of the calculated intermolecular potentials have been thoroughly studied. We have employed basis sets from the Slater type orbitals fitted with Gaussian functions (STO-nG, n = 3 approximately 6), Pople's medium size basis sets [up to 6-311++G(3df,3pd)], to Dunning's correlation consistent basis sets (cc-pVXZ and aug-cc-pVXZ, X = D, T, Q). With increasing basis size, the repulsion exponent and the equilibrium bond length converge at the 6-31G** basis set and the 6-311++G(3d,3p) basis set, respectively, while a large basis set (aug-cc-pVTZ) is required to converge the binding energy at a chemical accuracy ( approximately 0.05 kcal/mol). Up to the largest basis set used, the asymptotic dispersion coefficient has not converged to the expected C6 value from molecular polarizability calculations. We attribute the slow convergence partly to the inefficacy of using the MP2 calculations with Gaussian type functions to model the asymptotic behavior. Both the basis set superposition error (BSSE) corrected and uncorrected results are presented to emphasize the importance of including such corrections. Only the BSSE corrected results systematically converge to the expected potential curve with increasing basis size. The DFT calculations generate a wide range of interaction patterns, from purely unbound to strongly bound, underestimating or overestimating the binding energy. The binding energies calculated using the OPTXHCTH147, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals and the equilibrium bond lengths calculated using the MPWHCTH93, TPSSHCTH, PBEVP86, PBEP86, PW91TPSS, PW91PBE, and PW91PW91 functionals are close to the MP2 results using the 6-311++G(3df,3pd) basis set. A correlation between the calculated DFT potentials and the exchange and correlation enhancement factors at the low-density region has been elucidated. The asymptotic behaviors of the DFT potentials are also analyzed.