1
|
Burgess AC, Linscott E, O'Regan DD. Tilted-Plane Structure of the Energy of Finite Quantum Systems. PHYSICAL REVIEW LETTERS 2024; 133:026404. [PMID: 39073931 DOI: 10.1103/physrevlett.133.026404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/29/2023] [Revised: 04/18/2024] [Accepted: 05/17/2024] [Indexed: 07/31/2024]
Abstract
The piecewise linearity condition on the total energy with respect to the total magnetization of finite quantum systems is derived using the infinite-separation-limit technique. This generalizes the well-known constancy condition, related to static correlation error, in approximate density functional theory. The magnetic analog of Koopmans' theorem in density functional theory is also derived. Moving to fractional electron count, the tilted-plane condition is derived, lifting certain assumptions in previous works. This generalization of the flat-plane condition characterizes the total energy surface of a finite system for all values of electron count N and magnetization M. This result is used in combination with tabulated spectroscopic data to show the flat-plane structure of the oxygen atom, among others. We find that derivative discontinuities with respect to electron count sometimes occur at noninteger values. A diverse set of tilted-plane structures is shown to occur in d-orbital subspaces, depending on chemical coordination. General occupancy-based total-energy expressions are demonstrated thereby to be necessarily dependent on the symmetry-imposed degeneracies.
Collapse
|
2
|
Janesko BG. Multiconfigurational Correlation at DFT + U Cost: On-Site Electron-Electron Interactions Yield a Block-Localized Configuration Interaction Hamiltonian. J Phys Chem A 2024; 128:5077-5087. [PMID: 38878060 DOI: 10.1021/acs.jpca.4c02326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/28/2024]
Abstract
This work presents a first-principles wavefunction-in-DFT approach based on the Hubbard density functional theory (DFT) + U method. This approach begins with the standard DFT reference system of noninteracting electrons and introduces an electron-electron interaction projected onto DFT+U-type atomic states. The reference system's configuration interaction Hamiltonian is block-localized to these states and can be expressed in terms of state occupation numbers, state self-energies (which correspond to unscreened Hubbard U values), and the promotion energies of doubly excited Slater determinants. Simple approximations for the promotion energies provide multiconfigurational correlation energies without requiring explicit orbital localization/transform. Numerical results for fractionally occupied chromium atom, bonded chromium dimer, dissociating covalent bonds, and large active spaces show that the approach provides beyond-zero-sum accuracy at computational cost comparable to standard DFT+U.
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas 76129, United States
| |
Collapse
|
3
|
Macke E, Timrov I, Marzari N, Ciacchi LC. Orbital-Resolved DFT +U for Molecules and Solids. J Chem Theory Comput 2024; 20:4824-4843. [PMID: 38820347 PMCID: PMC11171274 DOI: 10.1021/acs.jctc.3c01403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Revised: 03/25/2024] [Accepted: 05/02/2024] [Indexed: 06/02/2024]
Abstract
We present an orbital-resolved extension of the Hubbard U correction to density-functional theory (DFT). Compared to the conventional shell-averaged approach, the prediction of energetic, electronic and structural properties is strongly improved, particularly for compounds characterized by both localized and hybridized states in the Hubbard manifold. The numerical values of all Hubbard parameters are readily obtained from linear-response calculations. The relevance of this more refined approach is showcased by its application to bulk solids pyrite (FeS2) and pyrolusite (β-MnO2), as well as to six Fe(II) molecular complexes. Our findings indicate that a careful definition of Hubbard manifolds is indispensable for extending the applicability of DFT+U beyond its current boundaries. The present orbital-resolved scheme aims to provide a computationally undemanding yet accurate tool for electronic structure calculations of charge-transfer insulators, transition-metal (TM) complexes and other compounds displaying significant orbital hybridization.
Collapse
Affiliation(s)
- Eric Macke
- Faculty
of Production Engineering, Bremen Center
for Computational Materials Science and MAPEX Center for Materials
and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, 28359 Bremen, Germany
| | - Iurii Timrov
- Theory
and Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Nicola Marzari
- Theory
and Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique
Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- University
of Bremen Excellence Chair, Bremen Center
for Computational Materials Science, Am Fallturm 1, 28359 Bremen, Germany
| | - Lucio Colombi Ciacchi
- Faculty
of Production Engineering, Bremen Center
for Computational Materials Science and MAPEX Center for Materials
and Processes, Hybrid Materials Interfaces Group, University of Bremen, Am Fallturm 1, 28359 Bremen, Germany
| |
Collapse
|
4
|
Daas K, Klute E, Seidl M, Gori-Giorgi P. Møller-Plesset Adiabatic Connection at Large Coupling Strengths for Open-Shell Systems. J Phys Chem A 2024; 128:4138-4149. [PMID: 38717868 PMCID: PMC11129316 DOI: 10.1021/acs.jpca.4c00788] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2024] [Revised: 03/29/2024] [Accepted: 04/05/2024] [Indexed: 05/24/2024]
Abstract
We study the adiabatic connection that has as weak-coupling expansion the Møller-Plesset perturbation series, generalizing to the open-shell case previous closed-shell results for the large-coupling limit. We first focus on the hydrogen atom with fractional spins, providing results along the adiabatic connection from small to large coupling strengths. We reveal an intriguing phase diagram and an equation for the large-coupling leading order that has closed-form solutions for specific choices of its relevant quantum numbers. We then show that the hydrogen atom results provide variational estimates for the large-coupling leading terms for the general many-electron open-shell case in terms of functionals of the Hartree-Fock α-spin and β-spin densities.
Collapse
Affiliation(s)
- Kimberly
J. Daas
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
| | - Eveline Klute
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
| | - Michael Seidl
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
| | - Paola Gori-Giorgi
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, Amsterdam 1081 HV, The Netherlands
- Microsoft
Research AI for Science, Evert van de Beekstraat 354, Schiphol 1118 CZ, The Netherlands
| |
Collapse
|
5
|
Cho Y, Kulik HJ. Improving gas adsorption modeling for MOFs by local calibration of Hubbard U parameters. J Chem Phys 2024; 160:154101. [PMID: 38624114 DOI: 10.1063/5.0201934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Accepted: 03/31/2024] [Indexed: 04/17/2024] Open
Abstract
While computational screening with density functional theory (DFT) is frequently employed for the screening of metal-organic frameworks (MOFs) for gas separation and storage, commonly applied generalized gradient approximations (GGAs) exhibit self-interaction errors, which hinder the predictions of adsorption energies. We investigate the Hubbard U parameter to augment DFT calculations for full periodic MOFs, targeting a more precise modeling of gas molecule-MOF interactions, specifically for N2, CO2, and O2. We introduce a calibration scheme for the U parameter, which is tailored for each MOF, by leveraging higher-level calculations on the secondary building unit (SBU) of the MOF. When applied to the full periodic MOF, the U parameter calibrated against hybrid HSE06 calculations of SBUs successfully reproduces hybrid-quality calculations of the adsorption energy of the periodic MOF. The mean absolute deviation of adsorption energies reduces from 0.13 eV for a standard GGA treatment to 0.06 eV with the calibrated U, demonstrating the utility of the calibration procedure when applied to the full MOF structure. Furthermore, attempting to use coupled cluster singles and doubles with perturbative triples calculations of isolated SBUs for this calibration procedure shows varying degrees of success in predicting the experimental heat of adsorption. It improves accuracy for N2 adsorption for cases of overbinding, whereas its impact on CO2 is minimal, and ambiguities in spin state assignment hinder consistent improvements of O2 adsorption. Our findings emphasize the limitations of cluster models and advocate the use of full periodic MOF systems with a calibrated U parameter, providing a more comprehensive understanding of gas adsorption in MOFs.
Collapse
Affiliation(s)
- Yeongsu Cho
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
6
|
Goshen Y, Kraisler E. Ensemble Ground State of a Many-Electron System with Fractional Electron Number and Spin: Piecewise-Linearity and Flat-Plane Condition Generalized. J Phys Chem Lett 2024; 15:2337-2343. [PMID: 38386920 PMCID: PMC10926161 DOI: 10.1021/acs.jpclett.3c03509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Revised: 02/05/2024] [Accepted: 02/09/2024] [Indexed: 02/24/2024]
Abstract
Description of many-electron systems with a fractional electron number (Ntot) and fractional spin (Mtot) is of great importance in physical chemistry, solid-state physics, and materials science. In this Letter, we provide an exact description of the zero-temperature ensemble ground state of a general, finite, many-electron system and characterize the dependence of the energy and the spin-densities on both Ntot and Mtot, when the total spin is at its equilibrium value. We generalize the piecewise-linearity principle and the flat-plane condition and determine which pure states contribute to the ground-state ensemble. We find a new derivative discontinuity, which manifests for spin variation at a constant Ntot, as a jump in the Kohn-Sham potential. We identify a previously unknown degeneracy of the ground state, such that the total energy and density are unique, but the spin-densities are not. Our findings serve as a basis for development of advanced approximations in density functional theory and other many-electron methods.
Collapse
Affiliation(s)
- Yuli Goshen
- Fritz Haber Research Center for Molecular
Dynamics and Institute of Chemistry, The
Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| | - Eli Kraisler
- Fritz Haber Research Center for Molecular
Dynamics and Institute of Chemistry, The
Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| |
Collapse
|
7
|
Jana S, Herbert JM. Fractional-Electron and Transition-Potential Methods for Core-to-Valence Excitation Energies Using Density Functional Theory. J Chem Theory Comput 2023; 19:4100-4113. [PMID: 37312236 DOI: 10.1021/acs.jctc.3c00202] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Methods for computing X-ray absorption spectra based on a constrained core hole (possibly containing a fractional electron) are examined. These methods are based on Slater's transition concept and its generalizations, wherein core-to-valence excitation energies are determined using Kohn-Sham orbital energies. Methods examined here avoid promoting electrons beyond the lowest unoccupied molecular orbital, facilitating robust convergence. Variants of these ideas are systematically tested, revealing a best-case accuracy of 0.3-0.4 eV (with respect to experiment) for K-edge transition energies. Absolute errors are much larger for higher-lying near-edge transitions but can be reduced below 1 eV by introducing an empirical shift based on a charge-neutral transition-potential method, in conjunction with functionals such as SCAN, SCAN0, or B3LYP. This procedure affords an entire excitation spectrum from a single fractional-electron calculation, at the cost of ground-state density functional theory and without the need for state-by-state calculations. This shifted transition-potential approach may be especially useful for simulating transient spectroscopies or in complex systems where excited-state Kohn-Sham calculations are challenging.
Collapse
Affiliation(s)
- Subrata Jana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
8
|
Jana S, Herbert JM. Slater transition methods for core-level electron binding energies. J Chem Phys 2023; 158:094111. [PMID: 36889976 DOI: 10.1063/5.0134459] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023] Open
Abstract
Methods for computing core-level ionization energies using self-consistent field (SCF) calculations are evaluated and benchmarked. These include a "full core hole" (or "ΔSCF") approach that fully accounts for orbital relaxation upon ionization, but also methods based on Slater's transition concept in which the binding energy is estimated from an orbital energy level that is obtained from a fractional-occupancy SCF calculation. A generalization that uses two different fractional-occupancy SCF calculations is also considered. The best of the Slater-type methods afford mean errors of 0.3-0.4 eV with respect to experiment for a dataset of K-shell ionization energies, a level of accuracy that is competitive with more expensive many-body techniques. An empirical shifting procedure with one adjustable parameter reduces the average error below 0.2 eV. This shifted Slater transition method is a simple and practical way to compute core-level binding energies using only initial-state Kohn-Sham eigenvalues. It requires no more computational effort than ΔSCF and may be especially useful for simulating transient x-ray experiments where core-level spectroscopy is used to probe an excited electronic state, for which the ΔSCF approach requires a tedious state-by-state calculation of the spectrum. As an example, we use Slater-type methods to model x-ray emission spectroscopy.
Collapse
Affiliation(s)
- Subrata Jana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
9
|
Rana B, Beran GJO, Herbert JM. Correcting π-delocalisation errors in conformational energies using density-corrected DFT, with application to crystal polymorphs. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2138789] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/31/2022]
Affiliation(s)
- Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| | | | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
10
|
Janesko BG. Unification of Perdew-Zunger self-interaction correction, DFT+U, and Rung 3.5 density functionals. J Chem Phys 2022; 157:151101. [PMID: 36272781 DOI: 10.1063/5.0109338] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This Communication presents a unified derivation of three different approximations used in density functional theory (DFT): the Perdew-Zunger self-interaction correction (PZSIC), the Hubbard correction DFT+U, and the Rung 3.5 density functionals. All three approximations can be derived by introducing electron self-interaction into the Kohn-Sham (KS) reference system of noninteracting electrons. The derivation uses the Adiabatic Projection formalism: one projects the electron-electron interaction operator onto certain states, introduces the projected operator into the reference system, and defines a density functional for the remainder. Projecting onto individual localized KS orbitals recovers our previous derivation of the PZSIC [B. G. Janesko, J. Phys. Chem. Lett. 13, 5698-5702 (2022)]. Projecting onto localized atom-centered orbitals recovers a variant of DFT+U. Projecting onto localized states at each point in space recovers Rung 3.5 approaches. New results include an "atomic state PZSIC" that does not require localizing the KS orbitals, a demonstration that typical Hubbard U parameters reproduce a scaled-down PZSIC, and a Rung 3.5 variant of DFT+U that does not require choosing atom-dependent states.
Collapse
Affiliation(s)
- Benjamin G Janesko
- Department of Chemistry and Biochemistry, Texas Christian University, 2800 S. University Dr., Fort Worth, Texas 76109, USA
| |
Collapse
|
11
|
Mahler A, Williams J, Su NQ, Yang W. Localized orbital scaling correction for periodic systems. PHYSICAL REVIEW. B 2022; 106:035147. [PMID: 37727592 PMCID: PMC10508887 DOI: 10.1103/physrevb.106.035147] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
Density functional theory offers accurate structure prediction at acceptable computational cost, but commonly used approximations suffer from delocalization error; this results in inaccurate predictions of quantities such as energy band gaps of finite and bulk systems, energy level alignments, and electron distributions at interfaces. The localized orbital scaling correction (LOSC) was developed to correct delocalization error by using orbitals localized in space and energy. These localized orbitals span both the occupied and unoccupied spaces and can have fractional occupations in order to correct both the total energy and the one-electron energy eigenvalues. We extend the LOSC method to periodic systems, in which the localized orbitals employed are dually localized Wannier functions. In light of the effect of the bulk environment on the electrostatic interaction between localized orbitals, we modify the LOSC energy correction to include a screened Coulomb kernel. For a test set of semiconductors and large-gap insulators, we show that the screened LOSC method consistently improves the band gap compared to the parent density functional approximation.
Collapse
Affiliation(s)
- Aaron Mahler
- Department of Physics, Duke University, Durham, North Carolina 27708, USA
| | - Jacob Williams
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Neil Qiang Su
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
- Department of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education) and Renewable Energy Conversion and Storage Center (RECAST), Nankai University, Tianjin 300071, China
| | - Weitao Yang
- Department of Physics, Duke University, Durham, North Carolina 27708, USA
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
12
|
Bajaj A, Duan C, Nandy A, Taylor MG, Kulik HJ. Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry. J Chem Phys 2022; 156:184112. [DOI: 10.1063/5.0089460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d1 MOPs to be suitable for nearly eliminating all energetic delocalization error and static correlation error. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization error and static correlation error within a non-empirical framework.
Collapse
Affiliation(s)
- Akash Bajaj
- Massachusetts Institute of Technology, United States of America
| | - Chenru Duan
- Massachusetts Institute of Technology, United States of America
| | - Aditya Nandy
- Massachusetts Institute of Technology, United States of America
| | | | - Heather J. Kulik
- Dept of Chemical Engineering, Massachusetts Institute of Technology, United States of America
| |
Collapse
|
13
|
Prokopiou G, Hartstein M, Govind N, Kronik L. Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals. J Chem Theory Comput 2022; 18:2331-2340. [PMID: 35369687 PMCID: PMC9009176 DOI: 10.1021/acs.jctc.2c00082] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Indexed: 11/29/2022]
Abstract
We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.
Collapse
Affiliation(s)
- Georgia Prokopiou
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Michal Hartstein
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Leeor Kronik
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
14
|
Nandy A, Duan C, Taylor MG, Liu F, Steeves AH, Kulik HJ. Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning. Chem Rev 2021; 121:9927-10000. [PMID: 34260198 DOI: 10.1021/acs.chemrev.1c00347] [Citation(s) in RCA: 70] [Impact Index Per Article: 23.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
15
|
Mei Y, Chen Z, Yang W. Exact Second-Order Corrections and Accurate Quasiparticle Energy Calculations in Density Functional Theory. J Phys Chem Lett 2021; 12:7236-7244. [PMID: 34310157 PMCID: PMC9367128 DOI: 10.1021/acs.jpclett.1c01962] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We develop a second-order correction to commonly used density functional approximations (DFAs) to eliminate the systematic delocalization error. The method, based on the previously developed global scaling correction (GSC), is an exact quadratic correction to the DFA for the fractional charge behavior and uses the analytical second derivatives of the total energy with respect to fractional occupation numbers of the canonical molecular orbitals. For small and medium-size molecules, this correction leads to ground-state orbital energies that are a highly accurate approximation to the corresponding quasiparticle energies. It provides excellent predictions of ionization potentials, electron affinities, photoemission spectrum, and photoexcitation energies beyond previous approximate second-order approaches, thus showing potential for broad applications in computational spectroscopy.
Collapse
Affiliation(s)
- Yuncai Mei
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Zehua Chen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
- Department of Physics, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
16
|
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
|
17
|
Bajaj A, Kulik HJ. Molecular DFT+U: A Transferable, Low-Cost Approach to Eliminate Delocalization Error. J Phys Chem Lett 2021; 12:3633-3640. [PMID: 33826346 DOI: 10.1021/acs.jpclett.1c00796] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
While density functional theory (DFT) is widely applied for its combination of cost and accuracy, corrections (e.g., DFT+U) that improve it are often needed to tackle correlated transition-metal chemistry. In principle, the functional form of DFT+U, consisting of a set of localized atomic orbitals (AOs) and a quadratic energy penalty for deviation from integer occupations of those AOs, enables the recovery of the exact conditions of piecewise linearity and the derivative discontinuity. Nevertheless, for practical transition-metal complexes, where both atomic states and ligand orbitals participate in bonding, standard DFT+U can fail to eliminate delocalization error (DE). Here, we show that by introducing an alternative valence-state (i.e., molecular orbital or MO) basis to the DFT+U approach, we recover exact conditions in cases for which standard DFT+U corrections have no error-reducing effect. This MO-based DFT+U also eliminates DE where standard AO-based DFT+U is already successful. We demonstrate the transferability of our approach on representative transition-metal complexes with a range of ligand field strengths, electron configurations (i.e., from Sc to Zn), and spin states.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
18
|
Affiliation(s)
- Heather J. Kulik
- Department of Chemical Engineering Massachusetts Institute of Technology 77 Massachusetts Ave Rm 66–464 Cambridge MA 02139 USA
| |
Collapse
|
19
|
Ramos C, Janesko BG. Nonlocal rung-3.5 correlation from the density matrix expansion: Flat-plane condition, thermochemistry, and kinetics. J Chem Phys 2020; 153:164116. [PMID: 33138396 DOI: 10.1063/5.0025160] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
The rung-3.5 approach to density functional theory constructs nonlocal approximate correlation from the expectation values of nonlocal one-electron operators. This offers an inexpensive solution to hybrid functionals' imbalance between exact nonlocal exchange and local approximate correlation. Our rung-3.5 correlation functionals also include a local complement to the nonlocal ingredient, analogous to the local exchange component of a hybrid functional. Here, we use the density matrix expansion (DME) to build rung-3.5 complements. We demonstrate how these provide a measure of local fractional occupancy and use them to approximate the flat-plane condition. We also use these complements in a three-parameter nonlocal correlation functional compatible with full nonlocal exchange. This functional approaches the accuracy of widely used hybrids for molecular thermochemistry and kinetics. The DME provides a foundation for practical, minimally empirical, nonlocal correlation functionals compatible with full nonlocal local exchange.
Collapse
Affiliation(s)
- Chloe Ramos
- Department of Chemistry and Biochemistry, Texas Christian University, 2800 S. University Dr., Fort Worth, Texas 76129, USA
| | - Benjamin G Janesko
- Department of Chemistry and Biochemistry, Texas Christian University, 2800 S. University Dr., Fort Worth, Texas 76129, USA
| |
Collapse
|
20
|
Dewo W, Łuczyńska K, Zorenko Y, Gorbenko V, Drużbicki K, Runka T. In silico Raman spectroscopy of YAlO 3 single-crystalline film. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2020; 231:118111. [PMID: 32035331 DOI: 10.1016/j.saa.2020.118111] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Revised: 01/23/2020] [Accepted: 01/24/2020] [Indexed: 06/10/2023]
Abstract
The Raman response of the YAlO3 (YAP) perovskite is modeled by means of periodic density functional theory. A number of different approximations to the exchange-correlation functional are benchmarked against the structural and spectroscopic data as imposing all-electron Gaussian-type basis sets. The WC1LYP functional was found to be superior, particularly outperforming other tested approaches in the prediction of the local structure of the AlO subunits, which reflects in the observed lattice-dynamics. The Raman response is further decomposed into the directional spectra, which are due to different components of the polarizability tensor, and confronted with the experimental Raman spectra, recorded in different scattering geometries of the single-crystalline film of YAP. The in silico lattice dynamics provides the unequivocal assignment of the observed bands with an excellent match to the experimental spectra, allowing for a complete analysis of the underlying phonon modes in terms of their energy, symmetry and the directional activity. The presented analysis serves as a high-quality reference, potentially useful in the future studies of other YAP materials, where Raman spectroscopy along with the X-Ray diffraction is the first method of choice.
Collapse
Affiliation(s)
- Wioletta Dewo
- Faculty of Technical Physics, Poznan University of Technology, Piotrowo 3, 60-965 Poznan, Poland
| | - Katarzyna Łuczyńska
- Department of Radiospectroscopy, Faculty of Physics, Adam Mickiewicz University, Uniwersytetu Poznanskiego 2, 61-614 Poznan, Poland; Frank Laboratory of Neutron Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russian Federation.
| | - Yuriy Zorenko
- Institute of Physics, Kazimierz Wielki University, Weyssenhoffa 11, 85-090 Bydgoszcz, Poland
| | - Vitaliy Gorbenko
- Institute of Physics, Kazimierz Wielki University, Weyssenhoffa 11, 85-090 Bydgoszcz, Poland
| | - Kacper Drużbicki
- Department of Radiospectroscopy, Faculty of Physics, Adam Mickiewicz University, Uniwersytetu Poznanskiego 2, 61-614 Poznan, Poland; Frank Laboratory of Neutron Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russian Federation
| | - Tomasz Runka
- Faculty of Technical Physics, Poznan University of Technology, Piotrowo 3, 60-965 Poznan, Poland.
| |
Collapse
|
21
|
Liu F, Kulik HJ. Impact of Approximate DFT Density Delocalization Error on Potential Energy Surfaces in Transition Metal Chemistry. J Chem Theory Comput 2019; 16:264-277. [DOI: 10.1021/acs.jctc.9b00842] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
22
|
Zhao Q, Kulik HJ. Stable Surfaces That Bind Too Tightly: Can Range-Separated Hybrids or DFT+U Improve Paradoxical Descriptions of Surface Chemistry? J Phys Chem Lett 2019; 10:5090-5098. [PMID: 31411023 PMCID: PMC6748670 DOI: 10.1021/acs.jpclett.9b01650] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2019] [Accepted: 08/14/2019] [Indexed: 05/25/2023]
Abstract
Approximate, semilocal density functional theory (DFT) suffers from delocalization error that can lead to a paradoxical model of catalytic surfaces that both overbind adsorbates yet are also too stable. We investigate the effect of two widely applied approaches for delocalization error correction, (i) affordable DFT+U (i.e., semilocal DFT augmented with a Hubbard U) and (ii) hybrid functionals with an admixture of Hartree-Fock (HF) exchange, on surface and adsorbate energies across a range of rutile transition metal oxides widely studied for their promise as water-splitting catalysts. We observe strongly row- and period-dependent trends with DFT+U, which increases surface formation energies only in early transition metals (e.g., Ti and V) and decreases adsorbate energies only in later transition metals (e.g., Ir and Pt). Both global and local hybrids destabilize surfaces and reduce adsorbate binding across the periodic table, in agreement with higher-level reference calculations. Density analysis reveals why hybrid functionals correct both quantities, whereas DFT+U does not. We recommend local, range-separated hybrids for the accurate modeling of catalysis in transition metal oxides at only a modest increase in computational cost over semilocal DFT.
Collapse
Affiliation(s)
- Qing Zhao
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Mechanical Engineering, Massachusetts
Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
23
|
Shishkin M, Sato H. DFT+U in Dudarev's formulation with corrected interactions between the electrons with opposite spins: The form of Hamiltonian, calculation of forces, and bandgap adjustments. J Chem Phys 2019; 151:024102. [PMID: 31301721 DOI: 10.1063/1.5090445] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Hubbard corrected density functional theory (DFT) methods, such as the DFT+U approach in Dudarev's approximation, are widely used for the description of energetics and electronic structure of strongly correlated materials, providing higher level of accuracy than local DFT calculations (e.g., local density approximation or generalized gradient approximation). However, the DFT+U method in Dudarev's formulation limits the introduced corrections to interactions between the electrons within the same spin channel, whereas interactions between the electrons with opposite spins are still treated using local DFT functional (e.g., Perdew-Burke-Ernzerhof). In recent years, the need for correction of these interactions between the electrons with opposite spins has been recognized and additional terms have been added to the Hubbard term to reflect it. Although such extended DFT+U functionals have been proposed, the form of respective Hamiltonian operator, defined as a total energy derivative over density with appropriate treatment of double counting corrections due to additional Hubbard terms, has not been explicitly presented. In this work, we provide an expression for such a type of Hamiltonian, which contains the respective double counting correction contributions. This formulation also allows evaluation of atomic forces, using computational settings discussed herein. In addition, we also introduce adjustments for too narrow theoretical bandgaps, using scissor operator technique. This allows for a greater level of corrections of energetics and magnetic properties of studied transition metal compounds, avoiding possible unphysical overlap between occupied and unoccupied electronic bands.
Collapse
Affiliation(s)
- M Shishkin
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Nishikyo-ku, Kyoto 615-8520, Japan
| | - H Sato
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Nishikyo-ku, Kyoto 615-8520, Japan
| |
Collapse
|
24
|
Bajaj A, Liu F, Kulik HJ. Non-empirical, low-cost recovery of exact conditions with model-Hamiltonian inspired expressions in jmDFT. J Chem Phys 2019; 150:154115. [PMID: 31005112 DOI: 10.1063/1.5091563] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Density functional theory (DFT) is widely applied to both molecules and materials, but well known energetic delocalization and static correlation errors in practical exchange-correlation approximations limit quantitative accuracy. Common methods that correct energetic delocalization errors, such as the Hubbard U correction in DFT+U or Hartree-Fock exchange in global hybrids, do so at the cost of worsening static correlation errors. We recently introduced an alternate approach [Bajaj et al., J. Chem. Phys. 147, 191101 (2017)] known as judiciously modified DFT (jmDFT), wherein the deviation from exact behavior of semilocal functionals over both fractional spin and charge, i.e., the so-called flat plane, was used to motivate functional forms of second order analytic corrections. In this work, we introduce fully nonempirical expressions for all four coefficients in a DFT+U+J-inspired form of jmDFT, where all coefficients are obtained only from energies and eigenvalues of the integer-electron systems. We show good agreement for U and J coefficients obtained nonempirically as compared with the results of numerical fitting in a jmDFT U+J/J' correction. Incorporating the fully nonempirical jmDFT correction reduces and even eliminates the fractional spin error at the same time as eliminating the energetic delocalization error. We show that this approach extends beyond s-electron systems to higher angular momentum cases including p- and d-electrons. Finally, we diagnose some shortcomings of the current jmDFT approach that limit its ability to improve upon DFT results for cases such as weakly bound anions due to poor underlying semilocal functional behavior.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| |
Collapse
|
25
|
Janet JP, Liu F, Nandy A, Duan C, Yang T, Lin S, Kulik HJ. Designing in the Face of Uncertainty: Exploiting Electronic Structure and Machine Learning Models for Discovery in Inorganic Chemistry. Inorg Chem 2019; 58:10592-10606. [PMID: 30834738 DOI: 10.1021/acs.inorgchem.9b00109] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Recent transformative advances in computing power and algorithms have made computational chemistry central to the discovery and design of new molecules and materials. First-principles simulations are increasingly accurate and applicable to large systems with the speed needed for high-throughput computational screening. Despite these strides, the combinatorial challenges associated with the vastness of chemical space mean that more than just fast and accurate computational tools are needed for accelerated chemical discovery. In transition-metal chemistry and catalysis, unique challenges arise. The variable spin, oxidation state, and coordination environments favored by elements with well-localized d or f electrons provide great opportunity for tailoring properties in catalytic or functional (e.g., magnetic) materials but also add layers of uncertainty to any design strategy. We outline five key mandates for realizing computationally driven accelerated discovery in inorganic chemistry: (i) fully automated simulation of new compounds, (ii) knowledge of prediction sensitivity or accuracy, (iii) faster-than-fast property prediction methods, (iv) maps for rapid chemical space traversal, and (v) a means to reveal design rules on the kilocompound scale. Through case studies in open-shell transition-metal chemistry, we describe how advances in methodology and software in each of these areas bring about new chemical insights. We conclude with our outlook on the next steps in this process toward realizing fully autonomous discovery in inorganic chemistry using computational chemistry.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Fang Liu
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Aditya Nandy
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Chenru Duan
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Sean Lin
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
26
|
Hait D, Head-Gordon M. Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number. J Phys Chem Lett 2018; 9:6280-6288. [PMID: 30339010 DOI: 10.1021/acs.jpclett.8b02417] [Citation(s) in RCA: 65] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C12H14 polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
27
|
Describing strong correlation with fractional-spin correction in density functional theory. Proc Natl Acad Sci U S A 2018; 115:9678-9683. [PMID: 30201706 DOI: 10.1073/pnas.1807095115] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
An effective fractional-spin correction is developed to describe static/strong correlation in density functional theory. Combined with the fractional-charge correction from recently developed localized orbital scaling correction (LOSC), a functional, the fractional-spin LOSC (FSLOSC), is proposed. FSLOSC, a correction to commonly used functional approximations, introduces the explicit derivative discontinuity and largely restores the flat-plane behavior of electronic energy at fractional charges and fractional spins. In addition to improving results from conventional functionals for the prediction of ionization potentials, electron affinities, quasiparticle spectra, and reaction barrier heights, FSLOSC properly describes the dissociation of ionic species, single bonds, and multiple bonds without breaking space or spin symmetry and corrects the spurious fractional-charge dissociation of heteroatom molecules of conventional functionals. Thus, FSLOSC demonstrates success in reducing delocalization error and including strong correlation, within low-cost density functional approximation.
Collapse
|
28
|
Withanage KPK, Trepte K, Peralta JE, Baruah T, Zope R, Jackson KA. On the Question of the Total Energy in the Fermi–Löwdin Orbital Self-Interaction Correction Method. J Chem Theory Comput 2018; 14:4122-4128. [DOI: 10.1021/acs.jctc.8b00344] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kushantha P. K. Withanage
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| | - Kai Trepte
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| | - Juan E. Peralta
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| | - Tunna Baruah
- Physics Department, University of Texas-El Paso, El Paso, Texas 79968, United States
| | - Rajendra Zope
- Physics Department, University of Texas-El Paso, El Paso, Texas 79968, United States
| | - Koblar A. Jackson
- Physics Department and Science of Advanced Materials Ph.D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, United States
| |
Collapse
|
29
|
Margraf JT, Bartlett R. Communication: Coupled cluster and many-body perturbation theory for fractional charges and spins. J Chem Phys 2018; 148:221103. [DOI: 10.1063/1.5040164] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Johannes T. Margraf
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Rodney Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|