Assessment of TD-DFT methods and of various spin scaled CIS(D) and CC2 versions for the treatment of low-lying valence excitations of large organic dyes

Excited states using the simplified Tamm-Dancoff-Approach for range-separated hybrid density functionals: development and application

The recently introduced sTDA methodology [S. Grimme, J. Chem. Phys., 2013, 138, 244104] to compute excitation spectra of huge molecular systems is extended to range-separated hybrid (RSH) density functionals. The three empirical parameters of the method which describe a screened two-electron interaction are obtained for some common RSH functionals (ωB97 family, CAM-B3LYP, LC-BLYP) from a fit to theoretical SCS-CC2 reference vertical excitation energies for a set of small to medium-sized chromophores. The method is cross-validated on a set of inter- and intramolecular charge transfer states and a set composed of typical valence transitions. Overall small deviations from reference data of only about 0.2-0.4 eV are found with best performance for CAM-B3LYP and ωB97X-D3. To demonstrate versatility and robustness of the new methodology, applications (the UV/Vis spectrum of the pyridine polymer and the ECD spectrum of (P)-[11]helicene) and frequently used charge transfer examples are discussed. In one case, 11 000+ excited electronic states of a system containing 330 atoms were calculated. We show that the asymptotically correct sTDA-RSH combination yields results often superior to those based on global hybrids and that it opens up new possibilities for the computation of excited states in materials science and bio-molecular systems.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

Taking up the cyanine challenge with quantum tools

Cyanine derivatives, named from the Greek word kyanos meaning dark-blue, were discovered more than 150 years ago and remain one of the most widely used classes of organic dyes with contemporary applications in photography (panchromatic emulsions), information storage (CD-R and DVD-R media) and biochemistry (DNA and protein labeling) fields. Cyanine chromogens consist of a charged π-conjugated segment containing an odd number of sp(2) carbon atoms with the chain capped at the extremities by two electronegative centers, typically nitrogen or oxygen atoms. Cyanines are characterized by a vanishing bond length alternation indicating nearly equal carbon-carbon bond lengths, as well as a very intense and sharp absorption band presenting a shoulder. This hallmark band undergoes a strong red shift when the chain is extended. This so-called vinyl shift is extremely large (ca. 100 nm for each pair of carbon atoms added in the π-conjugated path), making cyanines ideal building blocks for the design of devices with near-infrared applications. Numerous cyanines also exhibit emission bands with large quantum yields. These exceptional optical properties explain why both canonical cyanines and the corresponding fluoroborates (e.g., boron-dipyrromethene, BODIPY) remain the focus of an ever-growing body of experimental work. In turn, this popularity has stimulated quantum mechanical investigations aiming, on the one hand, at probing the specific electronic nature of cyanine dyes and, on the other hand, at helping to design new dyes. However, the adiabatic approximation to time-dependent density functional theory, the most widespread ab initio model for electronically excited states, fails to accurately reproduce the absorption spectra of cyanine derivatives: it yields a systematic and large underestimation of the experimental wavelengths irrespective of the details of the computational protocol. In contrast, highly correlated wave function approaches provide accurate transition energies for model systems but are hardly applicable to real-life cyanines and BODIPY. This indicates that setting up a computationally tractable theoretical protocol that provides both robust and accurate optical spectra for cyanine-based dyes is a major challenge that has only been taken up lately. In this Account, we compile the most recent advances in the field by considering both compact streptocyanines and large fluoroborates. For the former, we summarize the key results obtained with a large panel of theoretical approaches, allowing us not only to understand the origin of the cyanine challenge but also to pinpoint the schemes presenting the most promising accuracy/effort ratio. For the latter, we show via selected examples how theoretical models can be used to reproduce simultaneously experimental band shapes and transition energies, thus paving the way to an efficient in silico design of new compounds.

Photoisomerization action spectroscopy: flicking the protonated merocyanine–spiropyran switch in the gas phase

Molecular ions drifting through buffer gas are exposed to tunable laser light causing a detectable change in their collision cross section and arrival time at the ion detector. The effect can be exploited to obtain action electronic spectra for selected molecular isomers in the gas phase.

The relationship between the boron dipyrromethene (BODIPY) structure and the effectiveness of homogeneous and heterogeneous solar hydrogen-generating systems as well as DSSCs

BODIPY photosensitizers were used to investigate the relationship between structure and effectiveness of visible-light-driven hydrogen production as well as DSSCs.

Electronic couplings for molecular charge transfer: benchmarking CDFT, FODFT, and FODFTB against high-level ab initio calculations

We introduce a database (HAB11) of electronic coupling matrix elements (H(ab)) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H(ab) values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.

White emitters by tuning the excited-state intramolecular proton-transfer fluorescence emission in 2-(2'-hydroxybenzofuran)benzoxazole dyes

The synthesis, structural, and photophysical properties of a new series of original dyes based on 2-(2'-hydroxybenzofuran)benzoxazole (HBBO) is reported. Upon photoexcitation, these dyes exhibit intense dual fluorescence with contribution from the enol (E*) and the keto (K*) emission, with K* being formed through excited-state intramolecular proton transfer (ESIPT). We show that the ratio of emission intensity E*/K* can be fine-tuned by judiciously decorating the molecular core with electron-donating or -attracting substituents. Push-pull dyes 9 and 10 functionalized by a strong donor (nNBu2 ) and a strong acceptor group (CF3 and CN, respectively) exhibit intense dual emission, particularly in apolar solvents such as cyclohexane in which the maximum wavelength of the two bands is the more strongly separated. Moreover, all dyes exhibit strong solid-state dual emission in a KBr matrix and polymer films with enhanced quantum yields reaching up to 54 %. A wise selection of substituents led to white emission both in solution and in the solid state. Finally, these experimental results were analyzed by time-dependent density functional theory (TD-DFT) calculations, which confirm that, on the one hand, only E* and K* emission are present (no rotamer) and, on the other hand, the relative free energies of the two tautomers in the excited state guide the ratio of the E*/K* emission intensities.

Application of recent double-hybrid density functionals to low-lying singlet-singlet excitation energies of large organic compounds

The present work assesses some recently developed double-hybrid density functionals (B2π-PLYP, PBE0-DH, and PBE0-2) using linear-response Tamm-Dancoff Time-Dependent Density Functional Theory. This assessment is achieved against experimentally derived low-lying excitation energies of large organic dyes of recent interest, including some excitations dominated by charge-transfer transitions. Comparisons are made with some of the best-performing methods established from the literature, such as PBE0 or B3LYP hybrid or the recently proposed B2-PLYP and B2GP-PLYP double-hybrid models, to ascertain their quality and robustness on equal footing. The accuracy of parameter-free or empirical forms of double-hybrid functionals is also briefly discussed. Generally speaking, it turns out that double-hybrid expressions always provide more accurate estimates than corresponding hybrid methods. Double-hybrid functionals actually reach averaged accuracies of 0.2 eV, that can be admittedly considered close to any intended accuracy limit within the present theoretical framework.

Exploring molecular structures, orbital interactions, intramolecular proton-transfer reaction kinetics, electronic transitions and complexation of 3-hydroxycoumarin species using DFT methods

Optimal structures and electronic properties of various species of 3-hydroxycoumarin (3-HCou) have been explored using density functional theory (DFT) methods under polarizable continuum model (PCM) of solvation. Electron transfer from pyrone to benzene moieties is enhanced upon deprotonation. Anionic and radical species have similar orbital-interaction characteristics but the charges in the former are distributed more uniformly. The rate of intramolecular proton transfer for the neutral species increases many folds upon excitation. The HOMO-LUMO transition with π→π* character mainly accounts for the UV absorption of most 3-HCou species in solution. The wavelengths of maximal absorption predicted using TD-DFT method are in agreement with the previous experiment. For the charged species, calculations with the range-corrected functional yield better agreement with the previous experiment. Anionic 3-HCou species shows high degrees of complexation with chromium(III) and copper(II) compared with oxovanadium(IV) and zinc(II). Either oxovanadium(IV) or zinc(II) prefers forming two isomeric complexes with comparable degrees of formation.

Expanding the palette of phenanthridinium cations

5,6-Disubstituted phenanthridinium cations have a range of redox, fluorescence and biological properties. Some properties rely on phenanthridiniums intercalating into DNA, but the use of these cations as exomarkers for the reactive oxygen species (ROS), superoxide, and as inhibitors of acetylcholine esterase (AChE) do not require intercalation. A versatile modular synthesis of 5,6-disubstituted phenanthridiniums that introduces diversity by Suzuki–Miyaura coupling, imine formation and microwave-assisted cyclisation is presented. Computational modelling at the density functional theory (DFT) level reveals that the novel displacement of the aryl halide by an acyclic N-alkylimine proceeds by an S(N)Ar mechanism rather than electrocyclisation. It is found that the displacement of halide is concerted and there is no stable Meisenheimer intermediate, provided the calculations consistently use a polarisable solvent model and a diffuse basis set.

Performance of an Optimally Tuned Range-Separated Hybrid Functional for 0-0 Electronic Excitation Energies

Using a set of 40 conjugated molecules, we assess the performance of an "optimally tuned" range-separated hybrid functional in reproducing the experimental 0-0 energies. The selected protocol accounts for the impact of solvation using a corrected linear-response continuum approach and vibrational corrections through calculations of the zero-point energies of both ground and excited-states and provides basis set converged data thanks to the systematic use of diffuse-containing atomic basis sets at all computational steps. It turns out that an optimally tuned long-range corrected hybrid form of the Perdew-Burke-Ernzerhof functional, LC-PBE*, delivers both the smallest mean absolute error (0.20 eV) and standard deviation (0.15 eV) of all tested approaches, while the obtained correlation (0.93) is large but remains slightly smaller than its M06-2X counterpart (0.95). In addition, the efficiency of two other recently developed exchange-correlation functionals, namely SOGGA11-X and ωB97X-D, has been determined in order to allow more complete comparisons with previously published data.

A simplified Tamm-Dancoff density functional approach for the electronic excitation spectra of very large molecules

Two approximations in the Tamm-Dancoff density functional theory approach (TDA-DFT) to electronically excited states are proposed which allow routine computations for electronic ultraviolet (UV)- or circular dichroism (CD) spectra of molecules with 500-1000 atoms. Speed-ups compared to conventional time-dependent DFT (TD-DFT) treatments of about two to three orders of magnitude in the excited state part at only minor loss of accuracy are obtained. The method termed sTDA ("s" for simplified) employs atom-centered Löwdin-monopole based two-electron repulsion integrals with the asymptotically correct 1∕R behavior and perturbative single excitation configuration selection. It is formulated generally for any standard global hybrid density functional with given Fock-exchange mixing parameter ax. The method performs well for two standard benchmark sets of vertical singlet-singlet excitations for values of ax in the range 0.2-0.6. The mean absolute deviations from reference data are only 0.2-0.3 eV and similar to those from standard TD-DFT. In three cases (two dyes and one polypeptide), good mutual agreement between the electronic spectra (up to 10-11 eV excitation energy) from the sTDA method and those from TD(A)-DFT is obtained. The computed UV- and CD-spectra of a few typical systems (e.g., C60, two transition metal complexes, [7]helicene, polyalanine, a supramolecular aggregate with 483 atoms and about 7000 basis functions) compare well with corresponding experimental data. The method is proposed together with medium-sized double- or triple-zeta type atomic-orbital basis sets as a quantum chemical tool to investigate the spectra of huge molecular systems at a reliable DFT level.

Vibronic spectra of organic electronic chromophores

Vibronic effects in organic electronic building blocks are presented.

Intramolecular excimer formation in hexakis(pyrenyloxy)cyclotriphosphazene: photophysical properties, crystal structure, and theoretical investigation

A hexakis(pyrenyloxy)cyclotriphosphazene was synthesized by the reaction of N₃P₃Cl₆ with 2-hydroxypyrene, and its excimer emission has been investigated.

Noble-metal-free BODIPY–cobaloxime photocatalysts for visible-light-driven hydrogen production

Four BODIPY–cobaloxime systems Co-Bn (n = 1–4) were studied as single-component homogeneous photocatalytic systems for H₂ generation in aqueous media.

Excited-states of BODIPY–cyanines: ultimate TD-DFT challenges?

Several computational approaches are used to mimic the excited-state properties of twelve large BODIPY–cyanine dyes.

Ultrafast shock initiation of exothermic chemistry in hydrogen peroxide

We report observations of shock compressed, unreacted hydrogen peroxide at pressures up to the von Neumann pressure for a steady detonation wave, using ultrafast laser-driven shock wave methods. At higher laser drive energy we find evidence of exothermic chemical reactivity occurring in less than 100 ps after the arrival of the shock wave in the sample. The results are consistent with our MD simulations and analysis and suggest that reactivity in hydrogen peroxide is initiated on a sub-100 ps time scale under conditions found just subsequent to the lead shock in a steady detonation wave.

SpecDis: quantifying the comparison of calculated and experimental electronic circular dichroism spectra

This article outlines theory and practice of the comparison of calculated and experimental electronic circular dichroism (ECD) curves to determine the absolute configuration of chiral molecules. The focus is on the evaluation of excited-state calculations giving hints at the identification of the correct bandwidth and the application of the so-called "UV shift" as a correction factor. A similarity factor is introduced, which helps to quantify the degree of matching of curves. In addition, a few common errors are described that can be made during the measurements of ECD and UV spectra-and advice is given of how to avoid these mistakes. All equations mentioned in the article are implemented in our SpecDis software, which has been developed to rapidly compare calculated ECD and UV curves with experimental ones, and to produce graphics in publication quality.