Low-lying electronic states of M(3)O(9)(-) and M(3)O(9)(2-) (M = Mo, W)

Multiple low-lying electronic states of M(3)O(9)(-) and M(3)O(9)(2-) (M = Mo, W) arise from the occupation of the near-degenerate low-lying virtual orbitals in the neutral clusters. We used density functional theory (DFT) and coupled cluster theory (CCSD(T)) with correlation consistent basis sets to study the structures and energetics of the electronic states of these anions. The adiabatic and vertical electron detachment energies (ADEs and VDEs) of the anionic clusters were calculated with 27 exchange-correlation functionals including one local spin density approximation functional, 13 generalized gradient approximation (GGA) functionals, and 13 hybrid GGA functionals, as well as the CCSD(T) method. For M(3)O(9)(-), CCSD(T) and nearly all of the DFT exchange-correlation functionals studied predict the (2)A(1) state arising from the Jahn-Teller distortion due to singly occupying the degenerate e' orbital to be lower in energy than the (2)A(1)' state arising from singly occupying the nondegenerate a(1)' orbital. For W(3)O(9)(-), the (2)A(1) state was predicted to have essentially the same energy as the (2)A(1)' state at the CCSD(T) level with core-valence correlation corrections included and to be higher in energy or essentially isoenergetic with most DFT methods. The calculated VDEs from the CCSD(T) method are in reasonable agreement with the experimental values for both electronic states if estimates for the corrections due to basis set incompleteness are included. For M(3)O(9)(2-), the singlet state arising from doubly occupying the nondegenerate a(1)' orbital was predicted to be the most stable state for both M = Mo and W. However, whereas M(3)O(9)(2-) was predicted to be less stable than M(3)O(9)(-), W(3)O(9)(2-) was predicted to be more stable than W(3)O(9)(-).     

Structural and electronic properties of reduced transition metal oxide clusters, M4O10 and M4O10- (M = Cr, W), from photoelectron spectroscopy and quantum chemical calculations

Anion photoelectron spectroscopy and quantum chemical calculations at the density functional theory (DFT), coupled cluster theory (CCSD(T)), and complete active space self-consistent field (CASSCF) theory levels are employed to study the reduced transition metal oxide clusters M(4)O(10)(-) (M = Cr, W) and their neutrals. Photoelectron spectra are obtained at 193 and 157 nm photon energies, revealing very different electronic structures for the Cr versus W oxide clusters. The electron affinity and HOMO-LUMO gap are measured to be 3.68 ± 0.05 and 0.7 eV, respectively, for the Cr(4)O(10) neutral cluster, as compared to 4.41 ± 0.04 and 1.3 eV for W(4)O(10). A comprehensive search is performed to determine the ground-state structures for M(4)O(10) and M(4)O(10)(-), in terms of geometry and electronic states by carefully examining the calculated relative energies at the DFT, CCSD(T), and CASSCF levels. The ground states of Cr(4)O(10) and Cr(4)O(10)(-) have tetrahedral structures similar to that of P(4)O(10) with the anion having a lower symmetry due to a Jahn-Teller distortion. The ground states of W(4)O(10) and W(4)O(10)(-) have butterfly shape structures, featuring two fused five-member rings with a metal-metal multiple bond between the central metal atoms. The much stronger WW bonding than the CrCr bonding is found to be the primary cause for the different ground state structures of the reduced Cr(4)O(10)(0/-) versus W(4)O(10)(0/-) oxide clusters. The photoelectron spectra are assigned by comparing the experimental and theoretical adiabatic and vertical electron detachment energies, further confirming the determination of the ground electronic states of M(4)O(10) and M(4)O(10)(-). The time-dependent DFT method is used to calculate the excitation energies of M(4)O(10). The TD-DFT results in combination with the self-consistently calculated vertical detachment energies for some of the excited states at the DFT and CCSD(T) levels are used to assign the higher energy bands. Accurate clustering energies and heats of formation of M(4)O(10) are calculated and used to calculate accurate reaction energies for the reduction of M(4)O(12) to M(4)O(10) by CH(3)OH, as well as for the oxidation of M(4)O(10) to M(4)O(12) by O(2). The performance of the DFT method with the B3LYP and BP86 functionals in the calculations of the relative energies, electron detachment energies, and excitation energies are evaluated, and the BP86 functional is found to give superior results for most of these energetic properties.     

Electronic structures and electron detachment energies of halogen substituted acetate anions, XCH2COO- (X=F,Cl,Br)

The electronic structures and the halogen inductive effects on the acetate anion were investigated in XCH2COO- (X=F,Cl,Br) by photoelectron spectroscopy (PES) and ab initio calculations. The PES spectra indicated that the electron binding energies increased in the order of FCl>Br. These systematic changes of detachment energy and IPs were explained by examining the charge redistributions upon detaching electrons.     

Electron affinities of uracil: microsolvation effects and polarizable continuum model

We present adiabatic electron affinities (AEAs) and the vertical detachment energies (VDEs) of the uracil molecule interacting with one to five water molecules. Credibility of MP2 and DFT/B3LYP calculations is supported by comparison with available benchmark CCSD(T) data. AEAs and VDEs obtained by MP2 and DFT/B3LYP methods copy trends of benchmark CCSD(T) results for the free uracil and uracil-water complexes in the gas phase being by 0.20 - 0.28 eV higher than CCSD(T) values depending on the particular structure of the complex. AEAs and VDEs from MP2 are underestimated by 0.09-0.15 eV. For the free uracil and uracil-(H(2)O)(n) (n = 1,2,3,5) complexes, we also consider the polarizable continuum model (PCM) and discuss the importance of the microsolvation when combined with PCM. AEAs and VDEs of uracil and uracil-water complexes enhance rapidly with increasing relative dielectric constant (ε) of the solvent. Highest AEAs and VDEs of the U(H(2)O)(5) complexes from B3LYP with ε = 78.4 are 2.03 and 2.81 eV, respectively, utilizing the correction from CCSD(T). Specific structural features of the microsolvated uracil-(H(2)O)(n) complexes and their anions are preserved also upon considering PCM in calculations of AEAs and VDEs.     

Benchmark calculations on the adiabatic ionization potentials of M-NH(3) (M=Na,Al,Ga,In,Cu,Ag)

The ground states of the M-NH(3) (M=Na,Al,Ga,In,Cu,Ag) complexes and their cations have been studied with density functional theory and coupled cluster [CCSD(T)] methods. The adiabatic ionization potentials (AIPs) of these complexes are calculated, and these are compared to results from high-resolution zero-electron kinetic energy photoelectron spectroscopy. By extrapolating the CCSD(T) energies to the complete basis set (CBS) limit and including the core-valence, scalar relativistic, spin-orbit, and zero-point corrections, the CCSD(T) method is shown to be able to predict the AIPs of these complexes to better than 6 meV or 0.15 kcal/mol. 27 exchange-correlation functionals, including one in the local density approximation, 13 in the generalized gradient approximation (GGA), and 13 with hybrid GGAs, were benchmarked in the calculations of the AIPs. The B1B95, mPW1PW91, B98, B97-1, PBE1PBE, O3LYP, TPSSh, and HCTH93 functionals give an average error of 0.1 eV for all the complexes studied, with the B98 functional alone yielding a maximum error of 0.1 eV. In addition, the calculated metal-ammonia harmonic stretching frequencies with the CCSD(T) method are in excellent agreement with their experimental values, whereas the B3LYP method tends to underestimate these stretching frequencies. The metal-ammonia binding energies were also calculated at the CCSD(T)/CBS level, and are in excellent agreement with the available experimental values considering the error limits, except for Ag-NH(3) and Ag(+)-NH(3), where the calculations predict stronger bond energies than measured by about 4 kcal/mol, just outside the experimental error bars of +/-3 kcal/mol.     

Small clusters of aluminum and tin: highly correlated calculations and validation of density functional procedures

We present results of molecular electronic structure treatments of multireference configuration interaction (MRCI) type for clusters Al(n) and Sn(n) in the range up to n = 4, and of coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) type in the range up to n = 10. Basis sets of quadruple zeta size are employed, computed energy differences, such as cohesive energies, E(coh), or dissociation energies for the removal of a single atom, D(e), differ from the complete basis set limit by only a few 0.01 eV. MRCI and CCSD(T) results are then compared to those obtained from density functional theory (DFT) treatments, which show that all computational procedures agree with the general features of D(e) and E(coh). The best agreement of DFT with CCSD(T) is found for the meta-GGA (generalized gradient approximation) TPSS (Tao, Perdew, Staroverov, Scuseria) for which D(e) differs from CCSD(T) by at most 0.15 eV for Al(n) and 0.21 eV for Sn(n). The GGA PBE (Perdew, Burke, Ernzerhof) is slightly poorer with maximum deviations of 0.23 and 0.24 eV, whereas hybrid functionals are not competitive with GGA and meta-GGA functionals. A general conclusion is that errors of D(e) and/or energy differences of isomers computed with DFT procedures may easily reach 0.2 eV and errors for cohesive energies E(coh) 0.1 eV.     

Accuracy and limitations of second-order many-body perturbation theory for predicting vertical detachment energies of solvated-electron clusters

Vertical electron detachment energies (VDEs) are calculated for a variety of (H(2)O)(n)(-) and (HF)(n)(-) isomers, using different electronic structure methodologies but focusing in particular on a comparison between second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster theory with noniterative triples, CCSD(T). For the surface-bound electrons that characterize small (H(2)O)(n)(-) clusters (n< or = 7), the correlation energy associated with the unpaired electron grows linearly as a function of the VDE but is unrelated to the number of monomers, n. In every example considered here, including strongly-bound "cavity" isomers of (H(2)O)(24)(-), the correlation energy associated with the unpaired electron is significantly smaller than that associated with typical valence electrons. As a result, the error in the MP2 detachment energy, as a fraction of the CCSD(T) value, approaches a limit of about -7% for (H(2)O)(n)(-) clusters with VDEs larger than about 0.4 eV. CCSD(T) detachment energies are bounded from below by MP2 values and from above by VDEs calculated using second-order many-body perturbation theory with molecular orbitals obtained from density functional theory. For a variety of both strongly- and weakly-bound isomers of (H(2)O)(20)(-) and (H(2)O)(24)(-), including both surface states and cavity states, these bounds afford typical error bars of +/-0.1 eV. We have found only one case where the Hartree-Fock and density functional orbitals differ qualitatively; in this case the aforementioned bounds lie 0.4 eV apart, and second-order perturbation theory may not be reliable.     

Enhanced stability by three-dimensional aromaticity of endohedrally doped clusters X(10)M(0/-) with X = Ge, Sn, Pb and M = Cu, Ag, Au

The group 14 clusters encapsulated by coinage metals in neutral and anionic states X(10)M(0/-) (X = Ge, Sn, Pb and M = Cu, Ag, Au) are investigated using quantum chemical calculations with the DFT/B3LYP functional and coupled-cluster CCSD(T) theory. Addition of transition metals into the empty cages forms high symmetry endohedral structures, except for Ge(10)Ag(0/-). In agreement with experiments available for X(10)Cu, the D(4d) global minima of the anions are calculated to be magic clusters with large frontier orbital gaps, high vertical and adiabatic detachment energies, and large embedding energies and binding energies as compared to those of the empty cages X(10)(2-). The enhanced stability of these magic clusters can be rationalized by the three-dimensional aromaticity.     

Calculation of electron detachment energies for water cluster anions: an appraisal of electronic structure methods, with application to (H2O)20- AND (H2O)24-

We present benchmark calculations of vertical electron detachment energies (VDEs) for various conformers of (H2O)n-, using both wave function and density functional methods, in sequences of increasingly diffuse Gaussian basis sets. For small clusters (n < or = 6), a systematic examination of VDE convergence reveals that it is possible to converge this quantity to within approximately 0.01 eV of the complete-basis limit, using a highly diffuse but otherwise economical Pople-style basis set of double-zeta quality, with 28 atom-centered basis functions per water molecule. Floating-center basis functions can be useful but are not required to obtain accurate VDEs. Second-order M?ller-Plesset perturbation (MP2) theory suffices to obtain VDEs that are within 0.05 eV of the results from both experiment and coupled-cluster theory, and which always err toward underbinding the extra electron. In contrast to these consistent predictions, VDEs calculated using density functional theory (DFT) vary widely, according to the fraction of Hartree-Fock exchange in a given functional. Common functionals such as BLYP and B3LYP overestimate the VDE by 0.2-0.5 eV, whereas a variant of Becke's "half and half" functional is much closer to coupled-cluster predictions. Exploratory calculations for (H2O)20- and (H2O)24- cast considerable doubt on earlier calculations that were used to assign the photoelectron spectra of these species to particular cluster isomers.     

The performance of density functional theory for LnF (Ln=Nd,Eu, Gd,Yb) and YbH

Different density functional theory (DFT) functionals have been evaluated by studying geometries and bond strengths of YbH, YbF, EuF, GdF, and NdF and compared with accurate CCSD(T) results and, when available, experiment. The agreement between the CCSD(T) results and experiment, when available, is good. The agreement is also good between bond strengths calculated at the DFT level using relativistic effective core potentials and the CCSD(T) results. However, the all-electron ADF calculations systematically overestimate binding energies. The geometries obtained by both the all-electron and the effective-core-potential-based DFT calculations are generally in good agreement with the CCSD(T) results.Contribution to the Björn Roos Honorary Issue     

Probing the Properties of Polynuclear Superhalogens without Halogen Ligand via ab Initio Calculations: A Case Study on Double‐Bridged [Mg2(CN)5]−1 Anions

An ab initio study of the superhalogen properties of eighteen binuclear double‐bridged [Mg₂(CN)₅]^?1 clusters is reported herein by using various theoretical methods. High‐level CCSD(T) results indicate that all the clusters possess strong superhalogen properties owing to their high vertical electron detachment energies (VDEs), which exceed 6.8 eV (highest: 8.15 eV). The outer valence Green's function method provides inaccurate relative VDE values; hence, this method is not suitable for this kind of polynuclear superhalogens. Both the HF and MP2 results are generally consistent with the CCSD(T) level regarding the relative VDE values and—especially interesting—the average values of the HF and MP2 VDEs are extremely close to the CCSD(T) results. The distributions of the extra electrons of the anions are mainly aggregated into the terminal CN units. These distributions are apparently different from those of previously reported triple‐bridged isomers and may be the reason for the decreased VDE values of the clusters. In addition, comparisons of the VDEs of binuclear and mononuclear superhalogens as well as studies of the thermodynamic stabilities with respect to the detachment of various CN^?1 ligands are also performed. These results confirm that polynuclear structures with pseudohalogen ligands can be considered as probable new superhalogens with enhanced properties.     

Electronic structure of the hydroxo and methoxo oxometalate anions MO3(OH)- and MO3(OCH3)- (M = Cr, Mo, and W)

The electronic structure of the mononuclear hydroxo MO3(OH)- and methoxo MO3(OCH3)- Group 6 oxometalate anions (M = Cr, Mo, and W) were examined by photodetachment photoelectron spectroscopy and electronic structure calculations at the density functional and CCSD(T) levels of theory. All of the anions exhibited high electron binding energies (>4.9 eV), with the lowest-energy detachment features arising from oxygen 2p-based orbitals. The combined experimental and theoretical results allowed the change in molecular orbital energy levels to be investigated as a function of metal (Cr, Mo, or W) and ligand (-OH, -OCH3). A number of fundamental thermodynamic properties of the anions and corresponding neutrals were predicted on the basis of the theoretical calculations. The calculations indicate high O-H bond dissociation energies for MO2(OR)(O-H) (R = H, CH3) and MO3(O-H), consistent with their high Br?nsted acidities (just below that of H2SO4 in the gas phase) and the high ionization energies of their conjugate base anions. This suggests that the corresponding radicals should readily abstract H atoms from organic molecules.     

The peculiar trend of cyclic perfluoroalkane electron affinities with increasing ring size

The adiabatic electron affinities (AEAs), vertical electron affinities (VEAs), and vertical detachment energies (VDEs) of cyclic perfluoroalkanes, c-C(n)F(2n) (n = 3-7), and their monotrifluoromethyl derivatives were computed using various pure and hybrid density functionals with DZP++ (polarization and diffuse function augmented double-zeta) basis sets. The theoretical AEA of c-C(4)F(8) at KMLYP/DZP++ is 0.70 eV, which exhibits satisfactory agreement with the 0.63 +/- 0.05 eV experimental value. The nonzero-point-corrected AEA of c-C(4)F(8) is predicted to be 0.41 eV at the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ level of theory, which shows a slight deviation of 0.11 eV from the KMLYP estimated value of 0.52 eV for the same. With the zero-point correction from the MP2/6-311G(d) [Gallup, G. A. Chem. Phys. Lett. 2004, 399, 206] level of theory combined with the CCSD(T)/aug-cc-pVTZ//MP2/aug-cc-pVTZ result, the most reliable estimate of AEA of c-C(4)F(8) is 0.60 eV. c-C(3)F(6)(-), c-C(4)F(8)(-), and c-C(5)F(10)(-) are unusual in preferring planar to near planar ring structures. The ZPE-corrected AEAs of c-C(n)F(2n) increase from n = 3 (0.24 eV) to n = 5 (0.77 eV), but then dramatically fall off to 0.40 eV for both n = 6 and n = 7. All of the other functionals predict the same trend. This is due to a change in the structural preference: C(s)() c-C(6)F(12)(-) and C(1) c-C(7)F(14)(-) are predicted to favor nonplanar rings, each with an exceptionally long C-F bond. (There also is a second, higher energy D3d minimum for C(6)F(12)(-).) The SOMOs as well as the spin density plots of the c-PFA radical anions reveal that the "extra" electron is largely localized on the unique F atoms in the larger n = 6 and n = 7 rings but is delocalized in the multiatom SOMOs of the three- to five-membered ring radical anions. The computed AEAs are much larger than the corresponding VEAs; the latter are not consistent with different functionals. The AEAs are substantially larger when a c-C(n)()F(2)(n)() fluorine is replaced by a -CF(3) group. This behavior is general; PFAs with tertiary C-F bonds have large AEAs. The VDEs for all the anions are substantial, ranging from 1.89 to 3.64 eV at the KMLYP/DZP++ level.     

Assessment of theoretical procedures for hydrogen-atom abstraction by chlorine, and related reactions

We have examined a number of hydrogen-abstraction reactions and assessed various theoretical procedures with regard to their performance for geometry optimization and for calculating barriers and reaction energies. We find that the BH&H-LYP and M05-2X procedures with the 6-31+G(d,p) basis set provide reasonable predictions for the geometries of the transition structures and also yield reasonable imaginary frequencies when compared with our benchmark QCISD/6-31+G(d,p) and CCSD(T)/6-311+G(3df,2p) values. For the calculation of barriers and reaction energies, M05-2X appears to be the most accurate of the hybrid functionals. The double-hybrid functionals, B2K-PLYP, UB2-PLYP-09, ROB2-PLYP, and DSD-B-LYP-D3, when used in combination with an augmented triple-zeta basis set, give very good agreement with the benchmark URCCSD(T)/aug-cc-pVQZ energies. We find that for wavefunction procedures, use of CCSD(T) in combination with an augmented triple-zeta quality basis set is required for the accurate prediction of barriers and reaction energies for these reactions.     

Nickel-catalyzed alkyl coupling reactions: evaluation of computational methods

The B3LYP, M06, M06L, M062X, MPW1K, and PBE1PBE DFT methods were evaluated for modeling nickel-catalyzed coupling reactions. The reaction consists of a nucleophilic attack by a carbanion equivalent on the nickel complex, S(N)2 attack by the anionic nickel complex on an alkyl halide, and reductive elimination of the coupled alkane product, regenerating the nickel catalyst. On the basis of CCSD(T)//DFT single-point energies, the B3LYP, M06, and PBE1PBE functionals were judged to generate the best ground state geometries. M06 energies are generally comparable or superior to B3LYP and PBE1PBE energies for transition state calculations. The MP2 and CCSD methods were also evaluated for single-point energies at the M06 geometries. The rate-determining step of this reaction was found to be nucleophilic attack of a L(2)NiR anion on the alkyl halide.     

羰基镍簇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校正, 从而预测合理的解离能.     

Exploring the rich energy landscape of sulfate-water clusters SO4(2-) (H2O)(n=3-7): an electronic structure approach

We present a reinvestigation of sulfate-water clusters SO4(2-) (H2O)(n=3-7), which involves several new aspects. Using a joint molecular mechanics/first principles approach, we perform exhaustive searches for stable cluster geometries, showing that the sulfate-water landscape is much richer than anticipated previously. We check the compatibility of the new structures with experiment by comparing vertical detachment energies (VDEs) calculated at the B3LYP/6-311++G** level of theory and determine the energetic ordering of the isomers at the RI-MP2/aug-cc-pVTZ level. Our results are bench-marked carefully against reference energies of estimated CCSD(T)/aug-cc-VTZ quality and VDEs of CCSD(T)/aug-cc-pVDZ quality. Furthermore, we calculate anharmonic vibrational corrections for up to the n = 6 clusters, which are shown to be significant for isomer energy ordering. We use energy decomposition analysis (EDA) based on the absolutely localized fragment (ALMO) expansion to gain chemical insight into the binding motifs.     

Conventional strain energies of 1,2-dihydroazete, 2,3-dihydroazete, 1,2-dihydrophosphete, and 2,3-dihydrophosphete

The conventional strain energies of 1,2-dihydroazete, 2,3-dihydroazete, 1,2-dihydrophosphete, and 2,3-dihydrophosphete are determined within the isodesmic, homodesmotic, and hyperhomodesmotic models. Optimum equilibrium geometries, harmonic vibrational frequencies, and corresponding electronic energies and zero-point vibrational energies are computed for all pertinent molecular systems using SCF theory, second-order perturbation theory, and density functional theory and employing the correlation consistent basis sets cc-pVDZ, cc-pVTZ, and cc-pVQZ. Single-point fourth-order perturbation theory, CCSD, and CCSD(T) calculations employing the cc-pVTZ and the cc-pVQZ basis sets are computed using the MP2/cc-pVTZ and MP2/cc-pVQZ optimized geometries, respectfully, to ascertain the contribution of higher order correlation. Three DFT functionals, B3LYP, wB97XD, and M06-2X, are employed to determine whether they can yield results similar to those obtained at the CCSD(T) level.     

Frontier orbital consistent quantum capping potential (FOC-QCP) for bulky ligand of transition metal complexes

Chemically reasonable models of PR3 (R = Me, Et, iPr, and tBu) were constructed to apply the post Hartree-Fock method to large transition metal complexes. In this model, R is replaced by the H atom including the frontier orbital consistent quantum capping potential (FOC-QCP) which reproduces the frontier orbital energy of PR3. The steric effect is incorporated by the new procedure named steric repulsion correction (SRC). To examine the performance of this FOC-QCP method with the SRC, the activation barriers and reaction energies of the reductive elimination reactions of C2H6 and H2 from M(R1)2(PR2(3))2 (M = Ni, Pd, or Pt; R1 = Me for R2 = Me, Et, or iPr, or R1 = H for R2 = tBu) were evaluated with the DFT[B3PW91], MP4(SDQ), and CCSD(T) methods. The FOC-QCP method reproduced well the DFT[B3PW91]- and MP4(SDQ)-calculated energy changes of the real complexes with PMe3. For more bulky phosphine, the SRC is important to present correct energy change, in which the MP2 method presents reliable steric repulsion correction like the CCSD(T) method because the systems calculated in the SRC do not include a transition metal element. The monomerization energy of [RhCl(PiPr3)2]2 and the coordination energies of CO, H2, N2, and C2H4 with [RhCl(PiPr3)2]2 were theoretically calculated by the CCSD(T) method combined with the FOC-QCP and the SRC. The CCSD(T)-calculated energies agree well with the experimental ones, indicating the excellent performance of the combination of the FOC-QCP with the SRC. On the other hand, the DFT[B3PW91]-calculated energies of the real complexes considerably deviate from the experimental ones.