Effect Of Basis Set Superposition Error Research Articles

Systematic Evaluation of Counterpoise Correction in Density Functional Theory.

A widespread belief persists that the Boys-Bernardi function counterpoise (CP) procedure "overcorrects" supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ∼150 atoms and ligand-protein models with ∼300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. This is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP "overcorrects" for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

Behavior of counterpoise correction in many-body molecular clusters of organic compounds: Hartree-Fock interaction energy perspective.

The counterpoise (CP) correction by Boys and Bernardi has been well accepted as a reliable strategy to account for basis set superposition error (BSSE) in intermolecular complexes. The behavior of the CP correction was thoroughly studied in individual molecules of molecular complexes. This work studies the performance of the CP correction in many‐body clusters including three‐body clusters of organic compounds in the 3B‐69 dataset. Additionally, we used crystal structures of polymorphs of benzene, aspirin, and oxalyl dihydrazide (ODH) to construct a many‐body cluster dataset, abbreviated as the MBC‐36 dataset, consisting of two, four and eight molecules, and 16 molecules in the case of benzene. A series of Dunning's basis sets—cc‐pXZ and aug‐cc‐pXZ (X = D and T)—were used to predict CP‐corrected Hartree–Fock (HF) interaction energies of the 3B‐69 and MBC‐36 datasets. The CP‐corrected interaction energies were found to be basis‐set independent, whereas the non‐CP corrected interaction energies were found not to a follow a smooth exponential fitting as previously found for electronic energies of individual molecules. This observation was attributed to the presence of non‐additive induction forces in some clusters. Two 2 × 2 × 2 supercells of benzene polymorphs were constructed to explore the local nature of BSSE effects. A cut‐off radius of 10 Å was demonstrated to be sufficient to fully recover these effects. Although the behavior of CP correction was found to be non‐conventional in many‐body clusters of organic compounds, the use of a small basis set such as cc‐pVDZ showed excellent performance in the prediction of HF interaction energies.

How Do DFT-DCP, DFT-NL, and DFT-D3 Compare for the Description of London-Dispersion Effects in Conformers and General Thermochemistry?

The dispersion-core-potential corrected B3LYP-DCP method (Torres and DiLabio J. Phys. Chem. Lett. 2012, 3, 1738) is for the first time thoroughly assessed and compared with the B3LYP-NL (Hujo and Grimme J. Chem. Theory Comput. 2011, 7, 3866) and B3LYP-D3 (Grimme et al. J. Comput. Chem. 2011, 32, 1456) methods for a broad range of chemical problems that particularly shed light on intramolecular London-dispersion effects in conformers and general thermochemistry. The analysis is based on a compilation of 473 reference cases, the majority of which are taken from the GMTKN30 database (Goerigk and Grimme J. Chem. Theory Comput. 2010, 6, 107; 2011, 7, 291). The results confirm previous findings that B3LYP-DCP indeed predicts very good binding energies for noncovalently bound complexes, particularly with small basis sets. However, problems are identified for the description of intramolecular effects in some conformers and chemical reactions, for which B3LYP-DCP sometimes gives results similar or worse than uncorrected B3LYP. Surprisingly large errors for total atomization energies reveal an unwanted influence of the DCPs on the short-range electronic structure of the investigated systems. However, a recently modified carbon potential for B3LYP-DCP (DiLabio et al. Theor. Chem. Acc. 2013, 132, 1389) was additionally tested that seems to solve most of those problems and provides improved results. An overall comparison between all tested methods shows that B3LYP-NL is the most robust and accurate approach, closely followed by B3LYP-D3. This is also true when small basis sets of double-ζ quality are applied for which those methods have not been parametrized. However, binding energies of noncovalently bound complexes can be more strongly influenced by basis-set superposition-error effects than for B3LYP-DCP. Finally, it is noted that the DFT-D3 and DFT-NL schemes are readily applicable to a large range of chemical elements and they are therefore particularly recommended for more general applications.

Efficient Methods for the Quantum Chemical Treatment of Protein Structures: The Effects of London-Dispersion and Basis-Set Incompleteness on Peptide and Water-Cluster Geometries.

We demonstrate how quantum chemical Hartree-Fock (HF) or density functional theory (DFT) optimizations with small basis sets of peptide and water cluster structures are decisively improved if London-dispersion effects, the basis-set-superposition error (BSSE), and other basis-set incompleteness errors are addressed. We concentrate on three empirical corrections to these problems advanced by Grimme and co-workers that lead to computational strategies that are both accurate and efficient. Our analysis encompasses a reoptimized version of Hobza's P26 set of tripeptide structures, a new test set of conformers of cysteine dimers, and isomers of the water hexamer. These systems reflect features commonly found in protein crystal structures. In all cases, we recommend Grimme's DFT-D3 correction for London-dispersion. We recommend usage of large basis sets such as cc-pVTZ whenever possible to reduce any BSSE effects and, if this is not possible, to use Grimme's gCP correction to account for BSSE when small basis sets are used. We demonstrate that S-S and C-S bond lengths are very prone to basis-set incompleteness and that polarization functions should always be used on S atoms. At the double-ζ level, the PW6B95-D3-gCP DFT method combined with the SVP and 6-31G* basis sets yields accurate results. Alternatively, the HF-D3-gCP/SV method is recommended, with inclusion of polarization functions for S atoms only. Minimal basis sets offer an intriguing route to highly efficient calculations, but due to significant basis-set incompleteness effects, calculated bond lengths are seriously overestimated, making applications to large proteins very difficult, but we show that Grimme's newest HF-3c correction overcomes this problem and so makes this computational strategy very attractive. Our results provide a useful guideline for future applications to the optimization, quantum refinement, and dynamics of large proteins.

Spectroscopic properties of Arx–Zn and Arx–Ag+ (x = 1,2) van der Waals complexes

Potential energy curves have been constructed using coupled cluster with singles, doubles, and perturbative triple excitations (CCSD(T)) in combination with all-electron and pseudopotential-based multiply augmented correlation consistent basis sets [m-aug-cc-pV(n + d)Z; m = singly, doubly, triply, n = D,T,Q,5]. The effect of basis set superposition error on the spectroscopic properties of Ar-Zn, Ar2-Zn, Ar-Ag(+), and Ar2-Ag(+) van der Waals complexes was examined. The diffuse functions of the doubly and triply augmented basis sets have been constructed using the even-tempered expansion. The a posteriori counterpoise scheme of Boys and Bernardi and its generalized variant by Valiron and Mayer has been utilized to correct for basis set superposition error (BSSE) in the calculated spectroscopic properties for diatomic and triatomic species. It is found that even at the extrapolated complete basis set limit for the energetic properties, the pseudopotential-based calculations still suffer from significant BSSE effects unlike the all-electron basis sets. This indicates that the quality of the approximations used in the design of pseudopotentials could have major impact on a seemingly valence-exclusive effect like BSSE. We confirm the experimentally determined equilibrium internuclear distance (re), binding energy (De), harmonic vibrational frequency (ωe), and C(1)Π ← X(1)Σ transition energy for ArZn and also predict the spectroscopic properties for the low-lying excited states of linear Ar2-Zn (X(1)Σg, (3)Πg, (1)Πg), Ar-Ag(+) (X(1)Σ, (3)Σ, (3)Π, (3)Δ, (1)Σ, (1)Π, (1)Δ), and Ar2-Ag(+) (X(1)Σg, (3)Σg, (3)Πg, (3)Δg, (1)Σg, (1)Πg, (1)Δg) complexes, using the CCSD(T) and MR-CISD + Q methods, to aid in their experimental characterizations.

An investigation of the BSSE effect on the evaluation of Ab Initio interaction energies for cisplatin-water complexes

In this article the role played by the effective core potential (Pt atom) on the evaluation of basis set superposition error (BSSE) and interaction energies for the complex between cisplatin ([Pt(NH3)2Cl2]) and water was investigated. Eight distinct spatial monomer orientations in the cisplatin-water complex were considered with the counterpoise procedure being employed to calculate the BSSE. The BSSE effect on potential energy curves, required for the determination of Lennard-Jones parameters, is also discussed. The calculations were performed at the HF, DFT (BLYP, B3LYP, BP86, PBE, PW91 functional) and MP2 level with Pople's split valence and Dunning's correlated consistent basis sets.

Hydrogen bonding in the urea dimers and adenine–thymine DNA base pair: anharmonic effects in the intermolecular H-bond and intramolecular H-stretching vibrations

The equilibrium structures, binding energies, vibrational harmonic frequencies, and the anharmonic corrections for two different (cyclic and asymmetric) urea dimers and for the adenine–thymine DNA base pair system have been studied using the second-order Moller–Plesset perturbation theory (MP2) method and different density functional theory (DFT) exchange–correlation (XC) functionals (BLYP, B3LYP, PBE, HCTH407, KMLYP, and BH and HLYP) with the D95V, D95V**, and D95V++** basis sets. The widely used a posteriori Boys–Bernardi or counterpoise correction scheme for basis set superposition error (BSSE) has been included in the calculations to take into account the BSSE effects during geometry optimization (on structure), on binding energies and on the different levels of approximation used for calculating the vibrational frequencies. The results obtained with the ab initio MP2 method are compared with those calculated with different DFT XC functionals; and finally the suitability of these DFT XC functionals to describe intermolecular hydrogen bonds as well as harmonic frequencies and the anharmonic corrections is assessed and discussed.

Intramolecular Basis Set Superposition Error Effects on the Planarity of DNA and RNA Nucleobases.

Molecules of utmost importance like DNA and RNA nucleobases are predicted to be nonplanar by a typical ab initio method, such as second order Møller-Plesset perturbation theory (MP2) combined with standard Pople's basis sets. Similarly to the case of other planar aromatic systems, these pitfalls can be explained in terms of intramolecular basis set superposition error (BSSE) effects, induced by local basis set deficiencies. We demonstrate that conventional BSSE correction techniques such as the Counterpoise method can account for this wrong behavior and provide proper correction whenever spurious results occur, mainly in case of thymine, uracil and guanine but also to lower extent for adenine and cytosine. We also show that special care must be taken when assessing the BSSE by means of ghost-orbital calculations for strongly overlapping fragments. Often molecular orbitals in the isolated fragment calculation have a different orientation as in the ghost-orbital calculation. This can lead to bogus derivatives of the CP-correction term, essential to account for geometry and vibrational BSSE effects.

Quantitative assessment of the effect of basis set superposition error on the electron density of molecular complexes by means of quantum molecular similarity measures

AbstractThe Chemical Hamiltonian Approach (CHA) method is applied to obtain Basis Set Superposition Error (BSSE)‐free molecular orbitals at the Hartree–Fock (HF) and Density Functional Theory (DFT) levels of theory. To assess qualitatively the effect of the BSSE on the first‐order electron density, we had previously applied Bader's analysis of the intermolecular critical points located on the electron density, as well as density difference maps for several hydrogen bonded complexes. In this work, Quantum Molecular Similarity Measures are probed as an alternative avenue to properly quantify the electronic relaxation due to the BSSE removal by means of distance indices between the uncorrected and corrected charge densities. It is shown that BSSE contamination is more important at the DFT level of theory, and in some cases, changes on the topology of the electron density are observed upon BSSE correction. Inclusion of diffuse functions have been found to dramatically decrease the BSSE effect in both geometry and electron density. The CHA method represents a good compromise to obtain accurate results with small basis sets. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009

Intramolecular basis set superposition error effects on the planarity of benzene and other aromatic molecules: A solution to the problem

Recently, the surprising result that ab initio calculations on benzene and other planar arenes at correlated MP2, MP3, configuration interaction with singles and doubles (CISD), and coupled cluster with singles and doubles levels of theory using standard Pople's basis sets yield nonplanar minima has been reported. The planar optimized structures turn out to be transition states presenting one or more large imaginary frequencies, whereas single-determinant-based methods lead to the expected planar minima and no imaginary frequencies. It has been suggested that such anomalous behavior can be originated by two-electron basis set incompleteness error. In this work, we show that the reported pitfalls can be interpreted in terms of intramolecular basis set superposition error (BSSE) effects, mostly between the C-H moieties constituting the arenes. We have carried out counterpoise-corrected optimizations and frequency calculations at the Hartree-Fock, B3LYP, MP2, and CISD levels of theory with several basis sets for a number of arenes. In all cases, correcting for intramolecular BSSE fixes the anomalous behavior of the correlated methods, whereas no significant differences are observed in the single-determinant case. Consequently, all systems studied are planar at all levels of theory. The effect of different intramolecular fragment definitions and the particular case of charged species, namely, cyclopentadienyl and indenyl anions, respectively, are also discussed.

The structure and binding energies of the van der Waals complexes of Ar and N2 with phenol and its cation, studied by high level ab initio and density functional theory calculations

We have investigated, using both ab initio and density functional theory methods, the minimum energy structures and corresponding binding energies of the van der Waals complexes between phenol and argon or the nitrogen molecule, and the corresponding complexes involving the phenol cation. Structures were obtained at the MP2 level using a large basis, and the corresponding energies were corrected for basis set superposition error (BSSE), higher order electron correlation effects, and for basis set size. The structures of the global minima were further refined for the effects of BSSE and the corresponding binding energies were evaluated. For each neutral species, we find only a single true minimum, pi bonded for argon and OH bonded for nitrogen. For both cationic species, we find that the OH-bonded complex is preferred over other minima which we have identified as having Ar or N(2) between exogeneous atoms. The ab initio calculations are generally in excellent agreement with experimental binding energies and rotational constants. We find that the B3LYP functional is particularly poor at describing these complexes, while a density functional theory (DFT) method with an empirical correction for dispersive interactions (DFT-D) is very successful, as are some of the new functionals proposed by Zhao and Truhlar [J. Phys. Chem. A 109, 5656 (2005); J. Chem. Theory Comput. 2, 1009 (2006); Phys. Chem. Chem. Phys. 7, 2701 (2005); J. Phys. Chem. A 108, 6908 (2004)]. Both the ab initio and DFT-D methods accurately predict the intermolecular vibrational modes.

Comparison of ab Initio and DFT Electronic Structure Methods for Peptides Containing an Aromatic Ring: Effect of Dispersion and BSSE

We establish that routine B3LYP and MP2 methods give qualitatively wrong conformations for flexible organic systems containing pi systems and that recently developed methods can overcome the known inadequacies of these methods. This is illustrated for a molecule (a conformer of the Tyr-Gly dipeptide) for which B3LYP/6-31+G(d) and MP2/6-31+G(d) geometry optimizations yield strikingly different structures [Mol. Phys. 2006, 104, 559-570]: MP2 predicts a folded "closed-book" conformer with the glycine residue located above the tyrosine ring, whereas B3LYP predicts a more open conformation. By employing different levels of theory, including the local electron correlation methods LMP2 (local MP2) and LCCSD(T0) (local coupled cluster with single, double, and noniterative local triple excitations) and large basis sets (aug-cc-pVnZ, n=D, T, Q), it is shown that the folded MP2 minimum is an artifact caused by large intramolecular BSSE (basis set superposition error) effects in the MP2/6-31+G(d) calculations. The B3LYP functional gives the correct minimum, but the potential energy apparently rises too steeply when the glycine and tyrosine residues approach each other, presumably due to missing dispersion effects in the B3LYP calculations. The PWB6K and M05-2X functionals, designed to give good results for weak interactions, remedy this to some extent. The reduced BSSE in the LMP2 calculations leads to faster convergence with increasing basis set quality, and accurate results can be obtained with smaller basis sets as compared to canonical MP2. We propose LMP2 as a suitable method to study interactions with pi-electron clouds.

Is Spin-Component Scaled Second-Order Møller−Plesset Perturbation Theory an Appropriate Method for the Study of Noncovalent Interactions in Molecules?

Testing of the spin-component scaled second-order Møller-Plesset (SCS-MP2) method for the computation of noncovalent interaction energies is done with a database of 165 biologically relevant complexes. The effects of the spin-scaling procedure (i.e., MP2 vs SCS-MP2), the basis set size, and the corrections for basis set superposition error (BSSE) are systematically examined. When using two-point basis set extrapolations for the correlation energy, augmentation of the atomic orbital basis with computationally costly diffuse functions is found to be obsolete. In general, SCS-MP2 also improves results for noncovalent interactions statistically on MP2, and significant outliers are removed. Moreover, it is shown that effects of BSSE and one-particle basis set incompleteness almost cancel each other in the case of triple-zeta sets (SCS-MP2/TZVPP or SCS-MP2/cc-pVTZ without counterpoise correction), which opens a practical route to efficient computations for large systems. We recommend SCS-MP2 as the preferred quantum chemical wave function based method for the noncovalent interactions in large biologically relevant systems when reasonable coupled-cluster with single and double and perturbative triple excitations (CCSD(T)) calculations cannot be performed anymore. A comparison to MP2 and CCSD(T) interaction energies for n-alkane dimers, however, indicates (and this also holds to a lesser extent for hydrogen-bonded systems) limitations of SCS-MP2 when treating chemically "saturated" interactions. The different behavior of second-order perturbation theory for saturated and for stacked pi-systems is discussed.

Blue‐shifting hydrogen bond in the benzene–benzene and benzene–naphthalene complexes

Ab initio complete optimizations at MP2/6-31++G** level have been performed in the T-shaped geometry of the benzene-benzene and benzene-naphthalene complexes. To check the effect of the basis set superposition error (BSSE), optimizations have been done in the BSSE corrected and BSSE uncorrected potential energy surfaces. The BSSE effect in the calculation of the Hessian has also been evaluated to check its influence in the frequency values. Quantum theory atoms in molecules (QTAIM) calculations have also been performed on both dimers. Intermolecular energies differ around a 25% when the optimization is performed with or without counterpoise corrected gradients. The influence of BSSE is also noticeable in the distances. Frequency shifts show big changes because of the BSSE. Thus, uncorrected values are up 350% larger than corrected ones. The hypotheses given in the literature to explain the origin of the blue-shifting hydrogen bond do not seem to give a suitable explanation for all characteristics of the behavior found in the studied systems.

Calculation of smooth potential energy surfaces using local electron correlation methods

The geometry dependence of excitation domains in local correlation methods can lead to noncontinuous potential energy surfaces. We propose a simple domain merging procedure which eliminates this problem in many situations. The method is applied to heterolytic bond dissociations of ketene and propadienone, to SN2 reactions of Cl(-) with alkylchlorides, and in a quantum mechanical/molecular mechanical study of the chorismate mutase enzyme. It is demonstrated that smooth potentials are obtained in all cases. Furthermore, basis set superposition error effects are reduced in local calculations, and it is found that this leads to better basis set convergence when computing barrier heights or weak interactions. When the electronic structure strongly changes between reactants or products and the transition state, the domain merging procedure leads to a balanced description of all structures and accurate barrier heights.

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.

Test of a nonempirical density functional: Short-range part of the van der Waals interaction in rare-gas dimers

It is known that the nonempirical generalized gradient approximation (GGA) of Perdew, Burke, and Ernzerhof (PBE) provides a much more realistic description of the short-range part of the van der Waals (vdW) interaction than does the local spin density (LSD) approximation. In the present work, the ability of the higher-level nonempirical meta-GGA of Tao, Perdew, Staroverov, and Scuseria (TPSS) to describe vdW interaction is tested self-consistently in ten rare-gas dimers with Z< or =36. The one-parameter hybrid version (TPSSh) of the TPSS exchange-correlation functional is also included in this test. Calculations show that both TPSS and TPSSh functionals correctly yield vdW bonds in these dimers and significantly improve the prediction of bond lengths, binding energies, and harmonic vibrational frequencies over LSD. The rather close agreement of TPSS with PBE for these dimers confirms a principle of the TPSS construction: preservation of the PBE large-gradient behavior. More importantly, it suggests that TPSS can serve as a platform on which to construct a still-higher level of nonempirical functionals. Compared with the PBE GGA, TPSS, and TPSSh yield a slightly weaker binding. As for normally bonded molecules, TPSSh yields the most accurate vibrational frequencies. The typically too-long bond lengths and too-small binding energies of TPSS meta-GGA suggest the need for some long-range vdW interaction correction even in this class of systems. The effect of basis-set superposition error on the calculated properties of these vdW systems is investigated. We also show that the relatively strong anharmonic effects in the rare-gas dimers are described remarkably well by the Morse potential.

BSSE‐corrected geometry and harmonic and anharmonic vibrational frequencies of formamide–water and formamide–formamide dimers

AbstractThe basis set superposition error (BSSE) influence in the geometry structure, interaction energies, and intermolecular harmonic and anharmonic vibrational frequencies of cyclic formamide–formamide and formamide–water dimers have been studied using different basis sets (6‐31G, 6‐31G**, 6‐31++G**, D95V, D95V**, and D95V++**). The a posteriori “counterpoise” (CP) correction scheme has been compared with the a priori “chemical Hamiltonian approach” (CHA) both at the Hartree–Fock (HF) and second‐order Møller–Plesset many‐body perturbation (MP2) levels of theory. The effect of BSSE on geometrical parameters, interaction energies, and intermolecular harmonic vibrational frequencies are discussed and compared with the existing experimental data. As expected, the BSSE‐free CP and CHA interaction energies usually show less deep minima than those obtained from the uncorrected methods at both the HF and MP2 levels. Focusing on the correlated level, the amount of BSSE in the intermolecular interaction energies is much larger than that at the HF level, and this effect is also conserved in the values of the force constants and harmonic vibrational frequencies. All these results clearly indicate the importance of the proper BSSE‐free correlation treatment with the well‐defined basis functions. At the same time, the results show a good agreement between the a priori CHA and a posteriori CP correction scheme; this agreement is remarkable in the case of large and well‐balanced basis sets. The anharmonic frequency correction values also show an important BSSE dependence, especially for hydrogen bond stretching and for low frequencies belonging to the intermolecular normal modes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005

Theoretical studies of hydrogen bonding in water–cyanides and in the base pair Gu–Cy

Density-functional (DFT) and many-body-perturbation theories (MBPT/CC) are used to study the hydrogen bonding in the water–cyanide complexes H–CN⋯H 2O, H 3C–CN⋯H 2O and (CH 3) 3C–CN⋯H 2O. Structures, binding energies and changes in vibrational frequencies are analyzed. The calculated CN stretching frequency is found to shift to the blue upon complexation in H–CN⋯H 2O and H 3C–CN⋯H 2O. To investigate electron correlation effects on the binding energies of these complexes, single-point calculations are performed at the MBPT/CC (MP2, MP3, MP4, CCSD and CCSD(T)) levels using the optimized MP2 geometries. Binding energies are also obtained at different levels of DFT (B3LYP and PW91) and compared with the MBPT/CC results. All calculations include corrections for basis set superposition error (BSSE) and zero-point vibrational energies. Additionally, the triple hydrogen-bonded guanine–cytosine (Gu–Cy) base pair is analyzed. The binding energy of the Watson–Crick model for Gu–Cy is calculated using the Hartree–Fock calculations and DFT (B3LYP and BP86) methods. The results for the hydrogen bonding distances and binding energies are in good agreement with experimental and recent theoretical values. The calculated dipole moment of the Gu–Cy complex is compared with the direct vector sum of the isolated bases. After taking into account the BSSE effects we find that the electron polarization due to the hydrogen binding leads to an increase of ∼20% of the calculated dipole moment of the complex.

Infrared spectra of C2H4–HCl complex

We report the first rotationally resolved observation of the infrared spectrum of the molecular complex C2H4–HCl. The complex was produced by a supersonic expansion through a pulsed slit jet. By means of a high-resolution tunable diode laser spectrometer, we have recorded the spectrum of the HCl stretching vibration for the isotopes C2H4–H35Cl and C2H4–H37Cl. From the analysis of the spectra, we determined the rotational constants and vibrational frequencies of both isotopes. These experimental results have been interpreted in view of obtaining information on the intermolecular interaction. The experimental data have been complemented by standard coupled cluster singles and doubles model including connected triple excitations with a correlation consistent polarized valence triple zeta basis set ab initio followed by grid calculations, in order to study the anharmonicity, the coupling between intramolecular and intermolecular motions, and the basis set superposition error effects. The results obtained in this study are compared to those of our previous work on the C2H2–HCl [J. Chem. Phys. 113, 4876 (2000)].