TDDFT diagnostic testing and functional assessment for triazene chromophores

Visualizing and characterizing excited states from time-dependent density functional theory

Time-dependent density functional theory (TD-DFT) is the most widely-used electronic structure method for excited states, due to a favorable combination of low cost and semi-quantitative accuracy in many contexts, even if there are well recognized limitations. This Perspective describes various ways in which excited states from TD-DFT calculations can be visualized and analyzed, both qualitatively and quantitatively. This includes not just orbitals and densities but also well-defined statistical measures of electron-hole separation and of Frenkel-type exciton delocalization. Emphasis is placed on mathematical connections between methods that have often been discussed separately. Particular attention is paid to charge-transfer diagnostics, which provide indicators of when TD-DFT may not be trustworthy due to its categorical failure to describe long-range electron transfer. Measures of exciton size and charge separation that are directly connected to the underlying transition density are recommended over more ad hoc metrics for quantifying charge-transfer character.

Exploration of intramolecular charge transfer in para-substituted nitrobenzofurazan: Experimental and theoretical analyses

The present work aims at exploring the high electrophilic character of 4-chloro-7-nitrobenzo-2-oxa-1,3-diazole (NBD-Cl) toward the morpholine group by an S_NAr reaction in acetonitrile or water (thereafter referred to as NBD-Morph). The electron-donating ability of the morpholine causes intra-molecular charge transfer (ICT). In this report, we present a comprehensive study on the optical characteristics using UV-Vis, photoluminescence (cw-PL) and its time-resolved (TR-PL) to determine the properties of the emissive intramolecular charge transfer (ICT) in the NBD-Morph donor-acceptor system. An exhaustive theoretical investigation utilizing the density functional theory (DFT) and its extension TD-DFT methods is an essential complement of experiments to rationalize and understand the molecular structure and related properties. The findings from QTAIM, ELF, and RDG analyses establish that the bonding between morpholine and NBD moieties is of the electrostatic or hydrogen bond type. In addition, the Hirshfeld surfaces have been established to explore the types of interactions. Further, the non-linear optical (NLO) responses of the compound have been examined. The structure-property relationships obtained through the combined experimental and theoretical offer valuable insights for designing efficient NLO material.

A Computational Protocol for Vibrational Circular Dichroism Spectra of Cyclic Oligopeptides

Cyclic peptides are a promising class of compounds for next-generation antibiotics as they may provide new ways of limiting antibiotic resistance development. Although their cyclic structure will introduce some rigidity, their conformational space is large and they usually have multiple chiral centers that give rise to a wide range of possible stereoisomers. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) are used to assign stereochemistry and discriminate enantiomers of chiral molecules, often in combination with electronic structure methods. The reliable determination of the absolute configuration of cyclic peptides will require robust computational methods than can identify all significant conformers and their relative population and reliably assign their stereochemistry from their chiroptical spectra by comparison with ab initio calculated spectra. We here present a computational protocol for the accurate calculation of the VCD spectra of a series of flexible cyclic oligopeptides. The protocol builds on the Conformer-Rotamer Ensemble Sampling Tool (CREST) developed by Grimme and co-workers ( Phys. Chem. Chem. Phys. 2020, 22, 7169-7192 and J. Chem. Theory. Comput. 2019, 15, 2847-2862) in combination with postoptimizations using B3LYP and moderately sized basis sets. Our recommended computational protocol for the computation of VCD spectra of cyclic oligopeptides consists of three steps: (1) conformational sampling with CREST and tight-binding density functional theory (xTB); (2) energy ranking based on single-point energy calculations as well as geometry optimization and VCD calculations of conformers that are within 2.5 kcal/mol of the most stable conformer using B3LYP/6-31+G*/CPCM; and (3) VCD spectra generation based on Boltzmann weighting with Gibbs free energies. Our protocol provides a feasible basis for generating VCD spectra also for larger cyclic peptides of biological/pharmaceutical interest and can thus be used to investigate promising compounds for next-generation antibiotics.

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

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

Are all charge-transfer parameters created equally? A study of functional dependence and excited-state charge-transfer quantification across two dye families

Small molecule organic dyes have many potential uses in medicine, textiles, forensics, and light-harvesting technology. Being able to computationally predict the spectroscopic properties of these dyes could greatly expedite screening efforts, saving time and materials. Time-dependent density functional theory (TD-DFT) has been shown to be a good tool for this in many instances, but characterizing electronic excitations with charge-transfer (CT) character has historically been challenging and can be highly sensitive to the chosen exchange-correlation functional. Here we present a combined experimental and computational study of the excited-state electronic structure of twenty organic dyes obtained from the Max Weaver Dye Library at NCSU. Results of UV-vis spectra calculations on these dyes with six different exchange-correlation functionals, BP86, B3LYP, PBE0, M06, BH and HLYP, and CAM-B3LYP, were compared against their measured UV-vis spectra. It was found that hybrid functionals with modest amounts (20-30%) of included Hartree-Fock exchange are the most effective at matching the experimentally determined λ_max. The interplay between the observed error, the functional chosen, and the degree of CT was analyzed by quantifying the CT character of λ_max using four orbital and density-based metrics, Λ, Δr, S_C and D_CT, as well as the change in the dipole moment, Δμ. The results showed that the relationship between CT character and the functional dependence of error is not straightforward, with the observed behavior being dependent both on how CT was quantified and the functional groups present in the molecules themselves. It is concluded that this may be a result of the examined excitations having intermediate CT character. Ultimately it was found that the nature of the molecular "family" influenced how a given functional behaved as a function of CT character, with only two of the examined CT quantification methods, Δr and D_CT, showing consistent behavior between the different molecular families. This suggests that further work needs to be done to ensure that currently used CT quantification methods show the same general trends across large sets of multiple dye families.

Sulfur-centered hemi-bond radicals as active intermediates in S-DNA phosphorothioate oxidation

Phosphorothioate (PS) modifications naturally appear in bacteria and archaea genome and are widely used as antisense strategy in gene therapy. But the chemical effects of PS introduction as a redox active site into DNA (S-DNA) is still poorly understood. Herein, we perform time-resolved spectroscopy to examine the underlying mechanisms and dynamics of the PS oxidation by potent radicals in free model, in dinucleotide, and in S-oligomer. The crucial sulphur-centered hemi-bonded intermediates -P–S∴S–P- were observed and found to play critical roles leading to the stable adducts of -P–S–S–P-, which are backbone DNA lesion products. Moreover, the oxidation of the PS moiety in dinucleotides d[G_PSG], d[A_PSA], d[G_PSA], d[A_PSG] and in S-oligomers was monitored in real-time, showing that PS oxidation can compete with adenine but not with guanine. Significantly, hole transfer process from A^+• to PS and concomitant -P–S∴S–P- formation was observed, demonstrating the base-to-backbone hole transfer unique to S-DNA, which is different from the normally adopted backbone-to-base hole transfer in native DNA. These findings reveal the distinct backbone lesion pathway brought by the PS modification and also imply an alternative -P–S∴S–P-/-P–S–S–P- pathway accounting for the interesting protective role of PS as an oxidation sacrifice in bacterial genome.

Effects of defects in g-C3N4 on excited-state charge distribution and transfer: Potential for improved photocatalysis

Graphite phase carbon nitride (g-C₃N₄) with triazine ring structures is a polymeric metal-free semiconductor with a medium bandgap and two-dimensional layered structure. g-C₃N₄ has attracted attention because of its photocatalytic applications, such as the photodegradation of pollution and hydrogen production via water splitting. Defective elements and sites are two essential factors in rationally designing highly-efficient photocatalysts based on g-C₃N₄ at the nanoscale. When the molecule absorbs energy and enters an excited state, electrons migrate and the charge distribution changes accordingly. The properties of the excited states of g-C₃N₄ are related to its defect elements and sites. Therefore, it is necessary to understand the effects of defects on excited states in the design of g-C₃N₄ catalysts. In this paper, the excited-state characteristics of intrinsic g-C₃N₄ and g-C₃N₄ with C- and N-atom defects are analyzed by density functional theory. We apply quantum chemistry and wave function analysis to determine the hole-electron distributions and charge transfer directions. To measure and discuss the characteristics of electron excitation using quantitative numerical methods, the D, Sr, H, and t indices are calculated. Our results promote a deeper understanding of the roles of defective elements and sites in photocatalysis by g-C₃N₄.

Influence of nonmetal dopants on charge separation of graphitic carbon nitride by time-dependent density functional theory

We examined the charge separation performance of phosphorus-, oxygen-, and sulfur-doped g-C₃N₄ by using time dependent density functional theory and wave function analysis.

Optical excitations of chlorophyll a and b monomers and dimers

A necessary first step in the development of technologies such as artificial photosynthesis is understanding the photoexcitation process within the basic building blocks of naturally occurring light harvesting complexes (LHCs). The most important of these building blocks in biological LHCs such as LHC II from green plants are the chlorophyll a (Chl a) and chlorophyll b (Chl b) chromophores dispersed throughout the protein matrix. However, efforts to describe such systems are still hampered by the lack of computationally efficient and accurate methods that are able to describe optical absorption in large biomolecules. In this work, we employ a highly efficient linear combination of atomic orbitals (LCAOs) to represent the Kohn-Sham (KS) wave functions at the density functional theory (DFT) level and perform time-dependent density functional theory (TDDFT) calculations in either the reciprocal space and frequency domain (LCAO-TDDFT-k-ω) or real space and time domain (LCAO-TDDFT-r-t) of the optical absorption spectra of Chl a and b monomers and dimers. We find that our LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations reproduce results obtained with a plane-wave (PW) representation of the KS wave functions (PW-TDDFT-k-ω) but with a significant reduction in computational effort. Moreover, by applying the Gritsenko, van Leeuwen, van Lenthe, and Baerends solid and correlation derivative discontinuity correction Δ_x to the KS eigenenergies, with both LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t methods, we are able to semiquantitatively reproduce the experimentally measured photoinduced dissociation results. This work opens the path to first principles calculations of optical excitations in macromolecular systems.

Davydov-type excitonic effects on the absorption spectra of parallel-stacked and herringbone aggregates of pentacene: Time-dependent density-functional theory and time-dependent density-functional tight binding

Exciton formation leads to J-bands in solid pentacene. Describing these exciton bands represents a challenge for both time-dependent (TD) density-functional theory (DFT) and for its semi-empirical analog, namely, for TD density-functional tight binding (DFTB) for three reasons: (i) solid pentacene and pentacene aggregates are bound only by van der Waals forces which are notoriously difficult to describe with DFT and DFTB, (ii) the proper description of the long-range coupling between molecules, needed to describe Davydov splitting, is not easy to include in TD-DFT with traditional functionals and in TD-DFTB, and (iii) mixing may occur between local and charge transfer excitons, which may, in turn, require special functionals. We assess how far TD-DFTB has progressed toward a correct description of this type of exciton by including both a dispersion correction for the ground state and a range-separated hybrid functional for the excited state and comparing the results against corresponding TD-CAM-B3LYP/CAM-B3LYP+D3 results. Analytic results for parallel-stacked ethylene are derived which go beyond Kasha's exciton model [M. Kasha, H. R. Rawls, and A. El-Bayoumi, Pure Appl. Chem. 11, 371 (1965)] in that we are able to make a clear distinction between charge transfer and energy transfer excitons. This is further confirmed when it is shown that range-separated hybrids have a markedly greater effect on charge-transfer excitons than on energy-transfer excitons in the case of parallel-stacked pentacenes. TD-DFT calculations with the CAM-B3LYP functional and TD-lc-DFT calculations lead to negligible excitonic corrections for the herringbone crystal structure, possibly because of an overcorrection of charge-transfer effects (CAM refers to Coulomb attenuated method). In this case, TD-DFT calculations with the B3LYP functional or TD-DFTB calculations parameterized to B3LYP give the best results for excitonic corrections for the herringbone crystal structure as judged from comparison with experimental spectra and with Bethe-Salpeter equation calculations from the literature.

TD-DFT calculations of one- and two-photon absorption in Coumarin C153 and Prodan: attuning theory to experiment

We use TD-DFT to calculate the one-photon absorption (1PA) and two-photon absorption (2PA) properties of C153 and Prodan in toluene and DMSO, and benchmark different methods relative to accurate experimental data available from the literature on these particular systems. As the first step, we modify the range-separated TD-DFT to provide the best prediction for the peak 1PA wavelength, and then apply the optimized functionals to achieve quantitative predictions of the corresponding two-photon absorption cross section, σ_2PA, with an accuracy ∼10-20% in C153 and ∼20-30% in Prodan. To elucidate the origin of residual discrepancies between the theory and experimental observations, we invoked the two essential states model for σ_2PA, which allows us to verify not only the transition wavelength and the σ_2PA value, but also to quantitatively benchmark the calculation of key molecular parameters such as the transition dipole moment and the change of the permanent dipole moment. Such comprehensive cross-checking indicates that a larger discrepancy in Prodan is most likely caused by a noted failure of DFT to predict the relative intensity and relative ordering of closely lying excited states with different degrees of intramolecular charge transfer, which we further support by analyzing the predictions obtained by high-level coupled-cluster calculations in the gas phase. Our results highlight the utility of benchmarking the calculations not only relative to other theoretical methods, but also in comparison to the experimental measurements, wherever such data are available.

Rationalisation of the optical signatures of nor-dihydroxanthene-hemicyanine fused near-infrared fluorophores by first-principle tools

The perfect pair towards more effective nor-DXH-hemicyanine fused NIR fluoroohores.

Capturing the radical ion-pair intermediate in DNA guanine oxidation

Although the radical ion pair has been frequently invoked as a key intermediate in DNA oxidative damage reactions and photoinduced electron transfer processes, the unambiguous detection and characterization of this species remain formidable and unresolved due to its extremely unstable nature and low concentration. We use the strategy that, at cryogenic temperatures, the transient species could be sufficiently stabilized to be detectable spectroscopically. By coupling the two techniques (the cryogenic stabilization and the time-resolved laser flash photolysis spectroscopy) together, we are able to capture the ion-pair transient G+•⋯Cl- in the chlorine radical-initiated DNA guanine (G) oxidation reaction, and provide direct evidence to ascertain the intricate type of addition/charge separation mechanism underlying guanine oxidation. The unique spectral signature of the radical ion-pair G+•⋯Cl- is identified, revealing a markedly intense absorption feature peaking at 570 nm that is distinctive from G+• alone. Moreover, the ion-pair spectrum is found to be highly sensitive to the protonation equilibria within guanine-cytosine base pair (G:C), which splits into two resolved bands at 480 and 610 nm as the acidic proton transfers along the central hydrogen bond from G+• to C. We thus use this exquisite sensitivity to track the intrabase-pair proton transfer dynamics in the double-stranded DNA oligonucleotides, which is of critical importance for the description of the proton-coupled charge transfer mechanisms in DNA.

Exploring Redox States, Doping and Ordering of Electroactive Star-Shaped Oligo(aniline)s

We have prepared a simple star-shaped oligo(aniline) (TDPB) and characterised it in detail by MALDI-TOF MS, UV/Vis/NIR spectroscopy, time-dependent DFT, cyclic voltammetry and EPR spectroscopy. TDPB is part of an underdeveloped class of π-conjugated molecules with great potential for organic electronics, display and sensor applications. It is redox active and reacts with acids to form radical cations. Acid-doped TDPB shows behaviour similar to discotic liquid crystals, with X-ray scattering investigations revealing columnar self-assembled arrays. The combination of unpaired electrons and supramolecular stacking suggests that star-shaped oligo(aniline)s like TDPB have the potential to form conducting nanowires and organic magnetic materials.

Truncated Transition Densities for Analysis of (Nonlinear) Optical Properties of carbo-Chromophores

The optical properties of several quadrupolar carbo-benzene derivatives are investigated at various levels of calculation (TDDFT and CASPT2) and analyzed using a new theoretical tool here disclosed: The "visualization" of the transition dipole moment from the transition density truncated to the main monoexcitations involved in the electronic transition (TTD). The experimental or calculated one-photon UV-visible absorption spectra of the carbo-benzene derivatives fit with the Gouterman model originally proposed for porphyrins, where the first four excited states involve linear combinations of monoexcitations of the same four frontier molecular orbitals. The relative intensities of the absorption bands are analyzed from the transition dipole moments calculated from the TTDs and an analogy between porphyrins and carbo-benzenes is argued. The two-photon absorption (TPA) cross section related to the third-order nonlinear optical response is calculated for each two-photon-allowed excited state |f⟩ from the contribution of all possible intermediate excited states |i⟩ using the "sum-over-state" (SOS) scheme. The quadrupolar carbo-benzene derivatives fit into the three-level model, as their TPA cross section exhibits a dominant contribution of one of the intermediate excited states. The origin of TPA efficiency (enhancement) upon carbo-merisation of the C-C link to the para-substituents is discussed from the excitation energies of the intermediate and final excited states and from the two corresponding transition dipole moments (μ0i and μif). The latter may be calculated from the TTDs.

Tailoring of spectral response and intramolecular charge transfer in β-enaminones through band gap tuning: synthesis, spectroscopy and quantum chemical studies

In this paper, we investigate the synthetic tailoring of the spectral response and intramolecular charge transfer (ICT) of β-enaminones through bandgap modulation.

Rational design of cyclopenta[b]naphthalenes for better optoelectronic applications and their photophysical properties using DFT/TD-DFT methods

The optical properties of cyclopenta[b]naphthalenes (CPNs) can be fine-tuned by suitable substitutions and DFT calculations show that they can make efficient OLEDs.

Modeling the photosensitizing properties of thiolate-protected gold nanoclusters

An accurate computational strategy for studying the structural, redox and optical properties of thiolated gold nanoclusters (GNCs) using (time-dependent) density functional theory is proposed.

Chemical trends in the optical properties of rocksalt nanoparticles

The nature and magnitude of the optical gaps of rocksalt alkaline earth (MgO, CaO, SrO, MgS, MgSe) and transition metal chalcogenide (CdO, PbS) nanoparticles are studied using time-dependent density functional theory (TD-DFT) calculations on (MX)32 nanoparticles. We demonstrate, just as we previously showed for MgO, that TD-DFT calculations on rocksalt nanoparticles require the use of hybrid exchange-correlation (XC-)functionals with a high percentage of Hartree-Fock like exchange (e.g. BHLYP) or range-separated XC-functionals to circumvent problems related to the description of charge-transfer excitations. Concentrating on the results obtained with TD-BHLYP we show that the optical gap in rocksalt nanoparticles displays a wide range of behavior; ranging from large optical gaps stemming from a localized excitation involving corner atoms in alkaline earth oxides to a delocalized excitation and small optical gaps in the transition metal chalcogenides. Finally, we rationalize this wide range of behaviour in terms of differences in the degree to which the Coulombic interaction between the excited electron and hole is screened in the different nanoparticles, and relate it to the optical dielectric constants of the bulk materials the nanoparticles are made from.

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.

Enhancement of nonlinear optical (NLO) properties of indigo through modification of auxiliary donor, donor and acceptor

In this study, indigo based dyes with high non-linear optical response have been investigated. Density functional theory (DFT) was used to study non-linear optical properties of indigo and newly designed dyes (IM-Dye-0, IM-Dye-1, IM-Dye-2 and IM-Dye-3). The time dependant density functional theory (TDDFT) was used to calculate the excitation energies. The HOMO-LUMO energy gaps of newly designed dyes were smaller as compare with indigo dye. Absorption maxima of newly designed dyes strongly red shifted as compare with indigo dye. High non-linear optical (NLO) response of newly designed dyes revealed that these materials would be excellent for NLO applications. This theoretical approach of designing will pave the way for experimentalists to synthesize high response NLO compound.