1
|
Alikhani ME, Madebène B, Silvi B. Microsolvation of cobalt, nickel, and copper atoms with ammonia: a theoretical study of the solvated electron precursors. J Mol Model 2024; 30:220. [PMID: 38902588 DOI: 10.1007/s00894-024-06019-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Accepted: 06/10/2024] [Indexed: 06/22/2024]
Abstract
CONTEXT The s-block metals dissolved in ammonia form metal-ammonia complexes with diffuse electrons which could be used for redox catalysis. In this theoretical paper, we investigated the possibility of the d-bloc transition metals (Mn, Fe, Co, Ni, and Cu) solvated by ammonia. It has been demonstrated that both Mn and Fe atoms undergo into an oxidative reaction with NH3 forming an inserted species, HMNH2. On the contrary, the Co, Ni, and Cu atoms can accommodate four NH3, via the coordination bond, to form the first solvation sphere within C2v, D2d, and Td point groups, respectively. Addition of a fifth NH3 constitute the second solvation shell by forming hydrogen bond with the other NH3s. Interestingly, M(NH3)4 (M = Co, Ni, and Cu) is a so-called solvated electron precursor and should be considered as a monocation M(NH3)4+ kernel in tight contact with one electron distributed over its periphery. This nearly free electron could be used to capture a CO2 molecule and engages in a reduction reaction. METHODS Geometry optimization of the stationary points on the potential energy surface was performed using density functional theory - CAM-B3LYP functional including the GD3BJ dispersion contribution - in combination with the 6-311 + + G(2d, 2p) basis set for all the atoms. All first-principles calculations were performed using the Gaussian 09 quantum chemical packages. The natural electron configuration of transition atom engaged in the compounds has been found using the natural bond orbital (NBO) method. We used the EDR (electron delocalization range) approach to analyze the structure of solvated electrons in real space. We also used the electron localization function (ELF) to measure the degree of electronic localization within a chemical compound. The EDR and ELF analyses are done using the TopMod and Multiwfn packages, respectively.
Collapse
Affiliation(s)
| | - Bruno Madebène
- Sorbonne Université CNRS, MONARIS, UMR8233, F-75005, Paris, France
| | - Bernard Silvi
- Sorbonne Université CNRS, LCT, UMR7616, F-75005, Paris, France
| |
Collapse
|
2
|
Manukovsky N, Kamieniarz G, Kronik L. Spin state and magnetic coupling in polynuclear Ni(II) complexes from density functional theory: is there an optimal amount of Fock exchange? J Chem Phys 2023; 159:154103. [PMID: 37846951 DOI: 10.1063/5.0169105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2023] [Accepted: 09/25/2023] [Indexed: 10/18/2023] Open
Abstract
Reliable prediction of the ground-state spin and magnetic coupling constants in transition-metal complexes is a well-known challenge for density functional theory (DFT). One popular strategy for addressing this long-standing issue involves the modification of the fraction of Fock exchange in a hybrid functional. Here we explore the viability of this approach using three polynuclear metal-organic complexes based on a Ni4O4 cubane motif, having different ground state spin values (S = 0, 2, 4) owing to the use of different ligands. We systematically search for an optimum fraction of Fock exchange, across various global, range-separated, and double hybrid functionals. We find that for all functionals tested, at best there only exists a very narrow range of Fock exchange fractions which results in a correct prediction of the ground-state spin for all three complexes. The useful range is functional dependent, but general trends can be identified. Typically, at least two similar systems must be used in order to determine both an upper and lower limit of the optimal range. This is likely owing to conflicting demands of minimizing delocalization errors, which typically requires a higher percentage of Fock exchange, and addressing static correlation, which typically requires a lower one. Furthermore, we find that within the optimal range of Fock exchange, the sign and relative magnitude of Ni-Ni magnetic coupling constants are reasonably well reproduced, but there is still room for quantitative improvement in the prediction. Thus, the prediction of spin state and magnetic coupling in polynuclear complexes remains an ongoing challenge for DFT.
Collapse
Affiliation(s)
- Nurit Manukovsky
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 7610001, Israel
| | | | - Leeor Kronik
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 7610001, Israel
| |
Collapse
|
3
|
Terayama K, Osaki Y, Fujita T, Tamura R, Naito M, Tsuda K, Matsui T, Sumita M. Koopmans' Theorem-Compliant Long-Range Corrected (KTLC) Density Functional Mediated by Black-Box Optimization and Data-Driven Prediction for Organic Molecules. J Chem Theory Comput 2023; 19:6770-6781. [PMID: 37729470 DOI: 10.1021/acs.jctc.3c00764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/22/2023]
Abstract
Density functional theory (DFT) is a significant computational tool that has substantially influenced chemistry, physics, and materials science. DFT necessitates parametrized approximation for determining an expected value. Hence, to predict the properties of a given molecule using DFT, appropriate parameters of the functional should be set for each molecule. Herein, we optimize the parameters of range-separated functionals (LC-BLYP and CAM-B3LYP) via Bayesian optimization (BO) to satisfy Koopmans' theorem. Our results demonstrate the effectiveness of the BO in optimizing functional parameters. Particularly, Koopmans' theorem-compliant LC-BLYP (KTLC-BLYP) shows results comparable to the experimental UV-absorption values. Furthermore, we prepared an optimized parameter dataset of KTLC-BLYP for over 3000 molecules through BO for satisfying Koopmans' theorem. We have developed a machine learning model on this dataset to predict the parameters of the LC-BLYP functional for a given molecule. The prediction model automatically predicts the appropriate parameters for a given molecule and calculates the corresponding values. The approach in this paper would be useful to develop new functionals and to update the previously developed functionals.
Collapse
Affiliation(s)
- Kei Terayama
- Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro-cho, Tsurumi-ku Kanagawa 230-0045, Japan
- Center for Advanced Intelligence Project, RIKEN, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan
- MDX Research Center for Element Strategy, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8501, Japan
| | - Yamato Osaki
- Department of Chemistry, Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Takehiro Fujita
- Research Center for Macromolecules and Biomaterials, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan
| | - Ryo Tamura
- Center for Advanced Intelligence Project, RIKEN, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan
- Center for Basic Research on Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Masanobu Naito
- Research Center for Macromolecules and Biomaterials, National Institute for Materials Science, 1-2-1 Sengen, Tsukuba, Ibaraki 305-0047, Japan
| | - Koji Tsuda
- Center for Advanced Intelligence Project, RIKEN, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan
- Center for Basic Research on Materials, National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki 305-0044, Japan
- Graduate School of Frontier Sciences, The University of Tokyo, 5-1-5 Kashiwa-no-ha, Kashiwa, Chiba 277-8561, Japan
| | - Toru Matsui
- Department of Chemistry, Graduate School of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8577, Japan
| | - Masato Sumita
- Center for Advanced Intelligence Project, RIKEN, 1-4-1 Nihonbashi, Chuo-ku, Tokyo 103-0027, Japan
| |
Collapse
|
4
|
Cui H, Yang S, Xue KH, Huang J, Miao X. On the self-consistency of DFT-1/2. J Chem Phys 2023; 158:094103. [PMID: 36889979 DOI: 10.1063/5.0137534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
Density functional theory (DFT)-1/2 is an efficient bandgap rectification method for DFT under local density approximation (LDA) or generalized gradient approximation. It was suggested that non-self-consistent DFT-1/2 should be used for highly ionic insulators like LiF, whereas self-consistent DFT-1/2 should still be used for other compounds. Nevertheless, there is no quantitative criterion prescribed for which implementation should work for an arbitrary insulator, which leads to severe ambiguity in this method. In this work, we analyze the impact of self-consistency in DFT-1/2 and shell DFT-1/2 calculations in insulators or semiconductors with ionic bonds, covalent bonds, and intermediate cases and show that self-consistency is required even for highly ionic insulators for globally better electronic structure details. The self-energy correction renders electrons more localized around the anions in self-consistent LDA-1/2. The well-known delocalization error of LDA is rectified, but with strong overcorrection, due to the presence of additional self-energy potential. However, in non-self-consistent LDA-1/2 calculations, the electron wave functions indicate that such localization is much more severe and beyond a reasonable range because the strong Coulomb repulsion is not counted in the Hamiltonian. Another common drawback of non-self-consistent LDA-1/2 is that the ionicity of the bonding gets substantially enhanced, and the bandgap can be enormously high in mixed ionic-covalent compounds like TiO2.
Collapse
Affiliation(s)
- Hanli Cui
- School of Integrated Circuits, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Shengxin Yang
- School of Integrated Circuits, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Kan-Hao Xue
- School of Integrated Circuits, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Jinhai Huang
- School of Integrated Circuits, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Xiangshui Miao
- School of Integrated Circuits, Huazhong University of Science and Technology, Wuhan 430074, China
| |
Collapse
|
5
|
Zobel JP, Kruse A, Baig O, Lochbrunner S, Bokarev SI, Kühn O, González L, Bokareva OS. Can range-separated functionals be optimally tuned to predict spectra and excited state dynamics in photoactive iron complexes? Chem Sci 2023; 14:1491-1502. [PMID: 36794199 PMCID: PMC9906774 DOI: 10.1039/d2sc05839a] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Accepted: 12/27/2022] [Indexed: 01/14/2023] Open
Abstract
Density functional theory is an efficient computational tool to investigate photophysical and photochemical processes in transition metal complexes, giving invaluable assistance in interpreting spectroscopic and catalytic experiments. Optimally tuned range-separated functionals are particularly promising, as they were created to address some of the fundamental deficiencies present in approximate exchange-correlation functionals. In this paper, we scrutinize the selection of optimally tuned parameters and its influence on the excited state dynamics, using the example of the iron complex [Fe(cpmp)2]2+ with push-pull ligands. Various tuning strategies are contemplated based on pure self-consistent DFT protocols, as well as on the comparison with experimental spectra and multireference CASPT2 results. The two most promising sets of optimal parameters are then employed to carry out nonadiabatic surface-hopping dynamics simulations. Intriguingly, we find that the two sets lead to very different relaxation pathways and timescales. While the set of optimal parameters from one of the self-consistent DFT protocols predicts the formation of long-lived metal-to-ligand charge transfer triplet states, the set in better agreement with CASPT2 calculations leads to deactivation in the manifold of metal-centered states, in better agreement with the experimental reference data. These results showcase the complexity of iron-complex excited state landscapes and the difficulty of obtaining an unambiguous parametrization of long-range corrected functionals without experimental input.
Collapse
Affiliation(s)
- J. Patrick Zobel
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 191090 ViennaAustria
| | - Ayla Kruse
- Institute of Physics, University of Rostock, Albert-Einstein-Straße 23-24 18059 Rostock Germany .,Department of Life, Light and Matter, University of Rostock 18051 Rostock Germany
| | - Omar Baig
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19 1090 Vienna Austria
| | - Stefan Lochbrunner
- Institute of Physics, University of Rostock, Albert-Einstein-Straße 23-24 18059 Rostock Germany .,Department of Life, Light and Matter, University of Rostock 18051 Rostock Germany
| | - Sergey I. Bokarev
- Institute of Physics, University of Rostock, Albert-Einstein-Straße 23-2418059 RostockGermany,Chemistry Department, Technical University of Munich, Lichtenbergstr. 4Garching 85748Germany
| | - Oliver Kühn
- Institute of Physics, University of Rostock, Albert-Einstein-Straße 23-24 18059 Rostock Germany
| | - Leticia González
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währingerstr. 19 1090 Vienna Austria
| | - Olga S. Bokareva
- Institute of Physics, University of Rostock, Albert-Einstein-Straße 23-2418059 RostockGermany,Institute of Physics, University of KasselHeinrich-Plett-Straße 4034132 KasselGermany
| |
Collapse
|
6
|
Lonsdale DR, Goerigk L. One-electron self-interaction error and its relationship to geometry and higher orbital occupation. J Chem Phys 2023; 158:044102. [PMID: 36725505 DOI: 10.1063/5.0129820] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Density Functional Theory (DFT) sees prominent use in computational chemistry and physics; however, problems due to the self-interaction error (SIE) pose additional challenges to obtaining qualitatively correct results. As an unphysical energy an electron exerts on itself, the SIE impacts most practical DFT calculations. We conduct an in-depth analysis of the one-electron SIE in which we replicate delocalization effects for simple geometries. We present a simple visualization of such effects, which may help in future qualitative analysis of the one-electron SIE. By increasing the number of nuclei in a linear arrangement, the SIE increases dramatically. We also show how molecular shape impacts the SIE. Two- and three-dimensional shapes show an even greater SIE stemming mainly from the exchange functional with some error compensation from the one-electron error, which we previously defined [D. R. Lonsdale and L. Goerigk, Phys. Chem. Chem. Phys. 22, 15805 (2020)]. Most tested geometries are affected by the functional error, while some suffer from the density error. For the latter, we establish a potential connection with electrons being unequally delocalized by the DFT methods. We also show how the SIE increases if electrons occupy higher-lying atomic orbitals; seemingly one-electron SIE free methods in a ground are no longer SIE free in excited states, which is an important insight for some popular, non-empirical density functional approximations (DFAs). We conclude that the erratic behavior of the SIE in even the simplest geometries shows that robust DFAs are needed. Our test systems can be used as a future benchmark or contribute toward DFT development.
Collapse
Affiliation(s)
- Dale R Lonsdale
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
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
|
Ye F, Xue K, Yu H, Yang S, Yuan J, Gu R, Xu M, Miao X. In Quest of Low‐Leakage Dynamic Random Access Memory Enabled by Doped TiO
2
Dielectrics. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200614] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Fan Ye
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
| | - Kan‐Hao Xue
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
- Hubei Yangtze Memory Laboratories Wuhan 430205 China
| | - Heng Yu
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
| | - Shengxin Yang
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
| | - Jun‐Hui Yuan
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
| | - Rongchuan Gu
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
| | - Ming Xu
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
- Hubei Yangtze Memory Laboratories Wuhan 430205 China
| | - Xiangshui Miao
- School of Integrated Circuits School of Optical and Electronic Information Huazhong University of Science and Technology Wuhan 430074 China
- Hubei Yangtze Memory Laboratories Wuhan 430205 China
| |
Collapse
|
9
|
Falletta S, Pasquarello A. Many-Body Self-Interaction and Polarons. PHYSICAL REVIEW LETTERS 2022; 129:126401. [PMID: 36179212 DOI: 10.1103/physrevlett.129.126401] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Revised: 08/26/2022] [Accepted: 08/30/2022] [Indexed: 06/16/2023]
Abstract
We address the many-body self-interaction in relation to polarons in density functional theory. Our study provides (i) a unified theoretical framework encompassing many-body and one-body forms of self-interaction and (ii) an efficient semilocal scheme for charge localization. Our theoretical formulation establishes a quantitative connection between the many-body and one-body forms of self-interaction in terms of electron screening, thereby conferring superiority to the concept of many-body self-interaction. Our semilocal methodology involves the use of a weak localized potential and applies equally to electron and hole polarons. We find that polarons free from many-body self-interaction have formation energies that are robust with respect to the functional adopted.
Collapse
Affiliation(s)
- Stefano Falletta
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
10
|
Chaudhary J, Aarzoo, Roy R, Roy RK. Screening the Band Shape of Molecules by Optimal Tuning of Range-Separated Hybrid Functional with TD-DFT: A Molecular Designing Approach. J Phys Chem A 2022; 126:5252-5264. [PMID: 35939642 DOI: 10.1021/acs.jpca.2c02204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the present article we have demonstrated the effectiveness of optimally tuned range-separated hybrid (RSH) functional to determine the electronic transitions from two fluorophore moieties (blue and yellow/orange) within a single white light emitter (WLE). First, the optimally tuned range separation parameter (ω) is calculated for two white emitting fluorophores (W1 and W2) already reported in the literature. The success of the optimally tuned RSH functional ω*B97XD, used in the TDDFT study, encouraged the authors to design eight new single organic white light emitters with frustrated energy transfer between the two individual fluorophore moieties (blue and yellow/orange). The simulated spectra (the band shapes, to be more specific) generated by TDDFT study and outcomes through natural transition orbital (NTO) and natural bond orbital (NBO) studies clearly demonstrate that all the designed eight organic molecules are potential white light emitters and can be synthesized in future.
Collapse
Affiliation(s)
- Jagrity Chaudhary
- Department of Chemistry, Birla Institute of Technology and Science, Pilani 333031, Rajasthan, India
| | - Aarzoo
- Department of Chemistry, Birla Institute of Technology and Science, Pilani 333031, Rajasthan, India
| | - Ritaban Roy
- Department of Chemistry, Birla Institute of Technology and Science, Pilani 333031, Rajasthan, India
| | - Ram Kinkar Roy
- Department of Chemistry, Birla Institute of Technology and Science, Pilani 333031, Rajasthan, India
| |
Collapse
|
11
|
Kretz B, Egger D. Accurate Non-Adiabatic Couplings from Optimally-Tuned Range-Separated Hybrid Functionals. J Chem Phys 2022; 157:101104. [DOI: 10.1063/5.0099854] [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
Precise theoretical calculations of non-adiabatic couplings, which describe the interaction between two Born-Oppenheimer surfaces, are important for the modeling of radiationless decay mechanisms in photochemical processes. Here, we demonstrate that accurate non-adiabatic couplings can be calculated in the framework of linear-response time-dependent density functional theory by using non-empirical, optimally-tuned range-separated hybrid (OT-RSH) functionals. We focus on molecular radicals, in which ultrafast non-radiative decay plays a crucial role, to find that the OT-RSH functional compares well to wave-function based reference data and competes with the accuracy of semi-empirical CAM-B3LYP calculations. Our findings show that the OT-RSH approach provides very accurate non-adiabatic couplings and, therefore, provides a computationally efficient alternative to wave-function based techniques.
Collapse
|
12
|
Majumdar S, Roy AK. Recent Advances in Cartesian-Grid DFT in Atoms and Molecules. Front Chem 2022; 10:926916. [PMID: 35936092 PMCID: PMC9354079 DOI: 10.3389/fchem.2022.926916] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2022] [Accepted: 06/09/2022] [Indexed: 11/30/2022] Open
Abstract
In the past several decades, density functional theory (DFT) has evolved as a leading player across a dazzling variety of fields, from organic chemistry to condensed matter physics. The simple conceptual framework and computational elegance are the underlying driver for this. This article reviews some of the recent developments that have taken place in our laboratory in the past 5 years. Efforts are made to validate a viable alternative for DFT calculations for small to medium systems through a Cartesian coordinate grid- (CCG-) based pseudopotential Kohn–Sham (KS) DFT framework using LCAO-MO ansatz. In order to legitimize its suitability and efficacy, at first, electric response properties, such as dipole moment (μ), static dipole polarizability (α), and first hyperpolarizability (β), are calculated. Next, we present a purely numerical approach in CCG for proficient computation of exact exchange density contribution in certain types of orbital-dependent density functionals. A Fourier convolution theorem combined with a range-separated Coulomb interaction kernel is invoked. This takes motivation from a semi-numerical algorithm, where the rate-deciding factor is the evaluation of electrostatic potential. Its success further leads to a systematic self-consistent approach from first principles, which is desirable in the development of optimally tuned range-separated hybrid and hyper functionals. Next, we discuss a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, by means of “adiabatic connection theorem” and virial theorem. Optical gaps in organic chromophores, dyes, linear/non-linear PAHs, and charge transfer complexes are faithfully reproduced. In short, CCG-DFT is shown to be a successful route for various practical applications in electronic systems.
Collapse
|
13
|
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
|
14
|
Prokopiou G, Hartstein M, Govind N, Kronik L. Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals. J Chem Theory Comput 2022; 18:2331-2340. [PMID: 35369687 PMCID: PMC9009176 DOI: 10.1021/acs.jctc.2c00082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Indexed: 11/29/2022]
Abstract
We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.
Collapse
Affiliation(s)
- Georgia Prokopiou
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Michal Hartstein
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Leeor Kronik
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
15
|
Yang J, Falletta S, Pasquarello A. One-Shot Approach for Enforcing Piecewise Linearity on Hybrid Functionals: Application to Band Gap Predictions. J Phys Chem Lett 2022; 13:3066-3071. [PMID: 35352960 DOI: 10.1021/acs.jpclett.2c00414] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We present an efficient procedure for constructing nonempirical hybrid functionals to accurately predict band gaps of extended systems. We determine mixing parameters by enforcing the generalized Koopmans' condition on localized electron states, which are achieved by inserting an optimized potential probe. Application of this scheme to a large set of materials yields band gaps with a mean error of 0.30 eV with respect to experiment. Next, we consider a perturbative one-shot approach in which the single-particle eigenvalues are calculated with the wave functions obtained at the semilocal level. In this way, the computational cost is reduced by ∼85% without loss of accuracy. The scheme is found to be robust upon consideration of different defect species and functional forms.
Collapse
Affiliation(s)
- Jing Yang
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Stefano Falletta
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
16
|
Sim E, Song S, Vuckovic S, Burke K. Improving Results by Improving Densities: Density-Corrected Density Functional Theory. J Am Chem Soc 2022; 144:6625-6639. [PMID: 35380807 DOI: 10.1021/jacs.1c11506] [Citation(s) in RCA: 34] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Density functional theory (DFT) calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving π-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This Perspective explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.
Collapse
Affiliation(s)
- Eunji Sim
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Suhwan Song
- Department of Chemistry, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Korea
| | - Stefan Vuckovic
- Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni,Campus Unisalento, 73100 Lecce, Italy.,Department of Chemistry & Pharmaceutical Sciences and Amsterdam Institute of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Kieron Burke
- Departments of Chemistry and of Physics, University of California, Irvine, California 92697, United States
| |
Collapse
|
17
|
Brütting M, Bahmann H, Kümmel S. Hybrid functionals with local range separation: Accurate atomization energies and reaction barrier heights. J Chem Phys 2022; 156:104109. [PMID: 35291795 DOI: 10.1063/5.0082957] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Range-separated hybrid approximations to the exchange-correlation density functional mix exact and semi-local exchange in a position-dependent manner. In their conventional form, the range separation is controlled by a constant parameter. Turning this constant into a density functional leads to a locally space-dependent range-separation function and thus a more powerful and flexible range-separation approach. In this work, we explore the self-consistent implementation of a local range-separated hybrid, taking into account a one-electron self-interaction correction and the behavior under uniform density scaling. We discuss different forms of the local range-separation function that depend on the electron density, its gradient, and the kinetic energy density. For test sets of atomization energies, reaction barrier heights, and total energies of atoms, we demonstrate that our best model is a clear improvement over common global range-separated hybrid functionals and can compete with density functionals that contain multiple empirical parameters. Promising results for equilibrium bond lengths, harmonic vibrational frequencies, and vertical ionization potentials further underline the potential and flexibility of our approach.
Collapse
Affiliation(s)
- Moritz Brütting
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Hilke Bahmann
- Physical and Theoretical Chemistry, Saarland University, 66123 Saarbrücken, Germany
| | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| |
Collapse
|
18
|
Nazarov VU. Breakdown of the ionization potential theorem of density functional theory in mesoscopic systems. J Chem Phys 2021; 155:194105. [PMID: 34800939 DOI: 10.1063/5.0070429] [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
The ionization potential (IP)-theorem of Kohn-Sham (KS) density functional theory (DFT) states that the energy of the highest occupied molecular orbital (HOMO) ϵHOMO equals the negative of the first IP, thus ascribing a physical meaning to one of the eigenvalues of the KS Hamiltonian. We scrutinize the fact that the validity of the IP-theorem relies critically on the electron density n(r), far from the system, to be determined by HOMO only, behaving as n(r)∼r→∞e-2-2ϵHOMOr. While this behavior always holds for finite systems, it does not hold for mesoscopic ones, such as quasi-two-dimensional (Q2D) electron gas or Q2D crystals. We show that this leads to the violation of the IP-theorem for the latter class of systems. This finding has a strong bearing on the role of the KS valence band with respect to the work-function problem in the mesoscopic case. Based on our results, we introduce a concept of the IP band structure as an observable alternative to its unphysical KS counterpart. A practical method of the determination of the IP band structure in terms of DFT quantities is provided.
Collapse
Affiliation(s)
- Vladimir U Nazarov
- Moscow Institute of Physics and Technology (National Research University), Dolgoprudny, Russian Federation
| |
Collapse
|
19
|
Ghosal A, Roy AK. A self-consistent systematic optimisation of range-separated hybrid functionals from first principles. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1983056] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Abhisek Ghosal
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Nadia, Mohanpur, India
| | - Amlan K. Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Nadia, Mohanpur, India
| |
Collapse
|
20
|
Marzari N, Ferretti A, Wolverton C. Electronic-structure methods for materials design. NATURE MATERIALS 2021; 20:736-749. [PMID: 34045704 DOI: 10.1038/s41563-021-01013-3] [Citation(s) in RCA: 38] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Accepted: 04/19/2021] [Indexed: 05/24/2023]
Abstract
The accuracy and efficiency of electronic-structure methods to understand, predict and design the properties of materials has driven a new paradigm in research. Simulations can greatly accelerate the identification, characterization and optimization of materials, with this acceleration driven by continuous progress in theory, algorithms and hardware, and by adaptation of concepts and tools from computer science. Nevertheless, the capability to identify and characterize materials relies on the predictive accuracy of the underlying physical descriptions, and on the ability to capture the complexity of realistic systems. We provide here an overview of electronic-structure methods, of their application to the prediction of materials properties, and of the different strategies employed towards the broader goals of materials design and discovery.
Collapse
Affiliation(s)
- Nicola Marzari
- Theory and Simulation of Materials (THEOS), and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
| | | | - Chris Wolverton
- Department of Materials Science and Engineering, Northwestern University, Evanston, IL, USA
| |
Collapse
|
21
|
Bajaj A, Kulik HJ. Molecular DFT+U: A Transferable, Low-Cost Approach to Eliminate Delocalization Error. J Phys Chem Lett 2021; 12:3633-3640. [PMID: 33826346 DOI: 10.1021/acs.jpclett.1c00796] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
While density functional theory (DFT) is widely applied for its combination of cost and accuracy, corrections (e.g., DFT+U) that improve it are often needed to tackle correlated transition-metal chemistry. In principle, the functional form of DFT+U, consisting of a set of localized atomic orbitals (AOs) and a quadratic energy penalty for deviation from integer occupations of those AOs, enables the recovery of the exact conditions of piecewise linearity and the derivative discontinuity. Nevertheless, for practical transition-metal complexes, where both atomic states and ligand orbitals participate in bonding, standard DFT+U can fail to eliminate delocalization error (DE). Here, we show that by introducing an alternative valence-state (i.e., molecular orbital or MO) basis to the DFT+U approach, we recover exact conditions in cases for which standard DFT+U corrections have no error-reducing effect. This MO-based DFT+U also eliminates DE where standard AO-based DFT+U is already successful. We demonstrate the transferability of our approach on representative transition-metal complexes with a range of ligand field strengths, electron configurations (i.e., from Sc to Zn), and spin states.
Collapse
Affiliation(s)
- Akash Bajaj
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
22
|
Kretz B, Egger DA. Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory. J Chem Theory Comput 2021; 17:357-366. [PMID: 33284603 DOI: 10.1021/acs.jctc.0c00858] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.
Collapse
Affiliation(s)
- Bernhard Kretz
- Department of Physics, Technical University of Munich, James-Franck-Str. 1, 85748 Garching, Germany
| | - David A Egger
- Department of Physics, Technical University of Munich, James-Franck-Str. 1, 85748 Garching, Germany
| |
Collapse
|
23
|
Klinger A, Lindorfer D, Müh F, Renger T. Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer. J Chem Phys 2020; 153:215103. [PMID: 33291900 DOI: 10.1063/5.0027994] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The intermolecular contribution to the spectral density of the exciton-vibrational coupling of the homotrimeric Fenna-Matthews-Olson (FMO) light-harvesting protein of green sulfur bacteria P. aestuarii is analyzed by combining a normal mode analysis of the protein with the charge density coupling method for the calculation of local transition energies of the pigments. Correlations in site energy fluctuations across the whole FMO trimer are found at low vibrational frequencies. Including, additionally, the high-frequency intrapigment part of the spectral density, extracted from line-narrowing spectra, we study intra- and intermonomer exciton transfer. Whereas the intrapigment part of the spectral density is important for fast intramonomer exciton relaxation, the intermolecular contributions (due to pigment-environment coupling) determine the intermonomer exciton transfer. Neither the variations of the local Huang-Rhys factors nor the correlations in site energy fluctuations have a critical influence on energy transfer. At room temperature, the intermonomer transfer in the FMO protein occurs on a 10 ps time scale, whereas intramonomer exciton equilibration is roughly two orders of magnitude faster. At cryogenic temperatures, intermonomer transfer limits the lifetimes of the lowest exciton band. The lifetimes are found to increase between 20 ps in the center of this band up to 100 ps toward lower energies, which is in very good agreement with the estimates from hole burning data. Interestingly, exciton delocalization in the FMO monomers is found to slow down intermonomer energy transfer, at both physiological and cryogenic temperatures.
Collapse
Affiliation(s)
- Alexander Klinger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Dominik Lindorfer
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Frank Müh
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
24
|
Carter KJB, Grossman JC. Design Rules for Transparent Push-Pull Electron Acceptors: A Case Study on Perylenediimide Derivatives. J Phys Chem Lett 2020; 11:9265-9271. [PMID: 33079546 DOI: 10.1021/acs.jpclett.0c02485] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Transparent photovoltaics are receiving increased attention for their wide range of applicability, but there have been few attempts to systematically tune organic materials to achieve transparency. In this Letter, we study the influence of simple structural modifications on the photoabsorption spectrum of the nonfullerene electron acceptor perylenediimide. Motivated by push-pull design strategies, we explore the effects of applying electron-donating functional groups with varying strengths using three design motifs. We demonstrate that relative Mulliken electronegativity, which can be easily computed using an optimally tuned screened range-separated hybrid functional, is a useful metric for finding suitable donor groups. We also demonstrate that donor groups that include a conjugated spacer are crucial to obtaining a planar molecule with uniform conjugation and low-energy charge-transfer excitations. These simple design rules can be used to design near-infrared absorbing transparent electron acceptors based on perylenediimide and other promising molecular platforms.
Collapse
Affiliation(s)
- Ki-Jana B Carter
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Jeffrey C Grossman
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
25
|
Savin A. Models and corrections: Range separation for electronic interaction—Lessons from density functional theory. J Chem Phys 2020; 153:160901. [DOI: 10.1063/5.0028060] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Affiliation(s)
- Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, 4 Place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|