1
|
Kumar VB, Śmiga S, Grabowski I. A Critical Evaluation of the Hybrid KS DFT Functionals Based on the KS Exchange-Correlation Potentials. J Phys Chem Lett 2024:10219-10229. [PMID: 39356205 DOI: 10.1021/acs.jpclett.4c01979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/03/2024]
Abstract
We have developed a critical methodology for the evaluation of the quality of hybrid exchange-correlation (XC) density functional approximations (DFAs) based on very fundamental quantities, i.e., Kohn-Sham (KS) XC potentials, self-consistent electron densities, first ionization potentials (IPs), and total energies. Since the XC potentials, the primary objects in the current study, are not directly accessible for the hybrids, we calculate them by inverting the KS electron densities. Utilizing this methodology, we tested 155 hybrid DFAs available in the LIBXC library using FCI and CCSD(T) methods as a reference. We have found that a group of functionals produces very decent XC potentials, mainly those with a large mixture of Hartree-Fock exchange. Moreover, the value of IP strongly depends on the XC potential quality. On the other hand, we show that the XC energy is dominated by functional-driven error, which in some cases leads to substantial errors in electronic densities. The study shows new directions for constructing more accurate XC functionals within the KS-DFT framework.
Collapse
Affiliation(s)
- Vignesh Balaji Kumar
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| | - Ireneusz Grabowski
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudzia̧dzka 5, 87-100 Toruń, Poland
| |
Collapse
|
2
|
Mester D, Kállay M. Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level. J Chem Theory Comput 2023; 19:3982-3995. [PMID: 37326360 PMCID: PMC10339736 DOI: 10.1021/acs.jctc.3c00363] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Indexed: 06/17/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.
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
- ELKH-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
- ELKH-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
|
3
|
Śmiga S, Della Sala F, Gori-Giorgi P, Fabiano E. Self-Consistent Implementation of Kohn-Sham Adiabatic Connection Models with Improved Treatment of the Strong-Interaction Limit. J Chem Theory Comput 2022; 18:5936-5947. [PMID: 36094908 PMCID: PMC9558377 DOI: 10.1021/acs.jctc.2c00352] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Adiabatic connection
models (ACMs), which interpolate between the
limits of weak and strong interaction, are powerful tools to build
accurate exchange–correlation functionals. If the exact weak-interaction
expansion from the second-order perturbation theory is included, a
self-consistent implementation of these functionals is challenging
and still absent in the literature. In this work, we fill this gap
by presenting a fully self-consistent-field (SCF) implementation of
some popular ACM functionals. While using second-order perturbation
theory at weak interactions, we have also introduced new generalized
gradient approximations (GGAs), beyond the usual point-charge-plus-continuum
model, for the first two leading terms at strong interactions, which
are crucial to ensure robustness and reliability. We then assess the
SCF–ACM functionals for molecular systems and for prototypical
strong-correlation problems. We find that they perform well for both
the total energy and the electronic density and that the impact of
SCF orbitals is directly connected to the accuracy of the ACM functional
form. For the H2 dissociation, the SCF–ACM functionals
yield significant improvements with respect to standard functionals
also thanks to the use of the new GGAs for the strong-coupling functionals.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, ul. Grudziądzka 5, 87-100 Toruń, Poland
| | - Fabio Della Sala
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| | - Paola Gori-Giorgi
- Department of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Campus Unisalento, Lecce, Via Monteroni 73100, Italy.,Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano, Lecce 73010, Italy
| |
Collapse
|
4
|
Pernal K, Hapka M. Range‐separated multiconfigurational density functional theory methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1566] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Michał Hapka
- Lodz University of Technology Institute of Physics Lodz Poland
- Faculty of Chemistry University of Warsaw Warsaw Poland
| |
Collapse
|
5
|
Modrzejewski M, Yourdkhani S, Śmiga S, Klimeš J. Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage. J Chem Theory Comput 2021; 17:804-817. [PMID: 33445879 DOI: 10.1021/acs.jctc.0c00966] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.
Collapse
Affiliation(s)
- Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Pasteura 1, Poland.,Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Sirous Yourdkhani
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziądzka 5, 87-100 Toruń, Poland
| | - Jiří Klimeš
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, CZ-12116 Prague 2, Czech Republic
| |
Collapse
|
6
|
Jana S, Śmiga S, Constantin LA, Samal P. Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms. J Chem Theory Comput 2020; 16:7413-7430. [PMID: 33205659 PMCID: PMC7735712 DOI: 10.1021/acs.jctc.0c00823] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Indexed: 11/28/2022]
Abstract
Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.
Collapse
Affiliation(s)
- Subrata Jana
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Szymon Śmiga
- Institute
of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Consiglio
Nazionale delle Ricerche CNR-NANO, Istituto
di Nanoscienze, 41125 Modena, Italy
| | - Prasanjit Samal
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
7
|
Jana S, Patra A, Śmiga S, Constantin LA, Samal P. Insights from the density functional performance of water and water–solid interactions: SCAN in relation to other meta-GGAs. J Chem Phys 2020; 153:214116. [DOI: 10.1063/5.0028821] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Istituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, Italy
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
8
|
Śmiga S, Constantin LA. Modified Interaction-Strength Interpolation Method as an Important Step toward Self-Consistent Calculations. J Chem Theory Comput 2020; 16:4983-4992. [PMID: 32559078 PMCID: PMC7588043 DOI: 10.1021/acs.jctc.0c00328] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The modified point charge plus continuum (mPC) model [ConstantinL. A.; 2019, 99, 085117] solves
the important failures of the original counterpart, namely, the divergences
when the reduced gradient of the density is large, such as in the
tail of the density and in quasi-dimensional density regimes. The
mPC allows us to define a modified interaction-strength interpolation
(mISI) method inheriting these good features, which are important
steps toward the full self-consistent treatment. Here, we provide
an assessment of mISI for molecular systems (i.e.,
considering thermochemistry properties, correlation energies, vertical
ionization potentials, and several noncovalent interactions), harmonium
atoms, and functional derivatives in the strong-interaction limit.
For all our tests, mISI provides a systematic improvement over the
original ISI method. Semilocal approximations of the second-order
Görling–Levy (GL2) perturbation theory are also considered
in the mISI method, showing considerable worsening of the results.
Possible further development of mISI is briefly discussed.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, 87-100 Toruń, Poland
| | - Lucian A Constantin
- Consiglio Nazionale delle Ricerche CNR-NANO, Istituto di Nanoscienze, 41125 Modena, Italy
| |
Collapse
|
9
|
Śmiga S, Constantin LA. Unveiling the Physics Behind Hybrid Functionals. J Phys Chem A 2020; 124:5606-5614. [PMID: 32551627 PMCID: PMC7590981 DOI: 10.1021/acs.jpca.0c04156] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2020] [Revised: 06/16/2020] [Indexed: 11/30/2022]
Abstract
We show that accurate exchange-correlation hybrid functionals give very physically optimized effective-correlation potentials, capable of correctly describing the quantum oscillations of atoms and molecules. Based on this analysis and on understanding the error cancellation between semilocal exchange and correlation functionals, we propose a very simple, semilocal correlation potential model compatible with the exact exchange of density functional theory, which performs remarkably well for charge densities and orbital energies.
Collapse
Affiliation(s)
- Szymon Śmiga
- Institute
of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Center
for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano (LE), Italy
- Istituto
di Nanoscienze, Consiglio Nazionale delle
Ricerche CNR-NANO, 41125 Modena, Italy
| |
Collapse
|
10
|
Śmiga S, Marusiak V, Grabowski I, Fabiano E. The ab initio density functional theory applied for spin-polarized calculations. J Chem Phys 2020; 152:054109. [DOI: 10.1063/1.5128933] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, 87-100 Torun, Poland
| | - Volodymyr Marusiak
- Department of Theoretical Physics and Astronomy, I. I. Mechnikov Odessa National University, Dvoryanskaya 2, Odessa 65082, Ukraine
| | - Ireneusz Grabowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, 87-100 Torun, Poland
| | - Eduardo Fabiano
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| |
Collapse
|
11
|
Jana S, Patra A, Constantin LA, Samal P. Screened range-separated hybrid by balancing the compact and slowly varying density regimes: Satisfaction of local density linear response. J Chem Phys 2020; 152:044111. [PMID: 32007058 DOI: 10.1063/1.5131530] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Due to their quantitative accuracy and ability to solve several difficulties, screened range-separated hybrid exchange-correlation functionals are now a standard approach for ab initio simulation of condensed matter systems. However, the screened range-separated hybrid functionals proposed so far are biased either toward compact or slowly varying densities. In this paper, we propose a screened range-separated hybrid functional, named HSEint, which can well describe these density regimes, achieving good accuracy for both molecular and solid-state systems. The semilocal part of the proposed functional is based on the PBEint generalized gradient approximation [E. Fabiano et al., Phys. Rev. B 82, 113104 (2010)], constructed for hybrid interfaces. To improve the functional performance, we employ exact or nearly exact constraints in the construction of range-separated hybrid functional, such as recovering of the local density linear response and semiclassical atom linear response.
Collapse
Affiliation(s)
- Subrata Jana
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Abhilash Patra
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Lucian A Constantin
- Center for Biomolecular Nanotechnologies @UNILE, Istituto Italiano di Tecnologia, Via Barsanti, I-73010 Arnesano, Italy
| | - Prasanjit Samal
- School of Physical Sciences, National Institute of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|