1
|
Rodriguez-Mayorga M, Blase X, Duchemin I, D'Avino G. From Many-Body Ab Initio to Effective Excitonic Models: A Versatile Mapping Approach Including Environmental Embedding Effects. J Chem Theory Comput 2024; 20:8675-8688. [PMID: 39376072 DOI: 10.1021/acs.jctc.4c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/09/2024]
Abstract
We present an original multistate projective diabatization scheme based on Green's function formalisms that allows the systematic mapping of many-body ab initio calculations onto effective excitonic models. This method inherits the ability of the Bethe-Salpeter equation to describe Frenkel molecular excitons and intermolecular charge-transfer states equally well, as well as the possibility for an effective description of environmental effects in a QM/MM framework. The latter is found to be a crucial element in order to obtain accurate model parameters for condensed phases and to ensure their transferability to excitonic models for extended systems. The method is presented through a series of examples illustrating its quality, robustness, and internal consistency.
Collapse
Affiliation(s)
- Mauricio Rodriguez-Mayorga
- Grenoble Alpes University, CNRS, Grenoble INP, Institut Néel, 25 rue des Martyrs, Grenoble 38042, France
| | - Xavier Blase
- Grenoble Alpes University, CNRS, Grenoble INP, Institut Néel, 25 rue des Martyrs, Grenoble 38042, France
| | - Ivan Duchemin
- Grenoble Alpes University, CEA, IRIG-MEM-L Sim, Grenoble 38054, France
| | - Gabriele D'Avino
- Grenoble Alpes University, CNRS, Grenoble INP, Institut Néel, 25 rue des Martyrs, Grenoble 38042, France
| |
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
|
Wang X, Gao S, Luo Y, Liu X, Tom R, Zhao K, Chang V, Marom N. Computational Discovery of Intermolecular Singlet Fission Materials Using Many-Body Perturbation Theory. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:7841-7864. [PMID: 38774154 PMCID: PMC11103713 DOI: 10.1021/acs.jpcc.4c01340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/29/2024] [Revised: 04/16/2024] [Accepted: 04/17/2024] [Indexed: 05/24/2024]
Abstract
Intermolecular singlet fission (SF) is the conversion of a photogenerated singlet exciton into two triplet excitons residing on different molecules. SF has the potential to enhance the conversion efficiency of solar cells by harvesting two charge carriers from one high-energy photon, whose surplus energy would otherwise be lost to heat. The development of commercial SF-augmented modules is hindered by the limited selection of molecular crystals that exhibit intermolecular SF in the solid state. Computational exploration may accelerate the discovery of new SF materials. The GW approximation and Bethe-Salpeter equation (GW+BSE) within the framework of many-body perturbation theory is the current state-of-the-art method for calculating the excited-state properties of molecular crystals with periodic boundary conditions. In this Review, we discuss the usage of GW+BSE to assess candidate SF materials as well as its combination with low-cost physical or machine learned models in materials discovery workflows. We demonstrate three successful strategies for the discovery of new SF materials: (i) functionalization of known materials to tune their properties, (ii) finding potential polymorphs with improved crystal packing, and (iii) exploring new classes of materials. In addition, three new candidate SF materials are proposed here, which have not been published previously.
Collapse
Affiliation(s)
- Xiaopeng Wang
- School
of Foundational Education, University of
Health and Rehabilitation Sciences, Qingdao 266113, China
- Qingdao
Institute for Theoretical and Computational Sciences, Institute of
Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Siyu Gao
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Yiqun Luo
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Xingyu Liu
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Rithwik Tom
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Kaiji Zhao
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Vincent Chang
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Noa Marom
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
4
|
Knysh I, Raimbault D, Duchemin I, Blase X, Jacquemin D. Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels. J Chem Phys 2024; 160:144115. [PMID: 38602292 DOI: 10.1063/5.0203818] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 03/25/2024] [Indexed: 04/12/2024] Open
Abstract
We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Denez Raimbault
- 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
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
5
|
Cardia R, Dardenne N, Mula G, Pinna E, Rignanese GM, Charlier JC, Cappellini G. First-Principles Investigation of the Optical Properties of Eumelanin Protomolecules. J Phys Chem A 2023; 127:10797-10806. [PMID: 38109190 DOI: 10.1021/acs.jpca.3c04898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Using first-principles calculations, we investigate the absorption spectra (in the near-infrared, visible, and first UV range) of the two most probable eumelanin tetrameric molecules exhibiting either a linear open-chain or a cyclic porphyrine-like configuration. In order to simulate a realistic molecular system, an implicit solvent model is used in our calculations to mimic the effect of the solvated environment around the eumelanin molecule. Although the presence of solvent is found not to significantly affect the absorption pattern of both molecules, the onset of the spectra are shifted toward higher energies, especially for the linear tetramer. Interestingly, the absorption spectra and optical onsets of the two molecules differ significantly both in a vacuum and in ethanol. However, the two predicted spectra do not allow us to definitely discriminate between the two configurations when comparing the theoretical predictions with the available experimental spectrum. In addition, a mix of the two eumelanin configurations (close to fifty-fifty) leads to a maximum overlap between theoretical and experimental spectra. Consequently, this theoretical research shows that deeper insight can be gained using beyond DFT techniques on the real form of eumelanin protomolecules present in living systems as well as on their possible use in hybrid solar cells.
Collapse
Affiliation(s)
- Roberto Cardia
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
| | - Nicolas Dardenne
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain (UCLouvain), B-1348 Louvain-la-Neuve, Belgium
| | - Guido Mula
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
| | - Elisa Pinna
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
| | - Gian-Marco Rignanese
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain (UCLouvain), B-1348 Louvain-la-Neuve, Belgium
- European Theoretical Spectroscopy Facility (ETSF)
| | - Jean-Christophe Charlier
- Institute of Condensed Matter and Nanosciences, Université catholique de Louvain (UCLouvain), B-1348 Louvain-la-Neuve, Belgium
- European Theoretical Spectroscopy Facility (ETSF)
| | - Giancarlo Cappellini
- Department of Physics, Università degli Studi di Cagliari, Cittadella Universitaria I-09042 Monserrato, Cagliari, Italy
- European Theoretical Spectroscopy Facility (ETSF)
| |
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, 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
|
8
|
Rai V, Gerhard L, Balzer N, Valášek M, Holzer C, Yang L, Wegener M, Rockstuhl C, Mayor M, Wulfhekel W. Activating Electroluminescence of Charged Naphthalene Diimide Complexes Directly Adsorbed on a Metal Substrate. PHYSICAL REVIEW LETTERS 2023; 130:036201. [PMID: 36763403 DOI: 10.1103/physrevlett.130.036201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 09/02/2022] [Accepted: 12/13/2022] [Indexed: 06/18/2023]
Abstract
Electroluminescence from single molecules adsorbed on a conducting surface imposes conflicting demands for the molecule-electrode coupling. To conduct electrons, the molecular orbitals need to be hybridized with the electrodes. To emit light, they need to be decoupled from the electrodes to prevent fluorescence quenching. Here, we show that fully quenched 2,6-core-substituted naphthalene diimide derivative in a self-assembled monolayer directly deposited on a Au(111) surface can be activated with the tip of a scanning tunneling microscope to decouple the relevant frontier orbitals from the metallic substrate. In this way, individual molecules can be driven from a strongly hybridized state with quenched luminescence to a light-emitting state. The emission performance compares in terms of quantum efficiency, stability, and reproducibility to that of single molecules deposited on thin insulating layers. Quantum chemical calculations suggest that the emitted light originates from the singly charged cationic pair of the molecules.
Collapse
Affiliation(s)
- Vibhuti Rai
- Institute for Quantum Materials and Technologies, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
| | - Lukas Gerhard
- Institute for Quantum Materials and Technologies, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
| | - Nico Balzer
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
| | - Michal Valášek
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), D-76128 Karlsruhe, Germany
| | - Liang Yang
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
- Institute of Applied Physics, Karlsruhe Institute of Technology (KIT), D-76128 Karlsruhe, Germany
| | - Martin Wegener
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
- Institute of Applied Physics, Karlsruhe Institute of Technology (KIT), D-76128 Karlsruhe, Germany
| | - Carsten Rockstuhl
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), D-76128 Karlsruhe, Germany
| | - Marcel Mayor
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
- Department of Chemistry, University of Basel, St. Johannsring 19, CH-4056 Basel, Switzerland
- Lehn Institute of Functional Materials (LIFM), Sun Yat-Sen University (SYSU), Xingang West Road, Guangzhou, China
| | - Wulf Wulfhekel
- Institute for Quantum Materials and Technologies, Karlsruhe Institute of Technology (KIT), D-76344 Eggenstein-Leopoldshafen, Germany
- Physikalisches Institut, Karlsruhe Institute of Technology (KIT), D-76128 Karlsruhe, Germany
| |
Collapse
|
9
|
Duan X, Zhang J, Liu Y, Zhang M, Jiang YN, Ma Y. Crucial Role of Defect States in the Ultralong Phosphorescence of Organic Molecular Crystals. J Phys Chem Lett 2023; 14:230-236. [PMID: 36594617 DOI: 10.1021/acs.jpclett.2c03604] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Ultralong organic phosphorescence (UOP) in pure organic molecular crystals has attracted a lot of interest recently. There is much debate on the emission mechanism of this UOP. Two recent experimental works published in Nat. Photonics 2019, 13, 406-411 and Nat. Mater. 2021, 20, 175-180 attribute UOP in the 2,4,6-trimethoxy-1,3,5-triazine (TMOT) crystals and the carbazole crystals to H-aggregation of the TMOT molecules or the formation of charge-transfer excitons between the carbazole and impurity molecules. Our first-principles many-body Green's function theory calculations show that the lowest triplet states of these two crystals are in fact the localized defect states originating from the twisted TMOT molecules and the impurities, respectively. Energies of the H-aggregation-induced exciton and the charge-transfer exciton are too high to account for UOP. UOP should be mainly due to the little orbital overlap between the localized defect state and the delocalized band edges of the crystal. Strong intermolecular interactions suppress nonradiative decay of the triplet exciton localized on the defect.
Collapse
Affiliation(s)
- Xiujuan Duan
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Jie Zhang
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Yaru Liu
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Min Zhang
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Ya-Nan Jiang
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Yuchen Ma
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| |
Collapse
|
10
|
Li J, Golze D, Yang W. Combining Renormalized Singles GW Methods with the Bethe-Salpeter Equation for Accurate Neutral Excitation Energies. J Chem Theory Comput 2022; 18:6637-6645. [PMID: 36279250 PMCID: PMC9972216 DOI: 10.1021/acs.jctc.2c00686] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We apply the renormalized singles (RS) Green's function in the Bethe-Salpeter equation (BSE)/GW approach to predict accurate neutral excitation energies of molecular systems. The BSE calculations are performed on top of the GRSWRS method, which uses the RS Green's function also for the computation of the screened Coulomb interaction W. We show that the BSE/GRSWRS approach significantly outperforms BSE/G0W0 for predicting excitation energies of valence, Rydberg, and charge-transfer (CT) excitations by benchmarking the Truhlar-Gagliardi set, Stein CT set, and an atomic Rydberg test set. For the Truhlar-Gagliardi test set, BSE/GRSWRS provides comparable accuracy to time-dependent density functional theory (TDDFT) and is slightly better than BSE starting from eigenvalue self-consistent GW (evGW). For the Stein CT test set, BSE/GRSWRS significantly outperforms BSE/G0W0 and TDDFT with the accuracy comparable to BSE/evGW. We also show that BSE/GRSWRS predicts Rydberg excitation energies of atomic systems well. Besides the excellent accuracy, BSE/GRSWRS largely eliminates the dependence on the choice of the density functional approximation. This work demonstrates that the BSE/GRSWRS approach is accurate and efficient for predicting excitation energies for a broad range of systems, which expands the applicability of the BSE/GW approach.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Dorothea Golze
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01062 Dresden, Germany
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
11
|
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
|
12
|
Förster A, Visscher L. Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems. J Chem Theory Comput 2022; 18:6779-6793. [PMID: 36201788 PMCID: PMC9648197 DOI: 10.1021/acs.jctc.2c00531] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The GW-Bethe–Salpeter equation
(BSE) method
is promising for calculating the low-lying excitonic states of molecular
systems. However, so far it has only been applied to rather small
molecules and in the commonly implemented diagonal approximations
to the electronic self-energy, it depends on a mean-field starting
point. We describe here an implementation of the self-consistent and
starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on
large molecules. We herein show that eigenvalue-only self-consistency
can lead to an unfaithful description of some excitonic states for
chlorophyll dimers while the qsGW-BSE vertical excitation
energies (VEEs) are in excellent agreement with spectroscopic experiments
for chlorophyll monomers and dimers measured in the gas phase. Furthermore,
VEEs from time-dependent density functional theory calculations tend
to disagree with experimental values and using different range-separated
hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use
the new qsGW-BSE implementation to calculate the
lowest excitation energies of the six chromophores of the photosystem
II (PSII) reaction center (RC) with nearly 2000 correlated electrons.
Using more than 11,000 (6000) basis functions, the calculation could
be completed in less than 5 (2) days on a single modern compute node.
In agreement with previous TD-DFT calculations using RSH kernels on
models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters
in the low-energy spectrum of the hexameric complex. Earlier works
with RSH kernels have demonstrated that the protein environment facilitates
the experimentally observed interchromophoric charge transfer. Therefore,
future research will need to combine correlation effects beyond TD-DFT
with an explicit treatment of environmental electrostatics.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HVAmsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HVAmsterdam, The Netherlands
| |
Collapse
|
13
|
Li J, Jin Y, Su NQ, Yang W. Combining localized orbital scaling correction and Bethe-Salpeter equation for accurate excitation energies. J Chem Phys 2022; 156:154101. [PMID: 35459294 PMCID: PMC9033305 DOI: 10.1063/5.0087498] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2022] [Accepted: 03/29/2022] [Indexed: 11/14/2022] Open
Abstract
We applied localized orbital scaling correction (LOSC) in Bethe-Salpeter equation (BSE) to predict accurate excitation energies for molecules. LOSC systematically eliminates the delocalization error in the density functional approximation and is capable of approximating quasiparticle (QP) energies with accuracy similar to or better than GW Green's function approach and with much less computational cost. The QP energies from LOSC, instead of commonly used G0W0 and evGW, are directly used in BSE. We show that the BSE/LOSC approach greatly outperforms the commonly used BSE/G0W0 approach for predicting excitations with different characters. For the calculations of Truhlar-Gagliardi test set containing valence, charge transfer, and Rydberg excitations, BSE/LOSC with the Tamm-Dancoff approximation provides a comparable accuracy to time-dependent density functional theory (TDDFT) and BSE/evGW. For the calculations of Stein CT test set and Rydberg excitations of atoms, BSE/LOSC considerably outperforms both BSE/G0W0 and TDDFT approaches with a reduced starting point dependence. BSE/LOSC is, thus, a promising and efficient approach to calculate excitation energies for molecular systems.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Ye Jin
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Neil Qiang Su
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
14
|
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
|
15
|
Pham NNT, Han SH, Park JS, Lee SG. Optical and Electronic Properties of Organic NIR-II Fluorophores by Time-Dependent Density Functional Theory and Many-Body Perturbation Theory: GW-BSE Approaches. NANOMATERIALS (BASEL, SWITZERLAND) 2021; 11:2293. [PMID: 34578610 PMCID: PMC8466807 DOI: 10.3390/nano11092293] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Revised: 08/30/2021] [Accepted: 09/01/2021] [Indexed: 11/16/2022]
Abstract
Organic-molecule fluorophores with emission wavelengths in the second near-infrared window (NIR-II, 1000-1700 nm) have attracted substantial attention in the life sciences and in biomedical applications because of their excellent resolution and sensitivity. However, adequate theoretical levels to provide efficient and accurate estimations of the optical and electronic properties of organic NIR-II fluorophores are lacking. The standard approach for these calculations has been time-dependent density functional theory (TDDFT). However, the size and large excitonic energies of these compounds pose challenges with respect to computational cost and time. In this study, we used the GW approximation combined with the Bethe-Salpeter equation (GW-BSE) implemented in many-body perturbation theory approaches based on density functional theory. This method was used to perform calculations of the excited states of two NIR molecular fluorophores (BTC980 and BTC1070), going beyond TDDFT. In this study, the optical absorption spectra and frontier molecular orbitals of these compounds were compared using TDDFT and GW-BSE calculations. The GW-BSE estimates showed excellent agreement with previously reported experimental results.
Collapse
Affiliation(s)
- Nguyet N. T. Pham
- School of Chemical Engineering, Pusan National University, Busan 46241, Korea; (N.N.T.P.); (S.H.H.)
| | - Seong Hun Han
- School of Chemical Engineering, Pusan National University, Busan 46241, Korea; (N.N.T.P.); (S.H.H.)
| | - Jong S. Park
- School of Chemical Engineering, Pusan National University, Busan 46241, Korea; (N.N.T.P.); (S.H.H.)
- Department of Organic Material Science and Engineering, Pusan National University, Busan 46241, Korea
| | - Seung Geol Lee
- School of Chemical Engineering, Pusan National University, Busan 46241, Korea; (N.N.T.P.); (S.H.H.)
- Department of Organic Material Science and Engineering, Pusan National University, Busan 46241, Korea
| |
Collapse
|
16
|
Förster A, Visscher L. Low-Order Scaling Quasiparticle Self-Consistent GW for Molecules. Front Chem 2021; 9:736591. [PMID: 34540804 PMCID: PMC8446457 DOI: 10.3389/fchem.2021.736591] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2021] [Accepted: 08/13/2021] [Indexed: 11/13/2022] Open
Abstract
Low-order scaling GW implementations for molecules are usually restricted to approximations with diagonal self-energy. Here, we present an all-electron implementation of quasiparticle self-consistent GW for molecular systems. We use an efficient algorithm for the evaluation of the self-energy in imaginary time, from which a static non-local exchange-correlation potential is calculated via analytical continuation. By using a direct inversion of iterative subspace method, fast and stable convergence is achieved for almost all molecules in the GW100 database. Exceptions are systems which are associated with a breakdown of the single quasiparticle picture in the valence region. The implementation is proven to be starting point independent and good agreement of QP energies with other codes is observed. We demonstrate the computational efficiency of the new implementation by calculating the quasiparticle spectrum of a DNA oligomer with 1,220 electrons using a basis of 6,300 atomic orbitals in less than 4 days on a single compute node with 16 cores. We use then our implementation to study the dependence of quasiparticle energies of DNA oligomers consisting of adenine-thymine pairs on the oligomer size. The first ionization potential in vacuum decreases by nearly 1 electron volt and the electron affinity increases by 0.4 eV going from the smallest to the largest considered oligomer. This shows that the DNA environment stabilizes the hole/electron resulting from photoexcitation/photoattachment. Upon inclusion of the aqueous environment via a polarizable continuum model, the differences between the ionization potentials reduce to 130 meV, demonstrating that the solvent effectively compensates for the stabilizing effect of the DNA environment. The electron affinities of the different oligomers are almost identical in the aqueous environment.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, Amsterdam, Netherlands
| | | |
Collapse
|
17
|
Assessment of DFT methods for the prediction of detachment energies and electronic structures of complex and multiply charged anions. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113295] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
|
18
|
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
|
19
|
Duchemin I, Blase X. Cubic-Scaling All-Electron GW Calculations with a Separable Density-Fitting Space-Time Approach. J Chem Theory Comput 2021; 17:2383-2393. [PMID: 33797245 DOI: 10.1021/acs.jctc.1c00101] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present an implementation of the GW space-time approach that allows cubic-scaling all-electron calculations with standard Gaussian basis sets without exploiting any localization or sparsity considerations. The independent-electron susceptibility is constructed in a time representation over a nonuniform distribution of real-space locations {rk} optimized within a separable resolution-of-the-identity framework to reproduce standard Coulomb-fitting calculations with meV accuracy. The compactness of the obtained {rk} distribution leads to a crossover with the standard Coulomb-fitting scheme for system sizes below a few hundred electrons. The needed analytic continuation follows a recent approach that requires the continuation of the screened Coulomb potential rather than the much more structured self-energy. The present scheme is benchmarked over large molecular sets, and scaling properties are demonstrated on a family of defected hexagonal boron-nitride flakes containing up to 6000 electrons.
Collapse
Affiliation(s)
- Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L_Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Inst NEEL, F-38042 Grenoble, France
| |
Collapse
|
20
|
Hashemi Z, Leppert L. Assessment of the Ab Initio Bethe-Salpeter Equation Approach for the Low-Lying Excitation Energies of Bacteriochlorophylls and Chlorophylls. J Phys Chem A 2021; 125:2163-2172. [PMID: 33656894 PMCID: PMC8028335 DOI: 10.1021/acs.jpca.1c01240] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Bacteriochlorophyll and chlorophyll molecules are crucial building blocks of the photosynthetic apparatus in bacteria, algae, and plants. Embedded in transmembrane protein complexes, they are responsible for the primary processes of photosynthesis: excitation energy and charge transfer. Here, we use ab initio many-body perturbation theory within the GW approximation and Bethe-Salpeter equation (BSE) approach to calculate the electronic structure and optical excitations of bacteriochlorophylls a, b, c, d, and e and chlorophylls a and b. We systematically study the effects of the structure, basis set size, partial self-consistency in GW, and the underlying exchange-correlation approximation and compare our calculations with results from time-dependent density functional theory, multireference RASPT2, and experimental literature results. We find that optical excitations calculated with GW+BSE are in excellent agreement with experimental data, with an average deviation of less than 100 meV for the first three bright excitations of the entire family of (bacterio)chlorophylls. Contrary to state-of-the-art time-dependent density functional theory (TDDFT) with an optimally tuned range-separated hybrid functional, this accuracy is achieved in a parameter-free approach. Moreover, GW+BSE predicts the energy differences between the low-energy excitations correctly and eliminates spurious charge transfer states that TDDFT with (semi)local approximations is known to produce. Our study provides accurate reference results and highlights the potential of the GW+BSE approach for the simulation of larger pigment complexes.
Collapse
Affiliation(s)
- Zohreh Hashemi
- Institute of Physics, University of Bayreuth, Bayreuth 95440, Germany
| | - Linn Leppert
- Institute of Physics, University of Bayreuth, Bayreuth 95440, Germany.,MESA+ Institute for Nanotechnology, University of Twente, 7500 AE Enschede, The Netherlands
| |
Collapse
|
21
|
Osella S. Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead? NANOMATERIALS 2021; 11:nano11020299. [PMID: 33498961 PMCID: PMC7911014 DOI: 10.3390/nano11020299] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 01/19/2021] [Accepted: 01/20/2021] [Indexed: 12/13/2022]
Abstract
A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.
Collapse
Affiliation(s)
- Silvio Osella
- Chemical and Biological Systems Simulation Lab, Center of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland
| |
Collapse
|
22
|
Cortés-Mejía R, Höfener S, Klopper W. Effects of rotational conformation on electronic properties of 4,4′-bis(carbazol-9-yl)biphenyl (CBP): the single-molecule picture and beyond. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1876936] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Rodrigo Cortés-Mejía
- Theoretical Chemistry Group, Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe, Germany
| | - Sebastian Höfener
- 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
|
23
|
Tirimbò G, Baumeier B. Ab initio modeling of excitons: from perfect crystals to biomaterials. ADVANCES IN PHYSICS: X 2021. [DOI: 10.1080/23746149.2021.1912638] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022] Open
Affiliation(s)
- Gianluca Tirimbò
- Department of Mathematics and Computer Science, Eindhoven University of Technology, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Björn Baumeier
- Department of Mathematics and Computer Science, Eindhoven University of Technology, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, Eindhoven, The Netherlands
| |
Collapse
|
24
|
Adhikari S, Santra B, Ruan S, Bhattarai P, Nepal NK, Jackson KA, Ruzsinszky A. The Fermi–Löwdin self-interaction correction for ionization energies of organic molecules. J Chem Phys 2020; 153:184303. [DOI: 10.1063/5.0024776] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Santosh Adhikari
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Shiqi Ruan
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Puskar Bhattarai
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Niraj K. Nepal
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Koblar A. Jackson
- Department of Physics and Science of Advanced Materials Program, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Adrienn Ruzsinszky
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| |
Collapse
|
25
|
Blase X, Duchemin I, Jacquemin D, Loos PF. The Bethe-Salpeter Equation Formalism: From Physics to Chemistry. J Phys Chem Lett 2020; 11:7371-7382. [PMID: 32787315 DOI: 10.1021/acs.jpclett.0c01875] [Citation(s) in RCA: 71] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called GW approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/vis spectra, has shown to provide accurate singlet excitation energies with a typical error of 0.1-0.3 eV. With a similar computational cost as time-dependent density-functional theory (TD-DFT), BSE is able to provide an accuracy on par with the most accurate global and range-separated hybrid functionals without the unsettling choice of the exchange-correlation functional, resolving further known issues (e.g., charge-transfer excitations). In this Perspective, we provide a historical overview of BSE, with a particular focus on its condensed-matter roots. We also propose a critical review of its strengths and weaknesses in different chemical situations.
Collapse
Affiliation(s)
- Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Denis Jacquemin
- 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, Toulouse, France
| |
Collapse
|
26
|
Shee J, Head-Gordon M. Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile. J Chem Theory Comput 2020; 16:6244-6255. [DOI: 10.1021/acs.jctc.0c00635] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
27
|
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
|
28
|
Loos PF, Scemama A, Duchemin I, Jacquemin D, Blase X. Pros and Cons of the Bethe-Salpeter Formalism for Ground-State Energies. J Phys Chem Lett 2020; 11:3536-3545. [PMID: 32298578 DOI: 10.1021/acs.jpclett.0c00460] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The combination of the many-body Green's function GW approximation and the Bethe-Salpeter equation (BSE) formalism has shown to be a promising alternative to time-dependent density functional theory (TD-DFT) for computing vertical transition energies and oscillator strengths in molecular systems. The BSE formalism can also be employed to compute ground-state correlation energies thanks to the adiabatic-connection fluctuation-dissipation theorem (ACFDT). Here, we study the topology of the ground-state potential energy surfaces (PESs) of several diatomic molecules near their equilibrium bond length. Using comparisons with state-of-art computational approaches (CC3), we show that ACFDT@BSE is surprisingly accurate and can even compete with lower-order coupled cluster methods (CC2 and CCSD) in terms of total energies and equilibrium bond distances for the considered systems. However, we sometimes observe unphysical irregularities on the ground-state PES in relation with difficulties in the identification of a few GW quasiparticle energies.
Collapse
Affiliation(s)
- Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31077 Toulouse, France
| | - Anthony Scemama
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31077 Toulouse, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L Sim, 38054 Grenoble, France
| | - Denis Jacquemin
- Laboratoire CEISAM - UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut NEEL, F-38042 Grenoble, France
| |
Collapse
|
29
|
Valencia AM, Guerrini M, Cocchi C. Ab initio modelling of local interfaces in doped organic semiconductors. Phys Chem Chem Phys 2020; 22:3527-3538. [PMID: 31994551 DOI: 10.1039/c9cp06655a] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Doping in organic semiconductors remains a debated issue from both an experimental and ab initio perspective. Due to the complexity of these systems, which exhibit a low degree of crystallinity and high level of disorder, modelling doped organic semiconductors from first-principles calculations is not a trivial task, as their electronic and optical properties are sensitive to the choice of initial geometries. A crucial aspect to take into account, in view of rationalizing the electronic structure of these materials through ab initio calculations, is the role of local donor/acceptor interfaces. We address this problem in the framework of state-of-the-art density-functional theory and many-body perturbation theory, investigating the structural, electronic, and optical properties of quaterthiophene and sexithiophene oligomers doped by 2,3,5,6-tetrafluoro-7,7,8,8-tetracyano-quinodimethane (F4TCNQ). We consider different model structures ranging from isolated dimers and trimers, to periodic stacks. Our results demonstrate that the choice of the initial geometry critically impacts the resulting electronic structure and the degree of charge transfer in the materials, depending on the amount and on the nature of the local interfaces between donor and acceptor species. The optical spectra appear less sensitive to these parameters at least from a first glance, although a quantitative analysis of the excitations reveals that their Frenkel or charge-transfer character is affected by the characteristics of the donor/acceptor interfaces as well as by the donor length. Our findings represent an important step forward towards an insightful first-principles description of the microscopic properties of doped organic semiconductors complementary to experiments.
Collapse
Affiliation(s)
- Ana M Valencia
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, Zum Großen Windkanal 6, 12489 Berlin, Germany.
| | - Michele Guerrini
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, Zum Großen Windkanal 6, 12489 Berlin, Germany.
| | - Caterina Cocchi
- Humboldt-Universität zu Berlin, Physics Department and IRIS Adlershof, Zum Großen Windkanal 6, 12489 Berlin, Germany.
| |
Collapse
|
30
|
Fujita T, Noguchi Y, Hoshi T. Charge-transfer excited states in the donor/acceptor interface from large-scale GW calculations. J Chem Phys 2019; 151:114109. [PMID: 31542033 DOI: 10.1063/1.5113944] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Predicting the charge-transfer (CT) excited states across the donor/acceptor (D/A) interface is essential for understanding the charge photogeneration process in an organic solar cell. Here, we present a fragment-based GW implementation that can be applied to a D/A interface structure and thus enables accurate determination of the CT states. The implementation is based on the fragmentation approximation of the polarization function and the combined GW and Coulomb-hole plus screened exchange approximations for self-energies. The fragment-based GW is demonstrated by application to the pentacene/C60 interface structure containing more than 2000 atoms. The CT excitation energies were estimated from the quasiparticle energies and electron-hole screened Coulomb interactions; the computed energies are in reasonable agreement with experimental estimates from the external quantum efficiency measurements. We highlight the impact of the induced polarization effects on the electron-hole energetics. The proposed fragment-based GW method offers a first-principles tool to compute the quasiparticle energies and electronic excitation energies of organic materials.
Collapse
Affiliation(s)
| | - Yoshifumi Noguchi
- Department of Applied Chemistry and Biochemical Engineering, Graduate School of Engineering, Shizuoka University, Hamamatsu, Shizuoka 432-8561, Japan
| | - Takeo Hoshi
- Department of Applied Mathematics and Physics, Tottori University, Tottori 680-8550, Japan
| |
Collapse
|
31
|
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
|
32
|
Leng X, Jin F, Wei M, Ma H, Feng J, Ma Y. Electronic energy transfer studied by many-body Green’s function theory. J Chem Phys 2019; 150:164107. [DOI: 10.1063/1.5066290] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Xia Leng
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Fan Jin
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Min Wei
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Huizhong Ma
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Jin Feng
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| | - Yuchen Ma
- School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, China
| |
Collapse
|
33
|
Jin Y, Yang W. Excitation Energies from the Single-Particle Green's Function with the GW Approximation. J Phys Chem A 2019; 123:3199-3204. [PMID: 30920830 DOI: 10.1021/acs.jpca.9b02379] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quasi-particle energies are important in predicting molecular ionization energies and bulk band structures. The state-of-the-art method for quasi-particle energy calculations, particularly for bulk systems, is the GW approximation. For excited state calculations, one needs to go beyond the GW approximation. The Bethe-Salpeter equation (BSE) is the commonly used approach for bulk-system excited state calculations beyond the GW approximation, which is accurate but computationally cumbersome. In this Article, we develop a new method to extract excitation energies directly from the quasi-particle energies based on the GW approximation. Starting from the ( N - 1)-electron system, we are able to calculate molecular excitation energies with orbital energies at the GW level for HOMO excitations. Our calculations demonstrate that this method can accurately capture low-lying local excitations as well as charge transfer excitations in many molecular systems. Our method is shown to outperform the time-dependent density functional theory (TDDFT) and are comparable with higher level excited state calculations, including the equation-of-motion couple cluster (EOM-CC) theory and the BSE, but with less computational effort. This new approach provides an efficient alternative to the BSE method for accurate excited state calculations.
Collapse
Affiliation(s)
- Ye Jin
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States
| | - Weitao Yang
- Department of Chemistry , Duke University , Durham , North Carolina 27708 , United States.,Key Laboratory of Theoretical Chemistry of Environment, Ministry of Education, School of Chemistry and Environment , South China Normal University , Guangzhou 510006 , China
| |
Collapse
|
34
|
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
|
35
|
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
|
36
|
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
|
37
|
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
|
38
|
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
|
39
|
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
|
40
|
Holzer C, Klopper W. Communication: Symmetry-adapted perturbation theory with intermolecular induction and dispersion energies from the Bethe–Salpeter equation. J Chem Phys 2017; 147:181101. [DOI: 10.1063/1.5007929] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- 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
|
41
|
Coccia E, Varsano D, Guidoni L. Theoretical S1 ← S0 Absorption Energies of the Anionic Forms of Oxyluciferin by Variational Monte Carlo and Many-Body Green’s Function Theory. J Chem Theory Comput 2017; 13:4357-4367. [DOI: 10.1021/acs.jctc.7b00505] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Emanuele Coccia
- S3
Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | - Daniele Varsano
- S3
Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | - Leonardo Guidoni
- Dipartimento
di Scienze Fisiche e Chimiche, Universitá degli Studi dell’Aquila, via Vetoio, 67100, L’Aquila, Italy
| |
Collapse
|
42
|
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
|
43
|
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
|
44
|
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
|
45
|
Duchemin I, Jacquemin D, Blase X. Combining the GW formalism with the polarizable continuum model: A state-specific non-equilibrium approach. J Chem Phys 2017; 144:164106. [PMID: 27131530 DOI: 10.1063/1.4946778] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
We have implemented the polarizable continuum model within the framework of the many-body Green's function GW formalism for the calculation of electron addition and removal energies in solution. The present formalism includes both ground-state and non-equilibrium polarization effects. In addition, the polarization energies are state-specific, allowing to obtain the bath-induced renormalisation energy of all occupied and virtual energy levels. Our implementation is validated by comparisons with ΔSCF calculations performed at both the density functional theory and coupled-cluster single and double levels for solvated nucleobases. The present study opens the way to GW and Bethe-Salpeter calculations in disordered condensed phases of interest in organic optoelectronics, wet chemistry, and biology.
Collapse
Affiliation(s)
- Ivan Duchemin
- INAC, SP2M/L_Sim, CEA/UJF Cedex 09, 38054 Grenoble, France
| | - Denis Jacquemin
- Laboratoire CEISAM - UMR CNR 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208, 44322 Nantes Cedex 3, France
| | | |
Collapse
|
46
|
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
|
47
|
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
|
48
|
Varsano D, Caprasecca S, Coccia E. Theoretical description of protein field effects on electronic excitations of biological chromophores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:013002. [PMID: 27830666 DOI: 10.1088/0953-8984/29/1/013002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.
Collapse
Affiliation(s)
- Daniele Varsano
- S3 Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | | | | |
Collapse
|
49
|
Feng J, Liu G, Yuan S, Ma Y. Influence of functional groups on water splitting in carbon nanodot and graphitic carbon nitride composites: a theoretical mechanism study. Phys Chem Chem Phys 2017; 19:4997-5003. [DOI: 10.1039/c6cp08622e] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Charge transfer and separation mechanism in C–CHO/g-C3N4 under light irradiation through a two-step process (Route I) and direct excitation (Route II).
Collapse
Affiliation(s)
- Jin Feng
- School of Chemistry and Chemical Engineering
- Shandong University
- China
| | - Guokui Liu
- School of Chemistry and Chemical Engineering
- Shandong University
- China
| | - Shiling Yuan
- School of Chemistry and Chemical Engineering
- Shandong University
- China
| | - Yuchen Ma
- School of Chemistry and Chemical Engineering
- Shandong University
- China
| |
Collapse
|
50
|
Senn F, Park YC. Constricted variational density functional theory for spatially clearly separated charge-transfer excitations. J Chem Phys 2016; 145:244108. [PMID: 28049328 DOI: 10.1063/1.4972231] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Constricted Variational Density Functional Theory (CV-DFT) is known to be one of the successful methods in predicting charge-transfer excitation energies. In this paper, we apply the CV-DFT method to the well-known model systems ethylene-tetrafluoroethylene (C2H4 × C2F4) and the zincbacteriochlorin-bacteriochlorin complex (ZnBC-BC). The analysis of the CV-DFT energies enables us to understand the -1/R charge-transfer behaviour in CV-DFT for large separation distances R. With this we discuss the importance of orbital relaxations using the relaxed version of CV(∞)-DFT, the R-CV(∞)-DFT method. Possible effects of the optimization of the transition matrix for the relaxed self-consistent field version of CV(∞)-DFT, RSCF-CV(∞)-DFT in the case of large fragment separations are shown and we introduce two possible gradient restrictions to avoid the unwanted admixing of other transitions.
Collapse
Affiliation(s)
- Florian Senn
- Department of Chemistry, University of Calgary, 2500 University Drive NorthWest, Calgary, Alberta T2N 1N4, Canada
| | - Young Choon Park
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon 305-701, South Korea
| |
Collapse
|