1
|
Ohad G, Hartstein M, Gould T, Neaton JB, Kronik L. Nonempirical Prediction of the Length-Dependent Ionization Potential in Molecular Chains. J Chem Theory Comput 2024; 20. [PMID: 39137361 PMCID: PMC11360138 DOI: 10.1021/acs.jctc.4c00847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2024] [Accepted: 07/30/2024] [Indexed: 08/15/2024]
Abstract
The ionization potential of molecular chains is well-known to be a tunable nanoscale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nanosized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an ansatz that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nanoscale regime.
Collapse
Affiliation(s)
- Guy Ohad
- 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
| | - Tim Gould
- Queensland
Micro- and Nanotechnology Centre, Griffith
University, Nathan QLD 4111, Australia
| | - Jeffrey B. Neaton
- Department
of Physics, University of California, Berkeley, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kavli
Energy NanoSciences Institute at Berkeley, University of California, Berkeley, Berkeley, California 94720, United States
| | - Leeor Kronik
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
2
|
Brütting M, Bahmann H, Kümmel S. Combining Local Range Separation and Local Hybrids: A Step in the Quest for Obtaining Good Energies and Eigenvalues from One Functional. J Phys Chem A 2024; 128:5212-5223. [PMID: 38905018 DOI: 10.1021/acs.jpca.4c02787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/23/2024]
Abstract
Some of the most successful exchange-correlation approximations in density functional theory are "hybrids", i.e., they rely on combining semilocal density functionals with exact nonlocal Fock exchange. In recent years, two classes of hybrid functionals have emerged as particularly promising: range-separated hybrids on the one hand, and local hybrids on the other hand. These functionals offer the hope to overcome a long-standing "observable dilemma", i.e., the fact that density functionals typically yield either a good description of binding energies, as obtained, e.g., in global and local hybrids, or physically interpretable eigenvalues, as obtained, e.g., in optimally tuned range-separated hybrids. Obtaining both of these characteristics from one and the same functional with the same set of parameters has been a long-standing challenge. We here discuss combining the concepts of local range separation and local hybrids as part of a constraint-guided quest for functionals that overcome the observable dilemma.
Collapse
Affiliation(s)
- Moritz Brütting
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Hilke Bahmann
- Physical and Theoretical Chemistry, University of Wuppertal, 42097 Wuppertal, Germany
| | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| |
Collapse
|
3
|
Kuan KY, Yeh SH, Yang W, Hsu CP. Excited-State Charge Transfer Coupling from Quasiparticle Energy Density Functional Theory. J Phys Chem Lett 2024; 15:6126-6136. [PMID: 38830203 PMCID: PMC11181311 DOI: 10.1021/acs.jpclett.4c00850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 04/29/2024] [Accepted: 04/30/2024] [Indexed: 06/05/2024]
Abstract
The recently developed Quasiparticle Energy (QE) scheme, based on a DFT calculation with one more (or less) electron, offers a good description of excitation energies, even with charge transfer characters. In this work, QE is further extended to calculate electron transfer (ET) couplings involving two excited states. We tested it with a donor-acceptor complex, consisting of a furan and a 1,1-dicyanoethylene (DCNE), in which two low lying charge transfer and local excitation states are involved. With generalized Mülliken-Hush and fragment charge-difference schemes, couplings from the QE approach generally agree well with those obtained from TDDFT, except that QE couplings exhibit better exponential distance dependence. Couplings from half-energy gaps with an external field are also calculated and reported. Our results show that the QE scheme is robust in calculating ET couplings with greatly reduced computational time.
Collapse
Affiliation(s)
- Kai-Yuan Kuan
- Institute
of Chemistry, Academia Sinica, 128 Academia Road, Section 2, Nankang
District, Taipei 11529, Taiwan
| | - Shu-Hao Yeh
- Institute
of Chemistry, Academia Sinica, 128 Academia Road, Section 2, Nankang
District, Taipei 11529, Taiwan
- Department
of Chemistry, National Taiwan University, 1 Roosevelt Rd, Section 4, Da’an
District, Taipei City 10617, Taiwan
| | - Weitao Yang
- Department
of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Chao-Ping Hsu
- Institute
of Chemistry, Academia Sinica, 128 Academia Road, Section 2, Nankang
District, Taipei 11529, Taiwan
- Division
of Physics, National Center for Theoretical
Sciences, 1 Roosevelt
Road, Section 4, Taipei City 10617, Taiwan
| |
Collapse
|
4
|
Zope RR, Yamamoto Y, Baruah T. How well do one-electron self-interaction-correction methods perform for systems with fractional electrons? J Chem Phys 2024; 160:084102. [PMID: 38385511 DOI: 10.1063/5.0182773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2023] [Accepted: 01/28/2024] [Indexed: 02/23/2024] Open
Abstract
Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (zσ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew-Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of H2+ and He2+ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li+ to F- at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.
Collapse
Affiliation(s)
- Rajendra R Zope
- Department of Physics, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Yoh Yamamoto
- 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
| |
Collapse
|
5
|
Shi Y, Shi Y, Wasserman A. Stretching Bonds without Breaking Symmetries in Density Functional Theory. J Phys Chem Lett 2024; 15:826-833. [PMID: 38232318 DOI: 10.1021/acs.jpclett.3c03073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
Kohn-Sham density functional theory (KS-DFT) stands out among electronic structure methods due to its balance of accuracy and computational efficiency. However, to achieve chemically accurate energies, standard density functional approximations in KS-DFT often need to break underlying symmetries, a long-standing "symmetry dilemma". By employing fragment spin densities as the main variables in calculations (rather than total molecular densities, as in KS-DFT), we present an embedding framework in which this symmetry dilemma is understood and partially resolved. The spatial overlap between fragment densities is used as the main ingredient to construct a simple, physically motivated approximation to a universal functional of the fragment densities. This "overlap approximation" is shown to significantly improve semilocal KS-DFT binding energies of molecules without artificially breaking either charge or spin symmetries. The approach is shown to be applicable to covalently bonded molecules and to systems of the "strongly correlated" type.
Collapse
Affiliation(s)
- Yuming Shi
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, United States
| | - Yi Shi
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Adam Wasserman
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, United States
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
6
|
Graf D, Thom AJW. Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost. J Chem Phys 2023; 159:174106. [PMID: 37921249 DOI: 10.1063/5.0168569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| |
Collapse
|
7
|
Linscott EB, Colonna N, De Gennaro R, Nguyen NL, Borghi G, Ferretti A, Dabo I, Marzari N. koopmans: An Open-Source Package for Accurately and Efficiently Predicting Spectral Properties with Koopmans Functionals. J Chem Theory Comput 2023; 19:7097-7111. [PMID: 37610300 PMCID: PMC10601481 DOI: 10.1021/acs.jctc.3c00652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Indexed: 08/24/2023]
Abstract
Over the past decade we have developed Koopmans functionals, a computationally efficient approach for predicting spectral properties with an orbital-density-dependent functional framework. These functionals impose a generalized piecewise linearity condition to the entire electronic manifold, ensuring that orbital energies match the corresponding electron removal/addition energy differences (in contrast to semilocal DFT, where a mismatch between the two lies at the heart of the band gap problem and, more generally, the unreliability of Kohn-Sham orbital energies). This strategy has proven to be very powerful, yielding molecular orbital energies and solid-state band structures with comparable accuracy to many-body perturbation theory but at greatly reduced computational cost while preserving a functional formulation. This paper reviews the theory of Koopmans functionals, discusses the algorithms necessary for their implementation, and introduces koopmans, an open-source package that contains all of the code and workflows needed to perform Koopmans functional calculations and obtain reliable spectral properties of molecules and materials.
Collapse
Affiliation(s)
- Edward B. Linscott
- Theory
and Simulation of Materials (THEOS), École
Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Nicola Colonna
- Laboratory
for Neutron Scattering and Imaging, Paul
Scherrer Institut, 5232 Villigen, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
| | - Riccardo De Gennaro
- Theory
and Simulation of Materials (THEOS), École
Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Ngoc Linh Nguyen
- Faculty
of Materials Science and Engineering, Phenikaa
University, Hanoi 12116, Vietnam
- A&A
Green Phoenix Group JSC, Phenikaa Research
and Technology Institute (PRATI), No. 167 Hoang Ngan, Trung Hoa, Cau Giay, Hanoi 11313, Vietnam
| | - Giovanni Borghi
- Theory
and Simulation of Materials (THEOS), École
Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | | | - Ismaila Dabo
- Department
of Materials Science and Engineering, Materials Research Institute,
and Institutes of Energy and the Environment, The Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Nicola Marzari
- Theory
and Simulation of Materials (THEOS), École
Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, 1015 Lausanne, Switzerland
- Laboratory
for Materials Simulations, Paul Scherrer
Institut, 5232 Villigen, Switzerland
| |
Collapse
|
8
|
Kim JH, Kim D, Yang W, Baik MH. Fractional Charge Density Functional Theory and Its Application to the Electro-inductive Effect. J Phys Chem Lett 2023; 14:3329-3334. [PMID: 36989527 DOI: 10.1021/acs.jpclett.3c00323] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
We employed the chemical potential equalization principle to demonstrate that fractional electrons are involved in the electro-inductive effect as well as the vibrational Stark effect. By the chemical potential model, we were able to deduce that the frontier molecular orbitals of immobilized molecules can provide valuable insight into these effects. To further understand and quantify these findings, we introduced fractional charge density functional theory (FC-DFT), a canonical ensemble approach for open systems. This method allows for the calculation of electronic energies, nuclear gradients, and the Hessian matrix of fractional electronic systems. To correct the spurious delocalization error commonly found in approximate density functionals for small systems, we imposed the Perdew-Parr-Levy-Balduz (PPLB) condition through linear interpolation of two adjacent integer points (LI-FC-DFT). Although this approach is relatively simple in terms of molecular modeling, the results obtained through LI-FC-DFT calculations predict the same trend seen in experimental reactivity and the frequency change of immobilized molecules.
Collapse
Affiliation(s)
- Jun-Hyeong Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
- Center for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
| | - Dongju Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
- Center for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Mu-Hyun Baik
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
- Center for Catalytic Hydrocarbon Functionalizations, Institute for Basic Science (IBS), Daejeon 34141, Republic of Korea
| |
Collapse
|
9
|
Teale AM, Helgaker T, Savin A, Adamo C, Aradi B, Arbuznikov AV, Ayers PW, Baerends EJ, Barone V, Calaminici P, Cancès E, Carter EA, Chattaraj PK, Chermette H, Ciofini I, Crawford TD, De Proft F, Dobson JF, Draxl C, Frauenheim T, Fromager E, Fuentealba P, Gagliardi L, Galli G, Gao J, Geerlings P, Gidopoulos N, Gill PMW, Gori-Giorgi P, Görling A, Gould T, Grimme S, Gritsenko O, Jensen HJA, Johnson ER, Jones RO, Kaupp M, Köster AM, Kronik L, Krylov AI, Kvaal S, Laestadius A, Levy M, Lewin M, Liu S, Loos PF, Maitra NT, Neese F, Perdew JP, Pernal K, Pernot P, Piecuch P, Rebolini E, Reining L, Romaniello P, Ruzsinszky A, Salahub DR, Scheffler M, Schwerdtfeger P, Staroverov VN, Sun J, Tellgren E, Tozer DJ, Trickey SB, Ullrich CA, Vela A, Vignale G, Wesolowski TA, Xu X, Yang W. DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science. Phys Chem Chem Phys 2022; 24:28700-28781. [PMID: 36269074 PMCID: PMC9728646 DOI: 10.1039/d2cp02827a] [Citation(s) in RCA: 73] [Impact Index Per Article: 36.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 08/09/2022] [Indexed: 12/13/2022]
Abstract
In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.
Collapse
Affiliation(s)
- Andrew M. Teale
- School of Chemistry, University of Nottingham, University ParkNottinghamNG7 2RDUK
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, 4 Place Jussieu, CEDEX 05, 75252 Paris, France.
| | - Carlo Adamo
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany.
| | - Alexei V. Arbuznikov
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7Straße des 17. Juni 13510623Berlin
| | | | - Evert Jan Baerends
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56125 Pisa, Italy.
| | - Patrizia Calaminici
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Eric Cancès
- CERMICS, Ecole des Ponts and Inria Paris, 6 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France.
| | - Emily A. Carter
- Department of Mechanical and Aerospace Engineering and the Andlinger Center for Energy and the Environment, Princeton UniversityPrincetonNJ 08544-5263USA
| | | | - Henry Chermette
- Institut Sciences Analytiques, Université Claude Bernard Lyon1, CNRS UMR 5280, 69622 Villeurbanne, France.
| | - Ilaria Ciofini
- PSL University, CNRS, ChimieParisTech-PSL, Institute of Chemistry for Health and Life Sciences, i-CLeHS, 11 rue P. et M. Curie, 75005 Paris, France.
| | - T. Daniel Crawford
- Department of Chemistry, Virginia TechBlacksburgVA 24061USA,Molecular Sciences Software InstituteBlacksburgVA 24060USA
| | - Frank De Proft
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | | | - Claudia Draxl
- Institut für Physik and IRIS Adlershof, Humboldt-Universität zu Berlin, 12489 Berlin, Germany. .,Fritz-Haber-Institut der Max-Planck-Gesellschaft, 14195 Berlin, Germany
| | - Thomas Frauenheim
- Bremen Center for Computational Materials Science, University of Bremen, P.O. Box 330440, D-28334 Bremen, Germany. .,Beijing Computational Science Research Center (CSRC), 100193 Beijing, China.,Shenzhen JL Computational Science and Applied Research Institute, 518110 Shenzhen, China
| | - Emmanuel Fromager
- Laboratoire de Chimie Quantique, Institut de Chimie, CNRS/Université de Strasbourg, 4 rue Blaise Pascal, 67000 Strasbourg, France.
| | - Patricio Fuentealba
- Departamento de Física, Facultad de Ciencias, Universidad de Chile, Casilla 653, Santiago, Chile.
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, The James Franck Institute, and Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, USA.
| | - Giulia Galli
- Pritzker School of Molecular Engineering and Department of Chemistry, The University of Chicago, Chicago, IL, USA.
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China. .,Department of Chemistry, University of Minnesota, Minneapolis, MN 55455, USA
| | - Paul Geerlings
- Research Group of General Chemistry (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium.
| | - Nikitas Gidopoulos
- Department of Physics, Durham University, South Road, Durham DH1 3LE, UK.
| | - Peter M. W. Gill
- School of Chemistry, University of SydneyCamperdown NSW 2006Australia
| | - Paola Gori-Giorgi
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Andreas Görling
- Chair of Theoretical Chemistry, University of Erlangen-Nuremberg, Egerlandstrasse 3, 91058 Erlangen, Germany.
| | - Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Gold Coast, Qld 4222, Australia.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115 Bonn, Germany.
| | - Oleg Gritsenko
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands.
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark.
| | - Erin R. Johnson
- Department of Chemistry, Dalhousie UniversityHalifaxNova ScotiaB3H 4R2Canada
| | - Robert O. Jones
- Peter Grünberg Institut PGI-1, Forschungszentrum Jülich52425 JülichGermany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, 10623, Berlin.
| | - Andreas M. Köster
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav)CDMX07360Mexico
| | - Leeor Kronik
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth, 76100, Israel.
| | - Anna I. Krylov
- Department of Chemistry, University of Southern CaliforniaLos AngelesCalifornia 90089USA
| | - Simen Kvaal
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Andre Laestadius
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - Mel Levy
- Department of Chemistry, Tulane University, New Orleans, Louisiana, 70118, USA.
| | - Mathieu Lewin
- CNRS & CEREMADE, Université Paris-Dauphine, PSL Research University, Place de Lattre de Tassigny, 75016 Paris, France.
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC 27599-3420, USA. .,Department of Chemistry, University of North Carolina, Chapel Hill, NC 27599-3290, USA
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, France.
| | - Neepa T. Maitra
- Department of Physics, Rutgers University at Newark101 Warren StreetNewarkNJ 07102USA
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser Wilhelm Platz 1, D-45470 Mülheim an der Ruhr, Germany.
| | - John P. Perdew
- Departments of Physics and Chemistry, Temple UniversityPhiladelphiaPA 19122USA
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.
| | - Pascal Pernot
- Institut de Chimie Physique, UMR8000, CNRS and Université Paris-Saclay, Bât. 349, Campus d'Orsay, 91405 Orsay, France.
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA. .,Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| | - Elisa Rebolini
- Institut Laue Langevin, 71 avenue des Martyrs, 38000 Grenoble, France.
| | - Lucia Reining
- Laboratoire des Solides Irradiés, CNRS, CEA/DRF/IRAMIS, École Polytechnique, Institut Polytechnique de Paris, F-91120 Palaiseau, France. .,European Theoretical Spectroscopy Facility
| | - Pina Romaniello
- Laboratoire de Physique Théorique (UMR 5152), Université de Toulouse, CNRS, UPS, France.
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA.
| | - Dennis R. Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS – Centre for Molecular Simulation, IQST – Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary2500 University Drive NWCalgaryAlbertaT2N 1N4Canada
| | - Matthias Scheffler
- The NOMAD Laboratory at the FHI of the Max-Planck-Gesellschaft and IRIS-Adlershof of the Humboldt-Universität zu Berlin, Faradayweg 4-6, D-14195, Germany.
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand.
| | - Viktor N. Staroverov
- Department of Chemistry, The University of Western OntarioLondonOntario N6A 5B7Canada
| | - Jianwei Sun
- Department of Physics and Engineering Physics, Tulane University, New Orleans, LA 70118, USA.
| | - Erik Tellgren
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway.
| | - David J. Tozer
- Department of Chemistry, Durham UniversitySouth RoadDurhamDH1 3LEUK
| | - Samuel B. Trickey
- Quantum Theory Project, Deptartment of Physics, University of FloridaGainesvilleFL 32611USA
| | - Carsten A. Ullrich
- Department of Physics and Astronomy, University of MissouriColumbiaMO 65211USA
| | - Alberto Vela
- Departamento de Química, Centro de Investigación y de Estudios Avanzados (Cinvestav), CDMX, 07360, Mexico.
| | - Giovanni Vignale
- Department of Physics, University of Missouri, Columbia, MO 65203, USA.
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, Université de Genève30 Quai Ernest-Ansermet1211 GenèveSwitzerland
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry, Fudan University, Shanghai 200433, China.
| | - Weitao Yang
- Department of Chemistry and Physics, Duke University, Durham, NC 27516, USA.
| |
Collapse
|
10
|
Lemke Y, Graf D, Kussmann J, Ochsenfeld C. An assessment of orbital energy corrections for the direct random phase approximation and explicit σ-functionals. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2098862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Yannick Lemke
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - Jörg Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
- Max Planck Institute for Solid State Research, Stuttgart, Germany
| |
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
|
Li Y, Li C. Exact Analytical Form of Diatomic Molecular Orbitals. ACS OMEGA 2022; 7:22594-22600. [PMID: 35811915 PMCID: PMC9260940 DOI: 10.1021/acsomega.2c01905] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Accepted: 05/24/2022] [Indexed: 06/15/2023]
Abstract
We provide the exact analytical form of diatomic molecular orbitals, as given by the solutions of a single-electron diatomic molecule with arbitrary nuclear charges, using our recently developed method for solving Schrödinger equations. We claim that the best representation of the wave function is a factorized form including a power prefactor, an exponentially decaying term, a modulator function on the exponential, and additional factors accounting for nodal surfaces and the magnetic quantum number. Applying our method, we have identified unexpected extreme points along the potential energy curves, hence revealing the limitations of the well-known concepts of bonding and antibonding.
Collapse
Affiliation(s)
- Yunzhi Li
- Beijing
National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Chen Li
- Beijing
National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
13
|
Bryenton KR, Adeleke AA, Dale SG, Johnson ER. Delocalization error: The greatest outstanding challenge in density‐functional theory. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1631] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Kyle R. Bryenton
- Department of Physics and Atmospheric Science Dalhousie University Halifax Nova Scotia Canada
| | | | - Stephen G. Dale
- Queensland Micro‐ and Nanotechnology Centre Griffith University Nathan Queensland Australia
| | - Erin R. Johnson
- Department of Physics and Atmospheric Science Dalhousie University Halifax Nova Scotia Canada
- Department of Chemistry Dalhousie University Halifax Nova Scotia Canada
| |
Collapse
|
14
|
Yu J, Su NQ, Yang W. Describing Chemical Reactivity with Frontier Molecular Orbitalets. JACS AU 2022; 2:1383-1394. [PMID: 35783161 PMCID: PMC9241161 DOI: 10.1021/jacsau.2c00085] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Revised: 04/07/2022] [Accepted: 04/19/2022] [Indexed: 06/15/2023]
Abstract
Locality in physical space is critical in understanding chemical reactivity in the analysis of various phenomena and processes in chemistry, biology, and materials science, as exemplified in the concepts of reactive functional groups and active sites. Frontier molecular orbitals (FMOs) pinpoint the locality of chemical bonds that are chemically reactive because of the associated orbital energies and thus have achieved great success in describing chemical reactivity, mainly for small systems. For large systems, however, the delocalization nature of canonical molecular orbitals makes it difficult for FMOs to highlight the locality of the chemical reactivity. To obtain localized molecular orbitals that also reflect the frontier nature of the chemical processes, we develop the concept of frontier molecular orbitalets (FMOLs) for describing the reactivity of large systems. The concept of orbitalets was developed recently in the localized orbital scaling correction method, which aims for eliminating the delocalization error in common density functional approximations. Orbitalets are localized in both physical and energy spaces and thus contain both orbital locality and energy information. The FMOLs are thus the orbitalets with energies highest among occupied orbitalets and lowest among unoccupied ones. The applications of FMOLs to hexadeca-1,3,5,7,9,11,13,15-octaene in its equilibrium geometry, inter- and intra-molecular charge-transfer systems, and two transition states of a bifurcating reaction demonstrate that FMOLs can connect quantum mechanical treatments of chemical systems and chemical reactivities by locating the reactive region of large chemical systems. Therefore, FMOLs extend the role of FMOs for small systems and describe the chemical reactivity of large systems with energy and locality insight, with potentially broad applications.
Collapse
Affiliation(s)
- Jincheng Yu
- Department
of Chemistry, Duke University, Durham, North Carolina 27708, United States
- Taishan
College, Shandong University, Jinan 250100, China
| | - Neil Qiang Su
- Department
of Chemistry, Duke University, Durham, North Carolina 27708, United States
- 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 Chemistry, Duke University, Durham, North Carolina 27708, United States
- Department
of Physics, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
15
|
Qu X, Xu P, Jiang H, He L, Ren X. DFT+U within the framework of linear combination of numerical atomic orbitals. J Chem Phys 2022; 156:234104. [PMID: 35732533 DOI: 10.1063/5.0090122] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
We present a formulation and implementation of the density functional theory (DFT)+U method within the framework of linear combination of numerical atomic orbitals (NAO). Our implementation not only enables single-point total energy and electronic-structure calculations but also provides access to atomic forces and cell stresses, hence allowing for full structure relaxations of periodic systems. Furthermore, our implementation allows one to deal with non-collinear spin texture, with the spin-orbit coupling (SOC) effect treated self-consistently. The key aspect behind our implementation is a suitable definition of the correlated subspace when multiple atomic orbitals with the same angular momentum are used, and this is addressed via the "Mulliken charge projector" constructed in terms of the first (most localized) atomic orbital within the d/f angular momentum channel. The important Hubbard U and Hund J parameters can be estimated from a screened Coulomb potential of the Yukawa type, with the screening parameter either chosen semi-empirically or determined from the Thomas-Fermi screening model. Benchmark calculations are performed for four late transition metal monoxide bulk systems, i.e., MnO, FeO, CoO, and NiO, and for the 5d-electron compounds IrO2. For the former type of systems, we check the performance of our DFT+U implementation for calculating bandgaps, magnetic moments, electronic band structures, as well as forces and stresses; for the latter, the efficacy of our DFT+U+SOC implementation is assessed. Systematic comparisons with available experimental results, especially with the results from other implementation schemes, are carried out, which demonstrate the validity of our NAO-based DFT+U formalism and implementation.
Collapse
Affiliation(s)
- Xin Qu
- Rocket Force University of Engineering, Xi'an, 710025, Shaanxi, China
| | - Peng Xu
- Rocket Force University of Engineering, Xi'an, 710025, Shaanxi, China
| | - Hong Jiang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Lixin He
- CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei, 230026, Anhui, China
| | - Xinguo Ren
- Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
| |
Collapse
|
16
|
Ahart CS, Rosso KM, Blumberger J. Implementation and Validation of Constrained Density Functional Theory Forces in the CP2K Package. J Chem Theory Comput 2022; 18:4438-4446. [PMID: 35700315 PMCID: PMC9281399 DOI: 10.1021/acs.jctc.2c00284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Constrained density
functional theory (CDFT) is a powerful tool
for the prediction of electron transfer parameters in condensed phase
simulations at a reasonable computational cost. In this work we present
an extension to CDFT in the popular mixed Gaussian/plane wave electronic
structure package CP2K, implementing the additional force terms arising
from a constraint based on Hirshfeld charge partitioning. This improves
upon the existing Becke partitioning scheme, which is prone to give
unphysical atomic charges. We verify this implementation for a variety
of systems: electron transfer in (H2O)2+ in a vacuum, electron tunnelling
between oxygen vacancy centers in solid MgO, and electron self-exchange
in aqueous Ru2+–Ru3+. We find good agreement
with previous plane-wave CDFT results for the same systems, but at
a significantly lower computational cost, and we discuss the general
reliability of condensed phase CDFT calculations.
Collapse
Affiliation(s)
- Christian S Ahart
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| | - Kevin M Rosso
- Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Jochen Blumberger
- Department of Physics and Astronomy and Thomas Young Centre, University College London, London WC1E 6BT, United Kingdom
| |
Collapse
|
17
|
Ariyarathna IR, Duan C, Kulik HJ. Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations. J Chem Phys 2022; 156:184113. [PMID: 35568536 DOI: 10.1063/5.0090128] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf-O or Hf-B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X4Σ-) and HfO(X1Σ+) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (De) from HfO to HfPo and the singlet-triplet energy gap (ΔES-T) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange-correlation functionals spanning multiple rungs of "Jacob's ladder." When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.
Collapse
Affiliation(s)
- Isuru R Ariyarathna
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Chenru Duan
- 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
|
18
|
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
|
19
|
Li J, Jin Y, Su NQ, Yang W. Combining localized orbital scaling correction and Bethe-Salpeter equation for accurate excitation energies. J Chem Phys 2022; 156:154101. [PMID: 35459294 PMCID: PMC9033305 DOI: 10.1063/5.0087498] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2022] [Accepted: 03/29/2022] [Indexed: 11/14/2022] Open
Abstract
We applied localized orbital scaling correction (LOSC) in Bethe-Salpeter equation (BSE) to predict accurate excitation energies for molecules. LOSC systematically eliminates the delocalization error in the density functional approximation and is capable of approximating quasiparticle (QP) energies with accuracy similar to or better than GW Green's function approach and with much less computational cost. The QP energies from LOSC, instead of commonly used G0W0 and evGW, are directly used in BSE. We show that the BSE/LOSC approach greatly outperforms the commonly used BSE/G0W0 approach for predicting excitations with different characters. For the calculations of Truhlar-Gagliardi test set containing valence, charge transfer, and Rydberg excitations, BSE/LOSC with the Tamm-Dancoff approximation provides a comparable accuracy to time-dependent density functional theory (TDDFT) and BSE/evGW. For the calculations of Stein CT test set and Rydberg excitations of atoms, BSE/LOSC considerably outperforms both BSE/G0W0 and TDDFT approaches with a reduced starting point dependence. BSE/LOSC is, thus, a promising and efficient approach to calculate excitation energies for molecular systems.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Ye Jin
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Neil Qiang Su
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
20
|
Mei Y, Yu J, Chen Z, Su NQ, Yang W. LibSC: Library for Scaling Correction Methods in Density Functional Theory. J Chem Theory Comput 2022; 18:840-850. [DOI: 10.1021/acs.jctc.1c01058] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Yuncai Mei
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Jincheng Yu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Zehua Chen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Neil Qiang Su
- 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
|
21
|
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
|
22
|
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
|
23
|
Kocák J, Kraisler E, Schild A. Charge-Transfer Steps in Density Functional Theory from the Perspective of the Exact Electron Factorization. J Phys Chem Lett 2021; 12:3204-3209. [PMID: 33761257 DOI: 10.1021/acs.jpclett.1c00467] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
When a molecule dissociates, the exact Kohn-Sham (KS) and Pauli potentials may form step structures. Reproducing these steps correctly is central for the description of dissociation and charge-transfer processes in density functional theory (DFT): The steps align the KS eigenvalues of the dissociating subsystems relative to each other and determine where electrons localize. While the step height can be calculated from the asymptotic behavior of the KS orbitals, this provides limited insight into what causes the steps. We give an explanation of the steps with an exact mapping of the many-electron problem to a one-electron problem, the exact electron factorization (EEF). The potentials appearing in the EEF have a clear physical meaning that translates to the DFT potentials by replacing the interacting many-electron system with the KS system. With a simple model of a diatomic, we illustrate that the steps are a consequence of spatial electron entanglement and are the result of a charge transfer. From this mechanism, the step height can immediately be deduced. Moreover, two methods to approximately reproduce the potentials during dissociation are proposed. One is based on the states of the dissociated system, while the other one is based on an analogy to the Born-Oppenheimer treatment of a molecule. The latter method also shows that the steps connect adiabatic potential energy surfaces. The view of DFT from the EEF thus provides a better understanding of how many-electron effects are encoded in a one-electron theory and how they can be modeled.
Collapse
Affiliation(s)
- Jakub Kocák
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Eli Kraisler
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | - Axel Schild
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
24
|
Mei Y, Yang N, Yang W. Describing polymer polarizability with localized orbital scaling correction in density functional theory. J Chem Phys 2021; 154:054302. [PMID: 33557560 DOI: 10.1063/5.0035883] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Polarizability reflects the response of the molecular charge distribution to an applied external electric field and thus closely relates to the molecular electron density. For the calculation of polarizability within density functional theory (DFT), it is well known that conventional density functional approximations (DFAs) greatly overestimate the results for polymers with long chains and the π-conjugated system. This is a manifestation of the delocalization error of the commonly used DFAs-they normally produce too delocalized electron density and underestimate the total energy for systems with fractional charge character, which occurs for long molecules in a longitudinal electric field. Thus, to achieve an accurate description of polarizabilities for polymeric molecular systems from DFT, applying DFAs with minimal delocalization error is very important. In this work, we use the recently developed localized orbital scaling correction (LOSC) to the conventional DFAs, which has been shown to largely eliminate the delocalization error, to calculate and study the polarizabilities of three classic polymers, polyyne, polyacetylene, and hydrogen chain. The results from this work demonstrate that applying LOSC to conventional DFAs with self-consistent field calculations can largely improve the description of polarizability from DFT calculations and the improved quality of electron density in LOSC leads to the improved results of polarizability of the polymers. However, the improvement is not complete and adjustment of the parameters in the LOSC method can further improve the accuracy to reach the level similar to the MP2 method. This work also points to the direction for the further development of LOSC in self-consistent calculations.
Collapse
Affiliation(s)
- Yuncai Mei
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Nathan Yang
- North Carolina School of Science and Mathematics, 1219 Broad St., Durham, North Carolina 27705, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
25
|
Yang X, Zheng X, Yang W. Density Functional Prediction of Quasiparticle, Excitation, and Resonance Energies of Molecules With a Global Scaling Correction Approach. Front Chem 2020; 8:588808. [PMID: 33425848 PMCID: PMC7793789 DOI: 10.3389/fchem.2020.588808] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Accepted: 09/23/2020] [Indexed: 11/13/2022] Open
Abstract
Molecular quasiparticle and excitation energies determine essentially the spectral characteristics measured in various spectroscopic experiments. Accurate prediction of these energies has been rather challenging for ground-state density functional methods, because the commonly adopted density function approximations suffer from delocalization error. In this work, by presuming a quantitative correspondence between the quasiparticle energies and the generalized Kohn–Sham orbital energies, and employing a previously developed global scaling correction approach, we achieve substantially improved prediction of molecular quasiparticle and excitation energies. In addition, we also extend our previous study on temporary anions in resonant states, which are associated with negative molecular electron affinities. The proposed approach does not require any explicit self-consistent field calculation on the excited-state species, and is thus highly efficient and convenient for practical purposes.
Collapse
Affiliation(s)
- Xiaolong Yang
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, China
| | - Xiao Zheng
- Hefei National Laboratory for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, China
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, NC, United States.,Key Laboratory of Theoretical Chemistry of Environment, School of Chemistry and Environment, South China Normal University, Guangzhou, China
| |
Collapse
|
26
|
Jin Y, Su NQ, Chen Z, Yang W. Introductory lecture: when the density of the noninteracting reference system is not the density of the physical system in density functional theory. Faraday Discuss 2020; 224:9-26. [PMID: 33084699 PMCID: PMC7746600 DOI: 10.1039/d0fd00102c] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A major challenge in density functional theory (DFT) is the development of density functional approximations (DFAs) to overcome errors in existing DFAs, leading to more complex functionals. For such functionals, we consider roles of the noninteracting reference systems. The electron density of the Kohn-Sham (KS) reference with a local potential has been traditionally defined as being equal to the electron density of the physical system. This key idea has been applied in two ways: the inverse calculation of such a local KS potential for the reference from a given density and the direct calculation of density and energy based on given DFAs. By construction, the inverse calculation can yield a KS reference with the density equal to the input density of the physical system. In application of DFT, however, it is the direct calculation of density and energy from a DFA that plays a central role. For direct calculations, we find that the self-consistent density of the KS reference defined by the optimized effective potential (OEP), is not the density of the physical system, when the DFA is dependent on the external potential. This inequality holds also for the density of generalized KS (GKS) or generalized OEP reference, which allows a nonlocal potential, when the DFA is dependent on the external potential. Instead, the density of the physical system, consistent with a given DFA, is given by the linear response of the total energy with respect to the variation of the external potential. This is a paradigm shift in DFT on the use of noninteracting references: the noninteracting KS or GKS references represent the explicit computational variables for energy minimization, but not the density of the physical system for external potential-dependent DFAs. We develop the expressions for the electron density so defined through the linear response for general DFAs, demonstrate the results for orbital functionals and for many-body perturbation theory within the second-order and the random-phase approximation, and explore the connections to developments in DFT.
Collapse
Affiliation(s)
- Ye Jin
- Department of Chemistry, Duke University, Durham, NC 27708, USA.
| | - Neil Qiang Su
- Department of Chemistry, Duke University, Durham, NC 27708, USA.
| | - Zehua Chen
- Department of Chemistry, Duke University, Durham, NC 27708, USA.
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, NC 27708, USA. and Department of Physics, Duke University, Durham, NC 27708, USA
| |
Collapse
|
27
|
Mei Y, Chen Z, Yang W. Self-Consistent Calculation of the Localized Orbital Scaling Correction for Correct Electron Densities and Energy-Level Alignments in Density Functional Theory. J Phys Chem Lett 2020; 11:10269-10277. [PMID: 33215928 PMCID: PMC8341159 DOI: 10.1021/acs.jpclett.0c03133] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The recently developed localized orbital scaling correction (LOSC) method shows the ability to systematically and size-consistently reduce the delocalization error existing in conventional density functional approximations (DFAs). However, the application of LOSC to DFAs was mainly through a post self-consistent field (SCF) manner, and few results from applying LOSC to DFAs in an SCF manner have been reported. The reason is that the originally proposed SCF approach to SCF-LOSC calculation uses an approximate Hamiltonian and encounters convergence problems easily in practice. In this work, we develop a new SCF approach with a correct Hamiltonian and achieve reliable SCF-LOSC calculations. We demonstrate the capability of the new SCF approach for SCF-LOSC to correctly describe the electron densities, total energies, and energy-level alignment for the molecular dissociation process, while conventional DFAs or LOSC-DFAs with post-SCF calculations show large errors. This work demonstrates that the new SCF approach for SCF-LOSC would be a promising method for studying problems for correct electron densities and energy-level alignments in large systems.
Collapse
Affiliation(s)
| | | | - Weitao Yang
- Key Laboratory of Theoretical Chemistry of Environment, School of Chemistry and Environment, South China Normal University, Guangzhou 510006, China
| |
Collapse
|