1
|
Panchagnula K, Graf D, Johnson ER, Thom AJW. Targeting spectroscopic accuracy for dispersion bound systems from ab initio techniques: Translational eigenstates of Ne@C70. J Chem Phys 2024; 161:054308. [PMID: 39092939 DOI: 10.1063/5.0223298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2024] [Accepted: 07/17/2024] [Indexed: 08/04/2024] Open
Abstract
We investigate the endofullerene system Ne@C70 by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. This is constructed from electronic structure calculations from a plethora of methods, including MP2, SCS-MP2, SOS-MP2, RPA@PBE, and C(HF)-RPA, which were previously used for He@C60 in Panchagnula et al. [J. Chem. Phys. 160, 104303 (2024)], alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C60 to C70 introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the electronic structure methods. The nuclear Hamiltonian is diagonalized using a symmetrized double minimum basis set outlined in Panchagnula and Thom [J. Chem. Phys. 159, 164308 (2023)], with translational energies having error bars ±1 and ±2 cm-1. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well and different translational eigenspectra, which also differ from the Lennard-Jones (LJ) PES given in Mandziuk and Bačić [J. Chem. Phys. 101, 2126-2140 (1994)]. We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES nor recreate the eigenstates, and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C60, due to the lack of concordance between all electronic structure methods, we require more experimental data in order to properly validate the choice.
Collapse
Affiliation(s)
- K Panchagnula
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - D Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - E R Johnson
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Department of Chemistry, Dalhousie University, 6243 Alumni Crescent, Halifax, Nova Scotia B3H 4R2, Canada
| | - A J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
2
|
Kaplan AD, Shahi C, Sah RK, Bhetwal P, Kanungo B, Gavini V, Perdew JP. How Does HF-DFT Achieve Chemical Accuracy for Water Clusters? J Chem Theory Comput 2024; 20:5517-5527. [PMID: 38937987 DOI: 10.1021/acs.jctc.4c00560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024]
Abstract
Bolstered by recent calculations of exact functional-driven errors (FEs) and density-driven errors (DEs) of semilocal density functionals in the water dimer binding energy [Kanungo, B. J. Phys. Chem. Lett. 2024, 15, 323-328], we investigate approximate FEs and DEs in neutral water clusters containing up to 20 monomers, charged water clusters, and alkali- and halide-water clusters. Our proxy for the exact density is r2SCAN 50, a 50% global hybrid of exact exchange with r2SCAN, which may be less correct than r2SCAN for the compact water monomer but importantly more correct for long-range electron transfers in the noncompact water clusters. We show that SCAN makes substantially larger FEs for neutral water clusters than r2SCAN, while both make essentially the same DEs. Unlike the case for barrier heights, these FEs are small in a relative sense and become large in an absolute sense only due to an increase in cluster size. SCAN@HF, short for SCAN evaluated on the Hartree-Fock (HF) density, produces a cancellation of errors that makes it chemically accurate for predicting the absolute binding energies of water clusters. Likewise, adding a long-range dispersion correction to r2SCAN@HF, as in the composite method HF-r2SCAN-DC4, makes its FE more negative than in r2SCAN@HF, permitting a near-perfect cancellation of FE and DE. r2SCAN by itself (and even more so, r2SCAN evaluated on the r2SCAN 50 density), is almost perfect for the energy differences between water hexamers, and thus probably also for liquid water away from the boiling point. Thus, the accuracy of composite methods like SCAN@HF and HF-r2SCAN-DC4 is not due to the HF density being closer to the exact density, but to a compensation of errors from its greater degree of localization. We also give an argument for the approximate reliability of this unconventional error cancellation for diverse molecular properties. Finally, we confirm this unconventional error cancellation for the SCAN description of the water trimer via Kohn-Sham inversion of the CCSD(T) density.
Collapse
Affiliation(s)
- Aaron D Kaplan
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Chandra Shahi
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - Raj K Sah
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Pradeep Bhetwal
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Bikash Kanungo
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Vikram Gavini
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Materials Science and Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - John P Perdew
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
3
|
Singh Y, Peralta JE, Jackson KA. The rise and fall of stretched bond errors: Extending the analysis of Perdew-Zunger self-interaction corrections of reaction barrier heights beyond the LSDA. J Chem Phys 2024; 160:124105. [PMID: 38526103 DOI: 10.1063/5.0179261] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Accepted: 03/07/2024] [Indexed: 03/26/2024] Open
Abstract
Incorporating self-interaction corrections (SIC) significantly improves chemical reaction barrier height predictions made using density functional theory methods. We present a detailed orbital-by-orbital analysis of these corrections for three semi-local density functional approximations (DFAs) situated on the three lowest rungs of Jacob's ladder of approximations. The analysis is based on Fermi-Löwdin Orbital Self-Interaction Correction (FLOSIC) calculations performed at several steps along the reaction pathway from the reactants (R) to the transition state (TS) to the products (P) for four representative reactions selected from the BH76 benchmark set. For all three functionals, the major contribution to self-interaction corrections of the barrier heights can be traced to stretched bond orbitals that develop near the TS configuration. The magnitude of the ratio of the self-exchange-correlation energy to the self-Hartree energy (XC/H) for a given orbital is introduced as an indicator of one-electron self-interaction error. XC/H = 1.0 implies that an orbital's self-exchange-correlation energy exactly cancels its self-Hartree energy and that the orbital, therefore, makes no contribution to the SIC in the FLOSIC scheme. For the practical DFAs studied here, XC/H spans a range of values. The largest values are obtained for stretched or strongly lobed orbitals. We show that significant differences in XC/H for corresponding orbitals in the R, TS, and P configurations can be used to identify the major contributors to the SIC of barrier heights and reaction energies. Based on such comparisons, we suggest that barrier height predictions made using the strongly constrained and appropriately normed meta-generalized gradient approximation may have attained the best accuracy possible for a semi-local functional using the Perdew-Zunger SIC approach.
Collapse
Affiliation(s)
- Yashpal Singh
- Department of Physics, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Juan E Peralta
- Department of Physics and Science of Advanced Materials PhD Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Koblar A Jackson
- Department of Physics and Science of Advanced Materials PhD Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| |
Collapse
|
4
|
Kanungo B, Kaplan AD, Shahi C, Gavini V, Perdew JP. Unconventional Error Cancellation Explains the Success of Hartree-Fock Density Functional Theory for Barrier Heights. J Phys Chem Lett 2024; 15:323-328. [PMID: 38170179 DOI: 10.1021/acs.jpclett.3c03088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally efficient semilocal approximations for the exchange-correlation energy. The accuracy of a semilocal density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree-Fock electron densities, which has been known for ∼30 years. The conventional explanation is that the Hartree-Fock theory yields the more accurate density. This work presents a benchmark Kohn-Sham inversion of accurate coupled-cluster densities for the reaction H2 + F → HHF → H + HF and finds a strong, understandable cancellation between positive (excessively overcorrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree-Fock density functional theory. This confirms earlier conclusions (Kaplan, A. D., et al. J. Chem. Theory Comput. 2023, 19, 532-543) based on 76 barrier heights and three less reliable, but less expensive, fully nonlocal density functional proxies for the exact density.
Collapse
Affiliation(s)
- Bikash Kanungo
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Aaron D Kaplan
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Chandra Shahi
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - Vikram Gavini
- Department of Mechanical Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - John P Perdew
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| |
Collapse
|
5
|
Barneschi L, Kaliakin D, Huix-Rotllant M, Ferré N, Filatov Gulak M, Olivucci M. Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods. J Chem Theory Comput 2023; 19:8189-8200. [PMID: 37937990 DOI: 10.1021/acs.jctc.3c00879] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.
Collapse
Affiliation(s)
- Leonardo Barneschi
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, I-53100 Siena, Italy
| | - Danil Kaliakin
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| | - Miquel Huix-Rotllant
- Aix-Marseille Université, CNRS, Institut Chimie Radicalaire, 13013 Marseille, France
| | - Nicolas Ferré
- Aix-Marseille Université, CNRS, Institut Chimie Radicalaire, 13013 Marseille, France
| | - Michael Filatov Gulak
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Massimo Olivucci
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università di Siena, I-53100 Siena, Italy
- Department of Chemistry, Bowling Green State University, Bowling Green, Ohio 43403, United States
| |
Collapse
|
6
|
Mi W, Luo K, Trickey SB, Pavanello M. Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations. Chem Rev 2023; 123:12039-12104. [PMID: 37870767 DOI: 10.1021/acs.chemrev.2c00758] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.
Collapse
Affiliation(s)
- Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, PR China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, PR China
- International Center of Future Science, Jilin University, Changchun 130012, PR China
| | - Kai Luo
- Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, PR China
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Michele Pavanello
- Department of Physics and Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
7
|
Graf D, Thom AJW. Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost. J Chem Phys 2023; 159:174106. [PMID: 37921249 DOI: 10.1063/5.0168569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| |
Collapse
|
8
|
Graf D, Thom AJW. Simple and Efficient Route toward Improved Energetics within the Framework of Density-Corrected Density Functional Theory. J Chem Theory Comput 2023; 19:5427-5438. [PMID: 37525457 PMCID: PMC10448722 DOI: 10.1021/acs.jctc.3c00441] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2023] [Indexed: 08/02/2023]
Abstract
The crucial step in density-corrected Hartree-Fock density functional theory (DC(HF)-DFT) is to decide whether the density produced by the density functional for a specific calculation is erroneous and, hence, should be replaced by, in this case, the HF density. We introduce an indicator, based on the difference in noninteracting kinetic energies between DFT and HF calculations, to determine when the HF density is the better option. Our kinetic energy indicator directly compares the self-consistent density of the analyzed functional with the HF density, is size-intensive, reliable, and most importantly highly efficient. Moreover, we present a procedure that makes best use of the computed quantities necessary for DC(HF)-DFT by additionally evaluating a related hybrid functional and, in that way, not only "corrects" the density but also the functional itself; we call that procedure corrected Hartree-Fock density functional theory (C(HF)-DFT).
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| | - Alex J. W. Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
9
|
Morgante P, Autschbach J. Density-Corrected Density Functional Theory for Molecular Properties. J Phys Chem Lett 2023:4983-4989. [PMID: 37220345 DOI: 10.1021/acs.jpclett.3c00953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Density-corrected (DC) density functional theory (DFT) has been proposed to overcome difficulties related to the self-interaction error. The procedure uses the Hartree-Fock electron density (matrix) non-self-consistently in conjunction with an approximate functional. DC-DFT has so far mainly been tested for total energy differences, whereas other types of molecular properties have not been evaluated systematically. This work focuses on the performance of DC-DFT for molecular properties, namely, dipole moments, static polarizabilities, and electric field gradients (EFGs) at atomic nuclei. Accurate reference data were generated from coupled-cluster theory to assess the performance of DC and self-consistent DFT calculations for twelve molecules, including diatomics with transition metals. DC-DFT does no harm in dipole moment calculations, but it negatively impacts the polarizability in at least one case. DC-DFT performs well for EFGs, even for the difficult case of CuCl.
Collapse
Affiliation(s)
- Pierpaolo Morgante
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
10
|
Shukla PB, Mishra P, Baruah T, Zope RR, Jackson KA, Johnson JK. How Do Self-Interaction Errors Associated with Stretched Bonds Affect Barrier Height Predictions? J Phys Chem A 2023; 127:1750-1759. [PMID: 36787213 PMCID: PMC10032567 DOI: 10.1021/acs.jpca.2c07894] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
Density functional theory (DFT) suffers from self-interaction errors (SIEs) that generally result in the underestimation of chemical reaction barrier heights. This is commonly attributed to the tendency of density functional approximations to overstabilize delocalized densities that typically occur in the stretched bonds of transition state structures. The Perdew-Zunger self-interaction correction (PZSIC) and locally scaled self-interaction correction (LSIC) improve the prediction of barrier heights of chemical reactions, with LSIC giving better accuracy than PZSIC on average. These methods employ an orbital-by-orbital correction scheme to remove the one-electron SIE. In the context of barrier heights, this allows an analysis of how the self-interaction correction (SIC) for each orbital contributes to the calculated barriers using Fermi-Löwdin orbitals (FLOs). We hypothesize that the SIC contribution to the reaction barrier comes mainly from a limited number of orbitals that are directly involved in bond-breaking and bond-making in the reaction transition state. We call these participant orbitals (POs), in contrast to spectator orbitals (SOs) which are not directly involved in changes to the bonding. Our hypothesis is that ΔETotalSIC ≈ ΔEPOSIC, where ΔETotalSIC is the difference in the SIC corrections for the reactants or products and the transition state. We test this hypothesis for the reaction barriers of the BH76 benchmark set of reactions. We find that the stretched-bond orbitals indeed make the largest individual SIC contributions to the barriers. These contributions increase the barrier heights relative to LSDA, which underpredicts the barrier. However, the full stretched-bond hypothesis does not hold in all cases for either PZSIC or LSIC. There are many cases where the total SIC contribution from the SOs is significant and cannot be ignored. The size of the SIC contribution to the barrier height is a key indicator. A large SIC correction is correlated to a large LSDA error in the barrier, showing that PZSIC properly gives larger corrections when corrections are needed most. A comparison of the performance of PZSIC and LSIC shows that the two methods have similar accuracy for reactions with large LSDA errors, but LSIC is clearly better for reactions with small errors. We trace this to an improved description of reaction energies in LSIC.
Collapse
Affiliation(s)
- Priyanka B Shukla
- Department of Chemical & Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Prakash Mishra
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Tunna Baruah
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, United States
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Rajendra R Zope
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, United States
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Koblar A Jackson
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mount Pleasant, Michigan 48859, United States
| | - J Karl Johnson
- Department of Chemical & Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
11
|
Song S, Vuckovic S, Kim Y, Yu H, Sim E, Burke K. Extending density functional theory with near chemical accuracy beyond pure water. Nat Commun 2023; 14:799. [PMID: 36781855 PMCID: PMC9925738 DOI: 10.1038/s41467-023-36094-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 01/13/2023] [Indexed: 02/15/2023] Open
Abstract
Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r2SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.
Collapse
Affiliation(s)
- Suhwan Song
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea ,grid.266093.80000 0001 0668 7243Department of Chemistry, University of California, Irvine, CA 92697 USA
| | - Stefan Vuckovic
- grid.472716.10000 0004 1758 7362Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy ,grid.12380.380000 0004 1754 9227Departments 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
| | - Youngsam Kim
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea
| | - Hayoung Yu
- grid.15444.300000 0004 0470 5454Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722 Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722, Korea.
| | - Kieron Burke
- grid.266093.80000 0001 0668 7243Department of Chemistry, University of California, Irvine, CA 92697 USA ,grid.266093.80000 0001 0668 7243Departments of Physics & Astronomy, University of California, Irvine, CA 92697 USA
| |
Collapse
|
12
|
Kaplan AD, Shahi C, Bhetwal P, Sah RK, Perdew JP. Understanding Density-Driven Errors for Reaction Barrier Heights. J Chem Theory Comput 2023; 19:532-543. [PMID: 36599075 DOI: 10.1021/acs.jctc.2c00953] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Delocalization errors, such as charge-transfer and some self-interaction errors, plague computationally efficient and otherwise accurate density functional approximations (DFAs). Evaluating a semilocal DFA non-self-consistently on the Hartree-Fock (HF) density is often recommended as a computationally inexpensive remedy for delocalization errors. For sophisticated meta-GGAs like SCAN, this approach can achieve remarkable accuracy. This HF-DFT (also known as DFA@HF) is often presumed to work, when it significantly improves over the DFA, because the HF density is more accurate than the self-consistent DFA density in those cases. By applying the metrics of density-corrected density functional theory (DFT), we show that HF-DFT works for barrier heights by making a localizing charge-transfer error or density overcorrection, thereby producing a somewhat reliable cancellation of density- and functional-driven errors for the energy. A quantitative analysis of the charge-transfer errors in a few randomly selected transition states confirms this trend. We do not have the exact functional and electron densities that would be needed to evaluate the exact density- and functional-driven errors for the large BH76 database of barrier heights. Instead, we have identified and employed three fully nonlocal proxy functionals (SCAN 50% global hybrid, range-separated hybrid LC-ωPBE, and SCAN-FLOSIC) and their self-consistent proxy densities. These functionals are chosen because they yield reasonably accurate self-consistent barrier heights and because their self-consistent total energies are nearly piecewise linear in fractional electron number─two important points of similarity to the exact functional. We argue that density-driven errors of the energy in a self-consistent density functional calculation are second order in the density error and that large density-driven errors arise primarily from incorrect electron transfers over length scales larger than the diameter of an atom.
Collapse
Affiliation(s)
- Aaron D Kaplan
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Chandra Shahi
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Pradeep Bhetwal
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Raj K Sah
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
| | - John P Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania19122, United States
| |
Collapse
|
13
|
Santra G, Martin JM. Performance of Localized-Orbital Coupled-Cluster Approaches for the Conformational Energies of Longer n-Alkane Chains. J Phys Chem A 2022; 126:9375-9391. [PMID: 36508714 PMCID: PMC9791657 DOI: 10.1021/acs.jpca.2c06407] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We report an update and enhancement of the ACONFL (conformer energies of large alkanes [J. Phys. Chem. A2022,126, 3521-3535]) dataset. For the ACONF12 (n-dodecane) subset, we report basis set limit canonical coupled-cluster with singles, doubles, and perturbative triples [i.e., CCSD(T)] reference data obtained from the MP2-F12/cc-pV{T,Q}Z-F12 extrapolation, [CCSD(F12*)-MP2-F12]/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/aug-cc-pV{D,T}Z calculations. Then, we explored the performance of a variety of single and composite localized-orbital CCSD(T) approximations, ultimately finding an affordable localized natural orbital CCSD(T) [LNO-CCSD(T)]-based post-MP2 correction that agrees to 0.006 kcal/mol mean absolute deviation with the revised canonical reference data. In tandem with canonical MP2-F12 complete basis set extrapolation, this was then used to re-evaluate the ACONF16 and ACONF20 subsets for n-hexadecane and n-icosane, respectively. Combining those with the revised canonical reference data for the dodecane conformers (i.e., ACONF12 subset), a revised ACONFL set was obtained. It was then used to assess the performance of different localized-orbital coupled-cluster approaches, such as pair natural orbital localized CCSD(T) [PNO-LCCSD(T)] as implemented in MOLPRO, DLPNO-CCSD(T0) and DLPNO-CCSD(T1) as implemented in ORCA, and LNO-CCSD(T) as implemented in MRCC, at their respective "Normal", "Tight", "vTight", and "vvTight" accuracy settings. For a given accuracy threshold and basis set, DLPNO-CCSD(T1) and DLPNO-CCSD(T0) perform comparably. With "VeryTightPNO" cutoffs, explicitly correlated DLPNO-CCSD(T1)-F12/VDZ-F12 is the best pick among all the DLPNO-based methods tested. To isolate basis set incompleteness from localized-orbital-related truncation errors (domain, LNOs), we have also compared the localized coupled-cluster approaches with canonical DF-CCSD(T)/aug-cc-pVTZ for the ACONF12 set. We found that gradually tightening the cutoffs improves the performance of LNO-CCSD(T), and using a composite scheme such as vTight + 0.50[vTight - Tight] improves things further. For DLPNO-CCSD(T1), "TightPNO" and "VeryTightPNO" offer a statistically similar accuracy, which gets slightly better when TCutPNO is extrapolated to the complete PNO space limit. Similar to Brauer et al.'s [Phys. Chem. Chem. Phys.2016,18 (31), 20905-20925] previous report for the S66x8 noncovalent interactions, the dispersion-corrected direct random phase approximation (dRPA)-based double hybrids perform remarkably well for the ACONFL set. While the revised reference data do not affect any conclusions on the less accurate methods, they may upend orderings for more accurate methods with error statistics on the same order as the difference between reference datasets.
Collapse
|
14
|
Gao SS, Zhang J, Yin K, Chen YL, He YZ, Li T, Chen LS, Zeng YY, Zheng WX, Ni CL, Liu W, Zhang YM. Syntheses, crystal structures, optical properties, antibacterial activity and DFT calculations of two substituted benzyl triphenylphosphinium picrates. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
15
|
Santra G, Calinsky R, Martin JML. Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights. J Phys Chem A 2022; 126:5492-5505. [PMID: 35930677 PMCID: PMC9393870 DOI: 10.1021/acs.jpca.2c03922] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Revised: 07/28/2022] [Indexed: 11/28/2022]
Abstract
To better understand the thermochemical kinetics and mechanism of a specific chemical reaction, an accurate estimation of barrier heights (forward and reverse) and reaction energies is vital. Because of the large size of reactants and transition state structures involved in real-life mechanistic studies (e.g., enzymatically catalyzed reactions), density functional theory remains the workhorse for such calculations. In this paper, we have assessed the performance of 91 density functionals for modeling the reaction energies and barrier heights on a large and chemically diverse data set (BH9) composed of 449 organic chemistry reactions. We have shown that range-separated hybrid functionals perform better than the global hybrids for BH9 barrier heights and reaction energies. Except for the PBE-based range-separated nonempirical double hybrids, range separation of the exchange term helps improve the performance for barrier heights and reaction energies. The 16-parameter Berkeley double hybrid, ωB97M(2), performs remarkably well for both properties. However, our minimally empirical range-separated double hybrid functionals offer marginally better accuracy than ωB97M(2) for BH9 barrier heights and reaction energies.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| | - Rivka Calinsky
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Molecular Chemistry
and Materials Science, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
16
|
Bryenton KR, Adeleke AA, Dale SG, Johnson ER. Delocalization error: The greatest outstanding challenge in density‐functional theory. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1631] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Kyle R. Bryenton
- Department of Physics and Atmospheric Science Dalhousie University Halifax Nova Scotia Canada
| | | | - Stephen G. Dale
- Queensland Micro‐ and Nanotechnology Centre Griffith University Nathan Queensland Australia
| | - Erin R. Johnson
- Department of Physics and Atmospheric Science Dalhousie University Halifax Nova Scotia Canada
- Department of Chemistry Dalhousie University Halifax Nova Scotia Canada
| |
Collapse
|
17
|
Liang W, Pei Z, Mao Y, Shao Y. Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges. J Chem Phys 2022; 156:210901. [PMID: 35676148 PMCID: PMC9162785 DOI: 10.1063/5.0088271] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2022] [Accepted: 05/15/2022] [Indexed: 01/04/2023] Open
Abstract
Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.
Collapse
Affiliation(s)
- WanZhen Liang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Zheng Pei
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and Department of Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| |
Collapse
|
18
|
Palos E, Lambros E, Swee S, Hu J, Dasgupta S, Paesani F. Assessing the Interplay between Functional-Driven and Density-Driven Errors in DFT Models of Water. J Chem Theory Comput 2022; 18:3410-3426. [PMID: 35506889 DOI: 10.1021/acs.jctc.2c00050] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We investigate the interplay between functional-driven and density-driven errors in different density functional approximations within density functional theory (DFT) and the implications of these errors for simulations of water with DFT-based data-driven potentials. Specifically, we quantify density-driven errors in two widely used dispersion-corrected functionals derived within the generalized gradient approximation (GGA), namely BLYP-D3 and revPBE-D3, and two modern meta-GGA functionals, namely strongly constrained and appropriately normed (SCAN) and B97M-rV. The effects of functional-driven and density-driven errors on the interaction energies are first assessed for the water clusters of the BEGDB dataset. Further insights into the nature of functional-driven errors are gained from applying the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA) to the interaction energies, which demonstrates that functional-driven errors are strongly correlated with the nature of the interactions. We discuss cases where density-corrected DFT (DC-DFT) models display higher accuracy than the original DFT models and cases where reducing the density-driven errors leads to larger deviations from the reference energies due to the presence of large functional-driven errors. Finally, molecular dynamics simulations are performed with data-driven many-body potentials derived from DFT and DC-DFT data to determine the effect that minimizing density-driven errors has on the description of liquid water. Besides rationalizing the performance of widely used DFT models of water, we believe that our findings unveil fundamental relations between the shortcomings of some common DFT approximations and the requirements for accurate descriptions of molecular interactions, which will aid the development of a consistent, DFT-based framework for the development of data-driven and machine-learned potentials for simulations of condensed-phase systems.
Collapse
Affiliation(s)
- Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Steven Swee
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Jie Hu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
19
|
Song S, Vuckovic S, Sim E, Burke K. Density-Corrected DFT Explained: Questions and Answers. J Chem Theory Comput 2022; 18:817-827. [DOI: 10.1021/acs.jctc.1c01045] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722, Korea
| | - Stefan Vuckovic
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce, 73100, Italy
- Department of Chemistry&Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam, 1081HV, The Netherlands
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722, Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
20
|
Mishra P, Yamamoto Y, Johnson JK, Jackson KA, Zope RR, Baruah T. Study of self-interaction-errors in barrier heights using locally scaled and Perdew-Zunger self-interaction methods. J Chem Phys 2022; 156:014306. [PMID: 34998352 DOI: 10.1063/5.0070893] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose, we use the well-known Perdew-Zunger self-interaction-correction (PZSIC) [J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981)] as well as two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al., J. Chem. Phys. 151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation (LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor w based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. The LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughly equal size for the non-hydrogen transfer reactions.
Collapse
Affiliation(s)
- Prakash Mishra
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Yoh Yamamoto
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - J Karl Johnson
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Koblar A Jackson
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Rajendra R Zope
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
21
|
Landeros-Rivera B, Gallegos M, Munarriz J, Laplaza R, Contreras García J. New venues in electron density analysis. Phys Chem Chem Phys 2022; 24:21538-21548. [DOI: 10.1039/d2cp01517j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We provide a comprehensive overview of the chemical information within the electron density: how to extract information, but also how to obtain and how to assess the quality of the...
Collapse
|
22
|
Santra G, Martin JML. Pure and Hybrid SCAN, rSCAN, and r 2SCAN: Which One Is Preferred in KS- and HF-DFT Calculations, and How Does D4 Dispersion Correction Affect This Ranking? Molecules 2021; 27:141. [PMID: 35011372 PMCID: PMC8746565 DOI: 10.3390/molecules27010141] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2021] [Revised: 12/19/2021] [Accepted: 12/24/2021] [Indexed: 11/27/2022] Open
Abstract
Using the large and chemically diverse GMTKN55 dataset, we have tested the performance of pure and hybrid KS-DFT and HF-DFT functionals constructed from three variants of the SCAN meta-GGA exchange-correlation functional: original SCAN, rSCAN, and r2SCAN. Without any dispersion correction involved, HF-SCANn outperforms the two other HF-DFT functionals. In contrast, among the self-consistent variants, SCANn and r2SCANn offer essentially the same performance at lower percentages of HF-exchange, while at higher percentages, SCANn marginally outperforms r2SCANn and rSCANn. However, with D4 dispersion correction included, all three HF-DFT-D4 variants perform similarly, and among the self-consistent counterparts, r2SCANn-D4 outperforms the other two variants across the board. In view of the much milder grid dependence of r2SCAN vs. SCAN, r2SCAN is to be preferred across the board, also in HF-DFT and hybrid KS-DFT contexts.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovot 7610001, Israel;
| | | |
Collapse
|
23
|
Pure and Hybrid SCAN, rSCAN, and r2SCAN: Which One Is Preferred in KS- and HF-DFT Calculations, and How Does D4 Dispersion Correction Affect This Ranking? Molecules 2021. [DOI: 10.3390/molecules27010141 https://www.mdpi.com/1420-3049/27/1/141] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
Using the large and chemically diverse GMTKN55 dataset, we have tested the performance of pure and hybrid KS-DFT and HF-DFT functionals constructed from three variants of the SCAN meta-GGA exchange-correlation functional: original SCAN, rSCAN, and r2SCAN. Without any dispersion correction involved, HF-SCANn outperforms the two other HF-DFT functionals. In contrast, among the self-consistent variants, SCANn and r2SCANn offer essentially the same performance at lower percentages of HF-exchange, while at higher percentages, SCANn marginally outperforms r2SCANn and rSCANn. However, with D4 dispersion correction included, all three HF-DFT-D4 variants perform similarly, and among the self-consistent counterparts, r2SCANn-D4 outperforms the other two variants across the board. In view of the much milder grid dependence of r2SCAN vs. SCAN, r2SCAN is to be preferred across the board, also in HF-DFT and hybrid KS-DFT contexts.
Collapse
|
24
|
Efremenko I, Martin JML. Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls. J Phys Chem A 2021; 125:8987-8999. [PMID: 34586809 DOI: 10.1021/acs.jpca.1c05124] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency. Most tested (meta)-GGA and hybrid density functionals perform better for open-shell than for closed-shell complexes; this is not the case for the double hybrids considered.
Collapse
Affiliation(s)
- Irena Efremenko
- 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
|
25
|
Mosquera MA. Density Functional Calculations Based on the Exponential Ansatz. J Phys Chem A 2021; 125:8751-8763. [PMID: 34582684 DOI: 10.1021/acs.jpca.1c07113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
This work explores the application of the singles-based exponential ansatz to density functional calculations. In contrast to the standard approach where Kohn-Sham (KS) orbitals are determined prior to computing molecular quantities of interest, we consider the single-reference Hartree-Fock wave function as a starting point. Applying the exponential ansatz to this single reference gives an auxiliary wave function that is employed to calculate the electronic properties of the system. This wave function is determined self-consistently through the standard KS Hamiltonian but evaluated over the Hartree-Fock molecular orbital basis. By using spin-symmetry breaking, we recover size-consistent results free of unphysical fractional charges in the dissociation limit. Our method shows consistency with standard KS density functional calculations when the system geometry is similar to the equilibrium one or in repulsive configurations. For moderately long distances between atoms, not at dissociation, because of self-interaction the exponential ansatz may give instabilities in the form of large cluster amplitudes. To avoid these, this work introduces a relatively simple regularization method that preserves size-consistency and penalizes high amplitudes of the cluster operator, whereas the results remain physically meaningful. We also present the time-dependent extension of our theory and show that it can feature quantum states where multiple electrons are excited.
Collapse
Affiliation(s)
- Martín A Mosquera
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, Montana 59718, United States
| |
Collapse
|
26
|
Lambros E, Dasgupta S, Palos E, Swee S, Hu J, Paesani F. General Many-Body Framework for Data-Driven Potentials with Arbitrary Quantum Mechanical Accuracy: Water as a Case Study. J Chem Theory Comput 2021; 17:5635-5650. [PMID: 34370954 DOI: 10.1021/acs.jctc.1c00541] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a general framework for the development of data-driven many-body (MB) potential energy functions (MB-QM PEFs) that represent the interactions between small molecules at an arbitrary quantum-mechanical (QM) level of theory. As a demonstration, a family of MB-QM PEFs for water is rigorously derived from density functionals belonging to different rungs across Jacob's ladder of approximations within density functional theory (MB-DFT) and from Møller-Plesset perturbation theory (MB-MP2). Through a systematic analysis of individual MB contributions to the interaction energies of water clusters, we demonstrate that all MB-QM PEFs preserve the same accuracy as the corresponding ab initio calculations, with the exception of those derived from density functionals within the generalized gradient approximation (GGA). The differences between the DFT and MB-DFT results are traced back to density-driven errors that prevent GGA functionals from accurately representing the underlying molecular interactions for different cluster sizes and hydrogen-bonding arrangements. We show that this shortcoming may be overcome, within the MB formalism, by using density-corrected functionals (DC-DFT) that provide a more consistent representation of each individual MB contribution. This is demonstrated through the development of a MB-DFT PEF derived from DC-PBE-D3 data, which more accurately reproduce the corresponding ab initio results.
Collapse
Affiliation(s)
- Eleftherios Lambros
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Steven Swee
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Jie Hu
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, United States.,Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, United States.,San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
27
|
Ramos C, Muehlbrad J, Janesko BG. Density functionals with full nonlocal exchange, nonlocal rung-3.5 correlation, and D3 dispersion: Combined accuracy for general main-group thermochemistry, kinetics, and noncovalent interactions. J Comput Chem 2021; 42:1974-1981. [PMID: 34387364 DOI: 10.1002/jcc.26728] [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: 06/04/2021] [Accepted: 07/25/2021] [Indexed: 11/11/2022]
Abstract
We introduce the HF-R35-D3(BJ) functional combining full nonlocal exact (Hartree-Fock-like, HF) exchange, inexpensive rung-3.5 correlation constructed from nonlocal one-electron operators, and nonlocal D3 dispersion corrections. HF-R35-D3(BJ) is among the first full-exact-exchange functionals offering competitive accuracy for general main-group thermochemistry, kinetics, and noncovalent interactions. HF-R35-D3(BJ) gives weighted mean absolute deviation WTMAD-2 8.5 kcal/mol across the entire GMTKN55 dataset, outperforming most dispersion-corrected semilocal functionals and approaching the accuracy of dispersion-corrected global hybrids. This requires six fitted parameters, three each in the nonlocal correlation and dispersion corrections. Full nonlocal exchange appears to help give accurate binding energies and reasonable energy orderings for water hexamers. These results motivate continued exploration of inexpensive nonlocal correlation corrections to nonlocal exchange.
Collapse
Affiliation(s)
- Chloe Ramos
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas, USA
| | - Jeremiah Muehlbrad
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas, USA
| | - Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas, USA
| |
Collapse
|
28
|
Sen R, González-Espinoza CE, Zech A, Dreuw A, Wesolowski TA. Benchmark of the Extension of Frozen-Density Embedding Theory to Nonvariational Correlated Methods: The Embedded-MP2 Case. J Chem Theory Comput 2021; 17:4049-4062. [PMID: 34137597 DOI: 10.1021/acs.jctc.1c00228] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The extension of the frozen-density embedding theory for nonvariational methods [J. Chem. Theory Comput. 2020, 16, 6880] was utilized to evaluate intermolecular interaction energies for complexes in the Zhao-Truhlar basis set. In the applied method (FDET-MP2-FAT-LDA), the same auxiliary system is used to evaluate the correlation energy by means of the second-order Møller-Plesset perturbation theory (MP2), as in our previous work [J. Chem. Phys. 2019, 150, 121101]. Local density approximation is used for ExcTnad[ρA,ρB] in both cases. Additionally, the contribution to the energy due to the neglected correlation potential was evaluated and analyzed. The domain of applicability of the local density approximation for ExcTnad[ρA,ρB] was determined based on deviations from the interaction energies from the conventional MP2 calculations. The local density approximation for ExcTnad[ρA,ρB] performs well for hydrogen- or dipole-bound complexes. The relative errors in the interaction energy lie within 3-30%. While for charge-transfer complexes, this approximation fails consistently, and for other types of complexes, the performance of this approximation is not systematic. The sources of error are discussed in detail.
Collapse
Affiliation(s)
- Reena Sen
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Tomasz A Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
29
|
Martín-Fernández C, Harvey JN. On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT. J Phys Chem A 2021; 125:4639-4652. [PMID: 34018759 DOI: 10.1021/acs.jpca.1c01290] [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/27/2022]
Abstract
In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.
Collapse
Affiliation(s)
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan, 200F 3001 Leuven, Belgium
| |
Collapse
|
30
|
Janesko BG. Replacing hybrid density functional theory: motivation and recent advances. Chem Soc Rev 2021; 50:8470-8495. [PMID: 34060549 DOI: 10.1039/d0cs01074j] [Citation(s) in RCA: 59] [Impact Index Per Article: 19.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Density functional theory (DFT) is the most widely-used electronic structure approximation across chemistry, physics, and materials science. Every year, thousands of papers report hybrid DFT simulations of chemical structures, mechanisms, and spectra. Unfortunately, hybrid DFT's accuracy is ultimately limited by tradeoffs between over-delocalization and under-binding. This review summarizes these tradeoffs, and introduces six modern attempts to go beyond them while maintaining hybrid DFT's relatively low computational cost: DFT+U, self-interaction corrections, localized orbital scaling corrections, local hybrid functionals, real-space nondynamical correlation, and our rung-3.5 approach. The review concludes with practical suggestions for DFT users to identify and mitigate these tradeoffs' impact on their simulations.
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, 2800 S. University Dr, Fort Worth, TX 76129, USA.
| |
Collapse
|
31
|
Mariano LA, Vlaisavljevich B, Poloni R. Improved Spin-State Energy Differences of Fe(II) Molecular and Crystalline Complexes via the Hubbard U-Corrected Density. J Chem Theory Comput 2021; 17:2807-2816. [PMID: 33831303 DOI: 10.1021/acs.jctc.1c00034] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We recently showed that the DFT+U approach with a linear-response U yields adiabatic energy differences biased toward high spin [Mariano et al. J. Chem. Theory Comput. 2020, 16, 6755-6762]. Such bias is removed here by employing a density-corrected DFT approach where the PBE functional is evaluated on the Hubbard U-corrected density. The adiabatic energy differences of six Fe(II) molecular complexes computed using this approach, named PBE[U] here, are in excellent agreement with coupled cluster-corrected CASPT2 values for both weak- and strong-field ligands resulting in a mean absolute error (MAE) of 0.44 eV, smaller than that of the recently proposed Hartree-Fock density-corrected DFT (1.22 eV) and any other tested functional, including the best performer TPSSh (0.49 eV). We take advantage of the computational efficiency of this approach and compute the adiabatic energy differences of five molecular crystals using PBE[U] with periodic boundary conditions. The results show, again, an excellent agreement (MAE = 0.07 eV) with experimentally extracted values and a superior performance compared with the best performers M06-L (MAE = 0.08 eV) and TPSSh (MAE = 0.31 eV) computed on molecular fragments.
Collapse
Affiliation(s)
- Lorenzo A Mariano
- University Grenoble Alpes, CNRS, Grenoble-INP, SIMaP, F-38042 Grenoble, France
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Roberta Poloni
- University Grenoble Alpes, CNRS, Grenoble-INP, SIMaP, F-38042 Grenoble, France
| |
Collapse
|
32
|
Mehta N, Fellowes T, White JM, Goerigk L. CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions? J Chem Theory Comput 2021; 17:2783-2806. [PMID: 33881869 DOI: 10.1021/acs.jctc.1c00006] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| | - Thomas Fellowes
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia.,Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Victoria 3010, Australia
| | - Jonathan M White
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia.,Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
33
|
Nam S, Cho E, Sim E, Burke K. Explaining and Fixing DFT Failures for Torsional Barriers. J Phys Chem Lett 2021; 12:2796-2804. [PMID: 33710903 DOI: 10.1021/acs.jpclett.1c00426] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Most torsional barriers are predicted with high accuracies (about 1 kJ/mol) by standard semilocal functionals, but a small subset was found to have much larger errors. We created a database of almost 300 carbon-carbon torsional barriers, including 12 poorly behaved barriers, that stem from the Y═C-X group, where Y is O or S and X is a halide. Functionals with enhanced exchange mixing (about 50%) worked well for all barriers. We found that poor actors have delocalization errors caused by hyperconjugation. These problematic calculations are density-sensitive (i.e., DFT predictions change noticeably with the density), and using HF densities (HF-DFT) fixes these issues. For example, conventional B3LYP performs as accurately as exchange-enhanced functionals if the HF density is used. For long-chain conjugated molecules, HF-DFT can be much better than exchange-enhanced functionals. We suggest that HF-PBE0 has the best overall performance.
Collapse
Affiliation(s)
- Seungsoo Nam
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, South Korea
| | - Eunbyol Cho
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, South Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, South Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
34
|
Wagle K, Santra B, Bhattarai P, Shahi C, Pederson MR, Jackson KA, Perdew JP. Self-interaction correction in water–ion clusters. J Chem Phys 2021; 154:094302. [DOI: 10.1063/5.0041620] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kamal Wagle
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Puskar Bhattarai
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Chandra Shahi
- Department of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Mark R. Pederson
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Koblar A. Jackson
- Department of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - John P. Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
35
|
Santra G, Martin JML. What Types of Chemical Problems Benefit from Density-Corrected DFT? A Probe Using an Extensive and Chemically Diverse Test Suite. J Chem Theory Comput 2021; 17:1368-1379. [PMID: 33625863 PMCID: PMC8028055 DOI: 10.1021/acs.jctc.0c01055] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
![]()
For the large and
chemically diverse GMTKN55 benchmark suite, we
have studied the performance of density-corrected density functional
theory (HF-DFT), compared to self-consistent DFT, for several pure
and hybrid GGA and meta-GGA exchange–correlation (XC) functionals
(PBE, BLYP, TPSS, and SCAN) as a function of the percentage of HF
exchange in the hybrid. The D4 empirical dispersion correction has
been added throughout. For subsets dominated by dynamical correlation,
HF-DFT is highly beneficial, particularly at low HF exchange percentages.
This is especially true for noncovalent interactions where the electrostatic
component is dominant, such as hydrogen and halogen bonds: for π-stacking,
HF-DFT is detrimental. For subsets with significant nondynamical correlation
(i.e., where a Hartree–Fock determinant is not a good zero-order
wavefunction), HF-DFT may do more harm than good. While the self-consistent
series show optima at or near 37.5% (i.e., 3/8) for all four XC functionals—consistent
with Grimme’s proposal of the PBE38 functional—HF-BnLYP-D4, HF-PBEn-D4, and HF-TPSSn-D4 all exhibit minima nearer 25% (i.e., 1/4) as the use
of HF orbitals greatly mitigates the error at 25% for barrier heights.
Intriguingly, for HF-SCANn-D4, the minimum is near
10%, but the weighted mean absolute error (WTMAD2) for GMTKN55 is
only barely lower than that for HF-SCAN-D4 (i.e., where the post-HF
step is a pure meta-GGA). The latter becomes an attractive option,
only slightly more costly than pure Hartree–Fock, and devoid
of adjustable parameters other than the three in the dispersion correction.
Moreover, its WTMAD2 is only surpassed by the highly empirical M06-2X
and by the combinatorially optimized empirical range-separated hybrids
ωB97X-V and ωB97M-V.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
36
|
Mei Y, Yang N, Yang W. Describing polymer polarizability with localized orbital scaling correction in density functional theory. J Chem Phys 2021; 154:054302. [PMID: 33557560 DOI: 10.1063/5.0035883] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Polarizability reflects the response of the molecular charge distribution to an applied external electric field and thus closely relates to the molecular electron density. For the calculation of polarizability within density functional theory (DFT), it is well known that conventional density functional approximations (DFAs) greatly overestimate the results for polymers with long chains and the π-conjugated system. This is a manifestation of the delocalization error of the commonly used DFAs-they normally produce too delocalized electron density and underestimate the total energy for systems with fractional charge character, which occurs for long molecules in a longitudinal electric field. Thus, to achieve an accurate description of polarizabilities for polymeric molecular systems from DFT, applying DFAs with minimal delocalization error is very important. In this work, we use the recently developed localized orbital scaling correction (LOSC) to the conventional DFAs, which has been shown to largely eliminate the delocalization error, to calculate and study the polarizabilities of three classic polymers, polyyne, polyacetylene, and hydrogen chain. The results from this work demonstrate that applying LOSC to conventional DFAs with self-consistent field calculations can largely improve the description of polarizability from DFT calculations and the improved quality of electron density in LOSC leads to the improved results of polarizability of the polymers. However, the improvement is not complete and adjustment of the parameters in the LOSC method can further improve the accuracy to reach the level similar to the MP2 method. This work also points to the direction for the further development of LOSC in self-consistent calculations.
Collapse
Affiliation(s)
- Yuncai Mei
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Nathan Yang
- North Carolina School of Science and Mathematics, 1219 Broad St., Durham, North Carolina 27705, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
37
|
Abstract
Empirical fitting of parameters in approximate density functionals is common. Such fits conflate errors in the self-consistent density with errors in the energy functional, but density-corrected DFT (DC-DFT) separates these two. We illustrate with catastrophic failures of a toy functional applied to H2+ at varying bond lengths, where the standard fitting procedure misses the exact functional; Grimme's D3 fit to noncovalent interactions, which can be contaminated by large density errors such as in the WATER27 and B30 data sets; and double-hybrids trained on self-consistent densities, which can perform poorly on systems with density-driven errors. In these cases, more accurate results are found at no additional cost by using Hartree-Fock (HF) densities instead of self-consistent densities. For binding energies of small water clusters, errors are greatly reduced. Range-separated hybrids with 100% HF at large distances suffer much less from this effect.
Collapse
Affiliation(s)
- Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Stefan Vuckovic
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
38
|
Jin Y, Su NQ, Chen Z, Yang W. Introductory lecture: when the density of the noninteracting reference system is not the density of the physical system in density functional theory. Faraday Discuss 2020; 224:9-26. [PMID: 33084699 PMCID: PMC7746600 DOI: 10.1039/d0fd00102c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A major challenge in density functional theory (DFT) is the development of density functional approximations (DFAs) to overcome errors in existing DFAs, leading to more complex functionals. For such functionals, we consider roles of the noninteracting reference systems. The electron density of the Kohn-Sham (KS) reference with a local potential has been traditionally defined as being equal to the electron density of the physical system. This key idea has been applied in two ways: the inverse calculation of such a local KS potential for the reference from a given density and the direct calculation of density and energy based on given DFAs. By construction, the inverse calculation can yield a KS reference with the density equal to the input density of the physical system. In application of DFT, however, it is the direct calculation of density and energy from a DFA that plays a central role. For direct calculations, we find that the self-consistent density of the KS reference defined by the optimized effective potential (OEP), is not the density of the physical system, when the DFA is dependent on the external potential. This inequality holds also for the density of generalized KS (GKS) or generalized OEP reference, which allows a nonlocal potential, when the DFA is dependent on the external potential. Instead, the density of the physical system, consistent with a given DFA, is given by the linear response of the total energy with respect to the variation of the external potential. This is a paradigm shift in DFT on the use of noninteracting references: the noninteracting KS or GKS references represent the explicit computational variables for energy minimization, but not the density of the physical system for external potential-dependent DFAs. We develop the expressions for the electron density so defined through the linear response for general DFAs, demonstrate the results for orbital functionals and for many-body perturbation theory within the second-order and the random-phase approximation, and explore the connections to developments in DFT.
Collapse
Affiliation(s)
- Ye Jin
- Department of Chemistry, Duke University, Durham, NC 27708, USA.
| | - Neil Qiang Su
- Department of Chemistry, Duke University, Durham, NC 27708, USA.
| | - Zehua Chen
- Department of Chemistry, Duke University, Durham, NC 27708, USA.
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, NC 27708, USA. and Department of Physics, Duke University, Durham, NC 27708, USA
| |
Collapse
|
39
|
Callow TJ, Pearce BJ, Pitts T, Lathiotakis NN, Hodgson MJP, Gidopoulos NI. Improving the exchange and correlation potential in density-functional approximations through constraints. Faraday Discuss 2020; 224:126-144. [PMID: 32940317 DOI: 10.1039/d0fd00069h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review and expand on our work to impose constraints on the effective Kohn-Sham (KS) potential of local and semi-local density-functional approximations. Constraining the minimisation of the approximate total energy density-functional invariably leads to an optimised effective potential (OEP) equation, the solution of which yields the KS potential. We review briefly our previous work on this and demonstrate with numerous examples that despite the well-known mathematical issues of the OEP with finite basis sets, our OEP equations are numerically robust. We demonstrate that appropriately constraining the 'screening charge' which corresponds to the Hartree, exchange and correlation potential not only corrects its asymptotic behaviour but also allows the exchange and correlation potential to exhibit a non-zero derivative discontinuity, a feature of the exact KS potential that is necessary for the accurate prediction of band-gaps in solids but very hard to capture with semi-local approximations.
Collapse
Affiliation(s)
- Timothy J Callow
- Department of Physics, Durham University, South Road, Durham, DH1 3LE, UK.
| | | | | | | | | | | |
Collapse
|
40
|
Bogojeski M, Vogt-Maranto L, Tuckerman ME, Müller KR, Burke K. Quantum chemical accuracy from density functional approximations via machine learning. Nat Commun 2020; 11:5223. [PMID: 33067479 PMCID: PMC7567867 DOI: 10.1038/s41467-020-19093-1] [Citation(s) in RCA: 129] [Impact Index Per Article: 32.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Accepted: 09/24/2020] [Indexed: 12/21/2022] Open
Abstract
Kohn-Sham density functional theory (DFT) is a standard tool in most branches of chemistry, but accuracies for many molecules are limited to 2-3 kcal ⋅ mol-1 with presently-available functionals. Ab initio methods, such as coupled-cluster, routinely produce much higher accuracy, but computational costs limit their application to small molecules. In this paper, we leverage machine learning to calculate coupled-cluster energies from DFT densities, reaching quantum chemical accuracy (errors below 1 kcal ⋅ mol-1) on test data. Moreover, density-based Δ-learning (learning only the correction to a standard DFT calculation, termed Δ-DFT ) significantly reduces the amount of training data required, particularly when molecular symmetries are included. The robustness of Δ-DFT is highlighted by correcting "on the fly" DFT-based molecular dynamics (MD) simulations of resorcinol (C6H4(OH)2) to obtain MD trajectories with coupled-cluster accuracy. We conclude, therefore, that Δ-DFT facilitates running gas-phase MD simulations with quantum chemical accuracy, even for strained geometries and conformer changes where standard DFT fails.
Collapse
Affiliation(s)
- Mihail Bogojeski
- Machine Learning Group, Technische Universität Berlin, Marchstr. 23, 10587, Berlin, Germany
| | | | - Mark E Tuckerman
- Department of Chemistry, New York University, New York, NY, 10003, USA.
- Courant Institute of Mathematical Science, New York University, New York, NY, 10012, USA.
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Road North, Shanghai, 200062, China.
| | - Klaus-Robert Müller
- Machine Learning Group, Technische Universität Berlin, Marchstr. 23, 10587, Berlin, Germany.
- Department of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea.
- Max-Planck-Institut für Informatik, Stuhlsatzenhausweg, 66123, Saarbrücken, Germany.
| | - Kieron Burke
- Department of Physics and Astronomy, University of California, Irvine, CA, 92697, USA.
- Department of Chemistry, University of California, Irvine, CA, 92697, USA.
| |
Collapse
|
41
|
Nam S, Song S, Sim E, Burke K. Measuring Density-Driven Errors Using Kohn-Sham Inversion. J Chem Theory Comput 2020; 16:5014-5023. [PMID: 32667787 DOI: 10.1021/acs.jctc.0c00391] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Kohn-Sham (KS) inversion, that is, the finding of the exact KS potential for a given density, is difficult in localized basis sets. We study the precision and reliability of several inversion schemes, finding estimates of density-driven errors at a useful level of accuracy. In typical cases of substantial density-driven errors, Hartree-Fock density functional theory (HF-DFT) is almost as accurate as DFT evaluated on CCSD(T) densities. A simple approximation in practical HF-DFT also makes errors much smaller than the density-driven errors being calculated. Two paradigm examples, stretched NaCl and the HO·Cl- radical, illustrate just how accurate HF-DFT is.
Collapse
Affiliation(s)
- Seungsoo Nam
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
42
|
Moltved KA, Kepp KP. Using electronegativity and hardness to test density functionals. J Chem Phys 2020; 152:244113. [PMID: 32610960 DOI: 10.1063/5.0006189] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Density functional theory (DFT) is used in thousands of papers each year, yet lack of universality reduces DFT's predictive capacity, and functionals may produce energy-density imbalances. The absolute electronegativity (χ) and hardness (η) directly reflect the energy-density relationship via the chemical potential ∂E/∂N and we thus hypothesized that they probe universality. We studied χ and η for atoms Z = 1-36 using 50 diverse functionals covering all major classes. Very few functionals describe both χ and η well. η benefits from error cancellation, whereas χ is marred by error propagation from IP and EA; thus, almost all standard GGA and hybrid functionals display a plateau in the MAE at ∼0.2 eV-0.3 eV for η. In contrast, variable performance for χ indicates problems in describing the chemical potential by DFT. The accuracy and precision of a functional is far from linearly related, yet for a universal functional, we expect linearity. Popular functionals such as B3LYP, PBE, and revPBE perform poorly for both properties. Density sensitivity calculations indicate large density-derived errors as occupation of degenerate p- and d-orbitals causes "non-universality" and large dependency on exact exchange. Thus, we argue that performance for χ for the same systems is a hallmark of an important aspect of universality by probing ∂E/∂N. With this metric, B98, B97-1, PW6B95D3, MN-15, rev-TPSS, HSE06, and APFD are the most "universal" among the tested functionals. B98 and B97-1 are accurate for very diverse metal-ligand bonds, supporting that a balanced description of ∂E/∂N and ∂E2/∂N2, via χ and η, is probably a first simple probe of universality.
Collapse
Affiliation(s)
- Klaus A Moltved
- Technical University of Denmark, DTU Chemistry, Building 206, 2800 Kgs. Lyngby, Denmark
| | - Kasper P Kepp
- Technical University of Denmark, DTU Chemistry, Building 206, 2800 Kgs. Lyngby, Denmark
| |
Collapse
|
43
|
Skelton JM, Gunn DSD, Metz S, Parker SC. Accuracy of Hybrid Functionals with Non-Self-Consistent Kohn-Sham Orbitals for Predicting the Properties of Semiconductors. J Chem Theory Comput 2020; 16:3543-3557. [PMID: 32369352 DOI: 10.1021/acs.jctc.9b01218] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurately modeling the electronic structure of materials is a persistent challenge to high-throughput screening. A promising means of balancing accuracy against computational cost is non-self-consistent calculations with hybrid density-functional theory, where the electronic band energies are evaluated using a hybrid functional from orbitals obtained with a less demanding (semi)local functional. We have quantified the performance of this technique for predicting the physical properties of 16 tetrahedral semiconductors with bandgaps from 0.2 to 5.5 eV. Provided the base functional predicts a nonmetallic electronic structure, bandgaps within 5% of the PBE0 and HSE06 gaps can be obtained with an order of magnitude reduction in computing time. The positions of the valence and conduction band extrema and the Fermi level are well reproduced, enabling calculation of the band dispersion, density of states, and dielectric properties using Fermi's Golden Rule. While the error in the non-self-consistent total energies is ∼50 meV atom-1, the energy-volume curves are reproduced accurately enough to obtain the equilibrium volume and bulk modulus with minimal error. We also test the dielectric-dependent scPBE0 functional and obtain bandgaps and dielectric constants to within 2.5% of the self-consistent results, which amounts to a significant improvement over self-consistent PBE0 for a similar computational cost. We identify cases where the non-self-consistent approach is expected to perform poorly and demonstrate that partial self-consistency provides a practical and efficient workaround. Finally, we perform proof-of-concept calculations on CoO and NiO to demonstrate the applicability of the technique to strongly correlated open-shell transition-metal oxides.
Collapse
Affiliation(s)
- Jonathan M Skelton
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom.,Department of Chemistry, University of Bath, Bath BA2 7AY, United Kingdom
| | - David S D Gunn
- STFC Daresbury Laboratory, Keckwick Ln., Daresbury, Warrington WA4 4AD, United Kingdom
| | - Sebastian Metz
- STFC Daresbury Laboratory, Keckwick Ln., Daresbury, Warrington WA4 4AD, United Kingdom.,Fraunhofer Institute for Solar Energy Systems ISE, Heidenhofstraße 2, 79110 Freiburg, Germany
| | - Stephen C Parker
- Department of Chemistry, University of Bath, Bath BA2 7AY, United Kingdom
| |
Collapse
|
44
|
Abstract
This work examines the suitability of meta-GGA functionals for symmetry-adapted perturbation theory (SAPT) calculations. The assessment is based on the term-by-term comparison with the benchmark SAPT variant based on coupled-cluster singles and doubles description of monomers, SAPT(CCSD). Testing systems include molecular complexes ranging from strong to weak and the He dimer. The following nonempirical meta-GGAs are examined: TPSS, revTPSS, MVS, SCAN, and SCAN0 with and without the asymptotic correction (AC) of the exchange-correlation potential. One range-separated meta-GGA functional, LC-PBETPSS, is also included. The AC-corrected pure meta-GGAs (with the exception of MVS) represent a definite progress in SAPT(DFT) compared to pure GGA, such as PBEAC, with their more consistent predictions of energy components. However, none of the meta-GGAs is better than the hybrid GGA approach SAPT(PBE0AC). The SAPT(DFT) electrostatic energy offers the most sensitive probe of the quality of the underlying DFT density. Both SCAN- and TPSS-based electrostatic energies agree with reference to within 5% or better which is an excellent result. We find that SCAN0 can be used in SAPT without the AC correction. The long-range corrected LC-PBETPSS is a reliable performer both for the components and total interaction energies.
Collapse
|
45
|
Abstract
By invoking a divide-and-conquer strategy, subsystem DFT dramatically reduces the computational cost of large-scale, ab initio electronic structure simulations of molecules and materials. The central ingredient setting subsystem DFT apart from Kohn-Sham DFT is the nonadditive kinetic energy functional (NAKE). Currently employed NAKEs are at most semilocal (i.e., they only depend on the electron density and its gradient), and as a result of this approximation, so far large-scale simulations only included systems composed of weakly interacting subsystems. In this work, we advance the state-of-the-art by introducing fully nonlocal NAKEs in subsystem DFT simulations for the first time. A benchmark analysis based on the S22-5 test set shows that nonlocal NAKEs considerably improve the computed interaction energies and electron densities compared to commonly employed GGA NAKEs, especially when increasing intersubsystem electron density overlap is considered. Most importantly, we resolve the long-standing problem of too attractive interaction energy curves typically resulting from the use of GGA NAKEs.
Collapse
Affiliation(s)
- Wenhui Mi
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
- Department of Physics , Rutgers University , Newark , New Jersey 07102 , United States
| | - Michele Pavanello
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
- Department of Physics , Rutgers University , Newark , New Jersey 07102 , United States
| |
Collapse
|
46
|
Otero-de-la-Roza A, LeBlanc LM, Johnson ER. What is “many-body” dispersion and should I worry about it? Phys Chem Chem Phys 2020; 22:8266-8276. [DOI: 10.1039/d0cp01213k] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
“Many-body” dispersion can refer to two distinct phenomena, here termed electronic and atomic many-body effects, both of which cause the dispersion energy to be non-additive.
Collapse
Affiliation(s)
- A. Otero-de-la-Roza
- Departamento de Química Física y Analítica and MALTA-Consolider Team
- Facultad de Química
- Universidad de Oviedo
- 33006 Oviedo
- Spain
| | | | | |
Collapse
|
47
|
Mi W, Ramos P, Maranhao J, Pavanello M. Ab Initio Structure and Dynamics of CO 2 at Supercritical Conditions. J Phys Chem Lett 2019; 10:7554-7559. [PMID: 31738567 DOI: 10.1021/acs.jpclett.9b03054] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Green technologies rely on green solvents and fluids. Among them, supercritical CO2 already finds many important applications. The molecular-level understanding of the dynamics and structure of this supercritical fluid is a prerequisite for rational design of future green technologies. Unfortunately, the commonly employed Kohn-Sham density functional theory (DFT) is too computationally demanding to produce meaningfully converged dynamics within a reasonable time and with a reasonable computational effort. Thanks to subsystem DFT, we analyze finite-size effects by considering simulation cells of varying sizes (up to 256 independent molecules in the cell) and finite-time effects by running 100 ps trajectories. We find that the simulations are in reasonable and semiquantitative agreement with the available neutron diffraction experiments and that, as opposed to the gas phase, the CO2 molecules in the fluid are bent with an average OCO angle of 175.8°. Our simulations also confirm that the dimer T-shape is the most prevalent configuration. Our results further strengthen the experiment-simulation agreement for this fluid when comparing radial distribution functions and diffusion coefficient, confirming subsystem DFT as a viable tool for modeling structure and dynamics of condensed-phase systems.
Collapse
Affiliation(s)
- Wenhui Mi
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Pablo Ramos
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Jack Maranhao
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
| | - Michele Pavanello
- Department of Chemistry , Rutgers University , Newark , New Jersey 07102 , United States
- Department of Physics , Rutgers University , Newark , New Jersey 07102 , United States
| |
Collapse
|
48
|
Abstract
The attraction density functional theory (DFT) has for electronic structure theory is that it is easier to do computationally than ab initio, correlated wavefunction methods, due to its effective one-particle structure. On the contrary, ab initio theorists insist on the ability to converge to the right answer in appropriate limits, but this requires a treatment of the reduced two-particle density matrix. DFT avoids that by appealing to an "existence" theorem (not a constructive one) that all its effects are subsummed into a DFT functional of the one-particle density. However, the existence of thousands of DFT functionals emphasizes that there is no satisfactory way to systematically improve the Kohn-Sham (KS) version as most changes in parameterization or formulation seldom lead to a new functional that is genuinely better than others. Some researchers in the DFT community try to address this issue by imposing conditions rigorously derived from exact DFT considerations, but to date, no one has shown how this route will ever lead to converged results even for the ground state, much less for all the other electronic states obtained from time-dependent DFT that are critically important for chemistry. On the contrary, coupled-cluster (CC) theory and its equation-of-motion extensions provide rigorous results for both that KS-DFT methods are attempting to emulate. How to use them and their exact formal properties to tie CC theory to an effective one-particle form is the target of this perspective. This route addresses the devil's triangle of KS-DFT problems: the one-particle spectrum, self-interaction, and the integer discontinuity.
Collapse
Affiliation(s)
- Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
49
|
Vuckovic S, Song S, Kozlowski J, Sim E, Burke K. Density Functional Analysis: The Theory of Density-Corrected DFT. J Chem Theory Comput 2019; 15:6636-6646. [DOI: 10.1021/acs.jctc.9b00826] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Stefan Vuckovic
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| | - Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - John Kozlowski
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| | - Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
50
|
Yamamoto Y, Diaz CM, Basurto L, Jackson KA, Baruah T, Zope RR. Fermi-Löwdin orbital self-interaction correction using the strongly constrained and appropriately normed meta-GGA functional. J Chem Phys 2019; 151:154105. [PMID: 31640373 DOI: 10.1063/1.5120532] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Despite the success of density functional approximations (DFAs) in describing the electronic properties of many-electron systems, the most widely used approximations suffer from self-interaction errors (SIEs) that limit their predictive power. Here, we describe the effects of removing SIE from the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation using the Fermi-Löwdin Orbital Self-Interaction Correction (FLOSIC) method. FLOSIC is a size-extensive implementation of the Perdew-Zunger self-interaction correction (PZ-SIC) formalism. We find that FLOSIC-SCAN calculations require careful treatment of numerical details and an integration grid that yields reliable accuracy with this approach. We investigate the performance of FLOSIC-SCAN for predicting a wide array of properties and find that it provides better results than FLOSIC-LDA and FLOSIC-PBE in nearly all cases. It also gives better predictions than SCAN for orbital energies and dissociation energies where self-interaction effects are known to be important, but total energies and atomization energies are made worse. For these properties, we also investigate the use of the self-consistent FLOSIC-SCAN density in the SCAN functional and find that this DFA@FLOSIC-DFA approach yields improved results compared to pure, self-consistent SCAN calculations. Thus, FLOSIC-SCAN provides improved results over the parent SCAN functional in cases where SIEs are dominant, and even when they are not, if the SCAN@FLOSIC-SCAN method is used.
Collapse
Affiliation(s)
- Yoh Yamamoto
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Carlos M Diaz
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Luis Basurto
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Koblar A Jackson
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Tunna Baruah
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|