Rubrene: The Interplay between Intramolecular and Intermolecular Interactions Determines the Planarization of Its Tetracene Core in the Solid State

Rapid Estimation of the Intermolecular Electronic Couplings and Charge-Carrier Mobilities of Crystalline Molecular Organic Semiconductors through a Machine Learning Pipeline

Organic semiconductors (OSC) offer tremendous potential across a wide range of (opto)electronic applications. OSC development, however, is often limited by trial-and-error design, with computational modeling approaches deployed to evaluate and screen candidates through a suite of molecular and materials descriptors that generally require hours to days of computational time to accumulate. Such bottlenecks slow the pace and limit the exploration of the vast chemical space comprising OSC. When considering charge-carrier transport in OSC, a key parameter of interest is the intermolecular electronic coupling. Here, we introduce a machine learning (ML) model to predict intermolecular electronic couplings in organic crystalline materials from their three-dimensional (3D) molecular geometries. The ML predictions take only a few seconds of computing time compared to hours by density functional theory (DFT) methods. To demonstrate the utility of the ML predictions, we deploy the ML model in conjunction with mathematical formulations to rapidly screen the charge-carrier mobility anisotropy for more than 60,000 molecular crystal structures and compare the ML predictions to DFT benchmarks.

Electronic Couplings for Triplet-Triplet Annihilation Upconversion in Crystal Rubrene

Triplet-triplet annihilation photon upconversion (TTA-UC) is a process able to repackage two low-frequency photons into light of higher energy. This transformation is typically orchestrated by the electronic degrees of freedom within organic compounds possessing suitable singlet and triplet energies and electronic couplings. In this work, we propose a computational protocol for the assessment of electronic couplings crucial to TTA-UC in molecular materials and apply it to the study of crystal rubrene. Our methodology integrates sophisticated yet computationally affordable approaches to quantify couplings in singlet and triplet energy transfer, the binding of triplet pairs, and the fusion to the singlet exciton. Of particular significance is the role played by charge-transfer states along the b-axis of rubrene crystal, acting as both partial quenchers of singlet energy transfer and mediators of triplet fusion. Our calculations identify the π-stacking direction as holding notable triplet energy transfer couplings, consistent with the experimentally observed anisotropic exciton diffusion. Finally, we have characterized the impact of thermally induced structural distortions, revealing their key role in the viability of triplet fusion and singlet fission. We posit that our approaches are transferable to a broad spectrum of organic molecular materials, offering a feasible means to quantify electronic couplings.

Mechanism of Enhanced Triplet-Triplet Upconversion in Organic Molecules

We use time-dependent density functional theory (TDDFT) to investigate the mechanism of efficient triplet-triplet upconversion (TTU) in certain organic materials. In particular, we focus on materials where some singlets are generated in a two-step spin-nonconserving process (T1 + T1 → T2 → S1). For this mechanism to contribute significantly, the intersystem crossing (ISC) from the high-lying triplet to the singlet (T2 → S1) must outcompete the internal conversion (IC) to the low-lying triplet (T2 → T1). By considering multiple families of materials, we show that the T2 → S1 ISC can be enhanced in a number of ways: the substitution of electron-donating (ED) and electron-withdrawing (EW) groups at appropriate positions; the substitution of bulky groups that distort the molecular geometry; and the substitution of heavy atoms that enhance the spin-orbit coupling (SOC). In the first two cases, the enhancements are consistent with El-Sayed's rule in that rapid T2 → S1 ISC requires significant differences in the characters of the S1 and the T2 wavefunctions. Together, these effects enable a wide tunability of T2 → S1 ISC rates over at least 5 orders of magnitude. Meanwhile, the T2 → T1 IC is inhibited in these systems due to the large T2 - T1 energy gap >0.5 eV, which entails a high energy barrier to the T2 → T1 IC and the prediction of a slow rate regardless of the substituents or the presence of heavy atoms. In this way, tuning the T2 → S1 ISC appears to provide an effective strategy to achieve systematic improvement of TTU materials.

Highly Ordered Small Molecule Organic Semiconductor Thin-Films Enabling Complex, High-Performance Multi-Junction Devices

Organic semiconductors have opened up many new electronic applications, enabled by properties like flexibility, low-cost manufacturing, and biocompatibility, as well as improved ecological sustainability due to low energy use during manufacturing. Most current devices are made of highly disordered thin-films, leading to poor transport properties and, ultimately, reduced device performance as well. Here, we discuss techniques to prepare highly ordered thin-films of organic semiconductors to realize fast and highly efficient devices as well as novel device types. We discuss the various methods that can be implemented to achieve such highly ordered layers compatible with standard semiconductor manufacturing processes and suitable for complex devices. A special focus is put on approaches utilizing thermal treatment of amorphous layers of small molecules to create crystalline thin-films. This technique has first been demonstrated for rubrene─an organic semiconductor with excellent transport properties─and extended to some other molecular structures. We discuss recent experiments that show that these highly ordered layers show excellent lateral and vertical mobilities and can be electrically doped to achieve high n- and p-type conductivities. With these achievements, it is possible to integrate these highly ordered layers into specialized devices, such as high-frequency diodes or completely new device principles for organics, e.g., bipolar 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.

Complete Set of Diketopyrrolopyrrole Centrosymmetrical Cofacial Stacked Pairs

Stacked centrosymmetrical dimers and simultaneously H-bonded and stacked hexamers of thiophene-substituted diketopyrrolopyrrole (ThDPP) were studied using DFT as models for crystals with slipped-stacked molecules in 1D columns. Eight stacked dimer arrangements were found, six of which are driven by the minimisation of electron repulsion and realised by placing the partially negatively charged atoms of the diketopyrrolopyrrole rings below the centre of an adjancent thiophene ring. Four of these stacks are related to N,N'-diacylated derivative. An analogous set of eight stacks was discovered computationally for phenyl-substituted DPP (PhDPP), four of which are known among H-bonded DPP pigments, and one more among N,N'-dialkylated PhDPP derivatives. The results shed more light on the mechanisms that drive the formation of stacks between nonaromatic (DPP) and aromatic (Th, Ph) rings. The excitation energies of the lowest four singlet states computed by TD DFT enabled excitonic coupling and energy separation between Frenkel-resonsnce-type and charge-transfer states to be established, depending on the equilibrium stack geometry.

Functionalized Crystalline N-Trimethyltriindoles: Counterintuitive Influence of Peripheral Substituents on Their Semiconducting Properties

Three crystalline N-trimethyltriindoles endowed with different functionalities at 3, 8 and 13 positions (either unsubstituted or with three methoxy or three acetyl groups attached) are investigated, and clear correlations between the electronic nature of the substituents and their solid-state organization, electronic properties and semiconductor behavior are established. The three compounds give rise to similar columnar hexagonal crystalline structures; however, the insertion of electron-donor methoxy groups results in slightly shorter stacking distances when compared with the unsubstituted derivative, whereas the insertion of electron-withdrawing acetyl groups lowers the crystallinity of the system. Functionalization significantly affects hole mobilities with the triacetyl derivative showing the lowest mobility within the series in agreement with the lower degree of order. However, attaching three methoxy groups also results in lower hole mobility values in the OFETs (0.022 vs. 0.0014 cm² V^-1 s^-1) in spite of the shorter stacking distances. This counterintuitive behavior has been explained with the help of DFT calculations performed to rationalize the interplay between the intramolecular and intermolecular properties, which point to lower transfer integrals in the trimethoxy derivative due to the HOMO wave function extension over the peripheral methoxy groups. The results of this study provide useful insights into how peripheral substituents influence the fundamental charge transport parameters of chemically modified triindole platforms of fundamental importance to design new derivatives with improved semiconducting performance.

Molecular dynamics simulation of organic materials: structure, potentials and the MiCMoS computer platform

Molecular dynamics simulation shows atoms in action: liquids, crystals, mesophases, nanoparticles trace their paths in time, temperature and pressure in a computing platform designed for organic small molecules.

The crystalline state of rubrene materials: intermolecular recognition, isomorphism, polymorphism, and periodic bond-chain analysis of morphologies

Crystal structure analysis and lattice energy calculations of 33 rubrenes provide a background for new ideas on synthesis and planning of transport properties.

Influence of non-covalent interactions in dictating the polarity and mobility of charge carriers in a series of crystalline NDIs: a computational case study

Polycyclic aromatic compounds and their derivatives have emerged as potential molecular entities for air-stable n-type organic semiconductors. In particular, naphthalene diimide (NDI)-derived compounds stand out as one of the most promising classes of molecules that have been studied extensively. There have been a lot of debatable experimental reports on the OFET performance characteristics of some of these materials, which have not yet been resolved completely. Hence, the critical intrinsic aspect of the molecular materials during charge transport in a bulk crystalline state would be essential to categorise the potential candidates. As a case study, in this comprehensive computational approach, we investigated the structural and supramolecular organization in single crystals and the role of those aspects in the bulk carrier transport of a group of selected end-substituted NDI derivatives. A subtle alteration of the end group was observed to result in the modulation of the polarity of charge transport and the charge carrier mobility in the single crystalline state. The disparity is addressed by considering the electronic coupling of the transport states, symmetry of the frontier molecular orbitals and various non-covalent intermolecular interactions. We expect that the present study would benefit towards the rational designing of air-stable n-type organic molecular semiconductors for efficient electronic devices.

Bulk carrier transport properties of a group of selected N-substituted naphthalene diimide derivatives are investigated.

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 Packing: Another Key Point for the Performance of Organic and Polymeric Optoelectronic Materials

ConspectusOptoelectronic material properties are governed by the whole collective of organic moieties, and these aggregate states present the characteristic performance of extended assemblies with different molecular packing, not only of single molecules themselves. Thus, controlling molecular packing is an essential issue for obtaining the optimized optical and electronic properties. It is also a great challenge because of the unclear structures and complicated intermolecular interactions, including dispersion forces, electrostatic interactions and hydrogen bonding. Moreover, upon the introduction of some external force as the stimulus source, dynamic optical properties can be achieved with the transformation of molecular packing in some cases, such as the photoinduced room temperature phosphorescence (RTP) effect, mechanochromic luminescence, the thermal treatment-dependent mechanoluminescence effect, and the optimized nonlinear optical (NLO) property achieved after electric poling. Therefore, it is essential to understand the relation between characteristics of molecular packing and the resultant optoelectronic performance at the molecular level, which becomes increasingly demanding for the further development of functional materials for their applications in organic light-emitting diodes (OLEDs), chemo- and biosensors, organic solar cells, data storage, and anticounterfeiting devices.This Account gives a summary of our research on the molecular design of optoelectronic materials, with the consideration of a molecular uniting effect in different aggregated states, such as crystalline states, thin films, and nanoparticles. Through the systematical investigation of structure-packing-performance relationships, some strategies are afforded to partially control the molecular packing via the tunable size, shape, and configuration of aromatic moieties with different electronic and steric effects, together with different types of substituents as functional units to adjust the intermolecular interactions. The utilization of π-π interactions and hydrogen bonding by rational molecular design is considered as the key point to achieve the bright emission of organic materials, including the RTP and mechanoluminescence effects. Also, the dynamic optoelectronic properties are highlighted with different kinds of stimuli, including light irradiation, mechanical force, thermal treatment, and electric field, which are mainly related to the subtle molecular motions under external force and the changeable molecular packing as the metastable state. These selected examples will not only open a window for further development of organic and polymeric optoelectronic materials by the adjustable molecular packing and noncovalent interactions, but also prompt further advances for more interesting and exciting properties.

A DFT study on the photoelectric properties of rubrene and its derivatives

The way to obtain new materials was usually to introduce some groups to molecules. Correlations among the molecular structure and photoelectric properties of rubrene and its eight derivatives have been studied in this paper, and the influences of the introduction of different electron-donating and withdrawing substituents on the molecular orbital, reorganization energies, absorption spectra, and fluorescence spectra of rubrene and its derivatives were discussed. In the present work, density functional theory calculations were performed at the B3LYP/6-311G(d,p) level to optimize the structure, and TD-DFT was used to calculate the absorption and emission spectra. Quantum chemistry calculation results indicated that the maximum absorption wavelengths are redshifted due to the introduction of cyclopentadienyl or furan groups, the maximum absorption wavelengths are blueshifted while methoxy groups substituted on tetracene backbone. We also discussed the influence of substituents on the molecule structure, which suggests that introducing furan and cyclopentadienyl substituents on the tetracene backbone can increase the rigidity of rubrene and improve the fluorescence intensity. The results of reorganization show that the introduction of cyclopentadienyl or furan groups into rubrene is advantageous to the holes transportation, and the introduction of F groups contributes to the electrons transportation. This study provides an insight into the properties of rubrene and its derivatives and supplies an effective method to design new organic semiconductor materials.

Solid state photodimerization of 9-tert-butyl anthracene ester produces an exceptionally metastable polymorph according to first-principles calculations

Molecular crystal engineering seeks to tune the material properties by controlling the crystal packing.

Read between the Molecules: Computational Insights into Organic Semiconductors

The performance and key electronic properties of molecular organic semiconductors are dictated by the interplay between the chemistry of the molecular core and the intermolecular factors of which manipulation has inspired both experimentalists and theorists. This Perspective presents major computational challenges and modern methodological strategies to advance the field. The discussion ranges from insights and design principles at the quantum chemical level, in-depth atomistic modeling based on multiscale protocols, morphological prediction and characterization as well as energy-property maps involving data-driven analysis. A personal overview of the past achievements and future direction is also provided.

Organic 2D Optoelectronic Crystals: Charge Transport, Emerging Functions, and Their Design Perspective

2D organic semiconductor crystals are emerging as a fascinating platform with regard to their applications in organic field-effect transistors (OFETs), attributed to their enhanced charge transport efficiency and their new optoelectronic functions, based on their unique morphological features. Advances in material processing techniques have not only enabled easy fabrication of few-monolayered 2D nanostructures but also facilitated exploration of the interesting properties induced by characteristic 2D morphologies. However, to date, only a limited number of representative organic semiconductors have been utilized in organic 2D optoelectronics. Therefore, in order to further spur this research, an intuitive crystal engineering principle for realizing organic 2D crystals is required. In this regard, here, not only the important implications of applying 2D structures to OFET devices are discussed but also a crystal engineering protocol is provided that first predicts molecular arrangements depending on the molecular factors, which is followed by realizing 2D supramolecular synthon networks for different molecular packing motifs. It is expected that 2D organic semiconductor crystals developed by this approach will pave a promising way toward next-generation organic 2D optoelectronics.

Electronic structure of twisted and planar rubrene molecules: a density functional study

X-ray absorption spectra (XAS), the density of states (DOS) and the electron density distribution of the HOMO and LUMO for flat and twisted rubrene molecules have been calculated using density functional theory (DFT).

Triplet Separation Drives Singlet Fission after Femtosecond Correlated Triplet Pair Production in Rubrene

Singlet fission, a multistep molecular process in which one photon generates two triplet excitons, holds great technological promise. Here, by applying a combination of transient transmittance and two-dimensional electronic spectroscopy with 5 fs laser pulses, we resolve the full set of fission steps before the onset of spin dephasing. In addition to its role as a viable singlet fission material, single-crystalline rubrene is selected because its energetics and transition dipole alignment uniquely allow for the unambiguous identification of the various fission steps through their contributions to distinct spectroscopic features. The measurements reveal that the neighboring correlated triplet pair achieves its maximum population within 20 fs. Subsequent growth of the triplet signal on picosecond time scales is attributable to spatial separation of the triplets, proceeding nonadiabatically through weakly coupled but near-resonant states. As such, we provide evidence in crystalline rubrene for a singlet fission step that, until now, has not been convincingly observed.

Role of molecular conformations in rubrene polycrystalline films growth from vacuum deposition at various substrate temperatures

We report on the optical and structural characterization of rubrene polycrystalline films fabricated from vacuum deposition with various substrate temperatures (Tsub). Depending on Tsub, the role of twisted and planar rubrene conformational isomers on the properties of rubrene films is focused. The temperature (T)-dependent inverse optical transmission (IOT) and photoluminescence (PL) spectra were performed on these rubrene films. The origins of these IOT and PL peaks are explained in terms of the features from twisted and planar rubrene molecules and of the band characteristics from rubrene molecular solid films. Here, two rarely reported weak-peaks at 2.431 and 2.605 eV were observed from IOT spectra, which are associated with planar rubrene. Besides, the T-dependence of optical bandgap deduced from IOT spectra is discussed with respect to Tsub. Together with IOT and PL spectra, for Tsub > 170 °C, the changes in surface morphology and unit cell volume were observed for the first time, and are attributed to the isomeric transformation from twisted to planar rubrenes during the deposition processes. Furthermore, a unified schematic diagram in terms of Frenkel exciton recombination is suggested to explain the origins of the dominant PL peaks performed on these rubrene films at 15 K.

A theoretical description of charge reorganization energies in molecular organic P-type semiconductors

Charge transport properties of materials composed of small organic molecules are important for numerous optoelectronic applications. A material's ability to transport charges is considerably influenced by the charge reorganization energies of the composing molecules. Hence, predictions about charge-transport properties of organic materials deserve reliable statements about these charge reorganization energies. However, using density functional theory which is mostly used for the predictions, the computed reorganization energies depend strongly on the chosen functional. To gain insight, a benchmark of various density functionals for the accurate calculation of charge reorganization energies is presented. A correlation between the charge reorganization energies and the ionization potentials is found which suggests applying IP-tuning to obtain reliable values for charge reorganization energies. According to benchmark investigations with IP-EOM-CCSD single-point calculations, the tuned functionals provide indeed more reliable charge reorganization energies. Among the standard functionals, ωB97X-D and SOGGA11X yield accurate charge reorganization energies in comparison with IP-EOM-CCSD values. © 2016 Wiley Periodicals, Inc.

Controlling electronic effects and intermolecular packing in contorted polyaromatic hydrocarbons (c-PAHs): towards high mobility field effect transistors

Classical MD simulations followed by DFT calculations for computationally screened contorted polyaromatic hydrocarbons reveal encouraging OFET potential.