Basis set effects on the intermolecular interaction of the H2-H2 system obtained using ab initio molecular orbital calculations with the Møller-Plesset perturbation correction

The intermolecular interaction potential of the H₂-H₂ system was calculated by an ab initio molecular orbital method using several basis sets (up to 6-31 lG(3pd)) with inclusion of the electron correlation correction of the Møller-Plesset perturbation method and the basis set superposition error (BSSE) correction of the counterpoise method in order to evaluate the basis set effect. The calculated interaction energies depend strongly on the basis set used. Whereas the interaction energies of the repulsive and coulombic energy components calculated at the Hartree-Fock level are not affected by a change of basis set, the dispersion energy component depends strongly on the basis set used. Parameters of an exp-6-1 type non-bonding interaction potential were optimized on the basis of the MP4(SDTQ)/6-311G(3p) level intermolecular interaction energies of the H₂-H₂ system.     

Some comments on the counterpoise correction for the basis set superposition error at the correlated level

The use of the counterpoise method for the mitigation of basis set superposition error at the correlated level is discussed. Evidence is presented to show that the ghost basis plays a dual role in the counterpoise method: The orbitals of the system are improved by the ghost basis but at the expense of a nonphysical increase in the dimension of the virtual space. This second factor has no effect on application of the counterpoise method at the SCF level but it makes the use of the counterpoise method at the correlated level much less straightforward. © 1993 John Wiley & Sons, Inc.     

On basis set superposition error corrected stabilization energies for large n-body clusters

In this contribution, we propose an approximate basis set superposition error (BSSE) correction scheme for the site-site function counterpoise and for the Valiron-Mayer function counterpoise correction of second order to account for the basis set superposition error in clusters with a large number of subunits. The accuracy of the proposed scheme has been investigated for a water cluster series at the CCSD(T), CCSD, MP2, and self-consistent field levels of theory using Dunning's correlation consistent basis sets. The BSSE corrected stabilization energies for a series of water clusters are presented. A study regarding the possible savings with respect to computational resources has been carried out as well as a monitoring of the basis set dependence of the approximate BSSE corrections.     

An interaction potential between an alanine zwitterion and a water molecule based on ab initio calculations

Interaction energies between an alanine zwitterion and a water molecule at 150 different positions and orientations have been calculated using the ab initio method with the minimal basis set and employing the counterpoise method to eliminate the basis set superposition error. Dispersion energies are estimated using the Slater–Kirkwood formula. Out of a total of 150 computed interaction energies, 140 whose SCF interaction energies are below 5 kcal/mol have been fitted with a summation of atom-atom pair potentials in the form of the Lennard–Jones potential plus an electrostatic term. The standard deviation for this fitting is 0.49 kcal/mol. A sampling scheme regarding geometrical configurations is presented. Twenty rays are uniformly drawn from the origin of coordinates, a floatable division with equal ratios is made along each ray, and one of 60 orientations is randomly taken as the orientation of a water molecule. A nonlinear fitting method is used with a restriction on the sign change of fitting coefficients.     

Counterpoise corrected interaction energies are not systematically better than uncorrected ones: comparison with CCSD(T) CBS extrapolated values

The effect of the inclusion of counterpoise corrections (CP) on the accuracy of interaction energies has been studied for different systems accounting for (1) intermolecular interactions, (2) intramolecular interactions and (3) chemical reactions. To minimize the error associated with the method of choice, the energy calculations were performed using CCSD(T) in all the cases. The values obtained using aug-cc-pVXZ basis sets are compared to CBS-extrapolated values. It has been concluded that at least for the tested systems CP corrections systematically leads to results that differ from the CBS-extrapolated ones to a larger extension than the uncorrected ones. Accordingly, from a practical point of view, we do not recommend the inclusion of such corrections in the calculation of interaction energies, except for CBS extrapolations. The best way of dealing with basis set superposition error (BSSE) is not to use CP corrections, but to make a computational effort for increasing the basis set. This approach does not eliminate BSSE but significantly decreases it, and more importantly it proportionally decreases all the errors arising from the basis set truncation.     

Ne-CO2的从头算势能面及微波光谱

采用三重激发校正的耦合簇[CCSD(T)]方法和大基组计算了范德华复合物Ne-CO2的分子间势能面. 分子间相互作用能的计算采用考虑了基组重叠误差修正的超分子方法. 计算结果表明, 该势能面有两个极小值点, 分别对应T形构型和线性Ne-OCO构型. 采用离散变量表象(DVR)方法及Lanczos算法计算了Ne-CO2的振转能级. 计算结果表明, 体系势能面支持22个振动束缚态. 计算得到的微波光谱的跃迁频率与实验值吻合得很好.     

咪唑氟硼酸类离子液体熔点与分子内相互作用能的关系

运用密度泛函理论B3LYP方法及6-311++G(d,p)基组对11种咪唑氟硼酸离子液体进行了研究.选择相应化合物的离子体系{[XIM][BF4]n}(n-1)-(n=2,3)作为研究对象,即研究体系由一个烷基咪唑阳离子XIM+和2-3个BF4-阴离子构成,对其进行结构优化.在优化得到的最低能量构型的基础上计算了分子内阳离子与阴离子间的相互作用能,同时考虑了基组重叠误差的修正.结果表明所研究离子体系的离子间相互作用能与离子液体的实验熔点之间存在明确的线性关系,并且所得到的线性方程与氨基酸阳离子型离子液体中存在的线性关系相近.我们的工作为今后借助量子化学方法设计功能化离子液体提供了一定的理论基础.     

Ab initio intermolecular potential-energy surface and microwave spectra for the Ne-OCS complex

An ab initio potential-energy surface for the Ne-OCS complex was calculated using the coupled-cluster singles and doubles with noniterative inclusion of connected triples [CCSD(T)] with a large basis set containing bond functions. The interaction energies were obtained by the supermolecular approach with the full counterpoise correction for the basis set superposition error. The CCSD(T) potential was found to have three minima corresponding to the T-shaped and the linear Ne-SCO and Ne-OCS structures. The two-dimensional discrete variable representation method was employed to calculate the rovibrational energy levels for five isotopomers Ne-OCS, (22)Ne-OCS, Ne-OC(34)S, Ne-O(13)CS, and Ne-(18)OCS. The calculated pure rotational transition frequencies for the vibrational ground state of the five isotopomers are in good agreement with the observed values. The corresponding microwave spectra show that the b-type transitions (deltaK(a)=+/-1) are significantly stronger than the a-type transitions (deltaK(a)=0).     

Basis set superposition effects on properties of interacting systems. Dipole moments and polarizabilities

Basis set superposition effects which are not removed by the counterpoise correction are shown to modify the electric properties of interacting subsystems and influence indirectly the calculated interaction energies. The role of these higher-order basis set superposition effects is illustrated by the calculation of the water molecule dipole moment and polarizability in the basis set of the water dimer.     

Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods

Potential energy curves for the parallel-displaced, T-shaped and sandwich structures of the benzene dimer are computed with density fitted local second-order M?ller-Plesset perturbation theory (DF-LMP2) as well as with the spin-component scaled (SCS) variant of DF-LMP2. While DF-LMP2 strongly overestimates the dispersion interaction, in common with canonical MP2, the DF-SCS-LMP2 interaction energies are in excellent agreement with the best available literature values along the entire potential energy curves. The DF-SCS-LMP2 dissociation energies for the three structures are also compared with new complete basis set estimates of the interaction energies obtained from accurate coupled cluster (CCSD(T)) and DF-SCS-MP2 calculations. Since LMP2 is essentially free of basis set superposition errors, counterpoise corrections are not required. As a result, DF-SCS-LMP2 is computationally inexpensive and represents an attractive method for the study of larger pi-stacked systems such as truncated sections of DNA.     

Quantum mechanical study of the nonbonded forces in water-methanol complexes

The water-methanol dimer can adopt two possible configurations (WdM or MdW) depending on whether the water or the methanol acts as the hydrogen bond donor. The relative stability between the two configurations is less than 1 kcal/mol, and as a result, this dimer has been a challenging system to investigate using either theoretical or experimental techniques. In this paper, we present a systematic study of the dependence of the geometries, interaction energies, and harmonic frequencies on basis sets and treatment of electron correlation for the two configurations. At the highest theory level, MP2/aug-cc-pVQZ//MP2/aug-cc-pVTZ, interaction energies of -5.72 and -4.95 kcal/mol were determined for the WdM and MdW configurations, respectively, after correcting for basis set superposition error using the Boys-Bernardi counterpoise scheme. Extrapolating to the complete basis set limit resulted in interaction energies of -5.87 for WdM and -5.16 kcal/mol for MdW. The energy difference between the two configurations is larger than the majority of previously reported values, confirming that the WdM complex is preferred, in agreement with experimental observations. The effects that electron correlation have on the geometry were investigated by performing optimization at the MP2(full), MP4, and CCSD levels of theory. The approach trajectories for the formation of each dimer configuration are presented and the importance of these trajectories in the development of parameters for use in classical force fields is discussed.     

Conformational and energetic properties of the ammonia dimer—comparison of post-Hartree—Fock and density functional methods

The equilibrium structure of the ammonia dimer has been investigated with density functional and MP2 calculations. We used Slater- and Becke-exchange functionals combined with correlation functionals as recommended by Vosko-Wilk-Nusair, by Perdew, and by Lee-Yang-Parr, respectively. The potential energy surfaces was investigated. The asymmetric cyclic “microwave” structure could be identified as a minimum. Optimization of the intermolecular parameters showed that this structure has nearly the same energy as the centrosymmetric cyclic structure. Full optimization transformed the asymmetric cyclic structure into the linear structure. The interaction energies in the dimer were corrected for the basis set superposition error using the Boys-Bernardi counterpoise method and the a priori chemical Hamiltonian approach, respectively. © 1996 by John Wiley & Sons, Inc.     

Efficiency of numerical basis sets for predicting the binding energies of hydrogen bonded complexes: evidence of small basis set superposition error compared to Gaussian basis sets

Binding energies of selected hydrogen bonded complexes have been calculated within the framework of density functional theory (DFT) method to discuss the efficiency of numerical basis sets implemented in the DFT code DMol3 in comparison with Gaussian basis sets. The corrections of basis set superposition error (BSSE) are evaluated by means of counterpoise method. Two kinds of different numerical basis sets in size are examined; the size of the one is comparable to Gaussian double zeta plus polarization function basis set (DNP), and that of the other is comparable to triple zeta plus double polarization functions basis set (TNDP). We have confirmed that the magnitudes of BSSE in these numerical basis sets are comparative to or smaller than those in Gaussian basis sets whose sizes are much larger than the corresponding numerical basis sets; the BSSE corrections in DNP are less than those in the Gaussian 6-311+G(3df,2pd) basis set, and those in TNDP are comparable to those in the substantially large scale Gaussian basis set aug-cc-pVTZ. The differences in counterpoise corrected binding energies between calculated using DNP and calculated using aug-cc-pVTZ are less than 9 kJ/mol for all of the complexes studied in the present work. The present results have shown that the cost effectiveness in the numerical basis sets in DMol3 is superior to that in Gaussian basis sets in terms of accuracy per computational cost.     

Møller-Plesset perturbation energies and distances for HeC(20) extrapolated to the complete basis set limit

The potential energy surface for the C₂₀–He interaction is extrapolated for three representative cuts to the complete basis set limit using second‐order Møller–Plesset perturbation calculations with correlation consistent basis sets up to the doubly augmented variety. The results both with and without counterpoise correction show consistency with each other, supporting that extrapolation without such a correction provides a reliable scheme to elude the basis‐set‐superposition error. Converged attributes are obtained for the C₂₀– He interaction, which are used to predict the fullerene dimer ones. Time requirements show that the method can be drastically more economical than the counterpoise procedure and even competitive with Kohn‐Sham density functional theory for the title system. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Basis set dependence of solute-solvent interaction energy of benzene in water: a HF/DFT study

Solute-solvent interaction energies for the benzene molecule dissolved in water are computed using Hartree-Fock and B3LYP density functional theories. Explicit solvent molecules up to 14-A away from the dissolved benzene molecule are included in the calculation of interaction energies. Both basis set dependence and basis set superposition errors are carefully examined. It is found that the use of a larger basis set for the region near the solute together with a smaller basis set for the outer region gives results very close to what would have been obtained if the larger basis set had been used for the whole system. It is also shown that a correction for the basis set superposition error is a necessary component in this kind of calculations. With this correction, results obtained with different tested basis sets converge to within 1 kcal/mol.     

Theoretical study of the prion protein based on the fragment molecular orbital method

We performed fragment molecular orbital (FMO) calculations to examine the molecular interactions between the prion protein (PrP) and GN8, which is a potential curative agent for prion diseases. This study has the following novel aspects: we introduced the counterpoise method into the FMO scheme to eliminate the basis set superposition error and examined the influence of geometrical fluctuation on the interaction energies, thereby enabling rigorous analysis of the molecular interaction between PrP and GN8. This analysis could provide information on key amino acid residues of PrP as well as key units of GN8 involved in the molecular interaction between the two molecules. The present FMO calculations were performed using an original program developed in our laboratory, called “Parallelized ab initio calculation system based on FMO (PAICS)”. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

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

The 5-fluorocytosine (5-FC) is a fluorinated cytosine analog that is used as an antifungal agent. In this work, we present the hydrogen-bonding base pairs involving 5-FC 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, 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 decreases in the following order: 5-FC:G > 5-FC:C > 5-FC:A > 5-FC:T.     

Monomer basis-set truncation effects in calculations of interaction energies: A model study

Supermolecular interaction energies are analyzed in terms of the symmetry-adapted perturbation theory and operators defining the inaccuracy of the monomer wave functions. The basis set truncation effects are shown to be of first order in the monomer inaccuracy operators. On the contrary, the usual counterpoise correction schemes are of second order in these operators. Recognition of this difference is used to suggest an approach to corrections for basis-set truncation effects. Also earlier claims--that dimer-centered basis sets may lead to interaction energies free of basis-set superposition effects--are shown to be misleading. According to the present study the basis-set truncation contributions, evaluated by means of the symmetry-adapted perturbation theory with monomer-centered basis sets, provide physically and mathematically justified corrections to supermolecular results for interaction energies.     

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.