1
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
2
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
3
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
4
|
Mishra P, Yamamoto Y, Chang PH, Nguyen DB, Peralta JE, Baruah T, Zope RR. Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods. J Phys Chem A 2022; 126:1923-1935. [PMID: 35302373 DOI: 10.1021/acs.jpca.1c10354] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We examine the role of self-interaction error (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular, we analyze the effect of scaling down the self-interaction correction (SIC) for three nonempirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and the recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger SIC [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048.], the orbitalwise scaled SIC method [Vydrov, O. A. et al. J. Chem. Phys. 2006, 124, 094108.], and the recent local scaling method [Zope, R. R. et al. J. Chem. Phys. 2019, 151, 214108.]. We compute the magnetic exchange coupling constants using the spin projection and nonprojection approaches for sets of molecules composed of dinuclear and polynuclear H···He models, organic radical molecules, and chlorocuprate and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single-electron regions, PZSIC performs well, but for more complex organic systems and the chlorocuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases, LSIC with kinetic energy density ratio performs better than PZSIC. Analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.
Collapse
Affiliation(s)
- Prakash Mishra
- Computational Science Program, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Yoh Yamamoto
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Po-Hao Chang
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Duyen B Nguyen
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, United States
| | - Juan E Peralta
- Physics Department and Science of Advanced Materials Program, Central Michigan University, Mt. Pleasant, Michigan 48859, 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
| |
Collapse
|
5
|
Akter S, Vargas JA, Sharkas K, Peralta JE, Jackson KA, Baruah T, Zope RR. How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations? Phys Chem Chem Phys 2021; 23:18678-18685. [PMID: 34612405 DOI: 10.1039/d0cp06512a] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys., 2018, 20, 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew-Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi-Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys., 2019, 151, 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.
Collapse
Affiliation(s)
- Sharmin Akter
- Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA.
| | | | | | | | | | | | | |
Collapse
|
6
|
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
|
7
|
Diaz CM, Suryanarayana P, Xu Q, Baruah T, Pask JE, Zope RR. Implementation of Perdew–Zunger self-interaction correction in real space using Fermi–Löwdin orbitals. J Chem Phys 2021; 154:084112. [DOI: 10.1063/5.0031341] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Carlos M. Diaz
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Phanish Suryanarayana
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Qimen Xu
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Tunna Baruah
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - John E. Pask
- Physics Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Rajendra R. Zope
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
8
|
Abstract
The local self-interaction correction method with a simple scaling factor performs better than the Perdew-Zunger self-interaction correction method and also provides a good description of the binding energies of weakly bonded water clusters.
Collapse
Affiliation(s)
- Selim Romero
- Department of Physics
- University of Texas at El Paso
- El Paso
- USA
- Computational Science Program
| | - Yoh Yamamoto
- Department of Physics
- University of Texas at El Paso
- El Paso
- USA
| | - Tunna Baruah
- Department of Physics
- University of Texas at El Paso
- El Paso
- USA
- Computational Science Program
| | - Rajendra R. Zope
- Department of Physics
- University of Texas at El Paso
- El Paso
- USA
- Computational Science Program
| |
Collapse
|
9
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
10
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
11
|
Yamamoto Y, Salcedo A, Diaz CM, Alam MS, Baruah T, Zope RR. Assessing the effect of regularization on the molecular properties predicted by SCAN and self-interaction corrected SCAN meta-GGA. Phys Chem Chem Phys 2020; 22:18060-18070. [PMID: 32760934 DOI: 10.1039/d0cp02717k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent regularization of the SCAN meta-GGA functional (rSCAN) has simplified the numerical complexities of the SCAN functional, alleviating SCAN's stringent demand on the numerical integration grids to some extent. The regularization of rSCAN, however, results in the breaking of some constraints such as the uniform electron gas limit, the slowly varying density limit, and coordinate scaling of the iso-orbital indicator. Here, we assess the effects of regularization on the electronic, structural, vibrational, and magnetic properties of molecules by comparing the SCAN and rSCAN predictions. The properties studied include atomic energies, atomization energies, ionization potentials, electron affinities, barrier heights, infrared intensities, dissociation and reaction energies, spin moments of molecular magnets, and isomer ordering of water clusters. Our results show that rSCAN requires less dense numerical grids and gives very similar results to those of SCAN for all properties examined with the exception of atomization energies, which are worsened in rSCAN. We also examine the performance of self-interaction-corrected (SIC) rSCAN with respect to SIC-SCAN using the Perdew-Zunger (PZ) SIC method. The PZSIC method uses orbital densities to compute one-electron self-interaction errors and places an even more stringent demand on numerical grids. Our results show that SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work with numerical grids that are on average half or less as dense as that needed for SIC-SCAN.
Collapse
Affiliation(s)
- Yoh Yamamoto
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA.
| | - Alan Salcedo
- 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. and Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Md Shamsul Alam
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA. and Computational Science Program, 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. and Computational Science Program, 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. and Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
12
|
Yamamoto Y, Romero S, Baruah T, Zope RR. Improvements in the orbitalwise scaling down of Perdew–Zunger self-interaction correction in many-electron regions. J Chem Phys 2020; 152:174112. [DOI: 10.1063/5.0004738] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yoh Yamamoto
- 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
| | - Tunna Baruah
- 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
|
13
|
Schwalbe S, Trepte K, Fiedler L, Johnson AI, Kraus J, Hahn T, Peralta JE, Jackson KA, Kortus J. Interpretation and Automatic Generation of Fermi‐Orbital Descriptors. J Comput Chem 2019; 40:2843-2857. [DOI: 10.1002/jcc.26062] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Revised: 08/13/2019] [Accepted: 08/14/2019] [Indexed: 01/28/2023]
Affiliation(s)
- Sebastian Schwalbe
- TU Freiberg Institute of Theoretical Physics Leipziger Str. 23 D‐09599 Freiberg Germany
| | - Kai Trepte
- Department of Physics Central Michigan University Mount Pleasant Michigan 48859
| | - Lenz Fiedler
- TU Freiberg Institute of Theoretical Physics Leipziger Str. 23 D‐09599 Freiberg Germany
- Freiberg Instruments GmbH Delfter Str.6 D‐09599 Freiberg Germany
| | - Alex I. Johnson
- Department of Physics Central Michigan University Mount Pleasant Michigan 48859
| | - Jakob Kraus
- TU Freiberg Institute of Theoretical Physics Leipziger Str. 23 D‐09599 Freiberg Germany
| | - Torsten Hahn
- TU Freiberg Institute of Theoretical Physics Leipziger Str. 23 D‐09599 Freiberg Germany
| | - Juan E. Peralta
- Department of Physics Central Michigan University Mount Pleasant Michigan 48859
| | - Koblar A. Jackson
- Department of Physics Central Michigan University Mount Pleasant Michigan 48859
| | - Jens Kortus
- TU Freiberg Institute of Theoretical Physics Leipziger Str. 23 D‐09599 Freiberg Germany
| |
Collapse
|