Ab initio calculations on halogen-bonded complexes and comparison with density functional methods

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.     

Halogen bonds with benzene: An assessment of DFT functionals

The performance of an extensive set of density functional theory functionals has been tested against CCSD(T) and MP2 results, extrapolated to the complete basis set (CBS) limit, for the interaction of either DCl or DBr (D = H, HCC, F, and NC) with the aromatic system of benzene. It was found that double hybrid functionals explicitly including dispersion, that is, B2PLYPD and mPW2PLYPD, provide the better agreement with the CCSD(T)/CBS results on both energies and equilibrium geometry, indicating the importance of dispersive contributions in determining this interaction. Among the less expensive functionals, the better performance is provided by the ωB97X and M062X functionals, while the ωB97XD and B97D functionals are shown to work very well for bromine complexes but not so well for chlorine complexes. © 2013 Wiley Periodicals, Inc.     

The influence of peripheral ligand bulk on nitrogen activation by three‐coordinate molybdenum complexes—A theoretical study using the ONIOM method

Electronic structure methods have been combined with the ONIOM approach to carry out a comprehensive study of the effect of ligand bulk on the activation of dinitrogen with three‐coordinate molybdenum complexes. Calculations were performed with both density functional and CCSD(T) methods. Our results show that not only is there expected destabilization of the intermediate on the pathway due to direct steric interactions of the bulky groups, but also there is significant electronic destabilization as the size of the ligand increases. This latter destabilization is due to the inability of the molecule to accommodate a rotated amide group bound to the molybdenum once the amide reaches a certain size. This destabilization also leads to a clear preference for the triplet intermediate (rather than the singlet intermediate) for bulky substituents which is in agreement with experiment. Overall, the calculated reaction profile for the bulky substituents shows a good correlation with the available experimental data. © 2009 Wiley Periodicals, Inc. J Comput Chem 2009     

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.     

Coordination properties of the oxime analogue of glycine to Cu(II)

The coordination of Cu2+ by glyoxilic acid oxime (gao)--the oxime analogue of glycine amino acid--and its deprotonated (gao- and gao2-) species has been studied with different density functional methods. Single-point calculations have also been carried out at the single- and double- (triple) excitation coupled-cluster (CCSD(T)) level of theory. The isomers studied involve coordination of Cu2+ to electron-rich sites (O,N) of neutral, anionic, and dianionic gao species in different conformations. In contrast to Cu2+-glycine, for which the ground-state structure is bidentate with the CO2(-) terminus of zwitterionic glycine, for Cu2+-gao the most stable isomer shows monodentate binding of Cu2+ with the carbonylic oxygen of the neutral form. The most stable complexes of Cu2+ interacting with deprotonated gao species (gao- and gao2-) also take place through the carboxylic oxygens but in a bidentate manner. The results with different functionals show that, for these open shell (Cu2+-L) systems, the relative stability of complexes with different coordination environments (and so, different spin distribution) can be quite sensitive to the amount of "Hartree-Fock" exchange included in the functional. Among all the functionals tested in this work, the BHandHLYP is the one that better compares to CCSD(T) results.     

Reassessment of the Mechanisms of Thermal C−H Bond Activation of Methane by Cationic Magnesium Oxides: A Critical Evaluation of the Suitability of Different Density Functionals

The mechanisms of the thermal reactions of the two iconic magnesium oxide cations MgO^.+ and Mg₂O₂^.+ with methane have been re-evaluated at the CCSD(T)/CBS//CCSD/def2-TZVP level of theory. For the reaction of MgO^.+ with CH₄, only the classical hydrogen-atom transfer (HAT) was found; in contrast, for the Mg₂O₂^.+/CH₄ couple, both HAT and proton-coupled electron-transfer (PCET) exist as mechanistic variants. In order to evaluate the suitability of density functional theory (DFT) methods, the reactions were computed by using 27 density functionals. The results obtained demonstrate that the various DFT methods often deliver rather different results for both geometric and energetic features. As to the prediction of the apparent barriers, pure functionals give the largest mean absolute errors. BMK, ωB97XD, and the double-hybrid functional mPW2PLYP were confirmed to come closest to the results provided by CCSD(T)/CBS. Thus, mechanistic conclusions based on a single DFT method should be viewed with great caution. In summary, this study may assist in the selection of a suitable quantum chemical method to unravel the mechanistic details of C−H bond activation by charged metal oxides.     

Thermodynamics of titanium and vanadium reduction in non-aqueous environment calculated at various levels of theory

Reduction of titanium and vanadium compounds is a process accompanying the activation of coordinative olefin polymerization catalysts. Four density functional theory (DFT) functionals, coupled cluster with single, double, and perturbative triple excitations method CCSD(T) as well as complete active-space second-order perturbation theory method CASPT2 with a complete active-space self-consistent field CASSCF reference wave function were applied to investigate the thermodynamics of titanium and vanadium reduction. The performance of these theoretical methods was assessed and compared with experimental values. The calculations indicate that vanadium(IV) chloride is more easily reduced by trimethylaluminum than the corresponding titanium compound; the energies of reaction calculated at the CCSD(T) level are equal -57.21 and -33.10 kcal/mol, respectively. The calculations deal with the redox reactions of metal chlorides in the gas phase, rather than solvated ions in the aqueous solution. This approach may be more appropriate for olefin polymerization, usually carried out in nonpolar solvents.     

Assessment of the accuracy of density functionals for prediction of relative energies and geometries of low-lying isomers of water hexamers

Water hexamers provide a critical testing ground for validating potential energy surface predictions because they contain structural motifs not present in smaller clusters. We tested the ability of 11 density functionals (four of which are local and seven of which are nonlocal) to accurately predict the relative energies of a series of low-lying water hexamers, relative to the CCSD(T)/aug'-cc-pVTZ level of theory, where CCSD(T) denotes coupled cluster theory with an interative treatment of single and double excitations and a quasi-perturbative treatment of connected triple excitations. Five of the density functionals were tested with two different basis sets, making a total of 16 levels of density functional theory (DFT) tested. When single-point energy calculations are carried out on geometries obtained with second-order M?ller-Plesset perturbation theory (MP2), only three density functionals, M06-L, M05-2X, and M06-2X, are able to correctly predict the relative energy ordering of the hexamers. These three functionals predict that the range of energies spanned by the six isomers is 3.2-5.6 kcal/mol, whereas the other eight functionals predict ranges of 1.0-2.4 kcal/mol; the benchmark value for this range is 3.1 kcal/mol. When the hexamers are optimized at each level of theory, all methods are able to reproduce the MP2 geometries well for all isomers except the boat and bag isomers, and DFT optimization changes the energy ordering for seven of the 16 methods tested. The addition of zero-point energy changes the energy ordering for all of the density functionals studied except for M05-2X and M06-2X. The variation in relative energies predicted by the different methods highlights the necessity for exercising caution in the choice of density functionals used in future studies. Of the 11 density functionals tested, the most accurate results for energies were obtained with the PWB6K, MPWB1K, and M05-2X functionals.     

Oxidative addition of the fluoromethane C-F bond to Pd. An ab initio benchmark and DFT validation study

We have computed a state-of-the-art benchmark potential energy surface (PES) for two reaction pathways (oxidative insertion, OxIn, and S(N)2) for oxidative addition of the fluoromethane C-F bond to the palladium atom and have used this to evaluate the performance of 26 popular density functionals, covering LDA, GGA, meta-GGA, and hybrid density functionals, for describing these reactions. The ab initio benchmark is obtained by exploring the PES using a hierarchical series of ab initio methods (HF, MP2, CCSD, CCSD(T)) in combination with a hierarchical series of seven Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account through a full four-component all-electron approach. Our best estimate of kinetic and thermodynamic parameters is -5.3 (-6.1) kcal/mol for the formation of the reactant complex, 27.8 (25.4) kcal/mol for the activation energy for oxidative insertion (OxIn) relative to the separate reactants, 37.5 (31.8) kcal/mol for the activation energy for the alternative S(N)2 pathway, and -6.4 (-7.8) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). Our work highlights the importance of sufficient higher angular momentum polarization functions for correctly describing metal-d-electron correlation. Best overall agreement with our ab initio benchmark is obtained by functionals from all three categories, GGA, meta-GGA, and hybrid DFT, with mean absolute errors of 1.4-2.7 kcal/mol and errors in activation energies ranging from 0.3 to 2.8 kcal/mol. The B3LYP functional compares very well with a slight underestimation of the overall barrier for OxIn by -0.9 kcal/mol. For comparison, the well-known BLYP functional underestimates the overall barrier by -10.1 kcal/mol. The relative performance of these two functionals is inverted with respect to previous findings for the insertion of Pd into the C-H and C-C bonds. However, all major functionals yield correct trends and qualitative features of the PES, in particular, a clear preference for the OxIn over the alternative S(N)2 pathway.     

Testing the performance of density functionals for the calculation of energetic properties of complex-forming radical-molecule reactions

The performance of some popular and some more recent density functional methods for the calculation of energies of stationary points on the potential surfaces of radical-molecule reactions was examined. The functionals studied are B3-LYP, BH&H, BH&H-LYP, MPW1K, MPWB1K, TPSS, TPSSh, BB1K, M05 and M05-2X, in conjunction with nine different AO basis sets. The reaction energies, barrier heights and the relative energies of the pre-and post-reaction complexes were compared with those obtained at the CCSD(T)/CBS limit for the reactions of OH radicals with HOOH and CH₃OOH. Very poor barrier heights are provided by the B3-LYP, TPSS and TPSSh functionals. The best overall performance was obtained with the BB1K, MPW1K and MPWB1K functionals. In these reactions all of the studied functionals provide converged results only if they are used with large basis sets like aug-cc-pVTZ and def2-TZVP. The data show that before relying on a functional for a specific reaction, it is desirable to make some test calculations on the performance. The same functional can predict some relative energies very well and some others very poorly even in systems including chemically similar reactants.     

How accurate are electronic structure methods for actinoid chemistry?

The CASPT2, CCSD, and CCSD(T) levels of wave function theory and seven density functionals were tested against experiment for predicting the ionization potentials and bond dissociation energies of actinoid monoxides and dioxides with their cations. The goal is to guide future work by enabling the choice of an appropriate method when performing calculations on actinoid-containing systems. We found that four density functionals, namely MPW3LYP, B3LYP, M05, and M06, and three levels of wave function theory, namely CASPT2, CCSD, and CCSD(T), give similar mean unsigned errors for actinoid?Coxygen bond energies and for ionization potentials of actinoid oxides and their cations.     

Theoretical Aspects of Hydrolysis of Peptide Bonds by Zinc Metalloenzymes

Activation and reaction energies for four model systems capturing the essential physicochemical features of the hydrolysis of the peptide bond have been calculated at various level of theory, including the presumably accurate CCSD(T) calculations. The models studied covered a part of the spectrum encountered in biological systems: the hydrolysis in the absence of metal ions (represented by formamide and Ala–Ala) and the hydrolysis in the presence of one and two zinc(II) ions, mimicking the active sites of mono‐ and dizinc metallopeptidases, respectively (by using thermolysin and glutamate carboxypeptidase II as the model catalytic systems and formamide as the model substrate). The results obtained using CCSD(T)/def2‐TZVP and CCSD(T)/aug‐cc‐pVTZ calculations were used as the benchmark values to which the set of cheaper methods, such as (RI‐)DFT, (RI‐)MP2, and SCS‐MP2, were referenced. It was shown that deviations of 3–5 kcal mol^?1 (translating to 2–3 orders in reaction constants) with respect to the reference CCSD(T) barriers are frequently encountered for many correlated methods and most of studied DFT functionals. It has been concluded that from the set of wave‐function methods, both MP2 and SCS‐MP2 methods can be recommended for smaller models (measured by the mean absolute deviation of the activation barriers over the four systems studied), whereas among the popular DFT functionals, B3LYP and especially M06‐2X are likely to be reasonable choices for calculating the activation barriers of zinc metallopeptidases. Finally, with the model of glutamate carboxypeptidase II, issues related to the convergence of the calculated barriers with the size of the model system used as the representative of the enzyme active site were addressed. The intricacies related to system truncation are demonstrated, and suggest that the correlated wave‐function methods may suffer from problems, such as intramolecular BSSE, which make their usage for the larger system questionable. Altogether, the presented data should contribute to efforts to understand enzymatic catalysis more deeply and to gain control of the accuracy and deficiencies of the available theoretical methods and computational approaches.     

Assessment of density functional methods for reaction energetics: Iridium‐catalyzed water oxidation as case study

We investigate basis set convergence for a series of density functional theory (DFT) functionals (both hybrid and nonhybrid) and compare to coupled‐cluster with single and double excitations and perturbative triples [CCSD(T)] benchmark calculations. The case studied is the energetics of the water oxidation reaction by an iridium‐oxo complex. Complexation energies for the reactants and products complexes as well as the transition state (TS) energy are considered. Contrary to the expectation of relatively weak basis set dependence for DFT, the basis set effects are large, for example, more than 10 kcal mol^?1 difference from converged basis for the activation energy with “small” basis sets (DZ/6‐31G** for Ir/other atoms, or SVP) and still more than 6 kcal mol^?1 for def2‐TZVPP/6‐31G**. Inclusion of the dispersion correction in DFT‐D3 schemes affects the energies of reactant complex (RC), TS, and product complex (PC) by almost the same amount; it significantly improves the complexation energy (the formation of RC), but has little effect on the activation energy with respect to RC. With converged basis, some pure GGAs (PBE‐D3, BP86‐D3) as well as the hybrid functional B3LYP‐D3 are very accurate compared to benchmark CCSD(T) calculations. © 2012 Wiley Periodicals, Inc.     

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)(-).     

Electronic structure benchmark calculations of CO2 fixing elementary chemical steps in RuBisCO using the projector-based embedding approach

Ribulose 1,5-bisphosphate carboxylase-oxygenase (RuBisCO) is the main enzyme involved in atmospheric carbon dioxide (CO₂) fixation in the biosphere. This enzyme catalyzes a set of five chemical steps that take place in the same active-site within magnesium (II) coordination sphere. Here, a set of electronic structure benchmark calculations have been carried out on a reaction path proposed by Gready et al. by means of the projector-based embedding approach. Activation and reaction energies for all main steps catalyzed by RuBisCO have been calculated at the MP2, SCS-MP2, CCSD, and CCSD(T)/aug-cc-pVDZ and cc-pVDZ levels of theory. The treatment of the magnesium cation with post-HF methods is explored to determine the nature of its involvement in the mechanism. With the high-level ab initio values as a reference, we tested the performance of a set of density functional theory (DFT) exchange-correlation (xc) functionals in reproducing the reaction energetics of RuBisCO carboxylase activity on a set of model fragments. Different DFT xc-functionals show large variation in activation and reaction energies. Activation and reaction energies computed at the B3LYP level are close to the reference SCS-MP2 results for carboxylation, hydration and protonation reactions. However, for the carbon–carbon bond dissociation reaction, B3LYP and other functionals give results that differ significantly from the ab initio reference values. The results show the applicability of the projector-based embedding approach to metalloenzymes. This technique removes the uncertainty associated with the selection of different DFT xc-functionals and so can overcome some of inherent limitations of DFT calculations, complementing, and potentially adding to modeling of enzyme reaction mechanisms with DFT methods.     

Benchmark calculations on the electron detachment energies of MO3* and M2O6* (M = Cr, Mo, W)

Neutral and anionic molecules of the monomers and dimers of the group VIB transition metal oxides (MO3 and M2O6) were studied with density functional theory (DFT) and coupled cluster CCSD(T) theory. Franck-Condon simulations of the photoelectron spectra were carried out for the transition from the ground state of the anion to that of the neutral molecule. Molecular structures from the DFT and CCSD(T) methods are compared. Electron detachment energies reported in the literature were evaluated. The calculated adiabatic and vertical electron detachment energies (ADEs and VDEs) were compared with the experimental results. CCSD(T) gives results within 0.12 eV for the ADEs. CCSD(T) predicts VDEs that are in error by as much as 0.3 eV for M = Cr. DFT hybrid functionals were found to give poor results for the ADEs and VDEs for M = Cr due to the substantial amount of multireference character in the wavefunction, whereas the pure DFT functionals give superior results. For M = Mo and W, excellent agreement was found for both CCSD(T) and many DFT fucntionals. The BP86 functional yields the best overall results for the VDEs of all the metal oxide clusters considered. Heats of formation calculated at the CCSD(T) level extrapolated to the complete basis set limit are also in good agreement with available experimental data.     

CO oxidation catalyzed by a single gold atom: benchmark calculations and the performance of DFT methods

Quantum chemical calculations were carried out on CO oxidation catalyzed by a single gold atom. To investigate the performance of density functional theory (DFT) methods, 42 DFT functionals have been evaluated and compared with high-level wavefunction based methods. It was found that in order to obtain accurate results the functionals used must treat long range interaction well. The double-hybrid mPW2PLYP and B2PLYP functionals are the two functionals with best overall performance. CAM-B3LYP, a long range corrected hybrid GGA functional, also performs well. On the other hand, the popular B3LYP, PW91, and PBE functionals do not show good performance and the performance of the latter two are even at the bottom of the 42 functionals. Our accurate results calculated at the CCSD(T)/aug-cc-pVTZ//mPW2PLYP/aug-cc-pVTZ level of theory indicate that Au atom is a good catalysis for CO oxidation. The reaction follows the following mechanism where CO and O(2) adsorb on Au atom forming an Au(OCOO) intermediate and subsequently O(2) transfer one oxygen atom to CO to form CO(2) and AuO. Then AuO reacts with CO to form another CO(2) to complete the catalytic cycle. The overall energy barrier at 0 K for the first CO oxidation step (Au + CO + O(2)→ AuO + CO(2)) is just 4.8 kcal mol(-1), and that for the second CO oxidation step (AuO + CO → Au + CO(2)) is just 1.6 kcal mol(-1).     

Density functional theory and post-Hartree-Fock studies on molecular structure and harmonic vibrational spectrum of formaldehyde

Density functional theory (DFT) with the Becke's three-parameter exchange correlation functional and the functional of Lee, Yang and Parr, gradient-corrected functionals of Perdew, and Perdew and Wang [the DFT(B3LYP), DFT(B3P86) and DFT(B3PW91) methods, respectively], and several levels of conventional ab initio post-Hartree-Fock theory (second- and fourth-order perturbation theory M?ller-Plesset MP2 and MP4(SDTQ), coupled cluster with the single and double excitations (CCSD), and CCSD with perturbative triple excitation [CCSD(T)], configuration interaction with the single and double excitations [CISD], and quadratic configuration interaction method [QCISD(T)], using several basis sets [ranging from a simple 6-31G(d,p) basis set to a 6-311+ +G(3df, 2pd) one], were applied to study of the molecular structure (geometrical parameters, rotational constants, dipole moment) and harmonized infrared (IR) spectrum of formaldehyde (CH₂O). High-level ab initio methods CCSD(T) and QCISD(T) with the 6-311+ +G(3df, 2pd) predict correctly molecular parameters, vibrational harmonic wavenumbers and the shifts of the harmonic IR spectrum of ¹²CH₂ ¹⁶O upon isotopic substitution. Received: 30 January 1997 / Accepted: 7 May 1997     

Benchmark of density functional theory methods on the prediction of bond energies and bond distances of noble-gas containing molecules

We have tested three pure density functional theory (DFT) functionals, BLYP, MPWPW91, MPWB95, and ten hybrid DFT functionals, B3LYP, B3P86, B98, MPW1B95, MPW1PW91, BMK, M05-2X, M06-2X, B2GP-PLYP, and DSD-BLYP with a series of commonly used basis sets on the performance of predicting the bond energies and bond distances of 31 small neutral noble-gas containing molecules. The reference structures were obtained using the CCSD(T)∕aug-cc-pVTZ theory and the reference energies were based on the calculation at the CCSD(T)∕CBS level. While in general the hybrid functionals performed significantly better than the pure functionals, our tests showed a range of performance by these hybrid functionals. For the bond energies, the MPW1B95∕6-311+G(2df,2pd), BMK∕aug-cc-pVTZ, B2GP-PLYP∕aug-cc-pVTZ, and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 2.0-2.3 kcal∕mol per molecule. For the bond distances, the MPW1B95∕6-311+G(2df,2pd), MPW1PW91∕6-311+G(2df,2pd), and B3P86∕6-311+G(2df,2pd), DSD-BLYP∕6-311+G(2df,2pd), and DSD-BLYP∕aug-cc-pVTZ methods stood out with mean unsigned errors of 0.008-0.013 A? per bond. The current study showed that a careful selection of DFT functionals is very important in the study of noble-gas chemistry, and the most recommended methods are MPW1B95∕6-311+G(2df,2pd) and DSD-BLYP∕aug-cc-pVTZ.