1
|
Mester D, Kállay M. Basis Set Limit of CCSD(T) Energies: Explicit Correlation Versus Density-Based Basis-Set Correction. J Chem Theory Comput 2023; 19:8210-8222. [PMID: 37950703 PMCID: PMC10688194 DOI: 10.1021/acs.jctc.3c00979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 10/21/2023] [Accepted: 10/23/2023] [Indexed: 11/13/2023]
Abstract
A thorough comparison is carried out for explicitly correlated and density-based basis-set correction approaches, which were primarily developed to mitigate the basis-set incompleteness error of wave function methods. An efficient implementation of the density-based scheme is also presented, utilizing the density-fitting approximation. The performance of these approaches is comprehensively tested for the second-order Møller-Plesset (MP2), coupled-cluster singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] methods with respect to the corresponding complete basis set references. It is demonstrated that the density-based correction together with complementary auxiliary basis set (CABS)-corrected Hartree-Fock energies is highly robust and effectively reduces the error of the standard approaches; however, it does not outperform the corresponding explicitly correlated methods. Nevertheless, what still makes the density-corrected CCSD and CCSD(T) methods competitive is that their computational costs are roughly half of those of the corresponding explicitly correlated variants. Additionally, an incremental approach for standard CCSD and CCSD(T) is introduced. In this simple scheme, the total energies are corrected with the CABS correction and explicitly correlated MP2 contributions. As demonstrated, the resulting methods yield surprisingly good results, below 1 kcal/mol for thermochemical properties even with a double-ζ basis, while their computational expenses are practically identical to those of the density-based basis-set correction approaches.
Collapse
Affiliation(s)
- Dávid Mester
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN-BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
2
|
Hellmann L, Tölle J, Niemeyer N, Neugebauer J. Automated Generation of Optimized Auxiliary Basis Sets for Long-Range-Corrected TDDFT Using the Cholesky Decomposition. J Chem Theory Comput 2022; 18:2959-2974. [PMID: 35446029 DOI: 10.1021/acs.jctc.2c00131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Range-separated hybrid functionals making use of a smooth separation of the Coulomb operator in terms of the error function and its complement have proven to be a valuable tool for improving Kohn-Sham density functional theory (DFT) calculations. This holds in particular for obtaining accurate excitation energies from linear-response time-dependent DFT. Evaluating the long-range exchange contributions represents one of the most time-consuming tasks in such calculations. Prefitted auxiliary basis sets can be employed to speed up this step. Here, we present a way to generate auxiliary basis sets optimized to fit the long-range exchange contributions only, contrary to the common optimization strategies on the basis of the full Coulomb operator. For this purpose, we use the atomic Cholesky decomposition technique. The basis sets are generated on-the-fly using the specific range-separation parameter defined in the exchange-correlation functional. We obtain excitation energies and oscillator strengths which are of similar or better accuracy than those obtained with conventional resolution-of-the-identity auxiliary basis sets while drastically reducing the number of auxiliary functions required. This is demonstrated for the QUESTDB#5 benchmark set. In addition, we outline the benefits of this approach in sequences of calculations employing varying range-separation parameters, as is the case in the optimally tuned range-separation strategy. Finally, we illustrate the efficiency of this approach for real-world examples, namely, a chlorophyll tetramer from photosystem II and a carotenoid-porphyrin-C60 triad.
Collapse
Affiliation(s)
- Lars Hellmann
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany.,Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany.,Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Niklas Niemeyer
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany.,Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany.,Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
3
|
Stein F, Hutter J. Double-hybrid density functionals for the condensed phase: Gradients, stress tensor, and auxiliary-density matrix method acceleration. J Chem Phys 2022; 156:074107. [DOI: 10.1063/5.0082327] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich 8057, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, Zurich 8057, Switzerland
| |
Collapse
|
4
|
Mester D, Kállay M. Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach. J Chem Theory Comput 2022; 18:865-882. [PMID: 35023739 PMCID: PMC8830052 DOI: 10.1021/acs.jctc.1c01100] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Our second-order algebraic-diagrammatic construction [ADC(2)]-based double-hybrid (DH) ansatz (J. Chem. Theory Comput. 2019, 15, 4440. DOI: 10.1021/acs.jctc.9b00391) is combined with range-separation techniques. In the present scheme, both the exchange and the correlation contributions are range-separated, while spin-scaling approaches are also applied. The new methods are thoroughly tested for the most popular benchmark sets including 250 singlet and 156 triplet excitations, as well as 80 oscillator strengths. It is demonstrated that the range separation for the correlation contributions is highly recommended for both the genuine and the ADC(2)-based DH approaches. Our results show that the latter scheme slightly but consistently outperforms the former one for single excitation dominated transitions. Furthermore, states with larger fractions of double excitations are assessed as well, and challenging charge-transfer excitations are also discussed, where the recently proposed spin-scaled long-range corrected DHs fail. The suggested iterative fourth-power scaling RS-PBE-P86/SOS-ADC(2) method, using only three adjustable parameters, provides the most robust and accurate excitation energies within the DH theory. In addition, the relative error of the oscillator strengths is reduced by 65% compared to the best genuine DH functionals.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
5
|
Mester D, Kállay M. Spin-Scaled Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:4211-4224. [PMID: 34152771 PMCID: PMC8280718 DOI: 10.1021/acs.jctc.1c00422] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Indexed: 11/28/2022]
Abstract
Our recently presented range-separated (RS) double-hybrid (DH) time-dependent density functional approach [J. Chem. Theory Comput. 17, 927 (2021)] is combined with spin-scaling techniques. The proposed spin-component-scaled (SCS) and scaled-opposite-spin (SOS) variants are thoroughly tested for almost 500 excitations including the most challenging types. This comprehensive study provides useful information not only about the new approaches but also about the most prominent methods in the DH class. The benchmark calculations confirm the robustness of the RS-DH ansatz, while several tendencies and deficiencies are pointed out for the existing functionals. Our results show that the SCS variant consistently improves the results, while the SOS variant preserves the benefits of the original RS-DH method reducing its computational expenses. It is also demonstrated that, besides our approaches, only the nonempirical functionals provide balanced performance for general applications, while particular methods are only suggested for certain types of excitations.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
6
|
Mester D, Kállay M. A Simple Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:927-942. [PMID: 33400872 PMCID: PMC7884002 DOI: 10.1021/acs.jctc.0c01135] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Indexed: 01/12/2023]
Abstract
A simple and robust range-separated (RS) double-hybrid (DH) time-dependent density functional approach is presented for the accurate calculation of excitation energies of molecules within the Tamm-Dancoff approximation. The scheme can be considered as an excited-state extension of the ansatz proposed by Toulouse and co-workers [J. Chem. Phys. 2018, 148, 164105], which is based on the two-parameter decomposition of the Coulomb potential, for which both the exchange and correlation contributions are range-separated. A flexible and efficient implementation of the new scheme is also presented, which facilitates its extension to any combination of exchange and correlation functionals. The performance of the new approximation is tested for singlet excitations on several benchmark compilations and thoroughly compared to that of representative DH, RS hybrid, and RS DH functionals. The one-electron basis set dependence and computation times are also assessed. Our results show that the new approach improves on standard DHs in most cases, and it can provide a more robust and accurate alternative. In addition, on average, it noticeably surpasses the existing RS hybrid and RS DH functionals.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| |
Collapse
|
7
|
Stein F, Hutter J, Rybkin VV. Double-Hybrid DFT Functionals for the Condensed Phase: Gaussian and Plane Waves Implementation and Evaluation. Molecules 2020; 25:E5174. [PMID: 33172070 PMCID: PMC7664425 DOI: 10.3390/molecules25215174] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2020] [Revised: 10/23/2020] [Accepted: 10/29/2020] [Indexed: 11/16/2022] Open
Abstract
Intermolecular interactions play an important role for the understanding of catalysis, biochemistry and pharmacy. Double-hybrid density functionals (DHDFs) combine the proper treatment of short-range interactions of common density functionals with the correct description of long-range interactions of wave-function correlation methods. Up to now, there are only a few benchmark studies available examining the performance of DHDFs in condensed phase. We studied the performance of a small but diverse selection of DHDFs implemented within Gaussian and plane waves formalism on cohesive energies of four representative dispersion interaction dominated crystal structures. We found that the PWRB95 and ωB97X-2 functionals provide an excellent description of long-ranged interactions in solids. In addition, we identified numerical issues due to the extreme grid dependence of the underlying density functional for PWRB95. The basis set superposition error (BSSE) and convergence with respect to the super cell size are discussed for two different large basis sets.
Collapse
Affiliation(s)
- Frederick Stein
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| | - Vladimir V. Rybkin
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| |
Collapse
|
8
|
Kreppel 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. [PMID: 32329618 DOI: 10.1021/acs.jctc.9b01294] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A formulation of range-separated random phase approximation (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-separated RPA in general main group thermochemistry, kinetics, and noncovalent interactions. The results show that range-separated RPA performs stably over the broad range of molecular chemistry included in the GMTKN55 set. It improves significantly over semilocal DFT but it is still less accurate than modern dispersion corrected double-hybrid functionals. Furthermore, range-separated RPA shows a faster basis set convergence compared to standard full-range RPA making it a promising applicable approach with only one empirical parameter.
Collapse
Affiliation(s)
- Andrea Kreppel
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Henryk Laqua
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany.,Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
9
|
Kállay M, Nagy PR, Mester D, Rolik Z, Samu G, Csontos J, Csóka J, Szabó PB, Gyevi-Nagy L, Hégely B, Ladjánszki I, Szegedy L, Ladóczki B, Petrov K, Farkas M, Mezei PD, Ganyecz Á. The MRCC program system: Accurate quantum chemistry from water to proteins. J Chem Phys 2020; 152:074107. [DOI: 10.1063/1.5142048] [Citation(s) in RCA: 152] [Impact Index Per Article: 38.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Zoltán Rolik
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Gyula Samu
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csontos
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - P. Bernát Szabó
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - László Gyevi-Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Bence Hégely
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - István Ladjánszki
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Lóránt Szegedy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Bence Ladóczki
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Klára Petrov
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Máté Farkas
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Pál D. Mezei
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Ádám Ganyecz
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|