Benchmarking TD-DFT and Wave Function Methods for Oscillator Strengths and Excited-State Dipole Moments

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.

A DFT/MRCI Hamiltonian parameterized using only ab initio data: I. valence excited states

A new combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian parameterized solely using the benchmark ab initio vertical excitation energies obtained from the QUEST databases is presented. This new formulation differs from all previous versions of the method in that the choice of the underlying exchange-correlation (XC) functional employed to construct the one-particle (orbital) basis is considered, and a new XC functional, QTP17, is chosen for its ability to generate a balanced description of core and valence vertical excitation energies. The ability of the new DFT/MRCI Hamiltonian, termed QE8, to furnish accurate excitation energies is confirmed using benchmark quantum chemistry computations, and a mean absolute error of 0.16 eV is determined for the wide range of electronic excitations included in the validation dataset. In particular, the QE8 Hamiltonian dramatically improves the performance of DFT/MRCI for doubly excited states. The performance of fast approximate DFT/MRCI methods, p-DFT/MRCI and DFT/MRCI(2), is also evaluated using the QE8 Hamiltonian, and they are found to yield excitation energies in quantitative agreement with the parent DFT/MRCI method, with the two methods exhibiting a mean difference of 0.01 eV with respect to DFT/MRCI over the entire benchmark set.

Toward Reliable Dipole Moments without Single Excitations: The Role of Orbital Rotations and Dynamical Correlation

The dipole moment is a crucial molecular property linked to a molecular system's bond polarity and overall electronic structure. To that end, the electronic dipole moment, which results from the electron density of a system, is often used to assess the accuracy and reliability of new electronic structure methods. This work analyses electronic dipole moments computed with the pair coupled cluster doubles (pCCD) ansätze and its linearized coupled cluster (pCCD-LCC) corrections using the canonical Hartree-Fock and pCCD-optimized (localized) orbital bases. The accuracy of pCCD-based dipole moments is assessed against experimental and CCSD(T) reference values using relaxed and unrelaxed density matrices and different basis set sizes. Our test set comprises molecules of various bonding patterns and electronic structures, exposing pCCD-based methods to a wide range of electron correlation effects. Additionally, we investigate the performance of pCCD-in-DFT dipole moments of some model complexes. Finally, our work indicates the importance of orbital relaxation in the pCCD model and shows the limitations of the linearized couple cluster corrections in predicting electronic dipole moments of multiple-bonded systems. Most importantly, pCCD with a linearized CCD correction can reproduce the dipole moment surfaces in singly bonded molecules, which are comparable to the multireference ones.

State-Specific Coupled-Cluster Methods for Excited States

We reexamine ΔCCSD, a state-specific coupled-cluster (CC) with single and double excitations (CCSD) approach that targets excited states through the utilization of non-Aufbau determinants. This methodology is particularly efficient when dealing with doubly excited states, a domain in which the standard equation-of-motion CCSD (EOM-CCSD) formalism falls short. Our goal here to evaluate the effectiveness of ΔCCSD when applied to other types of excited states, comparing its consistency and accuracy with EOM-CCSD. To this end, we report a benchmark on excitation energies computed with the ΔCCSD and EOM-CCSD methods for a set of molecular excited-state energies that encompasses not only doubly excited states but also doublet-doublet transitions and (singlet and triplet) singly excited states of closed-shell systems. In the latter case, we rely on a minimalist version of multireference CC known as the two-determinant CCSD method to compute the excited states. Our data set, consisting of 276 excited states stemming from the quest database [Véril et al., WIREs Comput. Mol. Sci. 2021, 11, e1517], provides a significant base to draw general conclusions concerning the accuracy of ΔCCSD. Except for the doubly excited states, we found that ΔCCSD underperforms EOM-CCSD. For doublet-doublet transitions, the difference between the mean absolute errors (MAEs) of the two methodologies (of 0.10 and 0.07 eV) is less pronounced than that obtained for singly excited states of closed-shell systems (MAEs of 0.15 and 0.08 eV). This discrepancy is largely attributed to a greater number of excited states in the latter set exhibiting multiconfigurational characters, which are more challenging for ΔCCSD. We also found typically small improvements by employing state-specific optimized orbitals.

Beyond Explored Functionals: A Computational Journey of Two-Photon Absorption

We present a thorough investigation into the efficacy of 19 density functional theory (DFT) functionals, relative to RI-CC2 results, for computing two-photon absorption (2PA) cross sections (σ2PA) and key dipole moments (|μ00|, |μ11|, |Δμ|, |μ01|) for a series of coumarin dyes in the gas-phase. The functionals include different categories, including local density approximation (LDA), generalized gradient approximation (GGA), hybrid-GGA (H-GGA), range-separated hybrid-GGA (RSH-GGA), meta-GGA (M-GGA), and hybrid M-GGA (HM-GGA), with 14 of them being subjected to analysis for the first time with respect to predicting σ2PA values. Analysis reveals that functionals integrating both short-range (SR) and long-range (LR) corrections, particularly those within the RSH-GGA and HM-GGA classes, outperform the others. Furthermore, the range-separation approach was found more impactful compared to the varying percentages of Hartree-Fock exchange (HF Ex) within different functionals. The functionals traditionally recommended for 2PA do not appear among the top 9 in our study, which is particularly interesting, as these top-performing functionals have not been previously investigated in this context. This list is dominated by M11, QTP variants, ωB97X, ωB97X-V, and M06-2X, surpassing the performance of other functionals, including the commonly used CAM-B3LYP.

Modeling of High-Harmonic Generation in the C60 Fullerene Using Ab Initio, DFT-Based, and Semiempirical Methods

We report calculations of the high-harmonic generation spectra of the C60 fullerene molecule carried out by employing a diverse set of real-time time-dependent quantum chemical methods. All methodologies involve expanding the propagated electronic wave function in bases consisting of the ground and singly excited time-independent eigenstates obtained through the solution of the corresponding linear-response equations. We identify the correlation and exchange effect in the spectra by comparing the results from methods relying on the Hartree-Fock reference determinant with those obtained using approaches based on the density functional theory with different exchange-correlation functionals. The effect of the full random-phase approximation treatment of the excited electronic states is also analyzed and compared with the configuration interaction singles and the Tamm-Dancoff approximation. We also showcase the fact that the real-time extension of the semiempirical method INDO/S can be effectively applied for an approximate description of laser-driven dynamics in large systems.

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.

Azaindolizine proton cranes attached to 7-hydroxyquinoline and 3-hydroxypyridine: a comparative theoretical study

Theoretical design of several proton cranes, based on 7-hydroxyquinoline and 3-hydroxypyridine as proton-transfer frames, has been attempted using ground and excited-state density functional theory (DFT) calculations in various environments. Imidazo[1,2-a]pyridine, pyrazolo[1,5-a]pyridine and benzimidazole were considered as proton crane units. The proton crane action requires the existence of a single enol-like form in the ground state, which under excitation goes to the end keto-like one through a series of consecutive excited-state intramolecular proton transfers (ESIPT) and twisting steps with the participation of a crane unit, resulting in a long-range intramolecular proton transfer. The results suggest that 3-hydroxypyridine is not suitable for a proton-transfer frame and 8-(imidazo[1,2-a]pyridin-2-yl)quinolin-7-ol and 8-(pyrazolo[1,5-a]pyridin-2-yl)quinolin-7-ol behave as non-conjugated proton cranes, instead of tautomeric re-arrangement in the latter, which was thought to be possible.

Theoretical Understanding of the Long-Range Proton Transfer Mechanism in 7-Hydroxy Quinoline-Based Molecular Switches: Varma's Proton Crane and Its Analogues

Herein, the detailed mechanism of intramolecular proton transfer in molecular switches, constructed from 7-hydroxy quinoline substituted in the eight-position C-C single axle, connected to three different proton cranes (morpholine, piperidine, and 1,3,5-dioxazine), was investigated by means of theoretical chemistry. The theoretical interpretation of the rotational mechanism and its stable structures were proposed for the well-known Varma's proton crane, based on morpholine molecule. The reliability of the theoretical simulations was confirmed by the available literature data from time-dependent IR measurements.

Activating two-dimensional semiconductors for photocatalysis: a cross-dimensional strategy

The emerging two-dimensional (2D) semiconductors substantially extend materials bases for versatile applications such as semiconductor photocatalysis demanding semiconductive matrices and large surface areas. The dimensionality, while endowing 2D semiconductors the unique properties to host photocatalytic functionality of pollutant removal and hydrogen evolution, hurdles the activation paths to form heterogenous photocatalysts where the photochemical processes are normally superior over these on the mono-compositional counterparts. In this perspective, we present a cross-dimensional strategy to employ thenD (n= 0-2) clusters or nanomaterials as activation partners to boost the photocatalytic activities of the 2D semiconductors. The formation principles of heterogenous photocatalysts are illustrated specifically for the 2D matrices, followed by selection criteria of them among the vast 2D database. The computer investigations are illustrated in the density functional theory route and machine learning benefitted from the vast samples in the 2D library. Synthetic realizations and characterizations of the 2D heterogenous systems are introduced with an emphasis on chemical methods and advanced techniques to understand materials and mechanistic studies. The perspective outlooks cross-dimensional activation strategies of the 2D materials for other applications such as CO2removal, and materials matrices in other dimensions which may inspire incoming research within these fields.

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.

Effect of Benzothiadiazole-Based π-Spacers on Fine-Tuning of Optoelectronic Properties of Oligothiophene-Core Donor Materials for Efficient Organic Solar Cells: A DFT Study

In this work, five novel A-π-D-π-A type molecules D1-D5 were designed by adding unusual benzothiadiazole derivatives as π-spacer blocks to the efficient reference molecule DRCN5T for application as donor materials in organic solar cells (OSCs). Based on a density functional theory approach, a comprehensive theoretical study was performed with different functionals (B3LYP, B3LYP-GD3, B3LYP-GD3BJ, CAM-B3LYP, M06, M062X, and wB97XD) and with different solvent types (PCM and SMD) at the extended basis set 6-311+g(d,p) to evaluate the structural, optoelectronic, and intramolecular charge transfer properties of these molecules. The B3LYP-GD3BJ hybrid functional was used to optimize the studied molecules in CHCl3 solution with the SMD model solvent as it provided the best results compared to experimental data. Transition density matrix maps were simulated to examine the hole-electron localization and the electronic excitation processes in the excited state, and photovoltaic parameters including open-circuit photovoltage and fill factor were investigated to predict the efficiency of these materials. All the designed materials showed promising optoelectronic and photovoltaic characteristics, and for most of them, a red shift. Out of the proposed molecules, [1,2,5]thiadiazolo[3,4-d]pyridazine was selected as a promising π-spacer block to evaluate its interaction with PC61BM in a composite to understand the charge transfer between the donor and acceptor subparts. Overall, this study showed that adding π-spacer building blocks to the molecular structure is undoubtedly a potential strategy to further enhance the performance of donor materials for OSC applications.

Oscillator strengths and excited-state couplings for double excitations in time-dependent density functional theory

Although useful to extract excitation energies of states of double-excitation character in time-dependent density functional theory that are missing in the adiabatic approximation, the frequency-dependent kernel derived earlier [Maitra et al., J. Chem. Phys. 120, 5932 (2004)] was not designed to yield oscillator strengths. These are required to fully determine linear absorption spectra, and they also impact excited-to-excited-state couplings that appear in dynamics simulations and other quadratic response properties. Here, we derive a modified non-adiabatic kernel that yields both accurate excitation energies and oscillator strengths for these states. We demonstrate its performance on a model two-electron system, the Be atom, and on excited-state transition dipoles in the LiH molecule at stretched bond-lengths, in all cases producing significant improvements over the traditional approximations.

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.

Excess and excited-state dipole moments of real-life dyes: a comparison between wave-function, BSE/GW, and TD-DFT values

In this work, we assess the accuracy of the Bethe-Salpeter equation (BSE) many-body Green's function formalism, adopting the eigenvalue-self-consistent evGW exchange-correlation kernel, for the calculation of the excited-state (μES) and excess dipole moments (Δμ), the latter ones being the changes of dipole amplitude between the ground and excited states (ES), in organic dyes. We compare the results obtained with wave-function methods [ADC(2), CC2, and CCSD], time-dependent density functional theory (TD-DFT), and BSE/evGW levels of theory. First, we compute the evolution of the dipole moments of the two lowest singlet excited states of 4-(dimethylamino)benzonitrile (DMABN) upon twisting of the amino group. Next, we use a set of 25 dyes having ES characters ranging from locally excited to charge transfer to determine both μES and Δμ. For DMABN our results show that BSE/evGW provides Δμ values closer to the CCSD reference and more consistent trends than TD-DFT. Moreover, a statistical analysis of both Δμ and μES for the set of 25 dyes shows that the BSE/evGW accuracy is comparable or sometimes slightly better than that of TD-M06-2X and TD-CAM-B3LYP, BSE/evGW outperforming TD-DFT in challenging cases (zwitterionic and cyanine transitions). Finally, the starting point dependency of BSE/evGW seems to be larger for Δμ, ES dipoles, and oscillator strengths than for transition energies.

Single-atom infrared emission in doped silicon nanocrystals

Silicon luminescence, due to silicon being abundant, non-toxic and harmless, is a topic of pivotal importance in optoelectronics and biological imaging. However, a major challenge in developing high-efficiency silicon light sources is the relatively weak allowable transitions. This study focuses on single atom-doped silicon nanocrystals (Si NCs) and theoretically investigates the emission behavior of single atoms within a tetrahedral coordination field. Doping a single atom in Si NCs can result in a ∼102 times improvement at least in the squared transition dipole moment (TDM2), and induce a spectral shift towards near- and mid-infrared wavelengths. These findings offer a strong foundation for designing Si NCs for on-chip optical communication and single photon emitters.

Origin-Independent Dynamic Polarizability Density from Coupled Cluster Response Theory

The calculation of the origin-independent density of the dynamic electric dipole polarizability, previously presented for uncorrelated and density functional theory (DFT)-based methods, has been developed and implemented at the coupled cluster singles and doubles (CCSD) level of theory. A pointwise analysis of polarizability densities calculated for a number of molecules at Hartree-Fock (HF) and CCSD clearly shows that the electron correlation effect is much larger than one would argue considering the integrated dipole electric polarizability alone. Large error compensations occur during the integration process, which hide fairly large deviations mainly located in the internuclear regions. The same is observed when calculated CCSD and B3LYP polarizability densities are compared, with the remarkable feature that positive/negative deviations between CCSD and HF reverse sign, becoming negative/positive when comparing CCSD to B3LYP.

Full theoretical protocol for the design of metal-free organic electron donor-spacer-acceptor systems

A user-friendly (time-dependent) density functional theory based algorithm is proposed to design new donor-spacer-acceptor systems for electron transfer reactions. This algorithm is focused on metal-free organic compounds, most of which contain aromatic or alkene moieties. The oxidation and reduction potentials are calculated, together with the excited-state energy difference including the zero-point energy and the structural properties required to calculate an electron transfer Gibbs free energy change. The proposed algorithm has been tested on well-known systems, while two new compounds are suggested for photoinduced intramolecular electron transfer reactions using this scheme. The methodology here presented is intended to be a tool for synthetic physical-chemists, allowing them to evaluate the properties of hypothetical systems before the synthesis, enabling the study of limitless combinations of donor-spacer-acceptor arrangements.

A Theoretical Study on the Underlying Factors of the Difference in Performance of Organic Solar Cells Based on ITIC and Its Isomers

Recently, non-fullerene-based organic solar cells (OSCs) have made great breakthroughs, and small structural differences can have dramatic impacts on the power conversion efficiency (PCE). We take ITIC and its isomers as examples to study their effects on the performance of OSCs. ITIC and NFBDT only differed in the side chain position, and they were used as models with the same donor molecule, PBDB-T, to investigate the main reasons for the difference in their performance in terms of theoretical methods. In this work, a detailed comparative analysis of the electronic structure, absorption spectra, open circuit voltage and interfacial parameters of the ITIC and NFBDT systems was performed mainly by combining the density functional theory/time-dependent density functional theory and molecular dynamics simulations. The results showed that the lowest excited state of the ITIC molecule possessed a larger ∆q and more hybrid FE/CT states, and PBDB-T/ITIC had more charge separation paths as well as a larger kCS and smaller kCR. The reason for the performance difference between PBDB-T/ITIC and PBDB-T/NFBDT was elucidated, suggesting that ITIC is a superior acceptor based on a slight modulation of the side chain and providing a guiding direction for the design of superior-performing small molecule acceptor materials.

The Best Models of Bodipy's Electronic Excited State: Comparing Predictions from Various DFT Functionals with Measurements from Femtosecond Stimulated Raman Spectroscopy

Density functional theory (DFT) and time-dependent DFT (TD-DFT) are pivotal approaches for modeling electronically excited states of molecules. However, choosing a DFT exchange-correlation functional (XCF) among the myriad of alternatives is an overwhelming task that can affect the interpretation of results and lead to erroneous conclusions. The performance of these XCFs to describe the excited-state properties is often addressed by comparing them with high-level wave function methods or experimentally available vertical excitation energies; however, this is a limited analysis that relies on evaluation of a single point in the excited-state potential energy surface (PES). Different strategies have been proposed but are limited by the difficulty of experimentally accessing the electronic excited-state properties. In this work, we have tested the performance of 12 different XCFs and TD-DFT to describe the excited-state potential energy surface of Bodipy (2,6-diethyl-1,3,5,7-tetramethyl-8-phenyldipyrromethene difluoroborate). We compare those results with resonance Raman spectra collected by using femtosecond stimulated Raman spectroscopy (FSRS). By simultaneously fitting the absorption spectrum, fluorescence spectrum, and all of the resonance Raman excitation profiles within the independent mode displaced harmonic oscillator (IMDHO) formalism, we can describe the PES at the Franck-Condon (FC) region and determine the solvent and intramolecular reorganization energy after relaxation. This allows a direct comparison of the TD-DFT output with experimental observables. Our analysis reveals that using vertical absorption energies might not be a good criterion to determine the best XCF for a given molecular system and that FSRS opens up a new way to benchmark the excited-state performance of XCFs of fluorescent dyes.

Photo-Induced Active Lewis Acid-Base Pairs in a Metal-Organic Framework for H2 Activation

The establishment of active sites as the frustrated Lewis pair (FLP) has recently attracted much attention ranging from homogeneous to heterogeneous systems in the field of catalysis. Their unquenched reactivity of Lewis acid and base pairs in close proximity that are unable to form stable adducts has been shown to activate small molecules such as dihydrogen heterolytically. Herein, we show that grafted Ru metal-organic framework-based catalysts prepared via N-containing linkers are rather catalytically inactive for H2 activation despite the application of elevated temperatures. However, upon light illumination, charge polarization of the anchored Ru bipyridine complex can form a transient Lewis acid-base pair, Ru+-N- via metal-to-ligand charge transfer, as confirmed by time-dependent density functional theory (TDDFT) calculations to carry out effective H2-D2 exchange. FTIR and 2-D NMR endorse the formation of such reactive intermediate(s) upon light irradiation.

Tweaking the conjugation effects on a pair of new triazene compounds by targeted deprotonation: a spectroscopic and theoretical overview

CONTEXT

Triazene compounds (-NNN(H)-) exhibit versatility in biological, physical, and chemical applications. In their anionic form (-NNN-)(-), they can act as coordinating sites for metals, forming metallic complexes. In this study, two new isomeric triazene compounds with meta- and para-substituents in their neutral and anionic forms were investigated. A combination of detailed experimental spectroscopic characterization and computational chemistry analyses were employed. The new compounds, 1-(2-benzamide)-3-(3-nitrophenyl) triazene (m-TZN) and 1-(2-benzamide)-3-(4-nitrophenyl) triazene (p-TZN), were compared to 1,3-diphenyltriazene (dph-TZN) to understand the effects of functionalization and targeted triazene deprotonation. The anionic forms are stable, and our investigation suggests that these new compounds are suitable tridentate ligands that can act as chelating agents for metallic cations in stable complexes, similar to those found in vitamin B12.

METHODS

The absorption, vibrational, and electronic properties of the newly synthesized triazene compounds were extensively characterized using FT-IR/FT-Raman and UV-Vis spectroscopy. Their distinct molecular properties, intramolecular hydrogen bond effects, stability, and electronic transitions were investigated using the ORCA software. These analyses involved DFT and TD-DFT calculations at the ωB97X-D3/Def2-TZVP level of theory with THF CPCM implicit solvation to determine the molecular topology and electronic structure. The advanced STEOM-DLPNO-CCSD method for excited states was employed, enabling an in-depth analysis of ground and excited-state chemistry, accounting for precise electronic correlation and solvation effects. Explicit THF solvation was tested on the full TD-DFT ωB97X-D3/Def2-TZVP level and using ONIOM on the STEOM calculation. Reactivity was studied using Fukui functions, and action as chelating agents was investigated using GFN-xTB2 and DFT.

Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation

Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.

Multistate Energy Decomposition Analysis of Molecular Excited States

A multistate energy decomposition analysis (MS-EDA) method is described to dissect the energy components in molecular complexes in excited states. In MS-EDA, the total binding energy of an excimer or an exciplex is partitioned into a ground-state term, called local interaction energy, and excited-state contributions that include exciton excitation energy, superexchange stabilization, and orbital and configuration-state delocalization. An important feature of MS-EDA is that key intermediate states associated with different energy terms can be variationally optimized, providing quantitative insights into widely used physical concepts such as exciton delocalization and superexchange charge-transfer effects in excited states. By introducing structure-weighted adiabatic excitation energy as the minimum photoexcitation energy needed to produce an excited-state complex, the binding energy of an exciplex and excimer can be defined. On the basis of the nature of intermolecular forces through MS-EDA analysis, it was found that molecular complexes in the excited states can be classified into three main categories, including (1) encounter excited-state complex, (2) charge-transfer exciplex, and (3) intimate excimer or exciplex. The illustrative examples in this Perspective highlight the interplay of local excitation polarization, exciton resonance, and superexchange effects in molecular excited states. It is hoped that MS-EDA can be a useful tool for understanding photochemical and photobiological processes.

Organic compounds for solid state luminescence enhancement/aggregation induced emission: a theoretical perspective

Organic luminophores displaying one or more forms of luminescence enhancement in solid state are extremely promising for the development and performance optimization of functional materials essential to many modern key technologies. Yet, the effort to harness their huge potential is riddled with hurdles that ultimately come down to a limited understanding of the interactions that result in the diverse molecular environments responsible for the macroscopic response. In this context, the benefits of a theoretical framework able to provide mechanistic explanations to observations, supported by quantitative predictions of the phenomenon, are rather apparent. In this perspective, we review some of the established facts and recent developments about the current theoretical understanding of solid-state luminescence enhancement (SLE) with an accent on aggregation-induced emission (AIE). A description of the macroscopic phenomenon and the questions it raises is accompanied by a discussion of the approaches and quantum chemistry methods that are more apt to model these molecular systems with the inclusion of an accurate yet efficient simulation of the local environment. A sketch of a general framework, building from the current available knowledge, is then attempted via the analysis of a few varied SLE/AIE molecular systems from literature. A number of fundamental elements are identified offering the basis for outlining design rules for molecular architectures exhibiting SLE that involve specific structural features with the double role of modulating the optical response of the luminophores and defining the environment they experience in solid state.

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.

Automated Active Space Selection with Dipole Moments

Multireference calculations can provide accurate information of systems with strong correlation, which have increasing importance in the development of new molecules and materials. However, selecting a suitable active space for multireference calculations is nontrivial, and the selection of an unsuitable active space can sometimes lead to results that are not physically meaningful. Active space selection often requires significant human input, and the selection that leads to reasonable results often goes beyond chemical intuition. In this work, we have developed and evaluated two protocols for automated selection of the active space for multireference calculations based on a simple physical observable, the dipole moment, for molecules with nonzero ground-state dipole moments. One protocol is based on the ground-state dipole moment, and the other is based on the excited-state dipole moments. To evaluate the protocols, we constructed a dataset of 1275 active spaces from 25 molecules, each with 51 active space sizes considered, and have mapped out the relationship between the active space, dipole moments, and vertical excitation energies. We have demonstrated that, within this dataset, our protocols allow one to choose among a number of accessible active spaces one that is likely to give reasonable vertical excitation energies, especially for the first three excitations, with no parameters manually decided by the user. We show that, with large active spaces removed from consideration, the accuracy is similar and the time-to-solution can be reduced by more than 10 fold. We also show that the protocols can be applied to potential energy surface scans and determining the spin states of transition metal oxides.

Simulation of the VUV Absorption Spectra of Oxygenates and Hydrocarbons: A Joint Theoretical-Experimental Study

Vacuum UV absorption spectroscopy is regularly used to provide unambiguous identification of a target species, insight into the electronic structure of molecules, and quantitative species concentrations. As molecules of interest have become more complex, theoretical spectra have been used in tandem with laboratory spectroscopic analysis or as a replacement when experimental data is unavailable. However, it is difficult to determine which theoretical methodologies can best simulate experiment. This study examined the performance of EOM-CCSD and 10 TD-DFT functionals (B3LYP, BH&HLYP, BMK, CAM-B3LYP, HSE, M06-2X, M11, PBE0, ωB97X-D, and X3LYP) to produce reliable vacuum UV absorption spectra for 19 small oxygenates and hydrocarbons using vertical excitation energies. The simulated spectra were analyzed against experiment using both a qualitative analysis and quantitative metrics, including cosine similarity, relative integral change, mean signed error, and mean absolute error. Based on our ranking system, it was determined that M06-2X was consistently the top performing TD-DFT method with BMK, CAM-B3LYP, and ωB97X-D also producing reliable spectra for these small combustion species.

Scoring molecular wires subject to an ultrafast laser pulse for molecular electronic devices

A nonionizing ultrafast laser pulse of 20-fs duration with a peak amplitude electric-field ±E = 200 × 10^-4 a.u. was simulated. It was applied to the ethene molecule to consider its effect on the electron dynamics, both during the application of the laser pulse and for up to 100 fs after the pulse was switched off. Four laser pulse frequencies ω = 0.2692, 0.2808, 0.2830, and 0.2900 a.u. were chosen to correspond to excitation energies mid-way between the (S₁ ,S₂ ), (S₂ ,S₃ ), (S₃ ,S₄ ) and (S₄ ,S₅ ) electronic states, respectively. Scalar quantum theory of atoms in molecules (QTAIM) was used to quantify the shifts of the C1C2 bond critical points (BCPs). Depending on the frequencies ω selected, the C1C2 BCP shifts were up to 5.8 times higher after the pulse was switched off compared with a static E-field with the same magnitude. Next generation QTAIM (NG-QTAIM) was used to visualize and quantify the directional chemical character. In particular, polarization effects and bond strengths, in the form of bond-rigidity vs. bond-flexibility, were found, for some laser pulse frequencies, to increase after the laser pulse was switched off. Our analysis demonstrates that NG-QTAIM, in partnership with ultrafast laser irradiation, is useful as a tool in the emerging field of ultrafast electron dynamics, which will be essential for the design, and control of molecular electronic devices.

Exploring Bethe-Salpeter Excited-State Dipoles: The Challenging Case of Increasingly Long Push-Pull Oligomers

The change of molecular dipole moment induced by photon absorption is key to interpret the measured optical spectra. Except for compact molecules, time-dependent density functional theory (TD-DFT) remains the only theory allowing to quickly predict excited-state dipoles (μ^ES), albeit with a strong dependency on the selected exchange-correlation functional. This Letter presents the first assessment of the performances of the many-body Green's function Bethe-Salpeter equation (BSE) formalism for the evaluation of the μ^ES. We explore increasingly long push-pull oligomers as they present an excited-state nature evolving with system size. This work shows that BSE's μ^ES do present the same evolution with oligomeric length as their CC2 and CCSD counterparts, with a dependency on the starting exchange-correlation functional that is strongly decreased as compared to TD-DFT. This Letter demonstrates that BSE is a valuable alternative to TD-DFT for properties related to the excited-state density and not only for transition energies and oscillator strengths.

A direct observation of up-converted room-temperature phosphorescence in an anti-Kasha dopant-matrix system

It is common sense that emission maxima of phosphorescence spectra (λ_P) are longer than those of fluorescence spectra (λ_F). Here we report a serendipitous finding of up-converted room-temperature phosphorescence (RTP) with λ_P < λ_F and phosphorescence lifetime > 0.1 s upon doping benzophenone-containing difluoroboron β-diketonate (BPBF₂) into phenyl benzoate matrices. The up-converted RTP is originated from BPBF₂'s T_n (n ≥ 2) states which show typical ³n-π* characters from benzophenone moieties. Detailed studies reveal that, upon intersystem crossing from BPBF₂'s S₁ states of charge transfer characters, the resultant T₁ and T_n states build T₁-to-T_n equilibrium. Because of their ³n-π* characters, the T_n states possess large phosphorescence rates that can strongly compete RTP(T₁) to directly emit RTP(T_n) which violates Kasha's rule. The direct observation of up-converted RTP provides deep understanding of triplet excited state dynamics and opens an intriguing pathway to devise visible-light-excitable deep-blue afterglow emitters, as well as stimuli-responsive afterglow materials.

Cost-Effective Simulations of Vibrationally-Resolved Absorption Spectra of Fluorophores with Machine-Learning-Based Inhomogeneous Broadening

The results of electronic and vibrational structure simulations are an invaluable support for interpreting experimental absorption/emission spectra, which stimulates the development of reliable and cost-effective computational protocols. In this work, we contribute to these efforts and propose an efficient first-principle protocol for simulating vibrationally-resolved absorption spectra, including nonempirical estimations of the inhomogeneous broadening. To this end, we analyze three key aspects: (i) a metric-based selection of density functional approximation (DFA) so to benefit from the computational efficiency of time-dependent density function theory (TD-DFT) while safeguarding the accuracy of the vibrationally-resolved spectra, (ii) an assessment of two vibrational structure schemes (vertical gradient and adiabatic Hessian) to compute the Franck-Condon factors, and (iii) the use of machine learning to speed up nonempirical estimations of the inhomogeneous broadening. In more detail, we predict the absorption band shapes for a set of 20 medium-sized fluorescent dyes, focusing on the bright ππ^★ S₀ → S₁ transition and using experimental results as references. We demonstrate that, for the studied 20-dye set which includes structures with large structural variability, the preselection of DFAs based on an easily accessible metric ensures accurate band shapes with respect to the reference approach and that range-separated functionals show the best performance when combined with the vertical gradient model. As far as band widths are concerned, we propose a new machine-learning-based approach for determining the inhomogeneous broadening induced by the solvent microenvironment. This approach is shown to be very robust offering inhomogeneous broadenings with errors as small as 2 cm^-1 with respect to genuine electronic-structure calculations, with a total CPU time reduced by 98%.

Benchmarking DFT functionals for photophysics of pyranoflavylium cations

An intense absorption, phosphorescence, a long triplet excited state lifetime and singlet oxygen generation capabilities are characteristics of pyranoflavylium cations, analogues to pyranoanthocyanidins originated in the maturation process of red wine. Such properties make these compounds potential photosensitizers to be applied in photodynamic therapy. In this context, the photophysical processes underlying that treatment critically depend on the electronic structure of the pyranoflavylium molecules. When employing density functional theory to describe the electronic structure of molecules, the choice of the most suitable functional is not trivial, and benchmark studies are needed to orient practitioners in the field. In this work, a benchmark of seven of the most commonly used density functionals in addressing the photophysical properties of a set of eight pyranoflavylium cations is reported. Ground and excited state geometries, molecular orbitals, and absorption, fluorescence and phosphorescence transition energies were calculated using density functional theory approaches, and evaluated and compared to experimental data and monoreferential wave function-based methodologies. Statistical analysis of the results indicates that global-hybrid functionals allow an excellent description of absorption and emission energies, with errors around 0.05 eV, while range-separated variants led to somewhat larger errors in the range 0.1-0.2 eV. In contrast, range-separated functionals display excellent phosphorescence energies with errors close to 0.05 eV, in this case global-hybrids showing increased discrepancies around 0.5-0.1 eV.

Minimal Auxiliary Basis Set Approach for the Electronic Excitation Spectra of Organic Molecules

We report a minimal auxiliary basis model for time-dependent density functional theory (TDDFT) with hybrid density functionals that can accurately reproduce excitation energies and absorption spectra from TDDFT while reducing cost by about 2 orders of magnitude. Our method, dubbed TDDFT-ris, employs the resolution-of-the-identity technique with just one s-type auxiliary basis function per atom for the linear response operator, where the Gaussian exponents are parametrized across the periodic table using tabulated atomic radii with a single global scaling factor. By tuning on a small test set, we determine a single functional-independent scale factor that balances errors in excitation energies and absorption spectra. Benchmarked on organic molecules and compared to standard TDDFT, TDDFT-ris has an average energy error of only 0.06 eV and yields absorption spectra in close agreement with TDDFT. Thus, TDDFT-ris enables simulation of realistic absorption spectra in large molecules that would be inaccessible from standard TDDFT.

Characterization of excited states in time-dependent density functional theory using localized molecular orbitals

Localized molecular orbitals are often used for the analysis of chemical bonds, but they can also serve to efficiently and comprehensibly compute linear response properties. While conventional canonical molecular orbitals provide an adequate basis for the treatment of excited states, a chemically meaningful identification of the different excited-state processes is difficult within such a delocalized orbital basis. In this work, starting from an initial set of supermolecular canonical molecular orbitals, we provide a simple one-step top-down embedding procedure for generating a set of orbitals, which are localized in terms of the supermolecule but delocalized over each subsystem composing the supermolecule. Using an orbital partitioning scheme based on such sets of localized orbitals, we further present a procedure for the construction of local excitations and charge-transfer states within the linear response framework of time-dependent density functional theory (TDDFT). This procedure provides direct access to approximate diabatic excitation energies and, under the Tamm-Dancoff approximation, also their corresponding electronic couplings-quantities that are of primary importance in modeling energy transfer processes in complex biological systems. Our approach is compared with a recently developed diabatization procedure based on subsystem TDDFT using projection operators, which leads to a similar set of working equations. Although both of these methods differ in the general localization strategies adopted and the type of basis functions (Slaters vs Gaussians) employed, an overall decent agreement is obtained.

Structure of Geobacter cytochrome OmcZ identifies mechanism of nanowire assembly and conductivity

OmcZ nanowires produced by Geobacter species have high electron conductivity (>30 S cm-1). Of 111 cytochromes present in G. sulfurreducens, OmcZ is the only known nanowire-forming cytochrome essential for the formation of high-current-density biofilms that require long-distance (>10 µm) extracellular electron transport. However, the mechanisms underlying OmcZ nanowire assembly and high conductivity are unknown. Here we report a 3.5-Å-resolution cryogenic electron microscopy structure for OmcZ nanowires. Our structure reveals linear and closely stacked haems that may account for conductivity. Surface-exposed haems and charge interactions explain how OmcZ nanowires bind to diverse extracellular electron acceptors and how organization of nanowire network re-arranges in different biochemical environments. In vitro studies explain how G. sulfurreducens employ a serine protease to control the assembly of OmcZ monomers into nanowires. We find that both OmcZ and serine protease are widespread in environmentally important bacteria and archaea, thus establishing a prevalence of nanowire biogenesis across diverse species and environments.

Intramolecular Charge Transfer and Ultrafast Nonradiative Decay in DNA-Tethered Asymmetric Nitro- and Dimethylamino-Substituted Squaraines

Molecular (dye) aggregates are a materials platform of interest in light harvesting, organic optoelectronics, and nanoscale computing, including quantum information science (QIS). Strong excitonic interactions between dyes are key to their use in QIS; critically, properties of the individual dyes govern the extent of these interactions. In this work, the electronic structure and excited-state dynamics of a series of indolenine-based squaraine dyes incorporating dimethylamino (electron donating) and/or nitro (electron withdrawing) substituents, so-called asymmetric dyes, were characterized. The dyes were covalently tethered to DNA Holliday junctions to suppress aggregation and permit characterization of their monomer photophysics. A combination of density functional theory and steady-state absorption spectroscopy shows that the difference static dipole moment (Δd) successively increases with the addition of these substituents while simultaneously maintaining a large transition dipole moment (μ). Steady-state fluorescence and time-resolved absorption and fluorescence spectroscopies uncover a significant nonradiative decay pathway in the asymmetrically substituted dyes that drastically reduces their excited-state lifetime (τ). This work indicates that Δd can indeed be increased by functionalizing dyes with electron donating and withdrawing substituents and that, in certain classes of dyes such as these asymmetric squaraines, strategies may be needed to ensure long τ, e.g., by rigidifying the π-conjugated network.

Ground- and Excited-State Dipole Moments and Oscillator Strengths of Full Configuration Interaction Quality

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

Interexcited State Photophysics I: Benchmarking Density Functionals for Computing Nonadiabatic Couplings and Internal Conversion Rate Constants

We present the first benchmarking study of nonadiabatic matrix coupling elements (NACMEs) calculated using different density functionals. Using the S₁ → S₀ transition in perylene solvated in toluene as a case study, we calculate the photophysical properties and corresponding rate constants for a variety of density functionals from each rung of Jacob's ladder. The singlet photoluminescence quantum yield (sPLQY) is taken as a measure of accuracy, measured experimentally here as 0.955. Important quantum chemical parameters such as geometries, absorption, emission, and adiabatic energies, NACMEs, Hessians, and transition dipole moments were calculated for each density functional basis set combination (data set) using density functional theory based multireference configuration interaction (DFT/MRCI) and compared to experiment where possible. We were able to derive simple relations between the TDDFT and DFT/MRCI photophysical properties; with semiempirical damping factors of ∼0.843 ± 0.017 and ∼0.954 ± 0.064 for TDDFT transition dipole moments and energies to DFT/MRCI level approximations, respectively. NACMEs were dominated by out-of-plane derivative components belonging to the center-most ring atoms with weaker contributions from perturbations along the transverse and longitudinal axes. Calculated theoretical spectra compared well to both experiment and literature, with fluorescence lifetimes between 7.1 and 12.5 ns, agreeing within a factor of 2 with experiment. Internal conversion (IC) rates were then calculated and were found to vary wildly between 10⁶-10¹⁶ s^-1 compared with an experimental rate of the order 10⁷ s^-1. Following further testing by mixing data sets, we found a strong dependence on the method used to obtain the Hessian. The 5 characterized data sets ranked in order of most promising are PBE0/def2-TZVP, ωB97XD/def2-TZVP, HCTH407/TZVP, PBE/TZVP, and PBE/def2-TZVP.

Multiple Facets of Modeling Electronic Absorption Spectra of Systems in Solution

In this Perspective, we outline the essential physicochemical aspects that need to be considered when building a reliable approach to describe absorption properties of solvated systems. In particular, we focus on how to properly model the complexity of the solvation phenomenon, arising from dynamical aspects and specific, strong solute-solvent interactions. To this end, conformational and configurational sampling techniques, such as Molecular Dynamics, have to be coupled to accurate fully atomistic Quantum Mechanical/Molecular Mechanics (QM/MM) methodologies. By exploiting different illustrative applications, we show that an effective reproduction of experimental spectral signals can be achieved by delicately balancing exhaustive sampling, hydrogen bonding, mutual polarization, and nonelectrostatic effects.

Clustering Analysis, Structure Fingerprint Analysis, and Quantum Chemical Calculations of Compounds from Essential Oils of Sunflower (Helianthus annuus L.) Receptacles

Sunflower (Helianthus annuus L.) is an appropriate crop for current new patterns of green agriculture, so it is important to change sunflower receptacles from waste to useful resource. However, there is limited knowledge on the functions of compounds from the essential oils of sunflower receptacles. In this study, a new method was created for chemical space network analysis and classification of small samples, and applied to 104 compounds. Here, t-SNE (t-Distributed Stochastic Neighbor Embedding) dimensions were used to reduce coordinates as node locations and edge connections of chemical space networks, respectively, and molecules were grouped according to whether the edges were connected and the proximity of the node coordinates. Through detailed analysis of the structural characteristics and fingerprints of each classified group, our classification method attained good accuracy. Targets were then identified using reverse docking methods, and the active centers of the same types of compounds were determined by quantum chemical calculation. The results indicated that these compounds can be divided into nine groups, according to their mean within-group similarity (MWGS) values. The three families with the most members, i.e., the d-limonene group (18), α-pinene group (10), and γ-maaliene group (nine members) determined the protein targets, using PharmMapper. Structure fingerprint analysis was employed to predict the binding mode of the ligands of four families of the protein targets. Thence, quantum chemical calculations were applied to the active group of the representative compounds of the four families. This study provides further scientific information to support the use of sunflower receptacles.

Comparison of Length, Velocity, and Symmetric Gauges for the Calculation of Absorption and Electric Circular Dichroism Spectra with Real-Time Time-Dependent Density Functional Theory

A velocity and symmetric gauge implementation for real-time time-dependent density functional theory (RT-TDDFT) in the CP2K package using a Gaussian and plane wave approach is presented, including the explicit gauge-transformed contributions due to the nonlocal part of pseudopotentials. Absorption spectra of gas-phase α-pinene are calculated in length and velocity gauges in the long-wavelength approximation for the application of a δ pulse in linear and full order. The velocity gauge implementation is also applied to a solvated uracil molecule to showcase its use within periodic boundary conditions (PBC). For the calculation of the expectation value of the electric dipole moment in PBC, both the velocity representation and the modern theory of polarization give equivalent absorption spectra if a distributed reference point is used for the nonlocal term of the velocity operator. The discussion of linear response theory takes place in a unified framework in terms of linear response functions in propagator notation, distinguishing the parts of the linear response functions associated with perturbation and response. To further investigate gauge dependence, electric circular dichroism (ECD) spectra of α-pinene were calculated either as magnetic response to an electric field perturbation, in length or velocity gauge, or as electric response to a magnetic field perturbation in the symmetric gauge. Both approaches, electric and magnetic perturbations, have been found to yield equivalent ECD spectra.

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

Accurate Calculation of Excited-State Absorption for Small-to-Medium-Sized Conjugated Oligomers: Multiconfigurational Treatment vs Quadratic Response TD-DFT

Excited-state absorption (ESA) spectra of π-conjugated compounds are frequently calculated by (quadratic response) time-dependent density functional theory, (QR) TD-DFT, often giving a reasonable representation of the experimental results despite the (known) incomplete electronic description. To investigate whether this is inherent to the method, we calculate here the ESA spectra of small-to-medium-sized oligophenylenevinylenes (nPV) and oligothiophenes (nT) using QR TD-DFT as well as CASPT2 based on CASSCF geometries. CASPT2 gives indeed a reliable, theoretically correct description of the ESA features for all compounds; the computational effort can be reduced without significant loss of accuracy using TD-DFT geometries. QR TD-DFT, based on BHandHLYP and CAM-/B3LYP functionals, fails on short nTs but provides a reasonable description for spectral positions of nPVs and long nTs. The failure on short nTs is, however, only partly due to the incomplete configuration description but, in particular, related to an improper MO description, resulting in an asymmetric energy spacing of the occupied vs unoccupied MOs in the DFT scheme. Longer nTs, on the other side, adapt approximately the MO scheme for alternant hydrocarbons just like in nPVs, while contributions by two triplet excitations combined to a singlet (which inhibits an accurate treatment of polyenes with standard TD-DFT) do not play a relevant role in the current case. For such "well-behaved" systems, a reasonable representation of ESA spectra is found at the QR TD-DFT level due to the rather small energy shifts when including higher-order excitations.

Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids

The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.

Theoretical study on the mechanism of hot excitons combined with aggregation-induced emission in efficient red fluorescent molecules

Fluorescent emitters with the hot exciton mechanism combined with aggregation induced emission (AIE) character show prospective applications in organic light emitting diodes (OLEDs). However, theoretical studies on amorphous states are limited. In this work, a theoretical study is performed on the photophysical properties of the reported compound 4-(7-(10-ethyl-10H-phenothiazin-3-yl)benzo[c][1,2,5]thiadiazol-4-yl)-N,N-diphenylaniline (PBTPA), which possesses a hot exciton mechanism and AIE. The aggregation states of this molecule in a film are given by molecular dynamics (MD) simulations, and then the photophysical properties are studied by using the QM/MM method with the consideration of the solid-state effect (SSE). The results explain the hot exciton and AIE mechanism of the molecule. First, there is a hot exciton channel between the S₁ and T₂ state of the PBTPA. Second, the conformational changes of PBTPA between the ground state and the excited state are restricted in the aggregate state. Last, in the low frequency region, the rotation motion is suppressed, and then the reorganization energy and Huang-Rhys (HR) factor in the aggregate state are much smaller. Therefore, the molecules show strong fluorescence efficiency in the aggregated state.

Going beyond the Electric-Dipole Approximation in the Calculation of Absorption and (Magnetic) Circular Dichroism Spectra including Scalar Relativistic and Spin-Orbit Coupling Effects

In this work, a time-dependent density functional theory (TD-DFT) scheme for computing optical spectroscopic properties in the framework of linearly and circularly polarized light is presented. The scheme is based on a previously formulated theory for predicting optical absorption and magnetic circular dichroism (MCD) spectra. The scheme operates in the framework of the full semi-classical field-matter interaction operator, thus generating a powerful and general computational scheme capable of computing the absorption (ABS), circular dichroism (CD), and MCD spectra. In addition, our implementation includes the treatment of relativistic effects in the framework of quasidegenerate perturbation theory, which accounts for scalar relativistic effects (in the self-consistent field step) and spin-orbit coupling (in the TD-DFT step), as well as external magnetic field perturbations. Hence, this formalism is also able to probe spin-forbidden transitions. The random orientations of molecules are taken into account by a semi-numerical approach involving a Lebedev numerical quadrature alongside analytical integration. It is demonstrated the numerical quadrature requires as few as 14 points for satisfactory converged results thus leading to a highly efficient scheme, while the calculation of the exact transition moments creates no computational bottlenecks. It is demonstrated that at zero magnetic field, the CD spectrum is recovered while the sum of left and right circularly polarized light contributions provides the linear absorption spectrum. The virtues of this efficient and general protocol are demonstrated on a selected set of organic molecules where the various contributions to the spectral intensities have been analyzed in detail.

Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck–Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck–Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck–Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.

Comparative Study of Vibrational Raman Optical Activity with Different Time-Dependent Density Functional Approximations: The VROA36 Database

A new database, VROA36, is introduced to investigate the performance of computational approaches for vibrational Raman optical activity (VROA) calculations. The database is composed of 36 molecules with known experimental VROA spectra. It includes 93 conformers. Normal modes calculated with B3LYP-D3(BJ)/def2-TZVP are used to compute the VROA spectra with four functionals, B3LYP-D3(BJ), ωB97X-D, M11, and optimally tuned LC-PBE, as well as several basis sets. SimROA indices and frequency scaling factors are used to compare calculated spectra with each other and with experimental data. The four functionals perform equally well independently of the basis set and usually achieve good agreement with the experimental data. For molecules in near- or at-resonance conditions, the inclusion of a complex (damped) linear response approach is important to obtain physically meaningful VROA intensities. The use of any of the tested functional approximations with the def2-SVPD Gaussian-type basis set, or a basis of similar flexibility, can be recommended for efficient and reliable theoretical VROA studies.