1
|
Garcia-Alvarez JC, Gozem S. Absorption Intensities of Organic Molecules from Electronic Structure Calculations versus Experiments: the Effect of Solvation, Method, Basis Set, and Transition Moment Gauge. J Chem Theory Comput 2024; 20. [PMID: 39141425 PMCID: PMC11360136 DOI: 10.1021/acs.jctc.4c00642] [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/15/2024] [Revised: 07/24/2024] [Accepted: 07/25/2024] [Indexed: 08/15/2024]
Abstract
Recently, we derived experimental oscillator strengths (OSs) from well-defined UV-visible absorption spectral peaks of 100 molecules in solution. Here, we focus on a subset of transitions with the highest reliability to further benchmark the OSs from several wave function methods and density functionals. We consider multiple basis sets, transition moment gauges (length, velocity, and mixed), and solvent corrections. Most transitions in the comparison set come from conjugated molecules and have π → π* character. We use an automated algorithm to assign computed transitions to experimental bands. OSs computed using the Tamm-Dancoff approximation (TDA), CIS, or EOM-CCSD exhibited a strong gauge dependence, which is diminished in linear response theories (TD-DFT, TD-HF, and to a smaller degree LR-CCSD). OSs calculated from TD-DFT with PCM solvent models are systematically larger than apparent OSs derived from experimental spectra. For example, fcomp from hybrid functionals and PCM have mean absolute errors that are ∼10% of n·fexp, where n is a solvent refractive index factor that arises from the energy flux of the radiation field in a dielectric (solvent). Theoretical cavity field corrections considering spherical cavities do not improve the agreement between computed and experimental data. Corrections that account for the molecular shape and the direction of transition dipole moments, or that explicitly account for the effect of solvent molecules on the local field, should be more appropriate.
Collapse
Affiliation(s)
| | - Samer Gozem
- Department of Chemistry, Georgia
State University, Atlanta, Georgia 30302, United States
| |
Collapse
|
2
|
Holzer C, Franzke YJ. Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory. Chemphyschem 2024; 25:e202400120. [PMID: 38456204 DOI: 10.1002/cphc.202400120] [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: 02/02/2024] [Revised: 03/07/2024] [Accepted: 03/08/2024] [Indexed: 03/09/2024]
Abstract
Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.
Collapse
Affiliation(s)
- Christof Holzer
- Karlsruhe Institute of Technology (KIT), Institute of Theoretical Solid State Physics, Kaiserstraße 12, 76131, Karlsruhe, Germany
| | - Yannick J Franzke
- Friedrich Schiller University Jena, Otto Schott Institute of Materials Research, Löbdergraben 32, 07743, Jena, Germany
| |
Collapse
|
3
|
Rauwolf N, Klopper W, Holzer C. Non-linear light-matter interactions from the Bethe-Salpeter equation. J Chem Phys 2024; 160:061101. [PMID: 38341783 DOI: 10.1063/5.0191499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Accepted: 01/15/2024] [Indexed: 02/13/2024] Open
Abstract
A route to assess non-linear light-matter interactions from the increasingly popular GW-Bethe-Salpeter equation (GW-BSE) method is outlined. In the present work, the necessary analytic expressions within the static-screened exchange approximation of the BSE are derived. This enables a straightforward implementation of the computation of the first hyperpolarizability as well as two-photon absorption processes for molecular systems. Benchmark calculations on small molecular systems reveal that the GW-BSE method is intriguingly accurate for predicting both first hyperpolarizabilities and two-photon absorption strengths. Using state-of-the-art Kohn-Sham references as a starting point, the accuracy of the GW-BSE method rivals that of the coupled-cluster singles-and-doubles method, outperforming both second-order coupled-cluster and time-dependent density-functional theory.
Collapse
Affiliation(s)
- Nina Rauwolf
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology, Kaiserstraße 12, 76131 Karlsruhe, Germany
| |
Collapse
|
4
|
Dall’Osto G, Marsili M, Vanzan M, Toffoli D, Stener M, Corni S, Coccia E. Peeking into the Femtosecond Hot-Carrier Dynamics Reveals Unexpected Mechanisms in Plasmonic Photocatalysis. J Am Chem Soc 2024; 146:2208-2218. [PMID: 38199967 PMCID: PMC10811681 DOI: 10.1021/jacs.3c12470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Revised: 12/23/2023] [Accepted: 12/29/2023] [Indexed: 01/12/2024]
Abstract
Plasmonic-driven photocatalysis may lead to reaction selectivity that cannot be otherwise achieved. A fundamental role is played by hot carriers, i.e., electrons and holes generated upon plasmonic decay within the metal nanostructure interacting with molecular species. Understanding the elusive microscopic mechanism behind such selectivity is a key step in the rational design of hot-carrier reactions. To accomplish that, we present state-of-the-art multiscale simulations, going beyond density functional theory, of hot-carrier injections for the rate-determining step of a photocatalytic reaction. We focus on carbon dioxide reduction, for which it was experimentally shown that the presence of a rhodium nanocube under illumination leads to the selective production of methane against carbon monoxide. We show that selectivity is due to a (predominantly) direct hole injection from rhodium to the reaction intermediate CHO. Unexpectedly, such an injection does not promote the selective reaction path by favoring proper bond breaking but rather by promoting bonding of the proper molecular fragment to the surface.
Collapse
Affiliation(s)
- Giulia Dall’Osto
- Dipartimento
di Scienze Chimiche, Università di
Padova, via F. Marzolo 1, 35131 Padova, Italy
| | - Margherita Marsili
- Dipartimento
di Fisica e Astronomia “Augusto Righi”, University of Bologna, Viale Berti Pichat 6/2, 40127 Bologna, Italy
| | - Mirko Vanzan
- Dipartimento
di Scienze Chimiche, Università di
Padova, via F. Marzolo 1, 35131 Padova, Italy
- Dipartimento
di Fisica, University of Milan, Via Giovanni Celoria 16, 20133 Milano, Italy
| | - Daniele Toffoli
- Dipartimento
di Scienze Chimiche e Farmaceutiche, University
of Trieste, via L. Giorgieri 1, 34127 Trieste, Italy
| | - Mauro Stener
- Dipartimento
di Scienze Chimiche e Farmaceutiche, University
of Trieste, via L. Giorgieri 1, 34127 Trieste, Italy
| | - Stefano Corni
- Dipartimento
di Scienze Chimiche, Università di
Padova, via F. Marzolo 1, 35131 Padova, Italy
- Istituto
Nanoscienze-CNR, via
Campi 213/A, 41125 Modena, Italy
| | - Emanuele Coccia
- Dipartimento
di Scienze Chimiche e Farmaceutiche, University
of Trieste, via L. Giorgieri 1, 34127 Trieste, Italy
| |
Collapse
|
5
|
Soltani Nejad M, Alipour M. How does theory compare to experiment for oscillator strengths in electronic spectra? Proposing range-separated hybrids with reliable accountability. Phys Chem Chem Phys 2024; 26:879-894. [PMID: 38087910 DOI: 10.1039/d3cp04793h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2024]
Abstract
As an important quantity in atomic and molecular spectroscopy, oscillator strength should be mentioned. Oscillator strength is linked to the transition dipole moment and consequently to the transition probability between two states, where its magnitude is directly connected to the intensity of the peaks in ultraviolet-visible spectra. However, accurately accounting for oscillator strengths still remains one of the greatest challenges in theory and experiment. Given previous efforts in the context of investigations into oscillator strengths, the related theoretical treatments are relatively limited and have proven to be challenging. In this work, the oscillator strengths in the electronic spectra of organic compounds have thoroughly been investigated with the help of optimally tuned range-separated hybrids (OT-RSHs). In particular, variants of the OT-RSHs combined with the polarizable continuum model (PCM), OT-RSHs-PCM, as well as their screened versions accounting for the screening effects by the electron correlation through the dielectric constant, OT-SRSHs-PCM, are proposed for reliable prediction of the oscillator strengths. The role of the involved ingredients in the proposed methods, namely the underlying density functional approximations, short-range and long-range Hartree-Fock (HF) exchange, as well as the range-separation parameter, has been examined in detail. It is shown that any combination of the parameters in the proposed approximations does not render the reliable oscillator strengths, but a particular compromise among them is needed to describe the experimental data well. Perusing all the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM, coupled with the linear response theory in the non-equilibrium solvation regime, with the correct asymptotic behavior and incorporating no (low) HF exchange contributions in the short-range part. The best proposed approximations also reveal superior performances not only with respect to their standard counterparts with the default parameters but also as compared to earlier range-separated functionals. Finally, the applicability of the best approximation is also put into broader perspective, where it is used for predicting the oscillator strengths in other sets of compounds not included in the process of developing the approximations. Hopefully, our proposed method can function as an affordable alternative to the expensive wave function-based methods for both theoretical modeling and confirming the experimental observations in the field of electronic spectroscopy.
Collapse
Affiliation(s)
- Mahdi Soltani Nejad
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran.
| | - Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran.
| |
Collapse
|
6
|
Knysh I, Villalobos-Castro JDJ, Duchemin I, Blase X, Jacquemin D. Excess and excited-state dipole moments of real-life dyes: a comparison between wave-function, BSE/ GW, and TD-DFT values. Phys Chem Chem Phys 2023; 25:29993-30004. [PMID: 37905396 DOI: 10.1039/d3cp04467j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
In this work, we assess the accuracy of the Bethe-Salpeter equation (BSE) many-body Green's function formalism, adopting the eigenvalue-self-consistent evGW exchange-correlation kernel, for the calculation of the excited-state (μES) and excess dipole moments (Δμ), the latter ones being the changes of dipole amplitude between the ground and excited states (ES), in organic dyes. We compare the results obtained with wave-function methods [ADC(2), CC2, and CCSD], time-dependent density functional theory (TD-DFT), and BSE/evGW levels of theory. First, we compute the evolution of the dipole moments of the two lowest singlet excited states of 4-(dimethylamino)benzonitrile (DMABN) upon twisting of the amino group. Next, we use a set of 25 dyes having ES characters ranging from locally excited to charge transfer to determine both μES and Δμ. For DMABN our results show that BSE/evGW provides Δμ values closer to the CCSD reference and more consistent trends than TD-DFT. Moreover, a statistical analysis of both Δμ and μES for the set of 25 dyes shows that the BSE/evGW accuracy is comparable or sometimes slightly better than that of TD-M06-2X and TD-CAM-B3LYP, BSE/evGW outperforming TD-DFT in challenging cases (zwitterionic and cyanine transitions). Finally, the starting point dependency of BSE/evGW seems to be larger for Δμ, ES dipoles, and oscillator strengths than for transition energies.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France.
| | | | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L_Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut Néel, F-38042 Grenoble, France.
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France.
- Institut Universitaire de France, F-75005 Paris, France
| |
Collapse
|
7
|
Knysh I, Villalobos-Castro JDJ, Duchemin I, Blase X, Jacquemin D. Exploring Bethe-Salpeter Excited-State Dipoles: The Challenging Case of Increasingly Long Push-Pull Oligomers. J Phys Chem Lett 2023; 14:3727-3734. [PMID: 37042642 DOI: 10.1021/acs.jpclett.3c00699] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The change of molecular dipole moment induced by photon absorption is key to interpret the measured optical spectra. Except for compact molecules, time-dependent density functional theory (TD-DFT) remains the only theory allowing to quickly predict excited-state dipoles (μES), albeit with a strong dependency on the selected exchange-correlation functional. This Letter presents the first assessment of the performances of the many-body Green's function Bethe-Salpeter equation (BSE) formalism for the evaluation of the μES. We explore increasingly long push-pull oligomers as they present an excited-state nature evolving with system size. This work shows that BSE's μES do present the same evolution with oligomeric length as their CC2 and CCSD counterparts, with a dependency on the starting exchange-correlation functional that is strongly decreased as compared to TD-DFT. This Letter demonstrates that BSE is a valuable alternative to TD-DFT for properties related to the excited-state density and not only for transition energies and oscillator strengths.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | | | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
8
|
Knysh I, Letellier K, Duchemin I, Blase X, Jacquemin D. Excited state potential energy surfaces of N-phenylpyrrole upon twisting: reference values and comparison between BSE/ GW and TD-DFT. Phys Chem Chem Phys 2023; 25:8376-8385. [PMID: 36883347 DOI: 10.1039/d3cp00474k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
Abstract
The puzzling case of the mixing between the charge transfer (CT) and local excited (LE) characters upon twisting of the geometry of N-phenylpyrrole (N-PP) is investigated considering the six low-lying singlet excited states (ES). The theoretical calculations of the potential energy surfaces (PES) have been performed for these states using a Coupled Cluster method accounting for the impact of the contributions from the triples, many-body Green's function GW and Bethe-Salpeter equation (BSE) formalisms, as well as Time-Dependent Density Functional Theory (TD-DFT) using various exchange-correlation functionals. Our findings confirm that the BSE formalism is more reliable than TD-DFT for close-lying ES with mixed CT/LE nature. More specifically, BSE/GW yields a more accurate evolution of the excited state PES than TD-DFT when compared to the reference coupled cluster values. BSE/GW PES curves also show negligible exchange-correlation functional starting point dependency in sharp contrast with their TD-DFT counterparts.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000, Nantes, France.
| | | | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L_Sim, 38054, Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut Néel, F-38042, Grenoble, France.
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000, Nantes, France.
- Institut Universitaire de France, 75005, Paris, France
| |
Collapse
|
9
|
Damour Y, Quintero-Monsebaiz R, Caffarel M, Jacquemin D, Kossoski F, Scemama A, Loos PF. Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality. J Chem Theory Comput 2023; 19:221-234. [PMID: 36548519 DOI: 10.1021/acs.jctc.2c01111] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We report ground- and excited-state dipole moments and oscillator strengths (computed in different "gauges" or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as Configuration Interaction using a Perturbative Selection made Iteratively (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CC with singles and doubles (CCSD) to CC with singles, doubles, and triples (CCSDT) or from CCSDT to CC with singles, doubles, triples, and quadruples (CCSDTQ)) reduces the average error with respect to the near-FCI reference values by approximately 1 order of magnitude.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Raúl Quintero-Monsebaiz
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Michel Caffarel
- 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.,Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- 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
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| |
Collapse
|
10
|
McKeon CA, Hamed SM, Bruneval F, Neaton JB. An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules. J Chem Phys 2022; 157:074103. [DOI: 10.1063/5.0097582] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.
Collapse
Affiliation(s)
- Caroline A. McKeon
- Department of Physics, University of California, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Natural Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Samia M. Hamed
- Department of Physics, University of California, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Natural Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Fabien Bruneval
- Université Paris-Saclay, CEA, Service de Recherches de Métallurgie Physique, 91191 Gif-sur-Yvette, France
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Natural Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Kavli ENSI, University of California, Berkeley, California 94720, USA
| |
Collapse
|
11
|
Marsili M, Corni S. Electronic Dynamics of a Molecular System Coupled to a Plasmonic Nanoparticle Combining the Polarizable Continuum Model and Many-Body Perturbation Theory. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:8768-8776. [PMID: 35655939 PMCID: PMC9150096 DOI: 10.1021/acs.jpcc.2c02209] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Revised: 04/28/2022] [Indexed: 06/15/2023]
Abstract
The efficiency of plasmonic metallic nanoparticles in harvesting and concentrating light energy in their proximity triggers a wealth of important and intriguing phenomena. For example, spectroscopies are able to reach single-molecule and intramolecule sensitivities, and important chemical reactions can be effectively photocatalyzed. For the real-time description of the coupled dynamics of a molecule's electronic system and of a plasmonic nanoparticle, a methodology has been recently proposed (J. Phys. Chem. C. 120, 2016, 28774-28781) which combines the classical description of the nanoparticle as a polarizable continuum medium with a quantum-mechanical description of the molecule treated at the time-dependent configuration interaction (TDCI) level. In this work, we extend this methodology by describing the molecule using many-body perturbation theory: the molecule's excitation energies, transition dipoles, and potentials computed at the GW/Bethe-Salpeter equation (BSE) level. This allows us to overcome current limitations of TDCI in terms of achievable accuracy without compromising on the accessible molecular sizes. We illustrate the developed scheme by characterizing the coupled nanoparticle/molecule dynamics of two prototype molecules, LiCN and p-nitroaniline.
Collapse
Affiliation(s)
- Margherita Marsili
- Dipartimento
di Science Chimiche, Università di
Padova, via F. Marzolo 1, I-35131, Padova, Italy
| | - Stefano Corni
- Dipartimento
di Science Chimiche, Università di
Padova, via F. Marzolo 1, I-35131, Padova, Italy
- CNR
Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| |
Collapse
|
12
|
Yao Y, Golze D, Rinke P, Blum V, Kanai Y. All-Electron BSE@ GW Method for K-Edge Core Electron Excitation Energies. J Chem Theory Comput 2022; 18:1569-1583. [PMID: 35138865 DOI: 10.1021/acs.jctc.1c01180] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We present an accurate computational approach to calculate absolute K-edge core electron excitation energies as measured by X-ray absorption spectroscopy. Our approach employs an all-electron Bethe-Salpeter equation (BSE) formalism based on GW quasiparticle energies (BSE@GW) using numeric atom-centered orbitals (NAOs). The BSE@GW method has become an increasingly popular method for the computation of neutral valence excitation energies of molecules. However, it was so far not applied to molecular K-edge excitation energies. We discuss the influence of different numerical approximations on the BSE@GW calculation and employ in our final setup (i) exact numeric algorithms for the frequency integration of the GW self-energy, (ii) G0W0 and BSE starting points with ∼50% of exact exchange, (iii) the Tamm-Dancoff approximation and (iv) relativistic corrections. We study the basis set dependence and convergence with common Gaussian-type orbital and NAO basis sets. We identify the importance of additional spatially confined basis functions as well as of diffuse augmenting basis functions. The accuracy of our BSE@GW method is assessed for a benchmark set of small organic molecules, previously used for benchmarking the equation-of-motion coupled cluster method [Peng et al., J. Chem. Theory Comput., 2015, 11, 4146], as well as the medium-sized dibenzothiophene (DBT) molecule. Our BSE@GW results for absolute excitation energies are in excellent agreement with the experiment, with a mean average error of only 0.63 eV for the benchmark set and with errors <1 eV for the DBT molecule.
Collapse
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Dorothea Golze
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01062 Dresden, Germany.,Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Finland
| | - Patrick Rinke
- Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Finland
| | | | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
13
|
Bintrim SJ, Berkelbach TC. Full-frequency dynamical Bethe-Salpeter equation without frequency and a study of double excitations. J Chem Phys 2022; 156:044114. [PMID: 35105075 PMCID: PMC8807000 DOI: 10.1063/5.0074434] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Accepted: 12/31/2021] [Indexed: 01/30/2023] Open
Abstract
The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is O(N6), where N is the system size. To avoid the common static screening approximation, we show that the full-frequency dynamical BSE can be exactly reformulated as a frequency-independent eigenvalue problem in an expanded space of single and double excitations. When combined with an iterative eigensolver and the density fitting approximation to the electron repulsion integrals, this reformulation yields a dynamical BSE algorithm whose computational time is O(N5), which we verify numerically. Furthermore, the reformulation provides direct access to excited states with dominant double excitation character, which are completely absent in the spectrum of the statically screened BSE. We study the 21Ag state of butadiene, hexatriene, and octatetraene and find that GW/BSE overestimates the excitation energy by about 1.5-2 eV and significantly underestimates the double excitation character.
Collapse
Affiliation(s)
- Sylvia J. Bintrim
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
14
|
Tarleton AS, Garcia-Alvarez JC, Wynn A, Awbrey CM, Roberts TP, Gozem S. OS100: A Benchmark Set of 100 Digitized UV-Visible Spectra and Derived Experimental Oscillator Strengths. J Phys Chem A 2022; 126:435-443. [PMID: 35015532 DOI: 10.1021/acs.jpca.1c08988] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Excited-state quantum chemical calculations usually report excitation energies and oscillator strengths, f, for each electronic transition. On the other hand, UV-visible spectrophotometric experiments measure energy-dependent molar extinction/attenuation coefficients, ε(v), that give absorption band line shapes when plotted. ε(v) and f are related, but this relation is complicated by broadening and solvation effects. We fitted and integrated 100 experimental UV-visible spectra to obtain 164 fexp values for absorption bands appearing in these spectra. The 100 UV-visible spectra belong to solvated organic molecules ranging in size from 6-34 atoms. We estimated uncertainties in the fitting to indicate confidence level in the reported fexp values. The corresponding computed oscillator strengths (fcomp) were obtained with time-dependent density functional theory and a polarizable continuum solvent model. By expressing experimental and computed absorption strengths using a common quantity, we directly compared fcomp and fexp. Although fcomp and fexp are well correlated (linear regression R2 = 0.921), fcomp in most cases overestimated fexp (regression slope = 1.34). The agreement between absolute fcomp and fexp values was substantially improved by accounting for a solvent refractive index factor, as suggested in some derivations in the literature. The 100 digitized UV-visible spectra are included as plain text files in the Supporting Information to aid in benchmarking computational or machine learning methods that aim to simulate realistic UV-visible absorption spectra.
Collapse
Affiliation(s)
- Astrid S Tarleton
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Jorge C Garcia-Alvarez
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Anah Wynn
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Cade M Awbrey
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Tomas P Roberts
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| | - Samer Gozem
- Department of Chemistry, Georgia State University, Atlanta, Georgia 30302, United States
| |
Collapse
|
15
|
Franzke YJ, Holzer C, Mack F. NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation. J Chem Theory Comput 2022; 18:1030-1045. [PMID: 34981925 DOI: 10.1021/acs.jctc.1c00999] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
16
|
Wang X, Berkelbach TC. Absorption Spectra of Solids from Periodic Equation-of-Motion Coupled-Cluster Theory. J Chem Theory Comput 2021; 17:6387-6394. [PMID: 34559525 DOI: 10.1021/acs.jctc.1c00692] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present ab initio absorption spectra of six three-dimensional semiconductors and insulators calculated using Gaussian-based periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). The spectra are calculated efficiently by solving a system of linear equations at each frequency, giving access to an energy range of tens of electronvolts without explicit enumeration of excited states. We assess the impact of cost-saving approximations associated with Brillouin zone sampling, frozen orbitals, and the partitioned EOM-CCSD approximation. Although our most converged spectra exhibit line shapes that are in good agreement with experimental spectra, they are uniformly shifted to higher energies by about 1 eV, which is not explained by the remaining finite-size errors. We tentatively attribute this discrepancy to a combination of vibrational effects and the remaining electron correlation, i.e., triple excitations and above.
Collapse
Affiliation(s)
- Xiao Wang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Timothy C Berkelbach
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States.,Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
17
|
Grobas Illobre P, Marsili M, Corni S, Stener M, Toffoli D, Coccia E. Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms. J Chem Theory Comput 2021; 17:6314-6329. [PMID: 34486881 PMCID: PMC8515806 DOI: 10.1021/acs.jctc.1c00211] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2021] [Indexed: 12/16/2022]
Abstract
In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.
Collapse
Affiliation(s)
- P. Grobas Illobre
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - M. Marsili
- Dipartimento
di Scienze Chimiche, Universitá di
Padova, via Marzolo 1, Padova 35131, Italy
| | - S. Corni
- Dipartimento
di Scienze Chimiche, Universitá di
Padova, via Marzolo 1, Padova 35131, Italy
- CNR
Istituto di Nanoscienze, via Campi 213/A, Modena 41125, Italy
| | - M. Stener
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - D. Toffoli
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| | - E. Coccia
- Dipartimento
di Scienze Chimiche e Farmaceutiche, Universitá
di Trieste, via L. Giorgieri 1, Trieste 34127, Italy
| |
Collapse
|
18
|
Sarkar R, Boggio-Pasqua M, Loos PF, Jacquemin D. Benchmarking TD-DFT and Wave Function Methods for Oscillator Strengths and Excited-State Dipole Moments. J Chem Theory Comput 2021; 17:1117-1132. [DOI: 10.1021/acs.jctc.0c01228] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
19
|
Chrayteh A, Blondel A, Loos PF, Jacquemin D. Mountaineering Strategy to Excited States: Highly Accurate Oscillator Strengths and Dipole Moments of Small Molecules. J Chem Theory Comput 2020; 17:416-438. [DOI: 10.1021/acs.jctc.0c01111] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Amara Chrayteh
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Aymeric Blondel
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques, Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| |
Collapse
|
20
|
Kehry M, Franzke YJ, Holzer C, Klopper W. Quasirelativistic two-component core excitations and polarisabilities from a damped-response formulation of the Bethe–Salpeter equation. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1755064] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Max Kehry
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Yannick J. Franzke
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Christof Holzer
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Wim Klopper
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| |
Collapse
|
21
|
Cazzaniga M, Cargnoni F, Penconi M, Bossi A, Ceresoli D. Ab Initio Many-Body Perturbation Theory Calculations of the Electronic and Optical Properties of Cyclometalated Ir(III) Complexes. J Chem Theory Comput 2020; 16:1188-1199. [PMID: 31860292 DOI: 10.1021/acs.jctc.9b00763] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Cyclometalated Ir(III) compounds are the preferred choice as organic emitters in organic light-emitting diodes. In practice, the presence of the transition metal surrounded by carefully designed ligands allows fine-tuning of the emission frequency as well as good efficiency of the device. To support the development of new compounds, experimental measurements are generally compared with absorption and emission spectra obtained from ab initio calculations. The standard approach for these calculations is time-dependent density functional theory (TDDFT) with a hybrid exchange-correlation functional like B3LYP. Because of the size of these compounds, the application of more complex quantum chemistry approaches can be challenging. In this work, we used many-body perturbation theory approaches, in particular the GW approximation with the Bethe-Salpeter equation (BSE) implemented in Gaussian basis sets, to calculate the quasiparticle properties and the absorption spectra of six cyclometalated Ir(III) complexes, going beyond TDDFT. In the presented results, we compared standard TDDFT simulations with BSE calculations performed on top of perturbative G0W0 and accounting for eigenvalue self-consistency. Moreover, in order to investigate in detail the effect of the DFT starting point, we concentrated on Ir(ppy)3 and performed GW-BSE simulations starting from different DFT exchange-correlation potentials.
Collapse
Affiliation(s)
- Marco Cazzaniga
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Fausto Cargnoni
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Marta Penconi
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Alberto Bossi
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Davide Ceresoli
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| |
Collapse
|
22
|
Liu C, Kloppenburg J, Yao Y, Ren X, Appel H, Kanai Y, Blum V. All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals. J Chem Phys 2020; 152:044105. [DOI: 10.1063/1.5123290] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chi Liu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Jan Kloppenburg
- Institute of Condensed Matter and Nanoscience, Université Catholique de Louvain, Louvain-la-Neuve 1348, Belgium
| | - Yi Yao
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, USA
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| | - Xinguo Ren
- CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Heiko Appel
- Max Planck Institute for the Structure and Dynamics of Matter, Center for Free Electron Laser Science, 22761 Hamburg, Germany
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - Volker Blum
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
23
|
Abstract
The accuracy of three different complete active space (CAS) self-consistent field (CASSCF) methods is investigated for the electronically excited-state benchmark set of Schreiber , M. ; et al. J. Chem. Phys. 2008 , 128 , 134110 . Comparison of the CASSCF linear response (LR) methods MC-RPA and MC-TDA and the state-averaged (SA) CASSCF method is made for 122 singlet excitation energies and 69 oscillator strengths. Of all CASSCF methods, when considering the complete test set, MC-RPA performs best for both excitation energies and oscillator strengths with a mean absolute error (MAE) of 0.74 eV and 51%, respectively. MC-TDA and SA-CASSCF show a similar accuracy for the excitation energies with a MAE of ∼1 eV with respect to more accurate coupled cluster (CC3) excitation energies. The opposite trend is observed for the subset of n → π* excitation energies for which SA-CASSCF exhibits the least deviations (MAE 0.65 eV). By looking at s-tetrazine in more detail, we conclude that better performance for the n → π* SA-CASSCF excitation energies can be attributed to a fortunate error compensation. For oscillator strengths, SA-CASSCF performs worst for the complete test set (MAE 100%) as well as for the subsets of n → π* (MAE 192%) and π → π* excitations (MAE 84.9%). In general, CASSCF gives the worst performance for excitation energies of all excited-state ab initio methods considered so far due to lacking the major part of dynamic electron correlation, though MC-RPA and TD-DFT (BP86) show similar performance. Among all LR-type methods, LR-CASSCF oscillator strengths are the ones with the least accuracy for the same reason. As state-specific orbital relaxation effects are accounted for in LR-CASSCF, oscillator strengths are significantly more accurate than those of MS-CASPT2. Our findings should encourage further developments of response theory-based multireference methods with higher accuracy and feasibility.
Collapse
Affiliation(s)
- Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung , Kaiser-Wilhelm-Platz 1 , D-45470 Mülheim an der Ruhr , Germany
| |
Collapse
|
24
|
Holzer C, Klopper W. Ionized, electron-attached, and excited states of molecular systems with spin–orbit coupling: Two-component GW and Bethe–Salpeter implementations. J Chem Phys 2019; 150:204116. [DOI: 10.1063/1.5094244] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P. O. Box 6980, D-76049 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P. O. Box 6980, D-76049 Karlsruhe, Germany
- Centre for Advanced Study (CAS) at The Norwegian Academy of Science and Letters, Drammensveien 78, N-0271 Oslo, Norway
| |
Collapse
|
25
|
Ziaei V, Bredow T. Screening mixing GW/Bethe-Salpeter approach for triplet states of organic molecules. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:395501. [PMID: 30124435 DOI: 10.1088/1361-648x/aadb75] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We recently proposed a screening mixing many-body ansatz (Ziaei and Bredow 2017 Phys. Rev. B 96 195115) to decrease the typical overestimation of the Hartree-Fock based GW/Bethe-Salpeter equation (BSE) singlet-singlet excitation energies in molecular systems. Now we have evaluated the accuracy of the proposed scheme for triplet states of a set of 20 organic molecules known as the Thiel set. We show that by mixing different screenings into GW and BSE calculated within random phase approximation (in order to ensure best gap and an optimal exciton binding energy), the total mean absolute error of 0.59 eV in the standard Hartree-Fock based eigenvalue GW/BSE approach is reduced to 0.26 eV for 63 triplet states. We further demonstrate that the quasi-particle self-consistent GW/BSE approach in which orbitals and energies are updated in the Green's function and the dynamically screened interaction mostly and considerably underestimates the excitation energies as shown for a few molecules.
Collapse
Affiliation(s)
- Vafa Ziaei
- Mulliken Center for Theoretical Chemistry, University of Bonn, 53115 Bonn, Germany
| | | |
Collapse
|
26
|
Ghosh S, Verma P, Cramer CJ, Gagliardi L, Truhlar DG. Combining Wave Function Methods with Density Functional Theory for Excited States. Chem Rev 2018; 118:7249-7292. [PMID: 30044618 DOI: 10.1021/acs.chemrev.8b00193] [Citation(s) in RCA: 135] [Impact Index Per Article: 22.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.
Collapse
Affiliation(s)
- Soumen Ghosh
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Christopher J Cramer
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Laura Gagliardi
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455-0431 , United States
| |
Collapse
|
27
|
Duchemin I, Guido CA, Jacquemin D, Blase X. The Bethe-Salpeter formalism with polarisable continuum embedding: reconciling linear-response and state-specific features. Chem Sci 2018; 9:4430-4443. [PMID: 29896384 PMCID: PMC5956976 DOI: 10.1039/c8sc00529j] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2018] [Accepted: 04/02/2018] [Indexed: 11/25/2022] Open
Abstract
The Bethe-Salpeter equation (BSE) formalism has been recently shown to be a valuable alternative to time-dependent density functional theory (TD-DFT) with the same computing time scaling with system size. In particular, problematic transitions for TD-DFT such as charge-transfer, Rydberg and cyanine-like excitations were shown to be accurately described with BSE. We demonstrate here that combining the BSE formalism with the polarisable continuum model (PCM) allows us to include simultaneously linear-response and state-specific contributions to solvatochromism. This is confirmed by exploring transitions of various natures (local, charge-transfer, etc.) in a series of solvated molecules (acrolein, indigo, p-nitro-aniline, donor-acceptor complexes, etc.) for which we compare BSE solvatochromic shifts to those obtained by linear-response and state-specific TD-DFT implementations. Such a remarkable and unique feature is particularly valuable for the study of solvent effects on excitations presenting a hybrid localised/charge-transfer character.
Collapse
Affiliation(s)
- Ivan Duchemin
- Univ. Grenobles Alpes , CEA, INAC-MEM, L_Sim , F-38000 Grenoble , France . ;
| | - Ciro A Guido
- Laboratoire CEISAM - UMR CNR 6230 , Université de Nantes , 2 Rue de la Houssinière, BP 92208 , 44322 Nantes Cedex 3 , France
- Laboratoire MOLTECH - UMR CNRS 6200 , Université de Angers , 2 Bd Lavoisier , 49045 Angers Cedex , France
| | - Denis Jacquemin
- Laboratoire CEISAM - UMR CNR 6230 , Université de Nantes , 2 Rue de la Houssinière, BP 92208 , 44322 Nantes Cedex 3 , France
| | - Xavier Blase
- Univ. Grenobles Alpes , CNRS , Institut Néel , F-38042 Grenoble , France
| |
Collapse
|
28
|
Wang X, Liu X, Cook C, Schatschneider B, Marom N. On the possibility of singlet fission in crystalline quaterrylene. J Chem Phys 2018; 148:184101. [DOI: 10.1063/1.5027553] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Xiaopeng Wang
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Xingyu Liu
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Cameron Cook
- Department of Chemistry and Biochemistry, California State Polytechnic University at Pomona, Pomona, California 91768, USA
| | - Bohdan Schatschneider
- Department of Chemistry and Biochemistry, California State Polytechnic University at Pomona, Pomona, California 91768, USA
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
29
|
Sharifzadeh S. Many-body perturbation theory for understanding optical excitations in organic molecules and solids. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:153002. [PMID: 29460855 DOI: 10.1088/1361-648x/aab0d1] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Semiconductors composed of organic molecules are promising as components for flexible and inexpensive optoelectronic devices, with many recent studies aimed at understanding their electronic and optical properties. In particular, computational modeling of these complex materials has provided new understanding of the underlying properties which give rise to their excited-state phenomena. This article provides an overview of recent many-body perturbation theory (MBPT) studies of optical excitations within organic molecules and solids. We discuss the accuracy of MBPT within the GW/BSE approach in predicting excitation energies and absorption spectra, and assess the impact of two commonly used approximations, the DFT starting point and the Tamm-Dancoff approximation. Moreover, we summarize studies that elucidate the role of solid-state structure on the nature of excitons in organic crystals. These studies show that a rich physical understanding of organic materials can be obtained from GW/BSE.
Collapse
Affiliation(s)
- Sahar Sharifzadeh
- Division of Materials Science and Engineering, Department of Electrical and Computer Engineering, Department of Physics, Boston University, Boston, MA, United States of America
| |
Collapse
|
30
|
Gui X, Holzer C, Klopper W. Accuracy Assessment of GW Starting Points for Calculating Molecular Excitation Energies Using the Bethe–Salpeter Formalism. J Chem Theory Comput 2018; 14:2127-2136. [DOI: 10.1021/acs.jctc.8b00014] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Xin Gui
- Institute of Physical Chemistry, Theoretical Chemistry Group, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Physical Chemistry, Theoretical Chemistry Group, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, D-76049 Karlsruhe, Germany
| | - Wim Klopper
- Institute of Physical Chemistry, Theoretical Chemistry Group, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, D-76049 Karlsruhe, Germany
- Centre for Advanced Study (CAS) at The Norwegian Academy of Science and Letters, Drammensveien 78, N-0271 Oslo, Norway
| |
Collapse
|
31
|
Roma G, Bruneval F, Martin-Samos L. Optical Properties of Saturated and Unsaturated Carbonyl Defects in Polyethylene. J Phys Chem B 2018; 122:2023-2030. [PMID: 29360369 DOI: 10.1021/acs.jpcb.7b12172] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Polyethylene (PE), one of the simplest and most used aliphatic polymers, is generally provided with a number of additives, in particular antioxidants, because of its tendency to get oxidized. Carbonyl defects, a product of the oxidation of PE, are occurring in various forms, in particular saturated ones, known as ketones, where a C═O double bond substitutes a CH2 group, and various unsaturated ones, i.e., with further missing hydrogens. Many experimental investigations of the optical properties in the visible/UV range mainly attribute the photoluminescence of PE to one specific kind of unsaturated carbonyls, following analogies to the emission spectra of similar small molecules. However, the reason why saturated carbonyls should not be optically detected is not clear. We investigated the optical properties of PE with and without carbonyl defects using perturbative GW and the Bethe-Salpeter equation in order to take into account excitonic effects. We discuss the calculated excitonic states in comparison with experimental absorption/emission energies and the stability of both saturated and unsaturated carbonyl defects. We conclude that the unsaturated defects are indeed the best candidate for the luminescence of oxidized PE, and the reason is mainly due to oscillator strengths.
Collapse
Affiliation(s)
- Guido Roma
- DEN-Service de Recherches de Métallurgie Physique, CEA, Université Paris-Saclay , F-91191 Gif sur Yvette, France
| | - Fabien Bruneval
- DEN-Service de Recherches de Métallurgie Physique, CEA, Université Paris-Saclay , F-91191 Gif sur Yvette, France
| | - Layla Martin-Samos
- Materials Research Laboratory, University of Nova Gorica , SI-5000 Nova Gorica, Slovenia.,CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, c/o SISSA Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265, 34136 Trieste Italy
| |
Collapse
|
32
|
Blase X, Duchemin I, Jacquemin D. The Bethe–Salpeter equation in chemistry: relations with TD-DFT, applications and challenges. Chem Soc Rev 2018; 47:1022-1043. [DOI: 10.1039/c7cs00049a] [Citation(s) in RCA: 120] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
We review the Bethe–Salpeter formalism and analyze its performances for the calculation of the excited state properties of molecular systems.
Collapse
Affiliation(s)
- Xavier Blase
- Univ. Grenoble Alpes
- CNRS
- Inst NEEL
- F-38042 Grenoble
- France
| | - Ivan Duchemin
- Univ. Grenoble Alpes
- CEA
- INAC-MEM
- L-Sim
- F-38000 Grenoble
| | - Denis Jacquemin
- CEISAM UMR CNRS 6230
- Université de Nantes
- 44322 Nantes Cedex 3
- France
| |
Collapse
|
33
|
Azarias C, Ponce-Vargas M, Navizet I, Fleurat-Lessard P, Romieu A, Le Guennic B, Richard JA, Jacquemin D. Rationalisation of the optical signatures of nor-dihydroxanthene-hemicyanine fused near-infrared fluorophores by first-principle tools. Phys Chem Chem Phys 2018; 20:12120-12128. [DOI: 10.1039/c8cp01587b] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The perfect pair towards more effective nor-DXH-hemicyanine fused NIR fluoroohores.
Collapse
Affiliation(s)
- Cloé Azarias
- Chimie Et Interdisciplinarité, Synthèse, Analyse, Modélisation (CEISAM)
- UMR CNRS no. 6230
- BP 92208
- Université de Nantes
- 44322 Nantes Cedex 3
| | - Miguel Ponce-Vargas
- Institut de Chimie Moléculaire de l'Université de Bourgogne (ICMUB, UMR-CNRS 6302)
- Université de Bourgogne Franche-Comté
- 21078 Dijon
- France
- Université Paris-Est
| | - Isabelle Navizet
- Université Paris-Est
- Laboratoire Modélisation et Simulation Multi Echelle
- MSME
- UMR 8208 CNRS
- UPEM
| | - Paul Fleurat-Lessard
- Institut de Chimie Moléculaire de l'Université de Bourgogne (ICMUB, UMR-CNRS 6302)
- Université de Bourgogne Franche-Comté
- 21078 Dijon
- France
| | - Anthony Romieu
- Institut de Chimie Moléculaire de l'Université de Bourgogne (ICMUB, UMR-CNRS 6302)
- Université de Bourgogne Franche-Comté
- 21078 Dijon
- France
- Institut Universitaire de France
| | - Boris Le Guennic
- Univ Rennes
- CNRS
- ISCR (Institut des Sciences Chimiques de Rennes) – UMR 6226
- 35000 Rennes
- France
| | - Jean-Alexandre Richard
- Organic Chemistry
- Institute of Chemical and Engineering Sciences
- (ICES)
- Agency for Science, Technology and Research (A_STAR)
- Singapore
| | - Denis Jacquemin
- Chimie Et Interdisciplinarité, Synthèse, Analyse, Modélisation (CEISAM)
- UMR CNRS no. 6230
- BP 92208
- Université de Nantes
- 44322 Nantes Cedex 3
| |
Collapse
|
34
|
Azarias C, Habert C, Budzák Š, Blase X, Duchemin I, Jacquemin D. Calculations of n→π* Transition Energies: Comparisons Between TD-DFT, ADC, CC, CASPT2, and BSE/GW Descriptions. J Phys Chem A 2017; 121:6122-6134. [DOI: 10.1021/acs.jpca.7b05222] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Cloé Azarias
- CEISAM,
BP 92208, UMR CNRS 6230, Université de Nantes, 2, Rue de
la Houssiniere, 44322 Nantes, Cedex 3, France
| | - Chloé Habert
- CEISAM,
BP 92208, UMR CNRS 6230, Université de Nantes, 2, Rue de
la Houssiniere, 44322 Nantes, Cedex 3, France
| | - Šimon Budzák
- Department
of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, SK-97400 Banská Bystrica, Slovak Republic
| | - Xavier Blase
- CNRS, Inst NEEL, F-38042 Grenoble, France
- Univ.
Grenoble Alpes, Inst NEEL, F-38042 Grenoble, France
| | - Ivan Duchemin
- Univ.
Grenoble Alpes, Inst NEEL, F-38042 Grenoble, France
- Univ.
Grenobles Alpes, CEA, INAC-MEM, L_Sim, F-38000 Grenoble, France
| | - Denis Jacquemin
- CEISAM,
BP 92208, UMR CNRS 6230, Université de Nantes, 2, Rue de
la Houssiniere, 44322 Nantes, Cedex 3, France
- Institut Universitaire de France, 1, rue Descartes, F-75231 Paris Cedex 05, France
| |
Collapse
|
35
|
Jacquemin D, Duchemin I, Blase X. Is the Bethe-Salpeter Formalism Accurate for Excitation Energies? Comparisons with TD-DFT, CASPT2, and EOM-CCSD. J Phys Chem Lett 2017; 8:1524-1529. [PMID: 28301726 PMCID: PMC5385505 DOI: 10.1021/acs.jpclett.7b00381] [Citation(s) in RCA: 64] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2017] [Accepted: 03/16/2017] [Indexed: 05/25/2023]
Abstract
Developing ab initio approaches able to provide accurate excited-state energies at a reasonable computational cost is one of the biggest challenges in theoretical chemistry. In that framework, the Bethe-Salpeter equation approach, combined with the GW exchange-correlation self-energy, which maintains the same scaling with system size as TD-DFT, has recently been the focus of a rapidly increasing number of applications in molecular chemistry. Using a recently proposed set encompassing excitation energies of many kinds [J. Phys. Chem. Lett. 2016, 7, 586-591], we investigate here the performances of BSE/GW. We compare these results to CASPT2, EOM-CCSD, and TD-DFT data and show that BSE/GW provides an accuracy comparable to the two wave function methods. It is particularly remarkable that the BSE/GW is equally efficient for valence, Rydberg, and charge-transfer excitations. In contrast, it provides a poor description of triplet excited states, for which EOM-CCSD and CASPT2 clearly outperform BSE/GW. This contribution therefore supports the use of the Bethe-Salpeter approach for spin-conserving transitions.
Collapse
Affiliation(s)
- Denis Jacquemin
- Laboratoire
CEISAM - UMR CNRS 6230, Université
de Nantes, 2 Rue de la
Houssinière, BP 92208, 44322 Nantes Cedex 3, France
- Institut
Universitaire de France, 1 rue Descartes, 75231 Paris Cedex 5, France
| | - Ivan Duchemin
- Univ.
Grenobles Alpes, CEA, INAC-MEM, L_Sim, F-38000 Grenoble, France
- Univ.
Grenoble Alpes, Inst NEEL, F-38042 Grenoble, France
| | - Xavier Blase
- Univ.
Grenoble Alpes, Inst NEEL, F-38042 Grenoble, France
- CNRS, Inst
NEEL, F-38042 Grenoble, France
| |
Collapse
|
36
|
Escudero D, Duchemin I, Blase X, Jacquemin D. Modeling the Photochrome-TiO 2 Interface with Bethe-Salpeter and Time-Dependent Density Functional Theory Methods. J Phys Chem Lett 2017; 8:936-940. [PMID: 28178780 DOI: 10.1021/acs.jpclett.7b00015] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Hybrid organic-inorganic semiconductor systems have important applications in both molecular electronics and photoresponsive materials. The characterizations of the interface and of the electronic excited-states of these hybrid systems remain a challenge for state-of-the-art computational methods, as the systems of interest are large. In the present investigation, we present for the first time a many-body Green's function Bethe-Salpeter investigation of a series of photochromic molecules adsorbed onto TiO2 nanoclusters. On the basis of these studies, the performance of time-dependent density functional theory (TD-DFT) calculations is assessed. In addition, the photochromic properties of different hybrid systems are also evaluated. This work shows that qualitatively different conclusions can be reached with TD-DFT relying on various exchange-correlation functionals for such organic-inorganic interfaces and paves the way to more accurate simulation of many hybrid materials.
Collapse
Affiliation(s)
- Daniel Escudero
- CEISAM UMR CNRS 6230, Université de Nantes , 2 rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Ivan Duchemin
- Univ. Grenoble Alpes, CEA, INAC-MEM, L-Sim , F-38000 Grenoble, France
- Univ. Grenoble Alpes, CNRS, Inst NEEL , F-38042 Grenoble, France
| | - Xavier Blase
- Univ. Grenoble Alpes, CNRS, Inst NEEL , F-38042 Grenoble, France
| | - Denis Jacquemin
- CEISAM UMR CNRS 6230, Université de Nantes , 2 rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
- Institut Universitaire de France , 1 rue Descartes, F-75005 Paris Cedex 05, France
| |
Collapse
|
37
|
Jacquemin D, Duchemin I, Blondel A, Blase X. Benchmark of Bethe-Salpeter for Triplet Excited-States. J Chem Theory Comput 2017; 13:767-783. [PMID: 28107000 DOI: 10.1021/acs.jctc.6b01169] [Citation(s) in RCA: 51] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
We have evaluated the accuracy of the Bethe-Salpeter singlet-triplet transition energies as well as singlet-triplet and triplet-triplet splittings for 20 organic molecules, using as reference the CC3 values determined by Thiel and co-workers with both the TZVP and aug-cc-pVTZ atomic basis sets. Our excitation energies are obtained on the basis of GW quasiparticle energy levels that are self-consistently converged with respect to the starting DFT eigenvalues. In its current form, BSE/GW is often unable to provide a balanced description of both singlet and triplet excited-states. While the singlet-singlet and triplet-triplet energy separations are obtained accurately, triplets are located too close in energy from the ground-state, by typically -0.55 eV when using standard functionals to generate the starting eigenstates. Applying the Tamm-Dancoff approximation upshifts the BSE triplet energies and allows reducing this error to ca. -0.40 eV, while using M06-HF eigenstates allows a further increase and hence a reduction of the error for triplet states, but at the cost of larger errors for the singlet excited-states. At this stage, the most accurate TD-DFT estimates therefore remain competitive for computing singlet-triplet transition energies. Indeed, with M06-2X, irrespective of the application or not of the Tamm-Dancoff approximation and of the selected atomic basis set, the deviations obtained with TD-DFT are rather small.
Collapse
Affiliation(s)
- Denis Jacquemin
- Laboratoire CEISAM - UMR CNR 6230, Université de Nantes , 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France.,Institut Universitaire de France , 103 bd St. Michel, 75005 Paris Cedex 5, France
| | - Ivan Duchemin
- INAC, SP2M/L_Sim, CEA/UJF Cedex 09, 38054 Grenoble, France.,Institut NEEL, Université Grenoble Alpes , F-38042 Grenoble, France
| | - Aymeric Blondel
- Laboratoire CEISAM - UMR CNR 6230, Université de Nantes , 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Xavier Blase
- Institut NEEL, Université Grenoble Alpes , F-38042 Grenoble, France.,CNRS, Institut NEEL, F-38042 Grenoble, France
| |
Collapse
|
38
|
Ziaei V, Bredow T. Large-Scale Quantum Many-Body Perturbation on Spin and Charge Separation in the Excited States of the Synthesized Donor-Acceptor Hybrid PBI-Macrocycle Complex. Chemphyschem 2017; 18:579-583. [DOI: 10.1002/cphc.201601244] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2016] [Revised: 12/23/2016] [Indexed: 11/10/2022]
Affiliation(s)
- Vafa Ziaei
- Mulliken Center for Theoretical Chemistry; University of Bonn; Germany
| | - Thomas Bredow
- Mulliken Center for Theoretical Chemistry; University of Bonn; Germany
| |
Collapse
|
39
|
Azarias C, Duchemin I, Blase X, Jacquemin D. Bethe-Salpeter study of cationic dyes: Comparisons with ADC(2) and TD-DFT. J Chem Phys 2017; 146:034301. [PMID: 28109224 DOI: 10.1063/1.4974097] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a theoretical investigation of the excited-state properties of a large series of structurally diverse arylcarbonium derivatives that are known to be challenging for theoretical models. More specifically, we compare the pros and cons of TD-DFT (TD-M06-2X), ADC(2), and BSE/GW approaches for a large panel of compounds, using two different solvent models. Both 0-0 and vertical transition energies are considered and compared to the experimental values. All approaches reasonably reproduce the auxochromic and acidochromic shifts, although in most cases both TD-DFT and BSE/GW return larger correlation with experimental values than ADC(2) for these shifts. In contrast, the absolute transition energies obtained with ADC(2) tend to be closer to the measurements, TD-DFT using the M06-2X functional largely overestimating the experimental references (by ca. 0.5 eV), and BSE/GW providing intermediate values. In addition, we show that the selected solvent model has a significant impact on the results, the corrected linear-response approach providing larger transition energies than its linear-response counterpart.
Collapse
Affiliation(s)
- Cloé Azarias
- CEISAM, UMR CNRS 6230, BP 92208, Université de Nantes, 2, Rue de la Houssinière, 44322 Nantes Cedex 3, France
| | - Ivan Duchemin
- INAC, SP2M/L_Sim, CEA/UJF Cedex 09, 38054 Grenoble, France
| | - Xavier Blase
- Univ. Grenoble Alpes, Inst NEEL, F-38042 Grenoble, France
| | - Denis Jacquemin
- CEISAM, UMR CNRS 6230, BP 92208, Université de Nantes, 2, Rue de la Houssinière, 44322 Nantes Cedex 3, France
| |
Collapse
|