New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.     

Heats of formation for Fe(CO)n(n=1–4)

The dissociation energies of Fe(CO)_n (n=2–4) are computed using correlation consistent basis sets and the CCSD(T) approach. The dissociation energies are extrapolated to the CBS limit and are corrected for core–valence (CV), scalar relativistic, spin–orbit, zero-point, and thermal effects. Our iron carbonyl bond strengths agree with experiment within the respective error bars. We use our dissociations energies at 298 K to obtain the heats of formation of Fe(CO)_n (n=1–4).     

The effect of basis set superposition error on the convergence of interaction energies

 For the intermolecular interaction energies of ion-water clusters [OH⁻(H₂O) _n (n=1,2), F⁻(H₂O), Cl⁻(H₂O), H₃O⁺(H₂O) _n (n=1,2), and NH₄ ⁺(H₂O) _n (n=1,2)] calculated with correlation-consistent basis sets at MP2, MP4, QCISD(T), and CCSD(T) levels, the basis set superposition error is nearly zero in the complete basis set (CBS) limit. That is, the counterpoise-uncorrected intermolecular interaction energies are nearly equal to the counterpoise-corrected intermolecular interaction energies in the CBS limit. When the basis set is smaller, the counterpoise-uncorrected intermolecular interaction energies are more reliable than the counterpoise-corrected intermolecular interaction energies. The counterpoise-uncorrected intermolecular interaction energies evaluated using the MP2/aug-cc-pVDZ level is reliable. Received: 14 March 2001 / Accepted: 25 April 2001 / Published online: 9 August 2001     

On the basis set superposition error in supermolecular calculations of interaction-induced electric properties: many-body components

In the present paper we analyze basis set superposition error (BSSE) removal methods from many-body components of interaction-induced electric properties. The Valiron–Mayer function counterpoise (VMFC), site–site function counterpoise (SSFC) and TB methods have been employed in order to obtain the incremental optical components of linear hydrogen fluoride clusters (HF)n, where n = {3,4}. Following Mierzwicki and Latajka, who have performed similar calculations for the interaction energy, we compare those three methods of eliminating BSSE using several Dunning’s correlation consistent basis sets.     

Basis sets for the evaluation of van der Waals complex interaction energies: Ne–N2 intermolecular potential and microwave spectrum

In order to obtain efficient basis sets for the evaluation of van der Waals complex intermolecular potentials, we carry out systematic basis set studies. For this, interaction energies at representative geometries on the potential energy surfaces are evaluated using the CCSD(T) correlation method and large polarized LPol‐n and augmented polarization‐consistent aug‐pc‐2 basis sets extended with different sets of midbond functions. On the basis of the root mean square errors calculated with respect to the values for the most accurate potentials available, basis sets are selected for fitting the corresponding interaction energies and getting analytical potentials. In this work, we study the Ne–N₂ van der Waals complex and after the above procedure, the aug‐pc‐2–3321 and the LPol‐ds‐33221 basis set results are fitted. The obtained potentials are characterized by T‐shaped global minima at distances between the Ne atom and the N₂ center of mass of 3.39 Å, with interaction energies of ?49.36 cm^?1 for the aug‐pc‐2–3321 surface and ?50.28 cm^?1 for the LPol‐ds‐33221 surface. Both sets of results are in excellent agreement with the reference surface. To check the potentials further microwave transition frequencies are calculated that agree well with the experimental and the aV5Z‐33221 values. The success of this study suggests that it is feasible to carry out similar accurate calculations of interaction energies and ro‐vibrational spectra at reduced cost for larger complexes than has been possible hitherto. © 2013 Wiley Periodicals, Inc.     

Binding energies of hydrogen-bonded clusters from extrapolation-oriented basis sets

The MP2 and CCSD(T) basis set limit binding energies of various hydrogen-bonded clusters were estimated via basis set extrapolation employing the correlation consistent aug-cc-pVDZ and modified aug-cc-pVDZ set containing extra polarization functions from cc-pVTZ set. By adopting the optimal interval for the difference between the cardinal numbers (X) corresponding to two basis sets in the X ⁻³ type extrapolation scheme the estimated binding energies for (H₂O)_n and (HF)_n (n=3−5) are shown to be close to the reference basis set limit values within the error bounds in many cases, manifesting the significance of these basis sets in studying the structures and binding of large hydrogen-bonded clusters.     

Theoretical insight into structural and electronic properties of cationic Scn+ (n=2-13): A benchmark study

Geometric, thermodynamic and electronic properties of cationic scandium clusters are studied. Geometric optimizations and stable spin states of Sc₂⁺ are assessed on high level ab-initio coupled cluster method CCSD(T) with different dunning correlation consistent basis sets (aug-cc-pVDZ, aug-cc-pVTZ and aug-cc-pVQZ). Then, 23 DFT functionals belonging to different classes are evaluated at 6-31G (d), LANL2MB, LANL2DZ and Def2-SVP basis sets, and the results are compared with the benchmarked coupled cluster calculations. Due to excellent correlation, PBEPBE/LANL2DZ was chosen to perform calculation of higher scandium cationic clusters Sc_n⁺ (n = 3-13). In addition, we explored relative stability, binding energies, second order energy differences, vertical ionization energies, vertical electron affinities and HOMO-LUMO gaps. Moreover, these results are also compared with the neutral scandium clusters.     

Theoretical reference values for the AE6 and BH6 test sets from explicitly correlated coupled-cluster theory

We propose a new computational protocol to obtain highly accurate theoretical reference data. This protocol employs the explicitly correlated coupled-cluster method with iterative single and double excitations as well as perturbative triple excitations, CCSD(T)(F12), using quadruple-z\zeta basis sets. Higher excitations are accounted for by conventional CCSDT(Q) calculations using double-z\zeta basis sets, while core/core-valence correlation effects are estimated by conventional CCSD(T) calculations using quadruple-z\zeta basis sets. Finally, scalar-relativistic effects are accounted for by conventional CCSD(T) calculations using triple-z\zeta basis sets. In the present article, this protocol is applied to the popular test sets AE6 and BH6. An error analysis shows that the new reference values obtained by our computational protocol have an uncertainty of less than 1 kcal/mol (chemical accuracy). Furthermore, concerning the atomization energies, a cancellation of the basis set incompleteness error in the CCSD(T)(F12) perturbative triples contribution with the corresponding error in the contribution from higher excitations is observed. This error cancellation is diminished by the CCSD(T*)(F12) method. Thus, we recommend the use of the CCSD(T*)(F12) method only for small- and medium-sized basis sets, while the CCSD(T)(F12) approach is preferred for high-accuracy calculations in large basis sets.     

Structure, Electronic Properties and Interaction of MRn n + (n?=?1–3, M?=?Cu, Ag and Au) Clusters: Ab Initio Calculations

The structures and stabilities of MRn _n ⁺ (n?=?1–3, M?=?Cu, Ag and Au) series at the CCSD(T) theoretical level are performed. The n?=?2 systems are more stable than its neighbours. The role of the interaction is investigated using the natural bond orbital analysis, Laplacian, electron localization function and reduced density gradient analysis. The results show the intermediate character in the M–Rn interaction.     

Polarized basis sets for accurate calculations of static and dynamic electric properties of molecules

We report on the development and testing of large polarized basis sets (LPolX, where X is the element symbol) for accurate calculations of linear and nonlinear electric properties of molecules. The method used to generate LPolX sets is based on our studies of the analytic dependence of Gaussian functions on external time‐independent and time‐dependent electric fields. At variance with the earlier investigations of small, highly compact (ZPolX) basis sets for moderately accurate calculations of electric properties of large molecules, the present goal is to obtain basis sets that are nearly saturated with respect to the selected class of electric properties and can be used for accurate studies of interaction‐induced properties. This saturation makes the LPolX sets also useful in calculations of optical properties for chiral molecules. In this article, the LPolX sets are generated for X = H, C, N, O, and F, and examined in calculations of linear and nonlinear electric properties of four standard test systems: HF, N₂, CO, and HCN. The study of the performance of LPolX basis sets has been carried out at different levels of approximation ranging from the SCF HF method to highly correlated CCSD(T) approach. The results obtained in this study compare favorably with accurate reference data and show a high level of saturation of LPolX basis sets with respect to the polarization effect due to external electric fields. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

CCSDT and CCSD(T) calculations on model H-bonded and stacked complexes

The CCSD(T) and CCSDT interaction energies were determined for model planar H-bonded complexes (formamide…formamide, formamidine…formamidine) and stacked complexes (ethylene…ethylene, formaldehyde…formaldehyde). Various basis sets from the 6-31G*(0.25) to aug-cc-pVDZ were used. Difference between CCSD(T) and CCSDT interaction energies were small and become negligible (bellow 0.1 kcal/mol) if the aug-cc-pVDZ (or aug-cc-pVDZ/cc-pVDZ) basis set was applied. This result strongly supports the use of the CCSD(T) method for determination of true stabilization energies of extended complexes.     

Chemical bond between Cu(II) and Rn: ab initio study of CuRn n 2+ (n = 1–6) by coupled cluster method

Quantum chemical calculations of the structures, stabilities, and interactions of the title series at the CCSD(T) theoretical level are performed. The analysis of binding energies, average binding energies, and fragmentation energies indicate that the n = 2 system is more stable than its neighbors. Topological analysis of the natural bond orbital, electron density deformation, integrated charge transfer, bond critical point properties, reduced density gradient analysis are performed to explore the nature of the interaction. The results show that the present Rn–Cu²⁺ interactions fall into intermediate interaction type with a pronounced covalent character.     

Toward more efficient CCSD(T) calculations of intermolecular interactions in model hydrogen-bonded and stacked dimers

Interaction energies of the model H-bonded complexes, the formamide and formamidine dimers, as well as the stacked formaldehyde and ethylene dimers are calculated by the coupled cluster CCSD(T) method. These systems serve as a model for H-bonded and stacking interactions, typical in molecules participating in biological systems. We use the optimized virtual orbital space (OVOS) technique, by which the dimension of the space of virtual orbitals in coupled cluster CCSD(T) calculations can be significantly reduced. We demonstrate that when the space of virtual orbitals is reduced to 50% of the full space, which means reducing computational demands by 1 order of magnitude, the interaction energies for both H-bonded and stacked dimers are affected by no more than 0.1 kcal/mol. This error is much smaller than the error when interaction energies are calculated using limited basis sets.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

Highly Accurate CCSD(T) and DFT–SAPT Stabilization Energies of H‐Bonded and Stacked Structures of the Uracil Dimer

The CCSD(T) interaction energies for the H‐bonded and stacked structures of the uracil dimer are determined at the aug‐cc‐pVDZ and aug‐cc‐pVTZ levels. On the basis of these calculations we can construct the CCSD(T) interaction energies at the complete basis set (CBS) limit. The most accurate energies, based either on direct extrapolation of the CCSD(T) correlation energies obtained with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets or on the sum of extrapolated MP2 interaction energies (from aug‐cc‐pVTZ and aug‐cc‐pVQZ basis sets) and extrapolated ΔCCSD(T) correction terms [difference between CCSD(T) and MP2 interaction energies] differ only slightly, which demonstrates the reliability and robustness of both techniques. The latter values, which represent new standards for the H‐bonding and stacking structures of the uracil dimer, differ from the previously published data for the S22 set by a small amount. This suggests that interaction energies of the S22 set are generated with chemical accuracy. The most accurate CCSD(T)/CBS interaction energies are compared with interaction energies obtained from various computational procedures, namely the SCS–MP2 (SCS: spin‐component‐scaled), SCS(MI)–MP2 (MI: molecular interaction), MP3, dispersion‐augmented DFT (DFT–D), M06–2X, and DFT–SAPT (SAPT: symmetry‐adapted perturbation theory) methods. Among these techniques, the best results are obtained with the SCS(MI)–MP2 method. Remarkably good binding energies are also obtained with the DFT–SAPT method. Both DFT techniques tested yield similarly good interaction energies. The large magnitude of the stacking energy for the uracil dimer, compared to that of the benzene dimer, is explained by attractive electrostatic interactions present in the stacked uracil dimer. These interactions force both subsystems to approach each other and the dispersion energy benefits from a shorter intersystem separation.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Observable Structures of Small Neutral and Anionic Gold Clusters

Since gold clusters have mostly been studied theoretically by using DFT calculations, more accurate studies are of importance. Thus, small neutral and anionic gold clusters (Au_n and Au_n^?, n=4–7) were investigated by means of coupled cluster with singles, doubles, and perturbative triple excitations [CCSD(T)] calculations with large basis sets, and some differences between DFT and CCSD(T) results are discussed. Interesting isomeric structures that have dangling atoms were obtained. Structures having dangling atoms appear to be stable up to n=4 for neutral gold clusters and up to n=7 for anionic clusters. The relative stabilities and electronic properties of some isomers and major structures are discussed on the basis of the CCSD(T) calculations. This accurate structure prediction of small gold clusters corresponding to experimental photoelectron spectral peaks is valuable in the field of atom‐scale materials science including nanocatalysts.     

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.