General Density-Based Index to Analyze Charge Transfer Phenomena: From Models to Butterfly Molecules

Earth Mover's Charge Transfer Distance: A General and Robust Approach for Describing Excited State Locality

A novel approach for assessing the extent of electron displacement in optical transitions is proposed by implementing the Earth Mover's Distance (EMD) method, which quantifies the spatial dissimilarity between ground and excited state electron density distributions. In contrast to previous descriptors, this index provides a representative and intuitively understandable distance under a robust and computationally efficient scheme for all possible forms of locality, even in the most difficult to dissect topological cases. The theoretical differences among the existing indices and our method are first illustrated with the help of a simplified model system, followed by a benchmarking of several partial atomic charge models using experimentally relevant push-pull compounds with diverse symmetries. These same molecules are finally employed to further demonstrate the principal advantages of the EMD index and its capabilities in rationalizing charge transfer phenomena.

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.

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.

Metallo-dithiaporphyrin pigments for bulk-heterojunction solar cell applications: ab initio investigation of structural and optoelectronic properties

Metallo-dithiaporphyrin small molecules have been designed by substituting Ru(ii) with various transition metals at the same oxidation state (M = Mn, Fe, Ni, Cu) as donor materials for Bulk Heterojunction Organic Solar Cells (BHJ-OSCs). Density functional theory (DFT) and time-dependent DFT (TD-DFT) have been used to study the optoelectronic properties of metallo-dithiaporphyrin at various functionals and basis sets. We discovered that the open-circuit voltage (VOC) value increases when Ru(ii) in Ru(S2TTP)Cl2 (S2TTP = tetra-p-tolyldithiaporphyrin) is substituted. In addition, the light harvesting efficiency (LHE) of nickel, manganese, and iron complexes was found to be similar to that of ruthenium, and the iron complex furthermore presented a comparable charge transfer in the excited state corresponding to the Q-band, compared to Ru(S2TTP)Cl2. Hence M(S2TTP)Cl2 (M = Mn, Fe, Ni) appear to be potential low cost candidate donor molecules within a bulk heterojunction solar cell. We further propose suitable engineered acceptor pigments, fitted to provide a good overall solar cell efficiency.

Benchmarking DFT Functionals for Excited-State Calculations of Donor-Acceptor TADF Emitters: Insights on the Key Parameters Determining Reverse Inter-System Crossing

The importance of intermediate triplet states and the nature of excited states has gained interest in recent years for the thermally activated delayed fluorescence (TADF) mechanism. It is widely accepted that simple conversion between charge transfer (CT) triplet and singlet excited states is too crude, and a more complex route involving higher-lying locally excited triplet excited states has to be invoked to witness the magnitude of the rate of reverse inter-system crossing (RISC) rates. The increased complexity has challenged the reliability of computational methods to accurately predict the relative energy between excited states as well as their nature. Here, we compare the results of widely used density functional theory (DFT) functionals, CAM-B3LYP, LC-ωPBE, LC-ω*PBE, LC-ω*HPBE, B3LYP, PBE0, and M06-2X, against a wavefunction-based reference method, Spin-Component Scaling second-order approximate Coupled Cluster (SCS-CC2), in 14 known TADF emitters possessing a diversity of chemical structures. Overall, the use of the Tamm-Dancoff Approximation (TDA) together with CAM-B3LYP, M06-2X, and the two ω-tuned range-separated functionals LC-ω*PBE and LC-ω*HPBE demonstrated the best agreement with SCS-CC2 calculations in predicting the absolute energy of the singlet S₁, and triplet T₁ and T₂ excited states and their energy differences. However, consistently across the series and irrespective of the functional or the use of TDA, the nature of T₁ and T₂ is not as accurately captured as compared to S₁. We also investigated the impact of the optimization of S₁ and T₁ excited states on ΔE_ST and the nature of these states for three different functionals (PBE0, CAM-B3LYP, and M06-2X). We observed large changes in ΔE_ST using CAM-B3LYP and PBE0 functionals associated with a large stabilization of T₁ with CAM-B3LYP and a large stabilization of S₁ with PBE0, while ΔE_ST is much less affected considering the M06-2X functional. The nature of the S₁ state barely evolves after geometry optimization essentially because this state is CT by nature for the three functionals tested. However, the prediction of the T₁ nature is more problematic since these functionals for some compounds interpret the nature of T₁ very differently. SCS-CC2 calculations on top of the TDA-DFT optimized geometries lead to a large variation in terms of ΔE_ST and the excited-state nature depending on the chosen functionals, further stressing the large dependence of the excited-state features on the excited-state geometries. The presented work highlights that despite good agreement of energies, the description of the exact nature of the triplet states should be undertaken with caution.

Constructing high-performance TADF polymers from non-TADF monomers: a computational investigation

Thermally activated delayed fluorescence (TADF) polymers excelling in simple, low-cost and large-area solution process ability have attracted tremendous attention recently, but it remains a great challenge for the design of such materials due to the lack of reliable molecular construction guidelines. Here we perform a systematic computational investigation on the construction of TADF polymers from non-TADF monomers to elucidate the effects of polymerization sites, substituent positions and substituent types. The results indicate that the polymerization of 3,6-carbazole-based monomers with different substituents is efficient to build TADF polymers due to their facile π-conjugation extendability. Especially, polymers with para-phenyl-substituted monomers are promising in light of their separated frontier molecular orbitals for small ΔE_ST with favorable energy levels, bipolar charge transport properties and relatively strong absorption/emission intensity, which should be highly attractive for experimental investigations. These findings and insights are important in revealing the structure-property relation of TADF polymers made from non-TADF monomers with important clues for understanding the construction mechanism and molecular design principles of TADF polymers.

Following the density evolution using real time density functional theory and density based indexes: Application to model push-pull molecules

Considering as test case a family of organic rod like push-pull molecules, we derived and applied density based index enabling the description and diagnostic of the electronic density evolution in real time-time dependent density functional theory (RT-TDDFT) simulations. In particular, both the charge transfer (CT) distance and a diagnostic index, the D_CT and M_AC ^RT respectively, were computed on the fly from the density distribution obtained at a given time and the reference ground state density and their mean values were compared with what obtained at Linear Response-TDDFT level. Besides giving a way of analyzing the density redistribution occurring in time, these tools allowed to show how RT-TDDFT, which is definitely a powerful method to model the evolution of the density in CT or charge separation processes, can be affected by the same artifacts known for LR-TDDFT approaches and, particularly, to those related to the use of approximate exchange correlation functionals. The analysis here performed allowed to identify and discard on fly the electronic configurations corresponding to spurious situations.

Modeling the spectral properties of poly(x-phenylenediamine) conducting polymers using a combined TD-DFT and electrostatic embedding approach

The absorption spectra of polymers derived from ortho, meta and para phenylenediamines (o-PDA, m-PDA and p-PDA) have been simulated combining periodic density functional theory (DFT) calculations with time-dependent DFT simulations. These latter have been carried out on finite clusters embedded in a set of point charges devised to exactly reproduce the electrostatic potential of the periodic chains. The results are compared with those obtained for solvated o-PDA, m-PDA and p-PDA oligomers of increasing sizes extracted from the periodic structures. The electronic transitions involved have been investigated by a qualitative analysis based on isodensity maps completed by a quantitative analysis based on the density-based index (D_CT ). For poly-(o)- and poly-(p)- phenylenediamines the agreement with the experimental data is achieved already by modeling solvated dimers whereas the inclusion of long-range electrostatic effects is mandatory for poly-(m)-phenylenediamine highlighting the importance of an accurate treatment of the electrostatic environment when a finite cluster approach is considered.

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.

Unraveling the Symmetry Effects on the Second-Order Nonlinear Optical Responses of Molecular Switches: The Case of Ruthenium Complexes

Owing to their odd order, second-order nonlinear optical (NLO) responses are very sensitive to symmetry. Therefore, within hyper-Rayleigh scattering (HRS) technique, the symmetry impacts the amplitude of the molecular responses, the HRS first hyperpolarizability (β_HRS), and the depolarization ratio (DR). Starting from a challenging octupolar structure bearing six ruthenium(II) ammine centers π-conjugated via quaterpyridyl moieties to a tris-chelated zinc(II) core, together with its Λ shape and one-dimensional analogues built by replacing one or two Ru-quaterpyridyl moieties with bipyridine moieties, (time-dependent) density functional theory calculations have been enacted to unravel the symmetry-NLO response relationships as well as their Ru^II/III redox-triggered switching effects. The one-dimensional and Λ-shaped NLOphores present β_HRS values ∼3 times larger than those of the octupolar system, for both Ru oxidation states. However, using the few-state valence bond-charge transfer models demonstrates that the β_HRS response of the octupolar complex cation can become larger than those of its one-dimensional and Λ-shaped analogues provided stronger donor-acceptor groups are employed. In parallel, the DRs decrease from a strong dipolar character (DR ≈ 6) for the one-dimensional chromophore to a weaker dipolar character (DR ≈ 5) for the Λ-shaped one and to a clear octupolar character (DR ≈ 1.7) for the last one. In all cases, the β responses originate mostly from metal-to-ligand charge transfer excited states, as revealed using a new scheme for analyzing the variations in electron density upon excitation. The Ru^II/III oxidations lead to a strong decrease in the β_HRS responses, which is attributed to the loss of the donor character of the Ru centers and therefore to the reduction of the push-pull π-conjugation. These results demonstrate that the NLO contrast and the NLO switching behavior of these Ru cations are maintained for the different molecular symmetries. Finally, the character of the β responses of the oxidized species, as revealed by the DR values, further evidences a clear evolution from dipolar to octupolar NLOphores.

Assessing Effects of Different π bridges on Properties of Random Benzodithiophene-thienothiophene Donor and Non-fullerene Acceptor Based Active Layer

This report presents the effect of insertion of four different π bridges, furan, thienothiophene, thiophene, and thiazole, into a random benzodithiophene (BDT)-fluorinated-thienothiophene (TT-F) based donor. Starting from a structure of synthesized donor (D)-acceptor (A) random copolymer with 3:1 ratio, we have designed four D-π-A systems with four different π bridges. Structural, optoelectronic, and charge transport/transfer properties of these donors and donor/NDI (NDI = poly[N,N'-bis(2-hexyldecyl)-naphthalene-1,4,5,8-bis(dicarboximide)-2,6-diyl]-alt-5,5'-(2,2'-bithiophene)) blends are investigated using DFT and TD-DFT methodologies. Our results show that the thiazole based TzP1 oligomer has the deepest HOMO value resulting in the highest open circuit voltage among all systems. The maximum absorption wavelengths of π-linked systems are red-shifted compared to the parent molecule. Rates of charge transfer and charge recombination are the highest and smallest in case of the thiazole/NDI blend system. In addition, hole mobilities in thiophene, thienothiophene, and thiazole based systems are larger than in the parent system. The results indicate that the thiazole unit among the four π bridge units is the most suitable for active layer construction.

Dye-sensitized solar cells strike back

Dye-sensitized solar cells (DSCs) are celebrating their 30th birthday and they are attracting a wealth of research efforts aimed at unleashing their full potential. In recent years, DSCs and dye-sensitized photoelectrochemical cells (DSPECs) have experienced a renaissance as the best technology for several niche applications that take advantage of DSCs' unique combination of properties: at low cost, they are composed of non-toxic materials, are colorful, transparent, and very efficient in low light conditions. This review summarizes the advancements in the field over the last decade, encompassing all aspects of the DSC technology: theoretical studies, characterization techniques, materials, applications as solar cells and as drivers for the synthesis of solar fuels, and commercialization efforts from various companies.

Chasing unphysical TD-DFT excited states in transition metal complexes with a simple diagnostic tool

Transition Metal Complexes (TMCs) are known for the rich variety of their excited states showing different nature and degrees of locality. Describing the energies of these excited states with the same degree of accuracy is still problematic when using time-dependent density functional theory in conjunction with the most current density functional approximations. In particular, the presence of unphysically low lying excited states possessing a relevant Charge Transfer (CT) character may significantly affect the spectra computed at such a level of theory and, more relevantly, the interpretation of their photophysical behavior. In this work, we propose an improved version of the M_AC index, recently proposed by the authors and collaborators, as a simple and computationally inexpensive diagnostic tool that can be used for the detection and correction of the unphysically predicted low lying excited states. The analysis, performed on five prototype TMCs, shows that spurious and ghost states can appear in a wide spectral range and that it is difficult to detect them only on the basis of their CT extent. Indeed, both delocalization of the excited state and CT extent are criteria that must be combined, as in the M_AC index, to detect unphysical states.

Deciphering excited state properties utilizing algebraic diagrammatic construction schemes of decreasing order

Excited state properties are difficult to trace back to the common molecular orbital picture when the excited state wavefunction is a linear combination of two or more Slater determinants. Here, a theoretical methodology is introduced based on the algebraic diagrammatic construction scheme for the polarization propagator (ADC(n)) that allows to make this connection and to eventually derive structure-function relationships. The usefulness of this approach is demonstrated by an analysis of the transition dipole moments of the low-lying 1B_3u and 2B_3u states of anthracene and (1,4,5,8)-tetraazaanthracene.

Influence of (de)protonation on the photophysical properties of phenol-substituted diazine chromophores: experimental and theoretical studies

In this contribution, a series of seven push–pull systems has been designed by combining a protonable diazine heterocycle with a deprotonable phenol unit through various π-conjugated linkers (phenylene, thienylene, thienylenevinylene, and phenylenevinylene).