1
|
Nguyen DB, Jackson KA, Peralta JE. Bond length alternation of π-conjugated polymers predicted by the Fermi-Löwdin orbital self-interaction correction method. J Chem Phys 2024; 160:014101. [PMID: 38165094 DOI: 10.1063/5.0178251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 12/07/2023] [Indexed: 01/03/2024] Open
Abstract
π-conjugated polymers have been used in a wide range of practical applications, partly due to their unique properties that originate in the delocalization of electrons through the polymer backbone. The level of delocalization can be characterized by the induced bond length alternation (BLA), with shorter BLA connected with strong delocalization and vice versa. The accurate description of this structural parameter can be considered a benchmark for testing the capability of different electronic structure methods for self-interaction error (SIE) removal and electron correlation inclusion. Density functional theory (DFT), in its local or semi-local flavors, suffers from SIE and, thus, underestimates the BLA compared to self-interaction-free methods. In this work, we utilize the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for one-electron self-interaction removal to characterize the BLA of five oligomers with increasing length extrapolated to the polymeric limit. We compare the self-interaction-free BLA to several DFT approximations, Møller-Plesset second-order perturbation theory (MP2), and the BLA obtained with the domain based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] approximation. Our findings show that FLOSIC corrects for the small BLA given by (semi-)local DFT approximations, but it tends to overcorrect with respect to CAM-B3LYP, MP2, and DLPNO-CCSD(T).
Collapse
Affiliation(s)
- Duyen B Nguyen
- Department of Physics and Science of Advanced Materials, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Koblar A Jackson
- Department of Physics and Science of Advanced Materials, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Juan E Peralta
- Department of Physics and Science of Advanced Materials, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| |
Collapse
|
2
|
Pederson MR, Withanage KPK, Hooshmand Z, Johnson AI, Baruah T, Yamamoto Y, Zope RR, Kao DY, Shukla PB, Johnson JK, Peralta JE, Jackson KA. Use of FLOSIC for understanding anion-solvent interactions. J Chem Phys 2023; 159:154112. [PMID: 37861122 DOI: 10.1063/5.0172300] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 09/20/2023] [Indexed: 10/21/2023] Open
Abstract
An Achille's heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3-, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth's atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi-Löwdin-Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree-Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.
Collapse
Affiliation(s)
- Mark R Pederson
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | | | - Zahra Hooshmand
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Alex I Johnson
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Yoh Yamamoto
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Der-You Kao
- NASA Postdoctoral Program, NASA Goddard Space Flight Center, Greenbelt, Maryland 20771, USA
| | - Priyanka B Shukla
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - J Karl Johnson
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Juan E Peralta
- Physics Department, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Koblar A Jackson
- Physics Department, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| |
Collapse
|
3
|
Pederson MR, Johnson AI, Withanage KPK, Dolma S, Flores GB, Hooshmand Z, Khandal K, Lasode PO, Baruah T, Jackson KA. Downward quantum learning from element 118: Automated generation of Fermi-Löwdin orbitals for all atoms. J Chem Phys 2023; 158:084101. [PMID: 36859080 DOI: 10.1063/5.0135089] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi-Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew-Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N - 1 and a maximum of 2N - 2 initial starting points for systems composed of a smaller number of electrons. Applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.
Collapse
Affiliation(s)
- Mark R Pederson
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Alexander I Johnson
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | | | - Sherab Dolma
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Gustavo Bravo Flores
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Zahra Hooshmand
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Kusal Khandal
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Peter O Lasode
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Department of Physics, The 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
| |
Collapse
|
4
|
Ruan S, Jackson KA, Ruzsinszky A. Spin-crossover complexes: Self-interaction correction vs density correction. J Chem Phys 2023; 158:064303. [PMID: 36792493 DOI: 10.1063/5.0128950] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
Complexes containing a transition metal atom with a 3d4-3d7 electron configuration typically have two low-lying, high-spin (HS) and low-spin (LS) states. The adiabatic energy difference between these states, known as the spin-crossover energy, is small enough to pose a challenge even for electronic structure methods that are well known for their accuracy and reliability. In this work, we analyze the quality of electronic structure approximations for spin-crossover energies of iron complexes with four different ligands by comparing energies from self-consistent and post-self-consistent calculations for methods based on the random phase approximation and the Fermi-Löwdin self-interaction correction. Considering that Hartree-Fock densities were found by Song et al., J. Chem. Theory Comput. 14, 2304 (2018), to eliminate the density error to a large extent, and that the Hartree-Fock method and the Perdew-Zunger-type self-interaction correction share some physics, we compare the densities obtained with these methods to learn their resemblance. We find that evaluating non-empirical exchange-correlation energy functionals on the corresponding self-interaction-corrected densities can mitigate the strong density errors and improves the accuracy of the adiabatic energy differences between HS and LS states.
Collapse
Affiliation(s)
- Shiqi Ruan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Koblar A Jackson
- Physics Department and Science of Advanced Materials Ph. D. Program, Central Michigan University, Mount Pleasant, Michigan 48858, USA
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
5
|
Yamamoto Y, Baruah T, Chang PH, Romero S, Zope RR. Self-consistent implementation of locally scaled self-interaction-correction method. J Chem Phys 2023; 158:064114. [PMID: 36792502 DOI: 10.1063/5.0130436] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
Recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys. 151, 214108 (2019)] is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew-Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn-Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional [e.g., Perdew-Burke-Ernzerhof (PBE)] and barrier heights of the BH6 database better than some of the most widely used hybrid functionals (e.g., PBE0 and B3LYP). The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). This work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.
Collapse
Affiliation(s)
- Yoh Yamamoto
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Po-Hao Chang
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Selim Romero
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
6
|
Melo JI, Pederson MR, Peralta JE. Density Matrix Implementation of the Fermi-Löwdin Orbital Self-Interaction Correction Method. J Phys Chem A 2023; 127:527-534. [PMID: 36598275 DOI: 10.1021/acs.jpca.2c07646] [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/2023]
Abstract
The Fermi-Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi-Löwdin orbitals which are used to remove the self-interaction error in the Perdew-Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi-Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn-Sham with one additional term in the Kohn-Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn-Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. We analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville-von Neumann equation to calculate excitation energies.
Collapse
Affiliation(s)
- Juan I Melo
- Facultad de Ciencias Exactas y Naturales, Departamento de Física, Universidad de Buenos Aires, Buenos Aires1428, Argentina.,CONICET - Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Buenos Aires1428, Argentina
| | - Mark R Pederson
- Department of Physics, the University of Texas at El Paso, El Paso, Texas79968, United States
| | - Juan E Peralta
- Department of Physics, Central Michigan University, Mount Pleasant, Michigan48859, United States
| |
Collapse
|
7
|
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
|
8
|
Trepte K, Schwalbe S, Liebing S, Schulze WT, Kortus J, Myneni H, Ivanov AV, Lehtola S. Chemical bonding theories as guides for self-interaction corrected solutions: Multiple local minima and symmetry breaking. J Chem Phys 2021; 155:224109. [PMID: 34911315 DOI: 10.1063/5.0071796] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Fermi-Löwdin orbitals (FLOs) are a special set of localized orbitals, which have become commonly used in combination with the Perdew-Zunger self-interaction correction (SIC) in the FLO-SIC method. The FLOs are obtained for a set of occupied orbitals by specifying a classical position for each electron. These positions are known as Fermi-orbital descriptors (FODs), and they have a clear relation to chemical bonding. In this study, we show how FLOs and FODs can be used to initialize, interpret, and justify SIC solutions in a common chemical picture, both within FLO-SIC and in traditional variational SIC, and to locate distinct local minima in either of these approaches. We demonstrate that FLOs based on Lewis theory lead to symmetry breaking for benzene-the electron density is found to break symmetry already at the symmetric molecular structure-while ones from Linnett's double-quartet theory reproduce symmetric electron densities and molecular geometries. Introducing a benchmark set of 16 planar cyclic molecules, we show that using Lewis theory as the starting point can lead to artifactual dipole moments of up to 1 D, while Linnett SIC dipole moments are in better agreement with experimental values. We suggest using the dipole moment as a diagnostic of symmetry breaking in SIC and monitoring it in all SIC calculations. We show that Linnett structures can often be seen as superpositions of Lewis structures and propose Linnett structures as a simple way to describe aromatic systems in SIC with reduced symmetry breaking. The role of hovering FODs is also briefly discussed.
Collapse
Affiliation(s)
- Kai Trepte
- SUNCAT Center for Interface Science and Catalysis, Stanford University, Menlo Park, California 94025, USA
| | - Sebastian Schwalbe
- Institute of Theoretical Physics, TU Bergakademie Freiberg, D-09599 Freiberg, Germany
| | - Simon Liebing
- Joint Institute for Nuclear Research Dubna, Bogoliubov Laboratory of Theoretical Physics, 141980 Dubna, Russia
| | - Wanja T Schulze
- Institute of Theoretical Physics, TU Bergakademie Freiberg, D-09599 Freiberg, Germany
| | - Jens Kortus
- Institute of Theoretical Physics, TU Bergakademie Freiberg, D-09599 Freiberg, Germany
| | - Hemanadhan Myneni
- Science Institute and Faculty of Physical Sciences, VR-III, University of Iceland, 107 Reykjavík, Iceland
| | - Aleksei V Ivanov
- Science Institute and Faculty of Physical Sciences, VR-III, University of Iceland, 107 Reykjavík, Iceland
| | - Susi Lehtola
- Molecular Sciences Software Institute, Blacksburg, Virginia 24061, USA
| |
Collapse
|
9
|
Nguyen DB, Pederson MR, Perdew JP, Jackson KA, Peralta JE. Initial Fermi orbital descriptors for FLOSIC calculations: The quick-FOD method. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
10
|
Karanovich A, Yamamoto Y, Jackson KA, Park K. Electronic structure of mononuclear Cu-based molecule from density-functional theory with self-interaction correction. J Chem Phys 2021; 155:014106. [PMID: 34241401 DOI: 10.1063/5.0054439] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the electronic structure of a planar mononuclear Cu-based molecule [Cu(C6H4S2)2]z in two oxidation states (z = -2, -1) using density-functional theory (DFT) with Fermi-Löwdin orbital (FLO) self-interaction correction (SIC). The dianionic Cu-based molecule was proposed to be a promising qubit candidate. Self-interaction error within approximate DFT functionals renders severe delocalization of electron and spin densities arising from 3d orbitals. The FLO-SIC method relies on optimization of Fermi-Löwdin orbital descriptors (FODs) with which localized occupied orbitals are constructed to create SIC potentials. Starting with many initial sets of FODs, we employ a frozen-density loop algorithm within the FLO-SIC method to study the Cu-based molecule. We find that the electronic structure of the molecule remains unchanged despite somewhat different final FOD configurations. In the dianionic state (spin S = 1/2), FLO-SIC spin density originates from the Cu d and S p orbitals with an approximate ratio of 2:1, in quantitative agreement with multireference calculations, while in the case of SIC-free DFT, the orbital ratio is reversed. Overall, FLO-SIC lowers the energies of the occupied orbitals and, in particular, the 3d orbitals unhybridized with the ligands significantly, which substantially increases the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) compared to SIC-free DFT results. The FLO-SIC HOMO-LUMO gap of the dianionic state is larger than that of the monoanionic state, which is consistent with experiment. Our results suggest a positive outlook of the FLO-SIC method in the description of magnetic exchange coupling within 3d-element-based systems.
Collapse
Affiliation(s)
- Anri Karanovich
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Yoh Yamamoto
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Koblar Alan Jackson
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Kyungwha Park
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
11
|
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: 28] [Impact Index Per Article: 9.3] [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
|
12
|
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
|
13
|
Bhattarai P, Santra B, Wagle K, Yamamoto Y, Zope RR, Ruzsinszky A, Jackson KA, Perdew JP. Exploring and enhancing the accuracy of interior-scaled Perdew–Zunger self-interaction correction. J Chem Phys 2021; 154:094105. [DOI: 10.1063/5.0041646] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Puskar Bhattarai
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Kamal Wagle
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Yoh Yamamoto
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R. Zope
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Koblar A. Jackson
- Department of Physics and Science of Advanced Materials, 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
|
14
|
Adhikari S, Santra B, Ruan S, Bhattarai P, Nepal NK, Jackson KA, Ruzsinszky A. The Fermi–Löwdin self-interaction correction for ionization energies of organic molecules. J Chem Phys 2020; 153:184303. [DOI: 10.1063/5.0024776] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Santosh Adhikari
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Shiqi Ruan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Puskar Bhattarai
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Niraj K. Nepal
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Koblar A. Jackson
- Department of Physics and Science of Advanced Materials Program, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
15
|
Li L, Trepte K, Jackson KA, Johnson JK. Application of Self-Interaction Corrected Density Functional Theory to Early, Middle, and Late Transition States. J Phys Chem A 2020; 124:8223-8234. [DOI: 10.1021/acs.jpca.0c06485] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lin Li
- Department of Chemical and Petroleum Engineering University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Kai Trepte
- Department of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, United States
| | - Koblar A. Jackson
- Department of Physics and Science of Advanced Materials Program, Central Michigan University, Mount Pleasant, Michigan 48859, United States
| | - J. Karl Johnson
- Department of Chemical and Petroleum Engineering University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
16
|
Schwalbe S, Fiedler L, Kraus J, Kortus J, Trepte K, Lehtola S. PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction. J Chem Phys 2020; 153:084104. [DOI: 10.1063/5.0012519] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sebastian Schwalbe
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Lenz Fiedler
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Jakob Kraus
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Jens Kortus
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Kai Trepte
- Department of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Susi Lehtola
- Department of Chemistry, University of Helsinki, P.O. Box 55 (A. I. Virtasen Aukio 1), FI-00014 University of Helsinki, Finland
| |
Collapse
|
17
|
Johnson AI, Withanage KPK, Sharkas K, Yamamoto Y, Baruah T, Zope RR, Peralta JE, Jackson KA. The effect of self-interaction error on electrostatic dipoles calculated using density functional theory. J Chem Phys 2019; 151:174106. [PMID: 31703485 DOI: 10.1063/1.5125205] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Spurious electron self-interaction in density functional approximations (DFAs) can lead to inaccurate predictions of charge transfer in heteronuclear molecules that manifest as errors in calculated electrostatic dipoles. Here, we show the magnitude of these errors on dipoles computed for a diverse set of 47 molecules taken from the recent benchmark study of Hait and Head-Gordon [J. Chem. Theory Comput. 14, 1969 (2018)]. We compare the results of Perdew-Wang local spin density approximation (PW92), Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and strongly constrained and appropriately normed (SCAN) meta-GGA dipole calculations, along with those of their respective self-interaction-corrected (SIC) counterparts, to reference values from accurate wave function-based methods. The SIC calculations were carried out using the Fermi-Löwdin orbital (FLO-SIC) approach. We find that correcting for self-interaction generally increases the degree of charge transfer, thereby increasing the size of calculated dipole moments. The FLO-SIC-PW92 and FLO-SIC-PBE dipoles are in better agreement with reference values than their uncorrected DFA counterparts, particularly for strongly ionic molecules where significant improvement is seen. Applying FLO-SIC to SCAN does not improve dipole values overall. We also show that removing self-interaction improves the description of the dipole for stretched-bond geometries and recovers the physically correct separated atom limit of zero dipole. Finally, we find that the best agreement between the FLO-SIC-DFA and reference dipoles occurs when the molecular geometries are optimized using the FLO-SIC-DFA.
Collapse
Affiliation(s)
- Alexander I Johnson
- Physics Department, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Kushantha P K Withanage
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Kamal Sharkas
- Physics Department, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Yoh Yamamoto
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Juan E Peralta
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Koblar A Jackson
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| |
Collapse
|