1
|
Janesko BG. Multiconfigurational Correlation at DFT + U Cost: On-Site Electron-Electron Interactions Yield a Block-Localized Configuration Interaction Hamiltonian. J Phys Chem A 2024; 128:5077-5087. [PMID: 38878060 DOI: 10.1021/acs.jpca.4c02326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2024]
Abstract
This work presents a first-principles wavefunction-in-DFT approach based on the Hubbard density functional theory (DFT) + U method. This approach begins with the standard DFT reference system of noninteracting electrons and introduces an electron-electron interaction projected onto DFT+U-type atomic states. The reference system's configuration interaction Hamiltonian is block-localized to these states and can be expressed in terms of state occupation numbers, state self-energies (which correspond to unscreened Hubbard U values), and the promotion energies of doubly excited Slater determinants. Simple approximations for the promotion energies provide multiconfigurational correlation energies without requiring explicit orbital localization/transform. Numerical results for fractionally occupied chromium atom, bonded chromium dimer, dissociating covalent bonds, and large active spaces show that the approach provides beyond-zero-sum accuracy at computational cost comparable to standard DFT+U.
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas 76129, United States
| |
Collapse
|
2
|
Becke AD. A remarkably simple dispersion damping scheme and the DH24 double hybrid density functional. J Chem Phys 2024; 160:204118. [PMID: 38818895 DOI: 10.1063/5.0207682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2024] [Accepted: 04/25/2024] [Indexed: 06/01/2024] Open
Abstract
In recent papers, Becke et al. [J. Chem. Phys. 158, 151103 (2023)] and then Becke [J. Chem. Phys. 159, 241101 (2023)] have developed a novel double hybrid density functional, "DH23," whose terms are based on good local physics. Its 12 coefficients are trained on the GMTKN55 (general main-group thermochemistry, kinetics, and noncovalent interactions) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)]. The lowest GMTKN55 "WTMAD2" error to date for any hybrid or double hybrid density functional was obtained (1.73 kcal/mol for the revDH23 variant). Here, we simplify DH23 by introducing a dispersion damping scheme involving atomic numbers only and one global parameter. The resulting new functional, "DH24," performs as well as its predecessors.
Collapse
Affiliation(s)
- Axel D Becke
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, P.O. Box 15000, Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
3
|
Rask AE, Li L, Zimmerman PM. Kohn-Sham Density in a Slater Orbital Basis Set. J Phys Chem A 2024; 128:3194-3204. [PMID: 38602291 DOI: 10.1021/acs.jpca.3c08303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
Finite, atom-centered Slater basis sets are used to determine approximate Kohn-Sham molecular orbitals. This is achieved by minimizing the kinetic energy plus the sum-squared difference between the Kohn-Sham density and the full configuration interaction density. As a result of the finite basis, a weight factor is introduced to balance the two minimization components. Results herein show that this can be done systematically, without sensitive dependence on the choice of scaling factor. In addition, the algorithm is applied to the LiH diatomic for fractional electron counts, where stretching the bond introduces significant reorganization of the electron density. The analysis will show the correct KS orbital structure and reveal the effects of correlation and electron locality on the KS solutions.
Collapse
Affiliation(s)
- Alan E Rask
- SandboxAQ, 780 High Street, Palo Alto, California 94301, United States
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109, United States
| | - Liying Li
- Department of Mathematics, University of Michigan, 530 Church Street, Ann Arbor, Michigan 48109, United States
- Department of Mathematics, Southern University of Science and Technology, No. 1088 Xueyuan Road, Nanshan District, Shenzhen, 518055, Guangdong, China
| | - Paul M Zimmerman
- Department of Chemistry, University of Michigan, 930 N. University Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
4
|
Becke AD. Doubling down on density-functional theory. J Chem Phys 2023; 159:241101. [PMID: 38146827 DOI: 10.1063/5.0178236] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 11/30/2023] [Indexed: 12/27/2023] Open
Abstract
In a recent paper, Becke et al. [J. Chem. Phys. 158, 151103 (2023)] presented a novel double hybrid density functional, "DH23," whose terms are based on good physics. Its 12 coefficients were trained on the GMTKN55 (general main-group thermochemistry, kinetics, and noncovalent interactions) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)]. The lowest GMTKN55 "WTMAD2" error to date for any hybrid or double hybrid density functional was obtained (1.76 kcal/mol). Here, we make some revisions to DH23 and test its efficacy on reference data beyond GMTKN55, namely, organometallic reaction energies and barrier heights. The results confirm that DH23 is robust outside its training set. In the process, a slightly smaller GMTKN55 WTMAD2 of 1.73 kcal/mol is achieved.
Collapse
Affiliation(s)
- Axel D Becke
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, P.O. Box 15000, Halifax, Nova Scotia B3H 4R2, Canada
| |
Collapse
|
5
|
Fürst S, Kaupp M, Wodyński A. Range-Separated Local Hybrid Functionals with Small Fractional-Charge and Fractional-Spin Errors: Escaping the Zero-Sum Game of DFT Functionals. J Chem Theory Comput 2023. [PMID: 37972297 DOI: 10.1021/acs.jctc.3c00877] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Extending recent developments on strong-correlation (sc) corrections to local hybrid functionals to the recent accurate ωLH22t range-separated local hybrid, a series of highly flexible strong-correlation-corrected range-separated local hybrids (scRSLHs) has been constructed and evaluated. This has required the position-dependent reduction of both short- and long-range exact-exchange admixtures in regions of space characterized by strong static correlations. Using damping procedures provides scRSLHs that retain largely the excellent performance of ωLH22t for weakly correlated situations and, in particular, for accurate quasiparticle energies of a wide variety of systems while reducing dramatically static-correlation errors, e.g., in stretched-bond situations. An additional correction to the local mixing function to reduce delocalization errors in abnormal open-shell situations provides further improvements in thermochemical and kinetic parameters, making scRSLH functionals such as ωLH23tdE or ωLH23tdP promising tools for complex molecular or condensed-phase systems, where low fractional-charge and fractional-spin errors are simultaneously important. The proposed rung 4 functionals thereby largely escape the usual zero-sum game between these two quantities and are expected to open new areas of accurate computations by Kohn-Sham DFT. At the same time, they require essentially no extra computational effort over the underlying ωLH22t functional, which means that their use is only moderately more demanding than that of global, local, or range-separated hybrid functionals.
Collapse
Affiliation(s)
- Susanne Fürst
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Artur Wodyński
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
6
|
Grodner B, Żołek T, Kutner A. Nonaqueous Capillary Electrophoretic Separation of Analogs of (24 R)-1,24-Dihydroxyvitamin D 3 Derivative as Predicted by Quantum Chemical Calculations. Molecules 2023; 28:5055. [PMID: 37446718 DOI: 10.3390/molecules28135055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Revised: 06/16/2023] [Accepted: 06/25/2023] [Indexed: 07/15/2023] Open
Abstract
Nonaqueous capillary electrophoretic (NACE) separation was obtained of analogs of (24R)-1,24-dihydroxyvitamin D3 derivative (calcipotriol) as predicted by quantum chemical calculations supported by the density functional theory (DFT). Among the key electronic properties investigated, absolute values of the dipole polarizability and energy gap between HOMO and LUMO molecular orbitals of the analog molecules differ significantly for particular analogs, and there is a direct relationship with their electrophoretic migration time. These differences and relationships suggest that the structurally related analogs should be separable in the electrostatic field. Indeed, the robust, sensitive, and rapid NACE method was first developed for the identification and determination of the anticancer analog of calcipotriol (coded PRI-2205) and its process-related impurities (coded PRI-2201, PRI-2203, and PRI-2204) in organic and aqueous biological solutions. The direct relation between the calculated electronic properties of the analogs and the experimental electrophoretic migration time could be a promising prospect for theoretically predicting the electrophoretic separations.
Collapse
Affiliation(s)
- Błażej Grodner
- Department of Biochemistry and Pharmacogenomics, Faculty of Pharmacy, Medical University of Warsaw, 1 Banacha, 02-097 Warsaw, Poland
| | - Teresa Żołek
- Department of Organic and Physical Chemistry, Faculty of Pharmacy, Medical University of Warsaw, 1 Banacha, 02-097 Warsaw, Poland
| | - Andrzej Kutner
- Department of Drug Chemistry, Faculty of Pharmacy, Medical University of Warsaw, 1 Banacha, 02-097 Warsaw, Poland
| |
Collapse
|
7
|
Becke AD, Santra G, Martin JML. A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database. J Chem Phys 2023; 158:2882268. [PMID: 37094004 DOI: 10.1063/5.0141238] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Accepted: 03/20/2023] [Indexed: 04/26/2023] Open
Abstract
In two recent papers [A. D. Becke, J. Chem. Phys. 156, 214101 (2022) and A. D. Becke, J. Chem. Phys. 157, 234102 (2022)], we compared two Kohn-Sham density functionals based on physical modeling and theory with the best density-functional power-series fits in the literature. The best error statistics reported to date for a hybrid functional on the general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)] were obtained. In the present work, additional second-order perturbation-theory terms are considered. The result is a 12-parameter double-hybrid density functional with the lowest GMTKN55 WTMAD2 "weighted total mean absolute deviation" error (1.76 kcal/mol) yet seen for any hybrid or double-hybrid density-functional approximation. We call it "DH23."
Collapse
Affiliation(s)
- Axel D Becke
- Department of Chemistry, Dalhousie University, 6274 Coburg Road, P.O. Box 15000, Halifax, Nova Scotia B3H 4R2, Canada
| | - Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Rehovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Rehovot, Israel
| |
Collapse
|