A mountaineering strategy to excited states: Accurate vertical transition energies and benchmarks for substituted benzenes

In an effort to expand the existing QUEST database of accurate vertical transition energies [Véril et al. WIREs Comput. Mol. Sci. 2021, 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)  ⋆ , 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 very 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 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.

Spectroscopic Quantification of the Inverted Singlet-Triplet Gap in Pentaazaphenalene

We report the direct and accurate spectroscopic quantification of the inverted singlet-triplet gap in 1,3,4,6,9b-pentaazaphenalene. This measurement is achieved by directly probing the lowest singlet and triplet states via high-resolution cryogenic anion photoelectron spectroscopy. The assignment of the first excited singlet state is confirmed by visible absorption spectroscopy in an argon matrix at 20 K. Our measurements yield an inverted singlet-triplet gap with ΔEST= -0.047(7) eV. The accurate quantification of the singlet-triplet gap presented here allows for direct evaluation of various computational electronic structure methods and highlights the critical importance of the proper description of the double excitation character of these electronic states. Overall, this study validates the idea that despite Hund's multiplicity rule, useful organic chromophores can have inherently inverted singlet-triplet gaps.

Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels

We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.

Role of Singles Amplitudes in ADC(2) and CC2 for Low-Lying Electronically Excited States

The closely related second-order methods CC2 and ADC(2) usually perform very similarly for single excitations of organic molecules. However, as rationalized in this work, significant deviations between these two methods can arise if the ground state and a low-lying singly excited state arise from a strong coupling between their leading configurations. Such a configuration mixing is partially accounted for in CC2 through the ground-state singles amplitudes but is omitted in ADC(2). This can cause unusual deviations between the results obtained with these methods. In this work, we study how severe this effect can become at the example of two solvatochromic dyes: the negatively solvatochromic betaine dye N1-tBu and the positively solvatochromic bithiophene P1. These two dyes allow one to study the limits of both small and somewhat larger excitation energies and configuration mixing by tuning the S0 → S1 transition energy through the polarity of the environment. Higher-level calculations at the CC3 level provide information on the accuracy of ADC(2) and CC2 in these cases. The most extreme deviation between ADC(2) and CC2 is found for N1-tBu in vacuum, where the ADC(2) result is 0.45 eV below that of CC2. In this case, the methodical error of CC2 with respect to CC3 is only 0.05 eV. With increasing excitation energy in polar solvents, the CC2-ADC(2) deviation decreases and reaches a value of only 0.15 eV. For P1, which has larger excitation energies, these effects are reversed due to the opposite solvatochromism but also smaller in magnitude: the deviation increases from 0.08 eV in vacuum to 0.16 eV in the so-called conductor limit of the continuum solvation model. Although for these two dyes larger deviations are observed for smaller excitation energies, the extent of configuration mixing does not generally correlate with only the size of excitation energy. For example, s-triazine (0.15 eV), formamide (0.19 eV), and formaldehyde (0.23 eV) also show large deviations between CC2 and ADC(2) despite their much higher excitation energies compared to those of N1-tBu and P1.

DELFI: a computer oracle for recommending density functionals for excited states calculations

Density functional theory (DFT) is the workhorse of computational quantum chemistry. One of its main limitations is that choosing the right functional is a non-trivial task left for human experts. The choice is particularly hard for excited state calculations when using its time-dependent formulation (TD-DFT). This is due to the approximations of the method, but also because the photophysical properties of a molecule are defined by a manifold of states that all need to be properly described. This includes not only the relative energy of the states, but also capturing the correct character, order, and intensity of the transitions. In this work, we developed a neural network to recommend functionals to be used on molecules for TD-DFT calculations, by simultaneously considering all these properties for a manifold of states. This was possible by developing a scoring system to define the accuracy of an excited state's calculation against a higher-accuracy reference. The scoring system is generalizable to any level of theory; we here applied it to evaluate the performance of common functionals of different rungs against a higher accuracy method on a large set of organic molecules. The results are collected in a database that we released and made open, providing four million data points to the community for future applications. The scoring system assigns a value between zero and one hundred to each functional for each molecule, transforming the complicated task of learning photophysical properties into a simpler regression task. We used the dataset to train a graph attention neural network to predict the scores for unseen molecules. We call this oracle DELFI (Data-driven EvaLuation of Functionals by Inference), which can be used to quickly screen and predict the ranking of functionals to calculate the optical properties of organic molecules. We validated DELFI in two in silico experiments: choosing a common functional for a series of spiropyran-merocyanine isomers and a unique functional to screen a large dataset of over 50 000 organic photovoltaic molecules, for which an extensive benchmark would be unfeasible. A corresponding web application allows DELFI to be easily run and the results to be analyzed, alleviating the hurdle of choosing the right functional for TD-DFT calculations.

Cluster-Based Approach Utilizing Optimally Tuned TD-DFT to Calculate Absorption Spectra of Organic Semiconductor Thin Films

The photophysics of organic semiconductor (OSC) thin films or crystals has garnered significant attention in recent years since a comprehensive theoretical understanding of the various processes occurring upon photoexcitation is crucial for assessing the efficiency of OSC materials. To date, research in this area has relied on methods using Frenkel-Holstein Hamiltonians, calculations of the GW-Bethe-Salpeter equation with periodic boundaries, or cluster-based approaches using quantum chemical methods, with each of the three approaches having distinct advantages and disadvantages. In this work, we introduce an optimally tuned, range-separated time-dependent density functional theory approach to accurately reproduce the total and polarization-resolved absorption spectra of pentacene, tetracene, and perylene thin films, all representative OSC materials. Our approach achieves excellent agreement with experimental data (mostly ≤0.1 eV) when combined with the utilization of clusters comprising multiple monomers and a standard polarizable continuum model to simulate the thin-film environment. Our protocol therefore addresses a major drawback of cluster-based approaches and makes them attractive tools for OSC investigations. Its key advantages include its independence from external, system-specific fitting parameters and its straightforward application with well-known quantum chemical program codes. It demonstrates how chemical intuition can help to reduce computational cost and still arrive at chemically meaningful and almost quantitative results.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Heptazine, Cyclazine, and Related Compounds: Chemically-Accurate Estimates of the Inverted Singlet-Triplet Gap

Molecules that violate Hund's rule and exhibit an inverted gap between the lowest singlet S1 and triplet T1 excited states have attracted considerable attention due to their potential applications in optoelectronics. Among these molecules, the triangular-shaped heptazine, and its derivatives, have been in the limelight. However, conflicting reports have arisen regarding the relative energies of S1 and T1. 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 S1-T1 gap. We also determined the 0-0 energies to evaluate the significance of the vertical approximation. In addition, we compute reference S1-T1 gaps for a series of 10 related molecules. This enables us to benchmark lower-order methods for future applications in larger systems within the same family of compounds. This contribution can serve as a foundation for the design of triangular-shaped molecules with enhanced photophysical properties.

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.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

Performance of Common Density Functionals for Excited States of Tetraphenyldibenzoperiflanthene

Time-dependent density functional theory is the method of choice to efficiently calculate excitation spectra with the functional and basis set choice allowing one to compromise between accuracy and computational cost. In this work, the performance of different functionals as well as the second-order approximate coupled cluster singles and doubles model CC2 is evaluated by comparing the results to experimental results of the example molecule tetraphenyldibenzoperiflanthene (DBP). The choice of the functional has a significant impact on the calculated spectrum of DBP. The performance of a number of different functionals was evaluated, quantified, and, where possible, discussed. The best functional, tuned-CAM-B3LYP, is used to investigate DBP on a surface of hexagonal boron nitride (h-BN). The resulting spectrum shows excellent agreement with experimental results for a monolayer of DBP on h-BN.

Theoretical Study on the Photoacidity of Hydroxypyrene Derivatives in DMSO Using ADC(2) and CC2

This work applies the thermodynamic Förster cycle to theoretically investigate the pK_a^*, i.e., excited-state pK_a values of pyranine-derived superphotoacids developed by Jung and co-workers. The latter photoacids are strong enough to transfer a proton to the aprotic solvent dimethyl sulfoxide (DMSO). The Förster cycle provides access to pK_a^* via the ground-state pK_a and the electronic excitation energies. We use the conductor-like screening model for real solvents (COSMO-RS) to compute the ground-state pK_a and the correlated wavefunction-based methods ADC(2) and CC2 with the continuum solvation model COSMO to calculate the pK_a change upon excitation. A comparison of the calculated UV/Vis absorption and fluorescence emission energies to the experimental results leads us to infer that this approach allows for a proper description of the electronic excitations. In particular, implicit solvation by means of the COSMO model appears to be sufficient for the treatment of these photoacids in DMSO. The calculations confirm the presumption that a charge redistribution from the hydroxy group to the aromatic ring and the electron-withdrawing substituents is the origin of photoacidity for these photoacids. Moreover, the calculations with the continuum solvation model predict that the pK_a jump upon excitation decreases with increasing solvent polarity, as rationalized based on the Förster cycle.

On the origin of the inverted singlet-triplet gap of the 5th generation light-emitting molecules

Excitation energies of the lowest singlet and triplet state of molecules whose first excited singlet state lies energetically below the first triplet state have been studied computationally at (time-dependent) density functional theory, coupled-cluster, and second-order multiconfiguration perturbation theory levels. The calculations at the ab initio levels show that the singlet-triplet gap is inverted as compared to the one expected from Hund's rule, whereas all density functionals yield the triplet state as the lowest excited state. Double excitations responsible for the inverted singlet-triplet gap have been identified. Employing the spin-flip and ΔSCF methods, singlet-triplet inversion was obtained at the density functional theory level for some of the studied molecules.

Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models

Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for Time-Dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.

Charge-Transfer Excitations within Density Functional Theory: How Accurate Are the Most Recommended Approaches?

The performance of the most recent density functionals is assessed for charge-transfer (CT) excitations using comprehensive intra- and intermolecular CT benchmark sets with high-quality reference values. For this comparison, the state-of-the-art range-separated (RS) and long-range-corrected (LC) double hybrid (DH) approaches are selected, and global DH and LC hybrid functionals are also inspected. The correct long-range behavior of the exchange–correlation (XC) energy is extensively studied, and various CT descriptors are compared as well. Our results show that the most robust performance is attained by RS-PBE-P86/SOS-ADC(2), as it is suitable to describe both types of CT excitations with outstanding accuracy. Furthermore, concerning the intramolecular transitions, unexpectedly excellent results are obtained for most of the global DHs, but their limitations are also demonstrated for bimolecular complexes. Despite the outstanding performance of the LC-DH methods for common intramolecular excitations, serious deficiencies are pointed out for intermolecular CT transitions, and the wrong long-range behavior of the XC energy is revealed. The application of LC hybrids to such transitions is not recommended in any respect.

On the Photostability of Cyanuric Acid and Its Candidature as a Prebiotic Nucleobase

Cyanuric acid is a triazine derivative that has been identified from reactions performed under prebiotic conditions and has been proposed as a prospective precursor of ancestral RNA. For cyanuric acid to have played a key role during the prebiotic era, it would have needed to survive the harsh electromagnetic radiation conditions reaching the Earth’s surface during prebiotic times (≥200 nm). Therefore, the photostability of cyanuric acid would have been crucial for its accumulation during the prebiotic era. To evaluate the putative photostability of cyanuric acid in water, in this contribution, we employed density functional theory (DFT) and its time-dependent variant (TD-DFT) including implicit and explicit solvent effects. The calculations predict that cyanuric acid has an absorption maximum at ca. 160 nm (7.73 eV), with the lowest-energy absorption band extending to ca. 200 nm in an aqueous solution and exhibiting negligible absorption at longer wavelengths. Excitation of cyanuric acid at 160 nm or longer wavelengths leads to the population of S_5,6 singlet states, which have ππ* character and large oscillator strengths (0.8). The population reaching the S_5,6 states is expected to internally convert to the S_1,2 states in an ultrafast time scale. The S_1,2 states, which have nπ* character, are predicted to access a conical intersection with the ground state in a nearly barrierless fashion (ca. ≤ 0.13 eV), thus efficiently returning the population to the ground state. Furthermore, based on calculated spin–orbit coupling elements of ca. 6 to 8 cm⁻¹, the calculations predict that intersystem crossing to the triplet manifold should play a minor role in the electronic relaxation of cyanuric acid. We have also calculated the vertical ionization energy of cyanuric acid at 8.2 eV, which predicts that direct one-photon ionization of cyanuric acid should occur at ca. 150 nm. Collectively, the quantum-chemical calculations predict that cyanuric acid would have been highly photostable under the solar radiation conditions reaching the Earth’s surface during the prebiotic era in an aqueous solution. Of relevance to the chemical origin of life and RNA-first theories, these observations lend support to the idea that cyanuric acid could have accumulated in large quantities during the prebiotic era and thus strengthens its candidature as a relevant prebiotic nucleobase.

Impact of Charge-Resonance Excitations on CT-Mediated J-Type Aggregation in Singlet and Triplet Exciton States of Perylene Di-Imide Aggregates: A TDDFT Investigation

The modulation of intermolecular interactions upon aggregation induces changes in excited state properties of organic molecules that can be detrimental for some optoelectronic applications but can be exploited for others. The time-dependent density functional theory (TDDFT) is a cost-effective approach to determining the exciton states of molecular aggregates, and it has been shown to provide reliable results when coupled with the appropriate choice of the functional. Here we apply a general procedure to analyze the aggregates’ exciton states derived from TDDFT calculations in terms of diabatic states chosen to coincide with local (LE) and charge-transfer (CT) excitations within a restricted orbital space. We apply the approach to study energy profiles, interstate couplings, and the charge-transfer character of singlet and triplet exciton states of perylene di-imide aggregates (PDI). We focus on the intermolecular displacement along the longitudinal translation coordinate, which mimics different amounts of slip-stacking observed in PDI crystals. The analysis, in terms of symmetry-adapted Frenkel excitations (FE) and charge-resonance (CR) states and their interactions, discloses how the interchange of the H/J character for small longitudinal shifts, previously reported for singlet exciton states, also occurs for triplet excitons.

Excited States Computation of Models of Phenylalanine Protein Chains: TD-DFT and Composite CC2/TD-DFT Protocols

The present benchmark calculations testify to the validity of time-dependent density functional theory (TD-DFT) when exploring the low-lying excited states potential energy surfaces of models of phenylalanine protein chains. Among three functionals suitable for systems exhibiting charge-transfer excited states, LC-ωPBE, CAM-B3LYP, and ωB97X-D, which were tested on a reference peptide system, we selected the ωB97X-D functional, which gave the best results compared to the approximate coupled-cluster singles and doubles (CC2) method. A quantitative agreement for both the geometrical parameters and the vibrational frequencies was obtained for the lowest singlet excited state (a ππ* state) of the series of capped peptides. In contrast, only a qualitative agreement was met for the corresponding adiabatic zero-point vibrational energy (ZPVE)-corrected excitation energies. Two composite protocols combining CC2 and DFT/TD-DFT methods were then developed to improve these calculations. Both protocols substantially reduced the error compared to CC2 and experiment, and the best of both even led to results of CC2 quality at a lower cost, thus providing a reliable alternative to this method for very large systems.

Rotational analysis of naphthol-aromatic ring complexes stabilized by electrostatic and dispersion interactions

For complexes involving aromatic species, substitution effects can influence the preferred geometry. Using broadband rotational spectroscopy, we report the structures of three naphthol-aromatic ring complexes with different heteroatoms (furan and thiophene) and alkyl groups (2,5-dimethylfuran). The aim was to analyze the influence of the presence of heteroatoms or alkyl groups on the structure of the complex and the kind of intermolecular forces that control it. Face or edge arrangements can take place in these complexes via π-π or O-H⋯O/O-H⋯π interactions, respectively. All the experimentally observed complexes present O-H⋯O/O-H⋯π interactions with the hydroxyl group, with different structures and intermolecular interactions depending on the heteroatom present in the five-membered aromatic rings, yielding different symmetries in the experimental structure. Structures are experimentally identified through the use of planar moments of inertia. Further results from SAPT calculations show that dispersion and electrostatic interactions contribute similarly to the stabilization of all the studied complexes. These new spectroscopic results shed light on the influence of dispersion and hydrogen bonding in molecular aggregation of systems with substituted aromatic residues.

A Mountaineering Strategy to Excited States: Highly Accurate Energies and Benchmarks for Bicyclic Systems

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, furofuran, phthalazine, pyrrolopyrrole, quinoxaline, tetrathiafulvalene, and thienothiophene). In total, we provide 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), and 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.

Application of TD-DFT Theory to Studying Porphyrinoid-Based Photosensitizers for Photodynamic Therapy: A Review

An important focus for innovation in photodynamic therapy (PDT) is theoretical investigations. They employ mostly methods based on Time-Dependent Density Functional Theory (TD-DFT) to study the photochemical properties of photosensitizers. In the current article we review the existing state-of-the-art TD-DFT methods (and beyond) which are employed to study the properties of porphyrinoid-based systems. The review is organized in such a way that each paragraph is devoted to a separate aspect of the PDT mechanism, e.g., correct prediction of the absorption spectra, determination of the singlet-triplet intersystem crossing, and interaction with molecular oxygen. Aspects of the calculation schemes are discussed, such as the choice of the most suitable functional and inclusion of a solvent. Finally, quantitative structure-activity relationship (QSAR) methods used to explore the photochemistry of porphyrinoid-based systems are discussed.

Source of Bright Near-Infrared Luminescence in Gold Nanoclusters

Gold nanoclusters with near-infrared (NIR) photoluminescence (PL) have great potential as sensing and imaging materials in biomedical and bioimaging applications. In this work, Au₂₁(S-Adm)₁₅ and Au₃₈S₂(S-Adm)₂₀ are used to unravel the underlying mechanisms for the improved quantum yields (QY), large Stokes shifts, and long PL lifetimes in gold nanoclusters. Both nanoclusters show decent PL QY. In particular, the Au₃₈S₂(S-Adm)₂₀ nanocluster shows a bright NIR PL at 900 nm with QY up to 15% in normal solvents (such as toluene) at ambient conditions. The relatively lower QY for Au₂₁(S-Adm)₁₅ (4%) compared to that of Au₃₈S₂(S-Adm)₂₀ is attributed to the lowest-lying excited state being symmetry-disallowed, as evidenced by the pressure-dependent antispectral shift of the absorption spectra compared to PL, yet Au₂₁(S-Adm)₁₅ maintains some emissive properties due to a nearby symmetry-allowed excited state. Furthermore, our results show that suppression of nonradiative decay due to the surface "lock rings", which encircle the Au kernel and the surface "lock atoms" which bridge the fundamental Au kernel units (e.g., tetrahedra, icosahedra, etc.), is the key to obtaining high QYs in gold nanoclusters. The complicated excited-state processes and the small absorption coefficient of the band-edge transition lead to the large Stokes shifts and the long PL lifetimes that are widely observed in gold nanoclusters.

High-throughput virtual screening for organic electronics: a comparative study of alternative strategies

We present a review of the field of high-throughput virtual screening for organic electronics materials focusing on the sequence of methodological choices that determine each virtual screening protocol. These choices are present in all high-throughput virtual screenings and addressing them systematically will lead to optimised workflows and improve their applicability. We consider the range of properties that can be computed and illustrate how their accuracy can be determined depending on the quality and size of the experimental datasets. The approaches to generate candidates for virtual screening are also extremely varied and their relative strengths and weaknesses are discussed. The analysis of high-throughput virtual screening is almost never limited to the identification of top candidates and often new patterns and structure-property relations are the most interesting findings of such searches. The review reveals a very dynamic field constantly adapting to match an evolving landscape of applications, methodologies and datasets.

Simple Protocol for Capturing Both Linear-Response and State-Specific Effects in Excited-State Calculations with Continuum Solvation Models

We present an effective computational protocol (cLR²) to describe both solvatochromism and fluorosolvatochromism. This protocol, which couples the polarizable continuum model to time-dependent density functional theory, simultaneously accounts for both linear-response and state-specific solvation effects. A series of test cases, including solvatochromic and fluorosolvatochromic compounds and excited-state intramolecular proton transfers, are used to highlight that cLR² is especially beneficial for modeling bright excitations possessing a significant charge-transfer character, as well as cases in which an accurate balance between states of various polarities should be restored.

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies

The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.

Addressing the Frenkel and charge transfer character of exciton states with a model Hamiltonian based on dimer calculations: Application to large aggregates of perylene bisimide

To understand the influence of interchromophoric arrangements on photo-induced processes and optical properties of aggregates, it is fundamental to assess the contribution of local excitations [charge transfer (CT) and Frenkel (FE)] to exciton states. Here, we apply a general procedure to analyze the adiabatic exciton states derived from time-dependent density functional theory calculations, in terms of diabatic states chosen to coincide with local excitations within a restricted orbital space. In parallel, motivated by the need of cost-effective approaches to afford the study of larger aggregates, we propose to build a model Hamiltonian based on calculations carried out on dimers composing the aggregate. Both approaches are applied to study excitation energy profiles and CT character modulation induced by interchromophore rearrangements in perylene bisimide aggregates up to a tetramer. The dimer-based approach closely reproduces the results of full-aggregate calculations, and an analysis in terms of symmetry-adapted diabatic states discloses the effects of CT/FE interactions on the interchange of the H-/J-character for small longitudinal shifts of the chromophores.

Analytical Time-Dependent Long-Range Corrected Density Functional Tight Binding (TD-LC-DFTB) Gradients in DFTB+: Implementation and Benchmark for Excited-State Geometries and Transition Energies

The absorption and emission of light is a ubiquitous process in chemical and biological processes, making a theoretical description inevitable for understanding and predicting such properties. Although ab initio and DFT methods are capable of describing excited states with good accuracy in many cases, the investigation of dynamical processes and the need to sample the phase space in complex systems often requires methods with reduced computational costs but still sufficient accuracy. In the present work, we report the derivation and implementation of analytical nuclear gradients for time-dependent long-range corrected density functional tight binding (TD-LC-DFTB) in the DFTB+ program. The accuracy of the TD-LC-DFTB potential-energy surfaces is benchmarked for excited-state geometries and adiabatic as well as vertical transition energies. The benchmark set consists of more than 100 organic molecules taken as subsets from available benchmark sets. The reported method yields a mean deviation of 0.31 eV for adiabatic excitation energies with respect to CC2. In order to study more subtle effects, seminumerical second derivatives based on the analytical gradients are employed to simulate vibrationally resolved UV/vis spectra. This extensive test exhibits few problematic cases, which can be traced back to the parametrization of the repulsive potential.

Improved Description of Charge-Transfer Potential Energy Surfaces via Spin-Component-Scaled CC2 and ADC(2) Methods

The molecular level understanding of electronic transport properties depends on the reliable theoretical description of charge-transfer (CT)-type electronic states. In this paper, the performance of spin-component-scaled variants of the popular CC2 and ADC(2) methods is evaluated for CT states, following benchmark strategies of earlier studies that revealed a compromised accuracy of the unmodified models. In addition to statistics on the accuracy of vertical excitation energies at equilibrium and infinite separation of bimolecular complexes, potential energy surfaces of the ammonia-fluorine complex are also reported. The results show the capability of spin-component-scaled approaches to reduce the large errors of their regular counterparts to a significant extent, outperforming even the coupled-cluster single and double method in many cases. The cost-effective scaled-opposite-spin variants are found to provide a remarkably good agreement with the CCSDT-3 reference data, thereby being recommended methods of choice in the study of charge-transfer states.

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.

Density Functional Prediction of Quasiparticle, Excitation, and Resonance Energies of Molecules With a Global Scaling Correction Approach

Molecular quasiparticle and excitation energies determine essentially the spectral characteristics measured in various spectroscopic experiments. Accurate prediction of these energies has been rather challenging for ground-state density functional methods, because the commonly adopted density function approximations suffer from delocalization error. In this work, by presuming a quantitative correspondence between the quasiparticle energies and the generalized Kohn–Sham orbital energies, and employing a previously developed global scaling correction approach, we achieve substantially improved prediction of molecular quasiparticle and excitation energies. In addition, we also extend our previous study on temporary anions in resonant states, which are associated with negative molecular electron affinities. The proposed approach does not require any explicit self-consistent field calculation on the excited-state species, and is thus highly efficient and convenient for practical purposes.

A Potential Inhibition Process of Ricin Protein with the flavonoids Quercetin and Epigallocatechin Gallate. A Quantum-Chemical and Molecular Docking Study

Castor bean (Ricinus Communis) oil has been reported as one of the most important bio-based fuels; however, high amounts of toxic solid residue are generated in the production. This toxicity is due to several molecules, ricin protein being the most studied compound. The inhibition of the ricin protein is essential for eliminating its toxicity. The objective of this study is to predict the possible inhibition process via the interactions between the ricin protein and the flavonoids quercetin (Q) and epigallocatechin gallate (EGCG). The molecular structures of the complexes formed between the ricin protein and flavonoids were studied using quantum-chemical and molecular docking calculations to analyze the type of interaction, active site of the protein, binding energies, and different conformations in the inhibition process. Different methodologies were applied for the molecular structure determination; the best approximation was obtained with B3LYP/6-31G (d,p) theoretical methodology. Mappings of electrostatic potential (MEP) and frontier molecular orbitals were used for the identification of the probable sites of interaction, which were confirmed by molecular docking. The adjustment and alignment of flavonoid groups before and after the interaction, and charge transfer parameters, showed that Q and EGCG act as electron donors inside of the active site in ricin.

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.

Highly Polarized Coumarin Derivatives Revisited: Solvent-Controlled Competition Between Proton-Coupled Electron Transfer and Twisted Intramolecular Charge Transfer

Linking a polarized coumarin unit with an aromatic substituent via an amide bridge results in weak electronic coupling that affects the intramolecular electron-transfer (ET) process. As a result of this, interesting solvent-dependent photophysical properties can be observed. In polar solvents, electron transfer in coumarin derivatives of this type induces a mutual twist of the electron-donating and -accepting molecular units (TICT process) that facilitates radiationless decay processes (internal conversion). In the dyad with the strongest intramolecular hydrogen bond, the planar form is stabilized, such that twisting can only occur in highly polar solvents, whereas a fast proton-coupled electron-transfer (PCET process) occurs in nonpolar n-alkanes. The k_PCET rate constant decreases linearly with the energy of the fluorescence maximum in different solvents. This observation can be explained in terms of competition between electron- and proton-transfer from a highly polarized (ca. 15 D) and fluorescent locally excited (¹ LE) state to a much less polarized (ca. 4 D) charge-transfer (¹ CT) state, a unique occurrence. Photophysical measurements performed for a family of related coumarin dyads, together with results of quantum-chemical computations, give insight into the mechanism of the ET process, which is followed by either a TICT or a PCET process. Our results reveal that dielectric solvation of the excited state slows down the PCET process, even in nonpolar solvents.

Photoelimination of Nitrogen from Diazoalkanes: Involvement of Higher Excited Singlet States in the Carbene Formation

Although diazoalkanes are important carbene precursors in organic synthesis, a comprehensive mechanism of photochemical formation of carbenes from diazoalkanes has not been proposed. Synergies of experiments and computations demonstrate the involvement of higher excited singlet states in the photochemistry of diazoalkanes. In all investigated diazoalkanes, excitation to S₁ results in nonreactive internal conversion to S₀. On the contrary, excitation to higher-lying singlet states (S_n, n > 1) drives the reaction toward a different segment of the S₁/S₀ conical intersection seam and results in nitrogen elimination and formation of carbenes.

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.

Benchmark and performance of long-range corrected time-dependent density functional tight binding (LC-TD-DFTB) on rhodopsins and light-harvesting complexes

The chromophores of rhodopsins (Rh) and light-harvesting (LH) complexes still represent a major challenge for a quantum chemical description due to their size and complex electronic structure. Since gradient corrected and hybrid density functional approaches have been shown to fail for these systems, only range-separated functionals seem to be a promising alternative to the more time consuming post-Hartree-Fock approaches. For extended sampling of optical properties, however, even more approximate approaches are required. Recently, a long-range corrected (LC) functional has been implemented into the efficient density functional tight binding (DFTB) method, allowing to sample the excited states properties of chromophores embedded into proteins using quantum mechanical/molecular mechanical (QM/MM) with the time-dependent (TD) DFTB approach. In the present study, we assess the accuracy of LC-TD-DFT and LC-TD-DFTB for rhodopsins (bacteriorhodopsin (bR) and pharaonis phoborhodopsin (ppR)) and LH complexes (light-harvesting complex II (LH2) and Fenna-Matthews-Olson (FMO) complex). This benchmark study shows the improved description of the color tuning parameters compared to standard DFT functionals. In general, LC-TD-DFTB can exhibit a similar performance as the corresponding LC functionals, allowing a reliable description of excited states properties at significantly reduced cost. The two chromophores investigated here pose complementary challenges: while huge sensitivity to external field perturbation (color tuning) and charge transfer excitations are characteristic for the retinal chromophore, the multi-chromophoric character of the LH complexes emphasizes a correct description of inter-chromophore couplings, giving less importance to color tuning. None of the investigated functionals masters both systems simultaneously with satisfactory accuracy. LC-TD-DFTB, at the current stage, although showing a systematic improvement compared to TD-DFTB cannot be recommended for studying color tuning in retinal proteins, similar to some of the LC-DFT functionals, because the response to external fields is still too weak. For sampling of LH-spectra, however, LC-TD-DFTB is a viable tool, allowing to efficiently sample absorption energies, as shown for three different LH complexes. As the calculations indicate, geometry optimization may overestimate the importance of local minima, which may be averaged over when using trajectories. Fast quantum chemical approaches therefore may allow for a direct sampling of spectra in the near future.

Benchmarking the Performance of Time-Dependent Density Functional Theory Methods on Biochromophores

Quantum chemical calculations are important for elucidating light-capturing mechanisms in photobiological systems. The time-dependent density functional theory (TDDFT) has become a popular methodology because of its balance between accuracy and computational scaling, despite its problems in describing, for example, charge transfer states. As a step toward systematically understanding the performance of TDDFT calculations on biomolecular systems, we study here 17 commonly used density functionals, including seven long-range separated functionals, and compare the obtained results with excitation energies calculated at the approximate second order coupled-cluster theory level (CC2). The benchmarking set includes the first five singlet excited states of 11 chemical analogues of biochromophores from the green fluorescent protein, rhodopsin/bacteriorhodopsin (Rh/bR), and the photoactive yellow protein. We find that commonly used pure density functionals such as BP86, PBE, M11-L, and hybrid functionals with 20-25% of Hartree-Fock (HF) exchange (B3LYP, PBE0) have a tendency to consistently underestimate vertical excitation energies (VEEs) relative to the CC2 values, whereas hybrid density functionals with around 50% HF exchange such as BHLYP, PBE50, and M06-2X and long-range corrected functionals such as CAM-B3LYP, ωPBE, ωPBEh, ωB97X, ωB97XD, BNL, and M11 overestimate the VEEs. We observe that calculations using the CAM-B3LYP and ωPBEh functionals with 65% and 100% long-range HF exchange, respectively, lead to an overestimation of the VEEs by 0.2-0.3 eV for the benchmarking set. To reduce the systematic error, we introduce here two new empirical functionals, CAMh-B3LYP and ωhPBE0, for which we adjusted the long-range HF exchange to 50%. The introduced parameterization reduces the mean signed average (MSA) deviation to 0.07 eV and the root mean square (rms) deviation to 0.17 eV as compared to the CC2 values. In the present study, TDDFT calculations using the aug-def2-TZVP basis sets, the best performing functionals relative to CC2 are ωhPBE0 (rms = 0.17, MSA = 0.06 eV); CAMh-B3LYP (rms = 0.16, MSA = 0.07 eV); and PBE0 (rms = 0.23, MSA = -0.14 eV). For the popular range-separated CAM-B3LYP functional, we obtain an rms value of 0.31 eV and an MSA value of 0.25 eV, which can be compared with the rms and MSA values of 0.37 and -0.31 eV, respectively, as obtained at the B3LYP level.