1
|
Chaudhry I, Hu G, Ye H, Jensen L. Toward Modeling the Complexity of the Chemical Mechanism in SERS. ACS NANO 2024. [PMID: 39087679 DOI: 10.1021/acsnano.4c07198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/02/2024]
Abstract
Surface-enhanced Raman scattering (SERS) provides detailed information about the binding of molecules at interfaces and their interactions with the local environment due to the large enhancement of Raman scattering. This enhancement arises from a combination of the electromagnetic mechanism (EM) and chemical mechanism (CM). While it is commonly accepted that EM gives rise to most of the enhancement, large spectral changes originate from CM. To elucidate the rich information contained in SERS spectra about molecules at interfaces, a comprehensive understanding of the enhancement mechanisms is necessary. In this Perspective, we discuss the current understanding of the enhancement mechanisms and highlight their interplay in complex local environments. We will also discuss emerging areas where the development of computational and theoretical models is needed with specific attention given to how the CM contributes to the spectral changes. Future efforts in modeling should focus on overcoming the challenges presented in this review in order to capture the complexity of CM in SERS.
Collapse
Affiliation(s)
- Imran Chaudhry
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Gaohe Hu
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Hepeng Ye
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Lasse Jensen
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| |
Collapse
|
2
|
Neufeld O, Tancogne-Dejean N, Rubio A. Benchmarking Functionals for Strong-Field Light-Matter Interactions in Adiabatic Time-Dependent Density Functional Theory. J Phys Chem Lett 2024; 15:7254-7264. [PMID: 38976844 PMCID: PMC11261632 DOI: 10.1021/acs.jpclett.4c01383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2024] [Revised: 06/29/2024] [Accepted: 07/03/2024] [Indexed: 07/10/2024]
Abstract
In recent years, time-dependent density functional theory (TDDFT) has been extensively employed for highly nonlinear optics in molecules and solids, including high harmonic generation (HHG), photoemission, and more. TDDFT exhibits a relatively low numerical cost while still describing both light-matter and electron-electron interactions ab initio, making it highly appealing. However, the majority of implementations of the theory utilize the simplest possible approximations for the exchange-correlation (XC) functional-either the local density or generalized gradient approximations, which are traditionally considered to have rather poor chemical accuracy. We present the first systematic study of the XC functional effect on molecular HHG, testing various levels of theory. Our numerical results suggest justification for using simpler approximations for the XC functional, showing that hybrid and meta functionals (as well as Hartree-Fock) can, at times, lead to poor and unphysical results. The specific source of the failure in more elaborate functionals should be topic of future work, but we hypothesize that its origin might be connected to the adiabatic approximation of TDDFT.
Collapse
Affiliation(s)
- Ofer Neufeld
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Nicolas Tancogne-Dejean
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Angel Rubio
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
- Center
for Computational Quantum Physics (CCQ), The Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
3
|
Shepard C, Zhou R, Bost J, Carney TE, Yao Y, Kanai Y. Efficient exact exchange using Wannier functions and other related developments in planewave-pseudopotential implementation of RT-TDDFT. J Chem Phys 2024; 161:024111. [PMID: 38984957 DOI: 10.1063/5.0211238] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2024] [Accepted: 06/19/2024] [Indexed: 07/11/2024] Open
Abstract
The plane-wave pseudopotential (PW-PP) formalism is widely used for the first-principles electronic structure calculation of extended periodic systems. The PW-PP approach has also been adapted for real-time time-dependent density functional theory (RT-TDDFT) to investigate time-dependent electronic dynamical phenomena. In this work, we detail recent advances in the PW-PP formalism for RT-TDDFT, particularly how maximally localized Wannier functions (MLWFs) are used to accelerate simulations using the exact exchange. We also discuss several related developments, including an anti-Hermitian correction for the time-dependent MLWFs (TD-MLWFs) when a time-dependent electric field is applied, the refinement procedure for TD-MLWFs, comparison of the velocity and length gauge approaches for applying an electric field, and elimination of long-range electrostatic interaction, as well as usage of a complex absorbing potential for modeling isolated systems when using the PW-PP formalism.
Collapse
Affiliation(s)
- Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - John Bost
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Thomas E Carney
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
- Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
- Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
4
|
Taylor JT, Tozer DJ, Curchod BFE. On the Topological Phase around Conical Intersections with Tamm-Dancoff Linear-Response Time-Dependent Density Functional Theory. J Phys Chem A 2024; 128:5314-5320. [PMID: 38919046 PMCID: PMC11247484 DOI: 10.1021/acs.jpca.4c02503] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2024] [Revised: 05/30/2024] [Accepted: 06/17/2024] [Indexed: 06/27/2024]
Abstract
Regions of nuclear-configuration space away from the Franck-Condon geometry can prove problematic for some electronic structure methods, given the propensity of such regions to possess conical intersections, i.e., (highly connected) points of degeneracy between potential energy surfaces. With the likelihood (perhaps even inevitability) for nonadiabatic dynamics simulations to explore molecular geometries in close proximity to conical intersections, it is vital that the performance of electronic structure methods is routinely examined in this context. In a recent paper [Taylor, J. T. J. Chem. Phys. 2023, 159, 214115.], the ability of linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) to provide a proper description of conical intersections, in terms of their topology and topography, was investigated, with particular attention paid to conical intersections between two excited electronic states. For the same prototypical molecules, protonated formaldimine and pyrazine, we herein consider whether AA LR-TDDFT can correctly reproduce the topological phase accumulated by the adiabatic electronic wave function upon traversing a closed path around an excited-to-excited state conical intersection despite not using the appropriate quadratic-response nonadiabatic coupling vectors. Equally, we probe the ability of the ground-to-excited state intersection ring exhibited by AA LR-TDDFT in protonated formaldimine to give rise to a similar topological phase in spite of its incorrect dimensionality.
Collapse
Affiliation(s)
- Jack T. Taylor
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David J. Tozer
- Department
of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - Basile F. E. Curchod
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
5
|
Lieberherr AZ, Gori-Giorgi P, Giesbertz KJH. Optimal Transport Distances to Characterize Electronic Excitations. J Chem Theory Comput 2024; 20:5635-5642. [PMID: 38874479 PMCID: PMC11238536 DOI: 10.1021/acs.jctc.4c00289] [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/05/2024] [Revised: 05/31/2024] [Accepted: 06/03/2024] [Indexed: 06/15/2024]
Abstract
Understanding the character of electronic excitations is important in computational and reaction mechanistic studies, but their classification from simulations remains an open problem. Distances based on optimal transport have proven very useful in a plethora of classification problems and, therefore, seem a natural tool to try to tackle this challenge. We propose and investigate a new diagnostic Θ based on the Sinkhorn divergence from optimal transport. We evaluate a k-NN classification algorithm on Θ, the popular Λ diagnostic, and their combination, and assess their performance in labeling excitations, finding that (i) the combination only slightly improves the classification, (ii) Rydberg excitations are not separated well in any setting, and (iii) Θ breaks down for charge transfer in small molecules. We then define a length-scale-normalized version of Θ and show that the result correlates closely with Λ for results obtained with Gaussian basis functions. Finally, we discuss the orbital dependence of our approach and explore an orbital-independent version. Using an optimized combination of the optimal transport and overlap diagnostics together with a different metric is in our opinion the most promising for future classification studies.
Collapse
Affiliation(s)
- Annina Z. Lieberherr
- Department
of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QZ, U.K.
| | - Paola Gori-Giorgi
- Department
of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
- Microsoft
Research AI for Science, Evert van de Beekstraat 354, 1118CZ Schiphol, The Netherlands
| | - Klaas J. H. Giesbertz
- Department
of Chemistry and Pharmaceutical Sciences, Amsterdam Institute of Molecular
and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| |
Collapse
|
6
|
Lv ZP, Srivastava D, Conley K, Ruoko TP, Xu H, Lightowler M, Hong X, Cui X, Huang Z, Yang T, Wang HY, Karttunen AJ, Bergström L. Visualizing Noncovalent Interactions and Property Prediction of Submicron-Sized Charge-Transfer Crystals from ab-initio Determined Structures. SMALL METHODS 2024; 8:e2301229. [PMID: 38528393 DOI: 10.1002/smtd.202301229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Revised: 03/03/2024] [Indexed: 03/27/2024]
Abstract
The charge-transfer (CT) interactions between organic compounds are reflected in the (opto)electronic properties. Determining and visualizing crystal structures of CT complexes are essential for the design of functional materials with desirable properties. Complexes of pyranine (PYR), methyl viologen (MV), and their derivatives are the most studied water-based CT complexes. Nevertheless, very few crystal structures of CT complexes have been reported so far. In this study, the structures of two PYRs-MVs CT crystals and a map of the noncovalent interactions using 3D electron diffraction (3DED) are reported. Physical properties, e.g., band structure, conductivity, and electronic spectra of the CT complexes and their crystals are investigated and compared with a range of methods, including solid and liquid state spectroscopies and highly accurate quantum chemical calculations based on density functional theory (DFT). The combination of 3DED, spectroscopy, and DFT calculation can provide important insight into the structure-property relationship of crystalline CT materials, especially for submicrometer-sized crystals.
Collapse
Affiliation(s)
- Zhong-Peng Lv
- Department of Applied Physics, Aalto University, Espoo, FI 02150, Finland
| | - Divya Srivastava
- Department of Chemistry and Materials Science, Aalto University, Espoo, FI 02150, Finland
| | - Kevin Conley
- Department of Chemistry and Materials Science, Aalto University, Espoo, FI 02150, Finland
| | - Tero-Petri Ruoko
- Faculty of Engineering and Natural Sciences, Tampere University, Tampere, FI-33720, Finland
| | - Hongyi Xu
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, SE 10691, Sweden
| | - Molly Lightowler
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, SE 10691, Sweden
| | - Xiaodan Hong
- Department of Applied Physics, Aalto University, Espoo, FI 02150, Finland
| | - Xiaoqi Cui
- Department of Electronics and Nanoengineering, Aalto University, Espoo, FI 02150, Finland
| | - Zhehao Huang
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, SE 10691, Sweden
| | - Taimin Yang
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, SE 10691, Sweden
| | - Hai-Ying Wang
- School of Environmental Science, Nanjing Xiaozhuang University, Nanjing, 211171, P. R. China
| | - Antti J Karttunen
- Department of Chemistry and Materials Science, Aalto University, Espoo, FI 02150, Finland
| | - Lennart Bergström
- Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, SE 10691, Sweden
| |
Collapse
|
7
|
Osterloh WR, Desbois N, Conradie J, Gros CP, Kadish KM, Ghosh A. Inverse Hypercorroles. Inorg Chem 2024; 63:8739-8749. [PMID: 38696617 PMCID: PMC11094798 DOI: 10.1021/acs.inorgchem.4c00344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 04/08/2024] [Accepted: 04/10/2024] [Indexed: 05/04/2024]
Abstract
Ground-state and time-dependent density functional theory (TDDFT) calculations with the long-range-corrected, Coulomb-attenuating CAMY-B3LYP exchange-correlation functional and large, all-electron STO-TZ2P basis sets have been used to examine the potential "inverse hypercorrole" character of meso-p-nitrophenyl-appended dicyanidocobalt(III) corrole dianions. The effect is most dramatic for 5,15-bis(p-nitrophenyl) derivatives, where it manifests itself in intense NIR absorptions. The 10-aryl groups in these complexes play a modulatory role, as evinced by experimental UV-visible spectroscopic and electrochemical data for a series of 5,15-bis(p-nitrophenyl) dicyanidocobalt(III) corroles. TDDFT (CAMY-B3LYP) calculations ascribe these features clearly to a transition from the corrole's a2u-like HOMO (retaining the D4h irrep used for metalloporphyrins) to a nitrophenyl-based LUMO. The outward nature of this transition contrasts with the usual phenyl-to-macrocycle direction of charge transfer transitions in many hyperporphyrins and hypercorroles; thus, the complexes studied are aptly described as inverse hypercorroles.
Collapse
Affiliation(s)
- W. Ryan Osterloh
- ICMUB
(UMR CNRS 6302), Université de Bourgogne, 9, Avenue A. Savary, BP 47870, 21078 Dijon Cedex, France
- Department
of Chemistry, University of Houston, Houston, Texas 77204-5003, United
States
| | - Nicolas Desbois
- ICMUB
(UMR CNRS 6302), Université de Bourgogne, 9, Avenue A. Savary, BP 47870, 21078 Dijon Cedex, France
| | - Jeanet Conradie
- Department
of Chemistry, UiT − The Arctic University
of Norway, N-9037 Tromso̷, Norway
- Department
of Chemistry, University of the Free State, 9300 Bloemfontein, Republic of South Africa
| | - Claude P. Gros
- ICMUB
(UMR CNRS 6302), Université de Bourgogne, 9, Avenue A. Savary, BP 47870, 21078 Dijon Cedex, France
| | - Karl M. Kadish
- Department
of Chemistry, University of Houston, Houston, Texas 77204-5003, United
States
| | - Abhik Ghosh
- Department
of Chemistry, University of the Free State, 9300 Bloemfontein, Republic of South Africa
| |
Collapse
|
8
|
Avagliano D, Skreta M, Arellano-Rubach S, Aspuru-Guzik A. DELFI: a computer oracle for recommending density functionals for excited states calculations. Chem Sci 2024; 15:4489-4503. [PMID: 38516092 PMCID: PMC10952086 DOI: 10.1039/d3sc06440a] [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: 11/30/2023] [Accepted: 02/05/2024] [Indexed: 03/23/2024] Open
Abstract
Density functional theory (DFT) is the workhorse of computational quantum chemistry. One of its main limitations is that choosing the right functional is a non-trivial task left for human experts. The choice is particularly hard for excited state calculations when using its time-dependent formulation (TD-DFT). This is due to the approximations of the method, but also because the photophysical properties of a molecule are defined by a manifold of states that all need to be properly described. This includes not only the relative energy of the states, but also capturing the correct character, order, and intensity of the transitions. In this work, we developed a neural network to recommend functionals to be used on molecules for TD-DFT calculations, by simultaneously considering all these properties for a manifold of states. This was possible by developing a scoring system to define the accuracy of an excited state's calculation against a higher-accuracy reference. The scoring system is generalizable to any level of theory; we here applied it to evaluate the performance of common functionals of different rungs against a higher accuracy method on a large set of organic molecules. The results are collected in a database that we released and made open, providing four million data points to the community for future applications. The scoring system assigns a value between zero and one hundred to each functional for each molecule, transforming the complicated task of learning photophysical properties into a simpler regression task. We used the dataset to train a graph attention neural network to predict the scores for unseen molecules. We call this oracle DELFI (Data-driven EvaLuation of Functionals by Inference), which can be used to quickly screen and predict the ranking of functionals to calculate the optical properties of organic molecules. We validated DELFI in two in silico experiments: choosing a common functional for a series of spiropyran-merocyanine isomers and a unique functional to screen a large dataset of over 50 000 organic photovoltaic molecules, for which an extensive benchmark would be unfeasible. A corresponding web application allows DELFI to be easily run and the results to be analyzed, alleviating the hurdle of choosing the right functional for TD-DFT calculations.
Collapse
Affiliation(s)
- Davide Avagliano
- Department of Chemistry, University of Toronto 80 St. George Street Toronto ON M5S 3H6 Canada
- Department of Computer Science, University of Toronto 40 St. George Street Toronto ON M5S 2E4 Canada
| | - Marta Skreta
- Department of Computer Science, University of Toronto 40 St. George Street Toronto ON M5S 2E4 Canada
- Vector Institute for Artificial Intelligence 661 University Ave. Suite 710 ON M5G 1M1 Toronto Canada
| | | | - Alán Aspuru-Guzik
- Department of Chemistry, University of Toronto 80 St. George Street Toronto ON M5S 3H6 Canada
- Department of Computer Science, University of Toronto 40 St. George Street Toronto ON M5S 2E4 Canada
- Vector Institute for Artificial Intelligence 661 University Ave. Suite 710 ON M5G 1M1 Toronto Canada
- Department of Materials Science & Engineering, University of Toronto 184 College St Toronto M5S 3E4 Canada
- Department of Chemical Engineering & Applied Chemistry, University of Toronto 200 College St ON M5S 3E5 Toronto Canada
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR) 66118 University Ave. M5G 1M1 Toronto Canada
- Acceleration Consortium 80 St George St M5S 3H6 Toronto Canada
| |
Collapse
|
9
|
Crisci L, Coppola F, Petrone A, Rega N. Tuning ultrafast time-evolution of photo-induced charge-transfer states: A real-time electronic dynamics study in substituted indenotetracene derivatives. J Comput Chem 2024; 45:210-221. [PMID: 37706600 DOI: 10.1002/jcc.27231] [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: 07/20/2023] [Revised: 08/31/2023] [Accepted: 09/05/2023] [Indexed: 09/15/2023]
Abstract
Photo-induced charge transfer (CT) states are pivotal in many technological and biological processes. A deeper knowledge of such states is mandatory for modeling the charge migration dynamics. Real-time time-dependent density functional theory (RT-TD-DFT) electronic dynamics simulations are employed to explicitly observe the electronic density time-evolution upon photo-excitation. Asymmetrically substituted indenotetracene molecules, given their potential application as n-type semiconductors in organic photovoltaic materials, are here investigated. Effects of substituents with different electron-donating characters are analyzed in terms of the overall electronic energy spacing and resulting ultrafast CT dynamics through linear response (LR-)TD-DFT and RT-TD-DFT based approaches. The combination of the computational techniques here employed provided direct access to the electronic density reorganization in time and to its spatial and rational representation in terms of molecular orbital occupation time evolution. Such results can be exploited to design peculiar directional charge dynamics, crucial when photoactive materials are used for light-harvesting applications.
Collapse
Affiliation(s)
- Luigi Crisci
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Normale Superiore di Pisa, Pisa, Italy
| | | | - Alessio Petrone
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| | - Nadia Rega
- Department of Chemical Sciences, University of Napoli Federico II, Complesso Universitario di M.S. Angelo, Naples, Italy
- Scuola Superiore Meridionale, Naples, Italy
- Istituto Nazionale Di Fisica Nucleare, Sezione di Napoli, Complesso Universitario di M.S. Angelo ed. 6, Naples, Italy
| |
Collapse
|
10
|
Giarrusso S, Loos PF. Exact Excited-State Functionals of the Asymmetric Hubbard Dimer. J Phys Chem Lett 2023; 14:8780-8786. [PMID: 37739406 PMCID: PMC10561271 DOI: 10.1021/acs.jpclett.3c02052] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Accepted: 08/31/2023] [Indexed: 09/24/2023]
Abstract
The exact functionals associated with the (singlet) ground state and the two singlet excited states of the asymmetric Hubbard dimer at half-filling are calculated using both Levy's constrained search and Lieb's convex formulation. While the ground-state functional is, as is commonly known, a convex function with respect to the density, the functional associated with the doubly excited state is found to be concave. Also, because the density-potential mapping associated with the first excited state is noninvertible, its "functional" is a partial, multivalued function composed of one concave and one convex branch that correspond to two separate domains of the external potential. Remarkably, it is found that, although the one-to-one mapping between density and external potential may not apply (as in the case of the first excited state), each state-specific energy and corresponding universal functional are "functions" whose derivatives are each other's inverse, just as in the ground state formalism.
Collapse
Affiliation(s)
- Sara Giarrusso
- Laboratoire de Chimie et
Physique Quantiques (UMR 5626), Université
de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et
Physique Quantiques (UMR 5626), Université
de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
11
|
Ng BY, Zhou ZJ, Liu TT, Yoskamtorn T, Li G, Wu TS, Soo YL, Wu XP, Tsang SCE. Photo-Induced Active Lewis Acid-Base Pairs in a Metal-Organic Framework for H 2 Activation. J Am Chem Soc 2023; 145:19312-19320. [PMID: 37611205 PMCID: PMC10485891 DOI: 10.1021/jacs.3c05244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2023] [Indexed: 08/25/2023]
Abstract
The establishment of active sites as the frustrated Lewis pair (FLP) has recently attracted much attention ranging from homogeneous to heterogeneous systems in the field of catalysis. Their unquenched reactivity of Lewis acid and base pairs in close proximity that are unable to form stable adducts has been shown to activate small molecules such as dihydrogen heterolytically. Herein, we show that grafted Ru metal-organic framework-based catalysts prepared via N-containing linkers are rather catalytically inactive for H2 activation despite the application of elevated temperatures. However, upon light illumination, charge polarization of the anchored Ru bipyridine complex can form a transient Lewis acid-base pair, Ru+-N- via metal-to-ligand charge transfer, as confirmed by time-dependent density functional theory (TDDFT) calculations to carry out effective H2-D2 exchange. FTIR and 2-D NMR endorse the formation of such reactive intermediate(s) upon light irradiation.
Collapse
Affiliation(s)
- Bryan
Kit Yue Ng
- Department
of Chemistry, University of Oxford, Oxford OX1 3QR, U.K.
| | - Zi-Jian Zhou
- Key
Laboratory for Advanced Materials, Centre for Computational Chemistry
and Research Institute of Industrial Catalysis, East China University of Science and Technology, Shanghai 200237, People’s Republic of China
| | - Ting-Ting Liu
- Key
Laboratory for Advanced Materials, Centre for Computational Chemistry
and Research Institute of Industrial Catalysis, East China University of Science and Technology, Shanghai 200237, People’s Republic of China
| | | | - Guangchao Li
- Department
of Chemistry, University of Oxford, Oxford OX1 3QR, U.K.
| | - Tai-Sing Wu
- National
Synchrotron Radiation Research Center, 101 Hsin-Ann Road, Hsinchu 30076, Taiwan
| | - Yun-Liang Soo
- Department
of Physics, National Tsing Hua University, Hsin-chu 30013, Taiwan
| | - Xin-Ping Wu
- Key
Laboratory for Advanced Materials, Centre for Computational Chemistry
and Research Institute of Industrial Catalysis, East China University of Science and Technology, Shanghai 200237, People’s Republic of China
| | | |
Collapse
|
12
|
Herring C, Montemore MM. Recent Advances in Real-Time Time-Dependent Density Functional Theory Simulations of Plasmonic Nanostructures and Plasmonic Photocatalysis. ACS NANOSCIENCE AU 2023; 3:269-279. [PMID: 37601917 PMCID: PMC10436373 DOI: 10.1021/acsnanoscienceau.2c00061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Revised: 05/02/2023] [Accepted: 05/08/2023] [Indexed: 08/22/2023]
Abstract
Plasmonic catalysis provides a possible means for driving chemical reactions under relatively mild conditions. Rational design of these systems is impeded by the difficulty in understanding the electron dynamics and their interplay with reactions. Real-time, time-dependent density functional theory (RT-TDDFT) can provide dynamic information on excited states in plasmonic systems, including those relevant to plasmonic catalysis, at time scales and length scales that are otherwise out of reach of many experimental techniques. Here, we discuss previous RT-TDDFT studies of plasmonic systems, focusing on recent work that gains insight into plasmonic catalysis. These studies provide insight into plasmon dynamics, including size effects and the role of specific electronic states. Further, these studies provide significant insight into mechanisms underlying plasmonic catalysis, showing the importance of charge transfer between metal and adsorbate states, as well as local field enhancement, in different systems.
Collapse
Affiliation(s)
- Connor
J. Herring
- Department of Chemical and Biomolecular
Engineering, Tulane University, New Orleans, Louisiana 70115, United States
| | - Matthew M. Montemore
- Department of Chemical and Biomolecular
Engineering, Tulane University, New Orleans, Louisiana 70115, United States
| |
Collapse
|
13
|
Ranka K, Isborn CM. Size-dependent errors in real-time electron density propagation. J Chem Phys 2023; 158:2887545. [PMID: 37125706 DOI: 10.1063/5.0142515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Accepted: 04/14/2023] [Indexed: 05/02/2023] Open
Abstract
Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree-Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.
Collapse
Affiliation(s)
- Karnamohit Ranka
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| | - Christine M Isborn
- Chemistry and Biochemistry, University of California Merced, Merced, California 95343, USA
| |
Collapse
|
14
|
Brütting M, Foerster JM, Kümmel S. Understanding Primary Charge Separation in the Heliobacterial Reaction Center. J Phys Chem Lett 2023; 14:3092-3102. [PMID: 36951395 DOI: 10.1021/acs.jpclett.3c00377] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.
Collapse
|
15
|
Zhang J, Tang Z, Zhang X, Zhu H, Zhao R, Lu Y, Gao J. Target State Optimized Density Functional Theory for Electronic Excited and Diabatic States. J Chem Theory Comput 2023; 19:1777-1789. [PMID: 36917687 DOI: 10.1021/acs.jctc.2c01317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
A flexible self-consistent field method, called target state optimization (TSO), is presented for exploring electronic excited configurations and localized diabatic states. The key idea is to partition molecular orbitals into different subspaces according to the excitation or localization pattern for a target state. Because of the orbital-subspace constraint, orbitals belonging to different subspaces do not mix. Furthermore, the determinant wave function for such excited or diabatic configurations can be variationally optimized as a ground state procedure, unlike conventional ΔSCF methods, without the possibility of collapsing back to the ground state or other lower-energy configurations. The TSO method can be applied both in Hartree-Fock theory and in Kohn-Sham density functional theory (DFT). The density projection procedure and the working equations for implementing the TSO method are described along with several illustrative applications. For valence excited states of organic compounds, it was found that the computed excitation energies from TSO-DFT and time-dependent density functional theory (TD-DFT) are of similar quality with average errors of 0.5 and 0.4 eV, respectively. For core excitation, doubly excited states and charge-transfer states, the performance of TSO-DFT is clearly superior to that from conventional TD-DFT calculations. It is shown that variationally optimized charge-localized diabatic states can be defined using TSO-DFT in energy decomposition analysis to gain both qualitative and quantitative insights on intermolecular interactions. Alternatively, the variational diabatic states may be used in molecular dynamics simulation of charge transfer processes. The TSO method can also be used to define basis states in multistate density functional theory for excited states through nonorthogonal state interaction calculations. The software implementing TSO-DFT can be accessed from the authors.
Collapse
Affiliation(s)
- Jun Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Zhen Tang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Xiaoyong Zhang
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Hong Zhu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China
| | - Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,Institute of Theoretical Chemistry, Jilin University, Changchun, 130023 Jilin, P. R. China
| | - Yangyi Lu
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, P. R. China.,School of Chemical Biology and Biotechnology, Peking University Shenzhen Graduate School, Shenzhen 518055, P. R. China.,Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
16
|
Otis L, Neuscamman E. A promising intersection of excited‐state‐specific methods from quantum chemistry and quantum Monte Carlo. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2023. [DOI: 10.1002/wcms.1659] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Affiliation(s)
- Leon Otis
- Department of Physics University of California Berkeley Berkeley California USA
| | - Eric Neuscamman
- Department of Chemistry University of California Berkeley Berkeley California USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory Berkeley California USA
| |
Collapse
|
17
|
Kempfer-Robertson EM, Haase MN, Bersson JS, Avdic I, Thompson LM. Role of Exact Exchange in Difference Projected Double-Hybrid Density Functional Theory for Treatment of Local, Charge Transfer, and Rydberg Excitations. J Phys Chem A 2022; 126:8058-8069. [PMID: 36269072 DOI: 10.1021/acs.jpca.2c04338] [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/30/2022]
Abstract
Difference approaches to the study of excited states have undergone a renaissance in recent years, with the development of a plethora of algorithms for locating self-consistent field approximations to excited states. Density functional theory is likely to offer the best balance of cost and accuracy for difference approaches, and yet there has been little investigation of how the parametrization of density functional approximations affects performance. In this work, we aim to explore the role of the global Hartree-Fock exchange parameter in tuning accuracy of different excitation types within the framework of the recently introduced difference projected double-hybrid density functional theory approach and contrast the performance with conventional time-dependent double-hybrid density functional theory. Difference projected double-hybrid density functional theory was demonstrated to give vertical excitation energies with average error and standard deviation with respect to multireference perturbation theory comparable to more expensive linear-response coupled cluster approaches ( J. Chem. Phys.2020, 153, 074103). However, despite benchmarking of local excitations, there has been no investigation of the methods performance for charge transfer or Rydberg excitations. In this work we report a new benchmark of charge transfer, Rydberg, and local excited state vertical excitation energies and examine how the exact Hartree-Fock exchange affects the benchmark performance to provide a deeper understanding of how projection and nonlocal correlation balance differing sources of error in the ground and excited states.
Collapse
Affiliation(s)
| | - Meagan N Haase
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Jonathan S Bersson
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Irma Avdic
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| | - Lee M Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky40205, United States
| |
Collapse
|
18
|
Carter-Fenk K, Cunha LA, Arias-Martinez JE, Head-Gordon M. Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States. J Phys Chem Lett 2022; 13:9664-9672. [PMID: 36215404 DOI: 10.1021/acs.jpclett.2c02564] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| |
Collapse
|
19
|
Arias-Martinez JE, Cunha LA, Oosterbaan KJ, Lee J, Head-Gordon M. Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach. Phys Chem Chem Phys 2022; 24:20728-20741. [PMID: 36004629 DOI: 10.1039/d2cp01998a] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core-valence separated (FC-CVS) EOM-CCSD approach - a successful variant of EOM-CCSD tailored towards core excitations.
Collapse
Affiliation(s)
- Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Katherine J Oosterbaan
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA.
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York 10027, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA. .,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
20
|
Monino E, Boggio-Pasqua M, Scemama A, Jacquemin D, Loos PF. Reference Energies for Cyclobutadiene: Automerization and Excited States. J Phys Chem A 2022; 126:4664-4679. [PMID: 35820169 DOI: 10.1021/acs.jpca.2c02480] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Cyclobutadiene is a well-known playground for theoretical chemists and is particularly suitable to test ground- and excited-state methods. Indeed, due to its high spatial symmetry, especially at the D4h square geometry but also in the D2h rectangular arrangement, the ground and excited states of cyclobutadiene exhibit multiconfigurational characters and single-reference methods, such as standard adiabatic time-dependent density-functional theory (TD-DFT) or standard equation-of-motion coupled cluster (EOM-CC), are notoriously known to struggle in such situations. In this work, using a large panel of methods and basis sets, we provide an extensive computational study of the automerization barrier (defined as the difference between the square and rectangular ground-state energies) and the vertical excitation energies at D2h and D4h equilibrium structures. In particular, selected configuration interaction (SCI), multireference perturbation theory (CASSCF, CASPT2, and NEVPT2), and coupled-cluster (CCSD, CC3, CCSDT, CC4, and CCSDTQ) calculations are performed. The spin-flip formalism, which is known to provide a qualitatively correct description of these diradical states, is also tested within TD-DFT (combined with numerous exchange-correlation functionals) and the algebraic diagrammatic construction [ADC(2)-s, ADC(2)-x, and ADC(3)]. A theoretical best estimate is defined for the automerization barrier and for each vertical transition energy.
Collapse
Affiliation(s)
- Enzo Monino
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
21
|
Hancock AC, Goerigk L. Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models. RSC Adv 2022; 12:13014-13034. [PMID: 35520129 PMCID: PMC9062889 DOI: 10.1039/d2ra01703b] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Accepted: 04/12/2022] [Indexed: 01/21/2023] Open
Abstract
Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.
Collapse
Affiliation(s)
- Amy C Hancock
- School of Chemistry, The University of Melbourne Parkville Australia +61-3-8344-6784
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne Parkville Australia +61-3-8344-6784
| |
Collapse
|
22
|
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
|
23
|
Lebel M, Very T, Gloaguen E, Tardivel B, Mons M, Brenner V. Excited States Computation of Models of Phenylalanine Protein Chains: TD-DFT and Composite CC2/TD-DFT Protocols. Int J Mol Sci 2022; 23:621. [PMID: 35054802 PMCID: PMC8776158 DOI: 10.3390/ijms23020621] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 12/29/2021] [Accepted: 01/03/2022] [Indexed: 11/16/2022] Open
Abstract
The present benchmark calculations testify to the validity of time-dependent density functional theory (TD-DFT) when exploring the low-lying excited states potential energy surfaces of models of phenylalanine protein chains. Among three functionals suitable for systems exhibiting charge-transfer excited states, LC-ωPBE, CAM-B3LYP, and ωB97X-D, which were tested on a reference peptide system, we selected the ωB97X-D functional, which gave the best results compared to the approximate coupled-cluster singles and doubles (CC2) method. A quantitative agreement for both the geometrical parameters and the vibrational frequencies was obtained for the lowest singlet excited state (a ππ* state) of the series of capped peptides. In contrast, only a qualitative agreement was met for the corresponding adiabatic zero-point vibrational energy (ZPVE)-corrected excitation energies. Two composite protocols combining CC2 and DFT/TD-DFT methods were then developed to improve these calculations. Both protocols substantially reduced the error compared to CC2 and experiment, and the best of both even led to results of CC2 quality at a lower cost, thus providing a reliable alternative to this method for very large systems.
Collapse
Affiliation(s)
| | | | | | | | | | - Valérie Brenner
- LIDYL, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France; (M.L.); (T.V.); (E.G.); (B.T.); (M.M.)
| |
Collapse
|
24
|
Maitra NT. Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory. Annu Rev Phys Chem 2021; 73:117-140. [PMID: 34910562 DOI: 10.1146/annurev-physchem-082720-124933] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than otherwise possible. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. Here I review the fundamental challenges the approximate functionals have in describing double excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black-box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Neepa T Maitra
- Department of Physics, Rutgers University at Newark, Newark, New Jersey, USA;
| |
Collapse
|
25
|
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]
|
26
|
Prompt and Long-Lived Anti-Kasha Emission from Organic Dyes. Molecules 2021; 26:molecules26226999. [PMID: 34834093 PMCID: PMC8623836 DOI: 10.3390/molecules26226999] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 11/12/2021] [Accepted: 11/15/2021] [Indexed: 11/17/2022] Open
Abstract
Anti-Kasha behavior has been the subject of intense debate in the last few years, as demonstrated by the high number of papers appearing in the literature on this topic, dealing with both mechanistic and applicative aspects of this phenomenon. Examples of anomalous emitters reported in the last 10 years are collected in the present review, which is focused on strictly anti-Kasha organic molecules displaying radiative deactivation from Sn and/or Tn, with n greater than 1.
Collapse
|
27
|
Zhao R, Hettich CP, Chen X, Gao J. Minimal-active-space multistate density functional theory for excitation energy involving local and charge transfer states. NPJ COMPUTATIONAL MATERIALS 2021; 7:148. [PMID: 36713117 PMCID: PMC9881008 DOI: 10.1038/s41524-021-00624-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Accepted: 08/19/2021] [Indexed: 06/15/2023]
Abstract
Multistate density functional theory (MSDFT) employing a minimum active space (MAS) is presented to determine charge transfer (CT) and local excited states of bimolecular complexes. MSDFT is a hybrid wave function theory (WFT) and density functional theory, in which dynamic correlation is first incorporated in individual determinant configurations using a Kohn-Sham exchange-correlation functional. Then, nonorthogonal configuration-state interaction is performed to treat static correlation. Because molecular orbitals are optimized separately for each determinant by including Kohn-Sham dynamic correlation, a minimal number of configurations in the active space, essential to representing low-lying excited and CT states of interest, is sufficient to yield the adiabatic states. We found that the present MAS-MSDFT method provides a good description of covalent and CT excited states in comparison with experiments and high-level computational results. Because of the simplicity and interpretive capability through diabatic configuration weights, the method may be useful in dynamic simulations of CT and nonadiabatic processes.
Collapse
Affiliation(s)
- Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
- Institute of Theoretical Chemistry, Jilin University, Changchun, Jilin Province 130023, China
| | - Christian P. Hettich
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Xin Chen
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
- Beijing University Shenzhen Graduate School, Shenzhen 518055, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
- Beijing University Shenzhen Graduate School, Shenzhen 518055, China
| |
Collapse
|
28
|
Shepard C, Zhou R, Yost DC, Yao Y, Kanai Y. Simulating electronic excitation and dynamics with real-time propagation approach to TDDFT within plane-wave pseudopotential formulation. J Chem Phys 2021; 155:100901. [PMID: 34525811 DOI: 10.1063/5.0057587] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We give a perspective on simulating electronic excitation and dynamics using the real-time propagation approach to time-dependent density functional theory (RT-TDDFT) in the plane-wave pseudopotential formulation. RT-TDDFT is implemented in various numerical formalisms in recent years, and its practical application often dictates the most appropriate implementation of the theory. We discuss recent developments and challenges, emphasizing numerical aspects of studying real systems. Several applications of RT-TDDFT simulation are discussed to highlight how the approach is used to study interesting electronic excitation and dynamics phenomena in recent years.
Collapse
Affiliation(s)
- Christopher Shepard
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Ruiyi Zhou
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Dillon C Yost
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02142, USA
| | - Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, USA
| |
Collapse
|
29
|
Lacombe L, Maitra NT. Minimizing the Time-Dependent Density Functional Error in Ehrenfest Dynamics. J Phys Chem Lett 2021; 12:8554-8559. [PMID: 34464148 DOI: 10.1021/acs.jpclett.1c02020] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Simulating electron-ion dynamics using time-dependent density functional theory within an Ehrenfest dynamics scheme can be done in two ways that are in principle exact and identical: propagating time-dependent electronic Kohn-Sham equations or propagating electronic coefficients on surfaces obtained from linear-response. We show here that using an approximate functional leads to qualitatively different dynamics in the two approaches. We argue that the latter is more accurate because the functionals are evaluated on domains close to the ground state where currently used approximations perform better. We demonstrate this on an exactly solvable model of charge transfer and discuss implications for time-resolved spectroscopy.
Collapse
Affiliation(s)
- Lionel Lacombe
- Department of Physics, Rutgers University, Newark 07102, New Jersey United States
| | - Neepa T Maitra
- Department of Physics, Rutgers University, Newark 07102, New Jersey United States
| |
Collapse
|
30
|
Tao Z, Roy S, Schneider PE, Pavošević F, Hammes-Schiffer S. Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies. J Chem Theory Comput 2021; 17:5110-5122. [PMID: 34260237 DOI: 10.1021/acs.jctc.1c00454] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Saswata Roy
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Patrick E Schneider
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| |
Collapse
|
31
|
Loos PF, Comin M, Blase X, Jacquemin D. Reference Energies for Intramolecular Charge-Transfer Excitations. J Chem Theory Comput 2021; 17:3666-3686. [DOI: 10.1021/acs.jctc.1c00226] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, F-31400 Toulouse, France
| | | | - Xavier Blase
- Univ. Grenoble Alpes, CNRS, Inst NEEL, F-38042 Grenoble, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
32
|
Kocák J, Kraisler E, Schild A. Charge-Transfer Steps in Density Functional Theory from the Perspective of the Exact Electron Factorization. J Phys Chem Lett 2021; 12:3204-3209. [PMID: 33761257 DOI: 10.1021/acs.jpclett.1c00467] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
When a molecule dissociates, the exact Kohn-Sham (KS) and Pauli potentials may form step structures. Reproducing these steps correctly is central for the description of dissociation and charge-transfer processes in density functional theory (DFT): The steps align the KS eigenvalues of the dissociating subsystems relative to each other and determine where electrons localize. While the step height can be calculated from the asymptotic behavior of the KS orbitals, this provides limited insight into what causes the steps. We give an explanation of the steps with an exact mapping of the many-electron problem to a one-electron problem, the exact electron factorization (EEF). The potentials appearing in the EEF have a clear physical meaning that translates to the DFT potentials by replacing the interacting many-electron system with the KS system. With a simple model of a diatomic, we illustrate that the steps are a consequence of spatial electron entanglement and are the result of a charge transfer. From this mechanism, the step height can immediately be deduced. Moreover, two methods to approximately reproduce the potentials during dissociation are proposed. One is based on the states of the dissociated system, while the other one is based on an analogy to the Born-Oppenheimer treatment of a molecule. The latter method also shows that the steps connect adiabatic potential energy surfaces. The view of DFT from the EEF thus provides a better understanding of how many-electron effects are encoded in a one-electron theory and how they can be modeled.
Collapse
Affiliation(s)
- Jakub Kocák
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Eli Kraisler
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | - Axel Schild
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
33
|
Gould T, Kronik L. Ensemble generalized Kohn-Sham theory: The good, the bad, and the ugly. J Chem Phys 2021; 154:094125. [PMID: 33685152 DOI: 10.1063/5.0040447] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Two important extensions of Kohn-Sham (KS) theory are generalized: KS theory and ensemble KS theory. The former allows for non-multiplicative potential operators and greatly facilitates practical calculations with advanced, orbital-dependent functionals. The latter allows for quantum ensembles and enables the treatment of open systems and excited states. Here, we combine the two extensions, both formally and practically, first via an exact yet complicated formalism and then via a computationally tractable variant that involves a controlled approximation of ensemble "ghost interactions" by means of an iterative algorithm. The resulting formalism is illustrated using selected examples. This opens the door to the application of generalized KS theory in more challenging quantum scenarios and to the improvement of ensemble theories for the purpose of practical and accurate calculations.
Collapse
Affiliation(s)
- Tim Gould
- QLD Micro- and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia
| | - Leeor Kronik
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
34
|
Kraisler E, Hodgson MJP, Gross EKU. From Kohn-Sham to Many-Electron Energies via Step Structures in the Exchange-Correlation Potential. J Chem Theory Comput 2021; 17:1390-1407. [PMID: 33595312 PMCID: PMC8363072 DOI: 10.1021/acs.jctc.0c01093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Accurately
describing excited states within Kohn–Sham (KS)
density functional theory (DFT), particularly those which induce ionization
and charge transfer, remains a great challenge. Common exchange-correlation
(xc) approximations are unreliable for excited states owing, in part,
to the absence of a derivative discontinuity in the xc energy (Δ),
which relates a many-electron energy difference to the corresponding
KS energy difference. We demonstrate, analytically and numerically,
how the relationship between KS and many-electron energies leads to
the step structures observed in the exact xc potential in four scenarios:
electron addition, molecular dissociation, excitation of a finite
system, and charge transfer. We further show that steps in the potential
can be obtained also with common xc approximations, as simple as the
LDA, when addressed from the ensemble perspective. The article therefore
highlights how capturing the relationship between KS and many-electron
energies with advanced xc approximations is crucial for accurately
calculating excitations, as well as the ground-state density and energy
of systems which consist of distinct subsystems.
Collapse
Affiliation(s)
- Eli Kraisler
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| | - M J P Hodgson
- Department of Physics, Durham University, South Road, Durham DH1 3LE, United Kingdom.,Max-Planck-Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - E K U Gross
- Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, 9091401 Jerusalem, Israel
| |
Collapse
|
35
|
Gould T. Approximately Self-Consistent Ensemble Density Functional Theory: Toward Inclusion of All Correlations. J Phys Chem Lett 2020; 11:9907-9912. [PMID: 33170726 DOI: 10.1021/acs.jpclett.0c02894] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Recent theory developments in ensemble density functional theory (EDFT) promise to bring decades of work for ground states to the practical resolution of excited states, provided newly discovered "density-driven correlations" can be dealt with and adequate effective potentials can be found. This Letter introduces simple theories for both and shows that EDFT using these theories outperforms ΔSCF DFT and time-dependent DFT for low-lying gaps in most of the small atoms and molecules tested, even when all use the same density functional approximations. It thus establishes EDFT as a promising tool for low-cost studies of low-lying excited states and provides a clear route to practical EDFT implementation of arbitrary functional approximations.
Collapse
Affiliation(s)
- Tim Gould
- Qld Micro- and Nano-technology Centre, Griffith University, Nathan, Qld 4111, Australia
| |
Collapse
|
36
|
Conti I, Cerullo G, Nenov A, Garavelli M. Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going. J Am Chem Soc 2020; 142:16117-16139. [PMID: 32841559 PMCID: PMC7901644 DOI: 10.1021/jacs.0c04952] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Computational spectroscopy is becoming a mandatory tool for the interpretation of the
complex, and often congested, spectral maps delivered by modern non-linear multi-pulse
techniques. The fields of Electronic Structure Methods,
Non-Adiabatic Molecular Dynamics, and Theoretical
Spectroscopy represent the three pillars of the virtual ultrafast
optical spectrometer, able to deliver transient spectra in
silico from first principles. A successful simulation strategy requires a
synergistic approach that balances between the three fields, each one having its very
own challenges and bottlenecks. The aim of this Perspective is to demonstrate that,
despite these challenges, an impressive agreement between theory and experiment is
achievable now regarding the modeling of ultrafast photoinduced processes in complex
molecular architectures. Beyond that, some key recent developments in the three fields
are presented that we believe will have major impacts on spectroscopic simulations in
the very near future. Potential directions of development, pending challenges, and
rising opportunities are illustrated.
Collapse
Affiliation(s)
- Irene Conti
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Giulio Cerullo
- Dipartimento di Fisica, Politecnico di Milano, IFN-CNR, Piazza Leonardo da Vinci 32, I-20133 Milano, Italy
| | - Artur Nenov
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| | - Marco Garavelli
- Dipartimento di Chimica Industriale, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
| |
Collapse
|
37
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
38
|
Casanova-Páez M, Goerigk L. Assessing the Tamm–Dancoff approximation, singlet–singlet, and singlet–triplet excitations with the latest long-range corrected double-hybrid density functionals. J Chem Phys 2020; 153:064106. [DOI: 10.1063/5.0018354] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
39
|
Krumland J, Valencia AM, Pittalis S, Rozzi CA, Cocchi C. Understanding real-time time-dependent density-functional theory simulations of ultrafast laser-induced dynamics in organic molecules. J Chem Phys 2020; 153:054106. [PMID: 32770886 DOI: 10.1063/5.0008194] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Real-time time-dependent density functional theory, in conjunction with the Ehrenfest molecular dynamics scheme, is becoming a popular methodology to investigate ultrafast phenomena on the nanoscale. Thanks to recent developments, it is also possible to explicitly include in the simulations a time-dependent laser pulse, thereby accessing the transient excitation regime. However, the complexity entailed in these calculations calls for in-depth analysis of the accessible and yet approximate (either "dressed" or "bare") quantities in order to evaluate their ability to provide us with a realistic picture of the simulated processes. In this work, we analyze the ultrafast dynamics of three small molecules (ethylene, benzene, and thiophene) excited by a resonant laser pulse in the framework of the adiabatic local-density approximation. The electronic response to the laser perturbation in terms of induced dipole moment and excited-state population is compared to the results given by an exactly solvable two-level model. In this way, we can interpret the charge-carrier dynamics in terms of simple estimators, such as the number of excited electrons. From the computed transient absorption spectra, we unravel the appearance of nonlinear effects such as excited-state absorption and vibronic coupling. In this way, we observe that the laser excitation affects the vibrational spectrum by enhancing the anharmonicities therein, while the coherent vibrational motion contributes to stabilizing the electronic excitation already within a few tens of femtoseconds.
Collapse
Affiliation(s)
- Jannis Krumland
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, 12489 Berlin, Germany
| | - Ana M Valencia
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, 12489 Berlin, Germany
| | | | | | - Caterina Cocchi
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, 12489 Berlin, Germany
| |
Collapse
|
40
|
Holtrop F, Jupp AR, van Leest NP, Paradiz Dominguez M, Williams RM, Brouwer AM, de Bruin B, Ehlers AW, Slootweg JC. Photoinduced and Thermal Single-Electron Transfer to Generate Radicals from Frustrated Lewis Pairs. Chemistry 2020; 26:9005-9011. [PMID: 32259331 PMCID: PMC7496419 DOI: 10.1002/chem.202001494] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2020] [Revised: 04/06/2020] [Indexed: 12/17/2022]
Abstract
Archetypal phosphine/borane frustrated Lewis pairs (FLPs) are famed for their ability to activate small molecules. The mechanism is generally believed to involve two-electron processes. However, the detection of radical intermediates indicates that single-electron transfer (SET) generating frustrated radical pairs could also play an important role. These highly reactive radical species typically have significantly higher energy than the FLP, which prompted this investigation into their formation. Herein, we provide evidence that the classical phosphine/borane combinations PMes3 /B(C6 F5 )3 and PtBu3 /B(C6 F5 )3 both form an electron donor-acceptor (charge-transfer) complex that undergoes visible-light-induced SET to form the corresponding highly reactive radical-ion pairs. Subsequently, we show that by tuning the properties of the Lewis acid/base pair, the energy required for SET can be reduced to become thermally accessible.
Collapse
Affiliation(s)
- Flip Holtrop
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| | - Andrew R. Jupp
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| | - Nicolaas P. van Leest
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| | | | - René M. Williams
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| | - Albert M. Brouwer
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| | - Bas de Bruin
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| | - Andreas W. Ehlers
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
- Department of ChemistryScience FacultyUniversity of JohannesburgPO Box 254, Auckland ParkJohannesburgSouth Africa
| | - J. Chris Slootweg
- Van't Hoff Institute for Molecular SciencesUniversity of AmsterdamPO Box 941571090 GDAmsterdamThe Netherlands
| |
Collapse
|
41
|
Carter-Fenk K, Herbert JM. State-Targeted Energy Projection: A Simple and Robust Approach to Orbital Relaxation of Non-Aufbau Self-Consistent Field Solutions. J Chem Theory Comput 2020; 16:5067-5082. [DOI: 10.1021/acs.jctc.0c00502] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
42
|
Kozma B, Tajti A, Demoulin B, Izsák R, Nooijen M, Szalay PG. A New Benchmark Set for Excitation Energy of Charge Transfer States: Systematic Investigation of Coupled Cluster Type Methods. J Chem Theory Comput 2020; 16:4213-4225. [PMID: 32502351 PMCID: PMC7467641 DOI: 10.1021/acs.jctc.0c00154] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The
numerous existing publications on benchmarking quantum chemistry
methods for excited states rarely include Charge Transfer (CT) states,
although many interesting phenomena in, e.g., biochemistry and material
physics involve the transfer of electrons between fragments of the
system. Therefore, it is timely to test the accuracy of quantum chemical
methods for CT states, as well. In this study we first propose a new
benchmark set consisting of dimers having low-energy CT states. On
this set, the vertical excitation energy has been calculated with
Coupled Cluster methods including triple excitations (CC3, CCSDT-3,
CCSD(T)(a)*), as well as with methods including full or approximate
doubles (CCSD, STEOM-CCSD, CC2, ADC(2), EOM-CCSD(2)). The results
show that the popular CC2 and ADC(2) methods are much less accurate
for CT states than for valence states. On the other hand, EOM-CCSD
seems to have similar systematic overestimation of the excitation
energies for both types of states. Among the triples methods the novel
EOM-CCSD(T)(a)* method including noniterative triple excitations is
found to stand out with its consistently good performance for all
types of states, delivering essentially EOM-CCSDT quality results.
Collapse
Affiliation(s)
- Balázs Kozma
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| | - Attila Tajti
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| | - Baptiste Demoulin
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Péter G Szalay
- Institute of Chemistry, Laboratory of Theoretical Chemistry, ELTE Eötvös Loránd University, P.O. Box 32, H-1518, Budapest 112, Hungary
| |
Collapse
|
43
|
Ottochian A, Morgillo C, Ciofini I, Frisch MJ, Scalmani G, Adamo C. Double hybrids and time‐dependent density functional theory: An implementation and benchmark on charge transfer excited states. J Comput Chem 2020; 41:1242-1251. [DOI: 10.1002/jcc.26170] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2019] [Revised: 01/29/2020] [Accepted: 01/30/2020] [Indexed: 01/12/2023]
Affiliation(s)
- Alistar Ottochian
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences Paris France
| | - Carmela Morgillo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences Paris France
| | - Ilaria Ciofini
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences Paris France
| | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences Paris France
- Institut Universitaire de France Paris France
| |
Collapse
|
44
|
Callow TJ, Lathiotakis NN, Gidopoulos NI. Density-inversion method for the Kohn–Sham potential: Role of the screening density. J Chem Phys 2020; 152:164114. [DOI: 10.1063/5.0005781] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Timothy J. Callow
- Department of Physics, Durham University, South Road, Durham DH1 3LE, United Kingdom
- Max-Planck-Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - Nektarios N. Lathiotakis
- Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, Vass. Constantinou 48, 116 35 Athens, Greece
| | - Nikitas I. Gidopoulos
- Department of Physics, Durham University, South Road, Durham DH1 3LE, United Kingdom
| |
Collapse
|
45
|
Ramirez FF, Bustamente CM, González Lebrero MC, Scherlis DA. Transport and Spectroscopy in Conjugated Molecules: Two Properties and a Single Rationale. J Chem Theory Comput 2020; 16:2930-2940. [DOI: 10.1021/acs.jctc.9b01122] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Francisco F. Ramirez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| | - Carlos M. Bustamente
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, Buenos Aires C1428EHA, Argentina
| |
Collapse
|
46
|
Narsaria AK, Ruijter JD, Hamlin TA, Ehlers AW, Guerra CF, Lammertsma K, Bickelhaupt FM. Performance of TDDFT Vertical Excitation Energies of Core-Substituted Naphthalene Diimides. J Comput Chem 2020; 41:1448-1455. [PMID: 32142173 PMCID: PMC7317478 DOI: 10.1002/jcc.26188] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Revised: 02/10/2020] [Accepted: 02/22/2020] [Indexed: 01/15/2023]
Abstract
We have evaluated the performance of various density functionals, covering generalized gradient approximation (GGA), global hybrid (GH) and range‐separated hybrid (RSH), using time dependent density functional theory (TDDFT) for computing vertical excitation energies against experimental absorption maximum (λmax) for a set of 10 different core‐substituted naphthalene diimides (cNDI) recorded in dichloromethane. The computed excitation in case of GH PBE0 is most accurate while the trend is most systematic with RSH LCY‐BLYP compared to λmax. We highlight the importance of including solvent effects for optimal agreement with the λmax. Increasing the basis set size from TZ2P to QZ4P has a negligible influence on the computed excitation energies. Notably, RSH CAMY‐B3LYP gave the least error for charge‐transfer excitation. The poorest agreement with λmax is obtained with semi‐local GGA functionals. Use of the optimally‐tuned RSH LCY‐BLYP* is not recommended because of the high computational cost and marginal improvement in results.
Collapse
Affiliation(s)
- Ayush K Narsaria
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Julian D Ruijter
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Trevor A Hamlin
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Andreas W Ehlers
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands.,Department of Chemistry, University of Johannesburg, Johannesburg, South Africa
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, Leiden, The Netherlands
| | - Koop Lammertsma
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Department of Chemistry, University of Johannesburg, Johannesburg, South Africa
| | - F Matthias Bickelhaupt
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling (ACMM), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands.,Institute of Molecules and Materials, Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
47
|
Gould T, Pittalis S. Density-Driven Correlations in Ensemble Density Functional Theory: Insights from Simple Excitations in Atoms. Aust J Chem 2020. [DOI: 10.1071/ch19504] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Ensemble density functional theory extends the usual Kohn-Sham machinery to quantum state ensembles involving ground- and excited states. Recent work by the authors [Phys. Rev. Lett. 119, 243001 (2017); 123, 016401 (2019)] has shown that both the Hartree-exchange and correlation energies can attain unusual features in ensembles. Density-driven (DD) correlations – which account for the fact that pure-state densities in Kohn-Sham ensembles do not necessarily reproduce those of interacting pure states – are one such feature. Here we study atoms (specifically S–P and S–S transitions) and show that the magnitude and behaviour of DD correlations can vary greatly with the variation of the orbital angular momentum of the involved states. Such estimations are obtained through an approximation for DD correlations built from relevant exact conditions, Kohn-Sham inversion, and plausible assumptions for weakly correlated systems.
Collapse
|
48
|
Mandal I, Manna S, Venkatramani R. UV-Visible Lysine-Glutamate Dimer Excitations in Protein Charge Transfer Spectra: TDDFT Descriptions Using an Optimally Tuned CAM-B3LYP Functional. J Phys Chem B 2019; 123:10967-10979. [PMID: 31763835 DOI: 10.1021/acs.jpcb.9b07827] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Recent reports of distinctive UV-vis absorption profiles for monomeric proteins rich in charged amino acids that span 250-800 nm have opened up a new label-free optical spectral window for probing biomolecular structure and interactions. Combined experimental-computational studies have revealed that such broad absorption profiles of these proteins arise from photoexcited charge transfer (CT) transitions in spatially proximal charged amino acids such as lysine (Lys) and glutamate (Glu). Here, using time-dependent density functional theory (TDDFT) with an optimally tuned CAM-B3LYP functional, we refine the computed UV-vis spectra for Lys-Glu dimers within protein folds and quantify the percentage CT character of the constituent transitions. The optimally tuned functionals are derived through a careful analysis of the CAM-B3LYP parameter space for Lys-Glu dimers as a function of amino-acid conformation and side chain separation. Our studies reveal that the tuned Lys-Glu dimer spectrum spans 150-650 nm and possesses 5 specific types of CT excitations with diverse and large spatial charge separation length scales of 2-10 Å. These include inter-/intra-residue peptide backbone to peptide backbone (BB-CT) excitations spanning 160-210 nm, inter-/intra-residue peptide backbone to side chain (BS-CT) excitations spanning 160-260 nm, and side chain to side chain (SS-CT) excitations, which show the broadest absorption range spanning 260-650 nm.
Collapse
Affiliation(s)
- Imon Mandal
- Department of Chemical Sciences , Tata Institute of Fundamental Research , Homi Bhabha Road, Colaba , Mumbai 400 005 , India
| | - Sutanuka Manna
- Department of Chemical Sciences , Tata Institute of Fundamental Research , Homi Bhabha Road, Colaba , Mumbai 400 005 , India
| | - Ravindra Venkatramani
- Department of Chemical Sciences , Tata Institute of Fundamental Research , Homi Bhabha Road, Colaba , Mumbai 400 005 , India
| |
Collapse
|
49
|
Zhao L, Neuscamman E. Density Functional Extension to Excited-State Mean-Field Theory. J Chem Theory Comput 2019; 16:164-178. [DOI: 10.1021/acs.jctc.9b00530] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Luning Zhao
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
50
|
Partovi-Azar P, Sarap CS, Fyta M. In silico Complexes of Amino Acids and Diamondoids. Chemphyschem 2019; 20:2166-2170. [PMID: 31287610 DOI: 10.1002/cphc.201900394] [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: 04/17/2019] [Revised: 07/01/2019] [Indexed: 11/11/2022]
Abstract
We report on the specific interaction of a small diamond-like molecule, known as diamondoid, with single amino-acids forming nano/bio molecular complexes. Using time-dependent density-functional theory calculations we have studied two different relative configurations of three prototypical amino acids, phenylalanine, tyrosine, and tryptophan, with the diamondoid. The optical and charge-transfer properties of these complexes exhibit amino acid and topology specific features which can be directly utilized for in the direction of novel biomolecule detection schemes.
Collapse
Affiliation(s)
- Pouya Partovi-Azar
- Institute of Chemistry, Martin Luther University Halle-Wittenberg, Von-Danckelmann-Platz 4, 06120, Halle (Saale), Germany
| | - Chandra Shekar Sarap
- Institute for Computational Physics, Universität Stuttgart, Allmandring 3, 70569, Stuttgart, Germany
| | - Maria Fyta
- Institute for Computational Physics, Universität Stuttgart, Allmandring 3, 70569, Stuttgart, Germany
| |
Collapse
|