Theoretical 0-0 Energies with Chemical Accuracy

Absorption Intensities of Organic Molecules from Electronic Structure Calculations versus Experiments: the Effect of Solvation, Method, Basis Set, and Transition Moment Gauge

Recently, we derived experimental oscillator strengths (OSs) from well-defined UV-visible absorption spectral peaks of 100 molecules in solution. Here, we focus on a subset of transitions with the highest reliability to further benchmark the OSs from several wave function methods and density functionals. We consider multiple basis sets, transition moment gauges (length, velocity, and mixed), and solvent corrections. Most transitions in the comparison set come from conjugated molecules and have π → π* character. We use an automated algorithm to assign computed transitions to experimental bands. OSs computed using the Tamm-Dancoff approximation (TDA), CIS, or EOM-CCSD exhibited a strong gauge dependence, which is diminished in linear response theories (TD-DFT, TD-HF, and to a smaller degree LR-CCSD). OSs calculated from TD-DFT with PCM solvent models are systematically larger than apparent OSs derived from experimental spectra. For example, fcomp from hybrid functionals and PCM have mean absolute errors that are ∼10% of n·fexp, where n is a solvent refractive index factor that arises from the energy flux of the radiation field in a dielectric (solvent). Theoretical cavity field corrections considering spherical cavities do not improve the agreement between computed and experimental data. Corrections that account for the molecular shape and the direction of transition dipole moments, or that explicitly account for the effect of solvent molecules on the local field, should be more appropriate.

How does theory compare to experiment for oscillator strengths in electronic spectra? Proposing range-separated hybrids with reliable accountability

As an important quantity in atomic and molecular spectroscopy, oscillator strength should be mentioned. Oscillator strength is linked to the transition dipole moment and consequently to the transition probability between two states, where its magnitude is directly connected to the intensity of the peaks in ultraviolet-visible spectra. However, accurately accounting for oscillator strengths still remains one of the greatest challenges in theory and experiment. Given previous efforts in the context of investigations into oscillator strengths, the related theoretical treatments are relatively limited and have proven to be challenging. In this work, the oscillator strengths in the electronic spectra of organic compounds have thoroughly been investigated with the help of optimally tuned range-separated hybrids (OT-RSHs). In particular, variants of the OT-RSHs combined with the polarizable continuum model (PCM), OT-RSHs-PCM, as well as their screened versions accounting for the screening effects by the electron correlation through the dielectric constant, OT-SRSHs-PCM, are proposed for reliable prediction of the oscillator strengths. The role of the involved ingredients in the proposed methods, namely the underlying density functional approximations, short-range and long-range Hartree-Fock (HF) exchange, as well as the range-separation parameter, has been examined in detail. It is shown that any combination of the parameters in the proposed approximations does not render the reliable oscillator strengths, but a particular compromise among them is needed to describe the experimental data well. Perusing all the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM, coupled with the linear response theory in the non-equilibrium solvation regime, with the correct asymptotic behavior and incorporating no (low) HF exchange contributions in the short-range part. The best proposed approximations also reveal superior performances not only with respect to their standard counterparts with the default parameters but also as compared to earlier range-separated functionals. Finally, the applicability of the best approximation is also put into broader perspective, where it is used for predicting the oscillator strengths in other sets of compounds not included in the process of developing the approximations. Hopefully, our proposed method can function as an affordable alternative to the expensive wave function-based methods for both theoretical modeling and confirming the experimental observations in the field of electronic spectroscopy.

Reference Vertical Excitation Energies for Transition Metal Compounds

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.

Gauge-Invariant Excited-State Linear and Quadratic Response Properties within the Meta-Generalized Gradient Approximation

Gauge invariance is a fundamental symmetry connected to charge conservation and is widely accepted as indispensable for any electronic structure method. Hence, the gauge variance of the time-dependent kinetic energy density τ used in many meta-generalized gradient approximations (MGGAs) to the exchange-correlation (XC) functional presents a major obstacle for applying MGGAs within time-dependent density functional theory (TDDFT). Replacing τ by the gauge-invariant generalized kinetic energy density τ̂ significantly improves the accuracy of various functionals for vertical excitation energies [R. Grotjahn, F. Furche, and M. Kaupp. J. Chem. Phys. 2022, 157, 111102]. However, the dependence of the resulting current-MGGAs (cMGGAs) on the paramagnetic current density gives rise to new exchange-correlation kernels and hyper-kernels ignored in previous implementations of quadratic and higher-order response properties. Here we report the first implementation of cMGGAs and hybrid cMGGAs for excited-state gradients and dipole moments, as well as an extension to quadratic response properties including dynamic hyperpolarizabilities and two-photon absorption cross sections. In the first comprehensive benchmark study of MGGAs and cMGGAs for two-photon absorption cross sections, the M06-2X functional is found to be superior to the GGA hybrid PBE0. Additionally, two case studies from the literature for the practical prediction of nonlinear optical properties are revisited and potential advantages of hybrid (c)MGGAs compared to hybrid GGAs are discussed. The effect of restoring gauge invariance varies depending on the employed MGGA functional, the type of excitation, and the property under investigation: While some individual excited-state equilibrium structures are significantly affected, on average, these changes result in marginal improvements when compared against high-level reference data. Although the gauge-variant MGGA quadratic response properties are generally close to their gauge-invariant counterparts, the resulting errors are not bounded and significantly exceed typical method errors in some of the cases studied. Despite the limited effects seen in benchmark studies, gauge-invariant implementations of cMGGAs for excited-state properties are desirable from a fundamental perspective, entail little additional computational cost, and are necessary for response properties consistent with cMGGA linear response calculations such as excitation energies.

Double Hybrid Density Functionals for the Electronic Excitation Energies of Linear Cyanines

The lowest bright electronic excitations of seven model linear cyanines (CN3-CN15) using 28 double-hybrid (DH) density functionals are benchmarked against accurate and recent CC3 results. Some of these DH functionals are recently designed specifically for excited electronic state calculations. In addition, CIS, CIS(D), SCS-CIS(D), and SOS-CIS(D) were also tested. Four different basis sets were used for the vertical electronic excitation calculations: cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis. Augmented basis sets (e.g. aug-cc-pVDZ and aug-cc-pVTZ) are found to be required for accurate and consistent results using DH functionals. The DH functionals tested in this work are classified into four main groups: global double-hybrids (GDH), range-separated double-hybrids (RSDH), spin-component and spin-opposite scaling global double-hybrids (SCS/SOS-GDH), and spin-component and spin-opposite scaling range-separated double-hybrids (SCS/SOS-RSDH). Within these groups, the SCS/SOS-RSDH group of functionals is found to provide the lowest mean absolute error (MAE) values (in the range 0.020-0.148 eV) in comparison with the GDH group (0.195-0.441 eV), the RSDH group (0.186-0.511 eV), and the SCS/SOS-GDH group (0.079-0.461 eV). Of all the DH functionals and ab initio methods investigated in the present contribution, the following functionals are found to be the most accurate and consistent: SCS-ωB2GPPLYP (MAE = 0.036 eV), SOS-ωB2GPPLYP (MAE = 0.020 eV), SOS-ωB88PP86 (MAE = 0.035 eV), and SOS-ωPBEPP86 (MAE = 0.037 eV). In general, the ab initio methods tested here show mediocre performance as compared to many DH functionals.

Benchmarking CASPT3 Vertical Excitation Energies

Based on 280 reference vertical transition energies of various natures (singlet, triplet, valence, Rydberg, n → π∗, π → π∗, 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.

Generalized Formulation of the Density Functional Tight Binding-Based Restricted Ensemble Kohn-Sham Method with Onsite Correction to Long-Range Correction

We present a generalized formulation for the combination of the density functional tight binding (DFTB) approach and the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method by considering onsite correction (OC) as well as the long-range corrected (LC) functional. The OC contribution provides more accurate energies and analytic gradients for individual microstates, while the multireference character of the SSR provides the correct description for conical intersections. We benchmark the LC-OC-DFTB/SSR method against various DFTB calculation methods for excitation energies and conical intersection structures with π/π* or n/π* characters. Furthermore, we perform excited-state molecular dynamics simulations with a molecular rotary motor with variations of LC-OC-DFTB/SSR approaches. We show that the OC contribution to the LC functional is crucial to obtain the correct geometry of conical intersections.

A theoretical investigation into gallic acid pyrolysis

Thermodynamic and kinetic information on the first two steps of gallic acid pyrolysis, a decarboxylation followed by a dehydrogenation, is obtained based on density functional theory and quantum chemistry. For the kinetics, transition states are identified with the help of the climbing image nudged elastic band method. Both reactions exhibit two transition states. One of them is related to the rotation of OH groups, and the other one is related to the breaking and forming of bonds. The gallic acid pyrolysis as a whole is judged to be endothermal, and it changes from endergonic to exergonic between 500 and 750 K. The second reaction, the dehydrogenation of pyrogallol, is identified as the rate-determining step of gallic acid pyrolysis, with reaction rate constants below 1 s^-1 for temperatures below 1250 K.

OS100: A Benchmark Set of 100 Digitized UV-Visible Spectra and Derived Experimental Oscillator Strengths

Excited-state quantum chemical calculations usually report excitation energies and oscillator strengths, f, for each electronic transition. On the other hand, UV-visible spectrophotometric experiments measure energy-dependent molar extinction/attenuation coefficients, ε(v), that give absorption band line shapes when plotted. ε(v) and f are related, but this relation is complicated by broadening and solvation effects. We fitted and integrated 100 experimental UV-visible spectra to obtain 164 f_exp values for absorption bands appearing in these spectra. The 100 UV-visible spectra belong to solvated organic molecules ranging in size from 6-34 atoms. We estimated uncertainties in the fitting to indicate confidence level in the reported f_exp values. The corresponding computed oscillator strengths (f_comp) were obtained with time-dependent density functional theory and a polarizable continuum solvent model. By expressing experimental and computed absorption strengths using a common quantity, we directly compared f_comp and f_exp. Although f_comp and f_exp are well correlated (linear regression R² = 0.921), f_comp in most cases overestimated f_exp (regression slope = 1.34). The agreement between absolute f_comp and f_exp values was substantially improved by accounting for a solvent refractive index factor, as suggested in some derivations in the literature. The 100 digitized UV-visible spectra are included as plain text files in the Supporting Information to aid in benchmarking computational or machine learning methods that aim to simulate realistic UV-visible absorption spectra.

Laser Photodissocation, Action Spectroscopy and Mass Spectrometry Unite to Detect and Separate Isomers

The separation and detection of isomers remains a challenge for many areas of mass spectrometry. This article highlights laser photodissociation and ion mobility strategies that have been deployed to tackle...

Optimizing the Cosensitization Effect of SQ02 Dye on BP-2 Dye-Sensitized Solar Cells: A Computational Quantum Chemical Study

Cosensitization of the semiconducting electrode in dye-sensitized solar cells (DSCs), with two or more light-harvesting dyes, is a chemical fabrication method that aims to achieve a panchromatic absorption spectrum emulating that of the solar emission spectrum. In this paper, SQ02 and BP-2 cosensitizers have been investigated, as isolated monomers/dimer and adsorbed monomers/dimer on the TiO₂ (101) anatase surface, by employing density functional theory (DFT) and time-dependent DFT calculations. Computed results showed that the dominant electron injection pathway is direct injection from each dye into the conduction band of TiO₂. The almost complete spectral overlap between the simulated absorption spectrum of BP-2 and fluorescence emissions of SQ02 implies that excitation energy transfer occurs between cosensitizers via the trivial reabsorption mechanism. However, the results showed very limited unidirectional intermolecular charge transfer (CT) from SQ02 dye to BP-2 dye (0.04 |e^-|). Therefore, this study also presents a stepwise molecular engineering of BP-2 dye, aiming at optimizing the cosensitization functionality. First, 14 redesigned dye candidates are reported to identify dyes with photophysical properties matching the requirements for efficient DSCs. Second, the four most promising dyes are shortlisted for testing as cosensitizers with the SQ02 dye. The molecular design factors of cosensitization that need validation are chemical compatibility, availability of CT between cosensitizers, and complementarity of the absorption spectra. This screening suggests the judicious choice of the modeled difluorenyl amine donor-based dye (BP-D4) as a very promising cosensitizer. In particular, the SQ02/BP-D4 dimer showed 10 times larger (0.53 |e^-|) unidirectional CT than that of SQ02/BP-2 dimer, in addition to the maximum increased electron population of acceptor moieties upon photoexcitation.

Time-Dependent Long-Range-Corrected Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling: A Comprehensive Analysis of Singlet-Singlet and Singlet-Triplet Excitation Energies

Following the work on spin-component and spin-opposite scaled (SCS/SOS) global double hybrids for singlet-singlet excitations by Schwabe and Goerigk [ J. Chem. Theory Comput. 2017, 13, 4307-4323] and our own works on new long-range corrected (LC) double hybrids for singlet-singlet and singlet-triplet excitations [ J. Chem. Theory Comput. 2019, 15, 4735-4744 and J. Chem. Phys. 2020, 153, 064106], we present new LC double hybrids with SCS/SOS that demonstrate further improvement over previously published results and methods. We introduce new unscaled and scaled versions of different global and LC double hybrids based on Becke88 or PBE exchange combined with LYP, PBE, or P86 correlation. For singlet-singlet excitations, we cross-validate them on six benchmark sets that cover small to medium-sized chromophores with different excitation types (local-valence, Rydberg, and charge transfer). For singlet-triplet excitations, we perform the cross-validation on three different benchmark sets following the same analysis as in our previous work in 2020. In total, 203 excitations are analyzed. Our results confirm and extend those of Schwabe and Goerigk regarding the superior performance of SCS and SOS variants compared to their unscaled parents by decreasing mean absolute deviations, root-mean-square deviations, or error spans by more than half and bringing absolute mean deviations closer to zero. Our SCS/SOS variants are shown to be highly efficient and robust for the computation of vertical excitation energies, which even outperform specialized double hybrids that also contain an LC in their perturbative part. In particular, our new SCS/SOS-ωPBEPP86 and SCS/SOS-ωB88PP86 functionals are four of the most accurate and robust methods tested in this work, and we fully recommend them for future applications. However, if the relevant SCS and SOS algorithms are not available to the user, we suggest ωPBEPP86 as the best unscaled method in this work.

Pairing double hybrid functionals with a tailored basis set for an accurate thermochemistry of hydrocarbons

A collection of five challenging datasets, including noncovalent interactions, reaction barriers and electronic rearrangements of medium-sized hydrocarbons, has been selected to verify the robustness of double-hybrid functionals used in conjunction with the small DH-SVPD basis set, especially developed for noncovalent interactions. The analysis is completed by other, more standard functionals, for a total of 17 models, including also empirical corrections for dispersion. The obtained results show that the chemical accuracy threshold, that is an error lower than 1.0 kcal mol⁻¹, can be obtained by pairing the nonempirical PBE-QIDH functional with the DH-SVPD basis set, as well as by other semi-empirical functionals, such as DSD-PBEP86, using larger basis sets and empirical corrections. More in general, a significant improvement can be obtained using the DH-SVPD basis set with DHs, without resorting to any empirical corrections. This choice leads to a fast computational protocol that, avoiding any empirical potential, remains on a fully quantum ground.

The pairing of the PBE-QIDH double-hybrid functional with a tailored split-valence basis set leads to a fast computational protocol for the accurate evaluation of hydrocarbon thermochemistry, without resorting to any empirical correction.

Systematic improvement of molecular excited state calculations by inclusion of nuclear quantum motion: A mode-resolved picture and the effect of molecular size

The energies of molecular excited states arise as solutions to the electronic Schrödinger equation and are often compared to experiment. At the same time, nuclear quantum motion is known to be important and to induce a redshift of excited state energies. However, it is thus far unclear whether incorporating nuclear quantum motion in molecular excited state calculations leads to a systematic improvement of their predictive accuracy, making further investigation necessary. Here, we present such an investigation by employing two first-principles methods for capturing the effect of quantum fluctuations on excited state energies, which we apply to the Thiel set of organic molecules. We show that accounting for zero-point motion leads to much improved agreement with experiment, compared to "static" calculations that only account for electronic effects, and the magnitude of the redshift can become as large as 1.36 eV. Moreover, we show that the effect of nuclear quantum motion on excited state energies largely depends on the molecular size, with smaller molecules exhibiting larger redshifts. Our methodology also makes it possible to analyze the contribution of individual vibrational normal modes to the redshift of excited state energies, and in several molecules, we identify a limited number of modes dominating this effect. Overall, our study provides a foundation for systematically quantifying the shift of excited state energies due to nuclear quantum motion and for understanding this effect at a microscopic level.

Global double hybrids do not work for charge transfer: A comment on "Double hybrids and time-dependent density functional theory: An implementation and benchmark on charge transfer excited states"

We comment on the results published by Ottochian et al. in J. Comput Chem. 2020, 41, 1242. Therein, the authors claim that the second-order, perturbative correlation correction applied to the time-dependent version of the PBE-QIDH global double-hybrid functional approximation (DHDFA) enables the description of charge-transfer (CT) excitations. Herein, we point out some inadvertent oversights related to what had already been previously known and achieved in the field of time-dependent DHDFAs. Exemplified for the same four systems that Ottochian et al. have used to analyze intermolecular CT excitations, we show how a systematic and unacceptably large redshift in global DHDFAs is rectified when using the latest long-range corrected DHDFAs published earlier in J. Chem. Theory Comput. 2019, 15, 4735.

The Trip to the Density Functional Theory Zoo Continues: Making a Case for Time-Dependent Double Hybrids for Excited-State Problems

This account is written for general users of time-dependent density functional theory (TD-DFT) methods as well as chemists who are unfamiliar with the field. It includes a brief overview of conventional TD-DFT approaches and recommendations for applications to organic molecules based on our own experience. The main emphasis of this work, however, lies in providing the first in-depth review of time-dependent double-hybrid density functionals. They were first established in 2007 with very promising follow-up studies in the subsequent four years before developments or applications became scarce. The topic has regained more interest since 2017, and this account reviews those latest developments led by our group. These developments have shown unprecedented robustness for a variety of different types of electronic excitations when compared to more conventional TD-DFT methods. In particular, time-dependent double hybrids do not suffer from artificial ghost states and are able to reproduce exciton-coupled absorption spectra. Our latest methods include range separation and belong to the currently best TD-DFT methods for singlet-singlet excitations in organic molecules. While there is still room for improvement and further development in this space, we hope that this account encourages users to adjust their computational protocols to such new methods to provide more real-life testing and scenarios.

Dynamical correction to the Bethe-Salpeter equation beyond the plasmon-pole approximation

The Bethe-Salpeter equation (BSE) formalism is a computationally affordable method for the calculation of accurate optical excitation energies in molecular systems. Similar to the ubiquitous adiabatic approximation of time-dependent density-functional theory, the static approximation, which substitutes a dynamical (i.e., frequency-dependent) kernel by its static limit, is usually enforced in most implementations of the BSE formalism. Here, going beyond the static approximation, we compute the dynamical correction of the electron-hole screening for molecular excitation energies, thanks to a renormalized first-order perturbative correction to the static BSE excitation energies. The present dynamical correction goes beyond the plasmon-pole approximation as the dynamical screening of the Coulomb interaction is computed exactly within the random-phase approximation. Our calculations are benchmarked against high-level (coupled-cluster) calculations, allowing one to assess the clear improvement brought by the dynamical correction for both singlet and triplet optical transitions.

Relationships between Orbital Energies, Optical and Fundamental Gaps, and Exciton Shifts in Approximate Density Functional Theory and Quasiparticle Theory

The relationships between Kohn-Sham (KS) and generalized KS (GKS) density functional orbital energies and fundamental gaps or optical gaps raise many interesting questions including the physical meanings of KS and GKS orbital energies when computed with currently available approximate density functionals (ADFs). In this work, by examining three diverse databases with various ADFs, we examine such relations from the point of view of the exciton shift of quasiparticle theory. We start by calculating a large number of excitation energies by time-dependent density functional theory (TDDFT) with a large number of ADFs. To relate the exciton shift implicit in TDDFT to the exciton shift that is explicit in Green's function theory, we define the exciton shift in TDDFT as the difference of the response shift and the quasiparticle shift. We found a strong correlation between the response shift and the amount of Hartree-Fock exchange included in the density functional, with the response shift varying between -1 and 5 eV. This range is an order of magnitude larger than the mean errors of the TDDFT excitation energies. This result suggests that, with currently available functionals, the KS or GKS orbital energies should be treated as intermediate mathematical variables in the calculation of excitation energies rather than as the energies of independent-particle reference states for quasiparticle theory.

A New Benchmark Set for Excitation Energy of Charge Transfer States: Systematic Investigation of Coupled Cluster Type Methods

The numerous existing publications on benchmarking quantum chemistry methods for excited states rarely include Charge Transfer (CT) states, although many interesting phenomena in, e.g., biochemistry and material physics involve the transfer of electrons between fragments of the system. Therefore, it is timely to test the accuracy of quantum chemical methods for CT states, as well. In this study we first propose a new benchmark set consisting of dimers having low-energy CT states. On this set, the vertical excitation energy has been calculated with Coupled Cluster methods including triple excitations (CC3, CCSDT-3, CCSD(T)(a)*), as well as with methods including full or approximate doubles (CCSD, STEOM-CCSD, CC2, ADC(2), EOM-CCSD(2)). The results show that the popular CC2 and ADC(2) methods are much less accurate for CT states than for valence states. On the other hand, EOM-CCSD seems to have similar systematic overestimation of the excitation energies for both types of states. Among the triples methods the novel EOM-CCSD(T)(a)* method including noniterative triple excitations is found to stand out with its consistently good performance for all types of states, delivering essentially EOM-CCSDT quality results.

The Quest for Highly Accurate Excitation Energies: A Computational Perspective

We provide an overview of the successive steps that made it 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 the original 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 the fair and accurate assessment of the performance of computationally lighter methods. We conclude by discussing possible future theoretical and technological developments in the field.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Medium Sized Molecules

Following our previous work focusing on compounds containing up to 3 non-hydrogen atoms [J. Chem. Theory Comput. 2018, 14, 4360-4379], we present here highly accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms: acetone, acrolein, benzene, butadiene, cyanoacetylene, cyanoformaldehyde, cyanogen, cyclopentadiene, cyclopropenone, cyclopropenethione, diacetylene, furan, glyoxal, imidazole, isobutene, methylenecyclopropene, propynal, pyrazine, pyridazine, pyridine, pyrimidine, pyrrole, tetrazine, thioacetone, thiophene, thiopropynal, and triazine. To obtain these energies, we use equation-of-motion/linear-response coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ) and selected configuration interaction (SCI) calculations (with tens of millions of determinants in the reference space), as well as the multiconfigurational n-electron valence state perturbation theory (NEVPT2) method. All these approaches are applied in combination with diffuse-containing atomic basis sets. For all transitions, we report at least CC3/aug-cc-pVQZ vertical excitation energies as well as CC3/aug-cc-pVTZ oscillator strengths for each dipole-allowed transition. We show that CC3 almost systematically delivers transition energies in agreement with higher-level methods with a typical deviation of ±0.04 eV, except for transitions with a dominant double excitation character where the error is much larger. The present contribution gathers a large, diverse, and accurate set of more than 200 highly accurate transition energies for states of various natures (valence, Rydberg, singlet, triplet, n → π*, π → π*, ...). We use this series of theoretical best estimates to benchmark a series of popular methods for excited state calculations: CIS(D), ADC(2), CC2, STEOM-CCSD, EOM-CCSD, CCSDR(3), CCSDT-3, CC3, and NEVPT2. The results of these benchmarks are compared to the available literature data.

Is ADC(3) as Accurate as CC3 for Valence and Rydberg Transition Energies?

The search for new models rapidly delivering accurate excited-state energies and properties is one of the most active research lines of theoretical chemistry. Along with these developments, the performance of known methods is constantly reassessed on the basis of new benchmark values. In this Letter, we show that the third-order algebraic diagrammatic construction, ADC(3), does not yield transition energies of the same quality as the third-order coupled cluster method, CC3. This is demonstrated by extensive comparisons with several hundred high-quality vertical transition energies obtained with FCI, CCSDTQ, and CCSDT. Direct comparisons with experimental 0-0 energies of small- and medium-size molecules support the same conclusion, which holds for both valence and Rydberg transitions. Considering these results, we introduce a composite approach, ADC(2.5), which consists of averaging the ADC(2) and ADC(3) excitation energies. Although ADC(2.5) does not match the CC3 accuracy, it significantly improves the ADC(3) results, especially for vertical energies.

Excited States with Selected Configuration Interaction-Quantum Monte Carlo: Chemically Accurate Excitation Energies and Geometries

We employ quantum Monte Carlo to obtain chemically accurate vertical and adiabatic excitation energies, and equilibrium excited-state structures for the small, yet challenging, formaldehyde and thioformaldehyde molecules. A key ingredient is a robust protocol to obtain balanced ground- and excited-state Jastrow-Slater wave functions at a given geometry, and to maintain such a balanced description as we relax the structure in the excited state. We use determinantal components generated via a selected configuration interaction scheme which targets the same second-order perturbation energy correction for all states of interest at different geometries, and fully optimize all variational parameters in the resultant Jastrow-Slater wave functions. Importantly, the excitation energies as well as the structural parameters in the ground and excited states are converged with very compact wave functions comprising few thousand determinants in a minimally augmented double-ζ basis set. These results are obtained already at the variational Monte Carlo level, the more accurate diffusion Monte Carlo method yielding only a small improvement in the adiabatic excitation energies. We find that matching Jastrow-Slater wave functions with similar variances can yield excitation energies compatible with our best estimates; however, the variance-matching procedure requires somewhat larger determinantal expansions to achieve the same accuracy, and it is less straightforward to adapt during structural optimization in the excited state.

Application of the Consonance Solvent Concept for Accurate Prediction of Buckminster Solubility in 180 Net Solvents using COSMO-RS Approach

The default COSMO-RS (Conductor like Screening Model for Real Solvents) approach is incapable of accurate computation of C60 solubility in net solvents. Additionally, there is no adequate selection of single or multiple reference solvent, which can be applied to the whole population of 180 solvents for improving prediction of mole fraction at saturated conditions. This failure cannot be addressed to inaccurate data of the Buckminster fusion, although they pose a challenge for experimental measurement due to intense sublimation of C60 at elevated temperatures and the possibility of solvates precipitation. However, taking advantage of the richness of experimental data of fullerene solubility, it is possible to identify the source of errors expressed in terms of fluidization affinity. Classification of solvents according to the value of this fluidization term allowed for formulation of a consonance solvents approach, which enables accurate prediction of C60 solubility using the single reference solvent method.