1
|
Fauser 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. [PMID: 38466924 DOI: 10.1021/acs.jctc.3c01132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
σ-Functionals belong to the class of Kohn-Sham (KS) correlation functionals based on the adiabatic-connection fluctuation-dissipation theorem and are technically closely related to the random phase approximation (RPA). They have the same computational demand as the latter, with the computational effort of an energy evaluation for both methods being lower than that of a preceding hybrid DFT calculation for typical systems but yield much higher accuracy, reaching chemical accuracy of 1 kcal/mol for quantities such as reactions and transition energies in main group chemistry. In previous work on σ-functionals, rather large Gaussian basis sets have been used. Here, we investigate the actual basis set requirements of σ-functionals and present three setups that employ smaller Gaussian basis sets ranging from quadruple-ζ (QZ) to triple-ζ (TZ) quality and represent a good compromise between accuracy and computational efficiency. Furthermore, we introduce an implementation of σ-functionals based on Slater-type basis sets and present two setups of QZ and TZ quality for this implementation. We test the accuracy of these setups on a large database of various physical properties and types of reactions, as well as equilibrium geometries and vibrational frequencies. As expected, the accuracy of σ-functional calculations becomes somewhat lower with a decreasing basis set size. However, for all setups considered here, calculations with σ-functionals are clearly more accurate than those within the RPA and even more so than those of the conventional KS methods. For the smallest setup using Gaussian-type basis functions and Slater-type basis functions, we introduce a reparametrization that reduces the loss in accuracy due to the basis set error to some extent. A comparison with the range-separated hybrid ωB97X-V and the double hybrid DSD-BLYP-D3 shows that σ functionals outperform in accuracy both of these accurate and, for their class, representative functionals.
Collapse
Affiliation(s)
- Steffen Fauser
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Leon Redeker
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Christian Neiss
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Jannis Erhard
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Martensstr. 1, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
- Erlangen National High Performance Computing Center (NHR@FAU), Martensstr. 1, D-91058 Erlangen, Germany
| |
Collapse
|
2
|
Dhingra D, Shori A, Förster A. Chemically accurate singlet-triplet gaps of organic chromophores and linear acenes by the random phase approximation and σ-functionals. J Chem Phys 2023; 159:194105. [PMID: 37966004 DOI: 10.1063/5.0177528] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Accepted: 10/23/2023] [Indexed: 11/16/2023] Open
Abstract
Predicting the energy differences between different spin-states is challenging for many widely used ab initio electronic structure methods. We here assess the ability of the direct random phase approximation (dRPA), dRPA plus two different screened second-order exchange (SOX) corrections, and σ-functionals to predict adiabatic singlet-triplet gaps. With mean absolute deviations of below 0.1 eV to experimental reference values, independent of the Kohn-Sham starting point, dRPA and σ-functionals accurately predict singlet-triplet gaps of 18 organic chromophores. The addition of SOX corrections to dRPA considerably worsens agreement with experiment, adding to the mounting evidence that dRPA+SOX methods are not generally applicable beyond-RPA methods. Also for a series of linear acene chains with up to ten fused rings, dRPA, and σ-functionals are in excellent agreement with coupled-cluster single double triple reference data. In agreement with advanced multi-reference methods, dRPA@PBE and σ-functional@PBE predict a singlet ground state for all chain lengths, while dRPA@PBE0 and σ-functional@PBE0 predict a triplet ground state for longer acenes. Our work shows dRPA and σ-functionals to be reliable methods for calculating singlet-triplet gaps in aromatic molecules.
Collapse
Affiliation(s)
- Daniella Dhingra
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Arjun Shori
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| | - Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1108, 1081 HZ Amsterdam, The Netherlands
| |
Collapse
|
3
|
Cieśliński D, Tucholska AM, Modrzejewski M. Post-Kohn-Sham Random-Phase Approximation and Correction Terms in the Expectation-Value Coupled-Cluster Formulation. J Chem Theory Comput 2023; 19:6619-6631. [PMID: 37774375 PMCID: PMC10569055 DOI: 10.1021/acs.jctc.3c00496] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Indexed: 10/01/2023]
Abstract
Using expectation-value coupled-cluster theory and many-body perturbation theory (MBPT), we formulate a series of corrections to the post-Kohn-Sham (post-KS) random-phase approximation (RPA) energy. The beyond-RPA terms are of two types: those accounting for the non-Hartree-Fock reference and those introducing the coupled-cluster doubles non-ring contractions. The contributions of the former type, introduced via the semicanonical orbital basis, drastically reduce the binding strength in noncovalent systems. The good accuracy is recovered by the attractive third-order doubles correction referred to as Ec2g. The existing RPA approaches based on KS orbitals neglect most of the proposed corrections but can perform well thanks to error cancellation. The proposed method accounts for every contribution in the state-of-the-art renormalized second-order perturbation theory (rPT2) approach but adds additional terms which initially contribute in the third order of MBPT. The cost of energy evaluation scales as noniterative O ( N 4 ) in the implementation with low-rank tensor decomposition. The numerical tests of the proposed approach demonstrate accurate results for noncovalent dimers of polar molecules and for the challenging many-body noncovalent cluster of CH4···(H2O)20.
Collapse
Affiliation(s)
- Dominik Cieśliński
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| | | | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
4
|
Hummel F. On the Chemical Potential of Many-Body Perturbation Theory in Extended Systems. J Chem Theory Comput 2023; 19:1568-1581. [PMID: 36790901 PMCID: PMC10018744 DOI: 10.1021/acs.jctc.2c01043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
Finite-temperature many-body perturbation theory in the grand-canonical ensemble is fundamental to numerous methods for computing electronic properties at nonzero temperature, such as finite-temperature coupled-cluster. In most applications it is the average number of electrons that is known rather than the chemical potential. Expensive correlation calculations must be repeated iteratively in search for the interacting chemical potential that yields the desired average number of electrons. In extended systems with mobile charges the situation is particular, however. Long-ranged electrostatic forces drive the charges such that the average ratio of negative and positive charges is one for any finite chemical potential. All properties per electron are expected to be virtually independent of the chemical potential, as they are in an electric wire at different voltage potentials. This work shows that per electron, the exchange-correlation free energy and the exchange-correlation grand potential indeed agree in the infinite-size limit. Thus, only one expensive correlation calculation suffices for each system size, sparing the search for the interacting chemical potential. This work also demonstrates the importance of regularizing the Coulomb interaction such that each electron on average interacts only with as many electrons as there are electrons in the simulation, avoiding interactions with periodic images. Numerical calculations of the warm uniform electron gas have been conducted with the Spencer-Alavi regularization employing the finite-temperature Hartree approximation for the self-consistent field and linearized finite-temperature direct-ring coupled-cluster doubles for treating correlation.
Collapse
Affiliation(s)
- Felix Hummel
- Institute for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| |
Collapse
|
5
|
Spadetto E, Philipsen PHT, Förster A, Visscher L. Toward Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy. J Chem Theory Comput 2023; 19:1499-1516. [PMID: 36787494 PMCID: PMC10018742 DOI: 10.1021/acs.jctc.2c01201] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
Pair atomic density fitting (PADF) has been identified as a promising strategy to reduce the scaling with system size of quantum chemical methods for the calculation of the correlation energy like the direct random-phase approximation (RPA) or second-order Møller-Plesset perturbation theory (MP2). PADF can however introduce large errors in correlation energies as the two-electron interaction energy is not guaranteed to be bounded from below. This issue can be partially alleviated by using very large fit sets, but this comes at the price of reduced efficiency and having to deal with near-linear dependencies in the fit set. One posibility is to use global density fitting (DF), but in this work, we introduce an alternative methodology to overcome this problem that preserves the intrinsically favorable scaling of PADF. We first regularize the Fock matrix by projecting out parts of the basis set which gives rise to orbital products that are hard to describe by PADF. After having thus obtained a reliable self-consistent field solution, we then also apply this projector to the orbital coefficient matrix to improve the precision of PADF-MP2 and PADF-RPA. We systematically assess the accuracy of this new approach in a numerical atomic orbital framework using Slater type orbitals (STO) and correlation consistent Gaussian type basis sets up to quintuple-ζ quality for systems with more than 200 atoms. For the small and medium systems in the S66 database we show the maximum deviation of PADF-MP2 and PADF-RPA relative correlation energies to DF-MP2 and DF-RPA reference results to be 0.07 and 0.14 kcal/mol, respectively. When the new projector method is used, the errors only slightly increase for large molecules and also when moderately sized fit sets are used the resulting errors are well under control. Finally, we demonstrate the computational efficiency of our algorithm by calculating the interaction energies of large, non-covalently bound complexes with more than 1000 atoms and 20000 atomic orbitals at the RPA@PBE/CC-pVTZ level of theory.
Collapse
Affiliation(s)
- Edoardo Spadetto
- Software for Chemistry and Materials NV, NL-1081HV Amsterdam, The Netherlands
| | | | - Arno Förster
- Software for Chemistry and Materials NV, NL-1081HV Amsterdam, The Netherlands.,Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
6
|
Förster A. Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies. J Chem Theory Comput 2022; 18:5948-5965. [PMID: 36150190 PMCID: PMC9558381 DOI: 10.1021/acs.jctc.2c00366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
With increasing interelectronic distance, the screening
of the
electron–electron interaction by the presence of other electrons
becomes the dominant source of electron correlation. This effect is
described by the random phase approximation (RPA) which is therefore
a promising method for the calculation of weak interactions. The success
of the RPA relies on the cancellation of errors, which can be traced
back to the violation of the crossing symmetry of the 4-point vertex,
leading to strongly overestimated total correlation energies. By the
addition of second-order screened exchange (SOSEX) to the correlation
energy, this issue is substantially reduced. In the adiabatic connection
(AC) SOSEX formalism, one of the two electron–electron interaction
lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A
related SOSEX expression in which both electron–electron interaction
lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the
evaluation of this correlation energy expression does not require
an expensive numerical frequency integration and is therefore advantageous
from a computational perspective. We compare the accuracy of the statically
screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical
reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for
charged excitations and noncovalent interactions where they lead to
major improvements over RPA.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV, Amsterdam, The Netherlands
| |
Collapse
|
7
|
Schäfer T, Daelman N, López N. Cerium Oxides without U: The Role of Many-Electron Correlation. J Phys Chem Lett 2021; 12:6277-6283. [PMID: 34212726 PMCID: PMC8397342 DOI: 10.1021/acs.jpclett.1c01589] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 06/29/2021] [Indexed: 05/14/2023]
Abstract
Electron transfer with changing occupation in the 4f subshell poses a considerable challenge for quantitative predictions in quantum chemistry. Using the example of cerium oxide, we identify the main deficiencies of common parameter-dependent one-electron approaches, such as density functional theory (DFT) with a Hubbard correction, or hybrid functionals. As a response, we present the first benchmark of ab initio many-electron theory for electron transfer energies and lattice parameters under periodic boundary conditions. We show that the direct random phase approximation clearly outperforms all DFT variations. From this foundation, we, then, systematically improve even further. Periodic second-order Møller-Plesset perturbation theory meanwhile manages to recover standard hybrid functional values. Using these approaches to eliminate parameter bias allows for highly accurate benchmarks of strongly correlated materials, the reliable assessment of various density functionals, and functional fitting via machine-learning.
Collapse
Affiliation(s)
- Tobias Schäfer
- Institute
for Theoretical Physics, TU Wien, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Nathan Daelman
- Institute
of Chemical Research of Catalonia, The Barcelona
Institute of Science and Technology, 43007 Tarragona, Spain
| | - Núria López
- Institute
of Chemical Research of Catalonia, The Barcelona
Institute of Science and Technology, 43007 Tarragona, Spain
| |
Collapse
|
8
|
Rishi V, Perera A, Bartlett RJ. A route to improving RPA excitation energies through its connection to equation-of-motion coupled cluster theory. J Chem Phys 2020; 153:234101. [DOI: 10.1063/5.0023862] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Varun Rishi
- Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
9
|
Yu F, Wang Y. Dual‐hybrid direct random phase approximation and second‐order screened exchange with nonlocal van der Waals correlations for noncovalent interactions. J Comput Chem 2020; 41:1018-1025. [DOI: 10.1002/jcc.26149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Revised: 01/05/2020] [Accepted: 01/06/2020] [Indexed: 11/09/2022]
Affiliation(s)
- Feng Yu
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| | - Yaoting Wang
- Department of Physics, School of ScienceXi'an Technological University Xi'an Shaanxi China
| |
Collapse
|
10
|
Klimeš J, Tew DP. Efficient and accurate description of adsorption in zeolites. J Chem Phys 2019; 151:234108. [PMID: 31864262 DOI: 10.1063/1.5123425] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.
Collapse
Affiliation(s)
- Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - David P Tew
- Max-Planck-Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
11
|
Ramberger B, Sukurma Z, Schäfer T, Kresse G. RPA natural orbitals and their application to post-Hartree-Fock electronic structure methods. J Chem Phys 2019; 151:214106. [DOI: 10.1063/1.5128415] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Benjamin Ramberger
- Faculty of Physics and Center for Computational Materials Sciences, University of Vienna, Sensengasse 8/12, 1090 Vienna, Austria
| | - Zoran Sukurma
- Faculty of Physics and Center for Computational Materials Sciences, University of Vienna, Sensengasse 8/12, 1090 Vienna, Austria
| | - Tobias Schäfer
- Institute for Theoretical Physics, Technical University of Vienna, Wiedner Hauptstraße 8-10/136, 1040 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics and Center for Computational Materials Sciences, University of Vienna, Sensengasse 8/12, 1090 Vienna, Austria
| |
Collapse
|