1
|
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
|
2
|
Damour Y, Scemama A, Jacquemin D, Kossoski F, Loos PF. State-Specific Coupled-Cluster Methods for Excited States. J Chem Theory Comput 2024; 20:4129-4145. [PMID: 38749498 PMCID: PMC11137840 DOI: 10.1021/acs.jctc.4c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Revised: 03/06/2024] [Accepted: 03/06/2024] [Indexed: 05/29/2024]
Abstract
We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.
Collapse
Affiliation(s)
- Yann Damour
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Anthony Scemama
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Denis Jacquemin
- Nantes
Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut
Universitaire de France (IUF), F-75005 Paris, France
| | - Fábris Kossoski
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| | - Pierre-François Loos
- Laboratoire
de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, 31000 Toulouse, France
| |
Collapse
|
3
|
Chaussy L, Chilkuri VG, Humbel S, Nava P. Spectroscopy of End-On Copper(II) Superoxido Complexes: A Wave Function-Based Analysis. Inorg Chem 2024; 63:8038-8049. [PMID: 38659336 DOI: 10.1021/acs.inorgchem.3c04401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Wave function methods are employed to analyze the ground and low-lying excited states of bipyramid trigonal copper(II) superoxido complexes, up to their characteristic ligand to metal charge transfer band. Several multireference methods have been combined to provide new insights into the interpretation of their experimental absorption spectra. We show that the intraligand transition on the dioxygen leads to a dark state. Among the results, we shall highlight the finding of doubly excited states in the region of the d-d transitions and the subtle interplay between Cu(I) and Cu(II) in the ground and excited states. Some of these findings could be obtained only with multireference methods.
Collapse
Affiliation(s)
- Léo Chaussy
- Aix Marseille Univ, CNRS, Centrale Marseille, iSm2, Marseille 13397, France
| | | | - Stéphane Humbel
- Aix Marseille Univ, CNRS, Centrale Marseille, iSm2, Marseille 13397, France
| | - Paola Nava
- Aix Marseille Univ, CNRS, Centrale Marseille, iSm2, Marseille 13397, France
| |
Collapse
|
4
|
Schraivogel T, Kats D. Two determinant distinguishable cluster. J Chem Phys 2024; 160:124109. [PMID: 38526108 DOI: 10.1063/5.0199274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Accepted: 03/06/2024] [Indexed: 03/26/2024] Open
Abstract
A two reference determinant version of the distinguishable cluster with singles and doubles (DCSD) has been developed. We have implemented the two determinant distinguishable cluster (2D-DCSD) and the corresponding traditional 2D-CCSD method in a new open-source package written in Julia called ElemCo.jl. The methods were benchmarked on singlet and triplet excited states of valence and Rydberg character, as well as for singlet-triplet gaps of diradicals. It is demonstrated that the distinguishable cluster approximation improves the accuracy of 2D-CCSD.
Collapse
Affiliation(s)
- Thomas Schraivogel
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Daniel Kats
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| |
Collapse
|
5
|
Ponra A, Bakasa C, Etindele AJ, Casida ME. Diagrammatic multiplet sum method (MSM) density functional theory (DFT): Investigation of the transferability of integrals in "simple" DFT-based approaches to multideterminantal problems. J Chem Phys 2023; 159:244306. [PMID: 38149739 DOI: 10.1063/5.0173572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 11/30/2023] [Indexed: 12/28/2023] Open
Abstract
Kohn-Sham density functional theory (DFT) typically works well for describing dynamic correlation. Two other types of correlation, arising in the cases of degenerate (static) or quasidegenerate (nondynamic) zero-order states, represent a difficult problem for DFT. When symmetry is present, multiplet sum method (MSM) DFT [Ziegler et al., Theor. Chim. Acta 4, 877 (1977)] provides one of the earliest and simplest ways to include static correlation in DFT. MSM-DFT assumes that DFT provides a good description of single-determinant energies and uses symmetry and simple ansätze to include the effects of static correlation. This is equivalent to determining the off-diagonal matrix elements in a small configuration interaction (CI) eigenvalue problem. Our ultimate goal, however, is nondynamic correlation in cases where symmetry is inadequate for fixing the dynamic-correlation limitation of DFT. To this end, we have developed a diagrammatic approach to MSM-DFT, which does not, by itself, solve the nondynamic correlation problem in DFT but which facilitates comparison with wave function CI and so allows educated guesses of off-diagonal CI matrix elements even in the absence of symmetry. In every case, an additional exchange-only ansatz (EXAN) allows the MSM-DFT formulas to be transformed into wave function formulas. This EXAN also works for transforming time-dependent DFT into time-dependent Hartree-Fock. Although not enough to uniquely guess DFT formulas from wave function formulas, the diagrammatic approach and the EXAN provide important constraints on any guesses that might be used. We illustrate how diagrammatic MSM-DFT may be used to guess a nondynamic correlation correction for the dissociation of H2 and how diagrammatic MSM-DFT may be used to guess a nonsymmetry-based coupling element in the O2 multiplet problem, which is reasonably close to a previous symmetry-derived result.
Collapse
Affiliation(s)
- Abraham Ponra
- Department of Physics, Faculty of Science, University of Maroua, P.O. Box 814, Maroua, Cameroon
| | - Carolyne Bakasa
- Technical University of Kenya, P.O. Box 52428-00200, Haile Selassie Avenue, Nairobe, Kenya
| | - Anne Justine Etindele
- Higher Teachers Training College, University of Yaounde I, P.O. Box 47, Yaounde, Cameroon
| | - Mark E Casida
- Laboratoire de Spectrométrie, Interactions et Chimie théorique (SITh), Département de Chimie Moléculaire (DCM, UMR CNRS/UGA 5250), Institut de Chimie Moléculaire de Grenoble (ICMG, FR2607), Université Grenoble Alpes (UGA) 301 rue de la Chimie, BP 53, F-38041 Grenoble Cedex 9, France
| |
Collapse
|
6
|
Ravi M, Perera A, Park YC, Bartlett RJ. Excited states with pair coupled cluster doubles tailored coupled cluster theory. J Chem Phys 2023; 159:094101. [PMID: 37655762 DOI: 10.1063/5.0161368] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Accepted: 08/11/2023] [Indexed: 09/02/2023] Open
Abstract
It is known that some non-dynamic effects of electron correlation can be included in coupled cluster theory using a tailoring technique that separates the effects of non-dynamic and dynamic correlations. Recently, the simple pCCD (pair coupled cluster doubles) wavefunction was shown to provide good results for some non-dynamic correlation problems, such as bond-breaking, in a spin-adapted way with no active space selection. In this paper, we report a study of excited states using "tailored coupled cluster singles and doubles," to attempt to use pCCD as a kernel for more complete coupled-cluster singles and doubles (CCSD) results for excited states. Several excited states are explored from those primarily due to single excitations to those dominated by doubly excited states and from singlet-triplet splittings for some diradical states. For the first two situations, tailored pCCD-TCCSD offers no improvement over equation of motion-CCSD. However, when we explore the singlet-triplet gap of diradical molecules that are manifestly multi-reference, a pCCD kernel provides improved results, particularly with generalized valence bond orbitals.
Collapse
Affiliation(s)
- Moneesha Ravi
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Young Choon Park
- Korea Institute of Fusion Energy, 37 Dongjangsan-ro, Gunsan, Jeollabuk-do 54004, Republic of Korea
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
7
|
Rishi V, Ravi M, Perera A, Bartlett RJ. Dark Doubly Excited States with Modified Coupled Cluster Models: A Reliable Compromise between Cost and Accuracy? J Phys Chem A 2023; 127:828-834. [PMID: 36640093 DOI: 10.1021/acs.jpca.2c07697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
To treat doubly excited states, the treatment of triple excitations is considered necessary in the framework of equation-of-motion coupled cluster (EOM-CC) methods. We investigate models without explicit triples and seek quantitative measure for the performance of EOM based on CC with singles and doubles (CCSD) or modified CCSD (Distinguishable Cluster Approximation) approaches for states with predominant double excitation character. We also test the efficacy of including triples in perturbative manner through EOM-CCSD(T) and in an iterative way through EOM-CCSDT-3 method. Extended similarity transformed EOM-CCSD(EXT-STEOM-CCSD) method is also tested and provides superior quality results at comparatively low cost. We use the QUEST2 benchmark set of double excitations proposed by Loos et al. [ J. Chem. Theory Comput.2019, 15, 1939] to investigate the performance of methods such as EOM-CCSD, EOM-DCSD, EXT-STEOM-CCSD, ΔCCSD, and ΔDCSD. We also test a tailored CC approach, ΔpairCCD-TCCSD.
Collapse
Affiliation(s)
- Varun Rishi
- Quantum Theory Project, University of Florida, Gainesville, Florida32611, United States
| | - Moneesha Ravi
- Quantum Theory Project, University of Florida, Gainesville, Florida32611, United States
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida32611, United States
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida32611, United States
| |
Collapse
|
8
|
Park YC, Perera A, Bartlett RJ. Density functionals for core excitations. J Chem Phys 2022; 157:094107. [PMID: 36075732 DOI: 10.1063/5.0111095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The core excitation energies and related principal ionization energies are obtained for selected molecules using several density functionals and compared with benchmark equation-of-motion coupled cluster (EOM-CC) results. Both time-dependent and time-independent formulations of excitation spectra in the time-dependent density functional theory and the EOM-CC are employed to obtain excited states that are not always easily accessible with the time-independent method. Among those functionals, we find that the QTP(00) functional, which is only parameterized to reproduce the five IPs of water, provides excellent core IPs and core excitation energies, consistently yielding better excitation and ionization energies. We show that orbital eigenvalues of KS density functional theory play an important role in determining the accuracy of the excitation and photoelectron spectra.
Collapse
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
9
|
Windom ZW, Perera A, Bartlett RJ. Examining fundamental and excitation gaps at the thermodynamic limit: A combined (QTP) DFT and coupled cluster study on trans-polyacetylene and polyacene. J Chem Phys 2022; 156:204308. [DOI: 10.1063/5.0086158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Interest in ab initio property prediction of π-conjugated polymers for technological applications places significant demand on “cost-effective” and conceptual computational methods, particularly effective, one-particle theories. This is particularly relevant in the case of Kohn–Sham Density Functional Theory (KS-DFT) and its new competitors that arise from correlated orbital theory, the latter defining the QTP family of DFT functionals. This study presents large, ab initio equation of motion-coupled cluster calculations using the massively parallel ACESIII to target the fundamental bandgap of two prototypical organic polymers, trans-polyacetylene (tPA) and polyacene (Ac), and provides an assessment of the new quantum theory project (QTP) functionals for this problem. Further results focusing on the 1 A g (1 A g), 1 B u (1 B2 u), and 3 B u (3 B2 u) excited states of tPA (Ac) are also presented. By performing calculations on oligomers of increasing size, extrapolations to the thermodynamic limit for the fundamental and all excitation gaps, as well as estimations of the exciton binding energy, are made. Thermodynamic-limit results for a combination of “optimal” and model geometries are presented. Calculated results for excitations that are adequately described using a single-particle model illustrate the benefits of requiring a KS-DFT functional to satisfy the Bartlett ionization potential theorem.
Collapse
Affiliation(s)
- Zachary W. Windom
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|