Following our previous work focussing on compounds containing up to 3 non-hydrogen atoms [\emph{J. Chem. Theory Comput.} {\bfseries 14} (2018) 4360--4379], we present here highly-accurate vertical transition energies obtained for 27 molecules encompassing 4, 5, and 6 non-hydrogen atoms. To obtain these energies, we use equation-of-motion coupled cluster theory up to the highest technically possible excitation order for these systems (CC3, EOM-CCSDT, and EOM-CCSDTQ), 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/\emph{aug}-cc-pVQZ vertical excitation energies as well as CC3/\emph{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 $\pm 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 \rightarrow \pi^*$, $\pi \rightarrow \pi^*$, \ldots). 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, as well as NEVPT2. The results of these benchmarks are compared to the available literature data.

Lire le billet

Quantum Monte Carlo (QMC) methods have been applied very successfully to ground state properties but still remain generally less effective than other non-stochastic methods for electronically excited states. Nonetheless, we have recently reported accurate excitation energies for small organic molecules at the fixed-node diffusion Monte Carlo (FN-DMC) within a Jastrow-free QMC protocol relying on a deterministic and systematic construction of nodal surfaces using the selected configuration interaction (sCI) algorithm known as CIPSI (Configuration Interaction using a Perturbative Selection made Iteratively). Albeit highly accurate, these all-electron calculations are computationally expensive due to the presence of core electrons. One very popular approach to remove these chemically-inert electrons from the QMC simulation is to introduce pseudopotentials. However, such an approach inevitably introduces a bias due to the approximate nature of these pseudopotentials. Furthermore, an additional bias may be introduced in DMC due to the localization of nonlocal pseudopotentials. Taking the water molecule as an example, we investigate the influence of pseudopotentials on vertical excitation energies obtained with sCI and FN-DMC methods. Although pseudopotentials are known to be relatively safe for ground state properties, we evidence that special care is required if one strives for highly accurate vertical transition energies. Indeed, comparing all-electron and valence-only calculations, we show that the approximate nature of the pseudopotentials can induce errors as large as 0.1 eV on the excitation energies. While acceptable for most chemical applications, it might become unacceptable for benchmark studies.

Lire le billet

Interest in fullerenes has been renewed recently in astrophysics as a consequence of their detection in circumstellar environments. In particular, C60 + was detected in the diffuse interstellar medium and its presence has been related to some diffuse interstellar bands (DIBs) whose origin was previously unknown. A single recent laboratory experiment (J. Phys. Chem. A 2017, 121, 7356-7361) shows that upon laser excitation at 785 nm, C60 + in neon matrices exhibits a radiative decay at 965 nm, while UV photoexcitation does not lead to any significant luminescence. To rationalize this original dual photophysical behavior, we have performed time-dependent density functional theory (TD-DFT) calculations on C60 + to investigate the potential energy surfaces of the relevant electronic states, completed by the simulations of vibrationally-resolved absorption and emission spectra. The proposed photophysical pathways shed light on the experimental measurements: The near-IR laser excitation populates the eleventh doublet excited state (D11) that decays to the lowest first bright excited state D5, from which photoluminescence is predicted. Indeed, D5 is largely separated from the lower electronic states (D0-D4). Thus, D5 behaves effectively as the first excited state, while the D0-D4 set of states act as the electronic ground state. In addition, there are no low-lying conical intersections between D5 and lower excited states energetically accessible upon near-IR excitation that can provide efficient nonradiative decay channels for this state, leaving radiative decay as the most likely deactivation pathway. However, a sloped conical intersection between D5 and D4 was located around 2.9 eV above D0. While it is too high in energy to be accessible upon near-IR excitation, it provides a funnel for efficient nonradiative decay down to the ground state (D0) accessible upon UV light excitation. Thus, the photophysics of C60+ is controlled by the ability to access this funnel: Upon near-IR excitation, the system fluoresces because the funnel for nonradiative decay cannot be reached, while UV irradiation provides a different route by opening up a radiationless decay channel via this funnel accounting for the absence of fluorescence.

Lire le billet

The self-sensitized photo-oxygenation and singlet oxygen thermal release mechanisms in a pyridinium-substituted dimethyldihydropyrene (DHP) derivative have been investigated using quantum chemical calculations. First, the main photophysical pathway for intersystem crossing was identified, allowing the production of a DHP triplet state. Second, the energy transfer pathway between this triplet state and the oxygen triplet ground state was computed revealing a very efficient route for the photosensitized generation of singlet oxygen. Finally, the thermal pathway for the formation of a metacyclophanediene endoperoxide and the singlet oxygen release was characterized. A concerted and a stepwise mechanism were identified, the first one being lower in energy. All these results are consistent with recent experimental results and confirm that this type of system could be attractive oxygen carriers and singlet oxygen delivery agents.

Lire le billet

Excited states exhibiting double excitation character are notoriously difficult to model using conventional single-reference methods, such as adiabatic time-dependent density-functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC). 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 diffuse-containing atomic basis sets. Our set gathers 20 vertical transitions from 14 small- and medium-size molecule. Depending on the size of the molecule, selected configuration interaction (sCI) and/or multiconfigurational 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 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 \emph{trans}-butadiene) involving a significant amount of singles. As expected, CC approaches including quadruples yield highly accurate results for any type of transitions. The quality of the excitation energies obtained with multiconfigurational 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. The inclusion of the $\sigma$ and $\sigma^*$ orbitals in the active space, even for transitions involving mostly $\pi$ and $\pi^*$ orbitals, is mandatory in order to reach high accuracy. A theoretical best estimate (TBE) is reported for each transition.

Lire le billet