Screened Range-Separated Hybrid Functional with Polarizable Continuum Model Overcomes Challenges in Describing Triplet Excitations in the Condensed Phase Using TDDFT

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.

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.

Room Temperature Phosphorescence in Crystalline Iodinated Eumelanin Monomer

We report the room temperature phosphorescence upon iodination on a crystalline eumelanin monomer with shielded hydroxyl moieties, ethyl 5,6-dimethoxyindole-2-carboxylate (DMICE). Ultrafast intersystem crossing (ISC) is observed in the iodinated (IDMICE) as well as brominated (BDMICE) analogues of the eumelanin monomer derivative in solution. The triplet quantum yields (φT) and intersystem crossing rates (kISC) of the halogenated eumelanin derivatives areφ T B D M I C E ${{\phi{} }_{T}^{BDMICE}}$ =25.4±1.1 %;k I S C B D M I C E ${{k}_{ISC}^{BDMICE}}$ =1.95×109 s-1 andφ T I D M I C E ${{\phi{} }_{T}^{IDMICE}}$ =59.1±1.6 %;k I S C I D M I C E = ${{k}_{ISC}^{IDMICE}=}$ 1.36×1010 s-1, as monitored using transient absorption spectroscopy. Theoretical calculations based on nuclear ensemble method reveal that computed kISC and spin-orbit coupling matrix elements for eumelanin derivatives are larger for IDMICE relative to BDMICE. The halogen and π-π interactions, with distinct excitonic coupling and higher ISC rate promote phosphorescence in IDMICE molecular crystals. Accessing triplet excited states and resultant photoluminescence through structural modification of eumelanin scaffolds paves way for exploring the versatility of eumelanin-inspired molecules as bio-functional materials.

Antioxidative Triplet Excitation Energy Transfer in Bacterial Reaction Center Using a Screened Range Separated Hybrid Functional

Excess energy absorbed by photosystems (PSs) can result in photoinduced oxidative damage. Transfer of such energy within the core pigments of the reaction center in the form of triplet excitation is important in regulating and preserving the functionality of PSs. In the bacterial reaction center (BRC), the special pair (P) is understood to act as the electron donor in a photoinduced charge transfer process, triggering the charge separation process through the photoactive branch A pigments that experience a higher polarizing environment. At this work, triplet excitation energy transfer (TEET) in BRC is studied using a computational perspective to gain insights into the roles of the dielectric environment and interpigment orientations. We find in agreement with experimental observations that TEET proceeds through branch B. The TEET process toward branch B pigment is found to be significantly faster than the hypothetical process proceeding through branch A pigments with ps and ms time scales, respectively. Our calculations find that conformational differences play a major role in this branch asymmetry in TEET, where the dielectric environment asymmetry plays only a secondary role in directing the TEET to proceed through branch B. We also address TEET processes asserting the role of carotenoid as the final triplet energy acceptor and in a mutant form, where the branch pigments adjacent to P are replaced by bacteriopheophytins. The necessary electronic excitation energies and electronic state couplings are calculated by the recently developed polarization-consistent framework combining a screened range-separated hybrid functional and a polarizable continuum mode. The polarization-consistent potential energy surfaces are used to parametrize the quantum mechanical approach, implementing Fermi's golden rule expression of the TEET rate calculations.

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.

A Computational Study of the Electronic Energy and Charge Transfer Rates and Pathways in the Tetraphenyldibenzoperiflanthene/Fullerene Interfacial Dyad

The electronic transition rates and pathways underlying interfacial charge separation in tetraphenyldibenzoperiflanthene:fullerene (DBP:C70) blends are investigated computationally. The analysis is based on a polarization-consistent framework employing screened range-separated hybrid functional in a polarizable continuum model to parametrize Fermi's golden rule rate theory. The model considers the possible transitions within the 25 lowest excited states of a DBP:C70 dyad that are accessible by photoexcitation. The different identified pathways contributing to charge carrier generation include electron and hole transfer and backtransfer, exciton transfer, and internal relaxation steps. The larger density of states of C70 appears to explain the previously observed larger efficiency for charge separation through hole transfer mechanism. We also analyze the validity of the high-temperature and short-time semiclassical approximations of the FGR theory, where both overestimated and underestimated Marcus theory based constants can be affected.

Optimized attenuated interaction: Enabling stochastic Bethe-Salpeter spectra for large systems

We develop an improved stochastic formalism for the Bethe-Salpeter equation (BSE), based on an exact separation of the effective-interaction W into two parts, W = (W - vW) + vW, where the latter is formally any translationally invariant interaction, vW(r - r'). When optimizing the fit of the exchange kernel vW to W, using a stochastic sampling W, the difference W - vW becomes quite small. Then, in the main BSE routine, this small difference is stochastically sampled. The number of stochastic samples needed for an accurate spectrum is then largely independent of system size. While the method is formally cubic in scaling, the scaling prefactor is small due to the constant number of stochastic orbitals needed for sampling W.

Uncovering the Mechanism of Thermally Activated Delayed Fluorescence in Coplanar Emitters Using Potential Energy Surface Analysis

Planarized emitters exhibiting thermally activated delayed fluorescence (TADF) have attracted attention due to their narrow emission spectra, improved photostability, and high quantum yields, but with large singlet-triplet energy gaps (ΔE_ST) and no heavy atoms, the origin of their TADF remains a subject of debate. Here we prepare two isomeric, coplanar donor-acceptor compounds, with HMAT-2PYM performing dual TADF and room-temperature phosphorescence but with HMAT-4PYM exhibiting only prompt fluorescence. Although conventional TADF design principles suggest that neither isomer should exhibit TADF, we reveal differences in the excited state potential energy surfaces that enable spin-flip processes in only one isomer. We also find that hydrogen bonding is absent between the planar units of these emitters, despite earlier claims of intramolecular hydrogen bonding in similar compounds. Overall, this work demonstrates that potential energy surface analysis is a practical strategy for designing coplanar TADF materials that might otherwise be overlooked by conventional TADF design metrics.

Effects of Solvent Dielectric on Thermally Activated Delayed Fluorescence: A Predictive Computational Polarization Consistent Approach

We study computationally thermally activated delayed fluorescence (TADF) in donor-acceptor compounds. The relevant electronic excited states that are strongly affected by the dielectric environment are treated by a polarization consistent framework. The high fidelity potential energy surfaces are used following a quantum-mechanical Fermi's golden rule (FGR) picture to calculate rates of intersystem crossing (ISC) and reverse intersystem crossing (RISC). To demonstrate the potency of the approach, we consider isomers of benzonitrile functionalized tert-butyl-substituted dimethylacridine (DMAC-BN), which were recently found to perform well as TADF emitters. The calculated excited state energies that appear to reproduce well measured spectral trends with respect to the dielectric constant are used to parametrize ISC/RISC FGR rates. The calculated rates reproduce well measured rates, whereas semiclassical based rates are grossly underestimated. In particular, we find in agreement with the recent experimental study [Phys. Rev. Appl.2019, 12, 044021] that the ortho and meta isomers are significantly more effective as TADF emitters. The computational framework provides valuable insight at the molecular level into RISC rates and therefore can contribute to the design of materials of increased TADF efficiency.

Luminescent 1H-1,3-benzazaphospholes

2-R-1H-1,3-Benzazaphospholes (R-BAPs) are an interesting class of σ²P heterocycles containing P[double bond, length as m-dash]C bonds. While closely related 2-R-1,3-benzoxaphospholes (R-BOPs) have been shown to be highly photoluminescent materials depending on specific R substituents, photoluminescence of R-BAPs has been previously limited to an example having a fused carbazole ring system. Here we detail the synthesis and structural characterization of a new R-BAP (3c, R = 2,2'-dithiophene), and compare its photoluminescence against two previously reported R-BAPs (3a, R, R' = Me and 3b, R = 2-thiophene). The significant fluorescence displayed by the thiophene derivatives 3b (φ = 0.53) and 3c (φ = 0.12) stands in contrast to the weakly emissive methyl substituted analogue 3a (φ = 0.08). Comparative computational investigations of 3a-c offer insights into the interplay between structure-function relationships affecting excited state relaxation processes.

Computational Analysis of the Behavior of BODIPY Decorated Monofunctional Platinum(II) Complexes in the Dark and under Light Irradiation

Dual-action drugs are occupying an important place in the scientific landscape of cancer research owing to the possibility to combine different therapeutic strategies into a single molecule. In the present work, the behavior of two BODIPY-appended monofunctional Pt(II) complexes, one mononuclear and one binuclear, recently synthesized and tested for their cytotoxicity have been explored both in the dark and under light irradiation. Quantum mechanical DFT calculations have been used to carry out the exploration of the key steps, aquation and guanine attack, of the mechanism of action of Pt(II) complexes in the dark. Due to the presence of the BODIPY chromophore and the potential capability of the two investigated complexes to work as photosensitizers in PDT, time dependent DFT has been employed to calculate their photophysical properties and to inspect how the sensitizing properties of BODIPY are affected by the presence of the platinum "heavy atom". Furthermore, also the eventual influence on of the photophysical properties due to the displacement of chlorido ligands by water and of water by guanine has been taken into consideration.

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.

Solvent Dependent Nuclear Magnetic Resonance Molecular Parameters Based on a Polarization Consistent Screened Range Separated Hybrid Density Functional Theory Framework

Nuclear magnetic resonance (NMR) properties of solvated molecules are significantly affected by the solvent. We, therefore, employ a polarization consistent framework that efficiently addresses the solvent polarizing environment effects. Toward this goal a dielectric screened range separated hybrid (SRSH) functional is invoked with a polarizable continuum model (PCM) to properly represent the orbital gap in the condensed phase. We build on the success of range separated hybrid (RSH) functionals to address the erroneous tendency of traditional density functional theory (DFT) to collapse the orbital gap. Recently, the impact of RSH that properly opens up the orbital gap in gas-phase calculations on NMR properties has been assessed. Here, we report the use of SRSH-PCM that produces properly solute orbital gaps in calculating isotropic nuclear magnetic shielding and chemical shift parameters of molecular systems in the condensed phase. We show that in contrast to simpler DFT-PCM approaches, SRSH-PCM successfully follows expected dielectric constant trends.

Accurate Non-Adiabatic Couplings from Optimally-Tuned Range-Separated Hybrid Functionals

Precise theoretical calculations of non-adiabatic couplings, which describe the interaction between two Born-Oppenheimer surfaces, are important for the modeling of radiationless decay mechanisms in photochemical processes. Here, we demonstrate that accurate non-adiabatic couplings can be calculated in the framework of linear-response time-dependent density functional theory by using non-empirical, optimally-tuned range-separated hybrid (OT-RSH) functionals. We focus on molecular radicals, in which ultrafast non-radiative decay plays a crucial role, to find that the OT-RSH functional compares well to wave-function based reference data and competes with the accuracy of semi-empirical CAM-B3LYP calculations. Our findings show that the OT-RSH approach provides very accurate non-adiabatic couplings and, therefore, provides a computationally efficient alternative to wave-function based techniques.

Role of Dielectric Screening in Calculating Excited States of Solvated Azobenzene: A Benchmark Study Comparing Quantum Embedding and Polarizable Continuum Model for Representing the Solvent

The low energy excited states of the conformational isomers of solvated azobenzene are calculated with several DFT methods accounting for the solute-solvent interaction implicitly with the polarizable continuum model or explicitly with subsystem DFT. For the latter, embedding potentials are calculated for 21 sampled snapshots of the solvent molecules. First, we find that accounting for the solvent implicitly or explicitly has little effect on the predicted cis-trans S₁ excitation energy gap. Second, we find that azobenzene's S₁ cis and trans energies are accurate as long as a screened range-separated hybrid exchange-correlation functional is employed. Finally, we also tested a simplified workflow whereby a single, averaged, embedding potential is used. Unfortunately, we find larger deviations against the experiment for the simplified workflow. This highlights a basic flaw in the approach, where the time scale of solvent averaging is much longer than that of the solute's electronic polarization.

Excited-state properties of organic semiconductor dyes as electrically pumped lasing candidates from new optimally tuned range-separated models

Even though many efforts have been devoted to optical lasing in recent years, the realization of lasing by direct electrical excitation of organic semiconductors is hampered mainly due to optical losses from electrical contacts and electrical losses induced by triplets and polarons at high current densities. Hereby, accurately accounting for the electrically pumped organic semiconductor laser diodes (OSLDs) still remains one of the greatest challenges in optoelectronics. In this work, the excited-state characteristics of the organic semiconductor dyes used in the electrically pumped OSLDs have thoroughly been investigated using optimally tuned range-separated hybrids (OT-RSHs). Considering several experimentally known compounds of the electrically pumped OSLDs as working models, several variants of OT-RSHs, their combination forms with the polarizable continuum model (PCM), OT-RSH-PCM, as well as their screened versions accounting for the screening effects by the electron correlation through the scalar dielectric constant, OT-SRSHs, have been proposed for reliable prediction of their emission energies and oscillator strengths in both the gas and solvent phases. The role of involved ingredients in the models, namely, the underlying density functional approximations, short- and long-range exact-like exchange, as well as the range-separation parameter, has been examined in detail. It is shown that the newly designed OT-RSHs with the correct behavior of asymptotic exchange-correlation potential outperform the standard RSHs and other density functionals with both fixed and interelectronic distance-dependent exact-like exchange for describing the excite-state properties of compounds of the electrically pumped OSLDs. Concerning the computational cost of the models, it is unveiled that performing both the optimal tuning procedure and subsequent excited-state computations using OT-RSHs in the gas phase can be considered as a more reliable and affordable framework. Finally, the applicability of the proposed models is also put into a broader perspective for the computational design of several compounds as promising candidates to be used in the OSLD materials. Hopefully, our recommended OT-RSHs can function as efficient models for both the related theoretical modeling and confirming the experimental observations in the field of electrically pumped OSLDs.

Electronic Spectra of C60 Films Using Screened Range Separated Hybrid Functionals

We study computationally the electronic spectra of C₆₀ thin films using the recently developed density functional theory (DFT) framework combining a screened range separated hybrid (SRSH) functional with a polarizable continuum model (PCM). The SRSH-PCM approach achieves excellent correspondence between the frontier orbital's energy levels and the ionization potential and electron affinity of the molecular system at the condensed phase and consequently leads to high quality electronic excitation energies when used in time-dependent DFT calculations. Our calculated excited states reproduce the experimentally main reported spectral peaks at the 3.6-4.6 eV energy range and when addressing excitonic effects also reproduce the red-shifted spectral feature. Notably, we analyze the low-lying peak at 2.7 eV and associate it to an excitonic state.

CTRAMER: An open-source software package for correlating interfacial charge transfer rate constants with donor/acceptor geometries in organic photovoltaic materials

In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)-an open-source software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi's golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system.

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.

Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory

The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.

Regioselectivity of aminomethylation in 3-acetyl-7-hydroxycoumarins: Mannich bases and Betti bases

The atomic electron density distribution in 3-acetyl-7-hydroxycoumarin was analysed to elucidate the regioselectivity of aminomethylation.

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.

How Well Does a Solvated Octa-acid Capsule Shield the Embedded Chromophore? A Computational Analysis Based on an Anisotropic Dielectric Continuum Model

The optical properties of chromophores embedded in a water-solvated dimer of octa-acid that forms a molecular-shaped capsule are investigated. In particular, we address the anisotropic dielectric environment that appears to blue-shift excitation energies compared to the free aqueous chromophores. Recently we reported that using an effective scalar dielectric constant ε ≈ 3 appears to reproduce the measured spectra of the embedded coumarins, suggesting that the capsule provides a significant, albeit not perfect, screening of the aqueous dielectric environment. Here, we report absorption energies using a theoretical treatment that includes continuum solvation affected by an anisotropic dielectric function reflecting the high-dielectric environment outside of the capsule and the low-dielectric region within. We report time-dependent density functional theory calculations using a range-separated functional with the Poisson boundary conditions that model the anisotropic dielectric environment. Our calculations find that the anisotropic environment due to the water-solvated hydrophobic capsule is equivalent to a homogeneous effective dielectric constant of ≈3. The calculated values also appear to reproduce measured absorption of the embedded coumarin, where we study the effect of the hydrophobic capsule on the excited state.