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
|
Montero-Cabrera LA, Montero-Alejo AL, Aspuru-Guzik A, García de la Vega JM, Piris M, Díaz-Fernández LA, Pérez-Badell Y, Guerra-Barroso A, Alfonso-Ramos JE, Rodríguez J, Fuentes ME, de Armas CM. Alternative CNDOL Fockians for fast and accurate description of molecular exciton properties. J Chem Phys 2024; 160:214108. [PMID: 38828812 DOI: 10.1063/5.0208809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2024] [Accepted: 05/14/2024] [Indexed: 06/05/2024] Open
Abstract
CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb-exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.
Collapse
Affiliation(s)
- Luis A Montero-Cabrera
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
- Donostia International Physics Center (DIPC), 20018 Donostia - San Sebastián, Basque Country, Spain
| | - Ana L Montero-Alejo
- Departamento de Física, Facultad de Ciencias Naturales, Matemática y del Medio Ambiente (FCNMM), Universidad Tecnológica Metropolitana; Ñuñoa, Santiago 7800002, Chile
| | - Alan Aspuru-Guzik
- Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | | | - Mario Piris
- Donostia International Physics Center (DIPC), 20018 Donostia - San Sebastián, Basque Country, Spain
| | - Lourdes A Díaz-Fernández
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Yoana Pérez-Badell
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Alberto Guerra-Barroso
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Javier E Alfonso-Ramos
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - Javier Rodríguez
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| | - María E Fuentes
- Facultad de Ciencias Químicas, Universidad Autónoma de Chihuahua, Chihuahua, 31100 Chihuahua, Mexico
| | - Carlos M de Armas
- Laboratorio de Química Computacional y Teórica, Departamento de Química Física, Universidad de La Habana, 10400 Havana, Cuba
| |
Collapse
|
3
|
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
|
4
|
Selenius E, Sigurdarson AE, Schmerwitz YLA, Levi G. Orbital-Optimized Versus Time-Dependent Density Functional Calculations of Intramolecular Charge Transfer Excited States. J Chem Theory Comput 2024; 20:3809-3822. [PMID: 38695313 DOI: 10.1021/acs.jctc.3c01319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2024]
Abstract
The performance of time-independent, orbital-optimized calculations of excited states is assessed with respect to charge transfer excitations in organic molecules in comparison to the linear-response time-dependent density functional theory (TD-DFT) approach. A direct optimization method to converge on saddle points of the electronic energy surface is used to carry out calculations with the local density approximation (LDA) and the generalized gradient approximation (GGA) functionals PBE and BLYP for a set of 27 excitations in 15 molecules. The time-independent approach is fully variational and provides a relaxed excited state electron density from which the extent of charge transfer is quantified. The TD-DFT calculations are generally found to provide larger charge transfer distances compared to the orbital-optimized calculations, even when including orbital relaxation effects with the Z-vector method. While the error on the excitation energy relative to theoretical best estimates is found to increase with the extent of charge transfer up to ca. -2 eV for TD-DFT, no correlation is observed for the orbital-optimized approach. The orbital-optimized calculations with the LDA and the GGA functionals provide a mean absolute error of ∼0.7 eV, outperforming TD-DFT with both local and global hybrid functionals for excitations with a long-range charge transfer character. Orbital-optimized calculations with the global hybrid functional B3LYP and the range-separated hybrid functional CAM-B3LYP on a selection of states with short- and long-range charge transfer indicate that inclusion of exact exchange has a small effect on the charge transfer distance, while it significantly improves the excitation energy, with the best-performing functional CAM-B3LYP providing an absolute error typically around 0.15 eV.
Collapse
Affiliation(s)
- Elli Selenius
- Science Institute of the University of Iceland, Reykjavík 107, Iceland
| | | | | | - Gianluca Levi
- Science Institute of the University of Iceland, Reykjavík 107, Iceland
| |
Collapse
|
5
|
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
|
6
|
Silva RC, Canisares FSM, Saraiva LF, Pires AM, Lima SAM. Featuring long-lifetime deep-red emitting iridium III complexes with high colour purity: insights into the excited state dynamics from spectroscopic and theoretical perspectives. Dalton Trans 2024; 53:5466-5477. [PMID: 38414337 DOI: 10.1039/d3dt02531d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
The significant attention drawn to IrIII-complexes in recent years has boosted the development of new compounds with advantageous photophysical features. However, obtaining IrIII deep-red-emitting complexes with long lived excited state, high colour purity and high quantum yield (Φ) remains a challenging task. To address this issue, this study reports the synthesis and photophysical investigation of three novel zwitterionic complexes, [Ir(C^N)2bqdc] (C^N = ppy, phq, and bzq), with ppy = 2-phenylpyridine (Ir-pb), phq = 2-phenylquinoline (Ir-qb), bzq = benzo[h]quinoline (Ir-bb), and bqdc = potassium 2,2'-biquinoline-4,4'-dicarboxylate. These complexes exhibit high quantum yields and long emission lifetimes with high colour purity in the deep-red region. The structural characterisation carried out by usual spectroscopic measurements supports the proposed structures, while the photophysical study unveiled the high contribution of the 3MLCT state to the hybrid emitter state, as endorsed by theoretical investigations. The desired correspondence between the calculations and the experimental data set affirmed the accuracy of the theoretical analysis, which enabled us to establish a relationship between the ground-to-excited-state geometry distortion and the lifetime through the nonradiative decay (knr). Furthermore, these newly synthesized complexes exhibit quenching in the presence of molecular oxygen. In deoxygenated DMSO solution the knr values halve, increasing the quantum yields (34.0, 10.6, and 26.6%) and the lifetimes (1.13, 1.11, and 1.72 μs), while leading to quite pure deep-red emission - CIE coordinates: (0.67, 0.33), (0.60, 0.40;), (0.65, 0.35) for Ir-pb, Ir-qb, and Ir-bb, respectively. These complexes display considerable potential for a wide range of applications, such as photodynamic therapy, due to their attractive photophysical properties, and are among the deep-red-emitting complexes reported in the literature with longer lifetimes and higher Φ.
Collapse
Affiliation(s)
- Renan C Silva
- São Paulo State University (Unesp), School of Technology and Sciences, Presidente Prudente, SP, Brazil.
- São Paulo State University (Unesp), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, SP, Brazil
| | - Felipe S M Canisares
- São Paulo State University (Unesp), School of Technology and Sciences, Presidente Prudente, SP, Brazil.
- São Paulo State University (Unesp), Institute of Chemistry, Araraquara, SP, Brazil
| | - Leonardo F Saraiva
- São Paulo State University (Unesp), School of Technology and Sciences, Presidente Prudente, SP, Brazil.
- São Paulo State University (Unesp), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, SP, Brazil
| | - Ana M Pires
- São Paulo State University (Unesp), School of Technology and Sciences, Presidente Prudente, SP, Brazil.
- São Paulo State University (Unesp), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, SP, Brazil
- São Paulo State University (Unesp), Institute of Chemistry, Araraquara, SP, Brazil
| | - Sergio A M Lima
- São Paulo State University (Unesp), School of Technology and Sciences, Presidente Prudente, SP, Brazil.
- São Paulo State University (Unesp), Institute of Biosciences, Humanities and Exact Sciences, São José do Rio Preto, SP, Brazil
- São Paulo State University (Unesp), Institute of Chemistry, Araraquara, SP, Brazil
| |
Collapse
|
7
|
Lew-Yee JFH, Bonfil-Rivera IA, Piris M, M. del Campo J. Excited States by Coupling Piris Natural Orbital Functionals with the Extended Random-Phase Approximation. J Chem Theory Comput 2024; 20:2140-2151. [PMID: 38353418 PMCID: PMC10938499 DOI: 10.1021/acs.jctc.3c01194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 01/30/2024] [Accepted: 01/31/2024] [Indexed: 03/13/2024]
Abstract
In this work, we explore the use of Piris natural orbital functionals (PNOFs) to calculate excited-state energies by coupling their reconstructed second-order reduced density matrix with the extended random-phase approximation (ERPA). We have named the general method PNOF-ERPA, and specific approaches are referred to as PNOF-ERPA0, PNOF-ERPA1, and PNOF-ERPA2, according to the way the excitation operator is built. The implementation has been tested in the first excited states of H2, HeH+, LiH, Li2, and N2 showing good results compared to the configuration interaction (CI) method. As expected, an increase in accuracy is observed on going from ERPA0 to ERPA1 and ERPA2. We also studied the effect of electron correlation included by PNOF5, PNOF7, and the recently proposed global NOF (GNOF) on the predicted excited states. PNOF5 appears to be good and may even provide better results in very small systems, but including more electron correlation becomes important as the system size increases, where GNOF achieves better results. Overall, the extension of PNOF to excited states has been successful, making it a promising method for further applications.
Collapse
Affiliation(s)
- Juan Felipe Huan Lew-Yee
- Departamento
de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México City C.P.
04510, Mexico
- Donostia
International Physics Center (DIPC), 20018 Donostia, Spain
| | - Iván Alejandro Bonfil-Rivera
- Departamento
de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México City C.P.
04510, Mexico
| | - Mario Piris
- Donostia
International Physics Center (DIPC), 20018 Donostia, Spain
- Kimika
Fakultatea, Euskal Herriko Unibertsitatea
(UPV/EHU), 20018 Donostia, Spain
- IKERBASQUE,
Basque Foundation for Science, 48013 Bilbao, Spain
| | - Jorge M. del Campo
- Departamento
de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México City C.P.
04510, Mexico
| |
Collapse
|
8
|
Arias-Martinez JE, Wu H, Head-Gordon M. Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone. J Chem Theory Comput 2024; 20:752-766. [PMID: 38164934 DOI: 10.1021/acs.jctc.3c01139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.
Collapse
Affiliation(s)
- Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Hamlin Wu
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Oña OB, Massaccesi GE, Melo JI, Torre A, Lain L, Alcoba DR, Peralta JE. Generalized spin σ-SCF method. J Chem Phys 2023; 159:214104. [PMID: 38047509 DOI: 10.1063/5.0178264] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 11/07/2023] [Indexed: 12/05/2023] Open
Abstract
We introduce a generalization of the σ-SCF method to approximate noncollinear spin ground and excited single-reference electronic states by minimizing the Hamiltonian variance. The new method is based on the σ-SCF method, originally proposed by Ye et al. [J. Chem. Phys. 147, 214104 (2017)], and provides a prescription to determine ground and excited noncollinear spin states on an equal footing. Our implementation was carried out utilizing an initial simulated annealing stage followed by a mean-field iterative self-consistent approach to simplify the cumbersome search introduced by generalizing the spin degrees of freedom. The simulated annealing stage ensures a broad exploration of the Hilbert space spanned by the generalized spin single-reference states with random complex element-wise rotations of the generalized density matrix elements in the simulated annealing stage. The mean-field iterative self-consistent stage employs an effective Fockian derived from the variance, which is utilized to converge tightly to the solutions. This process helps us to easily find complex spin structures, avoiding manipulating the initial guess. As proof-of-concept tests, we present results for Hn (n = 3-7) planar rings and polyhedral clusters with geometrical spin frustration. We show that most of these systems have noncollinear spin excited states that can be interpreted in terms of geometric spin frustration. These states are not directly targeted by energy minimization methods, which are meant to converge to the ground state. This stresses the capability of the σ-SCF methodology to find approximate noncollinear spin structures as mean-field excited states.
Collapse
Affiliation(s)
- Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- Instituto de Investigaciones Matemáticas "Luis A. Santaló" (IMAS), Consejo Nacional de Investigaciones Científicas y Técnicas, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Juan I Melo
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- CONICET - Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Diego R Alcoba
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- CONICET - Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Juan E Peralta
- Department of Physics and Science of Advanced Materials, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| |
Collapse
|
10
|
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
|
11
|
Kossoski F, Loos PF. Seniority and Hierarchy Configuration Interaction for Radicals and Excited States. J Chem Theory Comput 2023. [PMID: 37965728 DOI: 10.1021/acs.jctc.3c00946] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023]
Abstract
Hierarchy configuration interaction (hCI) has recently been introduced as an alternative configuration interaction (CI) route combining excitation degree and seniority number and has been shown to efficiently recover both dynamic and static correlations for closed-shell molecular systems [ J. Phys. Chem. Lett. 2022, 13, 4342]. Here we generalize hCI for an arbitrary reference determinant, allowing calculations for radicals and excited states in a state-specific way. We gauge this route against excitation-based CI (eCI) and seniority-based CI (sCI) by evaluating how different ground-state properties of radicals converge to the full CI limit. We find that hCI outperforms or matches eCI, whereas sCI is far less accurate, in line with previous observations for closed-shell molecules. Employing second-order Epstein-Nesbet (EN2) perturbation theory as a correction significantly accelerates the convergence of hCI and eCI. We further explore various hCI and sCI models to calculate the excitation energies of closed- and open-shell systems. Our results underline that the choice of both the reference determinant and the set of orbitals drives the fine balance between correlation of ground and excited states. State-specific hCI2 and higher-order models perform similarly to their eCI counterparts, whereas lower orders of hCI deliver poor results unless supplemented by the EN2 correction, which substantially improves their accuracy. In turn, sCI1 produces decent excitation energies for radicals, encouraging the development of related seniority-based coupled-cluster methods.
Collapse
Affiliation(s)
- Fábris Kossoski
- 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
|
12
|
Tuckman H, Neuscamman E. Excited-State-Specific Pseudoprojected Coupled-Cluster Theory. J Chem Theory Comput 2023; 19:6160-6171. [PMID: 37676752 DOI: 10.1021/acs.jctc.3c00194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
We present an excited-state-specific coupled-cluster approach in which both the molecular orbitals and cluster amplitudes are optimized for an individual excited state. The theory is formulated via a pseudoprojection of the traditional coupled-cluster wavefunction that allows correlation effects to be introduced atop an excited-state mean field starting point. The approach shares much in common with ground-state CCSD, including size extensivity and an N6 cost scaling. Preliminary numerical tests show that, when augmented with N5 cost perturbative corrections for key terms, the method can improve over excited-state-specific second-order perturbation theory in valence, charge transfer, and Rydberg states.
Collapse
Affiliation(s)
- Harrison Tuckman
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
13
|
Marie A, Burton HGA. Excited States, Symmetry Breaking, and Unphysical Solutions in State-Specific CASSCF Theory. J Phys Chem A 2023; 127:4538-4552. [PMID: 37141564 DOI: 10.1021/acs.jpca.3c00603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/06/2023]
Abstract
State-specific electronic structure theory provides a route toward balanced excited-state wave functions by exploiting higher-energy stationary points of the electronic energy. Multiconfigurational wave function approximations can describe both closed- and open-shell excited states and avoid the issues associated with state-averaged approaches. We investigate the existence of higher-energy solutions in complete active space self-consistent field (CASSCF) theory and characterize their topological properties. We demonstrate that state-specific approximations can provide accurate higher-energy excited states in H2 (6-31G) with more compact active spaces than would be required in a state-averaged formalism. We then elucidate the unphysical stationary points, demonstrating that they arise from redundant orbitals when the active space is too large or symmetry breaking when the active space is too small. Furthermore, we investigate the singlet-triplet crossing in CH2 (6-31G) and the avoided crossing in LiF (6-31G), revealing the severity of root flipping and demonstrating that state-specific solutions can behave quasi-diabatically or adiabatically. These results elucidate the complexity of the CASSCF energy landscape, highlighting the advantages and challenges of practical state-specific calculations.
Collapse
Affiliation(s)
- Antoine Marie
- Physical and Theoretical Chemical Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, U.K
| | - Hugh G A Burton
- Physical and Theoretical Chemical Laboratory, Department of Chemistry, University of Oxford, Oxford OX1 3QZ, U.K
| |
Collapse
|