1
|
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
|
2
|
Chan B. Limiting factors in the accuracy of DFT calculation for redox potentials. J Comput Chem 2024; 45:1177-1186. [PMID: 38311976 DOI: 10.1002/jcc.27320] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 01/17/2024] [Accepted: 01/18/2024] [Indexed: 02/06/2024]
Abstract
In the present study, we have investigated factors affecting the accuracy of computational chemistry calculation of redox potentials, namely the gas-phase ionization energy (IE) and electron affinity (EA), and the continuum solvation effect. In general, double-hybrid density functional theory methods yield IEs and EAs that are on average within ~0.1 eV of our high-level W3X-L benchmark, with the best performing method being DSD-BLYP/ma-def2-QZVPP. For lower-cost methods, the average errors are ~0.2-0.3 eV, with ωB97X-3c being the most accurate (~0.15 eV). For the solvation component, essentially all methods have an average error of ~0.3 eV, which shows the limitation of the continuum solvation model. Curiously, the directly calculated redox potentials show errors of ~0.3 eV for all methods. These errors are notably smaller than what can be expected from error propagation with the two components (IE and EA, and solvation effect). Such a discrepancy can be attributed to the cancellation of errors, with the lowest-cost GFN2-xTB method benefiting the most, and the most accurate ωB97X-3c method benefiting the least. For organometallic species, the redox potentials show large deviations exceeding ~0.5 eV even for DSD-BLYP. The large errors are attributed to those for the gas-phase IEs and EAs, which represents a major barrier to the accurate calculation of redox potentials for such systems.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Nagasaki, Japan
- RIKEN Center for Computational Science, Kobe, Japan
| |
Collapse
|
3
|
Franzke YJ, Bruder F, Gillhuber S, Holzer C, Weigend F. Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods. J Phys Chem A 2024; 128:670-686. [PMID: 38195394 DOI: 10.1021/acs.jpca.3c07093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2024]
Abstract
An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Florian Bruder
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Sebastian Gillhuber
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology (KIT), Engesserstraße 15, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| |
Collapse
|
4
|
Fransson T, Pettersson LGM. Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange. J Chem Theory Comput 2023; 19:7333-7342. [PMID: 37787584 PMCID: PMC10601479 DOI: 10.1021/acs.jctc.3c00728] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Indexed: 10/04/2023]
Abstract
The structure and dynamics of liquid water continue to be debated, with insight provided by, among others, X-ray emission spectroscopy (XES), which shows a split in the high-energy 1b1 feature. This split is yet to be reproduced by theory, and it remains unclear if these difficulties are related to inaccuracies in dynamics simulations, spectrum calculations, or both. We investigate the performance of different methods for calculating XES of liquid water, focusing on the ability of time-dependent density functional theory (TDDFT) to reproduce reference spectra obtained by high-level coupled cluster and algebraic-diagrammatic construction scheme calculations. A metric for evaluating the agreement between theoretical spectra termed the integrated absolute difference (IAD), which considers the integral of shifted difference spectra, is introduced and used to investigate the performance of different exchange-correlation functionals. We find that computed spectra of symmetric and asymmetric model water structures are strongly and differently influenced by the amount of Hartree-Fock exchange, with best agreement to reference spectra for ∼40-50%. Lower percentages tend to yield high density of contributing states, resulting in too broad features. The method introduced here is useful also for other spectrum calculations, in particular where the performance for ensembles of structures are evaluated.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| |
Collapse
|
5
|
Abstract
Weakly bound non-valence anions are molecular systems where the excess electron stabilizes in a very diffuse orbital whose size, shape, and binding energy (∼1-100 meV) are governed by the long-range electrostatic potential of the molecule. Its binding energy comes mainly from charge-dipole or charge-multipole interactions or dispersion forces. While highly correlated methods, like coupled cluster methods, are considered to be the state of the art for describing anionic systems, especially when the electron lies in a very diffuse orbital, we consider here the possibility to use DFT-based calculations. In such molecular anions, the outer electron experiences long-range exchange and correlation interactions. We show that DFT can describe long-range bound states provided that a correct asymptotic exchange and correlation potential is used, namely, that from a range-separated hybrid functional. This opens an alternative to the computationally demanding highly correlated method calculations. It is also suggested that the study of weakly bound anions could help in the construction of new DFT potentials to study systems where nonlocal effects are significant.
Collapse
Affiliation(s)
- Guillaume Thiam
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, UMR5306, F-69622 Villeurbanne, France
| | - Franck Rabilloud
- Université de Lyon, Université Claude Bernard Lyon 1, CNRS, Institut Lumière Matière, UMR5306, F-69622 Villeurbanne, France
| |
Collapse
|
6
|
Mendes R, da Mata V, Haiduke R. An electronic structure investigation of excited state intramolecular proton transfer in amino-benzazole derivatives: relative energies and electron density descriptors. J Photochem Photobiol A Chem 2023. [DOI: 10.1016/j.jphotochem.2023.114738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
|
7
|
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: 62] [Impact Index Per Article: 31.0] [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
|
8
|
Holzer C, Franzke YJ. A Local Hybrid Exchange Functional Approximation from First Principles. J Chem Phys 2022; 157:034108. [DOI: 10.1063/5.0100439] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Local hybrid functionals are a more flexible class of density functional approximations allowing for a position-dependent admixture of exact exchange. This additional flexibility, however, comes with a more involved mathematical form and a more complicated design. A common denominator for previously constructed local hybrid funtionals is usage of thermochemical benchmark data to construct these functionals. Herein, we design a local hybrid functional without relying on benchmark data. Instead, we construct it in a more ab initio manner, following the principles of modern meta-generalized gradient approximations and considering theoretical constrains. To achieve this, we make use of the density matrix expansion and a local mixing function based on an approximate correlation length. The accuracy of the developed density functional approximation is assessed for thermochemistry, excitation energies, polarizabilities, magnetizabilities, NMR spin-spincoupling constants, NMR shieldings and shifts, as well as EPR g-tensors and hyperfine coupling constants. Here, the new exchange functional shows a robust performance and is especially well suited for atomization energies, barrier heights, excitation energies, NMR coupling constants, and EPR properties, whereas it looses some ground for the NMR shifts.Therefore, the designed functional is a major step forwards for functionals that have been designed from first principles.
Collapse
Affiliation(s)
- Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruher Institut für Technologie Fakultät für Physik, Germany
| | - Yannick J. Franzke
- Fachbereich Chemie, Philipps-Universität Marburg Fachbereich Chemie, Germany
| |
Collapse
|
9
|
Franzke YJ, Yu JM. Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance. J Chem Theory Comput 2022; 18:2246-2266. [PMID: 35354319 DOI: 10.1021/acs.jctc.1c01175] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C5Me5)2Y(μ-S)2Mo(μ-S)2Y(C5Me5)2]-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C6Cl5)4]- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Jason M Yu
- Department of Chemistry, University of California─Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
10
|
Franzke YJ, Holzer C, Mack F. NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation. J Chem Theory Comput 2022; 18:1030-1045. [PMID: 34981925 DOI: 10.1021/acs.jctc.1c00999] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
11
|
Franzke YJ, Yu JM. Hyperfine Coupling Constants in Local Exact Two-Component Theory. J Chem Theory Comput 2021; 18:323-343. [PMID: 34928142 DOI: 10.1021/acs.jctc.1c01027] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within the one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, that is, the derivatives of the unitary decoupling transformation are considered. This requires the solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened nuclear spin-orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12-17 transition-metal compounds. The error introduced by DLU is negligible, and the DLU-X2C Hamiltonian accurately reproduces its four-component "fully" relativistic parent results. Functionals with a large amount of Hartree-Fock exchange such as CAM-QTP-02 and ωB97X-D are generally favorable. The pure density functional r2SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-ζ bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C6Cl5)4]- with more than 4700 primitive basis functions and four rare-earth single-molecule magnets: [La(OAr*)3]-, [Lu(NR2)3]-, [Lu(OAr*)3]-, and [TbPc2]-. Here, the results with the spin-orbit DLU-X2C Hamiltonian are in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| |
Collapse
|
12
|
Chan B, Dawson W, Nakajima T, Hirao K. Taking Advantage of a Systematic Energy Non-linearity Error in Density Functional Theory for the Calculation of Electronic Energy Levels. J Phys Chem A 2021; 125:10507-10513. [PMID: 34874170 DOI: 10.1021/acs.jpca.1c06723] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present an approximate approach for the calculation of ionization potential (IP) and electron affinity (EA) by exploiting the complementary energy non-linearity errors for a species M and its one-electron-ionized counterpart (M+). Reasonable IPs and EAs are thus obtained by averaging the orbital energies of M and M+, even with a low-level method such as BLYP/6-31G(d). By combining the corrected IPs and EAs, we can further obtain reasonable excitation energies. The errors in uncorrected valence IPs and uncorrected virtual-orbital energies show systematic trends. These characteristics provide a convenient and computationally efficient avenue for qualitative estimation of these properties with single corrections for multiple IPs and excitation energies.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| | - William Dawson
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Kimihiko Hirao
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan.,Fukui Institute for Fundamental Chemistry, Kyoto University, Takano, Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
| |
Collapse
|
13
|
Steffler F, Haiduke RL. Investigating the electronic excitations in Polyoxoniobates: (Nb6O19)8−, (Nb10O28)6− and (XNb12O40)Y with (X=As, P, Si, Ge) and (Y=15− and 16−). J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2021.131156] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
14
|
Manzhos S, Chueh CC, Giorgi G, Kubo T, Saianand G, Lüder J, Sonar P, Ihara M. Materials Design and Optimization for Next-Generation Solar Cell and Light-Emitting Technologies. J Phys Chem Lett 2021; 12:4638-4657. [PMID: 33974435 DOI: 10.1021/acs.jpclett.1c00714] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We review some of the most potent directions in the design of materials for next-generation solar cell and light-emitting technologies that go beyond traditional solid-state inorganic semiconductor-based devices, from both the experimental and computational standpoints. We focus on selected recent conceptual advances in tackling issues which are expected to significantly impact applied literature in the coming years. Specifically, we consider solution processability, design of dopant-free charge transport materials, two-dimensional conjugated polymeric semiconductors, and colloidal quantum dot assemblies in the fields of experimental synthesis, characterization, and device fabrication. Key modeling issues that we consider are calculations of optical properties and of effects of aggregation, including recent advances in methods beyond linear-response time-dependent density functional theory and recent insights into the effects of correlation when going beyond the single-particle ansatz as well as in the context of modeling of thermally activated fluorescence.
Collapse
Affiliation(s)
- Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| | - Chu-Chen Chueh
- Department of Chemical Engineering, National Taiwan University, Taipei 10617, Taiwan
- Advanced Research Center for Green Materials Science and Technology, National Taiwan University, Taipei 10617, Taiwan
| | - Giacomo Giorgi
- Department of Civil & Environmental Engineering (DICA), Università degli Studi di Perugia, Via G. Duranti 93, 06125 Perugia, Italy
- CNR-SCITEC, 06123 Perugia, Italy
| | - Takaya Kubo
- Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Gopalan Saianand
- School of Chemistry and Physics, Queensland University of Technology (QUT), 2 George Street, 4001 Brisbane, Australia
- Global Center for Environmental Remediation (GCER), College of Engineering, Science and Environment, The University of Newcastle, Callaghan, NSW 2308, Australia
| | - Johann Lüder
- Department of Materials and Optoelectronic Science, National Sun Yat-sen University, 80424, No. 70, Lien-Hai Road, Kaohsiung, Taiwan R.O.C
- Center of Crystal Research, National Sun Yat-sen University, 80424, No. 70, Lien-Hai Road, Kaohsiung, Taiwan R.O.C
| | - Prashant Sonar
- School of Chemistry and Physics, Queensland University of Technology (QUT), 2 George Street, 4001 Brisbane, Australia
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552, Japan
| |
Collapse
|
15
|
Holzer C, Franzke YJ, Kehry M. Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties. J Chem Theory Comput 2021; 17:2928-2947. [PMID: 33914504 DOI: 10.1021/acs.jctc.1c00203] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.
Collapse
Affiliation(s)
- Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany.,Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Max Kehry
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
16
|
Hirao K, Bae HS, Song JW, Chan B. Koopmans'-Type Theorem in Kohn-Sham Theory with Optimally Tuned Long-Range-Corrected (LC) Functionals. J Phys Chem A 2021; 125:3489-3502. [PMID: 33874719 DOI: 10.1021/acs.jpca.1c01593] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the present study, we have investigated the applicability of long-range-corrected (LC) functionals to a Kohn-Sham (KS) Koopmans'-type theorem. Specifically, we have examined the performance of optimally tuned LCgau-core functionals (in combination with BOP and PW86-PW91 exchange-correlation functionals) by calculating the ionization potential (IP) within the context of Koopmans' prediction. In the LC scheme, the electron repulsion operator, 1/r12, is divided into short-range and long-range components using a standard error function, with a range separation parameter μ determining the weight of the two ranges. For each system that we have examined (H2O, CO, benzene, N2, HF, H2CO, C2H4, and five-membered ring compounds cyclo-C4H4X, with X = CH2, NH, O, and S, and pyridine), the value of μ is optimized to minimize the deviation of the negative HOMO energy from the experimental IP. Our results demonstrate the utility of optimally tuned LC functionals in predicting the IP of outer valence levels. The accuracy is comparable to that of highly accurate ab initio theory. However, our Koopmans' method is less accurate for the inner valence and core levels. Overall, our results support the notion that orbitals in KS-DFT, when obtained with the LC functional, provide an accurate one-electron energy spectrum. This method represents a one-electron orbital theory that is attractive in its simple formulation and effective in its practical application.
Collapse
Affiliation(s)
- Kimihiko Hirao
- Fukui Institute for Fundamental Chemistry, Kyoto University, Takano, Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan.,RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Han-Seok Bae
- Department of Chemistry Education, Daegu University, Gyeongsan 113-8656, Korea
| | - Jong-Won Song
- Department of Chemistry Education, Daegu University, Gyeongsan 113-8656, Korea
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
17
|
Fransson T, Brumboiu IE, Vidal ML, Norman P, Coriani S, Dreuw A. XABOOM: An X-ray Absorption Benchmark of Organic Molecules Based on Carbon, Nitrogen, and Oxygen 1s → π* Transitions. J Chem Theory Comput 2021; 17:1618-1637. [PMID: 33544612 PMCID: PMC8023667 DOI: 10.1021/acs.jctc.0c01082] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Indexed: 01/05/2023]
Abstract
The performance of several standard and popular approaches for calculating X-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges of 40 primarily organic molecules up to the size of guanine has been evaluated, focusing on the low-energy and intense 1s → π* transitions. Using results obtained with CVS-ADC(2)-x and fc-CVS-EOM-CCSD as benchmark references, we investigate the performance of CC2, ADC(2), ADC(3/2), and commonly adopted density functional theory (DFT)-based approaches. Here, focus is on precision rather than on accuracy of transition energies and intensities-in other words, we target relative energies and intensities and the spread thereof, rather than absolute values. The use of exchange-correlation functionals tailored for time-dependent DFT calculations of core excitations leads to error spreads similar to those seen for more standard functionals, despite yielding superior absolute energies. Long-range corrected functionals are shown to perform particularly well compared to our reference data, showing error spreads in energy and intensity of 0.2-0.3 eV and ∼10%, respectively, as compared to 0.3-0.6 eV and ∼20% for a typical pure hybrid. In comparing intensities, state mixing can complicate matters, and techniques to avoid this issue are discussed. Furthermore, the influence of basis sets in high-level ab initio calculations is investigated, showing that reasonably accurate results are obtained with the use of 6-311++G**. We name this benchmark suite as XABOOM (X-ray absorption benchmark of organic molecules) and provide molecular structures and ground-state self-consistent field energies and spectroscopic data. We believe that it provides a good assessment of electronic structure theory methods for calculating X-ray absorption spectra and will become useful for future developments in this field.
Collapse
Affiliation(s)
- Thomas Fransson
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
- Fysikum, Stockholm University, Albanova, 10691 Stockholm, Sweden
| | - Iulia E. Brumboiu
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology, 34141 Daejeon, Korea
| | - Marta L. Vidal
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Patrick Norman
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Sonia Coriani
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
- Department
of Chemistry, NTNU-Norwegian University
of Science and Technology, N-7991 Trondheim, Norway
| | - Andreas Dreuw
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
18
|
Mendes RA, Haiduke RLA, Bartlett RJ. The Devil's Triangle of Kohn-Sham density functional theory and excited states. J Chem Phys 2021; 154:074106. [PMID: 33607901 DOI: 10.1063/5.0035446] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Exchange-correlation (XC) functionals from Density Functional Theory (DFT) developed under the rigorous arguments of Correlated Orbital Theory (COT) address the Devil's Triangle of prominent errors in Kohn-Sham (KS) DFT. At the foundation of this triangle lie the incorrect one-particle spectrum, the lack of integer discontinuity, and the self-interaction error. At the top level, these failures manifest themselves in incorrect charge transfer and Rydberg excitation energies, along with poor activation barriers. Accordingly, the Quantum Theory Project (QTP) XC functionals have been created to address several of the long-term issues encountered in KS theory and its Time Dependent DFT (TDDFT) variant for electronic excitations. Recognizing that COT starts with a correct one-particle spectrum, a condition imposed on the QTP functionals by means of minimum parameterization, the question that arises is how does this affect the electronically excited states? Among up to 28 XC functionals considered, the QTP family provides one of the smallest mean absolute deviations for charge-transfer excitations while also showing excellent results for Rydberg states. However, there is some room for improvement in the case of excitation energies to valence states, which are systematically underestimated by all functionals investigated. An alternative path for better treatment of excitation energies, mainly for valence states, is offered by using orbital energies from QTP functionals, especially by CAM-QTP-02 and LC-QTP. In this case, the deviations from the reference data can be reduced approximately by half.
Collapse
Affiliation(s)
- Rodrigo A Mendes
- Departamento de Química e Física Molecular, Instituto de Química de São Carlos, Universidade de São Paulo, CP 780, 13560-970 São Carlos, SP, Brazil
| | - Roberto L A Haiduke
- Departamento de Química e Física Molecular, Instituto de Química de São Carlos, Universidade de São Paulo, CP 780, 13560-970 São Carlos, SP, Brazil
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
19
|
Ghosal A, Gupta T, Mahato K, Roy AK. Excitation energies through Becke’s exciton model within a Cartesian-grid KS DFT. Theor Chem Acc 2021. [DOI: 10.1007/s00214-020-02699-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
20
|
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
|
21
|
Hirao K, Nakajima T, Chan B, Song JW, Bae HS. Core-Level Excitation Energies of Nucleic Acid Bases Expressed as Orbital Energies of the Kohn–Sham Density Functional Theory with Long-Range Corrected Functionals. J Phys Chem A 2020; 124:10482-10494. [DOI: 10.1021/acs.jpca.0c07087] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kimihiko Hirao
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Takano, Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| | - Jong-Won Song
- Department of Chemistry Education, Daegu University, Gyeongsan 113-8656, Korea
| | - Han-Seok Bae
- Department of Chemistry Education, Daegu University, Gyeongsan 113-8656, Korea
| |
Collapse
|
22
|
Hirao K, Chan B, Song JW, Bae HS. Charge-Transfer Excitation Energies Expressed as Orbital Energies of Kohn–Sham Density Functional Theory with Long-Range Corrected Functionals. J Phys Chem A 2020; 124:8079-8087. [DOI: 10.1021/acs.jpca.0c05414] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kimihiko Hirao
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Takano, Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| | - Jong-Won Song
- Department of Chemistry Education, Daegu University, Gyeongsan 113-8656, Korea
| | - Han-Seok Bae
- Department of Chemistry Education, Daegu University, Gyeongsan 113-8656, Korea
| |
Collapse
|
23
|
Hirao K, Chan B, Song J, Bhattarai K, Tewary S. Excitation energies expressed as orbital energies of Kohn–Sham density functional theory with long‐range corrected functionals. J Comput Chem 2020; 41:1368-1383. [DOI: 10.1002/jcc.26181] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Revised: 02/05/2020] [Accepted: 02/13/2020] [Indexed: 01/28/2023]
Affiliation(s)
- Kimihiko Hirao
- RIKEN Center for Computational Science Chuo‐ku Kobe Japan
- Fukui Institute for Fundamental ChemistryKyoto University Sakyo‐ku Kyoto Japan
| | - Bun Chan
- Graduate School of EngineeringNagasaki University Nagasaki Japan
| | - Jong‐Won Song
- Department of Chemistry EducationDaegu University Gyeongsan South Korea
| | - Kamala Bhattarai
- Department of Chemistry EducationDaegu University Gyeongsan South Korea
| | - Subrata Tewary
- RIKEN Center for Computational Science Chuo‐ku Kobe Japan
| |
Collapse
|
24
|
Martin JML, Santra G. Empirical Double‐Hybrid Density Functional Theory: A ‘Third Way’ in Between WFT and DFT. Isr J Chem 2019. [DOI: 10.1002/ijch.201900114] [Citation(s) in RCA: 81] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Jan M. L. Martin
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| | - Golokesh Santra
- Department of Organic Chemistry Weizmann Institute of Science 76100 Rehovot Israel
| |
Collapse
|
25
|
Dittmer A, Izsák R, Neese F, Maganas D. Accurate Band Gap Predictions of Semiconductors in the Framework of the Similarity Transformed Equation of Motion Coupled Cluster Theory. Inorg Chem 2019; 58:9303-9315. [PMID: 31240911 PMCID: PMC6750750 DOI: 10.1021/acs.inorgchem.9b00994] [Citation(s) in RCA: 40] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
In
this work, we present a detailed comparison between wave-function-based
and particle/hole techniques for the prediction of band gap energies
of semiconductors. We focus on the comparison of the back-transformed
Pair Natural Orbital Similarity Transformed Equation of Motion Coupled-Cluster
(bt-PNO-STEOM-CCSD) method with Time Dependent Density Functional
Theory (TD-DFT) and Delta Self Consistent Field/DFT (Δ-SCF/DFT)
that are employed to calculate the band gap energies in a test set
of organic and inorganic semiconductors. Throughout, we have used
cluster models for the calculations that were calibrated by comparing
the results of the cluster calculations to periodic DFT calculations
with the same functional. These calibrations were run with cluster
models of increasing size until the results agreed closely with the
periodic calculation. It is demonstrated that bt-PNO-STEOM-CC yields
accurate results that are in better than 0.2 eV agreement with the
experiment. This holds for both organic and inorganic semiconductors.
The efficiency of the employed computational protocols is thoroughly
discussed. Overall, we believe that this study is an important contribution
that can aid future developments and applications of excited state
coupled cluster methods in the field of solid-state chemistry and
heterogeneous catalysis. In this work, it is shown
that a combination of the embedded cluster approach with wave-function-based
ab initio methods in the framework of the Similarity Transformed Equation
of Motion Coupled Cluster (bt-PNO STEOM-CC) provides an accurate protocol
for band gap energy predictions in classes of organic and inorganic
semiconductors.
Collapse
Affiliation(s)
- Anneke Dittmer
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , 45470 Mülheim an der Ruhr , Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , 45470 Mülheim an der Ruhr , Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , 45470 Mülheim an der Ruhr , Germany
| | - Dimitrios Maganas
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , 45470 Mülheim an der Ruhr , Germany
| |
Collapse
|
26
|
Mei Y, Yang W. Excited-State Potential Energy Surfaces, Conical Intersections, and Analytical Gradients from Ground-State Density Functional Theory. J Phys Chem Lett 2019; 10:2538-2545. [PMID: 31038964 PMCID: PMC7449530 DOI: 10.1021/acs.jpclett.9b00712] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/09/2023]
Abstract
Kohn-Sham density functional theory (KS-DFT) has been a well-established theoretical foundation for ground-state electronic structure and has achieved great success in practical calculations. Recently, utilizing the eigenvalues from KS or generalized KS (GKS) calculations as an approximation to the quasiparticle energies, our group demonstrated a method to calculate the excitation energies from (G)KS calculation on the ground-state ( N - 1)-electron system. This method is now called QE-DFT (quasiparticle energies from DFT). In this work, we extend this QE-DFT method to describe excited-state potential energy surfaces (PESs), conical intersections, and the analytical gradients of excited-state PESs. The analytical gradients were applied to perform geometry optimization for excited states. In conjunction with several commonly used density functional approximations, QE-DFT can yield PESs in the vicinity of the equilibrium structure with accuracy similar to that from time-dependent DFT (TD-DFT). Furthermore, it describes conical intersection well, in contrast to TD-DFT. Good results for geometry optimization, especially bond length, of low-lying excitations for 14 small molecules are presented. The capability of describing excited-state PESs, conical intersections, and analytical gradients from QE-DFT and its efficiency based on just ground-state DFT calculations should be of great interest for describing photochemical and photophysical processes in complex systems.
Collapse
Affiliation(s)
- Yuncai Mei
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
- Key Laboratory of Theoretical Chemistry of Environment, School of Chemistry and Environment, South China Normal University, Guangzhou 510006, China
- Corresponding Author:
| |
Collapse
|
27
|
Mei Y, Yang W. Charge transfer excitation energies from ground state density functional theory calculations. J Chem Phys 2019; 150:144109. [PMID: 30981264 DOI: 10.1063/1.5087883] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Calculating charge transfer (CT) excitation energies with high accuracy and low computational cost is a challenging task. Kohn-Sham density functional theory (KS-DFT), due to its efficiency and accuracy, has achieved great success in describing ground state problems. To extend to excited state problems, our group recently demonstrated an approach with good numerical results to calculate low-lying and Rydberg excitation energies of an N-electron system from a ground state KS or generalized KS calculations of an (N - 1)-electron system via its orbital energies. In the present work, we explore further the same methodology to describe CT excitations. Numerical results from this work show that performance of conventional density functional approximations (DFAs) is not as good for CT excitations as for other excitations due to the delocalization error. Applying localized orbital scaling correction (LOSC) to conventional DFAs, a recently developed method in our group to effectively reduce the delocalization error, can improve the results. Overall, the performance of this methodology is better than time dependent DFT (TDDFT) with conventional DFAs. In addition, it shows that results from LOSC-DFAs in this method reach similar accuracy to other methods, such as ΔSCF, G0W0 with Bethe-Salpeter equations, particle-particle random phase approximation, and even high-level wavefunction methods like CC2. Our analysis shows that the correct 1/R trend for CT excitation can be captured from LOSC-DFA calculations, stressing that the application of DFAs with the minimal delocalization error is essential within this methodology. This work provides an efficient way to calculate CT excitation energies from ground state DFT.
Collapse
Affiliation(s)
- Yuncai Mei
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
28
|
Jin Y, Yang W. Excitation Energies from the Single-Particle Green's Function with the GW Approximation. J Phys Chem A 2019; 123:3199-3204. [PMID: 30920830 DOI: 10.1021/acs.jpca.9b02379] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quasi-particle energies are important in predicting molecular ionization energies and bulk band structures. The state-of-the-art method for quasi-particle energy calculations, particularly for bulk systems, is the GW approximation. For excited state calculations, one needs to go beyond the GW approximation. The Bethe-Salpeter equation (BSE) is the commonly used approach for bulk-system excited state calculations beyond the GW approximation, which is accurate but computationally cumbersome. In this Article, we develop a new method to extract excitation energies directly from the quasi-particle energies based on the GW approximation. Starting from the ( N - 1)-electron system, we are able to calculate molecular excitation energies with orbital energies at the GW level for HOMO excitations. Our calculations demonstrate that this method can accurately capture low-lying local excitations as well as charge transfer excitations in many molecular systems. Our method is shown to outperform the time-dependent density functional theory (TDDFT) and are comparable with higher level excited state calculations, including the equation-of-motion couple cluster (EOM-CC) theory and the BSE, but with less computational effort. This new approach provides an efficient alternative to the BSE method for accurate excited state calculations.
Collapse
Affiliation(s)
- Ye Jin
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States
| | - Weitao Yang
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States.,Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment , South China Normal University , Guangzhou 510006 , China
| |
Collapse
|
29
|
Ranasinghe DS, Margraf JT, Perera A, Bartlett RJ. Vertical valence ionization potential benchmarks from equation-of-motion coupled cluster theory and QTP functionals. J Chem Phys 2019; 150:074108. [DOI: 10.1063/1.5084728] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Duminda S. Ranasinghe
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Johannes T. Margraf
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
- Chair for Theoretical Chemistry, Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
30
|
Mei Y, Li C, Su NQ, Yang W. Approximating Quasiparticle and Excitation Energies from Ground State Generalized Kohn-Sham Calculations. J Phys Chem A 2019; 123:666-673. [PMID: 30589546 DOI: 10.1021/acs.jpca.8b10380] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quasiparticle energies and fundamental band gaps in particular are critical properties of molecules and materials. It was rigorously established that the generalized Kohn-Sham HOMO and LUMO orbital energies are the chemical potentials of electron removal and addition and thus good approximations to band edges and fundamental gaps from a density functional approximation (DFA) with minimal delocalization error. For other quasiparticle energies, their connection to the generalized Kohn-Sham orbital energies has not been established but remains highly interesting. We provide the comparison of experimental quasiparticle energies for many finite systems with calculations from the GW Green function and localized orbitals scaling correction (LOSC), a recently developed correction to semilocal DFAs, which has minimal delocalization error. Extensive results with over 40 systems clearly show that LOSC orbital energies achieve slightly better accuracy than the GW calculations with little dependence on the semilocal DFA, supporting the use of LOSC DFA orbital energies to predict quasiparticle energies. This also leads to the calculations of excitation energies of the N-electron systems from the ground state DFA calculations of the ( N - 1)-electron systems. Results show good performance with accuracy similar to TDDFT and the delta SCF approach for valence excitations with commonly used DFAs with or without LOSC. For Rydberg states, good accuracy was obtained only with the use of LOSC DFA. This work highlights the pathway to quasiparticle and excitation energies from ground density functional calculations.
Collapse
Affiliation(s)
- Yuncai Mei
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States
| | - Chen Li
- 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.,Key Laboratory of Theoretical Chemistry of Environment, School of Chemistry and Environment , South China Normal University , Guangzhou 510006 , China
| |
Collapse
|