Next Article in Journal
The bHLH Transcription Factors in Neural Development and Therapeutic Applications for Neurodegenerative Diseases
Next Article in Special Issue
New O- and N-N-Bridging Complexes of Tc(V), the Role of the Nitrogen Atom Position in Aromatic Rings: Reaction Mechanism, Spectroscopy, DTA, XRD and Hirshfeld Surface Analysis
Previous Article in Journal
Mitochondrial Fatty Acid β-Oxidation Disorders: From Disease to Lipidomic Studies—A Critical Review
Previous Article in Special Issue
Effects of Lewis Basicity and Acidity on σ-Hole Interactions in Carbon-Bearing Complexes: A Comparative Ab Initio Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics

College of Chemistry, Jilin University, Qianjin Street 2699, Changchun 130012, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(22), 13944; https://doi.org/10.3390/ijms232213944
Submission received: 25 October 2022 / Revised: 8 November 2022 / Accepted: 9 November 2022 / Published: 11 November 2022
(This article belongs to the Special Issue Non-covalent Interaction)

Abstract

:
When the bond lengths of 11 molecules containing van der Waals bonds are optimized by 572 methods and 20 basis sets, it is found that the best mean absolute deviations (MADs) of density-functional theory (DFT) methods are 0.005 Å (shown by APFD/6-311++G**), 0.007 Å (B2PLYPD3(Full)/aug-cc-pVQZ), and 0.010 Å (revDSDPBEP86/aug-cc-pVQZ), while the best MADs of ab initio methods are 0.008 Å (BD(T)/aug-cc-pVTZ) and 0.016 Å (MP4/aug-cc-pVQZ). Moreover, the best MADs calculated by 54 selected methods in combination with 60 other basis sets (such as 6-311++G, 6-31++G(3d′f,3p′d), and UGBS1V++) are not better. Therefore, these bond lengths can be calculated with extremely high accuracy by some special methods and basis sets, and CCSD(T) is also not as good as expected because its best MAD is only 0.023 Å (CCSD(T)/aug-cc-pVQZ).

1. Introduction

It is well known that experiment, theory, and computation are the three pillars of modern chemistry [1]. Like natural science, where experiment is the driving force for the establishment, testing, and development of a science, the core of chemistry is also experiment, such as synthetic chemistry, which has reached the stage of synthesizing any complexity (e.g., palytoxin) and any number (e.g., 200 million) of molecules. From experimental chemistry has developed theoretical chemistry making chemistry from empirical to rational, and theoretical chemistry becomes more and more important. As the center and the basics of theoretical chemistry, quantum chemistry was born when Heitler and London used quantum mechanics to explain why two neutral hydrogen atoms could form a chemical bond in 1927. It can be said that it is a science of describing chemical bonds at the molecular level. Therefore, chemical bonds are the core, the essence, and the golden key of chemistry. The developed result of theoretical chemistry is computational chemistry, which is also the executor and extension of theoretical chemistry, since computation is the bridge to connect theoretical models and experimental results. The development of chemical theories is inseparable from chemical experiments, although they can also be used to guide experiments. Computational chemistry combines theory and experiment, unites computer and chemistry, and uses computers as experimental tools. Theoretical chemistry and computational chemistry make chemistry more rigorous, quantitative, and mature, so that chemistry is no longer a purely experimental science.
The essence of chemistry is chemical reaction or chemical change, which is the breaking and forming of chemical bonds. It is also generally believed that structure determines everything, that is, molecular structure determines molecular properties, and molecular properties determine the chemical properties of a substance, so the properties of the substance are mainly determined by its molecular structure (e.g., sometimes external factors such as solvents are also important). However, a molecule is held together by chemical bond(s), or its structure depends on its chemical bond(s). That is to say, properties are the representation mode of structures, and structures are determinative. On the other hand, chemical bonds are the essence, and structures are their manifestation. Therefore, chemical bonds determine structures and properties by reasoning.
Chemical bonds are the attractions that join atoms in a compound, resulting in lower energy and greater stability. Although Pauling believed that Xenon might take place chemical reactions in 1933, it was not until Bartlett synthesized the first rare gas compound XePtF6 in 1962 [2] that the previous statement that rare gases could not participate in reactions was overturned. In 1873, van der Waals first discovered noncovalent interactions that helped him to reformulate the equation of state for real gases [3,4]. London used quantum mechanics to calculate dispersion forces in 1930 [5]. Bonding is inherently quantum mechanical in nature, so only quantum mechanics can provide valid theoretical tools to understand the nature of chemical bonds. For instance, the binding force of chemical bonds is formed by the wave action of electrons, or due to the special wave-particle duality, it is no wonder that it goes beyond the boundary of “same charges repel, different charges attract”. It is now recognized that it is convenient to classify chemical bonds into six types: covalent bond, ionic bond, coordination bond, metallic bond, hydrogen bond, and van der Waals bond. Similarly, the interactions between atoms can also be classified into six categories:
Two closed-shell neutral atoms, e.g., He-He, or a van der Waals bond;
Two open-shell neutral atoms, e.g., H-H, or a covalent bond in which two atoms share electron pair(s);
One closed-shell positive ion and one closed-shell negative ion, e.g., Na+−Cl, or an ionic bond of an electrostatic interaction;
One closed-shell atom or molecule with lone pair(s) of electrons and one or more open-shell atoms (including partially positively charged hydrogen(s) when bonded to an electronegative element) or molecules, e.g., a N−H·O hydrogen bond between oxygen and hydrogen, and a dative B−N bond or a special kind of covalent bonds in which one atom donates two electrons;
One open-shell ion (usually a positive ion) and one or several closed-shell ions or molecules, e.g., Mn+(X)m, or a coordination bond;
Many metal atoms aggregate together in which outer shell electrons break away from the core and move around all metal solids, e.g., a metallic bond.
The properties of chemical bonds are characterized by bond parameters such as bond length, bond angle, bond energy, bond order, and dipole moment, which can be measured by experiments directly or indirectly, or obtained from theoretical computations. Although direct experimental assessment of van der Waals interactions is intrinsically difficult [6], the notion that van der Waals bonds do exist can be seen from their very precisely measured bond lengths and bond energies. For example, the measured equilibrium bond length (re) and bond dissociation energy of NeKr are, respectively, 3.645444(13) Å [7] and 4.31 kJ/mol [8]. It is then doubtful whether computational chemistry can obtain these values theoretically and accurately, since theories are now usually considered to have achieved experimental accuracy. However, some popular density-functional theory (DFT) methods such as B3LYP and PBE1PBE have been found to be very poor in obtaining van der Waals bond lengths, and the “gold standard” [9] CCSD(T) is also not as good as expected when calculating molecular structures [10]. Additionally, van der Waals diatomics have not been very systematically studied in the references [11,12,13,14], so a more in-depth study is necessary, which is the main focus of this work.

2. Results and Discussion

2.1. Computational Results of 572 Methods and 20 Basis Sets

First, 11 rare gas molecules (diatomics) containing van der Waals bonds, represented below by 11-RG-Mols, are optimized by 572 methods and 20 basis sets. These methods include almost all DFT and ab initio methods known to us. Among these 20 basis sets, there are four of Ahlrichs’ [15], eight of Pople’s [16,17,18], and eight of Dunning’s [19,20] as they perform relatively better, are more commonly used, and can be used to calculate all these 11-RG-Mols. These 20 basis sets are numbered from BS01 to BS20 and are listed in Table S1 of the Supporting Information. The relative reliability of a method/basis set is based on its mean absolute deviation (MAD) value when computing the 11-RG-Mols. All their MADs are listed in Table S2 (their mean deviations (MDs) are also listed in Table S3). It is then found that the best MADs (0.005 Å, 0.006 Å, 0.007 Å, 0.008 Å, and 0.008 Å) are shown by the same DFT methods as in ref. [10] (i.e., APFD [14], B2PLYPD3 [21,22,23], and B2PLYPD3(Full)) in combination with several basis sets, i.e., shown respectively by APFD/6-311++G**, APFD/6-311++G(df,pd), B2PLYPD3(Full)/aug-cc-pVQZ, B2PLYPD3/aug-cc-pVQZ, and B2PLYPD3/aug-cc-pV5Z. The same best result of 0.008 Å is also shown by BD(T) [24], i.e., BD(T)/aug-cc-pVTZ, while the other best MADs of ab initio methods are 0.012 Å, 0.013 Å, 0.016 Å, and 0.016 Å, shown respectively by BD(T)(Full)/aug-cc-pVQZ, BD(T)/aug-cc-pVQZ, BD(T)/aug-cc-pV5Z, and MP4/aug-cc-pVQZ (i.e., other Møller-Plesset perturbation methods such as MP2, MP3, MP3(Full), MP4(Full), and MP5, whose MADs in combination with aug-cc-pVQZ are, respectively, 0.061 Å, 0.034 Å, 0.022 Å, 0.025 Å, and 0.027 Å are poorer). These data, and some other better or special methods, are listed in Table 1. The worst MAD (as large as 4.074 Å) is shown by G96B95/aug-cc-pVDZ. Unless otherwise stated, results of PM7, PM6, and UFF are not counted.
The methods which may be considered to be “accurate” (i.e., with MAD(s) between 0.010 Å and 0.020 Å) [25] or better (i.e., with MAD(s) even smaller than 0.010 Å) are shown by APFD [14], B2PLYPD3 [21,22,23], B2PLYPD3(Full), DSDPBEP86 [26,27], DSDPBEP86(Full), revDSDPBEP86 [28], revDSDPBEP86(Full), BD [24], BD(Full), MP4 [29,30], BD(T) [24], and BD(T)(Full). All these 12 methods are listed in Table 1. Their total number of method/basis sets is 33, its percentage relative to these 12 × 20 or 240 method/basis sets is 13.75%, and its percentage relative to all calculated 11436 method/basis sets is as low as 0.29%. Moreover, if the number (i.e., 792) of method/basis sets with MADs less than 0.100 Å (the MAD of UFF is 0.094 Å, also less than 0.100 Å) is counted, its percentage of 6.93% is also small. That is, most method/basis sets are very unreliable in computing these 11-RG-Mols. Methods with at least 10 times (i.e., at least 10 basis sets combined with a method with MADs less than 0.100 Å) include: 10 (BD(T)(Full), B2PLYPD, and M06-D3), 11 (PBEB95, LC-XaB95, and PBE1PBE-D3), 12 (LC-BRxB95, B2PLYP-D3, B2PLYP-D3(Full), and PBE1PBE-D3BJ), 13 (LGKCIS), 14 (wB97X, M06, LGB95, LC-G96B95, LC-RevTPSSB95, LC-TPSSB95, B2PLYPD3, B2PLYPD3(Full), and B3LYP-D3), 15 (LC-BB95, LC-LGB95, LC-mPWB95, LC-OB95, LC-PBEB95, LC-PBEhB95, LC-PKZBB95, LC-PW91B95, and LC-SB95), 16 (APFD and MN15L), 17 (MN15, DSDPBEP86, DSDPBEP86(Full), revDSDPBEP86, and revDSDPBEP86(Full)), and 18 times (PW6B95D3), which are all DFT methods except BD(T)(Full), e.g., the times of MP4(Full), CCSD(T)(Full), and QCISD(T)(Full) are all only 8.
In Table 2, the best and next best methods for these 20 basis sets are listed. It can be seen that their best MADs vary widely between 0.005 Å and 0.064 Å. The number of methods whose best MADs do not exceed 0.020 Å is 11. For these best MADs, the largest time is 8, shown by APFD, while the next largest time is 3, shown by DSDPBEP86. The other methods with 2 times include DSDPBEP86(Full), B2PLYPD3(Full), M06, and BD(T). For the next best methods, the largest times is 6, shown by DSDPBEP86(Full). For the basis sets whose best and next best MADs are both accurate or even better which are shown by methods of APFD, B2PLYPD3, B2PLYPD3(Full), DSDPBEP86, DSDPBEP86(Full), revDSDPBEP86(Full), and/or BD(T), only five Dunning’s basis sets can satisfy this requirement, i.e., aug-cc-pVTZ, cc-pVQZ, aug-cc-pVQZ, cc-pV5Z, and aug-cc-pV5Z. To obtain accurate or better bond lengths of the 11-RG-Mols, some special methods and basis sets must be used. As found before, diffuse functions are often important for calculating these 11-RG-Mols, e.g., the best MADs for 6-311G**, cc-pVDZ, and cc-pVTZ are, respectively, 0.064 Å, 0.056 Å, and 0.051 Å, while those for 6-311++G**, aug-cc-pVDZ, and aug-cc-pVTZ are, respectively, 0.005 Å, 0.025 Å, and 0.008 Å, which are much better. However, the situation is somewhat different for cc-pVQZ (with the best MAD of 0.016 Å), cc-pV5Z (0.012 Å), and Def2QZVP (0.012 Å).
The frequencies of these 36 methods listed in Table 1 in combination with all the 20 basis sets are also computed with their corresponding optimized bond lengths to see whether there are imaginary frequencies. An imaginary or a negative frequency means that this molecule is unstable (e.g., a transition state) judged from this method/basis set. As can be seen from their detailed data listed in Table S4, imaginary frequencies appear for all these 11-RG-Mols, but the number for one molecule ranges from 1 for 20Ne84Kr to 62 for 4He2. The total number of method/basis sets with imaginary frequency is 166, and its percentage relative to 7920 is only 2.10%, which is negligible, so nearly all geometries optimized by the corresponding method/basis sets are judged to be energy minima or stable conformations, meaning that these calculations can be considered successful. However, as listed in the Supporting Information of ref. [10], the number with imaginary frequency of 289 problematic molecules (almost all non-aromatic) is 69 when calculated by 26 method/basis sets (excluding PM6 and UFF results), and that of the total calculated 7514 method/basis sets is 741. The percentages of 69 versus 289 and 741 versus 7514 are, respectively, as high as 23.88% and 9.86%. It has been found that the number of imaginary frequencies which will exclude the stability of these two molecules is 1125, or its percentage relative to the total number of computed method/basis sets is 18.18% when planar benzene and nonplanar toluene are computed by more than one hundred methods and more than one hundred basis sets [31]. Therefore, the reason may be similar to calculating vibrational frequencies of aromatic molecules, i.e., it is also impossible to judge the stability of some non-aromatic molecules, including these van der Waals ones when using many method/basis sets of computational chemistry.

2.2. Computational Results of 54 Methods and 60 Other Basis Sets

Next, 54 relatively better methods are selected to optimize these 11-RG-Mols in combination with 60 other basis sets, including nearly all types for which they can be computed (e.g., STO-6G, TZVP, LanL2DZ, 6-31+G, CBSB7, 6-31++G(2df,2pd), 6-31++G(3d’,3p’), and UGBS), and their MADs are listed in Table S5 (their MDs are listed in Table S6). These 60 basis sets are numbered from BS21 to BS80, also listed in Table S1. If these DFT methods combined with UGBS and UGBS1V++ which can calculate only three molecules (i.e., like Dreiding, they can only calculate 4He2, Ne2, and HeNe) are not considered, it is seen that methods with “accurate” MADs only include BD(T), BD(T)(Full), and APFD, i.e., shown by BD(T)/6-31++G(3d,3p) (with MAD of 0.018 Å), BD(T)/6-31++G(3d’,3p’) (0.020 Å), BD(T)(Full)/6-31++G(3d’f,3p’d) (0.017 Å), APFD/6-311++G (0.017 Å), APFD/6-311++G* (0.011 Å), and APFD/6-311++G(df) (0.011 Å). Then, when they are used to optimize the 11-RG-Mols with APFD, the results of 6-311++G, 6-311++G*, and 6-311++G(df) are slightly worse than those of 6-311++G** and 6-311++G(df,pd) in Section 2.1, and other similar basis sets such as 6-31G* (with MAD of 0.354 Å), 6-31+G* (0.149 Å), 6-31++G** (0.108 Å), 6-31++G(3df,3pd) (0.173 Å), 6-311+G (0.034 Å), 6-311+G* (0.030 Å), 6-311+G** (0.049 Å), and CBSB7++ (0.028 Å) can be much worse. The worst MAD is 1.196 Å shown by MP3/CEP-31G or MP3(Full)/CEP-31G. The number of MADs smaller than 0.100 Å is 588, or its percentage relative to 3084 is 19.07%, much better than the overall 6.93% in Section 2.1. However, if only these 54 methods whose number with MAD smaller than 0.100 Å is 514 in Section 2.1 are counted, their percentage of 47.87% is much larger than this 19.07%, so the basis sets in Section 2.1 are generally better. It is noted that the number with MAD smaller than 0.100 Å of PW6B95D3 (with a maximum number of 18 for the former 20 basis sets) or APFD in combination with these 60 basis sets is 26, which is smaller than the maximum of 39 shown by M06, but as listed above, none of their MADs does not exceed 0.020 Å, except for six basis sets combined with APFD, BD(T), and BD(T)(Full).
The results of CCSD(T) [32] should be pointed out because it is generally considered as the only widely accepted gold-standard method. In Section 2.1, its best MAD is 0.023 Å, shown by CCSD(T)/aug-cc-pVQZ, as shown in ref. [10], while in this Section its best MAD is 0.026 Å, shown by CCSD(T)(Full)/6-31++G(3d’,3p’) (similar results for QCISD(T)). Therefore, the two best MADs are not accurate, nor better than revDSDPBEP86, BD(T), MP4, not to mention such DFT methods as APFD and B2PLYPD3(Full). In ref. [10], it is also found that not only are the best MADs of bond lengths (i.e., 0.021 Å shown by DSDPBEP86(Full)/Def2QZVP and MP2(Full)/cc-pVQZ) and bond angles (i.e., 2.3° shown by QCISD(T)/cc-pVTZ) not shown by CCSD(T), but also the best SUM1 (5.0 unit) and SUM2 (2.5 unit) of CCSD(T) shown by CCSD(T)/cc-pVTZ and CCSD(T)/6-311++G(3df,3pd) (shown also by QCISD(T)/cc-pVTZ and CCSD(Full)/cc-pVQZ for these same SUM1 and SUM2, with SUM1 and SUM2 measuring its overall reliability of a method/basis set calculating both bond length and bond angle) are larger than those (4.9 unit and 2.4 unit) of DSDPBEP86/cc-pVQZ, DSDPBEP86(Full)/cc-pVQZ, DSDPBEP86/aug-cc-pVQZ, DSDPBEP86/6-311++G(3df,3pd), DSDPBEP86(Full)/6-311++G(3df,3pd), and MP2/aug-cc-pVQZ. Therefore, it is no longer appropriate to call CCSD(T) the “gold standard” for the purpose of calculating molecular structures (at least for the lengths and angles of covalent and ionic bonds, and the lengths of van der Waals bonds).
Then, from this study and ref. [10], it is best to use DFT methods such as APFD, B2PLYPD3(Full), DSDPBEP86(Full), and revDSDPBEP86 to obtain better structures. For example, in ref. [10], it was found that DSDPBEP86, DSDPBEP86(Full), and APFD all performed very well (meaning that their SUM1s do not exceed 5.4 units and their SUM2s do not exceed 2.6 units) when used to optimize 179 molecules in combination with 6 basis sets (i.e., cc-pVTZ, cc-pVQZ, aug-cc-pVQZ, 6-311++G(3df,3pd), Def2TZVPP, and Def2QZVP) (this criterion is also met by 6 other DFT methods, i.e., OHSE1PBE, SOGGA11X, PBE1PBE, PBE1PBE-D3, PBE1PBE-D3BJ, and B3PW91-D3BJ). Although the number of these 6 basis sets which behave very well is only 4 for B2PLYPD3(Full), it is 2 for MP2 and CCSD(T). Not only are their computations (such as consumed memories and CPU times) much easier than CCSD(T), but also they are not inferior to CCSD(T), and can even be better.

3. Materials and Methods

These 11-RG-Mols include 4He2, Ne2, 40Ar2, 84Kr2, HeNe, HeAr, HeKr, 20Ne40Ar, 22Ne36Ar, 20Ne84Kr, and 40Ar84Kr, whose spectroscopic states are all X 1 g + or X1+, and whose adopted experimental equilibrium bond lengths (re) [7,11,12,13,33] are, respectively, 2.970 Å, 3.091 Å, 3.758 Å, 4.017(12) Å, 3.05 Å, 3.50 Å, 3.69 Å, 3.46360(12) Å, 3.46539(55) Å, 3.645444(13) Å, and 3.894358(50) Å, also listed in ref. [10]. As in ref. [10], their experimental bond lengths are usually used as their initial parameters to be optimized. These theoretically computed bond lengths are then compared with their corresponding experimental ones and the MAD ( MAD = ( 1 / 11 ) i = 1 11 r i , calc r i , exp [34]) of a method/basis set is obtained. These calculations are carried out by using Gaussian 16 program [35].

4. Conclusions

After the computations of the 11-RG-Mols containing van der Waals bonds with 572 methods and 80 basis sets, it is found that the best or next best MADs are usually shown by such methods as APFD, B2PLYPD3(Full), DSDPBEP86, and DSDPBEP86(Full), which are also very good for calculating the 179 problematic molecules with covalent and/or ionic bonds [10]. That is to say, to obtain very good structures for these three bond types, it is best to choose the above DFT methods in combination with some large or very large basis sets. CCSD(T) results do not meet the criterion of accurate bond lengths, so it is not suitable to be used as a comparing standard, as other ab initio methods like BD(T) and MP4 can be better than CCSD(T) in calculating these van der Waals bond lengths.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232213944/s1. Reference [36] are cited in the supplementary materials.

Author Contributions

Conceptualization, methodology, software, investigation, resources, data processing, and writing, Y.-L.Z.; earlier preliminary exploration, B.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support this study are available in the supplementary material and from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Y.L. Reaction Engineering; Chemical Industry Press: Beijing, China, 2020; p. 281. [Google Scholar]
  2. Bartlett, N. Xenon Hexafluoroplatinate(v) Xe+[PtF6]-. Proc. Chem. Soc. 1962, 218, 197–236. [Google Scholar] [CrossRef]
  3. Müller-Dethlefs, K.; Hobza, P. Noncovalent Interactions: A Challenge for Experiment and Theory. Chem. Rev. 2000, 100, 143–167. [Google Scholar] [CrossRef]
  4. Mahadevi, A.S.; Sastry, G.N. Cooperativity in Noncovalent Interactions. Chem. Rev. 2016, 116, 2775–2825. [Google Scholar] [CrossRef] [PubMed]
  5. Hobza, P.; Zahradník, R.; Müller-Dethlefs, K. The World of Non-Covalent Interactions: 2006. Collect. Czech. Chem. Commun. 2006, 71, 443–531. [Google Scholar] [CrossRef] [Green Version]
  6. Stöhr, M.; Voorhis, T.V.; Tkatchenko, A. Theory and practice of modeling van der Waals interactions in electronic-structure calculations. Chem. Soc. Rev. 2019, 48, 4118–4154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Hüttner, W. (Ed.) Landolt-Börnstein, New Series, Vol. II/24A; Springer: Berlin, Germany, 1998. [Google Scholar]
  8. Cheng, J.P. Bond Dissociation Energies. In CRC Handbook of Chemistry and Physics, 101st ed.; Rumble, J.R., Jr., Ed.; Taylor & Francis Group, LLC: Boca Raton, FL, USA, 2020; pp. 9–77. [Google Scholar]
  9. Dubecký, M.; Mitas, L.; Jurečka, P. Noncovalent Interactions by Quantum Monte Carlo. Chem. Rev. 2016, 116, 5188–5215. [Google Scholar] [CrossRef] [PubMed]
  10. Zhang, Y.L.; Wang, F.L.; Ren, A.M. Reliability of computed molecular structures. J. Comput. Chem. 2022, 43, 465–476. [Google Scholar] [CrossRef]
  11. Kannemann, F.O.; Becke, A.D. Van der Waals Interactions in Density-Functional Theory: Rare-Gas Diatomics. J. Chem. Theory Comput. 2009, 5, 719–727. [Google Scholar] [CrossRef]
  12. Kristyán, S.; Pulay, P. Can (semi)local density functional theory account for the London dispersion forces? Chem. Phys. Lett. 1994, 229, 175–180. [Google Scholar] [CrossRef]
  13. Tang, K.T.; Toennies, J.P. The van der Waals potentials between all the rare gas atoms from He to Rn. J. Chem. Phys. 2003, 118, 4976–4983. [Google Scholar] [CrossRef]
  14. Austin, A.; Petersson, G.A.; Frisch, M.J.; Dobek, F.J.; Scalmani, G.; Throssell, K. A Density Functional with Spherical Atom Dispersion Terms. J. Chem. Theory Comput. 2012, 8, 4989–5007. [Google Scholar] [CrossRef] [PubMed]
  15. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef] [PubMed]
  16. Hehre, W.J.; Ditchfield, R.; Pople, J.A. Self-consistent molecular orbital methods. XII. Further extensions of Gaussian-type basis sets for use in molecular orbital studies of organic molecules. J. Chem. Phys. 1972, 56, 2257–2261. [Google Scholar] [CrossRef]
  17. Hariharan, P.C.; Pople, J.A. The Influence of Polarization Functions on Molecular Orbital Hydrogenation Energies. Theoret. Chim. Acta 1973, 28, 213–222. [Google Scholar] [CrossRef]
  18. Krishnan, R.; Binkley, J.S.; Seeger, R.; Pople, J.A. Self-consistent molecular orbital methods. XX. A basis set for correlated wave functions. J. Chem. Phys. 1980, 72, 650–654. [Google Scholar] [CrossRef]
  19. Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  20. Kendall, R.A.; Dunning, T.H., Jr.; Harrison, R.J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796–6806. [Google Scholar] [CrossRef] [Green Version]
  21. Grimme, S. Semiempirical hybrid density functional with perturbative second-order correlation. J. Chem. Phys. 2006, 124, 034108. [Google Scholar] [CrossRef] [Green Version]
  22. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef]
  23. Goerigk, L.; Grimme, S. Efficient and Accurate Double-Hybrid-Meta-GGA Density Functionals─Evaluation with the Extended GMTKN30 Database for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput. 2011, 7, 291–309. [Google Scholar] [CrossRef]
  24. Handy, N.C.; Pople, J.A.; Head-Gordon, M.; Raghavachari, K.; Trucks, G.W. Size-consistent Brueckner theory limited to double substitutions. Chem. Phys. Lett. 1989, 164, 185–192. [Google Scholar] [CrossRef]
  25. Foresman, J.B.; Frisch, A.E. Exploring Chemistry with Electronic Structure Methods, 2nd ed.; Gaussian, Inc.: Pittsburgh, PA, USA, 1996; p. 118. [Google Scholar]
  26. Kozuch, S.; Martin, J.M.L. DSD-PBEP86: In search of the best double-hybrid DFT with spin-component scaled MP2 and dispersion corrections. Phys. Chem. Chem. Phys. 2011, 13, 20104–20107. [Google Scholar] [CrossRef]
  27. Kozuch, S.; Martin, J.M.L. Spin-Component-Scaled Double Hybrids: An Extensive Search for the Best Fifth-Rung Functionals Blending DFT and Perturbation Theory. J. Comput. Chem. 2013, 34, 2327–2344. [Google Scholar] [CrossRef]
  28. Santra, G.; Sylvetsky, N.; Martin, J.M.L. Minimally Empirical Double-Hybrid Functionals Trained against the GMTKN55 Database: revDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4. J. Phys. Chem. A 2019, 123, 5129–5143. [Google Scholar] [CrossRef] [PubMed]
  29. Raghavachari, K.; Pople, J.A. Approximate 4th-order perturbation-theory of electron correlation energy. Int. J. Quantum Chem. 1978, 14, 91–100. [Google Scholar] [CrossRef]
  30. Raghavachari, K.; Frisch, M.J.; Pople, J.A. Contribution of triple substitutions to the electron correlation energy in fourth-order perturbation theory. J. Chem. Phys. 1980, 72, 4244–4245. [Google Scholar] [CrossRef]
  31. Zhang, Y.L.; Zhang, G. Vibrational frequencies and determination of stability of aromatic molecules. Chin. J. Org. Chem. 2014, 34, 178–189. [Google Scholar] [CrossRef]
  32. Raghavachari, K.; Trucks, G.W.; Pople, J.A.; Head-Gordon, M. A fifth-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 1989, 157, 479–483. [Google Scholar] [CrossRef]
  33. Huber, K.P.; Herzberg, G. Molecular Spectra and Molecular Structure, IV. Constants of Diatomic Molecules; Van Nostrand Reinhold Company: New York, NY, USA, 1979. [Google Scholar]
  34. Levine, I.N. Quantum Chemistry, 6th ed.; Pearson Education, Inc.: Upper Saddle River, NJ, USA, 2009; p. 700. [Google Scholar]
  35. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision A.03; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  36. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision A.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
Table 1. MADs (Å) of some methods calculating the 11-RG-Mols.
Table 1. MADs (Å) of some methods calculating the 11-RG-Mols.
D2TPD2QPP++dp++3d++df++3dfaug-Taug-Qcc-5aug-5Ave
HF0.7911.3090.9140.7260.9130.7261.1051.2591.3081.3081.036
MP20.2190.1980.2730.0910.2410.0760.0570.0610.1530.0630.143
B2PLYPD30.0380.0530.1830.0990.1550.0950.0140.0080.0240.0080.068
B2PLYPD3(Full)0.0370.0380.1790.0850.1540.0830.0200.0070.0120.0100.063
DSDPBEP860.0180.0350.0520.0580.0480.0600.0280.0260.0130.0290.037
DSDPBEP86(Full)0.0180.0250.0500.0610.0400.0660.0720.0410.0140.0560.044
revDSDPBEP860.0320.0780.0670.0550.0580.0550.0150.0100.0340.0170.042
revDSDPBEP86(Full)0.0230.0590.0640.0600.0540.0600.0490.0190.0220.0330.044
CCSD(Full)0.2450.2250.2740.0920.2540.0780.0380.0250.1420.0330.140
BD0.2490.2150.2800.0900.2410.0710.0300.0190.0970.0220.131
BD(Full)0.2490.2120.2800.0900.2410.0710.0300.0190.0970.0220.131
MP40.2440.2030.2620.0980.2250.0710.0310.0160.1150.0270.129
CCSD(T)0.2500.1670.2600.0980.2260.0720.0300.0230.1070.0260.126
CCSD(T)(Full)0.2010.1660.2440.0560.2040.0490.0600.0250.0830.0600.115
BD(T)0.2070.1500.2240.0480.1870.0380.0080.0130.0980.0160.099
BD(T)(Full)0.1800.1780.2350.0660.2710.0430.0380.0120.0740.0270.112
APFD0.0940.0120.0050.0540.0060.0560.0230.0120.0120.0130.029
wB97X0.1020.0650.0820.0760.0830.0770.0700.0510.0670.0780.075
B3LYP0.9821.9101.4161.0691.3801.0741.7721.9101.8281.9901.533
PBE1PBE0.1480.1550.1480.1220.1420.1210.1470.1540.1530.1530.144
M060.0650.1410.0720.0700.0720.0700.0680.1700.1400.1520.102
M06L0.1210.1400.1270.0740.1270.0730.0660.0990.1100.0760.101
MN150.0940.0810.0750.0660.0750.0660.0660.0720.0820.0880.076
MN15L0.0720.0860.0760.0700.0760.0620.0610.0810.0860.0890.076
PW6B95D30.0810.0890.0820.0690.0820.0690.0840.0880.0890.0900.082
PBEB950.1190.1090.0570.0650.0570.0650.0840.0920.1090.0970.086
LC-PBEB950.0600.0990.0860.0840.0860.0840.0860.0960.0990.1000.088
B3LYP-D20.1700.0970.0930.1100.0940.1110.0910.1000.1000.0970.106
B3LYP-D30.1060.0640.0680.0870.0690.0870.0430.0650.0710.0660.073
M06-D30.0680.1720.1360.0730.1360.0730.1230.1720.1710.1180.124
M06L-D30.1220.1410.1280.0980.1280.0740.0670.1010.1110.1280.110
LC-wPBE-D30.0490.1570.1620.1150.1610.1120.1400.1730.1550.1660.139
PBE1PBE-D30.1210.0890.0960.1110.0970.1120.0840.0900.0920.0880.098
PBE1PBE-D3BJ0.1120.0870.0840.1020.0840.1020.0820.0900.0860.0870.091
PBEPBE-D3BJ0.1500.0970.0980.1180.0980.1190.0900.0950.1040.0970.107
B3LYP-D3BJ0.0520.0830.1650.1500.1650.1540.1510.0770.0730.1640.123
Ave0.1640.2000.1970.1290.1870.1240.1400.1490.1700.1580.162
The D2TP, D2QPP, ++dp, ++3d, ++df, ++3df, aug-T, aug-Q, cc-5, and aug-5 denote, respectively, Def2TZVP, Def2QZVPP, 6-311++G**, 6-311++G(3d,3p), 6-311++G(df,pd), 6-311++G(3df,3pd), aug-cc-pVTZ, aug-cc-pVQZ, cc-pV5Z, and aug-cc-pV5Z. Ave is abbreviation of average for these considered 36 methods or 10 basis sets. A green or blue value means, respectively, that this MAD is smaller than 0.020 Å or between 0.020~0.100 Å, while a black one means that this value is larger than 0.100 Å.
Table 2. The best and next best MADs (Å) and methods for the 20 basis sets.
Table 2. The best and next best MADs (Å) and methods for the 20 basis sets.
Basis SetThe BestThe Next Best
MADMethodMADMethod
Def2TZVP0.018DSDPBEP86,
or DSDPBEP86(Full)
0.023revDSDPBEP86(Full)
Def2TZVPP0.027DSDPBEP86(Full)0.029DSDPBEP86
Def2QZVP0.012APFD0.029DSDPBEP86(Full)
Def2QZVPP0.012APFD0.025DSDPBEP86(Full)
6-311G**0.064M060.069M06-D3
6-311++G**0.005APFD0.049wPBEhB95
6-311++G(2d,2p)0.036DSDPBEP860.040DSDPBEP86(Full)
6-311++G(3d,3p)0.054APFD0.055revDSDPBEP86
6-311++G(df,pd)0.006APFD0.040DSDPBEP86(Full)
6-311++G(2df,2pd)0.036revDSDPBEP86(Full)0.038DSDPBEP86
6-311++G(3df,3pd)0.038BD(T)0.043BD(T)(Full)
6-311++G(3d2f,3p2d)0.012APFD0.061PBEB95
cc-pVDZ0.055M06HF0.056M06HF-D3
aug-cc-pVDZ0.025APFD0.036DSDPBEP86(Full)
cc-pVTZ0.051M060.053M06-D3
aug-cc-pVTZ0.008BD(T)0.014B2PLYPD3
cc-pVQZ0.016DSDPBEP860.018DSDPBEP86(Full),
or revDSDPBEP86(Full)
aug-cc-pVQZ0.007B2PLYPD3(Full)0.008B2PLYPD3
cc-pV5Z0.012APFD,
or B2PLYPD3(Full)
0.013DSDPBEP86
aug-cc-pV5Z0.008B2PLYPD30.010B2PLYPD3(Full)
As in Table 1, a green or blue value means, respectively, that this MAD is smaller than 0.020 Å or between 0.020~0.100 Å. , while a black one means that this value is larger than 0.100 Å.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.-L.; Li, B. Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics. Int. J. Mol. Sci. 2022, 23, 13944. https://doi.org/10.3390/ijms232213944

AMA Style

Zhang Y-L, Li B. Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics. International Journal of Molecular Sciences. 2022; 23(22):13944. https://doi.org/10.3390/ijms232213944

Chicago/Turabian Style

Zhang, Yi-Liang, and Bin Li. 2022. "Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics" International Journal of Molecular Sciences 23, no. 22: 13944. https://doi.org/10.3390/ijms232213944

APA Style

Zhang, Y. -L., & Li, B. (2022). Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics. International Journal of Molecular Sciences, 23(22), 13944. https://doi.org/10.3390/ijms232213944

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop