Ionization Energies, Electron Affinities, and Polarization Energies of Organic Molecular Crystals: Quantitative Estimations from a Polarizable Continuum Model (PCM)-Tuned Range-Separated Density Functional Approach

Chain Conformation and Exciton Delocalization in a Push-Pull Conjugated Polymer

Linear and nonlinear optical line shapes reveal details of excitonic structure in polymer semiconductors. We implement absorption, photoluminescence, and transient absorption spectroscopies in DPP-DTT, an electron push-pull copolymer, to explore the relationship between their spectral line shapes and chain conformation, deduced from resonance Raman spectroscopy and from ab initio calculations. The viscosity of precursor polymer solutions before film casting displays a transition that suggests gel formation above a critical concentration. Upon crossing this viscosity deflection concentration, the line shape analysis of the absorption spectra within a photophysical aggregate model reveals a gradual increase in interchain excitonic coupling. We also observe a red-shifted and line-narrowed steady-state photoluminescence spectrum along with increasing resonance Raman intensity in the stretching and torsional modes of the dithienothiophene unit, which suggests a longer exciton coherence length along the polymer-chain backbone. Furthermore, we observe a change of line shape in the photoinduced absorption component of the transient absorption spectrum. The derivative-like line shape may originate from two possibilities: a new excited-state absorption or Stark effect, both of which are consistent with the emergence of a high-energy shoulder as seen in both photoluminescence and absorption spectra. Therefore, we conclude that the exciton is more dispersed along the polymer chain backbone with increasing concentrations, leading to the hypothesis that polymer chain order is enhanced when the push-pull polymers are processed at higher concentrations. Thus, tuning the microscopic chain conformation by concentration would be another factor of interest when considering the polymer assembly pathways for pursuing large-area and high-performance organic optoelectronic devices.

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.

Understanding Photophysical Properties of Molecules Relevant in Organic Semiconductor Laser Diodes from Electron Localization Function-Tuned and Solvent-Tuned Range-Separated Functionals

Organic semiconductor laser diodes (OSLDs) are prevalent in optoelectronics because of their sustainable energy applications. Organic molecules used in such diodes are usually large; hence, their studies are computationally challenging with high-end benchmark methods. Computational methods with reliable accuracy and efficiency are always indispensable. In the present work, we have applied our computationally inexpensive, nonempirically tuned [electron localization function (ELF*) and solvent (Sol*)] range-separated (RS) functionals to study five molecules used in OSLDs. The emission energies in three different environments [toluene, CBP (4,4'-bis(n-carbazolyl)-1,1'-biphenyl) film, and gas] have been computed with the tuned functionals and compared with the experimental emission energies. ELF* and Sol* functionals can accurately reproduce emission energies in toluene and CBP film environments. On the other hand, both ELF* and IP-tuned functionals with excited-state geometry (IP*) perform better in the gas phase. In addition, a comparative study is performed between time-dependent density functional theory and the Tamm-Dancoff approximation. Along with the emission energy, oscillator strength values have also been reported. Different IP-tuned RS parameters were obtained with the ground- and excited-state geometries. Interestingly, it has been observed that the optimally tuned RS parameter with excited-state geometry (IP*) performs better compared to that with ground-state geometries (IP). Fractional occupation calculations show that the tuned functionals exhibit less localization and delocalization error. The study envisages that ELF* and Sol* functionals can be used to design future candidates for OSLDs.

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.

Protonated Organic Semiconductors: Origin of Water-Induced Charge-Trap Generation

Despite dramatic improvements in the electronic characteristics of organic semiconductors, the low operational stability of organic field-effect transistors (OFETs) hinders their direct use in practical applications. Although the literature contains numerous reports on the effects of water on the operational stability of OFETs, the underlying mechanisms of trap generation induced by water remain unclear. Here, a protonation-induced trap generation of organic semiconductors is proposed as a possible origin of the operational instability in organic field-effect transistors. Spectroscopic and electronic investigation techniques combined with simulations reveal that the direct protonation of organic semiconductors by water during operation may be responsible for the trap generation induced by bias stress; this phenomenon is independent of the trap generation at an insulator surface. In addition, the same feature occurred in small-bandgap polymers with fused thiophene rings irrespective of their crystalline ordering, implying the generality of protonation induced trap generation in various polymer semiconductors with a small bandgap. The finding of the trap-generation process provides new perspectives for achieving greater operational stability of organic field-effect transistors.

Computational Approaches for Organic Semiconductors: From Chemical and Physical Understanding to Predicting New Materials

While a complete understanding of organic semiconductor (OSC) design principles remains elusive, computational methods─ranging from techniques based in classical and quantum mechanics to more recent data-enabled models─can complement experimental observations and provide deep physicochemical insights into OSC structure-processing-property relationships, offering new capabilities for in silico OSC discovery and design. In this Review, we trace the evolution of these computational methods and their application to OSCs, beginning with early quantum-chemical methods to investigate resonance in benzene and building to recent machine-learning (ML) techniques and their application to ever more sophisticated OSC scientific and engineering challenges. Along the way, we highlight the limitations of the methods and how sophisticated physical and mathematical frameworks have been created to overcome those limitations. We illustrate applications of these methods to a range of specific challenges in OSCs derived from π-conjugated polymers and molecules, including predicting charge-carrier transport, modeling chain conformations and bulk morphology, estimating thermomechanical properties, and describing phonons and thermal transport, to name a few. Through these examples, we demonstrate how advances in computational methods accelerate the deployment of OSCsin wide-ranging technologies, such as organic photovoltaics (OPVs), organic light-emitting diodes (OLEDs), organic thermoelectrics, organic batteries, and organic (bio)sensors. We conclude by providing an outlook for the future development of computational techniques to discover and assess the properties of high-performing OSCs with greater accuracy.

TDDFT versus GW/BSE Methods for Prediction of Light Absorption and Emission in a TADF Emitter

Design concepts for organic light emitting diode (OLED) emitters, which exhibit thermally activated delayed fluorescence (TADF) and thereby achieve quantum yields exceeding 25%, depend on singlet-triplet splitting energies of order kT to allow reverse intersystem crossing at ambient temperatures. Simulation methods for these systems must be able to treat relatively large organic molecules, as well as predict their excited state energies, transition energies, singlet-triplet splittings, and absorption and emission cross sections with reasonable accuracy, in order to prove useful in the design process. Here we compare predictions of TDDFT with M06-2X and ωB97X-D exchange-correlation functionals and a G_oW_o@HF/BSE method for these quantities in the well-studied DPTZ-DBTO2 TADF emitter molecule. Geometry optimization is performed for ground state (GS) and lowest donor-acceptor charge transfer (CT) state for each functional. Optical absorption and emission cross sections and energies are calculated at these geometries. Relaxation energies are on the order of 0.5 eV, and the importance of obtaining excited state equilibrium geometries in predicting delayed fluorescence is demonstrated. There are clear trends in predictions of G_oW_o@HF/BSE, and TDDFT/ωB97X-D and M06-2X methods in which the former method favors local exciton (LE) states while the latter favors DA CT states and ωB97X-D makes intermediate predictions. G_oW_o@HF/BSE suffers from triplet instability for LE states but not CT states relevant for TADF. Shifts in HOMO and LUMO levels on adding a conductor-like polarizable continuum model dielectric background are used to estimate changes in excitation energies on going from the gas phase to a solvated molecule.

First Hyperpolarizabilities of Intramolecular Charge-Transfer Architectures Based on Acenaphthene Derivatives in Gas, Solution, and Solid States

Constructing charge transfer (CT) systems and packing arrangement are common and effective methods to control the efficiency of nonlinear optical (NLO) materials. Apart from the traditional through-bond CT (TBCT) systems, through-space CT (TSCT) also leads to distinctive optical and electronic properties. Meanwhile, corresponding theoretical investigations of the aggregation effect are highly desired. In this work, some TSCT and model compounds incorporating acenaphthene as a scaffold and triphenylamine (TPA) as the donor are theoretically performed to systematically reveal the effect of both solvent and solid environments on their static first hyperpolarizabilities (β_tot) by using the polarizable continuum model (PCM) and the combined quantum mechanics and molecular mechanics (QM/MM) method. Results indicate that the dichloromethane solvent effect within the PCM approach causes an almost 2 times increase of the β_tot values. Besides, the different packing modes and intermolecular interactions have remarkable influence on the second-order NLO properties. For the case of TPA-ace-CN in the crystal state, the parallel arrangement will lead to large NLO responses (4.9 × 10^-30 esu) compared to the correspondingly isolated molecule (3.4 × 10^-30 esu). However, for the TPA-ace-TRZ compound with the TSCT architecture, selection of the molecular arrangement may make the aggregate ineffective due to the offset of the through-space dipole and charge transfer between D-A groups, which lead to the β_tot values decreasing from 15.2 × 10^-30 to 7.7 × 10^-30 esu. We believe that our calculation will serve as a guide for the exploration of more efficient NLO materials wherein the molecules are oriented in their most favorable arrangements.

Virtual Screening for Organic Solar Cells and Light Emitting Diodes

The field of organic semiconductors is multifaceted and the potentially suitable molecular compounds are very diverse. Representative examples include discotic liquid crystals, dye-sensitized solar cells, conjugated polymers, and graphene-based low-dimensional materials. This huge variety not only represents enormous challenges for synthesis but also for theory, which aims at a comprehensive understanding and structuring of the plethora of possible compounds. Eventually computational methods should point to new, better materials, which have not yet been synthesized. In this perspective, it is shown that the answer to this question rests upon the delicate balance between computational efficiency and accuracy of the methods used in the virtual screening. To illustrate the fundamentals of virtual screening, chemical design of non-fullerene acceptors, thermally activated delayed fluorescence emitters, and nanographenes are discussed.

Free charge photogeneration in a single component high photovoltaic efficiency organic semiconductor

Organic photovoltaics (OPVs) promise cheap and flexible solar energy. Whereas light generates free charges in silicon photovoltaics, excitons are normally formed in organic semiconductors due to their low dielectric constants, and require molecular heterojunctions to split into charges. Recent record efficiency OPVs utilise the small molecule, Y6, and its analogues, which – unlike previous organic semiconductors – have low band-gaps and high dielectric constants. We show that, in Y6 films, these factors lead to intrinsic free charge generation without a heterojunction. Intensity-dependent spectroscopy reveals that 60–90% of excitons form free charges at AM1.5 light intensity. Bimolecular recombination, and hole traps constrain single component Y6 photovoltaics to low efficiencies, but recombination is reduced by small quantities of donor. Quantum-chemical calculations reveal strong coupling between exciton and CT states, and an intermolecular polarisation pattern that drives exciton dissociation. Our results challenge how current OPVs operate, and renew the possibility of efficient single-component OPVs.

When light hits organic semiconductors, bound charge pairs, called excitons, are usually produced. Here, the authors show that in the best performing organic solar material to date, free charges, rather than excitons, are directly created by light.

Poisoning density functional theory with benchmark sets of difficult systems

Large benchmark sets like GMTKN55 [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184] let us analyse the performance of density functional theory over a diverse range of systems and bonding types. However, assessing over a large and diverse set can miss cases where approaches fail badly, and can give a misleading sense of security. To this end we introduce a series of 'poison' benchmark sets, P30-5, P30-10 and P30-20, comprising systems with up to 5, 10 and 20 atoms, respectively. These sets represent the most difficult-to-model systems in GMTKN55. We expect them to be useful in developing new approximations, identifying weak points in existing ones, and to aid in selecting appropriate DFAs for computational studies involving difficult physics, e.g. catalysis.

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.

Virtual Screening of TADF Emitters for Single-Layer OLEDs

Thermally-activated delayed fluorescence (TADF) is a concept which helps to harvest triplet excitations, boosting the efficiency of an organic light-emitting diode. TADF can be observed in molecules with spatially separated donor and acceptor groups with a reduced triplet-singlet energy level splitting. TADF materials with balanced electron and hole transport are attractive for realizing efficient single-layer organic light emitting diodes, greatly simplifying their manufacturing and improving their stability. Our goal here is to computationally screen such materials and provide a comprehensive database of compounds with a range of emission wavelengths, ionization energies, and electron affinities.

Two-dimensional hole gas in organic semiconductors

A highly conductive metallic gas that is quantum mechanically confined at a solid-state interface is an ideal platform to explore non-trivial electronic states that are otherwise inaccessible in bulk materials. Although two-dimensional electron gases have been realized in conventional semiconductor interfaces, examples of two-dimensional hole gases, the counterpart to the two-dimensional electron gas, are still limited. Here we report the observation of a two-dimensional hole gas in solution-processed organic semiconductors in conjunction with an electric double layer using ionic liquids. A molecularly flat single crystal of high-mobility organic semiconductors serves as a defect-free interface that facilitates two-dimensional confinement of high-density holes. A remarkably low sheet resistance of 6 kΩ and high hole-gas density of 10¹⁴ cm^-2 result in a metal-insulator transition at ambient pressure. The measured degenerate holes in the organic semiconductors provide an opportunity to tailor low-dimensional electronic states using molecularly engineered heterointerfaces.

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.

Electrostatic gating of ion transport in carbon nanotube porins: A modeling study

Carbon nanotube porins (CNTPs) are biomimetic membrane channels that demonstrate excellent biocompatibility and unique water and ion transport properties. Gating transport in CNTPs with external voltage could increase control over ion flow and selectivity. Herein, we used continuum modeling to probe the parameters that enable and further affect CNTP gating efficiency, including the size and composition of the supporting lipid membrane, slip flow in the carbon nanotube, and the intrinsic electronic properties of the nanotube. Our results show that the optimal gated CNTP device consists of a semiconducting CNTP inserted into a small membrane patch containing an internally conductive layer. Moreover, we demonstrate that the ionic transport modulated by gate voltages is controlled by the charge distribution along the CNTP under the external gate electric potential. The theoretical understanding developed in this study offers valuable guidance for the design of gated CNTP devices for nanofluidic studies, novel biomimetic membranes, and cellular interfaces in the future.

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.

Tunable lifetimes and efficiencies of room temperature phosphorescent liquids by modulating the length and number of alkyl chains

Organic room temperature phosphorescence (RTP) liquid composites exhibit the potential to make innovative changes in large area flexible lighting applications, and it is extremely challenging to achieve high-efficiency RTP in pure organic solvent-free liquid systems. The excited state properties and inner lighting mechanisms of these composites are unclear; therefore, a theoretical perspective to design high efficiency RTP liquids with tunable lifetime is highly desired. Herein, we systematically investigate the photophysical properties of a series of long swallow-tailed bromonaphthalimide (BT unit) molecules by the newly proposed optimally tuned range-separated (RS) functional method, and a state-of-the-art RTP molecule with an absolute quantum yield (ΦRTP) of 57.1% and a lifetime (τ) of 160 ms in solvent-free liquid is obtained. Moreover, theoretical results show that the energy gap between the lowest singlet excited state (S1) and triplet excited state (T1) can be reduced and the non-radiative energy consumption process can be restricted by modulating the length and number of alkyl chains in organic RTP molecules. Thus, a wise molecular design strategy is proposed and five additional efficient RTP molecules with tunable lifetimes (43, 19, 136, 0.11 and 0.005 ms) and efficiencies (11.3%, 6.8%, 5.9%, 0.2% and 0.05%) are theoretically proposed. This study sheds light on the relationship among molecular structure, lifetime and efficiency, and can provide an important prototype to explore high-efficiency RTP by pure organic solvent-free liquid systems.

Quasiparticles and Band Structures in Organized Nanostructures of Donor-Acceptor Copolymers

The performance of organic semiconductor devices is linked to highly ordered nanostructures of self-assembled molecules and polymers. Many-body perturbation theory is employed to study the excited states in bulk copolymers. The results show that acceptors in the polymer scaffold introduce a, hitherto unrecognized, conduction impurity band that leads to electron localization. The donor states are responsible for the formation of conjugated bands, which are only mildly perturbed by the presence of the acceptors. Along the polymer axis, the nonlocal electronic correlations among copolymer strands hinder efficient band transport, which is, however, strongly enhanced across individual chains. Holes are most effectively transported along the π-π stacking, while electrons in the impurity band follow the edge-to-edge directions. The copolymers exhibit regions with inverted transport polarity, in which electrons and holes are efficiently transported in mutually orthogonal directions.

Impact of Nanomorphology on Surface Doping of Organic Semiconductors: The Pentacene-C60F48 Interface

Establishing the rather complex correlation between the structure and the charge transfer in organic-organic heterostructures is of utmost importance for organic electronics and requires spatially resolved structural, chemical, and electronic details. Insight into this issue is provided here by combining atomic force microscopy, Kelvin probe force microscopy, photoemission electron microscopy, and low-energy electron microscopy for investigating a case study. We select the interface formed by pentacene (PEN), benchmark among the donor organic semiconductors, and a p-type dopant from the family of fluorinated fullerenes. As for Buckminsterfullerene (C₆₀), the growth of its fluorinated derivative C₆₀F₄₈ is influenced by the thickness and crystallinity of the PEN buffer layer, but the behavior is markedly different. We provide a microscopic description of the C₆₀F₄₈/PEN interface formation and analyze the consequences in the electronic properties of the final heterostructure. For just one single layer of PEN, a laterally complete but noncompact C₆₀F₄₈/PEN interface is created, importantly affecting the surface work function. Nonetheless, from the very beginning of the second layer formation, the presence of epitaxial and nonepitaxial PEN domains dramatically influences the growth dynamics and extremely well packed two-dimensional C₆₀F₄₈ islands develop. Insightful elemental maps of the C₆₀F₄₈/PEN surface spatially resolve the nonuniform distribution of the dopant molecules, which leads to a heterogeneous work function landscape.

Molecular-Scale Characterization of Photoinduced Charge Separation in Mixed-Dimensional InSe-Organic van der Waals Heterostructures

Layered indium selenide (InSe) is an emerging two-dimensional semiconductor that has shown significant promise for high-performance transistors and photodetectors. The range of optoelectronic applications for InSe can potentially be broadened by forming mixed-dimensional van der Waals heterostructures with zero-dimensional molecular systems that are widely employed in organic electronics and photovoltaics. Here, we report the spatially resolved investigation of photoinduced charge separation between InSe and two molecules (C₇₀ and C₈-BTBT) using scanning tunneling microscopy combined with laser illumination. We experimentally and computationally show that InSe forms type-II and type-I heterojunctions with C₇₀ and C₈-BTBT, respectively, due to an interplay of charge transfer and dielectric screening at the interface. Laser-excited scanning tunneling spectroscopy reveals a ∼0.25 eV decrease in the energy of the lowest unoccupied molecular orbital of C₇₀ with optical illumination. Furthermore, photoluminescence spectroscopy and Kelvin probe force microscopy indicate that electron transfer from InSe to C₇₀ in the type-II heterojunction induces a photovoltage that quantitatively matches the observed downshift in the tunneling spectra. In contrast, no significant changes are observed upon optical illumination in the type-I heterojunction formed between InSe and C₈-BTBT. Density functional theory calculations further show that, despite the weak coupling between the molecular species and InSe, the band alignment of these mixed-dimensional heterostructures strongly differs from the one suggested by the ionization potential and electronic affinities of the isolated components. Self-energy-corrected density functional theory indicates that these effects are the result of the combination of charge redistribution at the interface and heterogeneous dielectric screening of the electron-electron interactions in the heterostructure. In addition to providing specific insight for mixed-dimensional InSe-organic van der Waals heterostructures, this work establishes a general experimental methodology for studying localized charge transfer at the molecular scale that is applicable to other photoactive nanoscale systems.

Molecular Semiconductors for Logic Operations: Dead-End or Bright Future?

The field of organic electronics has been prolific in the last couple of years, leading to the design and synthesis of several molecular semiconductors presenting a mobility in excess of 10 cm² V^-1 s^-1 . However, it is also started to recently falter, as a result of doubtful mobility extractions and reduced industrial interest. This critical review addresses the community of chemists and materials scientists to share with it a critical analysis of the best performing molecular semiconductors and of the inherent charge transport physics that takes place in them. The goal is to inspire chemists and materials scientists and to give them hope that the field of molecular semiconductors for logic operations is not engaged into a dead end. To the contrary, it offers plenty of research opportunities in materials chemistry.

A computational study of anisotropic charge transport in air-stable fluorinated benzobisbenzothiophene (FBBBT) derivatives

A computational study of anisotropical charge transport properties of fluorinated benzobisbenzothiohphene derivatives (FBBBT) is presented. The values of IP_adia of all FBBBTs are found in the range of 6.00-6.20 eV inferring the fact that the investigated compounds have ambient air-stability. In addition, the energy levels of FBBBT s are found to be lower than those of benzobisbenzothiophene (BBBT) compound indicating higher charge carrier stability in the former. Hirshfield surface analyses showed that, in all the studied compounds, the principal identifiable interaction were mostly due to F⋯H and H⋯H intermolecular couplings with no contribution from S⋯S bondings. The calculated maximum μ_hole(μ_elec) value of the compounds FBBBT-a and FBBBT-b was found to be 0.483 (0.794) cm²V^- 1s^- 1 and 0.688 (0.542) cm²V^- 1s^- 1 respectively in the direction of transistor channel (Φ = 93.39 ^∘(273.30^∘) for FBBBT-a and Φ = 92.24 ^∘/272.72 ^∘ for FBBBT-b). For FBBBT-c, the maximum μ_elec(μ_hole) value of 0.933 (0.233) cm²V^- 1s^- 1 appeared for Φ = 0 ^∘/179.90 ^∘. In addition, the compounds FBBBT-a and FBBBT-b possess two additional fluorine atoms attached at the X positions in the backbone, which result in an increment in μ_elec values (1.4 times and 0.78 times higher than μ_hole) in these two compounds at a particular crystal direction.

Rational Tuning of Molecular Interaction and Energy Level Alignment Enables High-Performance Organic Photovoltaics

The performance of organic photovoltaics (OPVs) has rapidly improved over the past years. Recent work in material design has primarily focused on developing near-infrared nonfullerene acceptors with broadening absorption that pair with commercialized donor polymers; in the meanwhile, the influence of the morphology of the blend film and the energy level alignment on the efficiency of charge separation needs to be synthetically considered. Herein, the selection rule of the donor/acceptor blend is demonstrated by rationally considering the molecular interaction and energy level alignment, and highly efficient OPV devices using both-fluorinated or both-nonfluorinated donor/acceptor blends are realized. With the enlarged absorption, ideal morphology, and efficient charge transfer, the devices based on the PBDB-T-F/Y1-4F blend and PBDB-T-F/Y6 exhibit champion power conversion efficiencies as high as 14.8% and 15.9%, respectively.

Concentration dependent energy levels shifts in donor-acceptor mixtures due to intermolecular electrostatic interaction

Recent progress in the improvement of organic solar cells lead to a power conversion efficiency to over 16%. One of the key factors for this improvement is a more favorable energy level alignment between donor and acceptor materials, which demonstrates that the properties of interfaces between donor and acceptor regions are of paramount importance. Recent investigations showed a significant dependence of the energy levels of organic semiconductors upon admixture of different materials, but its origin is presently not well understood. Here, we use multiscale simulation protocols to investigate the molecular origin of the mixing induced energy level shifts and show that electrostatic properties, in particular higher-order multipole moments and polarizability determine the strength of the effect. The findings of this study may guide future material-design efforts in order to improve device performance by systematic modification of molecular properties.

Prediction of Excited-State Properties of Oligoacene Crystals Using Fragment-Based Quantum Mechanical Method

A fundamental understanding of the excited-state properties of molecular crystals is of central importance for their optoelectronics applications. In this study, we developed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for the quantitative characterization of the excited-state properties of locally excited molecular clusters. The accuracy of the EE-GMF method is systematically assessed for oligoacene crystals. Our result demonstrates that the EE-GMF method is capable of providing the lowest vertical singlet (S₁) and triplet excitation energies (T₁), in excellent agreement with the full-system quantum mechanical calculations. Using this method, we also investigated the performance of different density functionals in predicting the excited-state properties of the oligoacene crystals. Among the 13 tested functionals, B3LYP and MN15 give the two lowest overall mean unsigned errors with reference to the experimental S₁ and T₁ excitation energies. The EE-GMF approach can be readily utilized for studying the excited-state properties of large-scale organic solids at diverse ab initio levels.

Toward Design of Novel Materials for Organic Electronics

Materials for organic electronics are presently used in prominent applications, such as displays in mobile devices, while being intensely researched for other purposes, such as organic photovoltaics, large-area devices, and thin-film transistors. Many of the challenges to improve and optimize these applications are material related and there is a nearly infinite chemical space that needs to be explored to identify the most suitable material candidates. Established experimental approaches struggle with the size and complexity of this chemical space. Herein, the development of simulation methods is addressed, with a particular emphasis on predictive multiscale protocols, to complement experimental research in the identification of novel materials and illustrate the potential of these methods with a few prominent recent applications. Finally, the potential of machine learning and methods based on artificial intelligence is discussed to further accelerate the search for new materials.

Transport and Optical Gaps in Amorphous Organic Molecular Materials

The standard procedure to identify the hole- or electron-acceptor character of amorphous organic materials used in OLEDs is to look at the values of a pair of basic parameters, namely, the ionization potential (IP) and the electron affinity (EA). Recently, using published experimental data, the present authors showed that only IP matters, i.e., materials with IP > 5.7 (<5.7) showing electron (hole) acceptor character. Only three materials fail to obey this rule. This work reports ab initio calculations of IP and EA of those materials plus two materials that behave according to that rule, following a route which describes the organic material by means of a single molecule embedded in a polarizable continuum medium (PCM) characterized by a dielectric constant ε . PCM allows to approximately describe the extended character of the system. This "compound" system was treated within density functional theory (DFT) using several combinations of the functional/basis set. In the preset work ε was derived by assuming Koopmans' theorem to hold. Optimal ε values are in the range 4.4⁻5.0, close to what is expected for this material family. It was assumed that the optical gap corresponds to the excited state with a large oscillator strength among those with the lowest energies, calculated with time-dependent DFT. Calculated exciton energies were in the range 0.76⁻1.06 eV, and optical gaps varied from 3.37 up to 4.50 eV. The results are compared with experimental data.

Beyond Koopmans' theorem: electron binding energies in disordered materials

The topical review focuses on calculating ionization energies (IE), or electronic polarons in quasi-particle terminology, in large disordered systems, e.g. for a solute dissolved in a molecular solvent. The simplest estimate of the ionization energy is provided by one-electron energies in the Hartree-Fock theory, but the calculated quantities are not accurate. Density functional theory as many-body theory provides a principal opportunity for calculating one-electron energies including correlation and relaxation effects, i.e. the true energies of electronic polarons. We argue that such a principal possibility materializes within the concept of optimally tuned range-separated hybrid functionals (OT-RSH). We describe various schemes for optimal tuning. Importantly, the OT-RSH scheme is investigated for systems capped with dielectric continuum, providing a consistent picture on the QM/dielectric boundary. Finally, some limitations and open issues of the OT-RSH approach are addressed.

The role of initial and final states in molecular spectroscopies

Interpreting experimental spectra of thin films of organic semiconductors is challenging, and understanding the relationship between experimental data obtained by different spectroscopic techniques requires a careful consideration of the initial and final states for each process.

Deconstructing the behavior of donor–acceptor copolymers in solution & the melt: the case of PTB7

Molecular dynamics simulations of the donor–acceptor copolymer PTB7 at near experimental scale reveal structure–dynamics correlations in the condensed phase.

Dielectric Screening Meets Optimally Tuned Density Functionals

A short overview of recent attempts at merging two independently developed methods is presented. These are the optimal tuning of a range-separated hybrid (OT-RSH) functional, developed to provide an accurate first-principles description of the electronic structure and optical properties of gas-phase molecules, and the polarizable continuum model (PCM), developed to provide an approximate but computationally tractable description of a solvent in terms of an effective dielectric medium. After a brief overview of the OT-RSH approach, its combination with the PCM as a potentially accurate yet low-cost approach to the study of molecular assemblies and solids, particularly in the context of photocatalysis and photovoltaics, is discussed. First, solvated molecules are considered, with an emphasis on the challenge of balancing eigenvalue and total energy trends. Then, it is shown that the same merging of methods can also be used to study the electronic and optical properties of molecular solids, with a similar discussion of the pros and cons. Tuning of the effective scalar dielectric constant as one recent approach that mitigates some of the difficulties in merging the two approaches is considered.

High-Throughput Screening Approach for the Optoelectronic Properties of Conjugated Polymers

We propose a general high-throughput virtual screening approach for the optical and electronic properties of conjugated polymers. This approach makes use of the recently developed xTB family of low-computational-cost density functional tight-binding methods from Grimme and co-workers, calibrated here to (Time-Dependent) Density Functional Theory ((TD)DFT) data computed for a representative diverse set of (co)polymers. Parameters drawn from the resulting calibration using a linear model can then be applied to the xTB derived results for new polymers, thus generating near DFT-quality data with orders of magnitude reduction in computational cost. As a result, after an initial computational investment for calibration, this approach can be used to quickly and accurately screen on the order of thousands of polymers for target applications. We also demonstrate that the (opto)electronic properties of the conjugated polymers show only a very minor variation when considering different conformers and that the results of high-throughput screening are therefore expected to be relatively insensitive with respect to the conformer search methodology applied.

Benchmark study of ionization potentials and electron affinities of armchair single-walled carbon nanotubes using density functional theory

The intrinsic parameters of carbon nanotubes (CNTs) such as ionization potential (IP) and electron affinity (EA) are closely related to their unique properties and associated applications. In this work, we demonstrated the success of optimal tuning method based on range-separated (RS) density functionals for both accurate and efficient prediction of vertical IPs and electron affinities (EAs) of a series of armchair single-walled carbon nanotubes C_20n H₂₀ (n  =  2-6) compared to the high-level IP/EA equation-of-motion coupled-cluster method with single and double substitutions (IP/EA-EOM-CCSD). Notably, the resulting frontier orbital energies (-ε _HOMO and -ε _LUMO) from the tuning method exhibit an excellent approximation to the corresponding IPs and EAs, that significantly outperform other conventional density functionals. In addition, it is suggested that the RS density functionals that possess both a fixed amount of exact exchange in the short-range and a correct long-range asymptotic behavior are suitable for calculating electronic structures of finite-sized CNTs. Next the performance of density functionals for description of various molecular properties such as chemical potential, hardness and electrophilicity are assessed as a function of tube length. Thanks to the efficiency and accuracy of this tuning method, the related behaviors of much longer armchair single-walled CNTs until C₂₀₀H₂₀ were studied. Lastly, the present work is proved to provide an efficient theoretical tool for future materials design and reliable characterization of other interesting properties of CNT-based systems.