Search

The Bethe-Salpeter equation has been extensively employed to compute the two-body electron-hole propagator and its poles which correspond to the neutral excitation energies of the system. Through a different time-ordering, the two-body Green's function can also describe the propagation of two electrons or two holes. The corresponding poles are the double ionization potentials and double electron affinities of the system. In this work, a Bethe-Salpeter equation for the two-body particle-particle propagator is derived within the linear-response formalism using a pairing field and anomalous propagators. This framework allows us to compute kernels corresponding to different self-energy approximations ($GW$, $T$-matrix, and second-Born) as in the usual electron-hole case. The performance of these various kernels is gauged for singlet and triplet valence double ionization potentials using a set of 23 small molecules. The description of double core hole states is also analyzed., Comment: 19 pages, 8 figures

Reduced density matrix functional theory (RDMFT) and coupled cluster theory restricted to paired double excitations (pCCD) are emerging as efficient methodologies for accounting for the so-called non-dynamic electronic correlation effects. Up to now, molecular calculations have been performed with real-valued orbitals. However, before extending the applicability of these methodologies to extended systems, where Bloch states are employed, the subtleties of working with complex-valued orbitals and the consequences of imposing time-reversal symmetry must be carefully addressed. In this work, we describe the theoretical and practical implications of adopting time-reversal symmetry in RDMFT and pCCD when allowing for complex-valued orbital coefficients. The theoretical considerations primarily affect the optimization algorithms, while the practical implications raise fundamental questions about the stability of solutions. Specifically, we find that complex solutions lower the energy when non-dynamic electronic correlation effects are pronounced. We present numerical examples to illustrate and discuss these instabilities and possible problems introduced by N-representability violations.

To expand the QUEST database of highly-accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly-accurate excitation energies. Subsequently, we employ these reference values to benchmark a series of lower-order wave function approaches, including the popular ADC(2) and CC2 schemes, as well as time-dependent density-functional theory (TD-DFT), both with and without applying the Tamm-Dancoff approximation (TDA). At the TD-DFT level, we evaluate a large panel of global, range-separated, local, and double hybrid functionals. Additionally, we assess the performance of the Bethe-Salpeter equation (BSE) formalism relying on both $G_0W_0$ and ev$GW$ quasiparticle energies evaluated from various starting points. It turns out that CC2 and ADC(2.5) are the most accurate models amongst those with respective $\mathcal{O}(N^5)$ and $\mathcal{O}(N^6)$ scalings with system size. In contrast, CCSD does not outperform CC2. The best performing exchange-correlation functionals include BMK, M06-2X, M06-SX, CAM-B3LYP, $\omega$B97X-D, and LH20t, with average deviations of approximately 0.20 eV or slightly below. Errors on vertical excitation energies can be further reduced by considering double hybrids. Both SOS-$\omega$B88PP86 and SOS-$\omega$PBEPP86 exhibit particularly attractive performances with overall quality on par with CC2, whereas PBE0-DH and PBE-QIDH are only slightly less efficient. BSE/ev$GW$ calculations based on Kohn-Sham starting points have been found to be particularly effective for singlet transitions, but much less for their triplet counterparts., Comment: 27 pages, 10 figures (Supporting Information available)

Electronic resonances are metastable states that can decay by electron loss. They are ubiquitous across various fields of science, such as chemistry, physics, and biology. However, current theoretical and computational models for resonances cannot yet rival the level of accuracy achieved by bound-state methodologies. Here, we generalize selected configuration interaction (SCI) to treat resonances using the complex absorbing potential (CAP) technique. By modifying the selection procedure and the extrapolation protocol of standard SCI, the resulting CAP-SCI method yields resonance positions and widths of full configuration interaction quality. Initial results for the shape resonances of \ce{N2-} and \ce{CO-} reveal the important effect of high-order correlation, which shifts the values obtained with CAP-augmented equation-of-motion coupled-cluster with singles and doubles by more than \SI{0.1}{\eV}. The present CAP-SCI approach represents a cornerstone in the development of highly-accurate methodologies for resonances., Comment: 10 pages, 2 figures (supporting information available)

ipie is a Python-based auxiliary-field quantum Monte Carlo (AFQMC) package that has undergone substantial improvements since its initial release [J. Chem. Theory Comput., 2023, 19(1): 109-121]. This paper outlines the improved modularity and new capabilities implemented in ipie. We highlight the ease of incorporating different trial and walker types and the seamless integration of ipie with external libraries. We enable distributed Hamiltonian simulations of large systems that otherwise would not fit on single CPU node or GPU card. This development enabled us to compute the interaction energy of a benzene dimer with 84 electrons and 1512 orbitals with multi-GPUs. Using CUDA and cupy for NVIDIA GPUs, ipie supports GPU-accelerated multi-slater determinant trial wavefunctions [arXiv:2406.08314] to enable efficient and highly accurate simulations of large-scale systems. This allows for near-exact ground state energies of multi-reference clusters, [Cu$_2$O$_2$]$^{2+}$ and [Fe$_2$S$_2$(SCH$_3$)$_4$]$^{2-}$. We also describe implementations of free projection AFQMC, finite temperature AFQMC, AFQMC for electron--phonon systems, and automatic differentiation in AFQMC for calculating physical properties. These advancements position ipie as a leading platform for AFQMC research in quantum chemistry, facilitating more complex and ambitious computational method development and their applications., Comment: 18 pages, 13 figures

Hedin's equations provide an elegant route to compute the exact one-body Green's function (or propagator) via the self-consistent iteration of a set of non-linear equations. Its first-order approximation, known as $GW$, corresponds to a resummation of ring diagrams and has shown to be extremely successful in physics and chemistry. Systematic improvement is possible, although challenging, via the introduction of vertex corrections. Considering anomalous propagators and an external pairing potential, we derive a new self-consistent set of closed equations equivalent to the famous Hedin equations but having as a first-order approximation the particle-particle (pp) $T$-matrix approximation where one performs a resummation of the ladder diagrams. This pp version of Hedin's equations offers a way to go systematically beyond the $T$-matrix approximation by accounting for low-order pp vertex corrections., Comment: 12 pages, 5 figures (Supp. Mat. available)

Although selected configuration interaction (SCI) algorithms can tackle much larger Hilbert spaces than the conventional full CI (FCI) method, the scaling of their computational cost with respect to the system size remains inherently exponential. Additionally, inaccuracies in describing the correlation hole at small interelectronic distances lead to the slow convergence of the electronic energy relative to the size of the one-electron basis set. To alleviate these effects, we show that the non-Hermitian, transcorrelated (TC) version of SCI significantly compactifies the determinant space, allowing to reach a given accuracy with a much smaller number of determinants. Furthermore, we note a significant acceleration in the convergence of the TC-SCI energy as the basis set size increases. The extent of this compression and the energy convergence rate are closely linked to the accuracy of the correlation factor used for the similarity transformation of the Coulombic Hamiltonian. Our systematic investigation of small molecular systems in increasingly large basis sets illustrates the magnitude of these effects., Comment: 18 pages, 10 figures (supporting information available)

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 [http://dx.doi.org/10.1021/acs.jctc.8b01205; 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 ($\%T_1$) 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., Comment: 25 pages, 3 figures (Supporting Information available)

The cumulant expansion of the Green's function is a computationally efficient beyond-$GW$ approach renowned for its significant enhancement of satellite features in materials. In contrast to the ubiquitous $GW$ approximation of many-body perturbation theory, \textit{ab initio} cumulant expansions performed on top of $GW$ ($GW$+C) have demonstrated the capability to handle multi-particle processes by incorporating higher-order correlation effects or vertex corrections, yielding better agreements between experiment and theory for satellite structures. While widely employed in condensed matter physics, very few applications of $GW$+C have been published on molecular systems. Here, we assess the performance of this scheme on a series of 10-electron molecular systems (\ce{Ne}, \ce{HF}, \ce{H2O}, \ce{NH3}, and \ce{CH4}) where full configuration interaction estimates of the outer-valence quasiparticle and satellite energies are available., Comment: 11 pages, 3 figures (supporting information available)

Upon ionization of an atom or a molecule, another electron (or more) can be simultaneously excited. These concurrently generated states are called "satellites" (or shake-up 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 of full configuration interaction (FCI) quality computed in small molecular systems. Following the protocol developed for the quest database [https://doi.org/10.1002/wcms.1517, V\'eril et al. 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 ionization potentials are also analyzed. Their limitations in correctly modeling satellite transitions are discussed., Comment: 22 pages, 3 figures (supporting information available)

Recently, a new distributed implementation of the full configuration interaction (FCI) method has been reported [Gao et al. J. Chem Theory Comput. 2024, 20, 1185]. Thanks to a hybrid parallelization scheme, the authors were able to compute the exact energy of propane (\ce{C3H8}) in the minimal basis STO-3G. This formidable task involves handling an active space of 26 electrons in 23 orbitals or a Hilbert space of \SI{1.3d12} determinants. This is, by far, the largest FCI calculation reported to date. Here, we illustrate how, from a general point of view, selected configuration interaction (SCI) can achieve microhartree accuracy at a fraction of the computational and memory cost, via a sparse exploration of the FCI space. The present SCI calculations are performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) algorithm, as implemented in a determinant-driven way in the \textsc{quantum package} software. The present study reinforces the common wisdom that among the exponentially large number of determinants in the FCI space, only a tiny fraction of them significantly contribute to the energy. More importantly, it demonstrates the feasibility of achieving comparable accuracy using more reasonable and sustainable computational resources, hence reducing the ever-growing carbon footprint of computational chemistry., Comment: 5 pages, 2 figures (supporting information available)

To expand the existing QUEST database of accurate vertical transition energies [\href{https://doi.org/10.1002/wcms.1517}{V\'eril et al.~\textit{WIREs Comput.~Mol.~Sci.} \textbf{2021}, \textit{11}, e1517}], we have modeled more than 100 electronic excited states of different natures (local, charge-transfer, Rydberg, singlet, and triplet) in a dozen of mono- and di-substituted benzenes, including aniline, benzonitrile, chlorobenzene, fluorobenzene, nitrobenzene, among others. To establish theoretical best estimates for these vertical excitation energies, we have employed advanced coupled-cluster methods including iterative triples (CC3 and CCSDT) and, when technically possible, iterative quadruples (CC4). These high-level computational approaches provide a robust foundation for benchmarking a series of popular wave function methods. The evaluated methods all include contributions from double excitations (ADC(2), CC2, CCSD, CIS(D), EOM-MP2, STEOM-CCSD), along with schemes that also incorporate perturbative or iterative triples (ADC(3), CCSDR(3), CCSD(T)(a)$^\star$, and CCSDT-3). This systematic exploration not only broadens the scope of the QUEST database but also facilitates a rigorous assessment of different theoretical approaches in the framework of a homologous chemical series, offering valuable insights into the accuracy and reliability of these methods in such cases. We found that both ADC(2.5) and CCSDT-3 can provide consistent estimates, whereas among less expensive methods SCS-CC2 is likely the most effective approach. Importantly, we show that some lower order methods may offer reasonable trends in the homologous series while providing quite large average errors, and \emph{vice versa}. Consequently, benchmarking the accuracy of a model based solely on absolute transition energies may not be meaningful for applications involving a series of similar compounds., Comment: 14 pages, 4 figures

We reexamine $\Delta$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 where the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here is to evaluate the effectiveness of $\Delta$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 $\Delta$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 dataset, consisting of 276 excited states stemming from the \textsc{quest} database [V\'eril \textit{et al.}, \textit{WIREs Comput. Mol. Sci.} \textbf{2021}, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of $\Delta$CCSD. Except for the doubly-excited states, we found that $\Delta$CCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of \SI{0.10}{\eV} and \SI{0.07}{\eV}) is less pronounced than that obtained for singly-excited states of closed-shell systems (MAEs of \SI{0.15}{\eV} and \SI{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 $\Delta$CCSD., Comment: 16 pages, 6 figures

This work merges two different types of many-electron ensembles, namely the Theophilou-Gross-Oliveira-Kohn ensembles of ground and neutrally-excited states, and the more recent $N$-centered ensembles of neutral and charged ground states. On that basis, an in-principle exact and general, so-called extended $N$-centered, ensemble density-functional theory of charged and neutral electronic excitations is derived. We revisit in this context the concept of density-functional derivative discontinuity for neutral excitations, without ever invoking nor using the asymptotic behavior of the ensemble electronic density. The present mathematical construction fully relies on the weight dependence of the ensemble Hartree-exchange-correlation density-functional energy, which makes the theory applicable to lattice models and opens new perspectives for the description of gaps in mesoscopic systems., Comment: 21 pages, second revision

Due to the infinite summation of bubble diagrams, the $GW$ approximation of Green's function perturbation theory has proven particularly effective in the weak correlation regime, where this family of Feynman diagrams is important. However, the performance of $GW$ in multireference molecular systems, characterized by strong electron correlation, remains relatively unexplored. In the present study, we investigate the ability of $GW$ to handle closed-shell multireference systems in their singlet ground state by examining four paradigmatic scenarios. Firstly, we analyze a prototypical example of a chemical reaction involving strong correlation: the potential energy curve of \ce{BeH2} during the insertion of a beryllium atom into a hydrogen molecule. Secondly, we compute the electron detachment and attachment energies of a set of molecules that exhibit a variable degree of multireference character at their respective equilibrium geometries: \ce{LiF}, \ce{BeO}, \ce{BN}, \ce{C2}, \ce{B2}, and \ce{O3}. Thirdly, we consider a \ce{H6} cluster with a triangular arrangement, which features a notable degree of spin frustration. Finally, the dissociation curve of the \ce{HF} molecule is studied as an example of single bond breaking. These investigations highlight a nuanced perspective on the performance of $GW$ for strong correlation, depending on the level of self-consistency, the choice of initial guess, and the presence of spin-symmetry breaking at the Hartree-Fock level., Comment: 11 pages, 4 figures

Selected configuration interaction (SCI) methods have emerged as state-of-the-art methodologies for achieving high accuracy and generating benchmark reference data for ground and excited states in small molecular systems. However, their precision relies heavily on extrapolation procedures to produce a final estimate of the exact result. Using the structure of the exact electronic energy landscape, we provide a rationale for the common linear extrapolation of the variational energy as a function of the second-order perturbative correction. In particular, we demonstrate that the energy gap and the coupling between the so-called internal and external spaces are the key factors determining the rate at which the linear regime is reached. Starting from first principles, we also derive a new non-linear extrapolation formula that improves the post-processing of data generated from SCI methods and can be applied to both ground- and excited-state energies., Comment: 9 pages, 8 figures

We provide an in-depth examination of the $GW$ approximation of Green's function many-body perturbation theory by detailing both its theoretical and practical aspects in the realm of quantum chemistry. First, the quasiparticle context is introduced before delving into the derivation of Hedin's equations. From these, we explain how to derive the well-known $GW$ approximation of the self-energy. In a second time, we meticulously explain each step involved in a $GW$ calculation and what type of physical quantities can be computed. To illustrate its versatility, we consider two contrasting systems: the water molecule, a weakly correlated system, and the carbon dimer, a strongly correlated system. Each stage of the process is thoroughly detailed and explained alongside numerical results and illustrative plots. We hope that the contribution will facilitate the dissemination and democratization of Green's function-based formalisms within the computational and theoretical quantum chemistry community., Comment: 14 pages, 7 figures (supporting information available)

Molecules that violate Hund's rule and exhibit an inverted gap between the lowest singlet $S_1$ and triplet $T_1$ excited states have attracted considerable attention due to their potential applications in optoelectronics. Amongst these molecules, the triangular-shaped heptazine, and its derivatives, have been in the limelight. However, conflicting reports have arisen regarding the relative energies of $S_1$ and $T_1$. Here, we employ highly accurate levels of theory, such as CC3, to not only resolve the debate concerning the sign but also quantify the magnitude of the $S_1$-$T_1$ gap. We also determine the 0-0 energies to evaluate the significance of the vertical approximation. In addition, we compute reference $S_1$-$T_1$ gaps for a series of 10 related molecules. This enables us to benchmark lower-order methods for future applications on larger systems within the same family of compounds. This contribution can serve as a reliable foundation for the design of triangular-shaped molecules with enhanced photophysical properties., Comment: 8 pages, 3 figures, supporting information available

To enrich and enhance the diversity of the \textsc{quest} database of highly-accurate excitation energies [\href{https://doi.org/10.1002/wcms.1517}{V\'eril \textit{et al.}, \textit{WIREs Comput.~Mol.~Sci.}~\textbf{11}, e1517 (2021)}], we report vertical transition energies in transition metal compounds. Eleven diatomic molecules with singlet or doublet ground state containing a fourth-row transition metal (\ce{CuCl}, \ce{CuF}, \ce{CuH}, \ce{ScF}, \ce{ScH}, \ce{ScO}, \ce{ScS}, \ce{TiN}, \ce{ZnH}, \ce{ZnO}, and \ce{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 some cases, to provide more comprehensive benchmark data, we also provide full configuration interaction estimates computed with the \textit{"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 CC and multiconfigurational methods for this specific set of challenging transitions. Furthermore, comparisons with experimental data and previous theoretical results are also reported., Comment: 17 pages, 3 figures

We derive the explicit expression of the three self-energies that one encounters in many-body perturbation theory: the well-known $GW$ self-energy, as well as the particle-particle and electron-hole $T$-matrix self-energies. Each of these can be easily computed via the eigenvalues and eigenvectors of a different random-phase approximation (RPA) linear eigenvalue problem that completely defines their corresponding response function. For illustrative and comparative purposes, we report the principal ionization potentials of a set of small molecules computed at each level of theory. The performance of these schemes on strongly correlated systems (\ce{B2} and \ce{C2}) is also discussed., Comment: 13 pages, 6 figures

We introduce \textsc{qcmath}, a user-friendly quantum chemistry software tailored for electronic structure calculations, implemented using the Wolfram Mathematica language and available at \url{https://github.com/LCPQ/qcmath}. This software, designed with accessibility in mind, takes advantage of the symbolic capabilities intrinsic to Mathematica. Its primary goal is to provide a supportive environment for newcomers to the field of quantum chemistry, enabling them to easily conceptualize, develop, and test their own ideas. The functionalities of \textsc{qcmath} encompass a broad spectrum of methods, catering to both ground- and excited-state calculations. We provide a comprehensive overview of these capabilities, complemented by essential theoretical insights. To facilitate ease of use, we offer an exhaustive blueprint of the software's architecture. Furthermore, we provide users with comprehensive guides, addressing both the operational aspects and the more intricate programming facets of \textsc{qcmath}., Comment: 9 pages

Hierarchy configuration interaction (hCI) has been recently introduced as an alternative configuration interaction (CI) route combining excitation degree and seniority number, which showed to efficiently recover both dynamic and static correlations for closed-shell molecular systems [\href{https://doi.org/10.1021/acs.jpclett.2c00730}{\textit{J.~Phys.~Chem.~Lett.}~\textbf{2022}, \textit{13}, 4342}]. Here, we generalize hCI for an arbitrary reference determinant, allowing calculations for radicals and for 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 the 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 excitation energies of closed- and open-shell systems. Our results underline that both the choice of the reference determinant and the set of orbitals drive 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., Comment: 18 pages, 5 figures

We present an equation generator algorithm that utilizes second-quantized operators in normal order with respect to a correlated or non-correlated reference and the corresponding Wick theorem. The algorithm proposed here, written with \textsc{mathematica}, enables the generation of non-redundant strings of second-quantized operators that, after classification, are directly assigned to intermediate quantities used to construct the coupled-cluster (CC) effective Hamiltonian. We demonstrate the capabilities of the algorithm by computing the CC amplitude equations and various blocks of the equation-of-motion CC effective Hamiltonian. A comprehensive description of this four-step algorithm is provided alongside concrete examples., Comment: 15 pages, 2 figures

In recent years, Green's function methods have garnered considerable interest due to their ability to target both charged and neutral excitations. Among them, the well-established $GW$ approximation provides accurate ionization potentials and electron affinities and can be extended to neutral excitations using the Bethe-Salpeter equation (BSE) formalism. Here, we investigate the connections between various Green's function methods and evaluate their performance for charged and neutral excitations. Comparisons with other widely-known second-order wave function methods are also reported. Additionally, we calculate the singlet-triplet gap of cycl[3,3,3]azine, a model molecular emitter for thermally activated delayed fluorescence, which has the particularity of having an inverted gap thanks to a substantial contribution from the double excitations. We demonstrate that, within the $GW$ approximation, a second-order BSE kernel with dynamical correction is required to predict this distinctive characteristic., Comment: 18 pages, 1 figure

The exact functionals associated with the (singlet) ground and the two singlet excited states of the asymmetric Hubbard dimer at half-filling are calculated using both Levy's constrained search and Lieb's convex formulation. While the ground-state functional is, as commonly known, a convex function with respect to the density (or, more precisely, the site occupation), the functional associated with the (highest) doubly-excited state is found to be concave. Also, because the density of the first-excited state is non-invertible, its ``functional'' is a partial, multi-valued function composed of one concave and one convex branch that correspond to two separate sets of values of the external potential. Remarkably, it is found that, although the one-to-one mapping between density and external potential may not apply (as in the case of the first excited state), each state-specific energy and corresponding universal functional are ``functions'' whose derivatives are each other's inverse, just as in the ground state formalism. These findings offer insight into the challenges of developing state-specific excited-state density functionals for general applications in electronic structure theory., Comment: 7 pages, 6 figures

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 \href{https://doi.org/10.1063/5.0089317}{[J. Chem. Phys. 156, 231101 (2022)]} 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$ (qs$GW$) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qs$GW$ calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code., Comment: 14 pages, 7 figures (supporting information available)

The Bethe-Salpeter equation (BSE) is the key equation in many-body perturbation theory based on Green's functions to access response properties. Within the $GW$ approximation to the exchange-correlation kernel, the BSE has been successfully applied to several finite and infinite systems. However, it also shows some failures, such as underestimated triplet excitation energies, lack of double excitations, ground-state energy instabilities in the dissociation limit, etc. In this work, we study the performance of the BSE within the $GW$ approximation as well as the $T$-matrix approximation for the excitation energies of the exactly solvable asymmetric Hubbard dimer. This model allows one to study various correlation regimes by varying the on-site Coulomb interaction $U$ as well as the degree of the asymmetry of the system by varying the difference of potential $\Delta v$ between the two sites. We show that, overall, the $GW$ approximation gives more accurate excitation energies than $GT$ over a wide range of $U$ and $\Delta v$. However, the strongly-correlated (i.e., large $U$) regime still remains a challenge., Comment: 11 pages, 9 figures

We report ground- and excited-state dipole moments and oscillator strengths (computed in different ``gauges'' or representations) of full configuration interaction (FCI) quality using the selected configuration interaction method known as \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI). Thanks to a set encompassing 35 ground- and excited-state properties computed in 11 small molecules, the present near-FCI estimates allow us to assess the accuracy of high-order coupled-cluster (CC) calculations including up to quadruple excitations. In particular, we show that incrementing the excitation degree of the CC expansion (from CCSD to CCSDT or from CCSDT to CCSDTQ) reduces the average error with respect to the near-FCI reference values by approximately one order of magnitude., Comment: 12 pages, 8 figures (supporting information available)

We introduce and benchmark a systematically improvable route for excited-state calculations, state-specific configuration interaction ($\Delta$CI), \alert{which is a particular realization of multiconfigurational self-consistent field and multireference configuration interaction.} Starting with a reference built from optimized configuration state functions, separate CI calculations are performed for each targeted state (hence state-specific orbitals and determinants). Accounting for single and double excitations produces the $\Delta$CISD model, which can be improved with second-order Epstein-Nesbet perturbation theory ($\Delta$CISD+EN2) or a posteriori Davidson corrections ($\Delta$CISD+Q). These models were gauged against a vast and diverse set of 294 reference excitation energies. We have found that $\Delta$CI is significantly more accurate than standard ground-state-based CI, whereas close performances were found between $\Delta$CISD and EOM-CC2, and between $\Delta$CISD+EN2 and EOM-CCSD. For larger systems, $\Delta$CISD+Q delivers more accurate results than EOM-CC2 and EOM-CCSD. The $\Delta$CI route can handle challenging multireference problems, singly- and doubly-excited states, from closed- and open-shell species, with overall comparable accuracy, and thus represents a promising alternative to more established methodologies. In its current form, however, it is only reliable for relatively low-lying excited states., Comment: 14 pages, 3 figures (supplementary information available)

The sampling of the configuration space in diffusion Monte Carlo (DMC) is done using walkers moving randomly. In a previous work on the Hubbard model [\href{https://doi.org/10.1103/PhysRevB.60.2299}{Assaraf et al.~Phys.~Rev.~B \textbf{60}, 2299 (1999)}], it was shown that the probability for a walker to stay a certain amount of time in the same state obeys a Poisson law and that the on-state dynamics can be integrated out exactly, leading to an effective dynamics connecting only different states. Here, we extend this idea to the general case of a walker trapped within domains of arbitrary shape and size. The equations of the resulting effective stochastic dynamics are derived. The larger the average (trapping) time spent by the walker within the domains, the greater the reduction in statistical fluctuations. A numerical application to the Hubbard model is presented. Although this work presents the method for finite linear spaces, it can be generalized without fundamental difficulties to continuous configuration spaces., Comment: 14 pages, 4 figures

Here, we build on the works of Scuseria (et al.) http://dx.doi.org/10.1063/1.3043729 and Berkelbach https://doi.org/10.1063/1.5032314 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 introduced by Nooijen are also put forward. The present work allows to easily transfer key developments and 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., Comment: 6 pages

Based on 280 reference vertical transition energies of various natures (singlet, triplet, valence, Rydberg, $n\to\pi^*$, $\pi\to\pi^*$, and double excitations) extracted from the QUEST database, we assess the accuracy of third-order multireference perturbation theory, CASPT3, in the context of molecular excited states. When one applies the disputable ionization-potential-electron-affinity (IPEA) shift, we show that CASPT3 provides a similar accuracy as its second-order counterpart, CASPT2, with the same mean absolute error of $0.11$ eV. However, as already reported, we also observe that the accuracy of CASPT3 is almost insensitive to the IPEA shift, irrespective of the transition type and system size, with a small reduction of the mean absolute error to $0.09$ eV when the IPEA shift is switched off., Comment: 12 pages, 3 figures (supplementary material available)

Cyclobutadiene is a well-known playground for theoretical chemists and is particularly suitable to test ground- and excited-state methods. Indeed, due to its high spatial symmetry, especially at the $D_{4h}$ square geometry but also in the $D_{2h}$ rectangular arrangement, the ground and excited states of cyclobutadiene exhibit multi-configurational characters and single-reference methods, such as adiabatic time-dependent density-functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC), are notoriously known to struggle in such situations. In this work, using a large panel of methods and basis sets, we provide an extensive computational study of the automerization barrier (defined as the difference between the square and rectangular ground-state energies) and the vertical excitation energies at $D_{2h}$ and $D_{4h}$ equilibrium structures. In particular, selected configuration interaction (SCI), multi-reference perturbation theory (CASSCF, CASPT2, and NEVPT2), and coupled-cluster (CCSD, CC3, CCSDT, CC4, and CCSDTQ) calculations are performed. The spin-flip formalism, which is known to provide a qualitatively correct description of these diradical states, is also tested within TD-DFT (combined with numerous exchange-correlation functionals) and the algebraic diagrammatic construction [ADC(2)-s, ADC(2)-x, and ADC(3)] schemes. A theoretical best estimate is defined for the automerization barrier and for each vertical transition energy., Comment: 15 pages, 4 figures (supporting information available)

We propose a novel partitioning of the Hilbert space, hierarchy configuration interaction (hCI), where the excitation degree (with respect to a given reference determinant) and the seniority number (i.e., the number of unpaired electrons) are combined in a single hierarchy parameter. The key appealing feature of hCI is that each hierarchy level accounts for all classes of determinants whose number share the same scaling with system size. By surveying the dissociation of multiple molecular systems, we found that the overall performance of hCI usually exceeds or, at least, parallels that of excitation-based CI. For higher orders of hCI and excitation-based CI, the additional computational burden related to orbital optimization usually do not compensate the marginal improvements compared with results obtained with Hartree-Fock orbitals. The exception is orbital-optimized CI with single excitations, a minimally correlated model displaying the qualitatively correct description of single bond breaking, at a very modest computational cost., Comment: 7 pages, 3 figures (supporting information available)

The uniform electron gas (UEG), a hypothetical system with finite homogenous electron density composed by an infinite number of electrons in a box of infinite volume, is the practical pillar of density-functional theory (DFT) and the foundation of the most acclaimed approximation of DFT, the local-density approximation (LDA). In the last thirty years, the knowledge of analytical parametrizations of the infinite UEG (IUEG) exchange-correlation energy has allowed researchers to perform a countless number of approximate electronic structure calculations for atoms, molecules, and solids. Recently, it has been shown that the traditional concept of the IUEG is not the unique example of UEGs, and systems, in their lowest-energy state, consisting of electrons that are confined to the surface of a sphere provide a new family of UEGs with more customizable properties. Here, we show that, some of the excited states associated with these systems can be classified as transient UEGs (TUEGs) as their electron density is only homogenous for very specific values of the radius of the sphere even though the electronic wave function is not rotationally invariant. Concrete examples are provided in the case of two-electron systems., Comment: 5 pages, 2 figures, submitted for Peter Gill Festschrift to appear in Mol. Phys

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\pm1)$-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 convergence of partially self-consistent $GW$ calculations., Comment: 8 pages, 5 figures (supplementary material available)

While the well-established $GW$ approximation corresponds to a resummation of the direct ring diagrams and is particularly well suited for weakly-correlated systems, the $T$-matrix approximation does sum ladder diagrams up to infinity and is supposedly more appropriate in the presence of strong correlation. Here, we derive and implement, for the first time, the static and dynamic Bethe-Salpeter equations when one considers $T$-matrix quasiparticle energies as well as a $T$-matrix-based kernel. The performance of the static scheme and its perturbative dynamical correction are assessed by computing the neutral excited states of molecular systems. Comparison with more conventional schemes as well as other wave function methods are also reported. Our results suggest that the $T$-matrix-based formalism performs best in few-electron systems where the electron density remains low., Comment: 10 pages, 6 figures

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events, but also to provide reference values for vertical transition energies, hence allowing to benchmark lower-order models. In this category, both CASPT2 and NEVPT2 are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches the present work relies on highly-accurate ($\pm 0.03$ eV) \emph{aug}-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 $n \to \pi^*$, 119 $\pi \to \pi^*$, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially-contracted (PC) and strongly-contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of $0.11$ and $0.13$ eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focussed on single-reference wave function methods (\textit{J.~Chem. Theory Comput.} \textbf{14}, 4360 (2018); \emph{ibid.}, \textbf{16}, 1711 (2020)), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character., Comment: 21 pages, 3 figure (supporting information available)

Pursuing our efforts to define highly-accurate estimates of the relative energies of excited states in organic molecules, we investigate, with coupled-cluster methods including iterative triples (CC3 and CCSDT), the vertical excitation energies of 10 bicyclic molecules (azulene, benzoxadiazole, benzothiadiazole, diketopyrrolopyrrole, fuofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide \emph{aug}-cc-pVTZ reference vertical excitation energies for 91 excited states of these relatively large systems. We use these reference values to benchmark various wave function methods, i.e., CIS(D), EOM-MP2, CC2, CCSD, STEOM-CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, ADC(2), ADC(2.5), ADC(3), as well as some spin-scaled variants of both CC2 and ADC(2). These results are compared to those obtained previously on smaller molecules. It turns out that while the accuracy of some methods is almost unaffected by system size, e.g., CIS(D) and CC3, the performance of others can significantly deteriorate as the systems grow, e.g., EOM-MP2 and CCSD, whereas others, e.g., ADC(2) and CC2, become more accurate for larger derivatives., Comment: 19 pages, 2 figures

Following our recent work on the benzene molecule [\href{https://doi.org/10.1063/5.0027617}{J.~Chem.~Phys.~\textbf{153}, 176101 (2020)}], itself motivated by the blind challenge of Eriksen \textit{et al.} [\href{https://doi.org/10.1021/acs.jpclett.0c02621}{J.~Phys.~Chem.~Lett.~\textbf{11}, 8922 (2020)}] on the same system, we report accurate full configuration interaction (FCI) frozen-core correlation energy estimates for twelve five- and six-membered ring molecules in the standard correlation-consistent double-$\zeta$ Dunning basis set (cc-pVDZ). Our FCI correlation energy estimates, with estimated error smaller than 1 millihartree, are based on energetically optimized-orbital selected configuration interaction (SCI) calculations performed with the \textit{Configuration Interaction using a Perturbative Selection made Iteratively} (CIPSI) algorithm. Having at our disposal these accurate reference energies, the respective performance and convergence properties of several popular and widely-used families of single-reference quantum chemistry methods are investigated. In particular, we study the convergence properties of i) the M{\o}ller-Plesset perturbation series up to fifth-order (MP2, MP3, MP4, and MP5), ii) the iterative approximate coupled-cluster series CC2, CC3, and CC4, and iii) the coupled-cluster series CCSD, CCSDT, and CCSDTQ. The performance of the ground-state gold standard CCSD(T) as well as the completely renormalized CC model, CR-CC(2,3), are also investigated., Comment: 13 pages, 5 figures

Using the simple (symmetric) Hubbard dimer, we analyze some important features of the $GW$ approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot $GW$ method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard $GW$ approximation and find, in particular, that i) some neutral excitation energies become complex when the electron-electron interaction $U$ increases, which can be traced back to the approximate nature of the $GW$ quasiparticle energies; ii) the BSE formalism yields accurate correlation energies over a wide range of $U$ when the trace (or plasmon) formula is employed; iii) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); iv) the trace formula has the correct behavior for weak (\ie, small $U$) interaction, unlike the ACFDT expression., Comment: 12 pages, 4 figures

In single-reference coupled-cluster (CC) methods, one has to solve a set of non-linear polynomial equations in order to determine the so-called amplitudes which are then used to compute the energy and other properties. Although it is of common practice to converge to the (lowest-energy) ground-state solution, it is also possible, thanks to tailored algorithms, to access higher-energy roots of these equations which may or may not correspond to genuine excited states. Here, we explore the structure of the energy landscape of variational CC (VCC) and we compare it with its (projected) traditional version (TCC) in the case where the excitation operator is restricted to paired double excitations (pCCD). By investigating two model systems (the symmetric stretching of the linear \ce{H4} molecule and the continuous deformation of the square \ce{H4} molecule into a rectangular arrangement) in the presence of weak and strong correlations, the performance of VpCCD and TpCCD are gauged against their configuration interaction (CI) equivalent, known as doubly-occupied CI (DOCI), for reference Slater determinants made of ground- or excited-state Hartree-Fock orbitals or state-specific orbitals optimized directly at the VpCCD level. The influence of spatial symmetry breaking is also investigated., Comment: 16 pages, 8 figures

Fractional-spin errors are inherent in all current approximate density functionals, including Hartree-Fock theory, and their origin has been related to strong static correlation effects. The conventional way to encode fractional-spin calculations is to construct an ensemble density that scales between the high-spin and low-spin densities. In this article, we explore the variation of the Hartree-Fock fractional-spin (or ghost-interaction) error in one-electron systems using restricted and unrestricted ensemble densities, and the exact generalized Hartree-Fock representation. By considering the hydrogen atom and H$_2^+$ cation, we analyze how the unrestricted and generalized Hartree-Fock schemes minimize this error by localizing the electrons or rotating the spin coordinates. We also reveal a clear similarity between the Coulomb hole of He-like ions and the density depletion near the nucleus induced by the fractional-spin error in the unpolarized hydrogen atom. Finally, we analyze the effect of the fractional-spin error on the M{\o}ller-Plesset adiabatic connection, excited states, and functional- and density-driven errors., Comment: 12 pages, 9 figures

We report the first investigation of the performance of EOM-CC4 -- an approximate equation-of-motion coupled-cluster model which includes iterative quadruple excitations -- for vertical excitation energies in molecular systems. By considering a set of 28 excited states in 10 small molecules for which we have computed CCSDTQP and FCI reference energies, we show that, in the case of excited states with a dominant contribution from the single excitations, CC4 yields excitation energies with sub-kJ~mol$^{-1}$ accuracy (i.e., error below $0.01$ eV), in very close agreement with its more expensive CCSDTQ parent. Therefore, if one aims at high accuracy, CC4 stands as a highly competitive approximate method to model molecular excited states, with a significant improvement over both CC3 and CCSDT. Our results also evidence that, although the same qualitative conclusions hold, one cannot reach the same level of accuracy for transitions with a dominant contribution from the double excitations., Comment: 6 pages, 2 figures

The pair coupled cluster doubles (pCCD) method (where the excitation manifold is restricted to electron pairs) has a series of interesting features. Among others, it provides ground-state energies very close to what is obtained with doubly-occupied configuration interaction (DOCI), but with polynomial cost (compared with the exponential cost of the latter). Here, we address whether this similarity holds for excited states, by exploring the symmetric dissociation of the linear \ce{H4} molecule. When ground-state Hartree-Fock (HF) orbitals are employed, pCCD and DOCI excited-state energies do not match, a feature that is assigned to the poor HF reference. In contrast, by optimizing the orbitals at the pCCD level (oo-pCCD) specifically for each excited state, the discrepancies between pCCD and DOCI decrease by one or two orders of magnitude. Therefore, the pCCD and DOCI methodologies still provide comparable energies for excited states, but only if suitable, state-specific orbitals are adopted. We also assessed whether a pCCD approach could be used to directly target doubly-excited states, without having to resort to the equation-of-motion (EOM) formalism. In our $\Delta$oo-pCCD model, excitation energies were extracted from the energy difference between separate oo-pCCD calculations for the ground state and the targeted excited state. For a set comprising the doubly-excited states of \ce{CH+}, \ce{BH}, nitroxyl, nitrosomethane, and formaldehyde, we found that $\Delta$oo-pCCD provides quite accurate excitation energies, with root mean square deviations (with respect to full configuration interaction results) lower than CC3 and comparable to EOM-CCSDT, two methods with much higher computational cost., Comment: 12 pages, 4 figures

In the aim of completing our previous efforts devoted to local and Rydberg transitions in organic compounds, we provide a series of highly-accurate vertical transition energies for intramolecular charge-transfer transitions occurring in ($\pi$-conjugated) molecular compounds. To this end we apply a composite protocol consisting of linear-response CCSDT excitation energies determined with Dunning's double-$\zeta$ basis set corrected by CC3/CCSDT-3 energies obtained with the corresponding triple-$\zeta$ basis. Further basis set corrections (up to \emph{aug}-cc-pVQZ) are obtained at the CCSD and CC2 level. We report 30 transitions obtained in 17 compounds. These reference values are then used to benchmark a series of wave function (CIS(D), SOPPA, RPA(D), EOM-MP2, CC2, CCSD, CCSD(T)(a)*, CCSDR(3), CCSDT-3, CC3, ADC(2), ADC(3), and ADC(2.5)), the Green's function-based Bethe-Salpeter equation (BSE) formalism performed on top of the partially self-consistent ev$GW$ scheme considering two different starting points (BSE/ev$GW$@HF and BSE/ev$GW$@PBE0), and TD-DFT combined with several exchange-correlation functionals (B3LYP, PBE0, M06-2X, CAM-B3LYP, LC-$\omega$HPBE, $\omega$B97X, $\omega$B97X-D, and M11)., Comment: 25 pages, 3 figures, 3 tables (SI available)

Like adiabatic time-dependent density-functional theory (TD-DFT), the Bethe-Salpeter equation (BSE) formalism of many-body perturbation theory, in its static approximation, is "blind" to double (and higher) excitations, which are ubiquitous, for example, in conjugated molecules like polyenes. Here, we apply the spin-flip \textit{ansatz} (which considers the lowest triplet state as the reference configuration instead of the singlet ground state) to the BSE formalism in order to access, in particular, double excitations. The present scheme is based on a spin-unrestricted version of the $GW$ approximation employed to compute the charged excitations and screened Coulomb potential required for the BSE calculations. Dynamical corrections to the static BSE optical excitations are taken into account via an unrestricted generalization of our recently developed (renormalized) perturbative treatment. The performance of the present spin-flip BSE formalism is illustrated by computing excited-state energies of the beryllium atom, the hydrogen molecule at various bond lengths, and cyclobutadiene in its rectangular and square-planar geometries., Comment: 14 pages, 3 figures and 3 tables

We explore the non-Hermitian extension of quantum chemistry in the complex plane and its link with perturbation theory. We observe that the physics of a quantum system is intimately connected to the position of complex-valued energy singularities, known as exceptional points. After presenting the fundamental concepts of non-Hermitian quantum chemistry in the complex plane, including the mean-field Hartree--Fock approximation and Rayleigh--Schr\"odinger perturbation theory, we provide a historical overview of the various research activities that have been performed on the physics of singularities. In particular, we highlight seminal work on the convergence behaviour of perturbative series obtained within M{\o}ller--Plesset perturbation theory, and its links with quantum phase transitions. We also discuss several resummation techniques (such as Pad\'e and quadratic approximants) that can improve the overall accuracy of the M{\o}ller--Plesset perturbative series in both convergent and divergent cases. Each of these points is illustrated using the Hubbard dimer at half filling, which proves to be a versatile model for understanding the subtlety of analytically-continued perturbation theory in the complex plane., Comment: 22 page, 14 figures, 4 tables

We describe our efforts of the past few years to create a large set of more than 500 highly-accurate vertical excitation energies of various natures ($\pi \to \pi^*$, $n \to \pi^*$, double excitation, Rydberg, singlet, doublet, triplet, etc) in small- and medium-sized molecules. These values have been obtained using an incremental strategy which consists in combining high-order coupled cluster and selected configuration interaction calculations using increasingly large diffuse basis sets in order to reach high accuracy. One of the key aspect of the so-called QUEST database of vertical excitations is that it does not rely on any experimental values, avoiding potential biases inherently linked to experiments and facilitating theoretical cross comparisons. Following this composite protocol, we have been able to produce theoretical best estimate (TBEs) with the aug-cc-pVTZ basis set for each of these transitions, as well as basis set corrected TBEs (i.e., near the complete basis set limit) for some of them. The TBEs/aug-cc-pVTZ have been employed to benchmark a large number of (lower-order) wave function methods such as CIS(D), ADC(2), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, ADC(3), CC3, NEVPT2, and others (including spin-scaled variants). In order to gather the huge amount of data produced during the QUEST project, we have created a website [https://lcpq.github.io/QUESTDB_website] where one can easily test and compare the accuracy of a given method with respect to various variables such as the molecule size or its family, the nature of the excited states, the type of basis set, etc. We hope that the present review will provide a useful summary of our effort so far and foster new developments around excited-state methods., Comment: 44 pages, 5 figures, 4 Tables, supplementary information available at https://doi.org/10.5281/zenodo.4297012

Using a set of oscillator strengths and excited-state dipole moments of near full configuration interaction (FCI) quality determined for small compounds, we benchmark the performances of several single-reference wave function methods (CC2, CCSD, CC3, CCSDT, ADC(2), and ADC(3/2)) and time-dependent density-functional theory (TD-DFT) with various functionals (B3LYP, PBE0, M06-2X, CAM-B3LYP, and $\omega$B97X-D). We consider the impact of various gauges (length, velocity, and mixed) and formalisms: equation of motion (EOM) \emph{vs} linear response (LR), relaxed \emph{vs} unrelaxed orbitals, etc. Beyond the expected accuracy improvements and a neat decrease of formalism sensitivy when using higher-order wave function methods, the present contribution shows that, for both ADC(2) and CC2, the choice of gauge impacts more significantly the magnitude of the oscillator strengths than the choice of formalism, and that CCSD yields a notable improvement on this transition property as compared to CC2. For the excited-state dipole moments, switching on orbital relaxation appreciably improves the accuracy of both ADC(2) and CC2, but has a rather small effect at the CCSD level. Going from ground to excited states, the typical errors on dipole moments for a given method tend to roughly triple. Interestingly, the ADC(3/2) oscillator strengths and dipoles are significantly more accurate than their ADC(2) counterparts, whereas the two models do deliver rather similar absolute errors for transition energies. Concerning TD-DFT, one finds: i) a rather negligible impact of the gauge on oscillator strengths for all tested functionals (except for M06-2X); ii) deviations of ca.~0.10 D on ground-state dipoles for all functionals; iii) the better overall performance of CAM-B3LYP for the two considered excited-state properties., Comment: 18 pages, 7 figures