Top-Down versus Bottom-Up Approaches for σ-Functionals Based on the Approximate Exchange KernelClick to copy article linkArticle link copied!
- Yannick LemkeYannick LemkeDepartment of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 5-13, D-81377 Munich, GermanyMore by Yannick Lemke
- Christian Ochsenfeld*Christian Ochsenfeld*Email: [email protected]Department of Chemistry, Ludwig-Maximilians-Universität München, Butenandtstr. 5-13, D-81377 Munich, GermanyMax-Planck-Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, GermanyMore by Christian Ochsenfeld
Abstract
Recently, we investigated a number of so-called σ- and τ-functionals based on the adiabatic-connection fluctuation–dissipation theorem (ACFDT); particularly, extensions of the random phase approximation (RPA) with inclusion of an exchange kernel in the form of an antisymmetrized Hartree kernel. One of these functionals, based upon the approximate exchange kernel (AXK) of Bates and Furche, leads to a nonlinear contribution of the spline function used within σ-functionals, which we previously avoided through the introduction of a simplified “top-down” approach in which the σ-functional modification is inserted a posteriori following the analytic coupling strength integration within the framework of the ACFDT and which was shown to provide excellent performance for the GMTKN55 database when using hybrid PBE0 reference orbitals. In this work, we examine the analytic “bottom-up” approach in which the spline function is inserted a priori, i.e., before evaluation of the analytic coupling strength integral. The new bottom-up functionals, denoted σ↑AXK, considerably improve upon their top-down counterparts for problems dominated by self-interaction and delocalization errors. Despite a small loss of accuracy for noncovalent interactions, the σ↑AXK@PBE0 functionals comprehensively outperform regular σ-functionals, scaled σ-functionals, and the previously derived σ+SOSEX- and τ-functionals in the WTMAD-1 and WTMAD-2 metrics of the GMTKN55 database.
This publication is licensed under
License Summary*
You are free to share(copy and redistribute) this article in any medium or format and to adapt(remix, transform, and build upon) the material for any purpose, even commercially within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format and to adapt(remix, transform, and build upon) the material for any purpose, even commercially within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format and to adapt(remix, transform, and build upon) the material for any purpose, even commercially within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
Special Issue
Published as part of The Journal of Physical Chemistry A special issue “Trygve Helgaker Festschrift”.
Introduction
Theory
Fluctuation–Dissipation Theorem
Direct Random Phase Approximation and σ-Functionals
Inclusion of Exchange Kernels and τ-Functionals
Approximate Exchange Kernel
Top-Down “σ↓AXK” Correlation Energy

Bottom-Up “σ↑AXK” Correlation Energy
Computational Details
Results and Discussion
Optimization of Bottom-Up σ-Functionals
Figure 1
Figure 1. MAE (solid) and wMAE (striped) for the ASCDB database using different σ-functionals with PBE and PBE0 reference orbitals.
Evaluation Compared to Top-Down σ-Functionals
Figure 2
Figure 2. (a) MAE, (b) WTMAD-1, and (c) WTMAD-2 for the GMTKN55 database and its subcategories using different σ-functionals with PBE reference orbitals.
Figure 3
Figure 3. (a) MAE, (b) WTMAD-1, and (c) WTMAD-2 for the GMTKN55 database and its subcategories using different σ-functionals with PBE0 reference orbitals.
Figure 4
Figure 4. Error distribution for the W4–17RE benchmark set using the top-down σ↓AXK(A2) and bottom-up σ↑AXK(A2) functionals for (a) PBE and (b) PBE0 reference orbitals.
Physical Interpretation of Top-Down versus Bottom-Up Approaches
Figure 5
Figure 5. Weight functions of the Hartree (a, c) and exchange (b, d) kernels for σ↓AXK-functionals (a, b) and σ↑AXK-functionals (c, d) for PBE0 reference orbitals.
@PBE | @PBE0 | |
---|---|---|
dRPA | –20.768 | –18.714 |
σ(W1) | –23.399 | –16.080 |
scσ(S1) | –2.110 | –1.512 |
dRPA+AXK | –2.219 | –1.737 |
σ↓AXK(A1) | –22.174 | –9.121 |
σ↓AXK(A2) | –24.143 | –10.391 |
σ↑AXK(A1) | –8.851 | –5.595 |
σ↑AXK(A2) | –7.602 | –5.540 |
refb | DFT | dRPA | σ | scσ | AXK | σ↓AXK | σ↑AXK | |||
---|---|---|---|---|---|---|---|---|---|---|
W1 | S1 | A1 | A2 | A1 | A2 | |||||
De (cm–1) | 28.868 | 44.050 | 11.087 | 11.669 | 18.607 | 11.851 | 13.296 | 10.695 | 19.054 | 16.614 |
Re (Å) | 3.1007 | 3.0748 | 3.2029 | 3.1598 | 3.1927 | 3.1895 | 3.1672 | 3.2045 | 3.0790 | 3.0875 |
Figure 6
Figure 6. CBS-extrapolated neon dimer interaction energies computed using PBE reference orbitals for (a) KS-DFT, (b) σ-functionals and scaled σ-functionals, (c) top-down σ↓AXK-functionals, and (d) bottom-up σ↑AXK-functionals. For reference, corrected CCSD(T) values from ref (97) are also included; data between points has been interpolated using cubic splines.
Conclusions
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpca.4c05289.
MAEs for all subsets of the ASCDB and GMTKN55 databases for σ↑AXK functionals; weight functions wH and wx for PBE reference orbitals; and spline coefficients for σ↑AXK functionals (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.
Acknowledgments
The authors thank Dr. Jörg Kussmann (LMU Munich) for providing access to a development version of the FermiONs++ program package. Financial support was provided by the Deutsche Forschungsgemeinschaft (DFG) through grant CRC 325 “Assembly Controlled Chemical Photocatalysis” (grant no. 444632635). CO acknowledges additional financial support as Max-Planck-Fellow at the Max Planck Institute for Solid State Research (MPI-FKF) Stuttgart.
References
This article references 102 other publications.
- 1Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864– B871, DOI: 10.1103/PhysRev.136.B864Google ScholarThere is no corresponding record for this reference.
- 2Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133– A1138, DOI: 10.1103/PhysRev.140.A1133Google ScholarThere is no corresponding record for this reference.
- 3Mardirossian, N.; Head-Gordon, M. Thirty years of density functional theory in computational chemistry: An overview and extensive assessment of 200 density functionals. Mol. Phys. 2017, 115, 2315– 2372, DOI: 10.1080/00268976.2017.1333644Google Scholar3https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVCltb3O&md5=ba27d707ee3f5fcdd949644d3d2cbd5eThirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionalsMardirossian, Narbe; Head-Gordon, MartinMolecular Physics (2017), 115 (19), 2315-2372CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)In the past 30 years, Kohn-Sham d. functional theory has emerged as the most popular electronic structure method in computational chem. To assess the ever-increasing no. of approx. exchange-correlation functionals, this review benchmarks a total of 200 d. functionals on a mol. database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochem., and barrier heights. In addn., the evolution of non-empirical and semi-empirical d. functional design is reviewed, and guidelines are provided for the proper and effective use of d. functionals. The most promising functional considered is ωB97M-V, a range-sepd. hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chem. applications, and the principal remaining limitations are assocd. with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.
- 4Perdew, J. P.; Schmidt, K. Jacob’s ladder of density functional approximations for the exchange-correlation energy. AIP Conf. Proc. 2001, 577, 1– 20, DOI: 10.1063/1.1390175Google Scholar4https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXmtVGnsb0%253D&md5=cac0386cbbd9df5b30216875abcd1d24Jacob's ladder of density functional approximations for the exchange-correlation energyPerdew, John P.; Schmidt, KarlaAIP Conference Proceedings (2001), 577 (Density Functional Theory and Its Application to Materials), 1-20CODEN: APCPCS; ISSN:0094-243X. (American Institute of Physics)The ground-state energy and d. of a many-electron system are often calcd. by Kohn-Sham d. functional theory. We describe a ladder of approxns. for the exchange-correlation energy as a functional of the electron d. At the lowest rung of this ladder, the contribution to the energy from a vol. element of 3-dimensional space is detd. by the local d. there. Higher rungs or levels incorporate increasingly complex ingredients constructed from the d. or the Kohn-Sham orbitals in or around this vol. element. We identify which addnl. exact conditions can be satisfied at each level, and discuss the extent to which the functionals at each level may be constructed without empirical input. We also discuss the research that remains to be done at the exact-exchange level, and present our "dreams of a final theory". "Jacob left Beer-sheba and went toward Haran. He came to a certain place and stayed there for the night, because the sun had set. Taking one of the stones of the place, he put it under his head and lay down in that place. And he dreamed that there was a ladder set up on the earth, the top of it reaching to heaven; and the angels of God were ascending and descending on it.".
- 5Becke, A. D.; Santra, G.; Martin, J. M. L. A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database. J. Chem. Phys. 2023, 158, 151103, DOI: 10.1063/5.0141238Google ScholarThere is no corresponding record for this reference.
- 6Kállay, M. Linear-scaling implementation of the direct random-phase approximation. J. Chem. Phys. 2015, 142, 204105, DOI: 10.1063/1.4921542Google Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXpt1Sjsb4%253D&md5=a66b1691368dd94596a8052c569de4d2Linear-scaling implementation of the direct random-phase approximationKallay, MihalyJournal of Chemical Physics (2015), 142 (20), 204105/1-204105/16CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We report the linear-scaling implementation of the direct RPA (dRPA) for closed-shell mol. systems. As a bonus, linear-scaling algorithms are also presented for the second-order screened exchange extension of dRPA as well as for the second-order Moller-Plesset (MP2) method and its spin-scaled variants. Our approach is based on an incremental scheme which is an extension of our previous local correlation method [Rolik et al., J. Chem. Phys. 139, 094105 (2013)]. The approach extensively uses local natural orbitals to reduce the size of the MO basis of local correlation domains. In addn., we also demonstrate that using natural auxiliary functions [M. Kallay, J. Chem. Phys. 141, 244113 (2014)], the size of the auxiliary basis of the domains and thus that of the three-center Coulomb integral lists can be reduced by an order of magnitude, which results in significant savings in computation time. The new approach is validated by extensive test calcns. for energies and energy differences. Our benchmark calcns. also demonstrate that the new method enables dRPA calcns. for mols. with more than 1000 atoms and 10 000 basis functions on a single processor. (c) 2015 American Institute of Physics.
- 7Schurkus, H. F.; Ochsenfeld, C. Communication: An effective linear-scaling atomic-orbital reformulation of the random-phase approximation using a contracted double-Laplace transformation. J. Chem. Phys. 2016, 144, 031101, DOI: 10.1063/1.4939841Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xht1Ciurw%253D&md5=b7f953b877f5fea31e2508e9f914be77Communication: An effective linear-scaling atomic-orbital reformulation of the random-phase approximation using a contracted double-Laplace transformationSchurkus, Henry F.; Ochsenfeld, ChristianJournal of Chemical Physics (2016), 144 (3), 031101/1-031101/5CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)An at.-orbital (AO) reformulation of the RPA correlation energy is presented allowing to reduce the steep computational scaling to linear, so that large systems can be studied on simple desktop computers with fully numerically controlled accuracy. Our AO-RPA formulation introduces a contracted double-Laplace transform and employs the overlap-metric resoln.-of-the-identity. First timings of our pilot code illustrate the reduced scaling with systems comprising up to 1262 atoms and 10 090 basis functions. (c) 2016 American Institute of Physics.
- 8Luenser, A.; Schurkus, H. F.; Ochsenfeld, C. Vanishing-Overhead Linear-Scaling Random Phase Approximation by Cholesky Decomposition and an Attenuated Coulomb-Metric. J. Chem. Theory Comput. 2017, 13, 1647, DOI: 10.1021/acs.jctc.6b01235Google Scholar8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXjs1yhsb0%253D&md5=a7a01e028f23cc0441ba2d752ebe200aVanishing-Overhead Linear-Scaling Random Phase Approximation by Cholesky Decomposition and an Attenuated Coulomb-MetricLuenser, Arne; Schurkus, Henry F.; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2017), 13 (4), 1647-1655CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A reformulation of the RPA within the Resoln.-of-the-Identity (RI) scheme is presented, that is competitive to canonical MO RI-RPA already for small- to medium-sized mols. For electronically sparse systems drastic speedups due to the reduced scaling behavior compared to the MO formulation are demonstrated. Our reformulation is based on two ideas which are independently useful: First, a Cholesky decompn. of d. matrixes that reduces the scaling with basis set size for a fixed-size mol. by one order, leading to massive performance improvements. Second, replacement of the overlap RI metric used in the original AO-RPA by an attenuated Coulomb metric. Accuracy is significantly improved compared to the overlap metric, while locality and sparsity of the integrals are retained, as is the effective linear scaling behavior.
- 9Graf, D.; Beuerle, M.; Schurkus, H. F.; Luenser, A.; Savasci, G.; Ochsenfeld, C. Accurate and Efficient Parallel Implementation of an Effective Linear-Scaling Direct Random Phase Approximation Method. J. Chem. Theory Comput. 2018, 14, 2505, DOI: 10.1021/acs.jctc.8b00177Google Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXnsFWisrg%253D&md5=b886c4aaf9c1b9e4a4f8482c453aa68bAccurate and Efficient Parallel Implementation of an Effective Linear-Scaling Direct Random Phase Approximation MethodGraf, Daniel; Beuerle, Matthias; Schurkus, Henry F.; Luenser, Arne; Savasci, Goekcen; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2018), 14 (5), 2505-2515CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)An efficient algorithm for calcg. the RPA correlation energy is presented that is as accurate as the canonical MO resoln.-of-the-identity RPA (RI-RPA) with the important advantage of an effective linear-scaling behavior (instead of quartic) for large systems due to a formulation in the local AO space. The high accuracy is achieved by utilizing optimized minimax integration schemes and the local Coulomb metric attenuated by the complementary error function for the RI approxn. The memory bottleneck of former AO (AO)-RI-RPA implementations is addressed by precontraction of the large 3-center integral matrix with the Cholesky factors of the ground state d. reducing the memory requirements of that matrix by a factor of (Nbasis/Nocc). Furthermore, we present a parallel implementation of our method, which not only leads to faster RPA correlation energy calcns. but also to a scalable decrease in memory requirements, opening the door for investigations of large mols. even on small- to medium-sized computing clusters. Although it is known that AO methods are highly efficient for extended systems, where sparsity allows for reaching the linear-scaling regime, we show that our work also extends the applicability when considering highly delocalized systems for which no linear scaling can be achieved. As an example, the interlayer distance of two covalent org. framework pore fragments (comprising 384 atoms in total) is analyzed.
- 10Beuerle, M.; Graf, D.; Schurkus, H. F.; Ochsenfeld, C. Efficient calculation of beyond RPA correlation energies in the dielectric matrix formalism. J. Chem. Phys. 2018, 148, 204104, DOI: 10.1063/1.5025938Google Scholar10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtVWlu7%252FM&md5=e7fa58c478a2a8532793dc946e50e57dEfficient calculation of beyond RPA correlation energies in the dielectric matrix formalismBeuerle, Matthias; Graf, Daniel; Schurkus, Henry F.; Ochsenfeld, ChristianJournal of Chemical Physics (2018), 148 (20), 204104/1-204104/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present efficient methods to calc. beyond RPA correlation energies for mol. systems with up to 500 atoms. To reduce the computational cost, we employ the resoln.-of-the-identity and a double-Laplace transform of the non-interacting polarization propagator in conjunction with an AO formalism. Further improvements are achieved using integral screening and the introduction of Cholesky decompd. densities. Our methods are applicable to the dielec. matrix formalism of RPA including second-order screened exchange (RPA-SOSEX), the RPA electron-hole time-dependent Hartree-Fock (RPA-eh-TDHF) approxn., and RPA renormalized perturbation theory using an approx. exchange kernel (RPA-AXK). We give an application of our methodol. by presenting RPA-SOSEX benchmark results for the L7 test set of large, dispersion dominated mols., yielding a mean abs. error below 1 kcal/mol. The present work enables calcg. beyond RPA correlation energies for significantly larger mols. than possible to date, thereby extending the applicability of these methods to a wider range of chem. systems. (c) 2018 American Institute of Physics.
- 11Drontschenko, V.; Graf, D.; Laqua, H.; Ochsenfeld, C. Lagrangian-Based Minimal-Overhead Batching Scheme for the Efficient Integral-Direct Evaluation of the RPA Correlation Energy. J. Chem. Theory Comput. 2021, 17, 5623, DOI: 10.1021/acs.jctc.1c00494Google Scholar11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXhvVGqt7rN&md5=f4d0ae3663784bc63f97e43c80a919f8Lagrangian-Based Minimal-Overhead Batching Scheme for the Efficient Integral-Direct Evaluation of the RPA Correlation EnergyDrontschenko, Viktoria; Graf, Daniel; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2021), 17 (9), 5623-5634CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A highly memory-efficient integral-direct RPA method based on our ω-CDGD-RI-RPA method is presented that completely alleviates the memory bottleneck of storing the multidimensional three-center integral tensor, which severely limited the tractable system sizes. Based on a Lagrangian formulation, we introduce an optimized batching scheme over the auxiliary and basis-function indexes, which allows to compute the optimal no. of batches for a given amt. of system memory, while minimizing the batching overhead. Thus, our optimized batching constitutes the best tradeoff between program runtime and memory demand. Within this batching scheme, the half-transformed three-center integral tensor BiμM is recomputed for each batch of auxiliary and basis functions. This allows the computation of systems that were out of reach before. The largest system within this work consists of a DNA fragment comprising 1052 atoms and 11 230 basis functions calcd. on a single node, which emphasizes the new possibilities of our integral-direct RPA method.
- 12Harl, J.; Kresse, G. Cohesive energy curves for noble gas solids calculated by adiabatic connection fluctuation-dissipation theory. Phys. Rev. B 2008, 77, 045136, DOI: 10.1103/PhysRevB.77.045136Google Scholar12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhslKrtLY%253D&md5=645e7ee23ef184cb921e0800974565fdCohesive energy curves for noble gas solids calculated by adiabatic connection fluctuation-dissipation theoryHarl, Judith; Kresse, GeorgPhysical Review B: Condensed Matter and Materials Physics (2008), 77 (4), 045136/1-045136/8CODEN: PRBMDO; ISSN:1098-0121. (American Physical Society)We present first-principles calcns. for the fcc noble gas solids Ne, Ar, and Kr applying the adiabatic connection fluctuation-dissipation theorem (ACFDT) to evaluate the correlation energy. The ACFDT allows us to describe long-range correlation effects including London dispersion or van der Waals interaction on top of conventional d. functional theory calcns. Even within the RPA, the typical 1/V2 vol. dependence for the cohesive energy of the noble gas solids is reproduced, and equil. cohesive energies and lattice consts. are improved compared to d. functional theory calcns. Furthermore, we present atomization energies for H2, N2, and O2 within the same post-d.-functional-theory framework, finding an excellent agreement with previously published data.
- 13Su, H.; Wu, Q.; Wang, H.; Wang, H. An assessment of the random-phase approximation functional and characteristics analysis for noncovalent cation−π interactions. Phys. Chem. Chem. Phys. 2017, 19, 26014– 26021, DOI: 10.1039/C7CP04504BGoogle ScholarThere is no corresponding record for this reference.
- 14Modrzejewski, M.; Yourdkhani, S.; Klimeš, J. Random Phase Approximation Applied to Many-Body Noncovalent Systems. J. Chem. Theory Comput. 2020, 16, 427– 442, DOI: 10.1021/acs.jctc.9b00979Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3MfgtlWgtw%253D%253D&md5=257b7378ea6502c717458e2ebf398a11Random Phase Approximation Applied to Many-Body Noncovalent SystemsModrzejewski Marcin; Yourdkhani Sirous; Klimes Jiri; Modrzejewski MarcinJournal of chemical theory and computation (2020), 16 (1), 427-442 ISSN:.The random phase approximation (RPA) has received considerable interest in the field of modeling systems where noncovalent interactions are important. Its advantages over widely used density functional theory (DFT) approximations are the exact treatment of exchange and the description of long-range correlation. In this work, we address two open questions related to RPA. First, we demonstrate how accurately RPA describes nonadditive interactions encountered in many-body expansion of a binding energy. We consider three-body nonadditive energies in molecular and atomic clusters. Second, we address how the accuracy of RPA depends on input provided by different DFT models, without resorting to self-consistent RPA procedure, which is currently impractical for calculations employing periodic boundary conditions. We find that RPA based on the SCAN0 and PBE0 models, that is, hybrid DFT, achieves an overall accuracy between CCSD and MP3 on a data set of molecular trimers from Rezac et al. ( J. Chem. Theory. Comput. 2015 , 11 , 3065 ). Finally, many-body expansion for molecular clusters and solids often leads to a large number of small contributions that need to be calculated with high precision. We therefore present a cubic-scaling (or self-consistent field (SCF)-like) implementation of RPA in atomic basis set, which is designed for calculations with high numerical precision.
- 15Harl, J.; Kresse, G. Accurate Bulk Properties from Approximate Many-Body Techniques. Phys. Rev. Lett. 2009, 103, 056401, DOI: 10.1103/PhysRevLett.103.056401Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXptl2jsbc%253D&md5=1a9a747d56ca24b35808a9aed92bb4dbAccurate Bulk Properties from Approximate Many-Body TechniquesHarl, Judith; Kresse, GeorgPhysical Review Letters (2009), 103 (5), 056401/1-056401/4CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)For ab initio electronic structure calcns., the RPA to the correlation energy is supposed to be a suitable complement to the exact exchange energy. We show that lattice consts., atomization energies of solids, and adsorption energies on metal surfaces evaluated using this approxn. are in very good agreement with expt. Since the method is fairly efficient and handles ionic, metallic, and van der Waals bonded systems equally well, it is a very promising choice to improve upon d. functional theory calcns., without resorting to more demanding diffusion Monte Carlo or quantum chem. methods.
- 16Kreppel, A.; Graf, D.; Laqua, H.; Ochsenfeld, C. Range-Separated Density-Functional Theory in Combination with the Random Phase Approximation: An Accuracy Benchmark. J. Chem. Theory Comput. 2020, 16, 2985– 2994, DOI: 10.1021/acs.jctc.9b01294Google Scholar16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXnslCit7c%253D&md5=e14a638f541b4fb9a2541961f2062ccbRange-Separated Density-Functional Theory in Combination with the Random Phase Approximation: An Accuracy BenchmarkKreppel, Andrea; Graf, Daniel; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2020), 16 (5), 2985-2994CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A formulation of range-sepd. RPA based on our efficient ω-CDGD-RI-RPA [J. Chem. Theory Comput.2018, 14, 2505] method and a large scale benchmark study are presented. By application to the GMTKN55 data set, we obtain a comprehensive picture of the performance of range-sepd. RPA in general main group thermochem., kinetics, and noncovalent interactions. The results show that range-sepd. RPA performs stably over the broad range of mol. chem. included in the GMTKN55 set. It improves significantly over semilocal DFT but it is still less accurate than modern dispersion cor. double-hybrid functionals. Furthermore, range-sepd. RPA shows a faster basis set convergence compared to std. full-range RPA making it a promising applicable approach with only one empirical parameter.
- 17Trushin, E.; Thierbach, A.; Görling, A. Toward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energy. J. Chem. Phys. 2021, 154, 014104, DOI: 10.1063/5.0026849Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXosFWrtg%253D%253D&md5=d3d9596431d9804901d50736957d01abToward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energyTrushin, Egor; Thierbach, Adrian; Goerling, AndreasJournal of Chemical Physics (2021), 154 (1), 014104CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We introduce new functionals for the Kohn-Sham correlation energy that are based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem and are named σ-functionals. Like in the well-established direct RPA (dRPA), σ-functionals require as input exclusively eigenvalues σ of the frequency-dependent KS response function. In the new functionals, functions of σ replace the σ-dependent dRPA expression in the coupling-const. and frequency integrations contained in the ACFD theorem. We optimize σ-functionals with the help of ref. sets for atomization, reaction, transition state, and non-covalent interaction energies. The optimized functionals are to be used in a post-self-consistent way using orbitals and eigenvalues from conventional Kohn-Sham calcns. employing the exchange-correlation functional of Perdew, Burke, and Ernzerhof. The accuracy of the presented approach is much higher than that of dRPA methods and is comparable to that of high-level wave function methods. Reaction and transition state energies from σ-functionals exhibit accuracies close to 1 kcal/mol and thus approach chem. accuracy. For the 10 966 reactions of the W4-11RE ref. set, the mean abs. deviation is 1.25 kcal/mol compared to 3.21 kcal/mol in the dRPA case. Non-covalent binding energies are accurate to a few tenths of a kcal/mol. The presented approach is highly efficient, and the post-self-consistent calcn. of the total energy requires less computational time than a d.-functional calcn. with a hybrid functional and thus can be easily carried out routinely. σ-Functionals can be implemented in any existing dRPA code with negligible programming effort. (c) 2021 American Institute of Physics.
- 18Fauser, S.; Trushin, E.; Neiss, C.; Görling, A. Chemical accuracy with σ-functionals for the Kohn–Sham correlation energy optimized for different input orbitals and eigenvalues. J. Chem. Phys. 2021, 155, 134111, DOI: 10.1063/5.0059641Google Scholar18https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXit1Sqtr3N&md5=7c717988178d285b3f1c00b20c750598Chemical accuracy with σ-functionals for the Kohn-Sham correlation energy optimized for different input orbitals and eigenvaluesFauser, Steffen; Trushin, Egor; Neiss, Christian; Goerling, AndreasJournal of Chemical Physics (2021), 155 (13), 134111CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recently, a new type of orbital-dependent functional for the Kohn-Sham (KS) correlation energy, σ-functionals, was introduced. Tech., σ-functionals are closely related to the well-known direct RPA (dRPA). Within the dRPA, a function of the eigenvalues σ of the frequency-dependent KS response function is integrated over purely imaginary frequencies. In σ-functionals, this function is replaced by one that is optimized with respect to ref. sets of atomization, reaction, transition state, and non-covalent interaction energies. The previously introduced σ-functional uses input orbitals and eigenvalues from KS calcns. with the generalized gradient approxn. (GGA) exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). Here, σ-functionals using input orbitals and eigenvalues from the meta-GGA TPSS and the hybrid-functionals PBE0 and B3LYP are presented and tested. The no. of ref. sets taken into account in the optimization of the σ-functionals is larger than in the first PBE based σ-functional and includes sets with 3d-transition metal compds. Therefore, also a reparameterized PBE based σ-functional is introduced. The σ-functionals based on PBE0 and B3LYP orbitals and eigenvalues reach chem. accuracy for main group chem. For the 10 966 reactions from the highly accurate W4-11RE ref. set, the B3LYP based σ-functional exhibits a mean av. deviation of 1.03 kcal/mol compared to 1.08 kcal/mol for the coupled cluster singles doubles perturbative triples method if the same valence quadruple zeta basis set is used. For 3d-transition metal chem., accuracies of about 2 kcal/mol are reached. The computational effort for the post-self-consistent evaluation of the σ-functional is lower than that of a preceding PBE0 or B3LYP calcn. for typical systems. (c) 2021 American Institute of Physics.
- 19Fauser, S.; Förster, A.; Redeker, L.; Neiss, C.; Erhard, J.; Trushin, E.; Görling, A. Basis Set Requirements of σ-Functionals for Gaussian- and Slater-Type Basis Functions and Comparison with Range-Separated Hybrid and Double Hybrid Functionals. J. Chem. Theory Comput. 2024, 20, 2404– 2422, DOI: 10.1021/acs.jctc.3c01132Google ScholarThere is no corresponding record for this reference.
- 20Drontschenko, V.; Graf, D.; Laqua, H.; Ochsenfeld, C. Efficient Method for the Computation of Frozen-Core Nuclear Gradients within the Random Phase Approximation. J. Chem. Theory Comput. 2022, 18, 7359– 7372, DOI: 10.1021/acs.jctc.2c00774Google Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB38Xisl2gtrzP&md5=01af0b01ddae622b53bd3164b285a9d9Efficient Method for the Computation of Frozen-Core Nuclear Gradients within the Random Phase ApproximationDrontschenko, Viktoria; Graf, Daniel; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2022), 18 (12), 7359-7372CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A method for the evaluation of anal. frozen-core gradients within the RPA is presented. We outline an efficient way to evaluate the response of the d. of active electrons arising only when introducing the frozen-core approxn. and constituting the main difficulty, together with the response of the std. Kohn-Sham d. The general framework allows to extend the outlined procedure to related electron correlation methods in the AO basis that require the evaluation of d. responses, such as second-order Moller-Plesset perturbation theory or coupled cluster variants. By using Cholesky decompd. densities - which reintroduce the occupied index in the time-detg. steps - we are able to achieve speedups of 20-30% (depending on the size of the basis set) by using the frozen-core approxn., which is of similar magnitude as for MO formulations. We further show that the errors introduced by the frozen-core approxn. are practically insignificant for mol. geometries.
- 21Neiss, C.; Fauser, S.; Görling, A. Geometries and vibrational frequencies with Kohn–Sham methods using σ-functionals for the correlation energy. J. Chem. Phys. 2023, 158, 044107, DOI: 10.1063/5.0129524Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3sXhvVGqsL4%253D&md5=f80b8233816be1e4e3ee4925fc0a5357Geometries and vibrational frequencies with Kohn-Sham methods using σ-functionals for the correlation energyNeiss, Christian; Fauser, Steffen; Goerling, AndreasJournal of Chemical Physics (2023), 158 (4), 044107CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recently, Kohn-Sham (KS) methods with new correlation functionals, called σ-functionals, have been introduced. Tech., σ-functionals are closely related to the well-known RPA; formally, σ-functionals are rooted in perturbation theory along the adiabatic connection. If employed in a post-SCF manner in a Gaussian basis set framework, then, σ-functional methods are computationally very efficient. Moreover, for main group chem., σ-functionals are highly accurate and can compete with high-level wave-function methods. For reaction and transition state energies, e.g., chem. accuracy of 1 kcal/mol is reached. Here, we show how to calc. first derivs. of the total energy with respect to nuclear coordinates for methods using σ-functionals and then carry out geometry optimizations for test sets of main group mols., transition metal compds., and non-covalently bonded systems. For main group mols., we addnl. calc. vibrational frequencies. σ-Functional methods are found to yield very accurate geometries and vibrational frequencies for main group mols. superior not only to those from conventional KS methods but also to those from RPA methods. For geometries of transition metal compds., not surprisingly, best geometries are found for RPA methods, while σ-functional methods yield somewhat less good results. This is attributed to the fact that in the optimization of σ-functionals, transition metal compds. could not be represented well due to the lack of reliable ref. data. For non-covalently bonded systems, σ-functionals yield geometries of the same quality as the RPA or as conventional KS schemes combined with dispersion corrections. (c) 2023 American Institute of Physics.
- 22Drontschenko, V.; Bangerter, F. H.; Ochsenfeld, C. Analytical Second-Order Properties for the Random Phase Approximation: Nuclear Magnetic Resonance Shieldings. J. Chem. Theory Comput. 2023, 19, 7542– 7554, DOI: 10.1021/acs.jctc.3c00542Google ScholarThere is no corresponding record for this reference.
- 23Fauser, S.; Drontschenko, V.; Ochsenfeld, C.; Görling, A. Accurate NMR Shieldings with σ-Functionals. J. Chem. Theory Comput. 2024, 20, 6028– 6036, DOI: 10.1021/acs.jctc.4c00512Google ScholarThere is no corresponding record for this reference.
- 24Drontschenko, V.; Ochsenfeld, C. Low-Scaling, Efficient and Memory Optimized Computation of Nuclear Magnetic Resonance Shieldings within the Random Phase Approximation using Cholesky-Decomposed Densities and an Attenuated Coulomb Metric. J. Phys. Chem. A 2024, 128, 7950– 7965, DOI: 10.1021/acs.jpca.4c02773Google ScholarThere is no corresponding record for this reference.
- 25Glasbrenner, M.; Graf, D.; Ochsenfeld, C. Benchmarking the Accuracy of the Direct Random Phase Approximation and σ-Functionals for NMR Shieldings. J. Chem. Theory Comput. 2022, 18, 192– 205, DOI: 10.1021/acs.jctc.1c00866Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXislWrtbvF&md5=25dffbe9eb6cff42e76bdf9a5d151a2eBenchmarking the Accuracy of the Direct Random Phase Approximation and σ-Functionals for NMR ShieldingsGlasbrenner, Michael; Graf, Daniel; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2022), 18 (1), 192-205CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A method for computing NMR shieldings with the direct RPA (RPA) and the closely related σ-functionals [Trushin, E.; Thierbach, A.; Gorling, A. Toward chem. accuracy at low computational cost: d. functional theory with σ-functionals for the correlation energy. J. Chem. Phys.2021,154, 014104] is presented, which is based on a finite-difference approach. The accuracy is evaluated in benchmark calcns. using high-quality coupled cluster values as a ref. Our results show that the accuracy of the computed NMR shieldings using direct RPA is strongly dependent on the d. functional theory ref. orbitals and improves with increasing amts. of exact Hartree-Fock exchange in the functional. NMR shieldings computed with the direct RPA with a Hartree-Fock ref. are significantly more accurate than MP2 shieldings and comparable to CCSD shieldings. Also, the basis set convergence is analyzed and it is shown that at least triple-zeta basis sets are required for reliable results.
- 26Görling, A. Hierarchies of methods towards the exact Kohn–Sham correlation energy based on the adiabatic-connection fluctuation-dissipation theorem. Phys. Rev. B 2019, 99, 235120, DOI: 10.1103/PhysRevB.99.235120Google ScholarThere is no corresponding record for this reference.
- 27Erhard, J.; Bleiziffer, P.; Görling, A. Power Series Approximation for the Correlation Kernel Leading to Kohn–Sham Methods Combining Accuracy, Computational Efficiency, and General Applicability. Phys. Rev. Lett. 2016, 117, 143002, DOI: 10.1103/PhysRevLett.117.143002Google ScholarThere is no corresponding record for this reference.
- 28Erhard, J.; Fauser, S.; Trushin, E.; Görling, A. Scaled σ-functionals for the Kohn–Sham correlation energy with scaling functions from the homogeneous electron gas. J. Chem. Phys. 2022, 157, 114105, DOI: 10.1063/5.0101641Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB38XisVWnsb%252FJ&md5=c3a5d9d86b9d5d799fbd3366b86f21d2Scaled σ-functionals for Kohn-Sham correlation energy with scaling function from homogeneous electron gasErhard, Jannis; Fauser, Steffen; Trushin, Egor; Goerling, AndreasJournal of Chemical Physics (2022), 157 (11), 114105CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. The σ-functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct RPA (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigen values σ of the Kohn-Sham response matrix that enters the coupling const. and frequency integration in the dRPA is replaced by another function optimized with the help of ref. sets of atomization, reaction, transition state, and non-covalent interaction energies, and σ-Functionals are highly accurate and yield chem. accuracy of 1 kcal/mol in reaction or transition state energies, in main group chem. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron no., such as ionization or electron attachment. This problem is attributed to unphys. self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are detected with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addn. yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts. (c) 2022 American Institute of Physics.
- 29Lemke, Y.; Ochsenfeld, C. Highly accurate σ- and τ-functionals for beyond-RPA methods with approximate exchange kernels. J. Chem. Phys. 2023, 159, 194104, DOI: 10.1063/5.0173042Google Scholar29https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3sXitlOgtr3N&md5=7cb6539f08f4d9e17b24a18fb1bfda82Highly accurate σ- and τ-functionals for beyond-RPA methods with approximate exchange kernelsLemke, Yannick; Ochsenfeld, ChristianJournal of Chemical Physics (2023), 159 (19), 194104CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)σ-Functionals are promising new developments for the Kohn-Sham correlation energy based upon the direct RPA (dRPA) within the adiabatic connection formalism, providing impressive improvements over dRPA for a broad range of benchmarks. However, σ-functionals exhibit a high amt. of self-interaction inherited from the approxns. made within dRPA. Inclusion of an exchange kernel in deriving the coupling-strength-dependent d.-d. response function leads to so-called τ-functionals, which - apart from a fourth-order Taylor series expansion - have only been realized in an approx. fashion so far to the best of our knowledge, most notably in the form of scaled σ-functionals. In this work, we derive, optimize, and benchmark three types of σ- and τ-functionals including approx. exchange effects in the form of an antisymmetrized Hartree kernel. These functionals, based on a second-order screened exchange type contribution in the adiabatic connection formalism, the electron-hole time-dependent Hartree-Fock kernel (eh-TDHF) otherwise known as RPA with exchange (RPAx), and an approxn. thereof known as approx. exchange kernel (AXK), are optimized on the ASCDB database using two new parametrizations named A1 and A2. In addn., we report a first full evaluation of σ- and τ-functionals on the GMTKN55 database, revealing our exchange-including functionals to considerably outperform existing σ-functionals while being highly competitive with some of the best double-hybrid functionals of the original GMTKN55 publication. In particular, the σ-functionals based on AXK and τ-functionals based on RPAx with PBE0 ref. stand out as highly accurate approaches for a wide variety of chem. relevant problems. (c) 2023 American Institute of Physics.
- 30Grüneis, A.; Marsman, M.; Harl, J.; Schimka, L.; Kresse, G. Making the random phase approximation to electronic correlation accurate. J. Chem. Phys. 2009, 131, 154115, DOI: 10.1063/1.3250347Google Scholar30https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhsVSjs7fJ&md5=d186251717a9fe10e70b360c9f350c53Making the random phase approximation to electronic correlation accurateGrueneis, Andreas; Marsman, Martijn; Harl, Judith; Schimka, Laurids; Kresse, GeorgJournal of Chemical Physics (2009), 131 (15), 154115/1-154115/5CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We show that the inclusion of second-order screened exchange to the RPA allows for an accurate description of electronic correlation in atoms and solids clearly surpassing the random phase approxn., but not yet approaching chem. accuracy. From a fundamental point of view, the method is self-correlation free for one-electron systems. From a practical point of view,the approach yields correlation energies for atoms, as well as for the jellium electron gas within a few kcal/mol of exact values, atomization energies within typically 2-3 kcal/mol of expt., and excellent lattice consts. for ionic and covalently bonded solids (0.2% error). The computational complexity is only O(N5), comparable to canonical second-order Moller-Plesset perturbation theory, which should allow for routine calcns. on many systems. (c) 2009 American Institute of Physics.
- 31Bates, J. E.; Furche, F. Communication: Random phase approximation renormalized many-body perturbation theory. J. Chem. Phys. 2013, 139, 171103, DOI: 10.1063/1.4827254Google Scholar31https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhsleiu73P&md5=bbfd81d02296e4acbc349e365c8b9b70Communication: Random phase approximation renormalized many-body perturbation theoryBates, Jefferson E.; Furche, FilippJournal of Chemical Physics (2013), 139 (17), 171103/1-171103/4CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We derive a renormalized many-body perturbation theory (MBPT) starting from the RPA. This RPA-renormalized perturbation theory extends the scope of single-ref. MBPT methods to small-gap systems without significantly increasing the computational cost. The leading correction to RPA, termed the approx. exchange kernel (AXK), substantially improves upon RPA atomization energies and ionization potentials without affecting other properties such as barrier heights where RPA is already accurate. Thus, AXK is more balanced than second-order screened exchange, which tends to overcorrect RPA for systems with stronger static correlation. Similarly, AXK avoids the divergence of second-order Moller-Plesset (MP2) theory for small gap systems and delivers a much more consistent performance than MP2 across the periodic table at comparable cost. RPA+AXK thus is an accurate, non-empirical, and robust tool to assess and improve semi-local d. functional theory for a wide range of systems previously inaccessible to first-principles electronic structure calcns. (c) 2013 American Institute of Physics.
- 32Goerigk, L.; Hansen, A.; Bauer, C.; Ehrlich, S.; Najibi, A.; Grimme, S. A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions. Phys. Chem. Chem. Phys. 2017, 19, 32184, DOI: 10.1039/C7CP04913GGoogle Scholar32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhslajtLnF&md5=f9393c9e3907336c4da053743797f8dfA look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactionsGoerigk, Lars; Hansen, Andreas; Bauer, Christoph; Ehrlich, Stephan; Najibi, Asim; Grimme, StefanPhysical Chemistry Chemical Physics (2017), 19 (48), 32184-32215CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We present the GMTKN55 benchmark database for general main group thermochem., kinetics and noncovalent interactions. Compared to its popular predecessor GMTKN30, it allows assessment across a larger variety of chem. problems - with 13 new benchmark sets being presented for the first time - and it also provides ref. values of significantly higher quality for most sets. GMTKN55 comprises 1505 relative energies based on 2462 single-point calcns. and it is accessible to the user community via a dedicated website. Herein, we demonstrate the importance of better ref. values, and we re-emphasize the need for London-dispersion corrections in d. functional theory (DFT) treatments of thermochem. problems, including Minnesota methods. We assessed 217 variations of dispersion-cor. and -uncorrected d. functional approxns., and carried out a detailed anal. of 83 of them to identify robust and reliable approaches. Double-hybrid functionals are the most reliable approaches for thermochem. and noncovalent interactions, and they should be used whenever tech. feasible. These are, in particular, DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), and B2GPPLYP-D3(BJ). The best hybrids are ωB97X-V, M052X-D3(0), and ωB97X-D3, but we also recommend PW6B95-D3(BJ) as the best conventional global hybrid. At the meta-generalized-gradient (meta-GGA) level, the SCAN-D3(BJ) method can be recommended. Other meta-GGAs are outperformed by the GGA functionals revPBE-D3(BJ), B97-D3(BJ), and OLYP-D3(BJ). We note that many popular methods, such as B3LYP, are not part of our recommendations. In fact, with our results we hope to inspire a change in the user community's perception of common DFT methods. We also encourage method developers to use GMTKN55 for cross-validation studies of new methodologies.
- 33Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865– 3868, DOI: 10.1103/PhysRevLett.77.3865Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XmsVCgsbs%253D&md5=55943538406ee74f93aabdf882cd4630Generalized gradient approximation made simplePerdew, John P.; Burke, Kieron; Ernzerhof, MatthiasPhysical Review Letters (1996), 77 (18), 3865-3868CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)Generalized gradient approxns. (GGA's) for the exchange-correlation energy improve upon the local spin d. (LSD) description of atoms, mols., and solids. We present a simple derivation of a simple GGA, in which all parameters (other than those in LSD) are fundamental consts. Only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked. Improvements over PW91 include an accurate description of the linear response of the uniform electron gas, correct behavior under uniform scaling, and a smoother potential.
- 34Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158– 6170, DOI: 10.1063/1.478522Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXitVCmt7Y%253D&md5=cad4185c69f9232753497f5203d6dc9fToward reliable density functional methods without adjustable parameters: the PBE0 modelAdamo, Carlo; Barone, VincenzoJournal of Chemical Physics (1999), 110 (13), 6158-6170CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present an anal. of the performances of a parameter free d. functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amt. of exact exchange. The results obtained for structural, thermodn., kinetic and spectroscopic (magnetic, IR and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chem. and condensed matter physics.
- 35Ernzerhof, M.; Scuseria, G. E. Assessment of the Perdew–Burke–Ernzerhof exchange-correlation functional. J. Chem. Phys. 1999, 110, 5029– 5036, DOI: 10.1063/1.478401Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXhsFWquro%253D&md5=6113a915baa66b0f9961a5e160cadaf1Assessment of the Perdew-Burke-Ernzerhof exchange-correlation functionalErnzerhof, Matthias; Scuseria, Gustavo E.Journal of Chemical Physics (1999), 110 (11), 5029-5036CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)In order to discriminate between approxns. to the exchange-correlation energy EXC[ρ↑,ρ↓], we employ the criterion of whether the functional is fitted to a certain exptl. data set or if it is constructed to satisfy phys. constraints. We present extensive test calcns. for atoms and mols., with the nonempirical local spin-d. (LSD) and the Perdew-Burke-Ernzerhof (PBE) functional and compare our results with results obtained with more empirical functionals. For the atomization energies of the G2 set, we find that the PBE functional shows systematic errors larger than those of commonly used empirical functionals. The PBE ionization potentials, electron affinities, and bond lengths are of accuracy similar to those obtained from empirical functionals. Furthermore, a recently proposed hybrid scheme using exact exchange together with PBE exchange and correlation is investigated. For all properties studied here, the PBE hybrid gives an accuracy comparable to the frequently used empirical B3LYP hybrid scheme. Phys. principles underlying the PBE and PBE hybrid scheme are examd. and the range of their validity is discussed.
- 36Langreth, D. C.; Perdew, J. P. The exchange-correlation energy of a metallic surface. Solid State Commun. 1975, 17, 1425, DOI: 10.1016/0038-1098(75)90618-3Google ScholarThere is no corresponding record for this reference.
- 37Langreth, D. C.; Perdew, J. P. Exchange-correlation energy of a metallic surface: Wave-vector analysis. Phys. Rev. B 1977, 15, 2884, DOI: 10.1103/PhysRevB.15.2884Google ScholarThere is no corresponding record for this reference.
- 38Petersilka, M.; Gossmann, U. J.; Gross, E. K. U. Excitation Energies from Time-Dependent Density-Functional Theory. Phys. Rev. Lett. 1996, 76, 1212– 1215, DOI: 10.1103/PhysRevLett.76.1212Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XhtFGjt7k%253D&md5=67ce109fced1909cc6dabf2f9aab0116Excitation energies from time-dependent density-functional theoryPetersilka, M.; Gossmann, U. J.; Gross, E. K. U.Physical Review Letters (1996), 76 (8), 1212-15CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)A new d.-functional approach to calc. the excitation spectrum of many-electron systems is proposed. It is shown that the full linear d. response of the interacting system, which has poles at the exact excitation energies, can rigorously be expressed in terms of the response function of the noninteracting (Kohn-Sham) system and a frequency-dependent exchange-correlation kernel. Using this expression, the poles of the full response function are obtained by systematic improvement upon the poles of the Kohn-Sham response function. Numerical results are presented for Be, Mg, Ca, Zn, Sr, and Cd atoms.
- 39Garrick, R.; Natan, A.; Gould, T.; Kronik, L. Exact Generalized Kohn–Sham Theory for Hybrid Functionals. Phys. Rev. X 2020, 10, 021040, DOI: 10.1103/PhysRevX.10.021040Google ScholarThere is no corresponding record for this reference.
- 40Whitten, J. L. Coulombic potential energy integrals and approximations. J. Chem. Phys. 1973, 58, 4496– 4501, DOI: 10.1063/1.1679012Google Scholar40https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE3sXktFSjtbo%253D&md5=05c510c8db660386b6fd3df789c22b10Coulombic potential energy integrals and approximationsWhitten, J. L.Journal of Chemical Physics (1973), 58 (10), 4496-501CODEN: JCPSA6; ISSN:0021-9606.Theorems are derived which establish a method of approxg. 2-particle Coulombic potential energy integrals, [.vphi.a(1)|r12-1|.vphi.b-(2)], in terms of approx. charge ds. .vphi.a' and .vphi.b'. Rigorous error bounds, |[.vphi.a(1)|r12-1|.vphi.b(2)] - [.vphi.a'(1)|r12-1|.vphi.b'(2)]| ≤ δ, are simply expressed in terms of information calcd. sep. for the pair of ds. .vphi.a and .vphi.b' and the pair .vphi.b and .vphi.b'. From the structure of the bound, a simple method of optimizing charge d. approxns. such that δ is minimized is derived. The framework of the theory appears to be well suited for application to the approxn. of electron repulsion integrals which occur in mol. structure theory, and applications to the approxn. of integrals over Slater orbitals or grouped Gaussian functions are discussed.
- 41Dunlap, B. I.; Connolly, J. W. D.; Sabin, J. R. On some approximations in applications of Xα theory. J. Chem. Phys. 1979, 71, 3396– 3402, DOI: 10.1063/1.438728Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE1MXlvF2hs7k%253D&md5=0d61bda3e94a1ef30d221972f51ea6cfOn some approximations in applications of X α theoryDunlap, B. I.; Connolly, J. W. D.; Sabin, J. R.Journal of Chemical Physics (1979), 71 (8), 3396-402CODEN: JCPSA6; ISSN:0021-9606.An approx. Xα functional is proposed from which the charge d. fitting equations follow variationally. LCAO Xα calcns. on at. Ni and H2 show the method independent of the fitting (auxiliary) bases to within 0.02 eV. Variational properties assocd. with both orbital and auxiliary basis set incompleteness are used to approach within 0.2 eV the Xα total energy limit for the N mol.
- 42Vahtras, O.; Almlöf, J.; Feyereisen, M. W. Integral approximations for LCAO-SCF calculations. Chem. Phys. Lett. 1993, 213, 514– 518, DOI: 10.1016/0009-2614(93)89151-7Google Scholar42https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3sXmsFKksL4%253D&md5=2c0f258baf7af4b891974d9b71f3bbb2Integral approximations for LCAO-SCF calculationsVahtras, O.; Almloef, J.; Feyereisen, M. W.Chemical Physics Letters (1993), 213 (5-6), 514-18CODEN: CHPLBC; ISSN:0009-2614.Three-center approxns. to the four-center integrals occurring in ab initio LCAO calcns. are investigated. Significant gains in computer time can be obtained without sacrificing accuracy, if a suitable expansion basis is chosen.
- 43Furche, F. Developing the random phase approximation into a practical post-Kohn–Sham correlation model. J. Chem. Phys. 2008, 129, 114105, DOI: 10.1063/1.2977789Google Scholar43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhtFOrtbbP&md5=b4d01f7faebebbfbbbf6ebc9c0077102Developing the random phase approximation into a practical post-Kohn-Sham correlation modelFurche, FilippJournal of Chemical Physics (2008), 129 (11), 114105/1-114105/8CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The RPA (RPA) to the d. functional correlation energy systematically improves upon many limitations of present semilocal functionals, but was considered too computationally expensive for widespread use in the past. Here a phys. appealing reformulation of the RPA correlation model is developed that substantially reduces its computational complexity. The d. functional RPA correlation energy is shown to equal one-half times the difference of all RPA electronic excitation energies computed at full and first order coupling. Thus, the RPA correlation energy may be considered as a difference of electronic zero point vibrational energies, where each eigenmode corresponds to an electronic excitation. This surprisingly simple result is intimately related to plasma theories of electron correlation. Differences to electron pair correlation models underlying popular correlated wave function methods are discussed. The RPA correlation energy is further transformed into an explicit functional of the Kohn-Sham orbitals. The only nontrivial ingredient to this functional is the sign function of the response operator. A stable iterative algorithm to evaluate this sign function based on the Newton-Schulz iteration is presented. Integral direct implementations scale as the fifth power of the system size, similar to second order Moeller-Plesset calcns. With these improvements, RPA may become the long-sought robust and efficient zero order post-Kohn-Sham correlation model. (c) 2008 American Institute of Physics.
- 44Eshuis, H.; Yarkony, J.; Furche, F. Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration. J. Chem. Phys. 2010, 132, 234114, DOI: 10.1063/1.3442749Google Scholar44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXnslSlsLw%253D&md5=90e538090088f1ded996c402ffc7b342Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integrationEshuis, Henk; Yarkony, Julian; Furche, FilippJournal of Chemical Physics (2010), 132 (23), 234114/1-234114/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The RPA is an increasingly popular post-Kohn-Sham correlation method, but its high computational cost has limited mol. applications to systems with few atoms. Here we present an efficient implementation of RPA correlation energies based on a combination of resoln. of the identity (RI) and imaginary frequency integration techniques. We show that the RI approxn. to four-index electron repulsion integrals leads to a variational upper bound to the exact RPA correlation energy if the Coulomb metric is used. Auxiliary basis sets optimized for second-order Moller-Plesset (MP2) calcns. are well suitable for RPA, as is demonstrated for the HEAT and MOLEKEL benchmark sets. Using imaginary frequency integration rather than diagonalization to compute the matrix square root necessary for RPA, evaluation of the RPA correlation energy requires O(N4logN) operations and O(N3) storage only; the price for this dramatic improvement over existing algorithms is a numerical quadrature. We propose a numerical integration scheme that is exact in the two-orbital case and converges exponentially with the no. of grid points. For most systems, 30-40 grid points yield μH accuracy in triple zeta basis sets, but much larger grids are necessary for small gap systems. The lowest-order approxn. to the present method is a post-Kohn-Sham frequency-domain version of opposite-spin Laplace-transform RI-MP2. Timings for polyacenes with up to 30 atoms show speed-ups of two orders of magnitude over previous implementations. The present approach makes it possible to routinely compute RPA correlation energies of systems well beyond 100 atoms, as is demonstrated for the octapeptide angiotensin II. (c) 2010 American Institute of Physics.
- 45Bleiziffer, P.; Heßelmann, A.; Görling, A. Resolution of identity approach for the Kohn–Sham correlation energy within the exact-exchange random-phase approximation. J. Chem. Phys. 2012, 136, 134102, DOI: 10.1063/1.3697845Google ScholarThere is no corresponding record for this reference.
- 46Lemke, Y.; Graf, D.; Kussmann, J.; Ochsenfeld, C. An assessment of orbital energy corrections for the direct random phase approximation and explicit σ-functionals. Mol. Phys. 2023, 121, e2098862, DOI: 10.1080/00268976.2022.2098862Google ScholarThere is no corresponding record for this reference.
- 47Heßelmann, A.; Görling, A. Random phase approximation correlation energies with exact Kohn–Sham exchange. Mol. Phys. 2010, 108, 359– 372, DOI: 10.1080/00268970903476662Google Scholar47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXjvFKksLw%253D&md5=099fc620e45aaffe5161baea6e7454acRandom phase approximation correlation energies with exact Kohn-Sham exchangeHesselmann, Andreas; Goerling, AndreasMolecular Physics (2010), 108 (3-4), 359-372CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)The RPA correlation energy is expressed in terms of the exact local Kohn-Sham (KS) exchange potential and corresponding adiabatic and nonadiabatic exchange kernels for d.-functional ref. determinants. The approach naturally extends the RPA method in which, conventionally, only the Coulomb kernel is included. By comparison with the coupled cluster singles doubles with perturbative triples method it is shown for a set of small mols. that the new RPA method based on KS exchange yields correlation energies more accurate than RPA on the basis of Hartree-Fock exchange.
- 48Ángyán, J. G.; Liu, R.-F.; Toulouse, J.; Jansen, G. Correlation Energy Expressions from the Adiabatic-Connection Fluctuation–Dissipation Theorem Approach. J. Chem. Theory Comput. 2011, 7, 3116– 3130, DOI: 10.1021/ct200501rGoogle Scholar48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXht1Sltb3P&md5=03529fc836d25b236e71c89e4170f1adCorrelation Energy Expressions from the Adiabatic-Connection Fluctuation-Dissipation Theorem ApproachAngyan, Janos G.; Liu, Ru-Fen; Toulouse, Julien; Jansen, GeorgJournal of Chemical Theory and Computation (2011), 7 (10), 3116-3130CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We explore several RPA correlation energy variants within the adiabatic-connection fluctuation-dissipation theorem approach. These variants differ in the way the exchange interactions are treated. One of these variants, named dRPA-II, is original to this work and closely resembles the second-order screened exchange (SOSEX) method. We discuss and clarify the connections among different RPA formulations. We derive the spin-adapted forms of all the variants for closed-shell systems and test them on a few at. and mol. systems with and without range sepn. of the electron-electron interaction.
- 49Heßelmann, A. Random-phase-approximation correlation method including exchange interactions. Phys. Rev. A 2012, 85, 012517, DOI: 10.1103/PhysRevA.85.012517Google Scholar49https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XjtFahu7g%253D&md5=1b15a12a2c98b16c7db2d9a9ced0d34eRandom-phase-approximation correlation method including exchange interactionsHesselmann, AndreasPhysical Review A: Atomic, Molecular, and Optical Physics (2012), 85 (1-A), 012517/1-012517/10CODEN: PLRAAN; ISSN:1050-2947. (American Physical Society)Two random-phase-approxn. correlation methods are introduced that take into account exchange interactions. The first one, termed RPAX, is obtained from a simple modification of the ring coupled-cluster doubles amplitude equation, while the second, termed RPAX2, is based on the first method using a slightly modified update equation for the amplitudes. It is shown that this second RPAX2 method can be implemented with a computational algorithm that scales only with the fifth power of the mol. size with the aid of d. fitting or the Cholesky decompn. of two-electron integrals. It is thus not much more costly than std. second-order perturbation theory methods and can be applied to quite large mol. systems. Moreover, numerical tests for chem. reaction energies and intermol. interaction energies have shown that the RPAX2 method, if based on a Perdew-Burke-Ernzerhof exchange Kohn-Sham ref. determinant, yields results which are very close to coupled-cluster with single, double, and perturbative triple excitations ref. results.
- 50Eshuis, H.; Bates, J. E.; Furche, F. Electron correlation methods based on the random phase approximation. Theor. Chem. Acc. 2012, 131, 1084, DOI: 10.1007/s00214-011-1084-8Google ScholarThere is no corresponding record for this reference.
- 51Mussard, B.; Rocca, D.; Jansen, G.; Ángyán, J. G. Dielectric Matrix Formulation of Correlation Energies in the Random Phase Approximation: Inclusion of Exchange Effects. J. Chem. Theory Comput. 2016, 12, 2191– 2202, DOI: 10.1021/acs.jctc.5b01129Google Scholar51https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XksVGksrk%253D&md5=6e7bf7138c73d863bdaa4353c728fc81Dielectric Matrix Formulation of Correlation Energies in the Random Phase Approximation: Inclusion of Exchange EffectsMussard, Bastien; Rocca, Dario; Jansen, Georg; Angyan, Janos G.Journal of Chemical Theory and Computation (2016), 12 (5), 2191-2202CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Starting from the general expression for the ground state correlation energy in the adiabatic-connection fluctuation-dissipation theorem (ACFDT) framework, it is shown that the dielec. matrix formulation, which is usually applied to calc. the direct RPA (dRPA) correlation energy, can be used for alternative RPA expressions including exchange effects. Within this famework, the ACFDT analog of the second order screened exchange (SOSEX) approxn. leads to a logarithmic formula for the correlation energy similar to the direct RPA expression. Alternatively, the contribution of the exchange can be included in the kernel used to evaluate the response functions. In this case, the use of an approx. kernel is crucial to simplify the formalism and to obtain a correlation energy in logarithmic form. Tech. details of the implementation of these methods are discussed, and it is shown that one can take advantage of d. fitting or Cholesky decompn. techniques to improve the computational efficiency; a discussion on the numerical quadrature made on the frequency variable is also provided. A series of test calcns. on at. correlation energies and mol. reaction energies shows that exchange effects are instrumental for improvement over direct RPA results.
- 52Dixit, A.; Ángyán, J. G.; Rocca, D. Improving the accuracy of ground-state correlation energies within a plane-wave basis set: The electron-hole exchange kernel. J. Chem. Phys. 2016, 145, 104105, DOI: 10.1063/1.4962352Google Scholar52https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFWnurnF&md5=d7cf5cda8cc53b00dd83016f7089b527Improving the accuracy of ground-state correlation energies within a plane-wave basis set: The electron-hole exchange kernelDixit, Anant; Angyan, Janos G.; Rocca, DarioJournal of Chemical Physics (2016), 145 (10), 104105/1-104105/12CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A new formalism was recently proposed to improve RPA correlation energies by including approx. exchange effects. Within this framework, by keeping only the electron-hole contributions to the exchange kernel, two approxns. can be obtained: An adiabatic connection analog of the second order screened exchange (AC-SOSEX) and an approx. electron-hole time-dependent Hartree-Fock (eh-TDHF). Here we show how this formalism is suitable for an efficient implementation within the plane-wave basis set. The response functions involved in the AC-SOSEX and eh-TDHF equations can indeed be compactly represented by an auxiliary basis set obtained from the diagonalization of an approx. dielec. matrix. Addnl., the explicit calcn. of unoccupied states can be avoided by using d. functional perturbation theory techniques and the matrix elements of dynamical response functions can be efficiently computed by applying the Lanczos algorithm. As shown by several applications to reaction energies and weakly bound dimers, the inclusion of the electron-hole kernel significantly improves the accuracy of ground-state correlation energies with respect to RPA and semi-local functionals. (c) 2016 American Institute of Physics.
- 53Chen, G. P.; Agee, M. M.; Furche, F. Performance and Scope of Perturbative Corrections to Random-Phase Approximation Energies. J. Chem. Theory Comput. 2018, 14, 5701– 5714, DOI: 10.1021/acs.jctc.8b00777Google Scholar53https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhslOlsrfM&md5=15027735d3a819dacfdeceef8ba62b99Performance and Scope of Perturbative Corrections to Random-Phase Approximation EnergiesChen, Guo P.; Agee, Matthew M.; Furche, FilippJournal of Chemical Theory and Computation (2018), 14 (11), 5701-5714CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)It has been suspected since the early days of the RPA that corrections to RPA correlation energies result mostly from short-range correlation effects and are thus amenable to perturbation theory. Here we test this hypothesis by analyzing formal and numerical results for the most common beyond-RPA perturbative corrections, including the bare second-order exchange (SOX), second-order screened exchange (SOSEX), and approx. exchange kernel (AXK) methods. Our anal. is facilitated by efficient and robust algorithms based on the resoln.-of-the-identity (RI) approxn. and numerical frequency integration, which enable benchmark beyond-RPA calcns. on medium- and large-size mols. with size-independent accuracy. The AXK method systematically improves upon RPA, SOX, and SOSEX for reaction barrier heights, reaction energies, and noncovalent interaction energies of main-group compds. The improved accuracy of AXK compared with SOX and SOSEX is attributed to stronger screening of bare SOX in AXK. For reactions involving transition-metal compds., particularly 3d transition-metal dimers, the AXK correction is too small and can even have the wrong sign. These observations are rationalized by a measure ‾α of the effective coupling strength for beyond-RPA correlation. When the effective coupling strength increases beyond a crit. ‾α value of approx. 0.5, the RPA errors increase rapidly and perturbative corrections become unreliable. Thus, perturbation theory can systematically correct RPA but only for systems and properties qual. well captured by RPA, as indicated by small ‾α values.
- 54Kussmann, J.; Ochsenfeld, C. Pre-selective screening for matrix elements in linear-scaling exact exchange calculations. J. Chem. Phys. 2013, 138, 134114, DOI: 10.1063/1.4796441Google Scholar54https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXlt1Glur0%253D&md5=ed3dfce7bf9beed628dae80664085777Pre-selective screening for matrix elements in linear-scaling exact exchange calculationsKussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Physics (2013), 138 (13), 134114/1-134114/7CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a simple but accurate preselection method based on Schwarz integral ests. to det. the significant elements of the exact exchange matrix before its evaluation, thus providing an asymptotical linear-scaling behavior for non-metallic systems. Our screening procedure proves to be highly suitable for exchange matrix calcns. on massively parallel computing architectures, such as graphical processing units, for which we present a first linear-scaling exchange matrix evaluation algorithm. (c) 2013 American Institute of Physics.
- 55Kussmann, J.; Ochsenfeld, C. Preselective Screening for Linear-Scaling Exact Exchange-Gradient Calculations for Graphics Processing Units and General Strong-Scaling Massively Parallel Calculations. J. Chem. Theory Comput. 2015, 11, 918, DOI: 10.1021/ct501189uGoogle Scholar55https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXjsFCnuro%253D&md5=b323b53aca62835e043f8ed31503c238Preselective Screening for Linear-Scaling Exact Exchange-Gradient Calculations for Graphics Processing Units and General Strong-Scaling Massively Parallel CalculationsKussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2015), 11 (3), 918-922CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present an extension of our recently presented PreLinK scheme (J. Chem. Phys. 2013, 138, 134114) for the exact exchange contribution to nuclear forces. The significant contributions to the exchange gradient are detd. by preselection based on accurate shell-pair contributions to the SCF exchange energy prior to the calcn. Therefore, our method is highly suitable for massively parallel electronic structure calcns. because of an efficient load balancing of the significant contributions only and an unhampered control flow. The efficiency of our method is shown for several illustrative calcns. on single GPU servers, as well as for hybrid MPI/CUDA parallel calcns. with the largest system comprising 3369 atoms and 26952 basis functions.
- 56Kussmann, J.; Laqua, H.; Ochsenfeld, C. Highly Efficient Resolution-of-Identity Density Functional Theory Calculations on Central and Graphics Processing Units. J. Chem. Theory Comput. 2021, 17, 1512, DOI: 10.1021/acs.jctc.0c01252Google Scholar56https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXksFWkurs%253D&md5=36c4001fa1dbe4ef9fdaeefa114cd851Highly Efficient Resolution-of-Identity Density Functional Theory Calculations on Central and Graphics Processing UnitsKussmann, Joerg; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2021), 17 (3), 1512-1521CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present an efficient method to evaluate Coulomb potential matrixes using the resoln. of identity approxn. and semilocal exchange-correlation potentials on central (CPU) and graphics processing units (GPU). The new GPU-based RI-algorithm shows a high performance and ensures the favorable scaling with increasing basis set size as the conventional CPU-based method. Furthermore, our method is based on the J-engine algorithm [White, Head-Gordon, J. Chem. Phys., 1996, 7, 2620], which allows for further optimizations that also provide a significant improvement of the corresponding CPU-based algorithm. Due to the increased performance for the Coulomb evaluation, the calcn. of the exchange-correlation potential of d. functional theory on CPUs quickly becomes a bottleneck to the overall computational time. Hence, we also present a GPU-based algorithm to evaluate the exchange-correlation terms, which results in an overall high-performance method for d. functional calcns. The algorithms to evaluate the potential and nuclear deriv. terms are discussed, and their performance on CPUs and GPUs is demonstrated for illustrative calcns.
- 57Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057, DOI: 10.1039/b515623hGoogle Scholar57https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xhs12ntrc%253D&md5=314690393f1e21096541a317a80e563cAccurate Coulomb-fitting basis sets for H to RnWeigend, FlorianPhysical Chemistry Chemical Physics (2006), 8 (9), 1057-1065CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)A series of auxiliary basis sets to fit Coulomb potentials for the elements H to Rn (except lanthanides) is presented. For each element only one auxiliary basis set is needed to approx. Coulomb energies in conjunction with orbital basis sets of split valence, triple zeta valence and quadruple zeta valence quality with errors of typically below ca. 0.15 kJ mol-1 per atom; this was demonstrated in conjunction with the recently developed orbital basis sets of types def2-SV(P), def2-TZVP and def2-QZVPP for a large set of small mols. representing (nearly) each element in all of its common oxidn. states. These auxiliary bases are slightly more than three times larger than orbital bases of split valence quality. Compared to non-approximated treatments, computation times for the Coulomb part are reduced by a factor of ca. 8 for def2-SV(P) orbital bases, ca. 25 for def2-TZVP and ca. 100 for def2-QZVPP orbital bases.
- 58Laqua, H.; Kussmann, J.; Ochsenfeld, C. Efficient and Linear-Scaling Seminumerical Method for Local Hybrid Density Functionals. J. Chem. Theory Comput. 2018, 14, 3451– 3458, DOI: 10.1021/acs.jctc.8b00062Google Scholar58https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtFarurvK&md5=a22f0d7c57ca427b744970b0a450034fEfficient and Linear-Scaling Seminumerical Method for Local Hybrid Density FunctionalsLaqua, Henryk; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2018), 14 (7), 3451-3458CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Local hybrid functionals, i.e., functionals with local dependence on the exact exchange energy d., generalize the popular class of global hybrid functionals and extend the applicability of d. functional theory to electronic structures that require an accurate description of static correlation. However, the higher computational cost compared to conventional Kohn-Sham d. functional theory restrained their widespread application. Here, we present a low-prefactor, linear-scaling method to evaluate the local hybrid exchange-correlation potential as well as the corresponding nuclear forces by employing a seminumerical integration scheme. In the seminumerical scheme, one integration in the electron repulsion integrals is carried out anal. and the other one is carried out numerically, employing an integration grid. A high computational efficiency is achieved by combining the preLinK method with explicit screening of integrals for batches of grid points to minimize the screening overhead. This new method, denoted as preLinX, provides an 8-fold performance increase for a DNA fragment contg. four base pairs as compared to existing S- and P-junction-based methods. In this way, our method allows for the evaluation of local hybrid functionals at a cost similar to that of global hybrid functionals. The linear-scaling behavior, efficiency, accuracy, and multi-node parallelization of the presented method is demonstrated for large systems contg. more than 1000 atoms.
- 59Laqua, H.; Thompson, T. H.; Kussmann, J.; Ochsenfeld, C. Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing Units. J. Chem. Theory Comput. 2020, 16, 1456, DOI: 10.1021/acs.jctc.9b00860Google Scholar59https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXivVejs7Y%253D&md5=606be25e20a70cc545ff71ebabf28152Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing UnitsLaqua, Henryk; Thompson, Travis H.; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2020), 16 (3), 1456-1468CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present a highly efficient and asymptotically linear-scaling graphic processing unit accelerated seminumerical exact-exchange method (sn-LinK). We go beyond our previous central processing unit-based method (H. Laqua et al., 2018) by employing our recently developed integral bounds (T.H. Thomson and C. Ochsenfeld, 2019) and high-accuracy numerical integration grid (H. Laqua et al., 2018). The accuracy is assessed for several established test sets, providing errors significantly below 1mEh for the smallest grid. Moreover, a comprehensive performance anal. for large mols. between 62 and 1347 atoms is provided, revealing the outstanding performance of our method, in particular, for large basis sets such as the polarized quadruple-zeta level with diffuse functions.
- 60Laqua, H.; Kussmann, J.; Ochsenfeld, C. An improved molecular partitioning scheme for numerical quadratures in density functional theory. J. Chem. Phys. 2018, 149, 204111, DOI: 10.1063/1.5049435Google Scholar60https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXitlylurvI&md5=9cc4db762351543f44c05b2c28868ccbAn improved molecular partitioning scheme for numerical quadratures in density functional theoryLaqua, Henryk; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Physics (2018), 149 (20), 204111/1-204111/8CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a modification to Becke's mol. partitioning scheme [A. D. Becke, J. Chem. Phys. 88, 2547 (1988)] that provides substantially better accuracy for weakly bound complexes and allows for a faster and linear scaling grid generation without introducing a cutoff error. We present the accuracy of our new partitioning scheme for atomization energies of small mols. and for interaction energies of van der Waals complexes. Furthermore, the efficiency and scaling behavior of the grid generation are demonstrated for large mol. systems with up to 1707 atoms. (c) 2018 American Institute of Physics.
- 61Lehtola, S.; Steigemann, C.; Oliveira, M. J. T.; Marques, M. A. L. Recent developments in libxc–A comprehensive library of functionals for density functional theory. Software X 2018, 7, 1– 5, DOI: 10.1016/j.softx.2017.11.002Google ScholarThere is no corresponding record for this reference.
- 62Weigend, F.; Furche, F.; Ahlrichs, R. Gaussian basis sets of quadruple zeta valence quality for atoms H–Kr. J. Chem. Phys. 2003, 119, 12753– 12762, DOI: 10.1063/1.1627293Google Scholar62https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXpvFOns7s%253D&md5=9d1664e0004c3a900b7efcb5da10dd6dGaussian basis sets of quadruple zeta valence quality for atoms H-KrWeigend, Florian; Furche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2003), 119 (24), 12753-12762CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present Gaussian basis sets of quadruple zeta valence quality with a segmented contraction scheme for atoms H to Kr. This extends earlier work on segmented contracted split valence (SV) and triple zeta valence (TZV) basis sets. Contraction coeffs. and orbital exponents are fully optimized in at. Hartree-Fock (HF) calcns. As opposed to other quadruple zeta basis sets, the basis set errors in at. ground-state HF energies are less than 1 mEh and increase smoothly across the Periodic Table, while the no. of primitives is comparably small. Polarization functions are taken partly from previous work, partly optimized in at. MP2 calcns., and for a few cases detd. at the HF level for excited at. states nearly degenerate with the ground state. This leads to basis sets denoted QZVP for HF and d. functional theory (DFT) calcns., and for some atoms to a larger basis recommended for correlated treatments, QZVPP. We assess the performance of the basis sets in mol. HF, DFT, and MP2 calcns. for a sample of diat. and small polyat. mols. by a comparison of energies, bond lengths, and dipole moments with results obtained numerically or using very large basis sets. It is shown that basis sets of quadruple zeta quality are necessary to achieve an accuracy of 1 kcal/mol per bond in HF and DFT atomization energies. For compds. contg. third row as well as alk. and earth alk. metals it is demonstrated that the inclusion of high-lying core orbitals in the active space can be necessary for accurate correlated treatments. The QZVPP basis sets provide sufficient flexibility to polarize the core in those cases. All test calcns. indicate that the new basis sets lead to consistent accuracies in HF, DFT, or correlated treatments even in crit. cases where other basis sets may show deficiencies.
- 63Rappoport, D.; Furche, F. Property-optimized Gaussian basis sets for molecular response calculations. J. Chem. Phys. 2010, 133, 134105, DOI: 10.1063/1.3484283Google Scholar63https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXht1Chs7nI&md5=17ee007503db923dedbb77f5e2214694Property-optimized Gaussian basis sets for molecular response calculationsRappoport, Dmitrij; Furche, FilippJournal of Chemical Physics (2010), 133 (13), 134105/1-134105/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)With recent advances in electronic structure methods, first-principles calcns. of electronic response properties, such as linear and nonlinear polarizabilities, have become possible for mols. with more than 100 atoms. Basis set incompleteness is typically the main source of error in such calcns. since traditional diffuse augmented basis sets are too costly to use or suffer from near linear dependence. To address this problem, we construct the first comprehensive set of property-optimized augmented basis sets for elements H-Rn except lanthanides. The new basis sets build on the Karlsruhe segmented contracted basis sets of split-valence to quadruple-zeta valence quality and add a small no. of moderately diffuse basis functions. The exponents are detd. variationally by maximization of at. Hartree-Fock polarizabilities using anal. deriv. methods. The performance of the resulting basis sets is assessed using a set of 313 mol. static Hartree-Fock polarizabilities. The mean abs. basis set errors are 3.6%, 1.1%, and 0.3% for property-optimized basis sets of split-valence, triple-zeta, and quadruple-zeta valence quality, resp. D. functional and second-order Moller-Plesset polarizabilities show similar basis set convergence. We demonstrate the efficiency of our basis sets by computing static polarizabilities of icosahedral fullerenes up to C720 using hybrid d. functional theory. (c) 2010 American Institute of Physics.
- 64Weigend, 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, DOI: 10.1039/b508541aGoogle Scholar64https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpsFWgu7o%253D&md5=a820fb6055c993b50c405ba0fc62b194Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracyWeigend, Florian; Ahlrichs, ReinhartPhysical Chemistry Chemical Physics (2005), 7 (18), 3297-3305CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 mols. representing (nearly) all elements-except lanthanides-in their common oxidn. states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, d. functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
- 65Hättig, C. Optimization of auxiliary basis sets for RI-MP2 and RI-CC2 calculations: Core–valence and quintuple-ζ basis sets for H to Ar and QZVPP basis sets for Li to Kr. Phys. Chem. Chem. Phys. 2005, 7, 59– 66, DOI: 10.1039/B415208EGoogle ScholarThere is no corresponding record for this reference.
- 66Hellweg, A.; Hättig, C.; Höfener, S.; Klopper, W. Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to Rn. Theor. Chem. Acc. 2007, 117, 587– 597, DOI: 10.1007/s00214-007-0250-5Google Scholar66https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXivVGmtrY%253D&md5=c140888c49b17c3671a51a45497e6284Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to RnHellweg, Arnim; Haettig, Christof; Hoefener, Sebastian; Klopper, WimTheoretical Chemistry Accounts (2007), 117 (4), 587-597CODEN: TCACFW; ISSN:1432-881X. (Springer GmbH)The introduction of the resoln.-of-the-identity (RI) approxn. for electron repulsion integrals in quantum chem. calcns. requires in addn. to the orbital basis so-called auxiliary or fitting basis sets. We report here such auxiliary basis sets optimized for second-order Moller-Plesset perturbation theory for the recently published (Weigend and Ahlrichs Phys. Chem. Chem. Phys., 2005, 7, 3297-3305) segmented contracted Gaussian basis sets of split, triple-ζ and quadruple-ζ valence quality for the atoms Rb-Rn (except lanthanides). These basis sets are designed for use in connection with small-core effective core potentials including scalar relativistic corrections. Hereby accurate resoln.-of-the-identity calcns. with second-order Moller-Plesset perturbation theory (MP2) and related methods can now be performed for mols. contg. elements from H to Rn. The error of the RI approxn. has been evaluated for a test set of 385 small and medium sized mols., which represent the common oxidn. states of each element, and is compared with the one-electron basis set error, estd. based on highly accurate explicitly correlated MP2-R12 calcns. With the reported auxiliary basis sets the RI error for MP2 correlation energies is typically two orders of magnitude smaller than the one-electron basis set error, independent on the position of the atoms in the periodic table.
- 67Hellweg, A.; Rappoport, D. Development of new auxiliary basis functions of the Karlsruhe segmented contracted basis sets including diffuse basis functions (def2-SVPD, def2-TZVPPD, and def2-QZVPPD) for RI-MP2 and RI-CC calculations. Phys. Chem. Chem. Phys. 2015, 17, 1010– 1017, DOI: 10.1039/C4CP04286GGoogle ScholarThere is no corresponding record for this reference.
- 68Kaltak, M.; Klimeš, J.; Kresse, G. Low Scaling Algorithms for the Random Phase Approximation: Imaginary Time and Laplace Transformations. J. Chem. Theory Comput. 2014, 10, 2498, DOI: 10.1021/ct5001268Google Scholar68https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXmslGlu70%253D&md5=032579aea73969146da834d532c29997Low Scaling Algorithms for the Random Phase Approximation: Imaginary Time and Laplace TransformationsKaltak, Merzuk; Klimes, Jiri; Kresse, GeorgJournal of Chemical Theory and Computation (2014), 10 (6), 2498-2507CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We det. efficient imaginary frequency and imaginary time grids for second-order Moller-Plesset (MP) perturbation theory. The least-squares and Minimax quadratures are compared for periodic systems, finding that the Minimax quadrature performs slightly better for the considered materials. We show that the imaginary frequency grids developed for second order also perform well for the correlation energy in the direct RPA. Furthermore, we show that the polarizabilities on the imaginary time axis can be Fourier-transformed to the imaginary frequency domain, since the time and frequency Minimax grids are dual to each other. The same duality is obsd. for the least-squares grids. The transformation from imaginary time to imaginary frequency allows one to reduce the time complexity to cubic (in system size), so that RPA correlation energies become accessible for large systems.
- 69Kaltak, M.; Kresse, G. Minimax isometry method: A compressive sensing approach for Matsubara summation in many-body perturbation theory. Phys. Rev. B 2020, 101, 205145, DOI: 10.1103/PhysRevB.101.205145Google Scholar69https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhsFChtLjJ&md5=0e8b7dcabc1bfd1ddbac581ec0b91706Minimax isometry method: A compressive sensing approach for Matsubara summation in many-body perturbation theoryKaltak, Merzuk; Kresse, GeorgPhysical Review B (2020), 101 (20), 205145CODEN: PRBHB7; ISSN:2469-9969. (American Physical Society)We present a compressive sensing approach for the long-standing problem of Matsubara summation in many-body perturbation theory. By constructing low-dimensional, almost isometric subspaces of the Hilbert space we obtain optimum imaginary time and frequency grids that allow for extreme data compression of fermionic and bosonic functions in a broad temp. regime. The method is applied to the random phase and self-consistent GW approxn. of the grand potential. Integration and transformation errors are investigated for Si and SrVO3.
- 70Andrae, D.; Häußermann, U.; Dolg, M.; Stoll, H.; Preuß, H. Energy-adjusted ab initio pseudopotentials for the second and third row transition elements. Theor. Chim. Acta 1990, 77, 123– 141, DOI: 10.1007/BF01114537Google Scholar70https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3cXkt12ntLo%253D&md5=8203c6bc6149924cbd4b23b3063715e1Energy-adjusted ab initio pseudopotentials for the second and third row transition elementsAndrae, D.; Haeussermann, U.; Dolg, M.; Stoll, H.; Preuss, H.Theoretica Chimica Acta (1990), 77 (2), 123-41CODEN: TCHAAM; ISSN:0040-5744.Nonrelativistic and quasirelativistic ab initio pseudopotentials substituting the M(z-28)+-core orbitals of the second row transition elements and M(z-60)+-core orbitals of the third row transition elements, resp., and optimized (8s7p6d)/[6s5p3d]-GTO valence basis sets for use in mol. calcns. have been generated. Addnl., corresponding spin-orbit operators have also been derived. At. excitation and ionization energies from numerical HF as well as from SCF pseudopotential calcns. using the derived basis sets differ in most cases by less than 0.1 eV from corresponding numerical all-electron results. Spin-orbit splittings for low-lying states are in reasonable agreement with corresponding all-electron Dirac-Fock (DF) results.
- 71Metz, B.; Stoll, H.; Dolg, M. Small-core multiconfiguration-Dirac–Hartree–Fock-adjusted pseudopotentials for post-d main group elements: Application to PbH and PbO. J. Chem. Phys. 2000, 113, 2563– 2569, DOI: 10.1063/1.1305880Google Scholar71https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXlsVWmsrg%253D&md5=3b605722931765b461713ea2fe318c7eSmall-core multiconfiguration-Dirac-Hartree-Fock-adjusted pseudopotentials for post-d main group elements: Application to PbH and PbOMetz, Bernhard; Stoll, Hermann; Dolg, MichaelJournal of Chemical Physics (2000), 113 (7), 2563-2569CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Relativistic pseudopotentials (PPs) of the energy-consistent variety have been generated for the post-d group 13-15 elements, by adjustment to multiconfiguration Dirac-Hartree-Fock data based on the Dirac-Coulomb-Breit Hamiltonian. The outer-core (n-1)spd shells are explicitly treated together with the nsp valence shell, with these PPs, and the implications of the small-core choice are discussed by comparison to a corresponding large-core PP, in the case of Pb. Results from valence ab initio one- and two-component calcns. using both PPs are presented for the fine-structure splitting of the ns2np2 ground-state configuration of the Pb atom, and for spectroscopic consts. of PbH (X 2Π1/2, 2Π3/2) and PbO (X 1Σ+). In addn., a combination of small-core and large-core PPs has been explored in spin-free-state shifted calcns. for the above mols.
- 72Peterson, K. A. Systematically convergent basis sets with relativistic pseudopotentials. I. Correlation consistent basis sets for the post-d group 13–15 elements. J. Chem. Phys. 2003, 119, 11099– 11112, DOI: 10.1063/1.1622923Google Scholar72https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXovFyms74%253D&md5=ab472eed340dd9b9ac5e2d7d2e76df7aSystematically convergent basis sets with relativistic pseudopotentials. I. Correlation consistent basis sets for the post-d group 13-15 elementsPeterson, Kirk A.Journal of Chemical Physics (2003), 119 (21), 11099-11112CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)New correlation consistent-like basis sets have been developed for the post-d group 13-15 elements (Ga-As, In-Sb, Tl-Bi) employing accurate, small-core relativistic pseudopotentials. The resulting basis sets, which are denoted cc-pVnZ-PP, are appropriate for valence electron correlation and range in size from (8s7p7d)/[4s3p2d] for the cc-pVDZ-PP to (16s13p12d3f2g1h)/[7s7p5d3f2g1h] for the cc-pV5Z-PP sets. Benchmark calcns. on selected diat. mols. (As2, Sb2, Bi2, AsN, SbN, BiN, GeO, SnO, PbO, GaCl, InCl, TlCl, GaH, InH, and TlH) are reported using these new basis sets at the coupled cluster level of theory. Much like their all-electron counterparts, the cc-pVnZ-PP basis sets yield systematic convergence of total energies and spectroscopic consts. In several cases all-electron benchmark calcns. were also carried out for comparison. The results from the pseudopotential and all-electron calcns. were nearly identical when scalar relativity was accurately included in the all-electron work. Diffuse-augmented basis sets, aug-cc-pVnZ-PP, have also been developed and have been used in calcns. of the at. electron affinities.
- 73Peterson, K. A.; Figgen, D.; Goll, E.; Stoll, H.; Dolg, M. Systematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis sets for the post-d group 16–18 elements. J. Chem. Phys. 2003, 119, 11113– 11123, DOI: 10.1063/1.1622924Google Scholar73https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXovFyms7w%253D&md5=2f5869fe14c68b766761d8d160acce6aSystematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis sets for the post-d group 16-18 elementsPeterson, Kirk A.; Figgen, Detlev; Goll, Erich; Stoll, Hermann; Dolg, MichaelJournal of Chemical Physics (2003), 119 (21), 11113-11123CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A series of correlation consistent basis sets have been developed for the post-d group 16-18 elements in conjunction with small-core relativistic pseudopotentials of the energy-consistent variety. The latter were adjusted to multiconfiguration Dirac-Hartree-Fock data based on the Dirac-Coulomb-Breit Hamiltonian. The outer-core (n-1)spd shells are explicitly treated together with the nsp valence shell with these PPs. The accompanying cc-pVnZ-PP and aug-cc-pVnZ-PP basis sets range in size from DZ to 5Z quality and yield systematic convergence of both Hartree-Fock and correlated total energies. In addn. to the calcn. of at. electron affinities and dipole polarizabilities of the rare gas atoms, numerous mol. benchmark calcns. (HBr, HI, HAt, Br2, I2, At2, SiSe, SiTe, SiPo, KrH+, XeH+, and RnH+) are also reported at the coupled cluster level of theory. For the purposes of comparison, all-electron calcns. using the Douglas-Kroll-Hess Hamiltonian have also been carried out for the halogen-contg. mols. using basis sets of 5Z quality.
- 74Morgante, P.; Peverati, R. Statistically representative databases for density functional theory via data science. Phys. Chem. Chem. Phys. 2019, 21, 19092– 19103, DOI: 10.1039/C9CP03211HGoogle Scholar74https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhsFygsL%252FO&md5=ed89b2e280632ec164cf0101e07f3fd5Statistically representative databases for density functional theory via data scienceMorgante, Pierpaolo; Peverati, RobertoPhysical Chemistry Chemical Physics (2019), 21 (35), 19092-19103CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)The amt. of data and no. of databases for the assessment and parameterization of d. functional theory methods has grown substantially in the past two decades. In this work, we introduce a novel cluster anal. technique for d. functional theory calcns. of the electronic structure of atoms and mols. with the goal of creating new statistically significant databases with broad chem. scope, and a manageable no. of data-points. By analyzing without a priori chem. assumptions a population of almost 350k data-points, we create a new database called ASCDB contg. only 200 data-points. This new database holds the same chem. information as the larger population of data from which it is obtained, but with a computational cost that is reduced by several orders of magnitude. The labeling of the significant chem. properties is performed a posteriori on the resulting 16 subsets, classifying them into four areas of chem. importance: non-covalent interactions, thermochem., non-local effects, and unbiased calcns. The anal. of the results and their transferability shows that ASCDB is capable of providing the same information as that of the larger collection of data-such as GMTKN55, MGCDB84, and Minnesota 2015B-for several d. functional theory methods and basis sets. In light of these results, we suggest the use of this new small database as a first inexpensive tool for the evaluation and parameterization of electronic structure theory methods.
- 75Morgante, P.; Peverati, R. ACCDB: A collection of chemistry databases for broad computational purposes. J. Comput. Chem. 2019, 40, 839– 848, DOI: 10.1002/jcc.25761Google Scholar75https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXisF2htLzI&md5=2e99af73352558db4926d2065f47dcdcACCDB: A collection of chemistry databases for broad computational purposesMorgante, Pierpaolo; Peverati, RobertoJournal of Computational Chemistry (2019), 40 (6), 839-848CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)The importance of databases of reliable and accurate data in chem. has substantially increased in the past two decades. Their main usage is to parametrize electronic structure theory methods, and to assess their capabilities and accuracy for a broad set of chem. problems. The collection we present here-ACCDB-includes data from 16 different research groups, for a total of 44,931 unique ref. data points, all at a level of theory significantly higher than d. functional theory, and covering most of the periodic table. It is composed of five databases taken from literature (GMTKN, MGCDB84, Minnesota2015, DP284, and W4-17), two newly developed reaction energy databases (W4-17-RE and MN-RE), and a new collection of databases contg. transition metals. A set of expandable software tools for its manipulation is also presented here for the first time, as well as a case study where ACCDB is used for benchmarking com. CPUs for chem. calcns. © 2018 Wiley Periodicals, Inc.
- 76Fritsch, F. N.; Butland, J. A Method for Constructing Local Monotone Piecewise Cubic Interpolants. SIAM J. Sci. Stat. Comp. 1984, 5, 300, DOI: 10.1137/0905021Google ScholarThere is no corresponding record for this reference.
- 77Peverati, R. Fitting elephants in the density functionals zoo: Statistical criteria for the evaluation of density functional theory methods as a suitable replacement for counting parameters. Int. J. Quantum Chem. 2021, 121, e26379, DOI: 10.1002/qua.26379Google Scholar77https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtlersrzK&md5=278ba83220ff86d50f5309e7f76151ecFitting elephants in the density functionals zoo: Statistical criteria for the evaluation of density functional theory methods as a suitable replacement for counting parametersPeverati, RobertoInternational Journal of Quantum Chemistry (2021), 121 (1), e26379CODEN: IJQCB2; ISSN:0020-7608. (John Wiley & Sons, Inc.)Counting parameters has become customary in the d. functional theory community as a way to infer the transferability of popular approxns. to the exchange-correlation functionals. Recent work in data science, however, has demonstrated that the no. of parameters of a fitted model is not related to the complexity of the model itself, nor to its eventual overfitting. Using similar arguments, here, we show that it is possible to represent every modern exchange-correlation functional approxns. using just one single parameter. This procedure proves the futility of the no. of parameters as a measure of transferability. To counteract this shortcoming, we introduce and analyze the performance of three statistical criteria for the evaluation of the transferability of exchange-correlation functionals. The three criteria are called Akaike information criterion, Vapnik-Chervonenkis criterion, and cross-validation criterion and are used in a preliminary assessment to rank 60 exchange-correlation functional approxns. using the ASCDB database of chem. data.
- 78Korth, M.; Grimme, S. Mindless DFT Benchmarking. J. Chem. Theory Comput. 2009, 5, 993– 1003, DOI: 10.1021/ct800511qGoogle Scholar78https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXis12itbY%253D&md5=d3a2772137b30cf4f0343825104798a5"Mindless" DFT BenchmarkingKorth, Martin; Grimme, StefanJournal of Chemical Theory and Computation (2009), 5 (4), 993-1003CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A diversity-oriented approach for the generation of thermochem. benchmark sets is presented. Test sets consisting of randomly generated "artificial mols." (AMs) are proposed that rely on systematic constraints rather than uncontrolled chem. biases. In this way, the narrow structural space of chem. intuition is opened up and electronically difficult cases can be produced in an unforeseeable manner. For the calcn. of chem. meaningful relative energies, AMs are systematically decompd. into small mols. (hydrides and diatomics). Two different example test sets contg. eight-atom, single-ref., main group AMs with chem. very diverse and unusual structures are generated. Highly accurate all-electron, estd. CCSD(T)/complete basis set ref. energies are also provided. They are used to benchmark the d. functionals S-VWN, BP86, B-LYP, B97-D, PBE, TPSS, PBEh, BH-LYP, B3-PW91, B3-LYP, B2-PLYP, B2GP-PLYP, BMK, MPW1B95, M05, M05-2X, PW6B95, M06, M06-L, and M06-2X. In selected cases, an empirical dispersion correction (DFT-D) has been applied. Due to the compn. of the sets, it is expected that a good performance indicates "robustness" in many different chem. applications. The results of a statistical anal. of the errors for the entire set with 165 entries (av. reaction energy of 117 kcal/mol, dubbed as the MB08-165 set) perfectly fit to the "Jacob's ladder" metaphor for the ordering of d. functionals according to their theor. complexity. The mean abs. deviation (MAD) decreases very strongly from LDA (20 kcal/mol) to GGAs (MAD of about 10 kcal/mol) but then was less pronounced to hybrid-GGAs (MAD of about 6-8 kcal/mol). The best performance (MAD of 4.1-4.2 kcal/mol) is found for the (fifth-rung) double-hybrid functionals B2-PLYP-D and B2GP-PLYP-D, followed by the M06-2X meta-hybrid (MAD of 4.8 kcal/mol). The significance of the proposed approach for thermodn. benchmarking is discussed and related to the obsd. performance ranking also regarding wave function based methods.
- 79Broyden, C. G. The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations. IMA J. Appl. Math. 1970, 6, 76, DOI: 10.1093/imamat/6.1.76Google ScholarThere is no corresponding record for this reference.
- 80Fletcher, R. A new approach to variable metric algorithms. Comput. J. 1970, 13, 317, DOI: 10.1093/comjnl/13.3.317Google ScholarThere is no corresponding record for this reference.
- 81Goldfarb, D. A family of variable-metric methods derived by variational methods. Math. Comput. 1970, 24, 23, DOI: 10.1090/S0025-5718-1970-0258249-6Google ScholarThere is no corresponding record for this reference.
- 82Shanno, D. F. Conditioning of quasi-Newton methods for function minimization. Math. Comput. 1970, 24, 647, DOI: 10.1090/S0025-5718-1970-0274029-XGoogle ScholarThere is no corresponding record for this reference.
- 83Parthiban, S.; Martin, J. M. L. Assessment of W1 and W2 theories for the computation of electron affinities, ionization potentials, heats of formation, and proton affinities. J. Chem. Phys. 2001, 114, 6014– 6029, DOI: 10.1063/1.1356014Google Scholar83https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXitlKls7w%253D&md5=c52c0d86eab1a2f42f893ecc19bac2adAssessment of W1 and W2 theories for the computation of electron affinities, ionization potentials, heats of formation, and proton affinitiesParthiban, Srinivasan; Martin, Jan M. L.Journal of Chemical Physics (2001), 114 (14), 6014-6029CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The performance of two recent ab initio computational thermochem. schemes, W1 and W2 theory [J. M. L. Martin and G. de Oliveira, J. Chem. Phys. 111, 1843 (1999)], is assessed for an enlarged sample of thermochem. data consisting of the ionization potentials and electron affinities in the G2-1 and G2-2 sets, as well as the heats of formation in the G2-1 and a subset of the G2-2 set. We find W1 theory to be several times more accurate for ionization potentials and electron affinities than commonly used (and less expensive) computational thermochem. schemes such as G2, G3, and CBS-QB3: W2 theory represents a slight improvement for electron affinities but no significant one for ionization potentials. The use of a two-point A+B/L5 rather than a three-point A+B/CL extrapolation for the SCF component greatly enhances the numerical stability of the W1 method for systems with slow basis set convergence. Inclusion of first-order spin-orbit coupling is essential for accurate ionization potentials and electron affinities involving degenerate electronic states: Inner-shell correlation is somewhat more important for ionization potentials than for electron affinities, while scalar relativistic effects are required for the highest accuracy. The mean deviation from expt. for the G2-1 heats of formation is within the av. exptl. uncertainty. W1 theory appears to be a valuable tool for obtaining benchmark quality proton affinities.
- 84Karton, A.; Sylvetsky, N.; Martin, J. M. L. W4–17: A diverse and high-confidence dataset of atomization energies for benchmarking high-level electronic structure methods. J. Comput. Chem. 2017, 38, 2063, DOI: 10.1002/jcc.24854Google Scholar84https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtV2gsL3K&md5=4e8a303a60ccb1e846019a1f69d61d19W4-17: A diverse and high-confidence dataset of atomization energies for benchmarking high-level electronic structure methodsKarton, Amir; Sylvetsky, Nitai; Martin, Jan M. L.Journal of Computational Chemistry (2017), 38 (24), 2063-2075CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Atomization reactions are among the most challenging tests for electronic structure methods. We use the first-principles Weizmann-4 (W4) computational thermochem. protocol to generate the W4-17 dataset of 200 total atomization energies (TAEs) with 3σ confidence intervals of 1 kJ mol-1. W4-17 is an extension of the earlier W4-11 dataset; it includes first- and second-row mols. and radicals with up to eight non-hydrogen atoms. These cover a broad spectrum of bonding situations and multireference character, and as such are an excellent benchmark for the parameterization and validation of highly accurate ab initio methods (e.g., CCSD(T) composite procedures) and double-hybrid d. functional theory (DHDFT) methods. The W4-17 dataset contains two subsets (i) a non-multireference subset of 183 systems characterized by dynamical or moderate nondynamical correlation effects (denoted W4-17-nonMR) and (ii) a highly multireference subset of 17 systems (W4-17-MR). We use these databases to evaluate the performance of a wide range of CCSD(T) composite procedures (e.g., G4, G4(MP2), G4(MP2)-6X, ROG4(MP2)-6X, CBS-QB3, ROCBS-QB3, CBS-APNO, ccCA-PS3, W1, W2, W1-F12, W2-F12, W1X-1, and W2X) and DHDFT methods (e.g., B2-PLYP, B2GP-PLYP, B2K-PLYP, DSD-BLYP, DSD-PBEP86, PWPB95, ωB97X-2(LP), and ωB97X-2(TQZ)). © 2017 Wiley Periodicals, Inc.
- 85Karton, A.; Daon, S.; Martin, J. M. W4–11: A high-confidence benchmark dataset for computational thermochemistry derived from first-principles W4 data. Chem. Phys. Lett. 2011, 510, 165, DOI: 10.1016/j.cplett.2011.05.007Google Scholar85https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXnvF2qsLs%253D&md5=41b53b92e3d9c4f450bc55a5e5df4eb7W4-11: A high-confidence benchmark dataset for computational thermochemistry derived from first-principles W4 dataKarton, Amir; Daon, Shauli; Martin, Jan M. L.Chemical Physics Letters (2011), 510 (4-6), 165-178CODEN: CHPLBC; ISSN:0009-2614. (Elsevier B.V.)We show that the purely first-principles Weizmann-4 (W4) computational thermochem. method developed in our group can reproduce available Active Thermochem. Tables atomization energies for 35 mols. with a 3σ uncertainty of under 1 kJ/mol. We then employ this method to generate the W4-11 dataset of 140 total atomization energies of small first-and second-row mols. and radicals. These cover a broad spectrum of bonding situations and multireference character, and as such are an excellent, quasi-automated benchmark (available electronically as Supporting Information) for parametrization and validation of more approx. methods (such as DFT functionals and composite methods). Secondary contributions such as relativity can be included or omitted at will, unlike with exptl. data. A broad variety of more approx. methods is assessed against the W4-11 benchmark and recommendations are made.
- 86Margraf, J. T.; Ranasinghe, D. S.; Bartlett, R. J. Automatic generation of reaction energy databases from highly accurate atomization energy benchmark sets. Phys. Chem. Chem. Phys. 2017, 19, 9798, DOI: 10.1039/C7CP00757DGoogle Scholar86https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXksVOqur4%253D&md5=4f425ab2aa394430e939871dcf8925b0Automatic generation of reaction energy databases from highly accurate atomization energy benchmark setsMargraf, Johannes T.; Ranasinghe, Duminda S.; Bartlett, Rodney J.Physical Chemistry Chemical Physics (2017), 19 (15), 9798-9805CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)In this contribution, we discuss how reaction energy benchmark sets can automatically be created from arbitrary atomization energy databases. As an example, over 11 000 reaction energies derived from the W4-11 database, as well as some relevant subsets are reported. Importantly, there is only very modest computational overhead involved in computing >11 000 reaction energies compared to 140 atomization energies, since the rate-detg. step for either benchmark is performing the same 140 quantum chem. calcns. The performance of commonly used electronic structure methods for the new database is analyzed. This allows investigating the relationship between the performances for atomization and reaction energy benchmarks based on an identical set of mols. The atomization energy is found to be a weak predictor for the overall usefulness of a method. The performance of d. functional approxns. in light of the no. of empirically optimized parameters used in their design is also discussed.
- 87Henderson, T. M.; Scuseria, G. E. The connection between self-interaction and static correlation: A random phase approximation perspective. Mol. Phys. 2010, 108, 2511, DOI: 10.1080/00268976.2010.507227Google Scholar87https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhtl2is7rJ&md5=62434dbb5a179a1b99b4dc07e0c25682The connection between self-interaction and static correlation: a random phase approximation perspectiveHenderson, Thomas M.; Scuseria, Gustavo E.Molecular Physics (2010), 108 (19-20), 2511-2517CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)Semi-local d. functional theory suggests a connection between static correlation and self-interaction. It is difficult to make such a connection from the wave function theory perspective, since few wave function methods permit self-interaction error. However, the RPA for ground-state correlation, which has a wave function derivation, does include self-interaction in its direct (Hartree) variant. This variant also describes left-right correlation. The self-interaction can be removed by means of second-order screened exchange; however, this also has neg. consequences for the description of static correlation. This paper discusses the connection between the two concepts (static correlation and self-interaction) from the perspective provided by the RPA.
- 88Grüning, M.; Gritsenko, O. V.; van Gisbergen, S. J. A.; Baerends, E. J. The Failure of Generalized Gradient Approximations (GGAs) and Meta-GGAs for the Two-Center Three-Electron Bonds in He2+, (H2O)2+, and (NH3)2+. J. Phys. Chem. A 2001, 105, 9211, DOI: 10.1021/jp011239kGoogle Scholar88https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXmvVWrsrc%253D&md5=e3944e590f1fc9aa215d593b809725e5The Failure of Generalized Gradient Approximations (GGAs) and Meta-GGAs for the Two-Center Three-Electron Bonds in He2+, (H2O)2+, and (NH3)2+Gruening, M.; Gritsenko, O. V.; van Gisbergen, S. J. A.; Baerends, E. J.Journal of Physical Chemistry A (2001), 105 (40), 9211-9218CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)The radical cations He2+, (H2O)2+, and (NH3)2+ with two-center three-electron A-A bonds are investigated at the CI, accurate Kohn-Sham (KS), generalized gradient approxn. (GGA), and meta-GGA levels. Assessment of seven different GGA and six meta-GGA methods shows that the A2+ systems remain a difficult case for d. functional theory (DFT). All methods tested consistently overestimate the stability of A2+: the corresponding De errors decrease for more diffuse valence densities in the series He2+ > (H2O)2+ > (NH3)2+. Upon comparison to the energy terms of the accurate Kohn-Sham solns., the approx. exchange functionals are found to be responsible for the errors of GGA-type methods, which characteristically overestimate the exchange in A2+. These so-called exchange functionals implicitly use localized holes. Such localized holes do occur if there is left-right correlation, i.e., the exchange functionals then also describe non-dynamical correlation. However, in the hemibonded A2+ systems the typical mol. (left-right, non-dynamical) correlation of the two-electron pair bond is absent. The non-dynamical correlation built into the exchange functionals is then spurious and yields too low energies.
- 89Peterson, K. A.; Woon, D. E.; Dunning, T. H. Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H+H2→H2+H reaction. J. Chem. Phys. 1994, 100, 7410– 7415, DOI: 10.1063/1.466884Google Scholar89https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2cXjtlWmurY%253D&md5=c8b0f6e0bd2b066d23bc77bddaf0be36Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H + H2 → H2 + H reactionPeterson, Kirk A.; Woon, David E.; Dunning, Thom H., Jr.Journal of Chemical Physics (1994), 100 (10), 7410-15CODEN: JCPSA6; ISSN:0021-9606.Using systematic sequences of correlation consistent Gaussian basis sets from double to sextuple zeta quality, the classical barrier height of the H + H2 exchange reaction has been calcd. by multireference CI (MRCI) methods. The MRCI calcns. for collinear H3 have also been calibrated against large basis set full CI (FCI) results, which demonstrate that the MRCI treatment leads to energies less than 1 μhartree (<0.001 kcal/mol) above the FCI energies. The dependence of both the H2 and H3 total energies on the basis set is found to be very regular, and this behavior has been used to extrapolate to the complete basis set (CBS) limits. The resulting est. of the H-H-H- CBS limit yields a classical barrier height, relative to exact H + H2, of 9.60 ± 0.02 kcal/mol; the best directly calcd. value for the barrier is equal to 9.62 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo calcns.
- 90Helgaker, T.; Klopper, W.; Koch, H.; Noga, J. Basis-set convergence of correlated calculations on water. J. Chem. Phys. 1997, 106, 9639– 9646, DOI: 10.1063/1.473863Google Scholar90https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXjvVCgu78%253D&md5=f4689c1b38fe30eb721e9cd7d607bdf7Basis-set convergence of correlated calculations on waterHelgaker, Trygve; Klopper, Wim; Koch, Henrik; Noga, JozefJournal of Chemical Physics (1997), 106 (23), 9639-9646CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The basis-set convergence of the electronic correlation energy in the water mol. is investigated at the second-order Moller-Plesset level and at the coupled-cluster singles-and-doubles level with and without perturbative triples corrections applied. The basis-set limits of the correlation energy are established to within 2mEh by means of (1) extrapolations from sequences of calcns. using correlation-consistent basis sets and (2) from explicitly correlated calcns. employing terms linear in the inter-electronic distances rij. For the extrapolations to the basis-set limit of the correlation energies, fits of the form a + bX-3 (where X is two for double-zeta sets, three for triple-zeta sets, etc.) are found to be useful. CCSD(T) calcns. involving as many as 492 AOs are reported.
- 91Halkier, A.; Helgaker, T.; Jo̷rgensen, P.; Klopper, W.; Koch, H.; Olsen, J.; Wilson, A. K. Basis-set convergence in correlated calculations on Ne, N2, and H2O. Chem. Phys. Lett. 1998, 286, 243– 252, DOI: 10.1016/S0009-2614(98)00111-0Google Scholar91https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXitVGqsLo%253D&md5=04274821d9c7fa664e9588855ed9a061Basis-set convergence in correlated calculations on Ne, N2, and H2OHalkier, Asger; Helgaker, Trygve; Jorgensen, Poul; Klopper, Wim; Koch, Henrik; Olsen, Jeppe; Wilson, Angela K.Chemical Physics Letters (1998), 286 (3,4), 243-252CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)Valence and all-electron correlation energies of Ne, N2, and H2O at fixed exptl. geometries are computed at the levels of second-order perturbation theory (MP2) and coupled cluster theory with singles and doubles excitations (CCSD), and singles and doubles excitations with a perturbative triples correction (CCSD(T)). Correlation-consistent polarized valence and core-valence basis sets up to sextuple zeta quality are employed. Guided by basis-set limits established by rij-dependent methods, a no. of extrapolation schemes for use with the correlation-consistent basis sets are investigated. Among the schemes considered here, a linear least-squares procedure applied to the quintuple and sextuple zeta results yields the most accurate extrapolations.
- 92Dunning, J.; Thom, H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007– 1023, DOI: 10.1063/1.456153Google Scholar92https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL1MXksVGmtrk%253D&md5=c6cd67a3748dc61692a9cb622d2694a0Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogenDunning, Thom H., Jr.Journal of Chemical Physics (1989), 90 (2), 1007-23CODEN: JCPSA6; ISSN:0021-9606.Guided by the calcns. on oxygen in the literature, basis sets for use in correlated at. and mol. calcns. were developed for all of the first row atoms from boron through neon, and for hydrogen. As in the oxygen atom calcns., the incremental energy lowerings, due to the addn. of correlating functions, fall into distinct groups. This leads to the concept of correlation-consistent basis sets, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation-consistent sets are given for all of the atoms considered. The most accurate sets detd. in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding at.-natural-orbital sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estd. that this set yields 94-97% of the total (HF + 1 + 2) correlation energy for the atoms neon through boron.
- 93Kendall, R. A.; Dunning, T. H.; Harrison, R. J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796– 6806, DOI: 10.1063/1.462569Google Scholar93https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK38XktFClurw%253D&md5=948a06eee10604a8fa37eae2b2ada4beElectron affinities of the first-row atoms revisited. Systematic basis sets and wave functionsKendall, Rick A.; Dunning, Thom H., Jr.; Harrison, Robert J.Journal of Chemical Physics (1992), 96 (9), 6796-806CODEN: JCPSA6; ISSN:0021-9606.The authors describe a reliable procedure for calcg. the electron affinity of an atom and present results for H, B, C, O, and F (H is included for completeness). This procedure involves the use of the recently proposed correlation-consistent basis sets augmented with functions to describe the more diffuse character of the at. anion coupled with a straightforward, uniform expansion of the ref. space for multireference singles and doubles configuration-interaction (MRSD-CI) calcns. A comparison is given with previous results and with corresponding full CI calcns. The most accurate EAs obtained from the MRSD-CI calcns. are (with exptl. values in parentheses): H 0.740 eV (0.754), B 0.258 (0.277), C 1.245 (1.263), O 1.384 (1.461), and F 3.337 (3.401). The EAs obtained from the MR-SDCI calcns. differ by less than 0.03 eV from those predicted by the full CI calcns.
- 94Peterson, K. A.; Dunning, T. H. Accurate correlation consistent basis sets for molecular core–valence correlation effects: The second row atoms Al–Ar, and the first row atoms B–Ne revisited. J. Chem. Phys. 2002, 117, 10548– 10560, DOI: 10.1063/1.1520138Google Scholar94https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XptFGgu70%253D&md5=146a124710aa803859c1a67e09b45b0aAccurate correlation consistent basis sets for molecular core-valence correlation effects: The second row atoms Al-Ar, and the first row atoms B-Ne revisitedPeterson, Kirk A.; Dunning, Thom H., Jr.Journal of Chemical Physics (2002), 117 (23), 10548-10560CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Correlation consistent basis sets for accurately describing core-core and core-valence correlation effects in atoms and mols. have been developed for the second row atoms Al-Ar. Two different optimization strategies were investigated, which led to two families of core-valence basis sets when the optimized functions were added to the std. correlation consistent basis sets (cc-pVnZ). In the first case, the exponents of the augmenting primitive Gaussian functions were optimized with respect to the difference between all-electron and valence-electron correlated calcns., i.e., for the core-core plus core-valence correlation energy. This yielded the cc-pCVnZ family of basis sets, which are analogous to the sets developed previously for the first row atoms [D. E. Woon and T. H. Dunning, Jr., J. Chem. Phys. 103, 4572 (1995)]. Although the cc-pCVnZ sets exhibit systematic convergence to the all-electron correlation energy at the complete basis set limit, the intershell (core-valence) correlation energy converges more slowly than the intrashell (core-core) correlation energy. Since the effect of including the core electrons on the calcn. of mol. properties tends to be dominated by core-valence correlation effects, a second scheme for detg. the augmenting functions was investigated. In this approach, the exponents of the functions to be added to the cc-pVnZ sets were optimized with respect to just the core-valence (intershell) correlation energy, except that a small amt. of core-core correlation energy was included in order to ensure systematic convergence to the complete basis set limit. These new sets, denoted weighted core-valence basis sets (cc-pwCVnZ), significantly improve the convergence of many mol. properties with n. Optimum cc-pwCVnZ sets for the first-row atoms were also developed and show similar advantages. Both the cc-pCVnZ and cc-pwCVnZ basis sets were benchmarked in coupled cluster [CCSD(T)] calcns. on a series of second row homonuclear diat. mols. (Al2, Si2, P2, S2, and Cl2), as well as on selected diat. mols. involving first row atoms (CO, SiO, PN, and BCl). For the calcn. of core correlation effects on energetic and spectroscopic properties, the cc-pwCVnZ basis sets are recommended over the cc-pCVnZ ones.
- 95Weigend, F.; Köhn, A.; Hättig, C. Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations. J. Chem. Phys. 2002, 116, 3175– 3183, DOI: 10.1063/1.1445115Google Scholar95https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XhtlSiu7k%253D&md5=0130fa656254a693e80d4be6b0f442b8Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculationsWeigend, Florian; Kohn, Andreas; Hattig, ChristofJournal of Chemical Physics (2002), 116 (8), 3175-3183CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The convergence of the second-order Moller-Plesset perturbation theory (MP2) correlation energy with the cardinal no. X is investigated for the correlation consistent basis-set series cc-pVXZ and cc-pV(X+d)Z. For the aug-cc-pVXZ and aug-cc-pV(X+d)Z series the convergence of the MP2 correlation contribution to the dipole moment is studied. It is found that, when d-shell electrons cannot be frozen, the cc-pVXZ and aug-cc-pVXZ basis sets converge much slower for third-row elements then they do for first- and second-row elements. Based on the results of these studies criteria are deduced for the accuracy of auxiliary basis sets used in the resoln. of the identity (RI) approxn. for electron repulsion integrals. Optimized auxiliary basis sets for RI-MP2 calcns. fulfilling these criteria are reported for the sets cc-pVXZ, cc-pV(X+d)Z, aug-cc-pVXZ, and aug-cc-pV(X+d)Z with X=D, T, and Q. For all basis sets the RI error in the MP2 correlation energy is more than two orders of magnitude smaller than the usual basis-set error. For the auxiliary aug-cc-pVXZ and aug-cc-pV(X+d)Z sets the RI error in the MP2 correlation contribution to the dipole moment is one order of magnitude smaller than the usual basis set error. Therefore extrapolations towards the basis-set limit are possible within the RI approxn. for both energies and properties. The redn. in CPU time obtained with the RI approxn. increases rapidly with basis set size. For the cc-pVQZ basis an acceleration by a factor of up to 170 is obsd.
- 96Boys, S. F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553– 566, DOI: 10.1080/00268977000101561Google Scholar96https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXht1alt7fM&md5=c4f847835a638b4fc449fae688c557e9The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errorsBoys, S. F.; Bernardi, F.Molecular Physics (1970), 19 (4), 553-566CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)A new direct difference method for the computation of mol. interactions has been based on a bivariational transcorrelated treatment, together with special methods for the balancing of other errors. It appears that these new features can give a strong redn. in the error of the interaction energy, and they seem to be particularly suitable for computations in the important region near the min. energy. It has been generally accepted that this problem is dominated by unresolved difficulties and the relation of the new methods of these apparent difficulties is analyzed here.
- 97Gdanitz, R. J. An accurate interaction potential for neon dimer (Ne2). Chem. Phys. Lett. 2001, 348, 67– 74, DOI: 10.1016/S0009-2614(01)01088-0Google ScholarThere is no corresponding record for this reference.
- 98Perdew, J. P. Proceedings of the 75. WE-Heraeus-Seminar and 21st Annual International Symposium on Electronic Structure of Solids; Akademie Verlag: Berlin, 1991; p 11.Google ScholarThere is no corresponding record for this reference.
- 99Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 1992, 46, 6671– 6687, DOI: 10.1103/PhysRevB.46.6671Google Scholar99https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK38XlvFyks7c%253D&md5=e52b617f523b1a8c545707544ce01d05Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlationPerdew, John P.; Chevary, J. A.; Vosko, S. H.; Jackson, Koblar A.; Pederson, Mark R.; Singh, D. J.; Fiolhais, CarlosPhysical Review B: Condensed Matter and Materials Physics (1992), 46 (11), 6671-87CODEN: PRBMDO; ISSN:0163-1829.Generalized gradient approxns. (GGA's) seek to improve upon the accuracy of the local-spin-d. (LSD) approxn. in electronic-structure calcns. Perdew and Wang have developed a GGA based on real-space cutoff of the spurious long-range components of the second-order gradient expansion for the exchange-correlation hole. Authors have found that this d. functional performs well in numerical tests for a variety of systems: Total energies of 30 atoms are highly accurate. Ionization energies and electron affinities are improved in a statistical sense, although significant interconfigurational and interterm errors remain. Accurate atomization energies are found for seven hydrocarbon mols., with a rms error per bond of 0.1 eV, compared with 0.7 eV for the LSD approxn. and 2.4 eV for the Hartree-Fock approxn. For atoms and mols., there is a cancellation of error between d. functionals for exchange and correlation, which is most striking whenever the Hartree-Fock result is furthest from expt. The surprising LSD underestimation of the lattice consts. of Li and Na by 3-4% is cor., and the magnetic ground state of solid Fe is restored. The work function, surface energy (neglecting the long-range contribution), and curvature energy of a metallic surface are all slightly reduced in comparison with LSD. Taking account of the pos. long-range contribution, authors find surface and curvature energies in good agreement with exptl. or exact values. Finally, a way is found to visualize and understand the nonlocality of exchange and correlation, its origins, and its phys. effects.
- 100Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098– 3100, DOI: 10.1103/PhysRevA.38.3098Google Scholar100https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL1cXmtlOhsLo%253D&md5=d4d219c134a5a90f689a8abed04d82ccDensity-functional exchange-energy approximation with correct asymptotic behaviorBecke, A. D.Physical Review A: Atomic, Molecular, and Optical Physics (1988), 38 (6), 3098-100CODEN: PLRAAN; ISSN:0556-2791.Current gradient-cor. d.-functional approxns. for the exchange energies of at. and mol. systems fail to reproduce the correct 1/r asymptotic behavior of the exchange-energy d. A gradient-cor. exchange-energy functional is given with the proper asymptotic limit. This functional, contg. only one parameter, fits the exact Hartree-Fock exchange energies of a wide variety of at. systems with remarkable accuracy, surpassing the performance of previous functionals contg. two parameters or more.
- 101Zhang, Y.; Pan, W.; Yang, W. Describing van der Waals interaction in diatomic molecules with generalized gradient approximations: The role of the exchange functional. J. Chem. Phys. 1997, 107, 7921– 7925, DOI: 10.1063/1.475105Google Scholar101https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXntFGisrY%253D&md5=ed2e0b29f00aa7dedbc4eb461bdbe0b1Describing van der Waals Interaction in diatomic molecules with generalized gradient approximations: the role of the exchange functionalZhang, Yingkai; Pan, Wei; Yang, WeitaoJournal of Chemical Physics (1997), 107 (19), 7921-7925CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Generalized gradient approxns. have been used to calc. the potential energy curves for six rare gas diat. mols. (He2, Ne2, Ar2, HeNe, ArHe, ArNe). Several generalized gradient approxns. are found to provide a good description of binding in these diat. mols. and show a significant improvement over the local d. approxn. in the prediction of bond lengths and dissocn. energies. It is shown here that the behavior of an exchange functional in the region of small d. and large d. gradient plays a very important role in the ability of the functional to describe this type of van der Waals attraction.
- 102Eshuis, H.; Furche, F. Basis set convergence of molecular correlation energy differences within the random phase approximation. J. Chem. Phys. 2012, 136, 084105, DOI: 10.1063/1.3687005Google Scholar102https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XivVeqtb0%253D&md5=d520f6c11aabd38dab2bc1d6b0ac6862Basis set convergence of molecular correlation energy differences within the random phase approximationEshuis, Henk; Furche, FilippJournal of Chemical Physics (2012), 136 (8), 084105/1-084105/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The basis set convergence of energy differences obtained from the RPA to the correlation energy is investigated for a wide range of mol. interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark. The use of very large basis sets (> quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable est. of the binding energy for these systems. Counterpoise cor. results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramol. dispersion, isomerization energies, and reaction energies of small org. mols. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency. (c) 2012 American Institute of Physics.
Cited By
This article has not yet been cited by other publications.
Article Views
Altmetric
Citations
Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.
Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.
The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.
Recommended Articles
Abstract
Figure 1
Figure 1. MAE (solid) and wMAE (striped) for the ASCDB database using different σ-functionals with PBE and PBE0 reference orbitals.
Figure 2
Figure 2. (a) MAE, (b) WTMAD-1, and (c) WTMAD-2 for the GMTKN55 database and its subcategories using different σ-functionals with PBE reference orbitals.
Figure 3
Figure 3. (a) MAE, (b) WTMAD-1, and (c) WTMAD-2 for the GMTKN55 database and its subcategories using different σ-functionals with PBE0 reference orbitals.
Figure 4
Figure 4. Error distribution for the W4–17RE benchmark set using the top-down σ↓AXK(A2) and bottom-up σ↑AXK(A2) functionals for (a) PBE and (b) PBE0 reference orbitals.
Figure 5
Figure 5. Weight functions of the Hartree (a, c) and exchange (b, d) kernels for σ↓AXK-functionals (a, b) and σ↑AXK-functionals (c, d) for PBE0 reference orbitals.
Figure 6
Figure 6. CBS-extrapolated neon dimer interaction energies computed using PBE reference orbitals for (a) KS-DFT, (b) σ-functionals and scaled σ-functionals, (c) top-down σ↓AXK-functionals, and (d) bottom-up σ↑AXK-functionals. For reference, corrected CCSD(T) values from ref (97) are also included; data between points has been interpolated using cubic splines.
References
This article references 102 other publications.
- 1Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864– B871, DOI: 10.1103/PhysRev.136.B864There is no corresponding record for this reference.
- 2Kohn, W.; Sham, L. J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133– A1138, DOI: 10.1103/PhysRev.140.A1133There is no corresponding record for this reference.
- 3Mardirossian, N.; Head-Gordon, M. Thirty years of density functional theory in computational chemistry: An overview and extensive assessment of 200 density functionals. Mol. Phys. 2017, 115, 2315– 2372, DOI: 10.1080/00268976.2017.13336443https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtVCltb3O&md5=ba27d707ee3f5fcdd949644d3d2cbd5eThirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionalsMardirossian, Narbe; Head-Gordon, MartinMolecular Physics (2017), 115 (19), 2315-2372CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)In the past 30 years, Kohn-Sham d. functional theory has emerged as the most popular electronic structure method in computational chem. To assess the ever-increasing no. of approx. exchange-correlation functionals, this review benchmarks a total of 200 d. functionals on a mol. database (MGCDB84) of nearly 5000 data points. The database employed, provided as Supplemental Data, is comprised of 84 data-sets and contains non-covalent interactions, isomerisation energies, thermochem., and barrier heights. In addn., the evolution of non-empirical and semi-empirical d. functional design is reviewed, and guidelines are provided for the proper and effective use of d. functionals. The most promising functional considered is ωB97M-V, a range-sepd. hybrid meta-GGA with VV10 nonlocal correlation, designed using a combinatorial approach. From the local GGAs, B97-D3, revPBE-D3, and BLYP-D3 are recommended, while from the local meta-GGAs, B97M-rV is the leading choice, followed by MS1-D3 and M06-L-D3. The best hybrid GGAs are ωB97X-V, ωB97X-D3, and ωB97X-D, while useful hybrid meta-GGAs (besides ωB97M-V) include ωM05-D, M06-2X-D3, and MN15. Ultimately, today's state-of-the-art functionals are close to achieving the level of accuracy desired for a broad range of chem. applications, and the principal remaining limitations are assocd. with systems that exhibit significant self-interaction/delocalisation errors and/or strong correlation effects.
- 4Perdew, J. P.; Schmidt, K. Jacob’s ladder of density functional approximations for the exchange-correlation energy. AIP Conf. Proc. 2001, 577, 1– 20, DOI: 10.1063/1.13901754https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXmtVGnsb0%253D&md5=cac0386cbbd9df5b30216875abcd1d24Jacob's ladder of density functional approximations for the exchange-correlation energyPerdew, John P.; Schmidt, KarlaAIP Conference Proceedings (2001), 577 (Density Functional Theory and Its Application to Materials), 1-20CODEN: APCPCS; ISSN:0094-243X. (American Institute of Physics)The ground-state energy and d. of a many-electron system are often calcd. by Kohn-Sham d. functional theory. We describe a ladder of approxns. for the exchange-correlation energy as a functional of the electron d. At the lowest rung of this ladder, the contribution to the energy from a vol. element of 3-dimensional space is detd. by the local d. there. Higher rungs or levels incorporate increasingly complex ingredients constructed from the d. or the Kohn-Sham orbitals in or around this vol. element. We identify which addnl. exact conditions can be satisfied at each level, and discuss the extent to which the functionals at each level may be constructed without empirical input. We also discuss the research that remains to be done at the exact-exchange level, and present our "dreams of a final theory". "Jacob left Beer-sheba and went toward Haran. He came to a certain place and stayed there for the night, because the sun had set. Taking one of the stones of the place, he put it under his head and lay down in that place. And he dreamed that there was a ladder set up on the earth, the top of it reaching to heaven; and the angels of God were ascending and descending on it.".
- 5Becke, A. D.; Santra, G.; Martin, J. M. L. A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database. J. Chem. Phys. 2023, 158, 151103, DOI: 10.1063/5.0141238There is no corresponding record for this reference.
- 6Kállay, M. Linear-scaling implementation of the direct random-phase approximation. J. Chem. Phys. 2015, 142, 204105, DOI: 10.1063/1.49215426https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXpt1Sjsb4%253D&md5=a66b1691368dd94596a8052c569de4d2Linear-scaling implementation of the direct random-phase approximationKallay, MihalyJournal of Chemical Physics (2015), 142 (20), 204105/1-204105/16CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We report the linear-scaling implementation of the direct RPA (dRPA) for closed-shell mol. systems. As a bonus, linear-scaling algorithms are also presented for the second-order screened exchange extension of dRPA as well as for the second-order Moller-Plesset (MP2) method and its spin-scaled variants. Our approach is based on an incremental scheme which is an extension of our previous local correlation method [Rolik et al., J. Chem. Phys. 139, 094105 (2013)]. The approach extensively uses local natural orbitals to reduce the size of the MO basis of local correlation domains. In addn., we also demonstrate that using natural auxiliary functions [M. Kallay, J. Chem. Phys. 141, 244113 (2014)], the size of the auxiliary basis of the domains and thus that of the three-center Coulomb integral lists can be reduced by an order of magnitude, which results in significant savings in computation time. The new approach is validated by extensive test calcns. for energies and energy differences. Our benchmark calcns. also demonstrate that the new method enables dRPA calcns. for mols. with more than 1000 atoms and 10 000 basis functions on a single processor. (c) 2015 American Institute of Physics.
- 7Schurkus, H. F.; Ochsenfeld, C. Communication: An effective linear-scaling atomic-orbital reformulation of the random-phase approximation using a contracted double-Laplace transformation. J. Chem. Phys. 2016, 144, 031101, DOI: 10.1063/1.49398417https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28Xht1Ciurw%253D&md5=b7f953b877f5fea31e2508e9f914be77Communication: An effective linear-scaling atomic-orbital reformulation of the random-phase approximation using a contracted double-Laplace transformationSchurkus, Henry F.; Ochsenfeld, ChristianJournal of Chemical Physics (2016), 144 (3), 031101/1-031101/5CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)An at.-orbital (AO) reformulation of the RPA correlation energy is presented allowing to reduce the steep computational scaling to linear, so that large systems can be studied on simple desktop computers with fully numerically controlled accuracy. Our AO-RPA formulation introduces a contracted double-Laplace transform and employs the overlap-metric resoln.-of-the-identity. First timings of our pilot code illustrate the reduced scaling with systems comprising up to 1262 atoms and 10 090 basis functions. (c) 2016 American Institute of Physics.
- 8Luenser, A.; Schurkus, H. F.; Ochsenfeld, C. Vanishing-Overhead Linear-Scaling Random Phase Approximation by Cholesky Decomposition and an Attenuated Coulomb-Metric. J. Chem. Theory Comput. 2017, 13, 1647, DOI: 10.1021/acs.jctc.6b012358https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXjs1yhsb0%253D&md5=a7a01e028f23cc0441ba2d752ebe200aVanishing-Overhead Linear-Scaling Random Phase Approximation by Cholesky Decomposition and an Attenuated Coulomb-MetricLuenser, Arne; Schurkus, Henry F.; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2017), 13 (4), 1647-1655CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A reformulation of the RPA within the Resoln.-of-the-Identity (RI) scheme is presented, that is competitive to canonical MO RI-RPA already for small- to medium-sized mols. For electronically sparse systems drastic speedups due to the reduced scaling behavior compared to the MO formulation are demonstrated. Our reformulation is based on two ideas which are independently useful: First, a Cholesky decompn. of d. matrixes that reduces the scaling with basis set size for a fixed-size mol. by one order, leading to massive performance improvements. Second, replacement of the overlap RI metric used in the original AO-RPA by an attenuated Coulomb metric. Accuracy is significantly improved compared to the overlap metric, while locality and sparsity of the integrals are retained, as is the effective linear scaling behavior.
- 9Graf, D.; Beuerle, M.; Schurkus, H. F.; Luenser, A.; Savasci, G.; Ochsenfeld, C. Accurate and Efficient Parallel Implementation of an Effective Linear-Scaling Direct Random Phase Approximation Method. J. Chem. Theory Comput. 2018, 14, 2505, DOI: 10.1021/acs.jctc.8b001779https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXnsFWisrg%253D&md5=b886c4aaf9c1b9e4a4f8482c453aa68bAccurate and Efficient Parallel Implementation of an Effective Linear-Scaling Direct Random Phase Approximation MethodGraf, Daniel; Beuerle, Matthias; Schurkus, Henry F.; Luenser, Arne; Savasci, Goekcen; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2018), 14 (5), 2505-2515CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)An efficient algorithm for calcg. the RPA correlation energy is presented that is as accurate as the canonical MO resoln.-of-the-identity RPA (RI-RPA) with the important advantage of an effective linear-scaling behavior (instead of quartic) for large systems due to a formulation in the local AO space. The high accuracy is achieved by utilizing optimized minimax integration schemes and the local Coulomb metric attenuated by the complementary error function for the RI approxn. The memory bottleneck of former AO (AO)-RI-RPA implementations is addressed by precontraction of the large 3-center integral matrix with the Cholesky factors of the ground state d. reducing the memory requirements of that matrix by a factor of (Nbasis/Nocc). Furthermore, we present a parallel implementation of our method, which not only leads to faster RPA correlation energy calcns. but also to a scalable decrease in memory requirements, opening the door for investigations of large mols. even on small- to medium-sized computing clusters. Although it is known that AO methods are highly efficient for extended systems, where sparsity allows for reaching the linear-scaling regime, we show that our work also extends the applicability when considering highly delocalized systems for which no linear scaling can be achieved. As an example, the interlayer distance of two covalent org. framework pore fragments (comprising 384 atoms in total) is analyzed.
- 10Beuerle, M.; Graf, D.; Schurkus, H. F.; Ochsenfeld, C. Efficient calculation of beyond RPA correlation energies in the dielectric matrix formalism. J. Chem. Phys. 2018, 148, 204104, DOI: 10.1063/1.502593810https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtVWlu7%252FM&md5=e7fa58c478a2a8532793dc946e50e57dEfficient calculation of beyond RPA correlation energies in the dielectric matrix formalismBeuerle, Matthias; Graf, Daniel; Schurkus, Henry F.; Ochsenfeld, ChristianJournal of Chemical Physics (2018), 148 (20), 204104/1-204104/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present efficient methods to calc. beyond RPA correlation energies for mol. systems with up to 500 atoms. To reduce the computational cost, we employ the resoln.-of-the-identity and a double-Laplace transform of the non-interacting polarization propagator in conjunction with an AO formalism. Further improvements are achieved using integral screening and the introduction of Cholesky decompd. densities. Our methods are applicable to the dielec. matrix formalism of RPA including second-order screened exchange (RPA-SOSEX), the RPA electron-hole time-dependent Hartree-Fock (RPA-eh-TDHF) approxn., and RPA renormalized perturbation theory using an approx. exchange kernel (RPA-AXK). We give an application of our methodol. by presenting RPA-SOSEX benchmark results for the L7 test set of large, dispersion dominated mols., yielding a mean abs. error below 1 kcal/mol. The present work enables calcg. beyond RPA correlation energies for significantly larger mols. than possible to date, thereby extending the applicability of these methods to a wider range of chem. systems. (c) 2018 American Institute of Physics.
- 11Drontschenko, V.; Graf, D.; Laqua, H.; Ochsenfeld, C. Lagrangian-Based Minimal-Overhead Batching Scheme for the Efficient Integral-Direct Evaluation of the RPA Correlation Energy. J. Chem. Theory Comput. 2021, 17, 5623, DOI: 10.1021/acs.jctc.1c0049411https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXhvVGqt7rN&md5=f4d0ae3663784bc63f97e43c80a919f8Lagrangian-Based Minimal-Overhead Batching Scheme for the Efficient Integral-Direct Evaluation of the RPA Correlation EnergyDrontschenko, Viktoria; Graf, Daniel; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2021), 17 (9), 5623-5634CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A highly memory-efficient integral-direct RPA method based on our ω-CDGD-RI-RPA method is presented that completely alleviates the memory bottleneck of storing the multidimensional three-center integral tensor, which severely limited the tractable system sizes. Based on a Lagrangian formulation, we introduce an optimized batching scheme over the auxiliary and basis-function indexes, which allows to compute the optimal no. of batches for a given amt. of system memory, while minimizing the batching overhead. Thus, our optimized batching constitutes the best tradeoff between program runtime and memory demand. Within this batching scheme, the half-transformed three-center integral tensor BiμM is recomputed for each batch of auxiliary and basis functions. This allows the computation of systems that were out of reach before. The largest system within this work consists of a DNA fragment comprising 1052 atoms and 11 230 basis functions calcd. on a single node, which emphasizes the new possibilities of our integral-direct RPA method.
- 12Harl, J.; Kresse, G. Cohesive energy curves for noble gas solids calculated by adiabatic connection fluctuation-dissipation theory. Phys. Rev. B 2008, 77, 045136, DOI: 10.1103/PhysRevB.77.04513612https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhslKrtLY%253D&md5=645e7ee23ef184cb921e0800974565fdCohesive energy curves for noble gas solids calculated by adiabatic connection fluctuation-dissipation theoryHarl, Judith; Kresse, GeorgPhysical Review B: Condensed Matter and Materials Physics (2008), 77 (4), 045136/1-045136/8CODEN: PRBMDO; ISSN:1098-0121. (American Physical Society)We present first-principles calcns. for the fcc noble gas solids Ne, Ar, and Kr applying the adiabatic connection fluctuation-dissipation theorem (ACFDT) to evaluate the correlation energy. The ACFDT allows us to describe long-range correlation effects including London dispersion or van der Waals interaction on top of conventional d. functional theory calcns. Even within the RPA, the typical 1/V2 vol. dependence for the cohesive energy of the noble gas solids is reproduced, and equil. cohesive energies and lattice consts. are improved compared to d. functional theory calcns. Furthermore, we present atomization energies for H2, N2, and O2 within the same post-d.-functional-theory framework, finding an excellent agreement with previously published data.
- 13Su, H.; Wu, Q.; Wang, H.; Wang, H. An assessment of the random-phase approximation functional and characteristics analysis for noncovalent cation−π interactions. Phys. Chem. Chem. Phys. 2017, 19, 26014– 26021, DOI: 10.1039/C7CP04504BThere is no corresponding record for this reference.
- 14Modrzejewski, M.; Yourdkhani, S.; Klimeš, J. Random Phase Approximation Applied to Many-Body Noncovalent Systems. J. Chem. Theory Comput. 2020, 16, 427– 442, DOI: 10.1021/acs.jctc.9b0097914https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3MfgtlWgtw%253D%253D&md5=257b7378ea6502c717458e2ebf398a11Random Phase Approximation Applied to Many-Body Noncovalent SystemsModrzejewski Marcin; Yourdkhani Sirous; Klimes Jiri; Modrzejewski MarcinJournal of chemical theory and computation (2020), 16 (1), 427-442 ISSN:.The random phase approximation (RPA) has received considerable interest in the field of modeling systems where noncovalent interactions are important. Its advantages over widely used density functional theory (DFT) approximations are the exact treatment of exchange and the description of long-range correlation. In this work, we address two open questions related to RPA. First, we demonstrate how accurately RPA describes nonadditive interactions encountered in many-body expansion of a binding energy. We consider three-body nonadditive energies in molecular and atomic clusters. Second, we address how the accuracy of RPA depends on input provided by different DFT models, without resorting to self-consistent RPA procedure, which is currently impractical for calculations employing periodic boundary conditions. We find that RPA based on the SCAN0 and PBE0 models, that is, hybrid DFT, achieves an overall accuracy between CCSD and MP3 on a data set of molecular trimers from Rezac et al. ( J. Chem. Theory. Comput. 2015 , 11 , 3065 ). Finally, many-body expansion for molecular clusters and solids often leads to a large number of small contributions that need to be calculated with high precision. We therefore present a cubic-scaling (or self-consistent field (SCF)-like) implementation of RPA in atomic basis set, which is designed for calculations with high numerical precision.
- 15Harl, J.; Kresse, G. Accurate Bulk Properties from Approximate Many-Body Techniques. Phys. Rev. Lett. 2009, 103, 056401, DOI: 10.1103/PhysRevLett.103.05640115https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXptl2jsbc%253D&md5=1a9a747d56ca24b35808a9aed92bb4dbAccurate Bulk Properties from Approximate Many-Body TechniquesHarl, Judith; Kresse, GeorgPhysical Review Letters (2009), 103 (5), 056401/1-056401/4CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)For ab initio electronic structure calcns., the RPA to the correlation energy is supposed to be a suitable complement to the exact exchange energy. We show that lattice consts., atomization energies of solids, and adsorption energies on metal surfaces evaluated using this approxn. are in very good agreement with expt. Since the method is fairly efficient and handles ionic, metallic, and van der Waals bonded systems equally well, it is a very promising choice to improve upon d. functional theory calcns., without resorting to more demanding diffusion Monte Carlo or quantum chem. methods.
- 16Kreppel, A.; Graf, D.; Laqua, H.; Ochsenfeld, C. Range-Separated Density-Functional Theory in Combination with the Random Phase Approximation: An Accuracy Benchmark. J. Chem. Theory Comput. 2020, 16, 2985– 2994, DOI: 10.1021/acs.jctc.9b0129416https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXnslCit7c%253D&md5=e14a638f541b4fb9a2541961f2062ccbRange-Separated Density-Functional Theory in Combination with the Random Phase Approximation: An Accuracy BenchmarkKreppel, Andrea; Graf, Daniel; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2020), 16 (5), 2985-2994CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A formulation of range-sepd. RPA based on our efficient ω-CDGD-RI-RPA [J. Chem. Theory Comput.2018, 14, 2505] method and a large scale benchmark study are presented. By application to the GMTKN55 data set, we obtain a comprehensive picture of the performance of range-sepd. RPA in general main group thermochem., kinetics, and noncovalent interactions. The results show that range-sepd. RPA performs stably over the broad range of mol. chem. included in the GMTKN55 set. It improves significantly over semilocal DFT but it is still less accurate than modern dispersion cor. double-hybrid functionals. Furthermore, range-sepd. RPA shows a faster basis set convergence compared to std. full-range RPA making it a promising applicable approach with only one empirical parameter.
- 17Trushin, E.; Thierbach, A.; Görling, A. Toward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energy. J. Chem. Phys. 2021, 154, 014104, DOI: 10.1063/5.002684917https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXosFWrtg%253D%253D&md5=d3d9596431d9804901d50736957d01abToward chemical accuracy at low computational cost: Density-functional theory with σ-functionals for the correlation energyTrushin, Egor; Thierbach, Adrian; Goerling, AndreasJournal of Chemical Physics (2021), 154 (1), 014104CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We introduce new functionals for the Kohn-Sham correlation energy that are based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem and are named σ-functionals. Like in the well-established direct RPA (dRPA), σ-functionals require as input exclusively eigenvalues σ of the frequency-dependent KS response function. In the new functionals, functions of σ replace the σ-dependent dRPA expression in the coupling-const. and frequency integrations contained in the ACFD theorem. We optimize σ-functionals with the help of ref. sets for atomization, reaction, transition state, and non-covalent interaction energies. The optimized functionals are to be used in a post-self-consistent way using orbitals and eigenvalues from conventional Kohn-Sham calcns. employing the exchange-correlation functional of Perdew, Burke, and Ernzerhof. The accuracy of the presented approach is much higher than that of dRPA methods and is comparable to that of high-level wave function methods. Reaction and transition state energies from σ-functionals exhibit accuracies close to 1 kcal/mol and thus approach chem. accuracy. For the 10 966 reactions of the W4-11RE ref. set, the mean abs. deviation is 1.25 kcal/mol compared to 3.21 kcal/mol in the dRPA case. Non-covalent binding energies are accurate to a few tenths of a kcal/mol. The presented approach is highly efficient, and the post-self-consistent calcn. of the total energy requires less computational time than a d.-functional calcn. with a hybrid functional and thus can be easily carried out routinely. σ-Functionals can be implemented in any existing dRPA code with negligible programming effort. (c) 2021 American Institute of Physics.
- 18Fauser, S.; Trushin, E.; Neiss, C.; Görling, A. Chemical accuracy with σ-functionals for the Kohn–Sham correlation energy optimized for different input orbitals and eigenvalues. J. Chem. Phys. 2021, 155, 134111, DOI: 10.1063/5.005964118https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXit1Sqtr3N&md5=7c717988178d285b3f1c00b20c750598Chemical accuracy with σ-functionals for the Kohn-Sham correlation energy optimized for different input orbitals and eigenvaluesFauser, Steffen; Trushin, Egor; Neiss, Christian; Goerling, AndreasJournal of Chemical Physics (2021), 155 (13), 134111CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recently, a new type of orbital-dependent functional for the Kohn-Sham (KS) correlation energy, σ-functionals, was introduced. Tech., σ-functionals are closely related to the well-known direct RPA (dRPA). Within the dRPA, a function of the eigenvalues σ of the frequency-dependent KS response function is integrated over purely imaginary frequencies. In σ-functionals, this function is replaced by one that is optimized with respect to ref. sets of atomization, reaction, transition state, and non-covalent interaction energies. The previously introduced σ-functional uses input orbitals and eigenvalues from KS calcns. with the generalized gradient approxn. (GGA) exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE). Here, σ-functionals using input orbitals and eigenvalues from the meta-GGA TPSS and the hybrid-functionals PBE0 and B3LYP are presented and tested. The no. of ref. sets taken into account in the optimization of the σ-functionals is larger than in the first PBE based σ-functional and includes sets with 3d-transition metal compds. Therefore, also a reparameterized PBE based σ-functional is introduced. The σ-functionals based on PBE0 and B3LYP orbitals and eigenvalues reach chem. accuracy for main group chem. For the 10 966 reactions from the highly accurate W4-11RE ref. set, the B3LYP based σ-functional exhibits a mean av. deviation of 1.03 kcal/mol compared to 1.08 kcal/mol for the coupled cluster singles doubles perturbative triples method if the same valence quadruple zeta basis set is used. For 3d-transition metal chem., accuracies of about 2 kcal/mol are reached. The computational effort for the post-self-consistent evaluation of the σ-functional is lower than that of a preceding PBE0 or B3LYP calcn. for typical systems. (c) 2021 American Institute of Physics.
- 19Fauser, S.; Förster, A.; Redeker, L.; Neiss, C.; Erhard, J.; Trushin, E.; Görling, A. Basis Set Requirements of σ-Functionals for Gaussian- and Slater-Type Basis Functions and Comparison with Range-Separated Hybrid and Double Hybrid Functionals. J. Chem. Theory Comput. 2024, 20, 2404– 2422, DOI: 10.1021/acs.jctc.3c01132There is no corresponding record for this reference.
- 20Drontschenko, V.; Graf, D.; Laqua, H.; Ochsenfeld, C. Efficient Method for the Computation of Frozen-Core Nuclear Gradients within the Random Phase Approximation. J. Chem. Theory Comput. 2022, 18, 7359– 7372, DOI: 10.1021/acs.jctc.2c0077420https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB38Xisl2gtrzP&md5=01af0b01ddae622b53bd3164b285a9d9Efficient Method for the Computation of Frozen-Core Nuclear Gradients within the Random Phase ApproximationDrontschenko, Viktoria; Graf, Daniel; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2022), 18 (12), 7359-7372CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A method for the evaluation of anal. frozen-core gradients within the RPA is presented. We outline an efficient way to evaluate the response of the d. of active electrons arising only when introducing the frozen-core approxn. and constituting the main difficulty, together with the response of the std. Kohn-Sham d. The general framework allows to extend the outlined procedure to related electron correlation methods in the AO basis that require the evaluation of d. responses, such as second-order Moller-Plesset perturbation theory or coupled cluster variants. By using Cholesky decompd. densities - which reintroduce the occupied index in the time-detg. steps - we are able to achieve speedups of 20-30% (depending on the size of the basis set) by using the frozen-core approxn., which is of similar magnitude as for MO formulations. We further show that the errors introduced by the frozen-core approxn. are practically insignificant for mol. geometries.
- 21Neiss, C.; Fauser, S.; Görling, A. Geometries and vibrational frequencies with Kohn–Sham methods using σ-functionals for the correlation energy. J. Chem. Phys. 2023, 158, 044107, DOI: 10.1063/5.012952421https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3sXhvVGqsL4%253D&md5=f80b8233816be1e4e3ee4925fc0a5357Geometries and vibrational frequencies with Kohn-Sham methods using σ-functionals for the correlation energyNeiss, Christian; Fauser, Steffen; Goerling, AndreasJournal of Chemical Physics (2023), 158 (4), 044107CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Recently, Kohn-Sham (KS) methods with new correlation functionals, called σ-functionals, have been introduced. Tech., σ-functionals are closely related to the well-known RPA; formally, σ-functionals are rooted in perturbation theory along the adiabatic connection. If employed in a post-SCF manner in a Gaussian basis set framework, then, σ-functional methods are computationally very efficient. Moreover, for main group chem., σ-functionals are highly accurate and can compete with high-level wave-function methods. For reaction and transition state energies, e.g., chem. accuracy of 1 kcal/mol is reached. Here, we show how to calc. first derivs. of the total energy with respect to nuclear coordinates for methods using σ-functionals and then carry out geometry optimizations for test sets of main group mols., transition metal compds., and non-covalently bonded systems. For main group mols., we addnl. calc. vibrational frequencies. σ-Functional methods are found to yield very accurate geometries and vibrational frequencies for main group mols. superior not only to those from conventional KS methods but also to those from RPA methods. For geometries of transition metal compds., not surprisingly, best geometries are found for RPA methods, while σ-functional methods yield somewhat less good results. This is attributed to the fact that in the optimization of σ-functionals, transition metal compds. could not be represented well due to the lack of reliable ref. data. For non-covalently bonded systems, σ-functionals yield geometries of the same quality as the RPA or as conventional KS schemes combined with dispersion corrections. (c) 2023 American Institute of Physics.
- 22Drontschenko, V.; Bangerter, F. H.; Ochsenfeld, C. Analytical Second-Order Properties for the Random Phase Approximation: Nuclear Magnetic Resonance Shieldings. J. Chem. Theory Comput. 2023, 19, 7542– 7554, DOI: 10.1021/acs.jctc.3c00542There is no corresponding record for this reference.
- 23Fauser, S.; Drontschenko, V.; Ochsenfeld, C.; Görling, A. Accurate NMR Shieldings with σ-Functionals. J. Chem. Theory Comput. 2024, 20, 6028– 6036, DOI: 10.1021/acs.jctc.4c00512There is no corresponding record for this reference.
- 24Drontschenko, V.; Ochsenfeld, C. Low-Scaling, Efficient and Memory Optimized Computation of Nuclear Magnetic Resonance Shieldings within the Random Phase Approximation using Cholesky-Decomposed Densities and an Attenuated Coulomb Metric. J. Phys. Chem. A 2024, 128, 7950– 7965, DOI: 10.1021/acs.jpca.4c02773There is no corresponding record for this reference.
- 25Glasbrenner, M.; Graf, D.; Ochsenfeld, C. Benchmarking the Accuracy of the Direct Random Phase Approximation and σ-Functionals for NMR Shieldings. J. Chem. Theory Comput. 2022, 18, 192– 205, DOI: 10.1021/acs.jctc.1c0086625https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXislWrtbvF&md5=25dffbe9eb6cff42e76bdf9a5d151a2eBenchmarking the Accuracy of the Direct Random Phase Approximation and σ-Functionals for NMR ShieldingsGlasbrenner, Michael; Graf, Daniel; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2022), 18 (1), 192-205CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A method for computing NMR shieldings with the direct RPA (RPA) and the closely related σ-functionals [Trushin, E.; Thierbach, A.; Gorling, A. Toward chem. accuracy at low computational cost: d. functional theory with σ-functionals for the correlation energy. J. Chem. Phys.2021,154, 014104] is presented, which is based on a finite-difference approach. The accuracy is evaluated in benchmark calcns. using high-quality coupled cluster values as a ref. Our results show that the accuracy of the computed NMR shieldings using direct RPA is strongly dependent on the d. functional theory ref. orbitals and improves with increasing amts. of exact Hartree-Fock exchange in the functional. NMR shieldings computed with the direct RPA with a Hartree-Fock ref. are significantly more accurate than MP2 shieldings and comparable to CCSD shieldings. Also, the basis set convergence is analyzed and it is shown that at least triple-zeta basis sets are required for reliable results.
- 26Görling, A. Hierarchies of methods towards the exact Kohn–Sham correlation energy based on the adiabatic-connection fluctuation-dissipation theorem. Phys. Rev. B 2019, 99, 235120, DOI: 10.1103/PhysRevB.99.235120There is no corresponding record for this reference.
- 27Erhard, J.; Bleiziffer, P.; Görling, A. Power Series Approximation for the Correlation Kernel Leading to Kohn–Sham Methods Combining Accuracy, Computational Efficiency, and General Applicability. Phys. Rev. Lett. 2016, 117, 143002, DOI: 10.1103/PhysRevLett.117.143002There is no corresponding record for this reference.
- 28Erhard, J.; Fauser, S.; Trushin, E.; Görling, A. Scaled σ-functionals for the Kohn–Sham correlation energy with scaling functions from the homogeneous electron gas. J. Chem. Phys. 2022, 157, 114105, DOI: 10.1063/5.010164128https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB38XisVWnsb%252FJ&md5=c3a5d9d86b9d5d799fbd3366b86f21d2Scaled σ-functionals for Kohn-Sham correlation energy with scaling function from homogeneous electron gasErhard, Jannis; Fauser, Steffen; Trushin, Egor; Goerling, AndreasJournal of Chemical Physics (2022), 157 (11), 114105CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. The σ-functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct RPA (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigen values σ of the Kohn-Sham response matrix that enters the coupling const. and frequency integration in the dRPA is replaced by another function optimized with the help of ref. sets of atomization, reaction, transition state, and non-covalent interaction energies, and σ-Functionals are highly accurate and yield chem. accuracy of 1 kcal/mol in reaction or transition state energies, in main group chem. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron no., such as ionization or electron attachment. This problem is attributed to unphys. self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are detected with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addn. yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts. (c) 2022 American Institute of Physics.
- 29Lemke, Y.; Ochsenfeld, C. Highly accurate σ- and τ-functionals for beyond-RPA methods with approximate exchange kernels. J. Chem. Phys. 2023, 159, 194104, DOI: 10.1063/5.017304229https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3sXitlOgtr3N&md5=7cb6539f08f4d9e17b24a18fb1bfda82Highly accurate σ- and τ-functionals for beyond-RPA methods with approximate exchange kernelsLemke, Yannick; Ochsenfeld, ChristianJournal of Chemical Physics (2023), 159 (19), 194104CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)σ-Functionals are promising new developments for the Kohn-Sham correlation energy based upon the direct RPA (dRPA) within the adiabatic connection formalism, providing impressive improvements over dRPA for a broad range of benchmarks. However, σ-functionals exhibit a high amt. of self-interaction inherited from the approxns. made within dRPA. Inclusion of an exchange kernel in deriving the coupling-strength-dependent d.-d. response function leads to so-called τ-functionals, which - apart from a fourth-order Taylor series expansion - have only been realized in an approx. fashion so far to the best of our knowledge, most notably in the form of scaled σ-functionals. In this work, we derive, optimize, and benchmark three types of σ- and τ-functionals including approx. exchange effects in the form of an antisymmetrized Hartree kernel. These functionals, based on a second-order screened exchange type contribution in the adiabatic connection formalism, the electron-hole time-dependent Hartree-Fock kernel (eh-TDHF) otherwise known as RPA with exchange (RPAx), and an approxn. thereof known as approx. exchange kernel (AXK), are optimized on the ASCDB database using two new parametrizations named A1 and A2. In addn., we report a first full evaluation of σ- and τ-functionals on the GMTKN55 database, revealing our exchange-including functionals to considerably outperform existing σ-functionals while being highly competitive with some of the best double-hybrid functionals of the original GMTKN55 publication. In particular, the σ-functionals based on AXK and τ-functionals based on RPAx with PBE0 ref. stand out as highly accurate approaches for a wide variety of chem. relevant problems. (c) 2023 American Institute of Physics.
- 30Grüneis, A.; Marsman, M.; Harl, J.; Schimka, L.; Kresse, G. Making the random phase approximation to electronic correlation accurate. J. Chem. Phys. 2009, 131, 154115, DOI: 10.1063/1.325034730https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhsVSjs7fJ&md5=d186251717a9fe10e70b360c9f350c53Making the random phase approximation to electronic correlation accurateGrueneis, Andreas; Marsman, Martijn; Harl, Judith; Schimka, Laurids; Kresse, GeorgJournal of Chemical Physics (2009), 131 (15), 154115/1-154115/5CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We show that the inclusion of second-order screened exchange to the RPA allows for an accurate description of electronic correlation in atoms and solids clearly surpassing the random phase approxn., but not yet approaching chem. accuracy. From a fundamental point of view, the method is self-correlation free for one-electron systems. From a practical point of view,the approach yields correlation energies for atoms, as well as for the jellium electron gas within a few kcal/mol of exact values, atomization energies within typically 2-3 kcal/mol of expt., and excellent lattice consts. for ionic and covalently bonded solids (0.2% error). The computational complexity is only O(N5), comparable to canonical second-order Moller-Plesset perturbation theory, which should allow for routine calcns. on many systems. (c) 2009 American Institute of Physics.
- 31Bates, J. E.; Furche, F. Communication: Random phase approximation renormalized many-body perturbation theory. J. Chem. Phys. 2013, 139, 171103, DOI: 10.1063/1.482725431https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhsleiu73P&md5=bbfd81d02296e4acbc349e365c8b9b70Communication: Random phase approximation renormalized many-body perturbation theoryBates, Jefferson E.; Furche, FilippJournal of Chemical Physics (2013), 139 (17), 171103/1-171103/4CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We derive a renormalized many-body perturbation theory (MBPT) starting from the RPA. This RPA-renormalized perturbation theory extends the scope of single-ref. MBPT methods to small-gap systems without significantly increasing the computational cost. The leading correction to RPA, termed the approx. exchange kernel (AXK), substantially improves upon RPA atomization energies and ionization potentials without affecting other properties such as barrier heights where RPA is already accurate. Thus, AXK is more balanced than second-order screened exchange, which tends to overcorrect RPA for systems with stronger static correlation. Similarly, AXK avoids the divergence of second-order Moller-Plesset (MP2) theory for small gap systems and delivers a much more consistent performance than MP2 across the periodic table at comparable cost. RPA+AXK thus is an accurate, non-empirical, and robust tool to assess and improve semi-local d. functional theory for a wide range of systems previously inaccessible to first-principles electronic structure calcns. (c) 2013 American Institute of Physics.
- 32Goerigk, L.; Hansen, A.; Bauer, C.; Ehrlich, S.; Najibi, A.; Grimme, S. A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions. Phys. Chem. Chem. Phys. 2017, 19, 32184, DOI: 10.1039/C7CP04913G32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhslajtLnF&md5=f9393c9e3907336c4da053743797f8dfA look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactionsGoerigk, Lars; Hansen, Andreas; Bauer, Christoph; Ehrlich, Stephan; Najibi, Asim; Grimme, StefanPhysical Chemistry Chemical Physics (2017), 19 (48), 32184-32215CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)We present the GMTKN55 benchmark database for general main group thermochem., kinetics and noncovalent interactions. Compared to its popular predecessor GMTKN30, it allows assessment across a larger variety of chem. problems - with 13 new benchmark sets being presented for the first time - and it also provides ref. values of significantly higher quality for most sets. GMTKN55 comprises 1505 relative energies based on 2462 single-point calcns. and it is accessible to the user community via a dedicated website. Herein, we demonstrate the importance of better ref. values, and we re-emphasize the need for London-dispersion corrections in d. functional theory (DFT) treatments of thermochem. problems, including Minnesota methods. We assessed 217 variations of dispersion-cor. and -uncorrected d. functional approxns., and carried out a detailed anal. of 83 of them to identify robust and reliable approaches. Double-hybrid functionals are the most reliable approaches for thermochem. and noncovalent interactions, and they should be used whenever tech. feasible. These are, in particular, DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), and B2GPPLYP-D3(BJ). The best hybrids are ωB97X-V, M052X-D3(0), and ωB97X-D3, but we also recommend PW6B95-D3(BJ) as the best conventional global hybrid. At the meta-generalized-gradient (meta-GGA) level, the SCAN-D3(BJ) method can be recommended. Other meta-GGAs are outperformed by the GGA functionals revPBE-D3(BJ), B97-D3(BJ), and OLYP-D3(BJ). We note that many popular methods, such as B3LYP, are not part of our recommendations. In fact, with our results we hope to inspire a change in the user community's perception of common DFT methods. We also encourage method developers to use GMTKN55 for cross-validation studies of new methodologies.
- 33Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865– 3868, DOI: 10.1103/PhysRevLett.77.386533https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XmsVCgsbs%253D&md5=55943538406ee74f93aabdf882cd4630Generalized gradient approximation made simplePerdew, John P.; Burke, Kieron; Ernzerhof, MatthiasPhysical Review Letters (1996), 77 (18), 3865-3868CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)Generalized gradient approxns. (GGA's) for the exchange-correlation energy improve upon the local spin d. (LSD) description of atoms, mols., and solids. We present a simple derivation of a simple GGA, in which all parameters (other than those in LSD) are fundamental consts. Only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked. Improvements over PW91 include an accurate description of the linear response of the uniform electron gas, correct behavior under uniform scaling, and a smoother potential.
- 34Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158– 6170, DOI: 10.1063/1.47852234https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXitVCmt7Y%253D&md5=cad4185c69f9232753497f5203d6dc9fToward reliable density functional methods without adjustable parameters: the PBE0 modelAdamo, Carlo; Barone, VincenzoJournal of Chemical Physics (1999), 110 (13), 6158-6170CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present an anal. of the performances of a parameter free d. functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amt. of exact exchange. The results obtained for structural, thermodn., kinetic and spectroscopic (magnetic, IR and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chem. and condensed matter physics.
- 35Ernzerhof, M.; Scuseria, G. E. Assessment of the Perdew–Burke–Ernzerhof exchange-correlation functional. J. Chem. Phys. 1999, 110, 5029– 5036, DOI: 10.1063/1.47840135https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1MXhsFWquro%253D&md5=6113a915baa66b0f9961a5e160cadaf1Assessment of the Perdew-Burke-Ernzerhof exchange-correlation functionalErnzerhof, Matthias; Scuseria, Gustavo E.Journal of Chemical Physics (1999), 110 (11), 5029-5036CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)In order to discriminate between approxns. to the exchange-correlation energy EXC[ρ↑,ρ↓], we employ the criterion of whether the functional is fitted to a certain exptl. data set or if it is constructed to satisfy phys. constraints. We present extensive test calcns. for atoms and mols., with the nonempirical local spin-d. (LSD) and the Perdew-Burke-Ernzerhof (PBE) functional and compare our results with results obtained with more empirical functionals. For the atomization energies of the G2 set, we find that the PBE functional shows systematic errors larger than those of commonly used empirical functionals. The PBE ionization potentials, electron affinities, and bond lengths are of accuracy similar to those obtained from empirical functionals. Furthermore, a recently proposed hybrid scheme using exact exchange together with PBE exchange and correlation is investigated. For all properties studied here, the PBE hybrid gives an accuracy comparable to the frequently used empirical B3LYP hybrid scheme. Phys. principles underlying the PBE and PBE hybrid scheme are examd. and the range of their validity is discussed.
- 36Langreth, D. C.; Perdew, J. P. The exchange-correlation energy of a metallic surface. Solid State Commun. 1975, 17, 1425, DOI: 10.1016/0038-1098(75)90618-3There is no corresponding record for this reference.
- 37Langreth, D. C.; Perdew, J. P. Exchange-correlation energy of a metallic surface: Wave-vector analysis. Phys. Rev. B 1977, 15, 2884, DOI: 10.1103/PhysRevB.15.2884There is no corresponding record for this reference.
- 38Petersilka, M.; Gossmann, U. J.; Gross, E. K. U. Excitation Energies from Time-Dependent Density-Functional Theory. Phys. Rev. Lett. 1996, 76, 1212– 1215, DOI: 10.1103/PhysRevLett.76.121238https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XhtFGjt7k%253D&md5=67ce109fced1909cc6dabf2f9aab0116Excitation energies from time-dependent density-functional theoryPetersilka, M.; Gossmann, U. J.; Gross, E. K. U.Physical Review Letters (1996), 76 (8), 1212-15CODEN: PRLTAO; ISSN:0031-9007. (American Physical Society)A new d.-functional approach to calc. the excitation spectrum of many-electron systems is proposed. It is shown that the full linear d. response of the interacting system, which has poles at the exact excitation energies, can rigorously be expressed in terms of the response function of the noninteracting (Kohn-Sham) system and a frequency-dependent exchange-correlation kernel. Using this expression, the poles of the full response function are obtained by systematic improvement upon the poles of the Kohn-Sham response function. Numerical results are presented for Be, Mg, Ca, Zn, Sr, and Cd atoms.
- 39Garrick, R.; Natan, A.; Gould, T.; Kronik, L. Exact Generalized Kohn–Sham Theory for Hybrid Functionals. Phys. Rev. X 2020, 10, 021040, DOI: 10.1103/PhysRevX.10.021040There is no corresponding record for this reference.
- 40Whitten, J. L. Coulombic potential energy integrals and approximations. J. Chem. Phys. 1973, 58, 4496– 4501, DOI: 10.1063/1.167901240https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE3sXktFSjtbo%253D&md5=05c510c8db660386b6fd3df789c22b10Coulombic potential energy integrals and approximationsWhitten, J. L.Journal of Chemical Physics (1973), 58 (10), 4496-501CODEN: JCPSA6; ISSN:0021-9606.Theorems are derived which establish a method of approxg. 2-particle Coulombic potential energy integrals, [.vphi.a(1)|r12-1|.vphi.b-(2)], in terms of approx. charge ds. .vphi.a' and .vphi.b'. Rigorous error bounds, |[.vphi.a(1)|r12-1|.vphi.b(2)] - [.vphi.a'(1)|r12-1|.vphi.b'(2)]| ≤ δ, are simply expressed in terms of information calcd. sep. for the pair of ds. .vphi.a and .vphi.b' and the pair .vphi.b and .vphi.b'. From the structure of the bound, a simple method of optimizing charge d. approxns. such that δ is minimized is derived. The framework of the theory appears to be well suited for application to the approxn. of electron repulsion integrals which occur in mol. structure theory, and applications to the approxn. of integrals over Slater orbitals or grouped Gaussian functions are discussed.
- 41Dunlap, B. I.; Connolly, J. W. D.; Sabin, J. R. On some approximations in applications of Xα theory. J. Chem. Phys. 1979, 71, 3396– 3402, DOI: 10.1063/1.43872841https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaE1MXlvF2hs7k%253D&md5=0d61bda3e94a1ef30d221972f51ea6cfOn some approximations in applications of X α theoryDunlap, B. I.; Connolly, J. W. D.; Sabin, J. R.Journal of Chemical Physics (1979), 71 (8), 3396-402CODEN: JCPSA6; ISSN:0021-9606.An approx. Xα functional is proposed from which the charge d. fitting equations follow variationally. LCAO Xα calcns. on at. Ni and H2 show the method independent of the fitting (auxiliary) bases to within 0.02 eV. Variational properties assocd. with both orbital and auxiliary basis set incompleteness are used to approach within 0.2 eV the Xα total energy limit for the N mol.
- 42Vahtras, O.; Almlöf, J.; Feyereisen, M. W. Integral approximations for LCAO-SCF calculations. Chem. Phys. Lett. 1993, 213, 514– 518, DOI: 10.1016/0009-2614(93)89151-742https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3sXmsFKksL4%253D&md5=2c0f258baf7af4b891974d9b71f3bbb2Integral approximations for LCAO-SCF calculationsVahtras, O.; Almloef, J.; Feyereisen, M. W.Chemical Physics Letters (1993), 213 (5-6), 514-18CODEN: CHPLBC; ISSN:0009-2614.Three-center approxns. to the four-center integrals occurring in ab initio LCAO calcns. are investigated. Significant gains in computer time can be obtained without sacrificing accuracy, if a suitable expansion basis is chosen.
- 43Furche, F. Developing the random phase approximation into a practical post-Kohn–Sham correlation model. J. Chem. Phys. 2008, 129, 114105, DOI: 10.1063/1.297778943https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXhtFOrtbbP&md5=b4d01f7faebebbfbbbf6ebc9c0077102Developing the random phase approximation into a practical post-Kohn-Sham correlation modelFurche, FilippJournal of Chemical Physics (2008), 129 (11), 114105/1-114105/8CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The RPA (RPA) to the d. functional correlation energy systematically improves upon many limitations of present semilocal functionals, but was considered too computationally expensive for widespread use in the past. Here a phys. appealing reformulation of the RPA correlation model is developed that substantially reduces its computational complexity. The d. functional RPA correlation energy is shown to equal one-half times the difference of all RPA electronic excitation energies computed at full and first order coupling. Thus, the RPA correlation energy may be considered as a difference of electronic zero point vibrational energies, where each eigenmode corresponds to an electronic excitation. This surprisingly simple result is intimately related to plasma theories of electron correlation. Differences to electron pair correlation models underlying popular correlated wave function methods are discussed. The RPA correlation energy is further transformed into an explicit functional of the Kohn-Sham orbitals. The only nontrivial ingredient to this functional is the sign function of the response operator. A stable iterative algorithm to evaluate this sign function based on the Newton-Schulz iteration is presented. Integral direct implementations scale as the fifth power of the system size, similar to second order Moeller-Plesset calcns. With these improvements, RPA may become the long-sought robust and efficient zero order post-Kohn-Sham correlation model. (c) 2008 American Institute of Physics.
- 44Eshuis, H.; Yarkony, J.; Furche, F. Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integration. J. Chem. Phys. 2010, 132, 234114, DOI: 10.1063/1.344274944https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXnslSlsLw%253D&md5=90e538090088f1ded996c402ffc7b342Fast computation of molecular random phase approximation correlation energies using resolution of the identity and imaginary frequency integrationEshuis, Henk; Yarkony, Julian; Furche, FilippJournal of Chemical Physics (2010), 132 (23), 234114/1-234114/9CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The RPA is an increasingly popular post-Kohn-Sham correlation method, but its high computational cost has limited mol. applications to systems with few atoms. Here we present an efficient implementation of RPA correlation energies based on a combination of resoln. of the identity (RI) and imaginary frequency integration techniques. We show that the RI approxn. to four-index electron repulsion integrals leads to a variational upper bound to the exact RPA correlation energy if the Coulomb metric is used. Auxiliary basis sets optimized for second-order Moller-Plesset (MP2) calcns. are well suitable for RPA, as is demonstrated for the HEAT and MOLEKEL benchmark sets. Using imaginary frequency integration rather than diagonalization to compute the matrix square root necessary for RPA, evaluation of the RPA correlation energy requires O(N4logN) operations and O(N3) storage only; the price for this dramatic improvement over existing algorithms is a numerical quadrature. We propose a numerical integration scheme that is exact in the two-orbital case and converges exponentially with the no. of grid points. For most systems, 30-40 grid points yield μH accuracy in triple zeta basis sets, but much larger grids are necessary for small gap systems. The lowest-order approxn. to the present method is a post-Kohn-Sham frequency-domain version of opposite-spin Laplace-transform RI-MP2. Timings for polyacenes with up to 30 atoms show speed-ups of two orders of magnitude over previous implementations. The present approach makes it possible to routinely compute RPA correlation energies of systems well beyond 100 atoms, as is demonstrated for the octapeptide angiotensin II. (c) 2010 American Institute of Physics.
- 45Bleiziffer, P.; Heßelmann, A.; Görling, A. Resolution of identity approach for the Kohn–Sham correlation energy within the exact-exchange random-phase approximation. J. Chem. Phys. 2012, 136, 134102, DOI: 10.1063/1.3697845There is no corresponding record for this reference.
- 46Lemke, Y.; Graf, D.; Kussmann, J.; Ochsenfeld, C. An assessment of orbital energy corrections for the direct random phase approximation and explicit σ-functionals. Mol. Phys. 2023, 121, e2098862, DOI: 10.1080/00268976.2022.2098862There is no corresponding record for this reference.
- 47Heßelmann, A.; Görling, A. Random phase approximation correlation energies with exact Kohn–Sham exchange. Mol. Phys. 2010, 108, 359– 372, DOI: 10.1080/0026897090347666247https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXjvFKksLw%253D&md5=099fc620e45aaffe5161baea6e7454acRandom phase approximation correlation energies with exact Kohn-Sham exchangeHesselmann, Andreas; Goerling, AndreasMolecular Physics (2010), 108 (3-4), 359-372CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)The RPA correlation energy is expressed in terms of the exact local Kohn-Sham (KS) exchange potential and corresponding adiabatic and nonadiabatic exchange kernels for d.-functional ref. determinants. The approach naturally extends the RPA method in which, conventionally, only the Coulomb kernel is included. By comparison with the coupled cluster singles doubles with perturbative triples method it is shown for a set of small mols. that the new RPA method based on KS exchange yields correlation energies more accurate than RPA on the basis of Hartree-Fock exchange.
- 48Ángyán, J. G.; Liu, R.-F.; Toulouse, J.; Jansen, G. Correlation Energy Expressions from the Adiabatic-Connection Fluctuation–Dissipation Theorem Approach. J. Chem. Theory Comput. 2011, 7, 3116– 3130, DOI: 10.1021/ct200501r48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXht1Sltb3P&md5=03529fc836d25b236e71c89e4170f1adCorrelation Energy Expressions from the Adiabatic-Connection Fluctuation-Dissipation Theorem ApproachAngyan, Janos G.; Liu, Ru-Fen; Toulouse, Julien; Jansen, GeorgJournal of Chemical Theory and Computation (2011), 7 (10), 3116-3130CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We explore several RPA correlation energy variants within the adiabatic-connection fluctuation-dissipation theorem approach. These variants differ in the way the exchange interactions are treated. One of these variants, named dRPA-II, is original to this work and closely resembles the second-order screened exchange (SOSEX) method. We discuss and clarify the connections among different RPA formulations. We derive the spin-adapted forms of all the variants for closed-shell systems and test them on a few at. and mol. systems with and without range sepn. of the electron-electron interaction.
- 49Heßelmann, A. Random-phase-approximation correlation method including exchange interactions. Phys. Rev. A 2012, 85, 012517, DOI: 10.1103/PhysRevA.85.01251749https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XjtFahu7g%253D&md5=1b15a12a2c98b16c7db2d9a9ced0d34eRandom-phase-approximation correlation method including exchange interactionsHesselmann, AndreasPhysical Review A: Atomic, Molecular, and Optical Physics (2012), 85 (1-A), 012517/1-012517/10CODEN: PLRAAN; ISSN:1050-2947. (American Physical Society)Two random-phase-approxn. correlation methods are introduced that take into account exchange interactions. The first one, termed RPAX, is obtained from a simple modification of the ring coupled-cluster doubles amplitude equation, while the second, termed RPAX2, is based on the first method using a slightly modified update equation for the amplitudes. It is shown that this second RPAX2 method can be implemented with a computational algorithm that scales only with the fifth power of the mol. size with the aid of d. fitting or the Cholesky decompn. of two-electron integrals. It is thus not much more costly than std. second-order perturbation theory methods and can be applied to quite large mol. systems. Moreover, numerical tests for chem. reaction energies and intermol. interaction energies have shown that the RPAX2 method, if based on a Perdew-Burke-Ernzerhof exchange Kohn-Sham ref. determinant, yields results which are very close to coupled-cluster with single, double, and perturbative triple excitations ref. results.
- 50Eshuis, H.; Bates, J. E.; Furche, F. Electron correlation methods based on the random phase approximation. Theor. Chem. Acc. 2012, 131, 1084, DOI: 10.1007/s00214-011-1084-8There is no corresponding record for this reference.
- 51Mussard, B.; Rocca, D.; Jansen, G.; Ángyán, J. G. Dielectric Matrix Formulation of Correlation Energies in the Random Phase Approximation: Inclusion of Exchange Effects. J. Chem. Theory Comput. 2016, 12, 2191– 2202, DOI: 10.1021/acs.jctc.5b0112951https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XksVGksrk%253D&md5=6e7bf7138c73d863bdaa4353c728fc81Dielectric Matrix Formulation of Correlation Energies in the Random Phase Approximation: Inclusion of Exchange EffectsMussard, Bastien; Rocca, Dario; Jansen, Georg; Angyan, Janos G.Journal of Chemical Theory and Computation (2016), 12 (5), 2191-2202CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Starting from the general expression for the ground state correlation energy in the adiabatic-connection fluctuation-dissipation theorem (ACFDT) framework, it is shown that the dielec. matrix formulation, which is usually applied to calc. the direct RPA (dRPA) correlation energy, can be used for alternative RPA expressions including exchange effects. Within this famework, the ACFDT analog of the second order screened exchange (SOSEX) approxn. leads to a logarithmic formula for the correlation energy similar to the direct RPA expression. Alternatively, the contribution of the exchange can be included in the kernel used to evaluate the response functions. In this case, the use of an approx. kernel is crucial to simplify the formalism and to obtain a correlation energy in logarithmic form. Tech. details of the implementation of these methods are discussed, and it is shown that one can take advantage of d. fitting or Cholesky decompn. techniques to improve the computational efficiency; a discussion on the numerical quadrature made on the frequency variable is also provided. A series of test calcns. on at. correlation energies and mol. reaction energies shows that exchange effects are instrumental for improvement over direct RPA results.
- 52Dixit, A.; Ángyán, J. G.; Rocca, D. Improving the accuracy of ground-state correlation energies within a plane-wave basis set: The electron-hole exchange kernel. J. Chem. Phys. 2016, 145, 104105, DOI: 10.1063/1.496235252https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFWnurnF&md5=d7cf5cda8cc53b00dd83016f7089b527Improving the accuracy of ground-state correlation energies within a plane-wave basis set: The electron-hole exchange kernelDixit, Anant; Angyan, Janos G.; Rocca, DarioJournal of Chemical Physics (2016), 145 (10), 104105/1-104105/12CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A new formalism was recently proposed to improve RPA correlation energies by including approx. exchange effects. Within this framework, by keeping only the electron-hole contributions to the exchange kernel, two approxns. can be obtained: An adiabatic connection analog of the second order screened exchange (AC-SOSEX) and an approx. electron-hole time-dependent Hartree-Fock (eh-TDHF). Here we show how this formalism is suitable for an efficient implementation within the plane-wave basis set. The response functions involved in the AC-SOSEX and eh-TDHF equations can indeed be compactly represented by an auxiliary basis set obtained from the diagonalization of an approx. dielec. matrix. Addnl., the explicit calcn. of unoccupied states can be avoided by using d. functional perturbation theory techniques and the matrix elements of dynamical response functions can be efficiently computed by applying the Lanczos algorithm. As shown by several applications to reaction energies and weakly bound dimers, the inclusion of the electron-hole kernel significantly improves the accuracy of ground-state correlation energies with respect to RPA and semi-local functionals. (c) 2016 American Institute of Physics.
- 53Chen, G. P.; Agee, M. M.; Furche, F. Performance and Scope of Perturbative Corrections to Random-Phase Approximation Energies. J. Chem. Theory Comput. 2018, 14, 5701– 5714, DOI: 10.1021/acs.jctc.8b0077753https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhslOlsrfM&md5=15027735d3a819dacfdeceef8ba62b99Performance and Scope of Perturbative Corrections to Random-Phase Approximation EnergiesChen, Guo P.; Agee, Matthew M.; Furche, FilippJournal of Chemical Theory and Computation (2018), 14 (11), 5701-5714CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)It has been suspected since the early days of the RPA that corrections to RPA correlation energies result mostly from short-range correlation effects and are thus amenable to perturbation theory. Here we test this hypothesis by analyzing formal and numerical results for the most common beyond-RPA perturbative corrections, including the bare second-order exchange (SOX), second-order screened exchange (SOSEX), and approx. exchange kernel (AXK) methods. Our anal. is facilitated by efficient and robust algorithms based on the resoln.-of-the-identity (RI) approxn. and numerical frequency integration, which enable benchmark beyond-RPA calcns. on medium- and large-size mols. with size-independent accuracy. The AXK method systematically improves upon RPA, SOX, and SOSEX for reaction barrier heights, reaction energies, and noncovalent interaction energies of main-group compds. The improved accuracy of AXK compared with SOX and SOSEX is attributed to stronger screening of bare SOX in AXK. For reactions involving transition-metal compds., particularly 3d transition-metal dimers, the AXK correction is too small and can even have the wrong sign. These observations are rationalized by a measure ‾α of the effective coupling strength for beyond-RPA correlation. When the effective coupling strength increases beyond a crit. ‾α value of approx. 0.5, the RPA errors increase rapidly and perturbative corrections become unreliable. Thus, perturbation theory can systematically correct RPA but only for systems and properties qual. well captured by RPA, as indicated by small ‾α values.
- 54Kussmann, J.; Ochsenfeld, C. Pre-selective screening for matrix elements in linear-scaling exact exchange calculations. J. Chem. Phys. 2013, 138, 134114, DOI: 10.1063/1.479644154https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXlt1Glur0%253D&md5=ed3dfce7bf9beed628dae80664085777Pre-selective screening for matrix elements in linear-scaling exact exchange calculationsKussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Physics (2013), 138 (13), 134114/1-134114/7CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a simple but accurate preselection method based on Schwarz integral ests. to det. the significant elements of the exact exchange matrix before its evaluation, thus providing an asymptotical linear-scaling behavior for non-metallic systems. Our screening procedure proves to be highly suitable for exchange matrix calcns. on massively parallel computing architectures, such as graphical processing units, for which we present a first linear-scaling exchange matrix evaluation algorithm. (c) 2013 American Institute of Physics.
- 55Kussmann, J.; Ochsenfeld, C. Preselective Screening for Linear-Scaling Exact Exchange-Gradient Calculations for Graphics Processing Units and General Strong-Scaling Massively Parallel Calculations. J. Chem. Theory Comput. 2015, 11, 918, DOI: 10.1021/ct501189u55https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXjsFCnuro%253D&md5=b323b53aca62835e043f8ed31503c238Preselective Screening for Linear-Scaling Exact Exchange-Gradient Calculations for Graphics Processing Units and General Strong-Scaling Massively Parallel CalculationsKussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2015), 11 (3), 918-922CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present an extension of our recently presented PreLinK scheme (J. Chem. Phys. 2013, 138, 134114) for the exact exchange contribution to nuclear forces. The significant contributions to the exchange gradient are detd. by preselection based on accurate shell-pair contributions to the SCF exchange energy prior to the calcn. Therefore, our method is highly suitable for massively parallel electronic structure calcns. because of an efficient load balancing of the significant contributions only and an unhampered control flow. The efficiency of our method is shown for several illustrative calcns. on single GPU servers, as well as for hybrid MPI/CUDA parallel calcns. with the largest system comprising 3369 atoms and 26952 basis functions.
- 56Kussmann, J.; Laqua, H.; Ochsenfeld, C. Highly Efficient Resolution-of-Identity Density Functional Theory Calculations on Central and Graphics Processing Units. J. Chem. Theory Comput. 2021, 17, 1512, DOI: 10.1021/acs.jctc.0c0125256https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXksFWkurs%253D&md5=36c4001fa1dbe4ef9fdaeefa114cd851Highly Efficient Resolution-of-Identity Density Functional Theory Calculations on Central and Graphics Processing UnitsKussmann, Joerg; Laqua, Henryk; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2021), 17 (3), 1512-1521CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present an efficient method to evaluate Coulomb potential matrixes using the resoln. of identity approxn. and semilocal exchange-correlation potentials on central (CPU) and graphics processing units (GPU). The new GPU-based RI-algorithm shows a high performance and ensures the favorable scaling with increasing basis set size as the conventional CPU-based method. Furthermore, our method is based on the J-engine algorithm [White, Head-Gordon, J. Chem. Phys., 1996, 7, 2620], which allows for further optimizations that also provide a significant improvement of the corresponding CPU-based algorithm. Due to the increased performance for the Coulomb evaluation, the calcn. of the exchange-correlation potential of d. functional theory on CPUs quickly becomes a bottleneck to the overall computational time. Hence, we also present a GPU-based algorithm to evaluate the exchange-correlation terms, which results in an overall high-performance method for d. functional calcns. The algorithms to evaluate the potential and nuclear deriv. terms are discussed, and their performance on CPUs and GPUs is demonstrated for illustrative calcns.
- 57Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057, DOI: 10.1039/b515623h57https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xhs12ntrc%253D&md5=314690393f1e21096541a317a80e563cAccurate Coulomb-fitting basis sets for H to RnWeigend, FlorianPhysical Chemistry Chemical Physics (2006), 8 (9), 1057-1065CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)A series of auxiliary basis sets to fit Coulomb potentials for the elements H to Rn (except lanthanides) is presented. For each element only one auxiliary basis set is needed to approx. Coulomb energies in conjunction with orbital basis sets of split valence, triple zeta valence and quadruple zeta valence quality with errors of typically below ca. 0.15 kJ mol-1 per atom; this was demonstrated in conjunction with the recently developed orbital basis sets of types def2-SV(P), def2-TZVP and def2-QZVPP for a large set of small mols. representing (nearly) each element in all of its common oxidn. states. These auxiliary bases are slightly more than three times larger than orbital bases of split valence quality. Compared to non-approximated treatments, computation times for the Coulomb part are reduced by a factor of ca. 8 for def2-SV(P) orbital bases, ca. 25 for def2-TZVP and ca. 100 for def2-QZVPP orbital bases.
- 58Laqua, H.; Kussmann, J.; Ochsenfeld, C. Efficient and Linear-Scaling Seminumerical Method for Local Hybrid Density Functionals. J. Chem. Theory Comput. 2018, 14, 3451– 3458, DOI: 10.1021/acs.jctc.8b0006258https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhtFarurvK&md5=a22f0d7c57ca427b744970b0a450034fEfficient and Linear-Scaling Seminumerical Method for Local Hybrid Density FunctionalsLaqua, Henryk; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2018), 14 (7), 3451-3458CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)Local hybrid functionals, i.e., functionals with local dependence on the exact exchange energy d., generalize the popular class of global hybrid functionals and extend the applicability of d. functional theory to electronic structures that require an accurate description of static correlation. However, the higher computational cost compared to conventional Kohn-Sham d. functional theory restrained their widespread application. Here, we present a low-prefactor, linear-scaling method to evaluate the local hybrid exchange-correlation potential as well as the corresponding nuclear forces by employing a seminumerical integration scheme. In the seminumerical scheme, one integration in the electron repulsion integrals is carried out anal. and the other one is carried out numerically, employing an integration grid. A high computational efficiency is achieved by combining the preLinK method with explicit screening of integrals for batches of grid points to minimize the screening overhead. This new method, denoted as preLinX, provides an 8-fold performance increase for a DNA fragment contg. four base pairs as compared to existing S- and P-junction-based methods. In this way, our method allows for the evaluation of local hybrid functionals at a cost similar to that of global hybrid functionals. The linear-scaling behavior, efficiency, accuracy, and multi-node parallelization of the presented method is demonstrated for large systems contg. more than 1000 atoms.
- 59Laqua, H.; Thompson, T. H.; Kussmann, J.; Ochsenfeld, C. Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing Units. J. Chem. Theory Comput. 2020, 16, 1456, DOI: 10.1021/acs.jctc.9b0086059https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXivVejs7Y%253D&md5=606be25e20a70cc545ff71ebabf28152Highly Efficient, Linear-Scaling Seminumerical Exact-Exchange Method for Graphic Processing UnitsLaqua, Henryk; Thompson, Travis H.; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Theory and Computation (2020), 16 (3), 1456-1468CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We present a highly efficient and asymptotically linear-scaling graphic processing unit accelerated seminumerical exact-exchange method (sn-LinK). We go beyond our previous central processing unit-based method (H. Laqua et al., 2018) by employing our recently developed integral bounds (T.H. Thomson and C. Ochsenfeld, 2019) and high-accuracy numerical integration grid (H. Laqua et al., 2018). The accuracy is assessed for several established test sets, providing errors significantly below 1mEh for the smallest grid. Moreover, a comprehensive performance anal. for large mols. between 62 and 1347 atoms is provided, revealing the outstanding performance of our method, in particular, for large basis sets such as the polarized quadruple-zeta level with diffuse functions.
- 60Laqua, H.; Kussmann, J.; Ochsenfeld, C. An improved molecular partitioning scheme for numerical quadratures in density functional theory. J. Chem. Phys. 2018, 149, 204111, DOI: 10.1063/1.504943560https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXitlylurvI&md5=9cc4db762351543f44c05b2c28868ccbAn improved molecular partitioning scheme for numerical quadratures in density functional theoryLaqua, Henryk; Kussmann, Joerg; Ochsenfeld, ChristianJournal of Chemical Physics (2018), 149 (20), 204111/1-204111/8CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present a modification to Becke's mol. partitioning scheme [A. D. Becke, J. Chem. Phys. 88, 2547 (1988)] that provides substantially better accuracy for weakly bound complexes and allows for a faster and linear scaling grid generation without introducing a cutoff error. We present the accuracy of our new partitioning scheme for atomization energies of small mols. and for interaction energies of van der Waals complexes. Furthermore, the efficiency and scaling behavior of the grid generation are demonstrated for large mol. systems with up to 1707 atoms. (c) 2018 American Institute of Physics.
- 61Lehtola, S.; Steigemann, C.; Oliveira, M. J. T.; Marques, M. A. L. Recent developments in libxc–A comprehensive library of functionals for density functional theory. Software X 2018, 7, 1– 5, DOI: 10.1016/j.softx.2017.11.002There is no corresponding record for this reference.
- 62Weigend, F.; Furche, F.; Ahlrichs, R. Gaussian basis sets of quadruple zeta valence quality for atoms H–Kr. J. Chem. Phys. 2003, 119, 12753– 12762, DOI: 10.1063/1.162729362https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXpvFOns7s%253D&md5=9d1664e0004c3a900b7efcb5da10dd6dGaussian basis sets of quadruple zeta valence quality for atoms H-KrWeigend, Florian; Furche, Filipp; Ahlrichs, ReinhartJournal of Chemical Physics (2003), 119 (24), 12753-12762CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)We present Gaussian basis sets of quadruple zeta valence quality with a segmented contraction scheme for atoms H to Kr. This extends earlier work on segmented contracted split valence (SV) and triple zeta valence (TZV) basis sets. Contraction coeffs. and orbital exponents are fully optimized in at. Hartree-Fock (HF) calcns. As opposed to other quadruple zeta basis sets, the basis set errors in at. ground-state HF energies are less than 1 mEh and increase smoothly across the Periodic Table, while the no. of primitives is comparably small. Polarization functions are taken partly from previous work, partly optimized in at. MP2 calcns., and for a few cases detd. at the HF level for excited at. states nearly degenerate with the ground state. This leads to basis sets denoted QZVP for HF and d. functional theory (DFT) calcns., and for some atoms to a larger basis recommended for correlated treatments, QZVPP. We assess the performance of the basis sets in mol. HF, DFT, and MP2 calcns. for a sample of diat. and small polyat. mols. by a comparison of energies, bond lengths, and dipole moments with results obtained numerically or using very large basis sets. It is shown that basis sets of quadruple zeta quality are necessary to achieve an accuracy of 1 kcal/mol per bond in HF and DFT atomization energies. For compds. contg. third row as well as alk. and earth alk. metals it is demonstrated that the inclusion of high-lying core orbitals in the active space can be necessary for accurate correlated treatments. The QZVPP basis sets provide sufficient flexibility to polarize the core in those cases. All test calcns. indicate that the new basis sets lead to consistent accuracies in HF, DFT, or correlated treatments even in crit. cases where other basis sets may show deficiencies.
- 63Rappoport, D.; Furche, F. Property-optimized Gaussian basis sets for molecular response calculations. J. Chem. Phys. 2010, 133, 134105, DOI: 10.1063/1.348428363https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXht1Chs7nI&md5=17ee007503db923dedbb77f5e2214694Property-optimized Gaussian basis sets for molecular response calculationsRappoport, Dmitrij; Furche, FilippJournal of Chemical Physics (2010), 133 (13), 134105/1-134105/11CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)With recent advances in electronic structure methods, first-principles calcns. of electronic response properties, such as linear and nonlinear polarizabilities, have become possible for mols. with more than 100 atoms. Basis set incompleteness is typically the main source of error in such calcns. since traditional diffuse augmented basis sets are too costly to use or suffer from near linear dependence. To address this problem, we construct the first comprehensive set of property-optimized augmented basis sets for elements H-Rn except lanthanides. The new basis sets build on the Karlsruhe segmented contracted basis sets of split-valence to quadruple-zeta valence quality and add a small no. of moderately diffuse basis functions. The exponents are detd. variationally by maximization of at. Hartree-Fock polarizabilities using anal. deriv. methods. The performance of the resulting basis sets is assessed using a set of 313 mol. static Hartree-Fock polarizabilities. The mean abs. basis set errors are 3.6%, 1.1%, and 0.3% for property-optimized basis sets of split-valence, triple-zeta, and quadruple-zeta valence quality, resp. D. functional and second-order Moller-Plesset polarizabilities show similar basis set convergence. We demonstrate the efficiency of our basis sets by computing static polarizabilities of icosahedral fullerenes up to C720 using hybrid d. functional theory. (c) 2010 American Institute of Physics.
- 64Weigend, 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, DOI: 10.1039/b508541a64https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXpsFWgu7o%253D&md5=a820fb6055c993b50c405ba0fc62b194Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracyWeigend, Florian; Ahlrichs, ReinhartPhysical Chemistry Chemical Physics (2005), 7 (18), 3297-3305CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)Gaussian basis sets of quadruple zeta valence quality for Rb-Rn are presented, as well as bases of split valence and triple zeta valence quality for H-Rn. The latter were obtained by (partly) modifying bases developed previously. A large set of more than 300 mols. representing (nearly) all elements-except lanthanides-in their common oxidn. states was used to assess the quality of the bases all across the periodic table. Quantities investigated were atomization energies, dipole moments and structure parameters for Hartree-Fock, d. functional theory and correlated methods, for which we had chosen Moller-Plesset perturbation theory as an example. Finally recommendations are given which type of basis set is used best for a certain level of theory and a desired quality of results.
- 65Hättig, C. Optimization of auxiliary basis sets for RI-MP2 and RI-CC2 calculations: Core–valence and quintuple-ζ basis sets for H to Ar and QZVPP basis sets for Li to Kr. Phys. Chem. Chem. Phys. 2005, 7, 59– 66, DOI: 10.1039/B415208EThere is no corresponding record for this reference.
- 66Hellweg, A.; Hättig, C.; Höfener, S.; Klopper, W. Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to Rn. Theor. Chem. Acc. 2007, 117, 587– 597, DOI: 10.1007/s00214-007-0250-566https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXivVGmtrY%253D&md5=c140888c49b17c3671a51a45497e6284Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to RnHellweg, Arnim; Haettig, Christof; Hoefener, Sebastian; Klopper, WimTheoretical Chemistry Accounts (2007), 117 (4), 587-597CODEN: TCACFW; ISSN:1432-881X. (Springer GmbH)The introduction of the resoln.-of-the-identity (RI) approxn. for electron repulsion integrals in quantum chem. calcns. requires in addn. to the orbital basis so-called auxiliary or fitting basis sets. We report here such auxiliary basis sets optimized for second-order Moller-Plesset perturbation theory for the recently published (Weigend and Ahlrichs Phys. Chem. Chem. Phys., 2005, 7, 3297-3305) segmented contracted Gaussian basis sets of split, triple-ζ and quadruple-ζ valence quality for the atoms Rb-Rn (except lanthanides). These basis sets are designed for use in connection with small-core effective core potentials including scalar relativistic corrections. Hereby accurate resoln.-of-the-identity calcns. with second-order Moller-Plesset perturbation theory (MP2) and related methods can now be performed for mols. contg. elements from H to Rn. The error of the RI approxn. has been evaluated for a test set of 385 small and medium sized mols., which represent the common oxidn. states of each element, and is compared with the one-electron basis set error, estd. based on highly accurate explicitly correlated MP2-R12 calcns. With the reported auxiliary basis sets the RI error for MP2 correlation energies is typically two orders of magnitude smaller than the one-electron basis set error, independent on the position of the atoms in the periodic table.
- 67Hellweg, A.; Rappoport, D. Development of new auxiliary basis functions of the Karlsruhe segmented contracted basis sets including diffuse basis functions (def2-SVPD, def2-TZVPPD, and def2-QZVPPD) for RI-MP2 and RI-CC calculations. Phys. Chem. Chem. Phys. 2015, 17, 1010– 1017, DOI: 10.1039/C4CP04286GThere is no corresponding record for this reference.
- 68Kaltak, M.; Klimeš, J.; Kresse, G. Low Scaling Algorithms for the Random Phase Approximation: Imaginary Time and Laplace Transformations. J. Chem. Theory Comput. 2014, 10, 2498, DOI: 10.1021/ct500126868https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXmslGlu70%253D&md5=032579aea73969146da834d532c29997Low Scaling Algorithms for the Random Phase Approximation: Imaginary Time and Laplace TransformationsKaltak, Merzuk; Klimes, Jiri; Kresse, GeorgJournal of Chemical Theory and Computation (2014), 10 (6), 2498-2507CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)We det. efficient imaginary frequency and imaginary time grids for second-order Moller-Plesset (MP) perturbation theory. The least-squares and Minimax quadratures are compared for periodic systems, finding that the Minimax quadrature performs slightly better for the considered materials. We show that the imaginary frequency grids developed for second order also perform well for the correlation energy in the direct RPA. Furthermore, we show that the polarizabilities on the imaginary time axis can be Fourier-transformed to the imaginary frequency domain, since the time and frequency Minimax grids are dual to each other. The same duality is obsd. for the least-squares grids. The transformation from imaginary time to imaginary frequency allows one to reduce the time complexity to cubic (in system size), so that RPA correlation energies become accessible for large systems.
- 69Kaltak, M.; Kresse, G. Minimax isometry method: A compressive sensing approach for Matsubara summation in many-body perturbation theory. Phys. Rev. B 2020, 101, 205145, DOI: 10.1103/PhysRevB.101.20514569https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhsFChtLjJ&md5=0e8b7dcabc1bfd1ddbac581ec0b91706Minimax isometry method: A compressive sensing approach for Matsubara summation in many-body perturbation theoryKaltak, Merzuk; Kresse, GeorgPhysical Review B (2020), 101 (20), 205145CODEN: PRBHB7; ISSN:2469-9969. (American Physical Society)We present a compressive sensing approach for the long-standing problem of Matsubara summation in many-body perturbation theory. By constructing low-dimensional, almost isometric subspaces of the Hilbert space we obtain optimum imaginary time and frequency grids that allow for extreme data compression of fermionic and bosonic functions in a broad temp. regime. The method is applied to the random phase and self-consistent GW approxn. of the grand potential. Integration and transformation errors are investigated for Si and SrVO3.
- 70Andrae, D.; Häußermann, U.; Dolg, M.; Stoll, H.; Preuß, H. Energy-adjusted ab initio pseudopotentials for the second and third row transition elements. Theor. Chim. Acta 1990, 77, 123– 141, DOI: 10.1007/BF0111453770https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3cXkt12ntLo%253D&md5=8203c6bc6149924cbd4b23b3063715e1Energy-adjusted ab initio pseudopotentials for the second and third row transition elementsAndrae, D.; Haeussermann, U.; Dolg, M.; Stoll, H.; Preuss, H.Theoretica Chimica Acta (1990), 77 (2), 123-41CODEN: TCHAAM; ISSN:0040-5744.Nonrelativistic and quasirelativistic ab initio pseudopotentials substituting the M(z-28)+-core orbitals of the second row transition elements and M(z-60)+-core orbitals of the third row transition elements, resp., and optimized (8s7p6d)/[6s5p3d]-GTO valence basis sets for use in mol. calcns. have been generated. Addnl., corresponding spin-orbit operators have also been derived. At. excitation and ionization energies from numerical HF as well as from SCF pseudopotential calcns. using the derived basis sets differ in most cases by less than 0.1 eV from corresponding numerical all-electron results. Spin-orbit splittings for low-lying states are in reasonable agreement with corresponding all-electron Dirac-Fock (DF) results.
- 71Metz, B.; Stoll, H.; Dolg, M. Small-core multiconfiguration-Dirac–Hartree–Fock-adjusted pseudopotentials for post-d main group elements: Application to PbH and PbO. J. Chem. Phys. 2000, 113, 2563– 2569, DOI: 10.1063/1.130588071https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3cXlsVWmsrg%253D&md5=3b605722931765b461713ea2fe318c7eSmall-core multiconfiguration-Dirac-Hartree-Fock-adjusted pseudopotentials for post-d main group elements: Application to PbH and PbOMetz, Bernhard; Stoll, Hermann; Dolg, MichaelJournal of Chemical Physics (2000), 113 (7), 2563-2569CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Relativistic pseudopotentials (PPs) of the energy-consistent variety have been generated for the post-d group 13-15 elements, by adjustment to multiconfiguration Dirac-Hartree-Fock data based on the Dirac-Coulomb-Breit Hamiltonian. The outer-core (n-1)spd shells are explicitly treated together with the nsp valence shell, with these PPs, and the implications of the small-core choice are discussed by comparison to a corresponding large-core PP, in the case of Pb. Results from valence ab initio one- and two-component calcns. using both PPs are presented for the fine-structure splitting of the ns2np2 ground-state configuration of the Pb atom, and for spectroscopic consts. of PbH (X 2Π1/2, 2Π3/2) and PbO (X 1Σ+). In addn., a combination of small-core and large-core PPs has been explored in spin-free-state shifted calcns. for the above mols.
- 72Peterson, K. A. Systematically convergent basis sets with relativistic pseudopotentials. I. Correlation consistent basis sets for the post-d group 13–15 elements. J. Chem. Phys. 2003, 119, 11099– 11112, DOI: 10.1063/1.162292372https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXovFyms74%253D&md5=ab472eed340dd9b9ac5e2d7d2e76df7aSystematically convergent basis sets with relativistic pseudopotentials. I. Correlation consistent basis sets for the post-d group 13-15 elementsPeterson, Kirk A.Journal of Chemical Physics (2003), 119 (21), 11099-11112CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)New correlation consistent-like basis sets have been developed for the post-d group 13-15 elements (Ga-As, In-Sb, Tl-Bi) employing accurate, small-core relativistic pseudopotentials. The resulting basis sets, which are denoted cc-pVnZ-PP, are appropriate for valence electron correlation and range in size from (8s7p7d)/[4s3p2d] for the cc-pVDZ-PP to (16s13p12d3f2g1h)/[7s7p5d3f2g1h] for the cc-pV5Z-PP sets. Benchmark calcns. on selected diat. mols. (As2, Sb2, Bi2, AsN, SbN, BiN, GeO, SnO, PbO, GaCl, InCl, TlCl, GaH, InH, and TlH) are reported using these new basis sets at the coupled cluster level of theory. Much like their all-electron counterparts, the cc-pVnZ-PP basis sets yield systematic convergence of total energies and spectroscopic consts. In several cases all-electron benchmark calcns. were also carried out for comparison. The results from the pseudopotential and all-electron calcns. were nearly identical when scalar relativity was accurately included in the all-electron work. Diffuse-augmented basis sets, aug-cc-pVnZ-PP, have also been developed and have been used in calcns. of the at. electron affinities.
- 73Peterson, K. A.; Figgen, D.; Goll, E.; Stoll, H.; Dolg, M. Systematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis sets for the post-d group 16–18 elements. J. Chem. Phys. 2003, 119, 11113– 11123, DOI: 10.1063/1.162292473https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3sXovFyms7w%253D&md5=2f5869fe14c68b766761d8d160acce6aSystematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis sets for the post-d group 16-18 elementsPeterson, Kirk A.; Figgen, Detlev; Goll, Erich; Stoll, Hermann; Dolg, MichaelJournal of Chemical Physics (2003), 119 (21), 11113-11123CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A series of correlation consistent basis sets have been developed for the post-d group 16-18 elements in conjunction with small-core relativistic pseudopotentials of the energy-consistent variety. The latter were adjusted to multiconfiguration Dirac-Hartree-Fock data based on the Dirac-Coulomb-Breit Hamiltonian. The outer-core (n-1)spd shells are explicitly treated together with the nsp valence shell with these PPs. The accompanying cc-pVnZ-PP and aug-cc-pVnZ-PP basis sets range in size from DZ to 5Z quality and yield systematic convergence of both Hartree-Fock and correlated total energies. In addn. to the calcn. of at. electron affinities and dipole polarizabilities of the rare gas atoms, numerous mol. benchmark calcns. (HBr, HI, HAt, Br2, I2, At2, SiSe, SiTe, SiPo, KrH+, XeH+, and RnH+) are also reported at the coupled cluster level of theory. For the purposes of comparison, all-electron calcns. using the Douglas-Kroll-Hess Hamiltonian have also been carried out for the halogen-contg. mols. using basis sets of 5Z quality.
- 74Morgante, P.; Peverati, R. Statistically representative databases for density functional theory via data science. Phys. Chem. Chem. Phys. 2019, 21, 19092– 19103, DOI: 10.1039/C9CP03211H74https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhsFygsL%252FO&md5=ed89b2e280632ec164cf0101e07f3fd5Statistically representative databases for density functional theory via data scienceMorgante, Pierpaolo; Peverati, RobertoPhysical Chemistry Chemical Physics (2019), 21 (35), 19092-19103CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)The amt. of data and no. of databases for the assessment and parameterization of d. functional theory methods has grown substantially in the past two decades. In this work, we introduce a novel cluster anal. technique for d. functional theory calcns. of the electronic structure of atoms and mols. with the goal of creating new statistically significant databases with broad chem. scope, and a manageable no. of data-points. By analyzing without a priori chem. assumptions a population of almost 350k data-points, we create a new database called ASCDB contg. only 200 data-points. This new database holds the same chem. information as the larger population of data from which it is obtained, but with a computational cost that is reduced by several orders of magnitude. The labeling of the significant chem. properties is performed a posteriori on the resulting 16 subsets, classifying them into four areas of chem. importance: non-covalent interactions, thermochem., non-local effects, and unbiased calcns. The anal. of the results and their transferability shows that ASCDB is capable of providing the same information as that of the larger collection of data-such as GMTKN55, MGCDB84, and Minnesota 2015B-for several d. functional theory methods and basis sets. In light of these results, we suggest the use of this new small database as a first inexpensive tool for the evaluation and parameterization of electronic structure theory methods.
- 75Morgante, P.; Peverati, R. ACCDB: A collection of chemistry databases for broad computational purposes. J. Comput. Chem. 2019, 40, 839– 848, DOI: 10.1002/jcc.2576175https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXisF2htLzI&md5=2e99af73352558db4926d2065f47dcdcACCDB: A collection of chemistry databases for broad computational purposesMorgante, Pierpaolo; Peverati, RobertoJournal of Computational Chemistry (2019), 40 (6), 839-848CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)The importance of databases of reliable and accurate data in chem. has substantially increased in the past two decades. Their main usage is to parametrize electronic structure theory methods, and to assess their capabilities and accuracy for a broad set of chem. problems. The collection we present here-ACCDB-includes data from 16 different research groups, for a total of 44,931 unique ref. data points, all at a level of theory significantly higher than d. functional theory, and covering most of the periodic table. It is composed of five databases taken from literature (GMTKN, MGCDB84, Minnesota2015, DP284, and W4-17), two newly developed reaction energy databases (W4-17-RE and MN-RE), and a new collection of databases contg. transition metals. A set of expandable software tools for its manipulation is also presented here for the first time, as well as a case study where ACCDB is used for benchmarking com. CPUs for chem. calcns. © 2018 Wiley Periodicals, Inc.
- 76Fritsch, F. N.; Butland, J. A Method for Constructing Local Monotone Piecewise Cubic Interpolants. SIAM J. Sci. Stat. Comp. 1984, 5, 300, DOI: 10.1137/0905021There is no corresponding record for this reference.
- 77Peverati, R. Fitting elephants in the density functionals zoo: Statistical criteria for the evaluation of density functional theory methods as a suitable replacement for counting parameters. Int. J. Quantum Chem. 2021, 121, e26379, DOI: 10.1002/qua.2637977https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhtlersrzK&md5=278ba83220ff86d50f5309e7f76151ecFitting elephants in the density functionals zoo: Statistical criteria for the evaluation of density functional theory methods as a suitable replacement for counting parametersPeverati, RobertoInternational Journal of Quantum Chemistry (2021), 121 (1), e26379CODEN: IJQCB2; ISSN:0020-7608. (John Wiley & Sons, Inc.)Counting parameters has become customary in the d. functional theory community as a way to infer the transferability of popular approxns. to the exchange-correlation functionals. Recent work in data science, however, has demonstrated that the no. of parameters of a fitted model is not related to the complexity of the model itself, nor to its eventual overfitting. Using similar arguments, here, we show that it is possible to represent every modern exchange-correlation functional approxns. using just one single parameter. This procedure proves the futility of the no. of parameters as a measure of transferability. To counteract this shortcoming, we introduce and analyze the performance of three statistical criteria for the evaluation of the transferability of exchange-correlation functionals. The three criteria are called Akaike information criterion, Vapnik-Chervonenkis criterion, and cross-validation criterion and are used in a preliminary assessment to rank 60 exchange-correlation functional approxns. using the ASCDB database of chem. data.
- 78Korth, M.; Grimme, S. Mindless DFT Benchmarking. J. Chem. Theory Comput. 2009, 5, 993– 1003, DOI: 10.1021/ct800511q78https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXis12itbY%253D&md5=d3a2772137b30cf4f0343825104798a5"Mindless" DFT BenchmarkingKorth, Martin; Grimme, StefanJournal of Chemical Theory and Computation (2009), 5 (4), 993-1003CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A diversity-oriented approach for the generation of thermochem. benchmark sets is presented. Test sets consisting of randomly generated "artificial mols." (AMs) are proposed that rely on systematic constraints rather than uncontrolled chem. biases. In this way, the narrow structural space of chem. intuition is opened up and electronically difficult cases can be produced in an unforeseeable manner. For the calcn. of chem. meaningful relative energies, AMs are systematically decompd. into small mols. (hydrides and diatomics). Two different example test sets contg. eight-atom, single-ref., main group AMs with chem. very diverse and unusual structures are generated. Highly accurate all-electron, estd. CCSD(T)/complete basis set ref. energies are also provided. They are used to benchmark the d. functionals S-VWN, BP86, B-LYP, B97-D, PBE, TPSS, PBEh, BH-LYP, B3-PW91, B3-LYP, B2-PLYP, B2GP-PLYP, BMK, MPW1B95, M05, M05-2X, PW6B95, M06, M06-L, and M06-2X. In selected cases, an empirical dispersion correction (DFT-D) has been applied. Due to the compn. of the sets, it is expected that a good performance indicates "robustness" in many different chem. applications. The results of a statistical anal. of the errors for the entire set with 165 entries (av. reaction energy of 117 kcal/mol, dubbed as the MB08-165 set) perfectly fit to the "Jacob's ladder" metaphor for the ordering of d. functionals according to their theor. complexity. The mean abs. deviation (MAD) decreases very strongly from LDA (20 kcal/mol) to GGAs (MAD of about 10 kcal/mol) but then was less pronounced to hybrid-GGAs (MAD of about 6-8 kcal/mol). The best performance (MAD of 4.1-4.2 kcal/mol) is found for the (fifth-rung) double-hybrid functionals B2-PLYP-D and B2GP-PLYP-D, followed by the M06-2X meta-hybrid (MAD of 4.8 kcal/mol). The significance of the proposed approach for thermodn. benchmarking is discussed and related to the obsd. performance ranking also regarding wave function based methods.
- 79Broyden, C. G. The Convergence of a Class of Double-rank Minimization Algorithms 1. General Considerations. IMA J. Appl. Math. 1970, 6, 76, DOI: 10.1093/imamat/6.1.76There is no corresponding record for this reference.
- 80Fletcher, R. A new approach to variable metric algorithms. Comput. J. 1970, 13, 317, DOI: 10.1093/comjnl/13.3.317There is no corresponding record for this reference.
- 81Goldfarb, D. A family of variable-metric methods derived by variational methods. Math. Comput. 1970, 24, 23, DOI: 10.1090/S0025-5718-1970-0258249-6There is no corresponding record for this reference.
- 82Shanno, D. F. Conditioning of quasi-Newton methods for function minimization. Math. Comput. 1970, 24, 647, DOI: 10.1090/S0025-5718-1970-0274029-XThere is no corresponding record for this reference.
- 83Parthiban, S.; Martin, J. M. L. Assessment of W1 and W2 theories for the computation of electron affinities, ionization potentials, heats of formation, and proton affinities. J. Chem. Phys. 2001, 114, 6014– 6029, DOI: 10.1063/1.135601483https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXitlKls7w%253D&md5=c52c0d86eab1a2f42f893ecc19bac2adAssessment of W1 and W2 theories for the computation of electron affinities, ionization potentials, heats of formation, and proton affinitiesParthiban, Srinivasan; Martin, Jan M. L.Journal of Chemical Physics (2001), 114 (14), 6014-6029CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The performance of two recent ab initio computational thermochem. schemes, W1 and W2 theory [J. M. L. Martin and G. de Oliveira, J. Chem. Phys. 111, 1843 (1999)], is assessed for an enlarged sample of thermochem. data consisting of the ionization potentials and electron affinities in the G2-1 and G2-2 sets, as well as the heats of formation in the G2-1 and a subset of the G2-2 set. We find W1 theory to be several times more accurate for ionization potentials and electron affinities than commonly used (and less expensive) computational thermochem. schemes such as G2, G3, and CBS-QB3: W2 theory represents a slight improvement for electron affinities but no significant one for ionization potentials. The use of a two-point A+B/L5 rather than a three-point A+B/CL extrapolation for the SCF component greatly enhances the numerical stability of the W1 method for systems with slow basis set convergence. Inclusion of first-order spin-orbit coupling is essential for accurate ionization potentials and electron affinities involving degenerate electronic states: Inner-shell correlation is somewhat more important for ionization potentials than for electron affinities, while scalar relativistic effects are required for the highest accuracy. The mean deviation from expt. for the G2-1 heats of formation is within the av. exptl. uncertainty. W1 theory appears to be a valuable tool for obtaining benchmark quality proton affinities.
- 84Karton, A.; Sylvetsky, N.; Martin, J. M. L. W4–17: A diverse and high-confidence dataset of atomization energies for benchmarking high-level electronic structure methods. J. Comput. Chem. 2017, 38, 2063, DOI: 10.1002/jcc.2485484https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXhtV2gsL3K&md5=4e8a303a60ccb1e846019a1f69d61d19W4-17: A diverse and high-confidence dataset of atomization energies for benchmarking high-level electronic structure methodsKarton, Amir; Sylvetsky, Nitai; Martin, Jan M. L.Journal of Computational Chemistry (2017), 38 (24), 2063-2075CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Atomization reactions are among the most challenging tests for electronic structure methods. We use the first-principles Weizmann-4 (W4) computational thermochem. protocol to generate the W4-17 dataset of 200 total atomization energies (TAEs) with 3σ confidence intervals of 1 kJ mol-1. W4-17 is an extension of the earlier W4-11 dataset; it includes first- and second-row mols. and radicals with up to eight non-hydrogen atoms. These cover a broad spectrum of bonding situations and multireference character, and as such are an excellent benchmark for the parameterization and validation of highly accurate ab initio methods (e.g., CCSD(T) composite procedures) and double-hybrid d. functional theory (DHDFT) methods. The W4-17 dataset contains two subsets (i) a non-multireference subset of 183 systems characterized by dynamical or moderate nondynamical correlation effects (denoted W4-17-nonMR) and (ii) a highly multireference subset of 17 systems (W4-17-MR). We use these databases to evaluate the performance of a wide range of CCSD(T) composite procedures (e.g., G4, G4(MP2), G4(MP2)-6X, ROG4(MP2)-6X, CBS-QB3, ROCBS-QB3, CBS-APNO, ccCA-PS3, W1, W2, W1-F12, W2-F12, W1X-1, and W2X) and DHDFT methods (e.g., B2-PLYP, B2GP-PLYP, B2K-PLYP, DSD-BLYP, DSD-PBEP86, PWPB95, ωB97X-2(LP), and ωB97X-2(TQZ)). © 2017 Wiley Periodicals, Inc.
- 85Karton, A.; Daon, S.; Martin, J. M. W4–11: A high-confidence benchmark dataset for computational thermochemistry derived from first-principles W4 data. Chem. Phys. Lett. 2011, 510, 165, DOI: 10.1016/j.cplett.2011.05.00785https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXnvF2qsLs%253D&md5=41b53b92e3d9c4f450bc55a5e5df4eb7W4-11: A high-confidence benchmark dataset for computational thermochemistry derived from first-principles W4 dataKarton, Amir; Daon, Shauli; Martin, Jan M. L.Chemical Physics Letters (2011), 510 (4-6), 165-178CODEN: CHPLBC; ISSN:0009-2614. (Elsevier B.V.)We show that the purely first-principles Weizmann-4 (W4) computational thermochem. method developed in our group can reproduce available Active Thermochem. Tables atomization energies for 35 mols. with a 3σ uncertainty of under 1 kJ/mol. We then employ this method to generate the W4-11 dataset of 140 total atomization energies of small first-and second-row mols. and radicals. These cover a broad spectrum of bonding situations and multireference character, and as such are an excellent, quasi-automated benchmark (available electronically as Supporting Information) for parametrization and validation of more approx. methods (such as DFT functionals and composite methods). Secondary contributions such as relativity can be included or omitted at will, unlike with exptl. data. A broad variety of more approx. methods is assessed against the W4-11 benchmark and recommendations are made.
- 86Margraf, J. T.; Ranasinghe, D. S.; Bartlett, R. J. Automatic generation of reaction energy databases from highly accurate atomization energy benchmark sets. Phys. Chem. Chem. Phys. 2017, 19, 9798, DOI: 10.1039/C7CP00757D86https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXksVOqur4%253D&md5=4f425ab2aa394430e939871dcf8925b0Automatic generation of reaction energy databases from highly accurate atomization energy benchmark setsMargraf, Johannes T.; Ranasinghe, Duminda S.; Bartlett, Rodney J.Physical Chemistry Chemical Physics (2017), 19 (15), 9798-9805CODEN: PPCPFQ; ISSN:1463-9076. (Royal Society of Chemistry)In this contribution, we discuss how reaction energy benchmark sets can automatically be created from arbitrary atomization energy databases. As an example, over 11 000 reaction energies derived from the W4-11 database, as well as some relevant subsets are reported. Importantly, there is only very modest computational overhead involved in computing >11 000 reaction energies compared to 140 atomization energies, since the rate-detg. step for either benchmark is performing the same 140 quantum chem. calcns. The performance of commonly used electronic structure methods for the new database is analyzed. This allows investigating the relationship between the performances for atomization and reaction energy benchmarks based on an identical set of mols. The atomization energy is found to be a weak predictor for the overall usefulness of a method. The performance of d. functional approxns. in light of the no. of empirically optimized parameters used in their design is also discussed.
- 87Henderson, T. M.; Scuseria, G. E. The connection between self-interaction and static correlation: A random phase approximation perspective. Mol. Phys. 2010, 108, 2511, DOI: 10.1080/00268976.2010.50722787https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXhtl2is7rJ&md5=62434dbb5a179a1b99b4dc07e0c25682The connection between self-interaction and static correlation: a random phase approximation perspectiveHenderson, Thomas M.; Scuseria, Gustavo E.Molecular Physics (2010), 108 (19-20), 2511-2517CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)Semi-local d. functional theory suggests a connection between static correlation and self-interaction. It is difficult to make such a connection from the wave function theory perspective, since few wave function methods permit self-interaction error. However, the RPA for ground-state correlation, which has a wave function derivation, does include self-interaction in its direct (Hartree) variant. This variant also describes left-right correlation. The self-interaction can be removed by means of second-order screened exchange; however, this also has neg. consequences for the description of static correlation. This paper discusses the connection between the two concepts (static correlation and self-interaction) from the perspective provided by the RPA.
- 88Grüning, M.; Gritsenko, O. V.; van Gisbergen, S. J. A.; Baerends, E. J. The Failure of Generalized Gradient Approximations (GGAs) and Meta-GGAs for the Two-Center Three-Electron Bonds in He2+, (H2O)2+, and (NH3)2+. J. Phys. Chem. A 2001, 105, 9211, DOI: 10.1021/jp011239k88https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXmvVWrsrc%253D&md5=e3944e590f1fc9aa215d593b809725e5The Failure of Generalized Gradient Approximations (GGAs) and Meta-GGAs for the Two-Center Three-Electron Bonds in He2+, (H2O)2+, and (NH3)2+Gruening, M.; Gritsenko, O. V.; van Gisbergen, S. J. A.; Baerends, E. J.Journal of Physical Chemistry A (2001), 105 (40), 9211-9218CODEN: JPCAFH; ISSN:1089-5639. (American Chemical Society)The radical cations He2+, (H2O)2+, and (NH3)2+ with two-center three-electron A-A bonds are investigated at the CI, accurate Kohn-Sham (KS), generalized gradient approxn. (GGA), and meta-GGA levels. Assessment of seven different GGA and six meta-GGA methods shows that the A2+ systems remain a difficult case for d. functional theory (DFT). All methods tested consistently overestimate the stability of A2+: the corresponding De errors decrease for more diffuse valence densities in the series He2+ > (H2O)2+ > (NH3)2+. Upon comparison to the energy terms of the accurate Kohn-Sham solns., the approx. exchange functionals are found to be responsible for the errors of GGA-type methods, which characteristically overestimate the exchange in A2+. These so-called exchange functionals implicitly use localized holes. Such localized holes do occur if there is left-right correlation, i.e., the exchange functionals then also describe non-dynamical correlation. However, in the hemibonded A2+ systems the typical mol. (left-right, non-dynamical) correlation of the two-electron pair bond is absent. The non-dynamical correlation built into the exchange functionals is then spurious and yields too low energies.
- 89Peterson, K. A.; Woon, D. E.; Dunning, T. H. Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H+H2→H2+H reaction. J. Chem. Phys. 1994, 100, 7410– 7415, DOI: 10.1063/1.46688489https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2cXjtlWmurY%253D&md5=c8b0f6e0bd2b066d23bc77bddaf0be36Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H + H2 → H2 + H reactionPeterson, Kirk A.; Woon, David E.; Dunning, Thom H., Jr.Journal of Chemical Physics (1994), 100 (10), 7410-15CODEN: JCPSA6; ISSN:0021-9606.Using systematic sequences of correlation consistent Gaussian basis sets from double to sextuple zeta quality, the classical barrier height of the H + H2 exchange reaction has been calcd. by multireference CI (MRCI) methods. The MRCI calcns. for collinear H3 have also been calibrated against large basis set full CI (FCI) results, which demonstrate that the MRCI treatment leads to energies less than 1 μhartree (<0.001 kcal/mol) above the FCI energies. The dependence of both the H2 and H3 total energies on the basis set is found to be very regular, and this behavior has been used to extrapolate to the complete basis set (CBS) limits. The resulting est. of the H-H-H- CBS limit yields a classical barrier height, relative to exact H + H2, of 9.60 ± 0.02 kcal/mol; the best directly calcd. value for the barrier is equal to 9.62 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo calcns.
- 90Helgaker, T.; Klopper, W.; Koch, H.; Noga, J. Basis-set convergence of correlated calculations on water. J. Chem. Phys. 1997, 106, 9639– 9646, DOI: 10.1063/1.47386390https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXjvVCgu78%253D&md5=f4689c1b38fe30eb721e9cd7d607bdf7Basis-set convergence of correlated calculations on waterHelgaker, Trygve; Klopper, Wim; Koch, Henrik; Noga, JozefJournal of Chemical Physics (1997), 106 (23), 9639-9646CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The basis-set convergence of the electronic correlation energy in the water mol. is investigated at the second-order Moller-Plesset level and at the coupled-cluster singles-and-doubles level with and without perturbative triples corrections applied. The basis-set limits of the correlation energy are established to within 2mEh by means of (1) extrapolations from sequences of calcns. using correlation-consistent basis sets and (2) from explicitly correlated calcns. employing terms linear in the inter-electronic distances rij. For the extrapolations to the basis-set limit of the correlation energies, fits of the form a + bX-3 (where X is two for double-zeta sets, three for triple-zeta sets, etc.) are found to be useful. CCSD(T) calcns. involving as many as 492 AOs are reported.
- 91Halkier, A.; Helgaker, T.; Jo̷rgensen, P.; Klopper, W.; Koch, H.; Olsen, J.; Wilson, A. K. Basis-set convergence in correlated calculations on Ne, N2, and H2O. Chem. Phys. Lett. 1998, 286, 243– 252, DOI: 10.1016/S0009-2614(98)00111-091https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK1cXitVGqsLo%253D&md5=04274821d9c7fa664e9588855ed9a061Basis-set convergence in correlated calculations on Ne, N2, and H2OHalkier, Asger; Helgaker, Trygve; Jorgensen, Poul; Klopper, Wim; Koch, Henrik; Olsen, Jeppe; Wilson, Angela K.Chemical Physics Letters (1998), 286 (3,4), 243-252CODEN: CHPLBC; ISSN:0009-2614. (Elsevier Science B.V.)Valence and all-electron correlation energies of Ne, N2, and H2O at fixed exptl. geometries are computed at the levels of second-order perturbation theory (MP2) and coupled cluster theory with singles and doubles excitations (CCSD), and singles and doubles excitations with a perturbative triples correction (CCSD(T)). Correlation-consistent polarized valence and core-valence basis sets up to sextuple zeta quality are employed. Guided by basis-set limits established by rij-dependent methods, a no. of extrapolation schemes for use with the correlation-consistent basis sets are investigated. Among the schemes considered here, a linear least-squares procedure applied to the quintuple and sextuple zeta results yields the most accurate extrapolations.
- 92Dunning, J.; Thom, H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007– 1023, DOI: 10.1063/1.45615392https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL1MXksVGmtrk%253D&md5=c6cd67a3748dc61692a9cb622d2694a0Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogenDunning, Thom H., Jr.Journal of Chemical Physics (1989), 90 (2), 1007-23CODEN: JCPSA6; ISSN:0021-9606.Guided by the calcns. on oxygen in the literature, basis sets for use in correlated at. and mol. calcns. were developed for all of the first row atoms from boron through neon, and for hydrogen. As in the oxygen atom calcns., the incremental energy lowerings, due to the addn. of correlating functions, fall into distinct groups. This leads to the concept of correlation-consistent basis sets, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation-consistent sets are given for all of the atoms considered. The most accurate sets detd. in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding at.-natural-orbital sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estd. that this set yields 94-97% of the total (HF + 1 + 2) correlation energy for the atoms neon through boron.
- 93Kendall, R. A.; Dunning, T. H.; Harrison, R. J. Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions. J. Chem. Phys. 1992, 96, 6796– 6806, DOI: 10.1063/1.46256993https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK38XktFClurw%253D&md5=948a06eee10604a8fa37eae2b2ada4beElectron affinities of the first-row atoms revisited. Systematic basis sets and wave functionsKendall, Rick A.; Dunning, Thom H., Jr.; Harrison, Robert J.Journal of Chemical Physics (1992), 96 (9), 6796-806CODEN: JCPSA6; ISSN:0021-9606.The authors describe a reliable procedure for calcg. the electron affinity of an atom and present results for H, B, C, O, and F (H is included for completeness). This procedure involves the use of the recently proposed correlation-consistent basis sets augmented with functions to describe the more diffuse character of the at. anion coupled with a straightforward, uniform expansion of the ref. space for multireference singles and doubles configuration-interaction (MRSD-CI) calcns. A comparison is given with previous results and with corresponding full CI calcns. The most accurate EAs obtained from the MRSD-CI calcns. are (with exptl. values in parentheses): H 0.740 eV (0.754), B 0.258 (0.277), C 1.245 (1.263), O 1.384 (1.461), and F 3.337 (3.401). The EAs obtained from the MR-SDCI calcns. differ by less than 0.03 eV from those predicted by the full CI calcns.
- 94Peterson, K. A.; Dunning, T. H. Accurate correlation consistent basis sets for molecular core–valence correlation effects: The second row atoms Al–Ar, and the first row atoms B–Ne revisited. J. Chem. Phys. 2002, 117, 10548– 10560, DOI: 10.1063/1.152013894https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XptFGgu70%253D&md5=146a124710aa803859c1a67e09b45b0aAccurate correlation consistent basis sets for molecular core-valence correlation effects: The second row atoms Al-Ar, and the first row atoms B-Ne revisitedPeterson, Kirk A.; Dunning, Thom H., Jr.Journal of Chemical Physics (2002), 117 (23), 10548-10560CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Correlation consistent basis sets for accurately describing core-core and core-valence correlation effects in atoms and mols. have been developed for the second row atoms Al-Ar. Two different optimization strategies were investigated, which led to two families of core-valence basis sets when the optimized functions were added to the std. correlation consistent basis sets (cc-pVnZ). In the first case, the exponents of the augmenting primitive Gaussian functions were optimized with respect to the difference between all-electron and valence-electron correlated calcns., i.e., for the core-core plus core-valence correlation energy. This yielded the cc-pCVnZ family of basis sets, which are analogous to the sets developed previously for the first row atoms [D. E. Woon and T. H. Dunning, Jr., J. Chem. Phys. 103, 4572 (1995)]. Although the cc-pCVnZ sets exhibit systematic convergence to the all-electron correlation energy at the complete basis set limit, the intershell (core-valence) correlation energy converges more slowly than the intrashell (core-core) correlation energy. Since the effect of including the core electrons on the calcn. of mol. properties tends to be dominated by core-valence correlation effects, a second scheme for detg. the augmenting functions was investigated. In this approach, the exponents of the functions to be added to the cc-pVnZ sets were optimized with respect to just the core-valence (intershell) correlation energy, except that a small amt. of core-core correlation energy was included in order to ensure systematic convergence to the complete basis set limit. These new sets, denoted weighted core-valence basis sets (cc-pwCVnZ), significantly improve the convergence of many mol. properties with n. Optimum cc-pwCVnZ sets for the first-row atoms were also developed and show similar advantages. Both the cc-pCVnZ and cc-pwCVnZ basis sets were benchmarked in coupled cluster [CCSD(T)] calcns. on a series of second row homonuclear diat. mols. (Al2, Si2, P2, S2, and Cl2), as well as on selected diat. mols. involving first row atoms (CO, SiO, PN, and BCl). For the calcn. of core correlation effects on energetic and spectroscopic properties, the cc-pwCVnZ basis sets are recommended over the cc-pCVnZ ones.
- 95Weigend, F.; Köhn, A.; Hättig, C. Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations. J. Chem. Phys. 2002, 116, 3175– 3183, DOI: 10.1063/1.144511595https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38XhtlSiu7k%253D&md5=0130fa656254a693e80d4be6b0f442b8Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculationsWeigend, Florian; Kohn, Andreas; Hattig, ChristofJournal of Chemical Physics (2002), 116 (8), 3175-3183CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The convergence of the second-order Moller-Plesset perturbation theory (MP2) correlation energy with the cardinal no. X is investigated for the correlation consistent basis-set series cc-pVXZ and cc-pV(X+d)Z. For the aug-cc-pVXZ and aug-cc-pV(X+d)Z series the convergence of the MP2 correlation contribution to the dipole moment is studied. It is found that, when d-shell electrons cannot be frozen, the cc-pVXZ and aug-cc-pVXZ basis sets converge much slower for third-row elements then they do for first- and second-row elements. Based on the results of these studies criteria are deduced for the accuracy of auxiliary basis sets used in the resoln. of the identity (RI) approxn. for electron repulsion integrals. Optimized auxiliary basis sets for RI-MP2 calcns. fulfilling these criteria are reported for the sets cc-pVXZ, cc-pV(X+d)Z, aug-cc-pVXZ, and aug-cc-pV(X+d)Z with X=D, T, and Q. For all basis sets the RI error in the MP2 correlation energy is more than two orders of magnitude smaller than the usual basis-set error. For the auxiliary aug-cc-pVXZ and aug-cc-pV(X+d)Z sets the RI error in the MP2 correlation contribution to the dipole moment is one order of magnitude smaller than the usual basis set error. Therefore extrapolations towards the basis-set limit are possible within the RI approxn. for both energies and properties. The redn. in CPU time obtained with the RI approxn. increases rapidly with basis set size. For the cc-pVQZ basis an acceleration by a factor of up to 170 is obsd.
- 96Boys, S. F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553– 566, DOI: 10.1080/0026897700010156196https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXht1alt7fM&md5=c4f847835a638b4fc449fae688c557e9The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errorsBoys, S. F.; Bernardi, F.Molecular Physics (1970), 19 (4), 553-566CODEN: MOPHAM; ISSN:0026-8976. (Taylor & Francis Ltd.)A new direct difference method for the computation of mol. interactions has been based on a bivariational transcorrelated treatment, together with special methods for the balancing of other errors. It appears that these new features can give a strong redn. in the error of the interaction energy, and they seem to be particularly suitable for computations in the important region near the min. energy. It has been generally accepted that this problem is dominated by unresolved difficulties and the relation of the new methods of these apparent difficulties is analyzed here.
- 97Gdanitz, R. J. An accurate interaction potential for neon dimer (Ne2). Chem. Phys. Lett. 2001, 348, 67– 74, DOI: 10.1016/S0009-2614(01)01088-0There is no corresponding record for this reference.
- 98Perdew, J. P. Proceedings of the 75. WE-Heraeus-Seminar and 21st Annual International Symposium on Electronic Structure of Solids; Akademie Verlag: Berlin, 1991; p 11.There is no corresponding record for this reference.
- 99Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 1992, 46, 6671– 6687, DOI: 10.1103/PhysRevB.46.667199https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK38XlvFyks7c%253D&md5=e52b617f523b1a8c545707544ce01d05Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlationPerdew, John P.; Chevary, J. A.; Vosko, S. H.; Jackson, Koblar A.; Pederson, Mark R.; Singh, D. J.; Fiolhais, CarlosPhysical Review B: Condensed Matter and Materials Physics (1992), 46 (11), 6671-87CODEN: PRBMDO; ISSN:0163-1829.Generalized gradient approxns. (GGA's) seek to improve upon the accuracy of the local-spin-d. (LSD) approxn. in electronic-structure calcns. Perdew and Wang have developed a GGA based on real-space cutoff of the spurious long-range components of the second-order gradient expansion for the exchange-correlation hole. Authors have found that this d. functional performs well in numerical tests for a variety of systems: Total energies of 30 atoms are highly accurate. Ionization energies and electron affinities are improved in a statistical sense, although significant interconfigurational and interterm errors remain. Accurate atomization energies are found for seven hydrocarbon mols., with a rms error per bond of 0.1 eV, compared with 0.7 eV for the LSD approxn. and 2.4 eV for the Hartree-Fock approxn. For atoms and mols., there is a cancellation of error between d. functionals for exchange and correlation, which is most striking whenever the Hartree-Fock result is furthest from expt. The surprising LSD underestimation of the lattice consts. of Li and Na by 3-4% is cor., and the magnetic ground state of solid Fe is restored. The work function, surface energy (neglecting the long-range contribution), and curvature energy of a metallic surface are all slightly reduced in comparison with LSD. Taking account of the pos. long-range contribution, authors find surface and curvature energies in good agreement with exptl. or exact values. Finally, a way is found to visualize and understand the nonlocality of exchange and correlation, its origins, and its phys. effects.
- 100Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098– 3100, DOI: 10.1103/PhysRevA.38.3098100https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL1cXmtlOhsLo%253D&md5=d4d219c134a5a90f689a8abed04d82ccDensity-functional exchange-energy approximation with correct asymptotic behaviorBecke, A. D.Physical Review A: Atomic, Molecular, and Optical Physics (1988), 38 (6), 3098-100CODEN: PLRAAN; ISSN:0556-2791.Current gradient-cor. d.-functional approxns. for the exchange energies of at. and mol. systems fail to reproduce the correct 1/r asymptotic behavior of the exchange-energy d. A gradient-cor. exchange-energy functional is given with the proper asymptotic limit. This functional, contg. only one parameter, fits the exact Hartree-Fock exchange energies of a wide variety of at. systems with remarkable accuracy, surpassing the performance of previous functionals contg. two parameters or more.
- 101Zhang, Y.; Pan, W.; Yang, W. Describing van der Waals interaction in diatomic molecules with generalized gradient approximations: The role of the exchange functional. J. Chem. Phys. 1997, 107, 7921– 7925, DOI: 10.1063/1.475105101https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXntFGisrY%253D&md5=ed2e0b29f00aa7dedbc4eb461bdbe0b1Describing van der Waals Interaction in diatomic molecules with generalized gradient approximations: the role of the exchange functionalZhang, Yingkai; Pan, Wei; Yang, WeitaoJournal of Chemical Physics (1997), 107 (19), 7921-7925CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)Generalized gradient approxns. have been used to calc. the potential energy curves for six rare gas diat. mols. (He2, Ne2, Ar2, HeNe, ArHe, ArNe). Several generalized gradient approxns. are found to provide a good description of binding in these diat. mols. and show a significant improvement over the local d. approxn. in the prediction of bond lengths and dissocn. energies. It is shown here that the behavior of an exchange functional in the region of small d. and large d. gradient plays a very important role in the ability of the functional to describe this type of van der Waals attraction.
- 102Eshuis, H.; Furche, F. Basis set convergence of molecular correlation energy differences within the random phase approximation. J. Chem. Phys. 2012, 136, 084105, DOI: 10.1063/1.3687005102https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XivVeqtb0%253D&md5=d520f6c11aabd38dab2bc1d6b0ac6862Basis set convergence of molecular correlation energy differences within the random phase approximationEshuis, Henk; Furche, FilippJournal of Chemical Physics (2012), 136 (8), 084105/1-084105/6CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The basis set convergence of energy differences obtained from the RPA to the correlation energy is investigated for a wide range of mol. interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark. The use of very large basis sets (> quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable est. of the binding energy for these systems. Counterpoise cor. results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramol. dispersion, isomerization energies, and reaction energies of small org. mols. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency. (c) 2012 American Institute of Physics.
Supporting Information
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpca.4c05289.
MAEs for all subsets of the ASCDB and GMTKN55 databases for σ↑AXK functionals; weight functions wH and wx for PBE reference orbitals; and spline coefficients for σ↑AXK functionals (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.