Dielectric Screening Meets Optimally Tuned Density Functionals

Optimal-tuning of range-separated density functionals to describe the optical and photophysical properties of rhodamine B dimers

Setting up computational approaches enabling the correct prediction of the photophysical properties of rhodamine B (RB) derivatives and their aggregates (referred to as dimers here) is of fundamental importance to rationally drive the design of novel systems of applicative relevance, such as artificial light-harvesting nanosystems. Currently, approaches rooted in time-dependent density functional theory (TD-DFT), which are appealing for their relatively low computational cost, nonetheless have limitations in terms of accuracy, especially while considering RB dimeric species. In this work, we investigated the performances of optimally tuned range-separated hybrid functionals for describing the excited states of RB and its H dimer, focusing on dimeric charge-transfer (CT) states. We compared different optimal-tuning (OT) procedures including or not solvent screening. The results show that the properties of CT states, such as their CT extent, brightness and relative energy ordering, are crucially affected by the Hartree-Fock exchange amount dictated by OT, which is in turn driven by the chosen asymptotic behavior. Finally, to understand the dissimilar pictures provided for the CT states by different tuning approaches, we performed an extensive analysis aimed at elucidating how CT states are affected by range-separation parameters. As a result, a simple procedure is finally provided for easily achieving a unique functional to be applied on both monomer and dimers, its current limitations are highlighted and possible perspectives for future development are envisaged.

Nonempirical Prediction of the Length-Dependent Ionization Potential in Molecular Chains

The ionization potential of molecular chains is well-known to be a tunable nanoscale property that exhibits clear quantum confinement effects. State-of-the-art methods can accurately predict the ionization potential in the small molecule limit and in the solid-state limit, but for intermediate, nanosized systems prediction of the evolution of the electronic structure between the two limits is more difficult. Recently, optimal tuning of range-separated hybrid functionals has emerged as a highly accurate method for predicting ionization potentials. This was first achieved for molecules using the ionization potential theorem (IPT) and more recently extended to solid-state systems, based on an ansatz that generalizes the IPT to the removal of charge from a localized Wannier function. Here, we study one-dimensional molecular chains of increasing size, from the monomer limit to the infinite polymer limit using this approach. By comparing our results with other localization-based methods and where available with experiment, we demonstrate that Wannier-localization-based optimal tuning is highly accurate in predicting ionization potentials for any chain length, including the nanoscale regime.

Thiocarbonyl-Bridged N-Heterotriangulenes for Energy Efficient Triplet Photosensitization: A Theoretical Perspective

Structurally-rigid metal-free organic molecules are of high demand for various triplet harvesting applications. However, inefficient intersystem crossing (ISC) due to large singlet-triplet gap (Δ E S - T ${\Delta {E}_{S-T}}$ ) and small spin-orbit coupling (SOC) between lowest excited singlet and triplet often limits their efficiency. Excited electronic states, fluorescence and ISC rates in several thiocarbonyl-bridged N-heterotriangulene ( m ${m}$ S-HTG) with systematically increased thione content (m = ${m=}$ 0-3) are investigated implementing polarization consistent time-dependent optimally-tuned range-separated hybrid. All m ${m}$ S-HTGs are dynamically stable and also thermodynamically feasible to synthesize. Relative energies of several low-lying singlets (S n ${{S}_{n}}$ ) and triplets (T n ${{T}_{n}}$ ), and their excitation nature (i. e.,n π * ${n{\pi }^{^{\ast}}}$ orπ π * ${\pi {\pi }^{^{\ast}}}$ ) and SOC are determined for these m ${m}$ S-HTGs in dichloromethane. Low-energy optical peak displays gradual red-shift with increasing thione content due to relatively smaller electronic gap resulted from greater degree of orbital delocalization. Significantly large SOC due to different orbital-symmetry and heavy-atom effect produces remarkably high ISC rates (k I S C ${{k}_{ISC}}$ ~1012 s-1) for enthalpically favouredS 1 n π * → T 2 ${{S}_{1}\left(n{\pi }^{^{\ast}}\right)\to {T}_{2}}$ (π π * ${\pi {\pi }^{^{\ast}}}$ ) channel in these m ${m}$ S-HTGs, which outcompete radiative fluorescence rates (~108 s-1) even directly from higher lying optically brightπ π * ${\pi {\pi }^{^{\ast}}}$ singlets. Importantly, high energy triplet excitons of ~1.7 eV resulting from such significantly large ISC rates from non-fluorescentS 1 n π * ${{S}_{1}\left(n{\pi }^{^{\ast}}\right)}$ make these thiocarbonylated HTGs ideal candidates for energy efficient triplet harvest including triplet-photosensitization.

Optical Properties of H-Bonded Heterotriangulene Supramolecular Polymers: Charge-Transfer Excitations Matter

H-bonded N-heterotriangulene (NHT) supramolecular polymers offer a nice playground to explore the nature and dynamics of electronic excitations in low-dimensional organic nanostructures. Here, we report on a comprehensive molecular modeling of the excited-state electronic structure and optical properties of model NHT stacks, highlighting the important role of intermolecular charge-transfer (CT) excitations in shaping their optical absorption and emission lineshapes. Most importantly, we show that the coupling between the local and CT excitations, modulated by the electric fields induced by the presence of polar amide groups forming H-bonded arrays along the stacks, significantly increases the resulting hybrid exciton bandwidth. We discuss these findings in the context of the efficient transport of singlet excitons over the μm length scale reported experimentally on individual self-assembled nanofibers with molecular-scale diameter.

Screened optimally tuned range separated hybrid functional for solvated low bandgap molecular systems

The description of electronic properties of low bandgap molecular system is often performed by using density functional theory (DFT) and time dependent (TD) DFT calculations with the optimally tuned range-separated hybrid (OT-RSH) functional, as it contains the necessary ingredients to reliably predict charge transfer excitations. However, the range separating parameter (ω) is system-dependent and its optimization, including the chemical environment, is intricate. Refaely-Abramson et al. demonstrated that the gap renormalization in molecular crystals, a ground state property, can be represented by an OT-RSH functional screened by ɛstatic [Phys. Rev. B 88, 081204(R) (2013)], the zero frequency scalar dielectric constant. In this study, we propose the use of an OT-RSH functional screened by the scalar dielectric constant in the high frequency limit (OT-sRSH), ɛ∞, an appropriate constraint for vertical ionization energies or excitations in a dielectric environment. We have performed calculations for S,N-heteroacene derivatives in tetrahydrofuran and dichloromethane. The "unscreened" OT-RSH functional tends to underestimate experimental ionization potentials (IPs) and optical gaps (Egs) by up to 1.5 and 0.5 eV, respectively. In contrast, OT-sRSH functional calculations underestimate IPs and Egs by only 0.4 and 0.2 eV. We also compared the OT-sRSH results to explicitly solvated OT-RSH functional calculations for oligothiophenes in dioxane, benzene in ammonia, and methylene blue in water. We observe that both the approaches perform similarly for weakly interacting intermolecular systems and deviate for solvent-solute interacting systems, as expected. In conclusion, the OT-sRSH functional can describe molecular systems with environmental polarization effects accurately, a step toward describing realistic molecular systems.

Predicting the Color Polymorphism of ROY from a Time-Dependent Optimally Tuned Screened Range-Separated Hybrid Functional

Polymorphism is a well-known property of molecular crystals, which allows the same molecule to form solids with several crystalline structures that can differ significantly in physical properties. Polymorphs that possess different optical absorption properties in the visible range may exhibit different perceived colors, a phenomenon known as color polymorphism. One striking example of color polymorphism is given by 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, known as ROY for its red-orange-yellow colors. First-principles prediction of color polymorphism may help in polymorph assignment and design but has proven to be challenging. Here, we predict the absorption spectra and simulate the colors of 12 ROY polymorphs using the general, nonempirical method of time-dependent (TD) optimally tuned screened range-separated hybrid (OT-SRSH) functional. For 5 ROY polymorphs with known experimental absorption spectra, we show that the TD-OT-SRSH approach predicts absorption spectra in quantitative agreement with experiment. For all polymorphs, we show that an accurate simulation of the colors is obtained, paving the way to a fully predictive, low-cost calculation of color polymorphism.

Combining Local Range Separation and Local Hybrids: A Step in the Quest for Obtaining Good Energies and Eigenvalues from One Functional

Some of the most successful exchange-correlation approximations in density functional theory are "hybrids", i.e., they rely on combining semilocal density functionals with exact nonlocal Fock exchange. In recent years, two classes of hybrid functionals have emerged as particularly promising: range-separated hybrids on the one hand, and local hybrids on the other hand. These functionals offer the hope to overcome a long-standing "observable dilemma", i.e., the fact that density functionals typically yield either a good description of binding energies, as obtained, e.g., in global and local hybrids, or physically interpretable eigenvalues, as obtained, e.g., in optimally tuned range-separated hybrids. Obtaining both of these characteristics from one and the same functional with the same set of parameters has been a long-standing challenge. We here discuss combining the concepts of local range separation and local hybrids as part of a constraint-guided quest for functionals that overcome the observable dilemma.

Nature and energetics of low-lying excited singlets/triplets and intersystem crossing rates in selone analogs of perylenediimide: A theoretical perspective

The structural rigidity and chemical diversity of the highly fluorescent perylenediimide (PDI) provide wide opportunities for developing triplet photosensitizers with sufficiently increased energy efficiency. Remarkably high intersystem crossing (ISC) rates with a complete fluorescence turn-off reported recently for several thione analogs of PDI due to substantially large spin-orbit coupling garners huge attention to develop other potential analogs. Here, several selone analogs of PDI, denoted as mSe-PDIs (m = 1-4) with varied Se content and positions, are investigated to provide a comprehensive and comparative picture down the group-16 using density functional theory (DFT) and time-dependent DFT implementing optimally tuned range-separated hybrid in toluene dielectric. All mSe-PDIs are confirmed to be dynamically stable and also thermodynamically feasible to synthesize from their oxygen and thione congeners. The first excited-state singlet (S1) of mSe-PDI with relatively low Se-content (m = 1, 2) is of nπ* character with an expected fluorescence turn-off. Whereas, the ππ* nature of the S1 for 3Se-PDI and 4Se-PDI suggests a possible fluorescence turn-on in the absence of any other active nonradiative deactivation pathways. However, ∼4-6 orders greater ISC rates (∼1012-1014 s-1) than the fluorescence ones (∼108 s-1) for all mSe-PDIs signify highly efficient triplet harvest. Importantly, significantly higher ISC rates for these mSe-PDIs than their thione congeners render them efficient triplet photosensitizers.

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.

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

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

Does the Intersystem Crossing Rate of β-Iodinated Phosphorus Corrole Depend on Iodine Numbers and/or Positions?

The heavy-atom effect is known to enhance the intersystem crossing (ISC) in organic molecular systems. Effects of iodine numbers and positions on the ISC rate of a few meso-difluorophenyl substituted β-iodinated phosphorus corroles (PCs) with axially ligated fluorine atoms (mI-FPC; m = 1-4) are studied using a time-dependent optimally tuned range-separated hybrid. Solvent effects are accounted for through a polarizable continuum model with a toluene dielectric. Calculations suggest similar thermodynamic stability for all mI-FPCs and also reproduce the experimentally measured 0-0 energies for some of the freebase phosphorus corrole (FPC) systems studied here. Importantly, our results reveal that all mI-FPCs display 10 times larger ISC rate (∼109 s-1) than the fluorescence rate (∼108 s-1), and the higher ISC rate stems from the improved spin-orbit coupling (SOC) introduced by lighter heteroatoms like central P and biaxial F rather than the I heavy-atom effect. However, an enhanced SOC is found with increasing I content for El-Sayed forbidden ISC channels. Research findings reported in this study unveil the impact of light heteroatoms and heavy atoms in promoting ISC in several iodinated PCs, which help in designing visible-light-driven efficient triplet photosensitizers.

Do Optimally Tuned Range-Separated Hybrid Functionals Require a Reparametrization of the Dispersion Correction? It Depends

For ground- and excited-state studies of large molecules, it is the state of the art to combine (time-dependent) DFT with dispersion-corrected range-separated hybrid functionals (RSHs), which ensures an asymptotically correct description of exchange effects and London dispersion. Specifically for studying excited states, it is common practice to tune the range-separation parameter ω (optimal tuning), which can further improve the accuracy. However, since optimal tuning essentially changes the functional, it is unclear if and how much the parameters used for the dispersion correction depend on the chosen ω value. To answer this question, we explore this interdependency by refitting the DFT-D4 dispersion model for six established RSHs over a wide range of ω values (0.05-0.45 a0-1) using a set of noncovalently bound molecular complexes. The results reveal some surprising differences among the investigated functionals: While PBE-based RSHs and ωB97M-D4 generally exhibit a weak interdependency and robust performance over a wide range of ω values, B88-based RSHs, specifically LC-BLYP, are strongly affected. For these, even a minor reduction of ω from the default value manifests in strong systematic overbinding and poor performance in the typical range of optimally tuned ω values. Finally, we discuss strategies to mitigate these issues and reflect the results in the context of the employed D4 parameter optimization algorithm and fit set, outlining strategies for future improvements.

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

Solvent effects on the photoinduced charge separation dynamics of directly linked zinc phthalocyanine-perylenediimide dyads: a nonadiabatic dynamics simulation with an optimally tuned screened range-separated hybrid functional

Herein, we have employed a combination of the optimally tuned screened range-separated hybrid (OT-SRSH) functional, the polarizable continuum model (PCM), and nonadiabatic dynamics (NAMD) simulations to investigate the photoinduced dynamics of directly linked donor-acceptor dyads formed using zinc phthalocyanine (ZnPc) and perylenediimide (PDI), in which ZnPc is the donor while PDI is the acceptor. Our simulations aim to analyze the behavior of these dyads upon local excitation of the ZnPc moiety in the gas phase and in benzonitrile. Our findings indicate that the presence of a solvent can significantly influence the excited state dynamics of ZnPc-PDI dyads. Specifically, the polar solvent benzonitrile effectively lowers the vertical excitation energies of the charge transfer (CT) state from ZnPc to PDI. As a result, the energetic order of the locally excited (LE) states of ZnPc and the CT states is reversed compared to the gas phase. Consequently, the photoinduced electron transfer (PET) dynamics from ZnPc to PDI, which is absent in the gas phase, takes place in benzonitrile with a time constant of 10.4 ps. Importantly, our present work not only qualitatively agrees with experimental results but also provides in-depth insights into the underlying mechanisms responsible for the photoinduced dynamics of ZnPc-PDI. Moreover, this study emphasizes the importance of appropriately considering solvent effects in NAMD simulation of organic donor-acceptor systems, taking into account the distinct excited state dynamics observed in the gas phase and benzonitrile. Furthermore, the combination of the OT-SRSH functional, the PCM solvent model, and nonadiabatic dynamics simulations shows promise as a strategy for investigating the complex excited state dynamics of organic donor-acceptor systems in solvents. These findings will be valuable for the future design of novel organic donor-acceptor structures with improved performance.

Theoretical Study of Structure and Photophysics of Homologous Series of Bis(arylydene)cycloalkanones

Photophysical properties of a series of bis(arylydene)cycloalkanone dyes with various donor substituents are studied using quantum chemistry. Their capacity for luminescence and nonradiative relaxation through trans-cis isomerization is related to their structure, in particular, to the donor capacity of the substituents and the degree of conjugation due to the central cycloalkanone moiety. It is shown that cyclohexanone central moiety introduces distortions and disrupts the conjugation, thus leading to a nonmonotonic change in their properties. The increasing donor capacity of the substituents causes increase in the HOMO energy (rise in the oxidation potential) and decrease in the HOMO-LUMO gap, which results in the red shift of the absorption spectra. The ability of the excited dye to relax through fluorescence or through trans-cis isomerization is governed by the height of the barrier between the Franck-Condon and S1-S0 conical intersection regions on the potential energy surface of the lowest π-π* excited state. This barrier also correlates with the donor capacity of the substituents and the degree of conjugation between the central and donor moieties. The calculated fluorescence and trans-cis isomerization rates are in good agreement with the observed fluorescence quantum yields.

Embedding vertex corrections in GW self-energy: Theory, implementation, and outlook

The vertex function (Γ) within the Green's function formalism encapsulates information about all higher-order electron-electron interaction beyond those mediated by density fluctuations. Herein, we present an efficient approach that embeds vertex corrections in the one-shot GW correlation self-energy for isolated and periodic systems. The vertex-corrected self-energy is constructed through the proposed separation-propagation-recombination procedure: the electronic Hilbert space is separated into an active space and its orthogonal complement denoted as the "rest;" the active component is propagated by a space-specific effective Hamiltonian different from the rest. The vertex corrections are introduced by a rescaled time-dependent nonlocal exchange interaction. The direct Γ correction to the self-energy is further updated by adjusting the rescaling factor in a self-consistent post-processing cycle. Our embedding method is tested mainly on donor-acceptor charge-transfer systems. The embedded vertex effects consistently and significantly correct the quasiparticle energies of the gap-edge states. The fundamental gap is generally improved by 1-3 eV upon the one-shot GW approximation. Furthermore, we provide an outlook for applications of (embedded) vertex corrections in calculations of extended solids.

Understanding Primary Charge Separation in the Heliobacterial Reaction Center

The homodimeric reaction center of heliobacteria retains features of the ancestral reaction center and can thus provide insights into the evolution of photosynthesis. Primary charge separation is expected to proceed in a two-step mechanism along either of the two reaction center branches. We reveal the first charge-separation step from first-principles calculations based on time-dependent density functional theory with an optimally tuned range-separated hybrid and ab initio Born-Oppenheimer molecular dynamics: the electron is most likely localized on the electron transfer cofactor 3 (EC3, OH-chlorophyll a), and the hole on the adjacent EC2. Including substantial parts of the surrounding protein environment into the calculations shows that a distinct structural mechanism is decisive for the relative energetic positioning of the electronic excitations: specific charged amino acids in the vicinity of EC3 lower the energy of charge-transfer excitations and thus facilitate efficient charge separation. These results are discussed considering recent experimental insights.

Either Accurate Singlet-Triplet Gaps or Excited-State Structures: Testing and Understanding the Performance of TD-DFT for TADF Emitters

The energy gap between the lowest singlet and triplet excited states (ΔE_ST) is a key property of thermally activated delayed fluorescence (TADF) emitters, where these states are dominated by charge-transfer (CT) character. Despite its well-known shortcomings concerning CT states, time-dependent density functional theory (TD-DFT) is widely used to predict this gap and study TADF. Moreover, polar CT states exhibit a strong interaction with their molecular environment, which further complicates their computational description. Addressing these two major challenges, this work studies the performance of Tamm-Dancoff-approximated TD-DFT (TDA-DFT) on the recent STGABS27 benchmark set,¹ exploring different strategies to include orbital and structural relaxation, as well as dielectric embedding. The results show that the best-performing strategy is to calculate ΔE_ST at the ground-state structure using functionals with a surprisingly small amount of Fock exchange of ≈10% and without a (complete) solvent model. However, as this approach heavily relies on error cancellation to mimic dielectric relaxation, it is not robust and exhibits large systematic deviations in excited state energies, state characters, and structures. More rigorous approaches, including state-specific solvation, do not share these systematic deviations, but their predicted ΔE_ST values exhibit larger statistical errors. We thus conclude that for the description of CT states in dielectric environments, none of the tested TDA-DFT methods is competitive with the recently presented ROKS/PCM approach regarding robustness, accuracy, and computational efficiency.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Decomposition analysis on the excitation behaviors of thiazolothiazole (TTz)-based dyes via the time-dependent dielectric density functional theory approach

Thiazolothiazole (TTz)-based materials have been attracting much attention because of their widespread applications. In this paper, we discuss the excited electronic behaviors of asymmetric TTz dyes in solvents based on the time-dependent dielectric density functional theory method. Based on dipole moment and charge distribution (population) analyses, we discuss large intramolecular electron transfers, which are triggered by photon excitations, toward the acceptor part of dyes. In addition, we explore the contributions of geometrical changes and solvent reorientations (reorganizations) to the solvatofluorochromic phenomena based on a decomposition technique. The decomposition analysis shows that the solvent reorientation effect mainly contributes to changes in the fluorescent spectra in response to solvents.

Toward highly efficient hyperfluorescence-based emitters through excited-states alignment using novel optimally tuned range-separated models

Hyperfluorescence has recently been introduced as a promising strategy to achieve organic light-emitting diodes (OLEDs) with high color purity and enhanced stability. In this approach, fluorescent emitters (FEs) with strong and narrow band fluorescence are integrated in thin films containing sensitizers exhibiting thermally activated delayed fluorescence (TADF). Toward highly efficient hyperfluorescence-based emitters, the excited-states ordering of the FEs should be well-aligned. Given some recent endeavors in this context, the related theoretical explorations are relatively limited and have proven to be challenging. In this work, alignments of the corresponding excited-states, crucial for both the fast Förster resonance energy transfer and suppression of the Dexter energy transfer from TADF sensitizers to FEs, have theoretically been investigated using optimally tuned range-separated hybrid functionals (OT-RSHs). We have proposed and validated several variants of the models including OT-RSHs, their coupled versions with the polarizable continuum model, OT-RSHs-PCM, as well as the screened versions accounting for the screening effects by the electron correlation through the scalar dielectric constant, OT-SRSHs, for a reliable description of the excited-states ordering in the FEs of the hyperfluorescence-based materials. Particular attention is paid to the influence of the underlying density functional approximations as well as the short- and long-range Hartree-Fock (HF) exchange contributions and the range-separation parameter. Considering a series of experimentally known hyperfluorescence-based emitters as working models, it is unveiled that any combination of the ingredients in the proposed models does not render the correct order of the excited-states of the FEs, but a particular compromise among the involved parameters is needed to more accurately account for the relevant excited-states alignment. Perusing the results of our developed methods, the best ones are found to be the generalized gradient approximation-based OT-RSHs-PCM with the correct asymptotic behavior and incorporating no (low) HF exchange contribution at the short-range regime. The proposed models show superior performances not only with respect to their standard counterparts with the default parameters but also as compared to other range-separated approximations. Accountability of the best-proposed model is also put into broader perspective, where it has been employed for the computational design of several molecules as promising FE candidates prone to be utilized in hyperfluorescence-based materials. Summing up, the proposed models in this study can be recommended for both the theoretical modeling and confirming the experimental observations in the field of hyperfluorescence-based OLEDs.

Screening the Band Shape of Molecules by Optimal Tuning of Range-Separated Hybrid Functional with TD-DFT: A Molecular Designing Approach

In the present article we have demonstrated the effectiveness of optimally tuned range-separated hybrid (RSH) functional to determine the electronic transitions from two fluorophore moieties (blue and yellow/orange) within a single white light emitter (WLE). First, the optimally tuned range separation parameter (ω) is calculated for two white emitting fluorophores (W1 and W2) already reported in the literature. The success of the optimally tuned RSH functional ω*B97XD, used in the TDDFT study, encouraged the authors to design eight new single organic white light emitters with frustrated energy transfer between the two individual fluorophore moieties (blue and yellow/orange). The simulated spectra (the band shapes, to be more specific) generated by TDDFT study and outcomes through natural transition orbital (NTO) and natural bond orbital (NBO) studies clearly demonstrate that all the designed eight organic molecules are potential white light emitters and can be synthesized in future.

Recent Advances in Cartesian-Grid DFT in Atoms and Molecules

In the past several decades, density functional theory (DFT) has evolved as a leading player across a dazzling variety of fields, from organic chemistry to condensed matter physics. The simple conceptual framework and computational elegance are the underlying driver for this. This article reviews some of the recent developments that have taken place in our laboratory in the past 5 years. Efforts are made to validate a viable alternative for DFT calculations for small to medium systems through a Cartesian coordinate grid- (CCG-) based pseudopotential Kohn-Sham (KS) DFT framework using LCAO-MO ansatz. In order to legitimize its suitability and efficacy, at first, electric response properties, such as dipole moment ( μ ), static dipole polarizability ( α ), and first hyperpolarizability ( β ), are calculated. Next, we present a purely numerical approach in CCG for proficient computation of exact exchange density contribution in certain types of orbital-dependent density functionals. A Fourier convolution theorem combined with a range-separated Coulomb interaction kernel is invoked. This takes motivation from a semi-numerical algorithm, where the rate-deciding factor is the evaluation of electrostatic potential. Its success further leads to a systematic self-consistent approach from first principles, which is desirable in the development of optimally tuned range-separated hybrid and hyper functionals. Next, we discuss a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, by means of "adiabatic connection theorem" and virial theorem. Optical gaps in organic chromophores, dyes, linear/non-linear PAHs, and charge transfer complexes are faithfully reproduced. In short, CCG-DFT is shown to be a successful route for various practical applications in electronic systems.

Excited-State Properties of Some Thermally Activated Delayed Fluorescence Emitters: Quest for an Accurate and Reliable Computational Method

Thermally activated delayed fluorescence (TADF) finds application in organic light-emitting diodes. The molecules exhibiting TADF are characterized by small singlet-triplet energy gaps that help reverse intersystem crossing. Recently, ionization potential (IP)-tuned range-separated (RS) density functionals have been well accepted for studying excited-state properties. In the present work, two efficient descriptor-based tuning schemes [electron localization function (ELF) and Sol] of RS density functionals have been used to accurately reproduce the excited-state properties of TADF emitters by performing a single self-consistent field calculation. The lowest singlet vertical excitation energies (E_VA(S₁)) and the vertical singlet-triplet energy gaps (ΔE_VST) are computed with ELF-, Sol-, and IP-tuned RS functionals (LC-BLYP, ωB97, ωB97X, and ωB97XD). Encouraging mean absolute deviations from the experimental values with ELF*-, Sol*-, and IP-tuned functionals are observed. Consistent performance of the non-empirical tuned functionals is noted in different solvent dielectrics. In addition to these, fractional occupation calculations have shown that our tuned functionals almost satisfy the energy linearity curve. Thus, ELF*- and Sol*-tuned functionals are promising and reliable alternatives in computing the excited-state properties. Considering the small experimental singlet-triplet gap, we recommend ELF* to calculate E_VA(S₁) and Sol* to calculate ΔE_VST.

Accurate Polarization-Resolved Absorption Spectra of Organic Semiconductor Thin Films Using First-Principles Quantum-Chemical Methods: Pentacene as a Case Study

Theoretical studies using clusters as model systems have been extremely successful in explaining various photophysical phenomena in organic semiconductor (OSC) thin films. But they have not been able to satisfactorily simulate total and polarization-resolved absorption spectra of OSCs so far. In this work, we demonstrate that accurate spectra are predicted by time-dependent density functional theory (TD-DFT) when the employed cluster reflects the symmetry of the crystal structure and all monomers feel the same environment. Additionally, long-range corrected optimal tuned functionals are mandatory. For pentacene thin films, the computed electronic spectra for thin films then reach an impressive accuracy compared with experimental data with a deviation of less than 0.1 eV. This allows for accurate peak assignments and mechanistic studies, which paves the way for a comprehensive understanding of OSCs using an affordable and easy-to-use cluster approach.

Density functionals with asymptotic-potential corrections are required for the simulation of spectroscopic properties of materials

Five effects of correction of the asymptotic potential error in density functionals are identified that significantly improve calculated properties of molecular excited states involving charge-transfer character. Newly developed materials-science computational methods are used to demonstrate how these effects manifest in materials spectroscopy. Connection is made considering chlorophyll-a as a paradigm for molecular spectroscopy, 22 iconic materials as paradigms for 3D materials spectroscopy, and the VN - defect in hexagonal boron nitride as an example of the spectroscopy of defects in 2D materials pertaining to nanophotonics. Defects can equally be thought of as being "molecular" and "materials" in nature and hence bridge the relms of molecular and materials spectroscopies. It is concluded that the density functional HSE06, currently considered as the standard for accurate calculations of materials spectroscopy, should be replaced, in most instances, by the computationally similar but asymptotically corrected CAM-B3LYP functional, with some specific functionals for materials-use only providing further improvements.

Linear fractional charge behavior in density functional theory through dielectric tuning of conductor-like polarizable continuum model

A property of exact density functional theory is linear fractional charge behavior as electrons are added or removed from a molecule. Typical density functional approximations (DFAs) exhibit delocalization error, which overstabilizes this fractional charge. Conversely, solvent corrections have been shown to erroneously destabilize this fractional charge. This work will show that an implicit solvent correction with a tuned dielectric can be used as an ad hoc correction to offset the delocalizing character of DFAs and achieve linear fractional charge behavior. While desirable, in principle, we find that this linear charge behavior degrades the vertical ionization energies reported by DFAs. Our results reveal that the localizing character of the solvent correction and the Hartree-Fock (HF) exchange offset each other. This helps explain the decreased ratios of HF exchange to DFA exchange in long-range hybrid tuning studies that use a solvent correction.

Ultrafast charge transfer in a nonfullerene all-small-molecule organic solar cell: a nonadiabatic dynamics simulation with optimally tuned range-separated functional

The combination of nonadiabatic dynamics simulation and optimally tuned range-separated functional might be a powerful tool for elucidating the ultrafast charge transfer in nonfullerene all-small-molecule organic solar cells.

A theoretical study on solvatofluorochromic asymmetric thiazolothiazole (TTz) dyes using dielectric-dependent density functional theory

In this work, the excitation energies of asymmetric thiazolothizaole (TTz) dye molecules have been theoretically studied using dielectric-dependent density functional theory (DFT). In the dielectric-dependent DFT approach, the ratio (fraction) of the nonlocal Hartree exchange term incorporated into the DFT exchange-correlation functional is a system-dependent parameter, which is inversely proportional to the dielectric constant of the target material. The dielectric-dependent DFT method is closely related to the Coulomb hole and screened exchange (COHSEX) approximation in the GW method and therefore has been applied to crystalline systems with periodic boundary conditions, such as semiconductors and inorganic materials. By focusing on the solvatofluorochromic phenomena of asymmetric TTz dyes, we show that excitation energy calculations obtained from the dielectric-dependent DFT method can reproduce the corresponding experimental UV-vis absorption and emission spectra of dyes in solvents.

Band gaps of crystalline solids from Wannier-localization-based optimal tuning of a screened range-separated hybrid functional

Accurate prediction of fundamental band gaps of crystalline solid-state systems entirely within density functional theory is a long-standing challenge. Here, we present a simple and inexpensive method that achieves this by means of nonempirical optimal tuning of the parameters of a screened range-separated hybrid functional. The tuning involves the enforcement of an ansatz that generalizes the ionization potential theorem to the removal of an electron from an occupied state described by a localized Wannier function in a modestly sized supercell calculation. The method is benchmarked against experiment for a set of systems ranging from narrow band-gap semiconductors to large band-gap insulators, spanning a range of fundamental band gaps from 0.2 to 14.2 electronvolts (eV), and is found to yield quantitative accuracy across the board, with a mean absolute error of ∼0.1 eV and a maximal error of ∼0.2 eV.

Unified Framework for Photophysical Rate Calculations in TADF Molecules

One of the challenges in organic light-emitting diodes research is finding ways to increase device efficiency by making use of the triplet excitons that are inevitably generated in the process of electroluminescence. One way to do so is by thermally activated delayed fluorescence (TADF), a process in which triplet excitons undergo upconversion to singlet states, allowing them to relax radiatively. The discovery of this phenomenon has ensued a quest for new materials that are able to effectively take advantage of this mechanism. From a theoretical standpoint, this requires the capacity to estimate the rates of the various processes involved in the photophysics of candidate molecules, such as intersystem crossing, reverse intersystem crossing, fluorescence, and phosphorescence. Here, we present a method that is able to, within a single framework, compute all of these rates and predict the photophysics of new molecules. We apply the method to two TADF molecules and show that results compare favorably with other theoretical approaches and experimental results. Finally, we use a kinetic model to show how the calculated rates act in concert to produce different photophysical behavior.

Efficient treatment of molecular excitations in the liquid phase environment via stochastic many-body theory

Accurate predictions of charge excitation energies of molecules in the disordered condensed phase are central to the chemical reactivity, stability, and optoelectronic properties of molecules and critically depend on the specific environment. Herein, we develop a stochastic GW method for calculating these charge excitation energies. The approach employs maximally localized electronic states to define the electronic subspace of a molecule and the rest of the system, both of which are randomly sampled. We test the method on three solute-solvent systems: phenol, thymine, and phenylalanine in water. The results are in excellent agreement with the previous high-level calculations and available experimental data. The stochastic calculations for supercells containing up to 1000 electrons representing the solvated systems are inexpensive and require ≤1000 central processing unit hrs. We find that the coupling with the environment accounts for ∼40% of the total correlation energy. The solvent-to-solute feedback mechanism incorporated in the molecular correlation term causes up to 0.6 eV destabilization of the quasiparticle energy. Simulated photo-emission spectra exhibit red shifts, state-degeneracy lifting, and lifetime shortening. Our method provides an efficient approach for an accurate study of excitations of large molecules in realistic condensed phase environments.

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low-energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. Intensity borrowing by these spurious states can affect intensities of the valence excitations, altering electronic bandshapes. At higher excitation energies, it is difficult to distinguish spurious charge-transfer states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces "natural charge-transfer orbitals" that provide a means to isolate orbitals that are most likely to participate in a CTTS excitation. Projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation and be used as an automated subspace selection algorithm for projection-based embedding of a high-level description of the CTTS state in a lower-level description of its environment. We apply this method to an ab initio molecular dynamics trajectory of I^-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with the experiment, and only one-third of the water molecules in the I^-(H₂O)₉₆ simulation cell need to be described with LR-TDDFT to obtain excitation energies that are converged to <0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments.

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso-meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4-phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi's golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

From Kohn-Sham to Many-Electron Energies via Step Structures in the Exchange-Correlation Potential

Accurately describing excited states within Kohn–Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy (Δ), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.

First principles theoretical spectroscopy of methylene blue: Between limitations of time-dependent density functional theory approximations and its realistic description in the solvent

Methylene blue [3,7-Bis(di-methylamino) phenothiazin-5-ium chloride] is a phenothiazine dye with applications as a sensitizer for photodynamic therapy, photoantimicrobials, and dye-sensitized solar cells. Time-dependent density functional theory (TDDFT), based on (semi)local and global hybrid exchange-correlation functionals, fails to correctly describe its spectral features due to known limitations for describing optical excitations of π-conjugated systems. Here, we use TDDFT with a non-empirical optimally tuned range-separated hybrid functional to explore the optical excitations of gas phase and solvated methylene blue. We compute solvated configurations using molecular dynamics and an iterative procedure to account for explicit solute polarization. We rationalize and validate that by extrapolating the optimized range separation parameter to an infinite amount of solvating molecules, the optical gap of methylene blue is well described. Moreover, this method allows us to resolve contributions from solvent-solute intermolecular interactions and dielectric screening. We validate our results by comparing them to first-principles calculations based on the GW+Bethe-Salpeter equation approach and experiment. Vibronic calculations using TDDFT and the generating function method account for the spectra's subbands and bring the computed transition energies to within 0.15 eV of the experimental data. This methodology is expected to perform equivalently well for describing solvated spectra of π-conjugated systems.

Morphology dependent interaction between Co(II)-tetraphenylporphyrin and the MgO(100) surface

Porphyrins are key elements in organic–inorganic hybrid systems for a wide range of applications. Understanding their interaction with the substrate gives a handle on structural and electronic device properties. Here we investigate a single transition-metal porphyrin, namely Co(ii)-tetraphenylporphyrin (CoTPP), on the MgO(100) surface and the effect of multilayer film formation within hybrid density-functional theory and many-body perturbation theory. We focus on the relevant adsorption sites, simulate their photoemission spectra as a key fingerprint and compare with experiments on MgO(100) films on Ag(100). While we find only weak interaction between the cobalt centre and terrace sites on the MgO(100) surface, a strong interaction manifests itself with the low-coordinated sites. This leads to distinct features in both the valence and core-level regions of the electronic structure, as observed in the ultraviolet and X-ray photoemission spectra, corroborated by simulated spectra and calculated cobalt core-level shifts. Our work thus demonstrates the relevance of morphology-related low-coordinated sites and their properties in the adsorption of CoTPP on the MgO(100) surface.

The adsorption of Co-tetraphenylporphyrin at relevant low-coordinated sites on MgO(100) shows distinct features from terrace-site and multilayer films in the near-valence and corelevel regions of the electronic structure.

Charge-transfer electronic states in organic solar cells: a TDDFT study

The prediction of new organic photovoltaic materials in organic solar cells (OSCs) must include a precise description of charge-transfer states because they are involved in electron-transfer processes such as charge separation and charge recombination which govern the device efficiency. Also, as the experimental performance of an optoelectronic device is measured for nonequilibrium nanostructures, computational approaches need models that can incorporate morphology effects. Usually, this aspect is treated by molecular dynamics simulation (MDS) methodologies; however, methodologies and formalisms to calculate the electron-transfer processes are still controversial and sometimes do not connect their information with the phase morphologies. In this work we propose a simple and fast characterization of electron-transfer processes to find the rate constants by analysing the distribution of vertical excitation energies of both local excitation (LE) and charge-transfer (CT) states using TD-DFT calculations in the donor-acceptor pair structures which were extracted from MDS. This proposal assumes that conformational changes are prevented and equilibria are not achieved while the electron-transfer events take effect, and thus the only pathway that connects the LE and CT states is their surface crossing point where an ideal distribution might exist. Different density functionals and dialectric models were tested. The results indicate a close relationship between the proposal and experimental data for electron-transfer events, suggesting the application of this method in the rational design of new photovoltaic materials.

Singlet fission relevant energetics from optimally tuned range-separated hybrids

As a promising idea to design high-efficiency organic photovoltaics, singlet fission (SF) mechanism, i.e., generating two triplet excitons out of a single photon absorption, has recently come into the spotlight. Even though much effort has been devoted to this arena, accurately accounting for the SF process from the theoretical perspective has proven to be challenging. Herein, the SF energetics have thoroughly been investigated with the help of optimally tuned range-separated hybrid functionals (OT-RSHs) in both gas and solvent phases. Taking a series of experimentally known SF chromophores as working models, we have proposed and validated several variants of OT-RSH approximations for the reliable prediction of the energy levels which match the crucial criteria for the SF process, namely, the negative singlet-triplet and triplet-triplet energy gaps. We scrutinize the role of the OT-RSH ingredients, i.e., the underlying density functional approximations, short- and long-range exact-like exchange, as well as the range-separation parameter, for our purpose. The newly designed OT-RSHs outperform the standard RSHs and other related schemes such as screened-exchange approximations as well as other density functionals from different rungs for describing the SF energetics. More importantly, it is unveiled that although the OT-RSH coupled with the polarizable continuum model, OT-RSH-PCM, as well as the screened versions, OT-SRSHs, which account for the screening effect by the electron correlation through the scalar dielectric constant have some advantages over gas-phase computations using OT-RSHs, the energetics criteria of the SF process may not necessarily be satisfied. This in turn corroborates the idea of performing both the optimal tuning procedure and subsequent computations of the SF relevant energetics using OT-RSHs as a more reliable and affordable framework, at least for the present purpose. The applicability of the proposed models is also put into broader perspective, where they are used for the computational design of several chromophores as promising candidates prone to utilization in the SF-based materials. Hopefully, our recommended OT-RSHs can function as efficient models for both the theoretical modeling of SF chromophores and confirming the experimental observations in the field.

Molecular excitations from meta-generalized gradient approximations in the Kohn-Sham scheme

Meta-Generalized Gradient Approximations (meta-GGAs) can, in principle, include spatial and temporal nonlocality in time-dependent density functional theory at a much lower computational cost than functionals that use exact exchange. We here test whether a meta-GGA that has recently been developed with a focus on capturing nonlocal response properties and the particle number discontinuity can realize such features in practice. To this end, we extended the frequency-dependent Sternheimer formalism to the meta-GGA case. Using the Krieger-Li-Iafrate (KLI) approximation, we calculate the optical response for the selected paradigm molecular systems and compare the meta-GGA Kohn-Sham response to the one found with exact exchange and conventional (semi-)local functionals. We find that the new meta-GGA captures important properties of the nonlocal exchange response. The KLI approximation, however, emerges as a limiting factor in the evaluation of charge-transfer excitations.

Computational screen-out strategy for electrically pumped organic laser materials

Electrically pumped organic lasing is one of the most challenging issues in organic optoelectronics. We present a systematic theoretical investigation to screen out electrical pumping lasing molecules over a wide range of organic materials. With the electronic structure information obtained from time-dependent density functional theory, we calculate multiple photophysical parameters of a set of optical pumping organic laser molecules in our self-developed molecular material property prediction package (MOMAP) to judge whether the electrically pumped lasing conditions can be satisfied, namely, to avoid reabsorption from excitons and/or polarons, and the accumulation of triplet excitons. In addition, a large oscillator strength of S₁ and weak intermolecular π–π interaction are preferred. With these criteria, we are able to conclude that BP3T, BSBCz, and CzPVSBF compounds are promising candidates for electrically pumped lasing, and the proposed computational strategy could serve as a general protocol for molecular design of organic lasing materials.

Though the goal of current organic solid-state laser research remains the realization of electrically pumped lasing, identifying organic semiconductors with ideal properties remains a challenge. Here, the authors report a computational strategy for screening electrical pumping lasing molecules.