1
|
Bhattacharya S, Li J, Yang W, Kanai Y. BSE@ GW Prediction of Charge Transfer Exciton in Molecular Complexes: Assessment of Self-Energy and Exchange-Correlation Dependence. J Phys Chem A 2024; 128:6072-6083. [PMID: 39011742 DOI: 10.1021/acs.jpca.4c02898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/17/2024]
Abstract
The Bethe-Salpeter equation using the GW approximation to the self-energy (BSE@GW) is a computationally attractive method for studying electronic excitation from first principles within the many-body Green's function theory framework. We examine its dependence on the underlying exchange-correlation (XC) approximation as well as on the GW approximation for predicting the charge transfer exciton formation at representative type-II interfaces between molecular systems of tetrachloro-1,2-benzoquinone (TCBQ) and acene derivatives. For the XC approximation, we consider several widely used generalized gradient approximation (GGA) and hybrid GGA functionals. For the GW self-energy approximation, we examine the recently proposed renormalized singles approach by Yang and coauthors [J. Phys. Chem. Lett. 2019, 10 (3), 447-452; J. Chem. Theory Comput. 2022, 18, 7570-7585] in addition to other commonly employed approximated GW schemes. We demonstrate a reliable prediction of the charge transfer exciton within the BSE@GW level of theory.
Collapse
Affiliation(s)
- Sampreeti Bhattacharya
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, United States
| | - Jiachen Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, United States
- Department of Physics and Astronomy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27514, United States
| |
Collapse
|
2
|
Kossoski F, Boggio-Pasqua M, Loos PF, Jacquemin D. Reference Energies for Double Excitations: Improvement and Extension. J Chem Theory Comput 2024; 20:5655-5678. [PMID: 38885174 DOI: 10.1021/acs.jctc.4c00410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
In the realm of photochemistry, the significance of double excitations (also known as doubly excited states), where two electrons are concurrently elevated to higher energy levels, lies in their involvement in key electronic transitions essential in light-induced chemical reactions as well as their challenging nature from the computational theoretical chemistry point of view. Based on state-of-the-art electronic structure methods (such as high-order coupled-cluster, selected configuration interaction, and multiconfigurational methods), we improve and expand our prior set of accurate reference excitation energies for electronic states exhibiting a substantial amount of double excitations [Loos et al. J. Chem. Theory Comput. 2019, 15, 1939]. This extended collection encompasses 47 electronic transitions across 26 molecular systems that we separate into two distinct subsets: (i) 28 "genuine" doubly excited states where the transitions almost exclusively involve doubly excited configurations and (ii) 19 "partial" doubly excited states which exhibit a more balanced character between singly and doubly excited configurations. For each subset, we assess the performance of high-order coupled-cluster (CC3, CCSDT, CC4, and CCSDTQ) and multiconfigurational methods (CASPT2, CASPT3, PC-NEVPT2, and SC-NEVPT2). Using as a probe the percentage of single excitations involved in a given transition (%T1) computed at the CC3 level, we also propose a simple correction that reduces the errors of CC3 by a factor of 3, for both sets of excitations. We hope that this more complete and diverse compilation of double excitations will help future developments of electronic excited-state methodologies.
Collapse
Affiliation(s)
- Fábris Kossoski
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Martial Boggio-Pasqua
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31062 Toulouse, France
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France (IUF), F-75005 Paris, France
| |
Collapse
|
3
|
Marie A, Loos PF. Reference Energies for Valence Ionizations and Satellite Transitions. J Chem Theory Comput 2024; 20:4751-4777. [PMID: 38776293 PMCID: PMC11171335 DOI: 10.1021/acs.jctc.4c00216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 04/10/2024] [Accepted: 04/11/2024] [Indexed: 05/24/2024]
Abstract
Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shakeup transitions) as they appear in ionization spectra as higher-energy peaks with weaker intensity and larger width than the main peaks associated with single-particle ionizations. Satellites, which correspond to electronically excited states of the cationic species, are notoriously challenging to model using conventional single-reference methods due to their high excitation degree compared to the neutral reference state. This work reports 42 satellite transition energies and 58 valence ionization potentials (IPs) of full configuration interaction quality computed in small molecular systems. Following the protocol developed for the quest database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; and Loos, P.-F. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2021, 11, e1517], these reference energies are computed using the configuration interaction using a perturbative selection made iteratively (CIPSI) method. In addition, the accuracy of the well-known coupled-cluster (CC) hierarchy (CC2, CCSD, CC3, CCSDT, CC4, and CCSDTQ) is gauged against these new accurate references. The performances of various approximations based on many-body Green's functions (GW, GF2, and T-matrix) for IPs are also analyzed. Their limitations in correctly modeling satellite transitions are discussed.
Collapse
Affiliation(s)
- Antoine Marie
- Laboratoire de Chimie et Physique
Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse 31062, France
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique
Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse 31062, France
| |
Collapse
|
4
|
Li J, Jin Y, Yu J, Yang W, Zhu T. Accurate Excitation Energies of Point Defects from Fast Particle-Particle Random Phase Approximation Calculations. J Phys Chem Lett 2024:2757-2764. [PMID: 38436573 DOI: 10.1021/acs.jpclett.4c00184] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2024]
Abstract
We present an efficient particle-particle random phase approximation (ppRPA) approach that predicts accurate excitation energies of point defects, including the nitrogen-vacancy (NV-) and silicon-vacancy (SiV0) centers in diamond and the divacancy center (VV0) in 4H silicon carbide, with errors of ±0.2 eV compared with experimental values. Starting from the (N + 2)-electron ground state calculated with density functional theory (DFT), the ppRPA excitation energies of the N-electron system are calculated as the differences between the two-electron removal energies of the (N + 2)-electron system. We demonstrate that the ppRPA excitation energies converge rapidly with a few hundred canonical active-space orbitals. We also show that active-space ppRPA has weak DFT starting-point dependence and is significantly cheaper than the corresponding ground-state DFT calculation. This work establishes ppRPA as an accurate and low-cost tool for investigating excited-state properties of point defects and opens up new opportunities for applications of ppRPA to periodic bulk materials.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Yu Jin
- Department of Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Jincheng Yu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Tianyu Zhu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| |
Collapse
|
5
|
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
|
6
|
Li J, Yu J, Chen Z, Yang W. Linear Scaling Calculations of Excitation Energies with Active-Space Particle-Particle Random-Phase Approximation. J Phys Chem A 2023; 127:7811-7822. [PMID: 37695567 DOI: 10.1021/acs.jpca.3c02834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/12/2023]
Abstract
We developed an efficient active-space particle-particle random-phase approximation (ppRPA) approach to calculate accurate charge-neutral excitation energies of molecular systems. The active-space ppRPA approach constrains both indexes in particle and hole pairs in the ppRPA matrix, which only selects frontier orbitals with dominant contributions to low-lying excitation energies. It employs the truncation in both orbital indexes in the particle-particle and the hole-hole spaces. The resulting matrix, whose eigenvalues are excitation energies, has a dimension that is independent of the size of the systems. The computational effort for the excitation energy calculation, therefore, scales linearly with system size and is negligible compared with the ground-state calculation of the (N - 2)-electron system, where N is the electron number of the molecule. With the active space consisting of 30 occupied and 30 virtual orbitals, the active-space ppRPA approach predicts the excitation energies of valence, charge-transfer, Rydberg, double, and diradical excitations with the mean absolute errors (MAEs) smaller than 0.03 eV compared with the full-space ppRPA results. As a side product, we also applied the active-space ppRPA approach in the renormalized singles (RS) T-matrix approach. Combining the non-interacting pair approximation that approximates the contribution to the self-energy outside the active space, the active-space GRSTRS@PBE approach predicts accurate absolute and relative core-level binding energies with the MAEs around 1.58 and 0.3 eV, respectively. The developed linear scaling calculation of excitation energies is promising for applications to large and complex systems.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Jincheng Yu
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Zehua Chen
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
7
|
Marie A, Loos PF. A Similarity Renormalization Group Approach to Green's Function Methods. J Chem Theory Comput 2023. [PMID: 37311565 DOI: 10.1021/acs.jctc.3c00281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [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
|
8
|
Orlando R, Romaniello P, Loos PF. Exploring new exchange-correlation kernels in the Bethe–Salpeter equation: A study of the asymmetric Hubbard dimer. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.02.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
|
9
|
Quintero-Monsebaiz R, Monino E, Marie A, Loos PF. Connections between many-body perturbation and coupled-cluster theories. J Chem Phys 2022; 157:231102. [PMID: 36550046 DOI: 10.1063/5.0130837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
Here, we build on the works of Scuseria et al. [J. Chem. Phys. 129, 231101 (2008)] and Berkelbach [J. Chem. Phys. 149, 041103 (2018)] to show connections between the Bethe-Salpeter equation (BSE) formalism combined with the GW approximation from many-body perturbation theory and coupled-cluster (CC) theory at the ground- and excited-state levels. In particular, we show how to recast the GW and Bethe-Salpeter equations as non-linear CC-like equations. Similitudes between BSE@GW and the similarity-transformed equation-of-motion CC method are also put forward. The present work allows us to easily transfer key developments and the general knowledge gathered in CC theory to many-body perturbation theory. In particular, it may provide a path for the computation of ground- and excited-state properties (such as nuclear gradients) within the GW and BSE frameworks.
Collapse
Affiliation(s)
- Raúl Quintero-Monsebaiz
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Enzo Monino
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| | - Antoine Marie
- 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
|
Li J, Golze D, Yang W. Combining Renormalized Singles GW Methods with the Bethe-Salpeter Equation for Accurate Neutral Excitation Energies. J Chem Theory Comput 2022; 18:6637-6645. [PMID: 36279250 PMCID: PMC9972216 DOI: 10.1021/acs.jctc.2c00686] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We apply the renormalized singles (RS) Green's function in the Bethe-Salpeter equation (BSE)/GW approach to predict accurate neutral excitation energies of molecular systems. The BSE calculations are performed on top of the GRSWRS method, which uses the RS Green's function also for the computation of the screened Coulomb interaction W. We show that the BSE/GRSWRS approach significantly outperforms BSE/G0W0 for predicting excitation energies of valence, Rydberg, and charge-transfer (CT) excitations by benchmarking the Truhlar-Gagliardi set, Stein CT set, and an atomic Rydberg test set. For the Truhlar-Gagliardi test set, BSE/GRSWRS provides comparable accuracy to time-dependent density functional theory (TDDFT) and is slightly better than BSE starting from eigenvalue self-consistent GW (evGW). For the Stein CT test set, BSE/GRSWRS significantly outperforms BSE/G0W0 and TDDFT with the accuracy comparable to BSE/evGW. We also show that BSE/GRSWRS predicts Rydberg excitation energies of atomic systems well. Besides the excellent accuracy, BSE/GRSWRS largely eliminates the dependence on the choice of the density functional approximation. This work demonstrates that the BSE/GRSWRS approach is accurate and efficient for predicting excitation energies for a broad range of systems, which expands the applicability of the BSE/GW approach.
Collapse
Affiliation(s)
- Jiachen Li
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| | - Dorothea Golze
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01062 Dresden, Germany
| | - Weitao Yang
- Department of Chemistry, Duke University, Durham, North Carolina 27708, United States
| |
Collapse
|
11
|
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
|