Striving to define very accurate vertical transition energies, we perform both high-level coupled cluster (CC) calculations (up to CCSDTQP) and selected configuration interaction (sCI) calculations (up to several millions of determinants) for 18 small compounds (water, hydrogen sulfide, ammonia, hydrogen chloride, dinitrogen, carbon monoxide, acetylene, ethylene, formaldehyde, methanimine, thioformaldehyde, acetaldehyde, cyclopropene, diazomethane, formamide, ketene, nitrosomethane, and the smallest streptocyanine). By systematically increasing the order of the CC expansion, the number of determinants in the CI expansion as well as the size of the one-electron basis set, we have been able to reach near full CI (FCI) quality transition energies. These calculations are carried out on CC3/ aug-cc-pVTZ geometries, using a series of increasingly large atomic basis sets systematically including diffuse functions. In this way, we define a list of 110 transition energies for states of various characters (valence, Rydberg, n → π, π → π*, singlet, triplet, etc.) to be used as references for further calculations. Benchmark transition energies are provided at the aug-cc-pVTZ level as well as with additional basis set corrections, in order to obtain results close to the complete basis set limit. These reference data are used to benchmark a series of 12 excited-state wave function methods accounting for double and triple contributions, namely ADC(2), ADC(3), CIS(D), CIS(D), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, CC3, CCSDT., and CCSDTQ. It turns out that CCSDTQ yields a negligible difference with the extrapolated CI values with a mean absolute error as small as 0.01 eV, whereas the coupled cluster approaches including iterative triples are also very accurate (mean absolute error of 0.03 eV). Consequently, CCSDT-3 and CC3 can be used to define reliable benchmarks. This observation does not hold for ADC(3) that delivers quite large errors for this set of small compounds, with a clear tendency to overcorrect its second-order version, ADC(2). Finally, we discuss the possibility to use basis set extrapolation approaches so as to tackle more easily larger compounds.

Excited states exhibiting double excitation character are notoriously di cult to model using conventional singlereference methods, such as adiabatic time-dependent density-functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC). In addition, these states are typical experimentally "dark" making their detection in photo-absorption spectra very challenging. Nonetheless, they play a key role in the faithful description of many physical, chemical, and biological processes. In the present work, we provide accurate reference excitation energies for transitions involving a substantial amount of double excitation using a series of increasingly large di use-containing atomic basis sets. Our set gathers 20 vertical transitions from 14 small-and medium-size molecules (acrolein, benzene, beryllium atom, butadiene, carbon dimer and trimer, ethylene, formaldehyde, glyoxal, hexatriene, nitrosomethane, nitroxyl, pyrazine, and tetrazine). Depending on the size of the molecule, selected con guration interaction (sCI) and/or multicon gurational (CASSCF, CASPT2, (X)MS-CASPT2 and NEVPT2) calculations are performed in order to obtain reliable estimates of the vertical transition energies. In addition, coupled cluster approaches including at least contributions from iterative triples (such as CC3, CCSDT, CCSDTQ, and CCSDTQP) are assessed. Our results clearly evidence that the error in CC methods is intimately related to the amount of double excitation character of the transition. For "pure" double excitations (i.e. for transitions which do not mix with single excitations), the error in CC3 can easily reach 1 eV, while it goes down to few tenths of an eV for more common transitions (like in trans-butadiene) involving a signi cant amount of singles. As expected, CC approaches including quadruples yield highly accurate results for any type of transitions. e quality of the excitation energies obtained with multicon gurational methods is harder to predict. We have found that the overall accuracy of these methods is highly dependent of both the system and the selected active space. e inclusion of the σ and σ orbitals in the active space, even for transitions involving mostly π and π orbitals, is mandatory in order to reach high accuracy. A theoretical best estimate (TBE) is reported for each transition. We believe that these reference data will be valuable for future methodological developments aiming at accurately describing double excitations. TOC graphical abstract arXiv:1811.12861v1 [physics.chem-ph] 30 Nov 2018

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 (π → π*, n → π*, 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 aspects 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 estimates (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 so on (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, and so on. We hope that the present review will provide a useful summary of our effort so far and foster new developments around excited‐state methods. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods

Aiming at completing the sets of FCI-quality transition energies that we recently developed (J. Chem. Theory Comput

A hybrid stochastic-deterministic approach for computing the second-order perturbative contribution E within multireference perturbation theory (MRPT) is presented. The idea at the heart of our hybrid scheme-based on a reformulation of E as a sum of elementary contributions associated with each determinant of the MR wave function-is to split E into a stochastic and a deterministic part. During the simulation, the stochastic part is gradually reduced by dynamically increasing the deterministic part until one reaches the desired accuracy. In sharp contrast with a purely stochastic Monte Carlo scheme where the error decreases indefinitely as t (where t is the computational time), the statistical error in our hybrid algorithm displays a polynomial decay ∼t with n = 3-4 in the examples considered here. If desired, the calculation can be carried on until the stochastic part entirely vanishes. In that case, the exact result is obtained with no error bar and no noticeable computational overhead compared to the fully deterministic calculation. The method is illustrated on the F and Cr molecules. Even for the largest case corresponding to the Cr molecule treated with the cc-pVQZ basis set, very accurate results are obtained for E for an active space of (28e, 176o) and a MR wave function including up to 2×10 determinants.

Aiming at completing the sets of FCI-quality transition energies that we recently developed (J. Chem. Theory Comput. 2018, 14, 4360–4379, ibid. 2019, 15, 1939–1956, and ibid. 2020, 16, 1711–1741), we provide, in the present contribution, ultra-accurate vertical excitation energies for a series of “exotic” closed-shell molecules containing F, Cl, P, and Si atoms and small radicals, such as CON and its variants, that were not considered to date in such investigations. This represents a total of 81 high-quality transitions obtained with a series of diffuse-containing basis sets of various sizes. For the exotic compounds, these transitions are used to perform benchmarks with a vast array of lower level models, i.e., CIS­(D), EOM-MP2, (SOS/SCS)-CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, (SOS-)­ADC(2), and ADC(3). Additional comparisons are made with literature data. For the open-shell compounds, we compared the performance of both the unrestricted and the restricted open-shell CCSD and CC3 formalisms.

We provide an overview of the successive steps that made possible to obtain increasingly accurate excitation energies with computational chemistry tools, eventually leading to chemically accurate vertical transition energies for small-and medium-size molecules. First, we describe the evolution of ab initio methods employed to define benchmark values, with originally Roos' CASPT2 method, then the CC3 method as in the renowned Thiel set, and more recently the resurgence of selected configuration interaction methods. The latter method has been able to deliver consistently, for both single and double excitations, highly accurate excitation energies for small molecules, as well as medium-size molecules with compact basis sets. Second, we describe how these high-level methods and the creation of representative benchmark sets of excitation energies have allowed to assess fairly and accurately the performance of computationally lighter methods. We conclude by discussing the future theoretical and technological developments in the field.

Selected configuration interaction (sCI) methods including second-order perturbative corrections provide near full CI (FCI) quality energies with only a small fraction of the determinants of the FCI space. Here, we introduce both a state-specific and a multi-state sCI method based on the configuration interaction using a perturbative selection made iteratively (CIPSI) algorithm. The present method revises the reference (internal) space under the effect of its interaction with the outer space via the construction of an effective Hamiltonian, following the shifted-Bk philosophy of Davidson and co-workers. In particular, the multi-state algorithm removes the storage bottleneck of the effective Hamiltonian via a low-rank factorization of the dressing matrix. Illustrative examples are reported for the state-specific and multi-state versions.