1
|
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
|
2
|
Abraham V, Harsha G, Zgid D. Relativistic Fully Self-Consistent GW for Molecules: Total Energies and Ionization Potentials. J Chem Theory Comput 2024; 20:4579-4590. [PMID: 38778459 DOI: 10.1021/acs.jctc.4c00075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The fully self-consistent GW (scGW) method with an iterative solution of the Dyson equation provides a consistent approach for describing the ground and excited states without any dependence on the mean-field reference. In this work, we present a relativistic version of scGW for molecules containing heavy elements using the exact two-component (X2C) Coulomb approximation. We benchmark SOC-81 data set containing closed shell heavy elements for the first ionization potential using the fully self-consistent GW as well as one-shot GW. The self-consistent GW provides superior results compared to G0W0 with PBE reference and comparable results to G0W0 with PBE0 while also removing the starting point dependence. The photoelectron spectra obtained at the X2C level demonstrate very good agreement with the experimental spectra. We also observe that scGW provides very good estimation of ionization potential for the inner d-shell orbitals. Additionally, using the well-conserved total energy, we investigate the equilibrium bond length and harmonic frequencies of a few halogen dimers using scGW. Overall, our findings demonstrate the applicability of the fully self-consistent GW method for accurate ionization potential, photoelectron spectra, and total energies in finite systems with heavy elements with a reasonable computational scaling.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Gaurav Harsha
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics and Astronomy, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
3
|
Jacquemin D, Kossoski F, Gam F, Boggio-Pasqua M, Loos PF. Reference Vertical Excitation Energies for Transition Metal Compounds. J Chem Theory Comput 2023. [PMID: 37965941 DOI: 10.1021/acs.jctc.3c01080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
To enrich and enhance the diversity of the quest database of highly accurate excitation energies [Véril, M.; et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with a singlet or doublet ground state containing a fourth-row transition metal (CuCl, CuF, CuH, ScF, ScH, ScO, ScS, TiN, ZnH, ZnO, and ZnS) are considered, and the corresponding excitation energies are computed using high-level coupled-cluster (CC) methods, namely, CC3, CCSDT, CC4, and CCSDTQ, as well as multiconfigurational methods such as CASPT2 and NEVPT2. In many cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. Based on these calculations, theoretical best estimates of the transition energies are established in both the aug-cc-pVDZ and aug-cc-pVTZ basis sets. This allows us to accurately assess the performance of the CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported.
Collapse
Affiliation(s)
- Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Franck Gam
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, F-31062 Toulouse, France
| |
Collapse
|
4
|
Förster A, van Lenthe E, Spadetto E, Visscher L. Two-Component GW Calculations: Cubic Scaling Implementation and Comparison of Vertex-Corrected and Partially Self-Consistent GW Variants. J Chem Theory Comput 2023; 19:5958-5976. [PMID: 37594901 PMCID: PMC10501001 DOI: 10.1021/acs.jctc.3c00512] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Indexed: 08/20/2023]
Abstract
We report an all-electron, atomic orbital (AO)-based, two-component (2C) implementation of the GW approximation (GWA) for closed-shell molecules. Our algorithm is based on the space-time formulation of the GWA and uses analytical continuation (AC) of the self-energy, and pair-atomic density fitting (PADF) to switch between AO and auxiliary basis. By calculating the dynamical contribution to the GW self-energy at a quasi-one-component level, our 2C-GW algorithm is only about a factor of 2-3 slower than in the scalar relativistic case. Additionally, we present a 2C implementation of the simplest vertex correction to the self-energy, the statically screened G3W2 correction. Comparison of first ionization potentials (IPs) of a set of 67 molecules with heavy elements (a subset of the SOC81 set) calculated with our implementation against results from the WEST code reveals mean absolute deviations (MAD) of around 70 meV for G0W0@PBE and G0W0@PBE0. We check the accuracy of our AC treatment by comparison to full-frequency GW calculations, which shows that in the absence of multisolution cases, the errors due to AC are only minor. This implies that the main sources of the observed deviations between both implementations are the different single-particle bases and the pseudopotential approximation in the WEST code. Finally, we assess the performance of some (partially self-consistent) variants of the GWA for the calculation of first IPs by comparison to vertical experimental reference values. G0W0@PBE0 (25% exact exchange) and G0W0@BHLYP (50% exact exchange) perform best with mean absolute deviations (MAD) of about 200 meV. Explicit treatment of spin-orbit effects at the 2C level is crucial for systematic agreement with experiment. On the other hand, eigenvalue-only self-consistent GW (evGW) and quasi-particle self-consistent GW (qsGW) significantly overestimate the IPs. Perturbative G3W2 corrections increase the IPs and therefore improve the agreement with experiment in cases where G0W0 alone underestimates the IPs. With a MAD of only 140 meV, 2C-G0W0@PBE0 + G3W2 is in best agreement with the experimental reference values.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical
Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Erik van Lenthe
- Software
for Chemistry and Materials NV, 1081 HV Amsterdam, The Netherlands
| | - Edoardo Spadetto
- Software
for Chemistry and Materials NV, 1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical
Chemistry, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
5
|
Marie A, Loos PF. A Similarity Renormalization Group Approach to Green's Function Methods. J Chem Theory Comput 2023; 19:3943-3957. [PMID: 37311565 PMCID: PMC10339683 DOI: 10.1021/acs.jctc.3c00281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Indexed: 06/15/2023]
Abstract
The family of Green's function methods based on the GW approximation has gained popularity in the electronic structure theory thanks to its accuracy in weakly correlated systems combined with its cost-effectiveness. Despite this, self-consistent versions still pose challenges in terms of convergence. A recent study [Monino and Loos J. Chem. Phys. 2022, 156, 231101.] has linked these convergence issues to the intruder-state problem. In this work, a perturbative analysis of the similarity renormalization group (SRG) approach is performed on Green's function methods. The SRG formalism enables us to derive, from first-principles, the expression of a naturally static and Hermitian form of the self-energy that can be employed in quasiparticle self-consistent GW (qsGW) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qsGW calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code.
Collapse
Affiliation(s)
- Antoine Marie
- Laboratoire de Chimie et
Physique Quantiques (UMR 5626), Université
de Toulouse, CNRS, UPS, 31400 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et
Physique Quantiques (UMR 5626), Université
de Toulouse, CNRS, UPS, 31400 Toulouse, France
| |
Collapse
|
6
|
Holzer C. Practical Post-Kohn-Sham Methods for Time-Reversal Symmetry Breaking References. J Chem Theory Comput 2023. [PMID: 37183702 DOI: 10.1021/acs.jctc.3c00156] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
The applicability of reduced scaling algorithms based on auxiliary subspace methods for the correlation energy from the random phase approximation (RPA) as well as the correlation part of the self-energy obtained from the GW method is outlined for time-reversal symmetry breaking Kohn-Sham (KS) references. The updated algorithms allow for an efficient evaluation of RPA energies and GW quasiparticle energies for molecular systems with KS references that break time-reversal symmetry. The latter occur, for example, in magnetic fields. Furthermore, KS references for relativistic open-shell molecules also break time-reversal symmetry due to the single determinant ansatz used. Errors of the updated reduced-scaling algorithms are shown to be negligible compared to reference implementations, while the overall computational scaling is reduced by 2 orders of magnitude. Ionization energies obtained from the GW approximation are shown to be robust even for the electronically complicated group of trivalent lanthanoid ions. Starting from GW quasiparticle energies, it is subsequently shown that light-matter interactions of these systems can be calculated using the Bethe-Salpeter equation (BSE). Using the combined GW-BSE method, the absorption and emission spectra of a molecular europium(III) complex can be obtained including spin-orbit coupling.
Collapse
Affiliation(s)
- Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| |
Collapse
|
7
|
Weng G, Mallarapu R, Vlček V. Embedding vertex corrections in GW self-energy: Theory, implementation, and outlook. J Chem Phys 2023; 158:144105. [PMID: 37061461 DOI: 10.1063/5.0139117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
The vertex function (Γ) within the Green's function formalism encapsulates information about all higher-order electron-electron interaction beyond those mediated by density fluctuations. Herein, we present an efficient approach that embeds vertex corrections in the one-shot GW correlation self-energy for isolated and periodic systems. The vertex-corrected self-energy is constructed through the proposed separation-propagation-recombination procedure: the electronic Hilbert space is separated into an active space and its orthogonal complement denoted as the "rest;" the active component is propagated by a space-specific effective Hamiltonian different from the rest. The vertex corrections are introduced by a rescaled time-dependent nonlocal exchange interaction. The direct Γ correction to the self-energy is further updated by adjusting the rescaling factor in a self-consistent post-processing cycle. Our embedding method is tested mainly on donor-acceptor charge-transfer systems. The embedded vertex effects consistently and significantly correct the quasiparticle energies of the gap-edge states. The fundamental gap is generally improved by 1-3 eV upon the one-shot GW approximation. Furthermore, we provide an outlook for applications of (embedded) vertex corrections in calculations of extended solids.
Collapse
Affiliation(s)
- Guorong Weng
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, USA
| | - Rushil Mallarapu
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, USA
| | - Vojtěch Vlček
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106-9510, USA
| |
Collapse
|
8
|
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
|
9
|
Monino E, Loos PF. Unphysical discontinuities, intruder states and regularization in GW methods. J Chem Phys 2022; 156:231101. [PMID: 35732525 DOI: 10.1063/5.0089317] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
By recasting the non-linear frequency-dependent GW quasiparticle equation into a linear eigenvalue problem, we explain the appearance of multiple solutions and unphysical discontinuities in various physical quantities computed within the GW approximation. Considering the GW self-energy as an effective Hamiltonian, it is shown that these issues are key signatures of strong correlation in the (N ± 1)-electron states and can be directly related to the intruder state problem. A simple and efficient regularization procedure inspired by the similarity renormalization group is proposed to avoid such issues and speed up the convergence of partially self-consistent GW calculations.
Collapse
Affiliation(s)
- Enzo Monino
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
10
|
Cho Y, Bintrim SJ, Berkelbach TC. Simplified GW/BSE Approach for Charged and Neutral Excitation Energies of Large Molecules and Nanomaterials. J Chem Theory Comput 2022; 18:3438-3446. [PMID: 35544591 DOI: 10.1021/acs.jctc.2c00087] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Inspired by Grimme's simplified Tamm-Dancoff density functional theory approach [Grimme, S. J. Chem. Phys. 2013, 138, 244104], we describe a simplified approach to excited-state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from ab initio GW calculations by only 0.2 eV. The performance of sBSE (based on the sGW input) is tested on the excitation energies of molecules in Thiel's set, for which it differs from ab initio GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms.
Collapse
Affiliation(s)
- Yeongsu Cho
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Sylvia J Bintrim
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States.,Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
11
|
Yao Y, Golze D, Rinke P, Blum V, Kanai Y. All-Electron BSE@ GW Method for K-Edge Core Electron Excitation Energies. J Chem Theory Comput 2022; 18:1569-1583. [PMID: 35138865 DOI: 10.1021/acs.jctc.1c01180] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We present an accurate computational approach to calculate absolute K-edge core electron excitation energies as measured by X-ray absorption spectroscopy. Our approach employs an all-electron Bethe-Salpeter equation (BSE) formalism based on GW quasiparticle energies (BSE@GW) using numeric atom-centered orbitals (NAOs). The BSE@GW method has become an increasingly popular method for the computation of neutral valence excitation energies of molecules. However, it was so far not applied to molecular K-edge excitation energies. We discuss the influence of different numerical approximations on the BSE@GW calculation and employ in our final setup (i) exact numeric algorithms for the frequency integration of the GW self-energy, (ii) G0W0 and BSE starting points with ∼50% of exact exchange, (iii) the Tamm-Dancoff approximation and (iv) relativistic corrections. We study the basis set dependence and convergence with common Gaussian-type orbital and NAO basis sets. We identify the importance of additional spatially confined basis functions as well as of diffuse augmenting basis functions. The accuracy of our BSE@GW method is assessed for a benchmark set of small organic molecules, previously used for benchmarking the equation-of-motion coupled cluster method [Peng et al., J. Chem. Theory Comput., 2015, 11, 4146], as well as the medium-sized dibenzothiophene (DBT) molecule. Our BSE@GW results for absolute excitation energies are in excellent agreement with the experiment, with a mean average error of only 0.63 eV for the benchmark set and with errors <1 eV for the DBT molecule.
Collapse
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Dorothea Golze
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01062 Dresden, Germany.,Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Finland
| | - Patrick Rinke
- Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Finland
| | | | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
12
|
Bintrim SJ, Berkelbach TC. Full-frequency dynamical Bethe-Salpeter equation without frequency and a study of double excitations. J Chem Phys 2022; 156:044114. [PMID: 35105075 PMCID: PMC8807000 DOI: 10.1063/5.0074434] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Accepted: 12/31/2021] [Indexed: 01/30/2023] Open
Abstract
The Bethe-Salpeter equation (BSE) that results from the GW approximation to the self-energy is a frequency-dependent (nonlinear) eigenvalue problem due to the dynamically screened Coulomb interaction between electrons and holes. The computational time required for a numerically exact treatment of this frequency dependence is O(N6), where N is the system size. To avoid the common static screening approximation, we show that the full-frequency dynamical BSE can be exactly reformulated as a frequency-independent eigenvalue problem in an expanded space of single and double excitations. When combined with an iterative eigensolver and the density fitting approximation to the electron repulsion integrals, this reformulation yields a dynamical BSE algorithm whose computational time is O(N5), which we verify numerically. Furthermore, the reformulation provides direct access to excited states with dominant double excitation character, which are completely absent in the spectrum of the statically screened BSE. We study the 21Ag state of butadiene, hexatriene, and octatetraene and find that GW/BSE overestimates the excitation energy by about 1.5-2 eV and significantly underestimates the double excitation character.
Collapse
Affiliation(s)
- Sylvia J. Bintrim
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
13
|
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
|
14
|
Armleder J, Strunk T, Symalla F, Friederich P, Enrique Olivares Peña J, Neumann T, Wenzel W, Fediai A. Computing Charging and Polarization Energies of Small Organic Molecules Embedded into Amorphous Materials with Quantum Accuracy. J Chem Theory Comput 2021; 17:3727-3738. [PMID: 34038113 DOI: 10.1021/acs.jctc.1c00036] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The ionization potential, electron affinity, and cation/anion polarization energies (IP, EA, P(+), P(-)) of organic molecules determine injection barriers, charge carriers balance, doping efficiency, and light outcoupling in organic electronics devices, such as organic light-emitting diodes (OLEDs). Computing IP and EA of isolated molecules is a common task for quantum chemistry methods. However, once molecules are embedded in an amorphous organic matrix, IP and EA values change, and accurate predictions become challenging. Here, we present a revised quantum embedding method [Friederich et al. J. Chem. Theory Comput. 2014, 10 (9), 3720-3725] that accurately predicts the dielectric permittivity and ionization potentials in three test materials, NPB, TCTA, and C60, and allows straightforward interpretation of their nature. The method paves the way toward reliable virtual screening of amorphous organic semiconductors with targeted IP/EA, polarization energies, and relative dielectric permittivity.
Collapse
Affiliation(s)
- Jonas Armleder
- Institute of Nanotechnology, Karlsruhe Institute of Technology, 76344 Karlsruhe, Germany
| | - Timo Strunk
- Nanomatch GmbH, Griesbachstraße 5, 76185 Karlsruhe, Germany
| | - Franz Symalla
- Nanomatch GmbH, Griesbachstraße 5, 76185 Karlsruhe, Germany
| | - Pascal Friederich
- Nanomatch GmbH, Griesbachstraße 5, 76185 Karlsruhe, Germany.,Institute of Theoretical Informatics, Karlsruhe Institute of Technology, 76344 Karlsruhe, Germany
| | | | - Tobias Neumann
- Nanomatch GmbH, Griesbachstraße 5, 76185 Karlsruhe, Germany
| | - Wolfgang Wenzel
- Institute of Nanotechnology, Karlsruhe Institute of Technology, 76344 Karlsruhe, Germany
| | - Artem Fediai
- Institute of Nanotechnology, Karlsruhe Institute of Technology, 76344 Karlsruhe, Germany
| |
Collapse
|
15
|
Diez-Cabanes V, Morales-García Á, Illas F, Pastore M. Understanding the Structural and Electronic Properties of Photoactive Tungsten Oxide Nanoparticles from Density Functional Theory and GW Approaches. J Chem Theory Comput 2021; 17:3462-3470. [PMID: 33966374 DOI: 10.1021/acs.jctc.1c00293] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Tungsten trioxide (WO3)-derived nanostructures have emerged recently as feasible semiconductors for photocatalytic purposes due to their visible-light harvesting that overcomes the drawbacks presented by TiO2-derived nanoparticles (NPs). However, applications are still limited by the lack of fundamental knowledge at the nanoscale due to poor understanding of the physical processes that affect their photoactivity. To fill this gap, we report here a detailed computational study using a combined density functional theory (DFT)-GW scheme to investigate the electronic structure of realistic WO3 NPs containing up to 1680 atoms. Different phases and morphologies are considered to provide reliable structure-property relationships. Upon proper benchmark of our DFT-GW methodology, we use this highly accurate approach to establish relevant rules for the design of photoactive WO3 nanostructures by pointing out the most stable morphologies at the nanoscale and the appropriate size regime at which the photoactive efficiency is enhanced.
Collapse
Affiliation(s)
- Valentin Diez-Cabanes
- Laboratoire de Physique et Chimie Théoriques (LPCT), Université de Lorraine & CNRS, UMR 7019, F-54000 Nancy, France
| | - Ángel Morales-García
- Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, c/Martí i Franqués 1-11, 08028 Barcelona, Spain
| | - Francesc Illas
- Departament de Ciència de Materials i Química Física & Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, c/Martí i Franqués 1-11, 08028 Barcelona, Spain
| | - Mariachiara Pastore
- Laboratoire de Physique et Chimie Théoriques (LPCT), Université de Lorraine & CNRS, UMR 7019, F-54000 Nancy, France
| |
Collapse
|
16
|
Holzer C, Franzke YJ, Kehry M. Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties. J Chem Theory Comput 2021; 17:2928-2947. [PMID: 33914504 DOI: 10.1021/acs.jctc.1c00203] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.
Collapse
Affiliation(s)
- Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany.,Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| | - Max Kehry
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), 76131 Karlsruhe, Germany
| |
Collapse
|
17
|
Advanced First-Principle Modeling of Relativistic Ruddlesden—Popper Strontium Iridates. APPLIED SCIENCES-BASEL 2021. [DOI: 10.3390/app11062527] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
In this review, we provide a survey of the application of advanced first-principle methods on the theoretical modeling and understanding of novel electronic, optical, and magnetic properties of the spin-orbit coupled Ruddlesden–Popper series of iridates Srn+1IrnO3n+1 (n = 1, 2, and ∞). After a brief description of the basic aspects of the adopted methods (noncollinear local spin density approximation plus an on-site Coulomb interaction (LSDA+U), constrained random phase approximation (cRPA), GW, and Bethe–Salpeter equation (BSE)), we present and discuss select results. We show that a detailed phase diagrams of the metal–insulator transition and magnetic phase transition can be constructed by inspecting the evolution of electronic and magnetic properties as a function of Hubbard U, spin–orbit coupling (SOC) strength, and dimensionality n, which provide clear evidence for the crucial role played by SOC and U in establishing a relativistic (Dirac) Mott–Hubbard insulating state in Sr2IrO4 and Sr3Ir2O7. To characterize the ground-state phases, we quantify the most relevant energy scales fully ab initio—crystal field energy, Hubbard U, and SOC constant of three compounds—and discuss the quasiparticle band structures in detail by comparing GW and LSDA+U data. We examine the different magnetic ground states of structurally similar n = 1 and n = 2 compounds and clarify that the origin of the in-plane canted antiferromagnetic (AFM) state of Sr2IrO4 arises from competition between isotropic exchange and Dzyaloshinskii–Moriya (DM) interactions whereas the collinear AFM state of Sr3Ir2O7 is due to strong interlayer magnetic coupling. Finally, we report the dimensionality controlled metal–insulator transition across the series by computing their optical transitions and conductivity spectra at the GW+BSE level from the the quasi two-dimensional insulating n = 1 and 2 phases to the three-dimensional metallic n=∞ phase.
Collapse
|
18
|
Rezaei M, Öğüt S. Photoelectron spectra of early 3d-transition metal dioxide molecular anions from GW calculations. J Chem Phys 2021; 154:094307. [PMID: 33685151 DOI: 10.1063/5.0042106] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Photoelectron spectra of early 3d-transition metal dioxide anions, ScO2 -, TiO2 -, VO2 -, CrO2 -, and MnO2 -, are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point (G0W0@PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV-0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew-Burke-Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G0W0@PBE0, the G-only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3d-transition metal dioxide anions. Adding eigenvalue self-consistency also in W, on the other hand, worsens the agreement with experiment overall. Our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.
Collapse
Affiliation(s)
- Meisam Rezaei
- Department of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, USA
| | - Serdar Öğüt
- Department of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, USA
| |
Collapse
|
19
|
Head-Marsden K, Flick J, Ciccarino CJ, Narang P. Quantum Information and Algorithms for Correlated Quantum Matter. Chem Rev 2020; 121:3061-3120. [PMID: 33326218 DOI: 10.1021/acs.chemrev.0c00620] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Johannes Flick
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Christopher J Ciccarino
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Prineha Narang
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
20
|
Förster A, Visscher L. Low-Order Scaling G0W0 by Pair Atomic Density Fitting. J Chem Theory Comput 2020; 16:7381-7399. [PMID: 33174743 PMCID: PMC7726916 DOI: 10.1021/acs.jctc.0c00693] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2020] [Indexed: 12/18/2022]
Abstract
We derive a low-scaling G0W0 algorithm for molecules using pair atomic density fitting (PADF) and an imaginary time representation of the Green's function and describe its implementation in the Slater type orbital (STO)-based Amsterdam density functional (ADF) electronic structure code. We demonstrate the scalability of our algorithm on a series of water clusters with up to 432 atoms and 7776 basis functions and observe asymptotic quadratic scaling with realistic threshold qualities controlling distance effects and basis sets of triple-ζ (TZ) plus double polarization quality. Also owing to a very small prefactor, a G0W0 calculation for the largest of these clusters takes only 240 CPU hours with these settings. We assess the accuracy of our algorithm for HOMO and LUMO energies in the GW100 database. With errors of 0.24 eV for HOMO energies on the quadruple-ζ level, our implementation is less accurate than canonical all-electron implementations using the larger def2-QZVP GTO-type basis set. Apart from basis set errors, this is related to the well-known shortcomings of the GW space-time method using analytical continuation techniques as well as to numerical issues of the PADF approach of accurately representing diffuse atomic orbital (AO) products. We speculate that these difficulties might be overcome by using optimized auxiliary fit sets with more diffuse functions of higher angular momenta. Despite these shortcomings, for subsets of medium and large molecules from the GW5000 database, the error of our approach using basis sets of TZ and augmented double-ζ (DZ) quality is decreasing with system size. On the augmented DZ level, we reproduce canonical, complete basis set limit extrapolated reference values with an accuracy of 80 meV on average for a set of 20 large organic molecules. We anticipate our algorithm, in its current form, to be very useful in the study of single-particle properties of large organic systems such as chromophores and acceptor molecules.
Collapse
Affiliation(s)
- Arno Förster
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
21
|
Cappellini G, Bosin A, Serra G, Furthmüller J, Bechstedt F, Botti S. Electronic and Optical Properties of Small Metal Fluoride Clusters. ACS OMEGA 2020; 5:13268-13277. [PMID: 32548513 PMCID: PMC7288709 DOI: 10.1021/acsomega.0c01317] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/24/2020] [Accepted: 05/07/2020] [Indexed: 06/11/2023]
Abstract
We report a systematic investigation on the electronic and optical properties of the smallest stable clusters of alkaline-earth metal fluorides, namely, MgF2, CaF2, SrF2, and BaF2. For these clusters, we perform density functional theory (DFT) and time-dependent DFT (TDDFT) calculations with a localized Gaussian basis set. For each molecule ((MF2) n , n = 1-3, M = Mg, Ca, Sr, Ba), we determine a series of molecular properties, namely, ground-state energies, fragmentation energies, electron affinities, ionization energies, fundamental energy gaps, optical absorption spectra, and exciton binding energies. We compare electronic and optical properties between clusters of different sizes with the same metal atom and between clusters of the same size with different metal atoms. From this analysis, it turns out that MgF2 clusters have distinguished ground-state and excited-state properties with respect to the other fluoride molecules. Sizeable reductions of the optical onset energies and a consistent increase of excitonic effects are observed for all clusters under study with respect to the corresponding bulk systems. Possible consequences of the present results are discussed with respect to applied and fundamental research.
Collapse
Affiliation(s)
- Giancarlo Cappellini
- Department
of Physics, University of Cagliari, Cittadella
Universitaria di Monserrato, Strada Prov.le Monserrato-Sestu, Km 0.700, 09042, Monserrato (Ca), Italy
| | - Andrea Bosin
- Department
of Physics, University of Cagliari, Cittadella
Universitaria di Monserrato, Strada Prov.le Monserrato-Sestu, Km 0.700, 09042, Monserrato (Ca), Italy
| | - Giovanni Serra
- Department
of Physics, University of Cagliari, Cittadella
Universitaria di Monserrato, Strada Prov.le Monserrato-Sestu, Km 0.700, 09042, Monserrato (Ca), Italy
| | | | | | - Silvana Botti
- ETSF
and IFTO, FSU-Jena, Max Wien Platz 1, D-07743 Jena, Germany
| |
Collapse
|
22
|
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
|
23
|
Balasubramani SG, Chen GP, Coriani S, Diedenhofen M, Frank MS, Franzke YJ, Furche F, Grotjahn R, Harding ME, Hättig C, Hellweg A, Helmich-Paris B, Holzer C, Huniar U, Kaupp M, Marefat Khah A, Karbalaei Khani S, Müller T, Mack F, Nguyen BD, Parker SM, Perlt E, Rappoport D, Reiter K, Roy S, Rückert M, Schmitz G, Sierka M, Tapavicza E, Tew DP, van Wüllen C, Voora VK, Weigend F, Wodyński A, Yu JM. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J Chem Phys 2020; 152:184107. [PMID: 32414256 PMCID: PMC7228783 DOI: 10.1063/5.0004635] [Citation(s) in RCA: 563] [Impact Index Per Article: 140.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Accepted: 04/07/2020] [Indexed: 01/30/2023] Open
Abstract
TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.
Collapse
Affiliation(s)
- Sree Ganesh Balasubramani
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Guo P Chen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| | - Michael Diedenhofen
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Yannick J Franzke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | | | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Arnim Hellweg
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Uwe Huniar
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Alireza Marefat Khah
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | | | - Thomas Müller
- Forschungszentrum Jülich, Jülich Supercomputer Centre, Wilhelm-Jonen Straße, 52425 Jülich, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Brian D Nguyen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, USA
| | - Eva Perlt
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Kevin Reiter
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Saswata Roy
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Matthias Rückert
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Marek Sierka
- TURBOMOLE GmbH, Litzenhardtstraße 19, 76135 Karlsruhe, Germany
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, USA
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstaße 1, 70569 Stuttgart, Germany
| | - Christoph van Wüllen
- Fachbereich Chemie and Forschungszentrum OPTIMAS, Technische Universität Kaiserslautern, Erwin-Schrödinger-Staße 52, 67663 Kaiserslautern, Germany
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Florian Weigend
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Artur Wodyński
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
24
|
Cazzaniga M, Cargnoni F, Penconi M, Bossi A, Ceresoli D. Ab Initio Many-Body Perturbation Theory Calculations of the Electronic and Optical Properties of Cyclometalated Ir(III) Complexes. J Chem Theory Comput 2020; 16:1188-1199. [PMID: 31860292 DOI: 10.1021/acs.jctc.9b00763] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Cyclometalated Ir(III) compounds are the preferred choice as organic emitters in organic light-emitting diodes. In practice, the presence of the transition metal surrounded by carefully designed ligands allows fine-tuning of the emission frequency as well as good efficiency of the device. To support the development of new compounds, experimental measurements are generally compared with absorption and emission spectra obtained from ab initio calculations. The standard approach for these calculations is time-dependent density functional theory (TDDFT) with a hybrid exchange-correlation functional like B3LYP. Because of the size of these compounds, the application of more complex quantum chemistry approaches can be challenging. In this work, we used many-body perturbation theory approaches, in particular the GW approximation with the Bethe-Salpeter equation (BSE) implemented in Gaussian basis sets, to calculate the quasiparticle properties and the absorption spectra of six cyclometalated Ir(III) complexes, going beyond TDDFT. In the presented results, we compared standard TDDFT simulations with BSE calculations performed on top of perturbative G0W0 and accounting for eigenvalue self-consistency. Moreover, in order to investigate in detail the effect of the DFT starting point, we concentrated on Ir(ppy)3 and performed GW-BSE simulations starting from different DFT exchange-correlation potentials.
Collapse
Affiliation(s)
- Marco Cazzaniga
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Fausto Cargnoni
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Marta Penconi
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Alberto Bossi
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| | - Davide Ceresoli
- Consiglio Nazionale delle Ricerche , Istituto di Scienze e Tecnologie Chimiche (CNR-SCITEC) , 20133 Milano , Italy
| |
Collapse
|
25
|
Byun YM, Öğüt S. Practical GW scheme for electronic structure of 3d-transition-metal monoxide anions: ScO -, TiO -, CuO -, and ZnO . J Chem Phys 2019; 151:134305. [PMID: 31594362 DOI: 10.1063/1.5118671] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The GW approximation to many-body perturbation theory is a reliable tool for describing charged electronic excitations, and it has been successfully applied to a wide range of extended systems for several decades using a plane-wave basis. However, the GW approximation has been used to test limited spectral properties of a limited set of finite systems (e.g., frontier orbital energies of closed-shell sp molecules) only for about a decade using a local-orbital basis. Here, we calculate the quasiparticle spectra of closed- and open-shell molecular anions with partially and completely filled 3d shells (shallow and deep 3d states, respectively), ScO-, TiO-, CuO-, and ZnO-, using various levels of GW theory, and compare them to experiments to evaluate the performance of the GW approximation on the electronic structure of small molecules containing 3d transition metals. We find that the G-only eigenvalue self-consistent GW scheme with W fixed to the PBE level (GnW0@PBE), which gives the best compromise between accuracy and efficiency for solids, also gives good results for both localized (d) and delocalized (sp) states of 3d-transition-metal oxide molecules. The success of GnW0@PBE in predicting electronic excitations in these systems reasonably well is likely due to the fortuitous cancellation effect between the overscreening of the Coulomb interaction by PBE and the underscreening by the neglect of vertex corrections. Together with the absence of the self-consistent field convergence error (e.g., spin contamination in open-shell systems) and the GW multisolution issue, the GnW0@PBE scheme gives the possibility to predict the electronic structure of complex real systems (e.g., molecule-solid and sp-d hybrid systems) accurately and efficiently.
Collapse
Affiliation(s)
- Young-Moo Byun
- Department of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, USA
| | - Serdar Öğüt
- Department of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, USA
| |
Collapse
|
26
|
Koval P, Ljungberg MP, Müller M, Sánchez-Portal D. Toward Efficient GW Calculations Using Numerical Atomic Orbitals: Benchmarking and Application to Molecular Dynamics Simulations. J Chem Theory Comput 2019; 15:4564-4580. [PMID: 31318555 DOI: 10.1021/acs.jctc.9b00436] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The use of atomic orbitals in Hedin's GW approximation provides, in principle, an inexpensive alternative to plane-wave basis sets, especially when modeling large molecules. However, benchmarking of the algorithms and basis sets is essential for a careful balance between cost and accuracy. In this paper, we present an implementation of the GW approximation using numerical atomic orbitals and a pseudopotential treatment of core electrons. The combination of a contour deformation technique with a one-shot extraction of quasiparticle energies provides an efficient scheme for many applications. The performance of the implementation with respect to the basis set convergence and the effect of the use of pseudopotentials has been tested for the 117 closed-shell molecules from the G2/97 test set and 24 larger acceptor molecules from another recently proposed test set. Moreover, to demonstrate the potential of our method, we compute the thermally averaged GW density of states of a large photochromic compound by sampling ab initio molecular dynamics trajectories at different temperatures. The computed thermal line widths indicate approximately twice as large electron-phonon couplings with GW than with standard DFT-GGA calculations. This is further confirmed using frozen-phonon calculations.
Collapse
Affiliation(s)
- Peter Koval
- Donostia International Physics Center , Paseo Manuel de Lardizabal 4 , 20018 Donostia-San Sebastián , Spain
| | - Mathias Per Ljungberg
- Donostia International Physics Center , Paseo Manuel de Lardizabal 4 , 20018 Donostia-San Sebastián , Spain
| | - Moritz Müller
- Donostia International Physics Center , Paseo Manuel de Lardizabal 4 , 20018 Donostia-San Sebastián , Spain
| | - Daniel Sánchez-Portal
- Donostia International Physics Center , Paseo Manuel de Lardizabal 4 , 20018 Donostia-San Sebastián , Spain.,Centro de Física de Materiales , Centro Mixto CSIC-UPV/EHU , Paseo Manuel de Lardizabal 5 , 20018 Donostia-San Sebastián , Spain
| |
Collapse
|
27
|
Golze D, Dvorak M, Rinke P. The GW Compendium: A Practical Guide to Theoretical Photoemission Spectroscopy. Front Chem 2019; 7:377. [PMID: 31355177 PMCID: PMC6633269 DOI: 10.3389/fchem.2019.00377] [Citation(s) in RCA: 153] [Impact Index Per Article: 30.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2018] [Accepted: 05/08/2019] [Indexed: 12/22/2022] Open
Abstract
The GW approximation in electronic structure theory has become a widespread tool for predicting electronic excitations in chemical compounds and materials. In the realm of theoretical spectroscopy, the GW method provides access to charged excitations as measured in direct or inverse photoemission spectroscopy. The number of GW calculations in the past two decades has exploded with increased computing power and modern codes. The success of GW can be attributed to many factors: favorable scaling with respect to system size, a formal interpretation for charged excitation energies, the importance of dynamical screening in real systems, and its practical combination with other theories. In this review, we provide an overview of these formal and practical considerations. We expand, in detail, on the choices presented to the scientist performing GW calculations for the first time. We also give an introduction to the many-body theory behind GW, a review of modern applications like molecules and surfaces, and a perspective on methods which go beyond conventional GW calculations. This review addresses chemists, physicists and material scientists with an interest in theoretical spectroscopy. It is intended for newcomers to GW calculations but can also serve as an alternative perspective for experts and an up-to-date source of computational techniques.
Collapse
Affiliation(s)
- Dorothea Golze
- Department of Applied Physics, Aalto University, School of Science, Espoo, Finland
| | | | | |
Collapse
|
28
|
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
|
29
|
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
|
30
|
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
|
31
|
Vogiatzis KD, Polynski MV, Kirkland JK, Townsend J, Hashemi A, Liu C, Pidko EA. Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities. Chem Rev 2019; 119:2453-2523. [PMID: 30376310 PMCID: PMC6396130 DOI: 10.1021/acs.chemrev.8b00361] [Citation(s) in RCA: 222] [Impact Index Per Article: 44.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Indexed: 12/28/2022]
Abstract
Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.
Collapse
Affiliation(s)
| | | | - Justin K. Kirkland
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Jacob Townsend
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ali Hashemi
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Chong Liu
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Evgeny A. Pidko
- TheoMAT
group, ITMO University, Lomonosova 9, St. Petersburg 191002, Russia
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| |
Collapse
|
32
|
Ziaei V, Bredow T. Screening mixing GW/Bethe-Salpeter approach for triplet states of organic molecules. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:395501. [PMID: 30124435 DOI: 10.1088/1361-648x/aadb75] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
We recently proposed a screening mixing many-body ansatz (Ziaei and Bredow 2017 Phys. Rev. B 96 195115) to decrease the typical overestimation of the Hartree-Fock based GW/Bethe-Salpeter equation (BSE) singlet-singlet excitation energies in molecular systems. Now we have evaluated the accuracy of the proposed scheme for triplet states of a set of 20 organic molecules known as the Thiel set. We show that by mixing different screenings into GW and BSE calculated within random phase approximation (in order to ensure best gap and an optimal exciton binding energy), the total mean absolute error of 0.59 eV in the standard Hartree-Fock based eigenvalue GW/BSE approach is reduced to 0.26 eV for 63 triplet states. We further demonstrate that the quasi-particle self-consistent GW/BSE approach in which orbitals and energies are updated in the Green's function and the dynamically screened interaction mostly and considerably underestimates the excitation energies as shown for a few molecules.
Collapse
Affiliation(s)
- Vafa Ziaei
- Mulliken Center for Theoretical Chemistry, University of Bonn, 53115 Bonn, Germany
| | | |
Collapse
|
33
|
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
|
34
|
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
|
35
|
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
|
36
|
Hashemi Z, Rafiezadeh S, Hafizi R, Hashemifar SJ, Akbarzadeh H. First-principles study of MoS2 and MoSe2 nanoclusters in the framework of evolutionary algorithm and density functional theory. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.03.008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
37
|
Colonna N, Nguyen NL, Ferretti A, Marzari N. Screening in Orbital-Density-Dependent Functionals. J Chem Theory Comput 2018; 14:2549-2557. [DOI: 10.1021/acs.jctc.7b01116] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Nicola Colonna
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Ngoc Linh Nguyen
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | | | - Nicola Marzari
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
38
|
Wilhelm J, Golze D, Talirz L, Hutter J, Pignedoli CA. Toward GW Calculations on Thousands of Atoms. J Phys Chem Lett 2018; 9:306-312. [PMID: 29280376 DOI: 10.1021/acs.jpclett.7b02740] [Citation(s) in RCA: 71] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
The GW approximation of many-body perturbation theory is an accurate method for computing electron addition and removal energies of molecules and solids. In a canonical implementation, however, its computational cost is [Formula: see text] in the system size N, which prohibits its application to many systems of interest. We present a full-frequency GW algorithm in a Gaussian-type basis, whose computational cost scales with N2 to N3. The implementation is optimized for massively parallel execution on state-of-the-art supercomputers and is suitable for nanostructures and molecules in the gas, liquid or condensed phase, using either pseudopotentials or all electrons. We validate the accuracy of the algorithm on the GW100 molecular test set, finding mean absolute deviations of 35 meV for ionization potentials and 27 meV for electron affinities. Furthermore, we study the length-dependence of quasiparticle energies in armchair graphene nanoribbons of up to 1734 atoms in size, and compute the local density of states across a nanoscale heterojunction.
Collapse
Affiliation(s)
- Jan Wilhelm
- Department of Chemistry, University of Zurich , Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Dorothea Golze
- COMP/Department of Applied Physics, Aalto University , P.O. Box 11100, FI-00076 Aalto, Finland
| | - Leopold Talirz
- Laboratory of Molecular Simulation, École Polytechnique Fédérale de Lausanne , Rue de l'Industrie 17, CH-1951 Sion, Switzerland
- Theory and Simulation of Materials, École Polytechnique Fédérale de Lausanne , Station 9, CH-1015 Lausanne, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich , Winterthurerstrasse 190, CH-8057 Zurich, Switzerland
| | - Carlo A Pignedoli
- Swiss Federal Laboratories for Materials Science and Technology (Empa) , Überlandstrasse 129, CH-8600 Dübendorf, Switzerland
| |
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
|
Reining L. The GW approximation: content, successes and limitations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1344] [Citation(s) in RCA: 85] [Impact Index Per Article: 12.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Lucia Reining
- Laboratoire des Solides Irradiés, École Polytechnique, CNRS, CEAUniversité Paris‐Saclay Palaiseau France
- European Theoretical Spectroscopy Facility (ETSF), http:///www.etsf.eu
| |
Collapse
|
41
|
Maggio E, Kresse G. GW Vertex Corrected Calculations for Molecular Systems. J Chem Theory Comput 2017; 13:4765-4778. [DOI: 10.1021/acs.jctc.7b00586] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Emanuele Maggio
- Faculty of Physics
and Center
for Computational Materials Science, University of Vienna, Sensengasse
8/12, A-1090 Vienna, Austria
| | - Georg Kresse
- Faculty of Physics
and Center
for Computational Materials Science, University of Vienna, Sensengasse
8/12, A-1090 Vienna, Austria
| |
Collapse
|
42
|
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
|
43
|
Atkins AJ, Talotta F, Freitag L, Boggio-Pasqua M, González L. Assessing Excited State Energy Gaps with Time-Dependent Density Functional Theory on Ru(II) Complexes. J Chem Theory Comput 2017; 13:4123-4145. [DOI: 10.1021/acs.jctc.7b00379] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Andrew J. Atkins
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Francesco Talotta
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
- Laboratoire de
Chimie et Physique Quantiques (UMR5626), CNRS et Université
de Toulouse 3, 31062 Toulouse, France
| | - Leon Freitag
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Martial Boggio-Pasqua
- Laboratoire de
Chimie et Physique Quantiques (UMR5626), CNRS et Université
de Toulouse 3, 31062 Toulouse, France
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| |
Collapse
|
44
|
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
|
45
|
Bruneval F. Optimized virtual orbital subspace for faster GW calculations in localized basis. J Chem Phys 2017; 145:234110. [PMID: 28010103 DOI: 10.1063/1.4972003] [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/14/2022] Open
Abstract
The popularity of the GW approximation to the self-energy to access the quasiparticle energies of molecules is constantly increasing. As the other methods addressing the electronic correlation, the GW self-energy unfortunately shows a very slow convergence with respect to the basis complexity, which precludes the calculation of accurate quasiparticle energies for large molecules. Here we propose a method to mitigate this issue that relies on two steps: (i) the definition of a reduced virtual orbital subspace, thanks to a much smaller basis set; (ii) the account of the remainder through the simpler one-ring approximation to the self-energy. We assess the quality of the corrected quasiparticle energies for simple molecules, and finally we show an application to large graphene chunks to demonstrate the numerical efficiency of the scheme.
Collapse
Affiliation(s)
- Fabien Bruneval
- CEA, DEN, Service de Recherches de Métallurgie Physique, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France
| |
Collapse
|
46
|
Hung L, Bruneval F, Baishya K, Öğüt S. Benchmarking the GW Approximation and Bethe–Salpeter Equation for Groups IB and IIB Atoms and Monoxides. J Chem Theory Comput 2017; 13:2135-2146. [DOI: 10.1021/acs.jctc.7b00123] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Linda Hung
- NIST
Center for Neutron Research, National Institute of Standard and Technology, Gaithersburg, Maryland 20899, United States
| | - Fabien Bruneval
- CEA,
DEN, Service de Recherches de Métallurgie Physique, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France
| | - Kopinjol Baishya
- Department
of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, United States
| | - Serdar Öğüt
- Department
of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, United States
| |
Collapse
|
47
|
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
|
48
|
Marom N. Accurate description of the electronic structure of organic semiconductors by GW methods. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:103003. [PMID: 28145283 DOI: 10.1088/1361-648x/29/10/103003] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Electronic properties associated with charged excitations, such as the ionization potential (IP), the electron affinity (EA), and the energy level alignment at interfaces, are critical parameters for the performance of organic electronic devices. To computationally design organic semiconductors and functional interfaces with tailored properties for target applications it is necessary to accurately predict these properties from first principles. Many-body perturbation theory is often used for this purpose within the GW approximation, where G is the one particle Green's function and W is the dynamically screened Coulomb interaction. Here, the formalism of GW methods at different levels of self-consistency is briefly introduced and some recent applications to organic semiconductors and interfaces are reviewed.
Collapse
Affiliation(s)
- Noa Marom
- Department of Materials Science and Engineering, Department of Chemistry, and Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, United States of America
| |
Collapse
|
49
|
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
|
50
|
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
|