Promises and Perils of Big Data: Philosophical Constraints on Chemical Ontologies

Chemistry is experiencing a paradigm shift in the way it interacts with data. So-called "big data" are collected and used at unprecedented scales with the idea that algorithms can be designed to aid in chemical discovery. As data-enabled practices become ever more ubiquitous, chemists must consider the organization and curation of their data, especially as it is presented to both humans and increasingly intelligent algorithms. One of the most promising organizational schemes for big data is a construct termed an ontology. In data science, ontologies are systems that represent relations among objects and properties in a domain of discourse. As chemistry encounters larger and larger data sets, the ontologies that support chemical research will likewise increase in complexity, and the future of chemistry will be shaped by the choices made in developing big data chemical ontologies. How such ontologies will work should therefore be a subject of significant attention in the chemical community. Now is the time for chemists to ask questions about ontology design and use: How should chemical data be organized? What can be reasonably expected from an organizational structure? Is a universal ontology tenable? As some of these questions may be new to chemists, we recommend an interdisciplinary approach that draws on the long history of philosophers of science asking questions about the organization of scientific concepts, constructs, models, and theories. This Perspective presents insights from these long-standing studies and initiates new conversations between chemists and philosophers.

Combination of Counterion Size and Doping Concentration Determines the Electronic and Thermoelectric Properties of Semiconducting Polymers

In both chemical and electrochemical doping of organic semiconductors (OSCs), a counterion, either from the electrolyte or ionized dopant, balances the charge introduced to the OSC. Despite the large influence of this counterion on OSC optical and electronic response, there remains substantial debate on how a fundamental parameter, ion size, impacts these properties. This work resolves much of this debate by accounting for two doping regimes. In the low-doping regime, the Coulomb binding energies between charge carriers on the OSC and the counterions are significant, and larger counterions lead to decreased Coulomb interactions, more delocalized charge carriers, and higher electrical conductivities. In the high-doping regime, the Coulomb binding energies become negligible due to the increased dielectric constant of the films and a smoothing of the energy landscape; thereby, the electrical conductivities depend primarily on the extent of morphological disorder in the OSC. Moreover, in regioregular poly(3-hexylthiophene), rr-P3HT, smaller counterions lead to greater bipolaron concentrations in the low-doping regime due to the increased Coulomb interactions. Emphasizing the impact of the counterion size, it is shown that larger counterions can lead to increased thermoelectric power factors for rr-P3HT.

Mo3S13 Chalcogel: A High-Capacity Electrode for Conversion-Based Li-Ion Batteries

Despite large theoretical energy densities, metal-sulfide electrodes for energy storage systems face several limitations that impact the practical realization. Here, we present the solution-processable, room temperature (RT) synthesis, local structures, and application of a sulfur-rich Mo3S13 chalcogel as a conversion-based electrode for lithium-sulfide batteries (LiSBs). The structure of the amorphous Mo3S13 chalcogel is derived through operando Raman spectroscopy, synchrotron X-ray pair distribution function (PDF), X-ray absorption near edge structure (XANES), and extended X-ray absorption fine structure (EXAFS) analysis, along with ab initio molecular dynamics (AIMD) simulations. A key feature of the three-dimensional (3D) network is the connection of Mo3S13 units through S-S bonds. Li/Mo3S13 half-cells deliver initial capacity of 1013 mAh g-1 during the first discharge. After the activation cycles, the capacity stabilizes and maintains 312 mAh g-1 at a C/3 rate after 140 cycles, demonstrating sustained performance over subsequent cycling. Such high-capacity and stability are attributed to the high density of (poly)sulfide bonds and the stable Mo-S coordination in Mo3S13 chalcogel. These findings showcase the potential of Mo3S13 chalcogels as metal-sulfide electrode materials for LiSBs.

Beyond n-dopants for organic semiconductors: use of bibenzo[d]imidazoles in UV-promoted dehalogenation reactions of organic halides

2,2'-Bis(4-dimethylaminophenyl)- and 2,2'-dicyclohexyl-1,1',3,3'-tetramethyl-2,2',3,3'-tetrahydro-2,2'-bibenzo[d]imidazole ((N-DMBI)2 and (Cyc-DMBI)2) are quite strong reductants with effective potentials of ca. -2 V vs ferrocenium/ferrocene, yet are relatively stable to air due to the coupling of redox and bond-breaking processes. Here, we examine their use in accomplishing electron transfer-induced bond-cleavage reactions, specifically dehalogenations. The dimers reduce halides that have reduction potentials less cathodic than ca. -2 V vs ferrocenium/ferrocene, especially under UV photoexcitation (using a 365 nm LED). In the case of benzyl halides, the products are bibenzyl derivatives, whereas aryl halides are reduced to the corresponding arenes. The potentials of the halides that can be reduced in this way, quantum-chemical calculations, and steady-state and transient absorption spectroscopy suggest that UV irradiation accelerates the reactions via cleavage of the dimers to the corresponding radical monomers.

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.

Multiexciton quintet state populations in a rigid pyrene-bridged parallel tetracene dimer

The multiexciton quintet state, 5TT, generated as a singlet fission intermediate in pairs of molecular chromophores, is a promising candidate as a qubit or qudit in future quantum information science schemes. In this work, we synthesize a pyrene-bridged parallel tetracene dimer, TPT, with an optimized interchromophore coupling strength to prevent the dissociation of 5TT to two decorrelated triplet (T1) states, which would contaminate the spin-state mixture. Long-lived and strongly spin-polarized pure 5TT state population is observed via transient absorption spectroscopy and transient/pulsed electron paramagnetic resonance spectroscopy, and its lifetime is estimated to be >35 µs, with the dephasing time (T2) for the 5TT-based qubit measured to be 726 ns at 10 K. Direct relaxation from 1TT to the ground state does diminish the overall excited state population, but the exclusive 5TT population at large enough persistent density for pulsed echo determination of spin coherence time is consistent with recent theoretical models that predict such behavior for strict parallel chromophore alignment and large exchange coupling.

Thin-Film Organic Heteroepitaxy

Incorporating crystalline organic semiconductors into electronic devices requires understanding of heteroepitaxy given the ubiquity of heterojunctions in these devices. However, while rules for commensurate epitaxy of covalent or ionic inorganic material systems are known to be dictated by lattice matching constraints, rules for heteroepitaxy of molecular systems are still being written. Here, it is found that lattice matching alone is insufficient to achieve heteroepitaxy in molecular systems, owing to weak intermolecular forces that describe molecular crystals. It is found that, in addition, the lattice matched plane also must be the lowest energy surface of the adcrystal to achieve one-to-one commensurate molecular heteroepitaxy over a large area. Ultraviolet photoelectron spectroscopy demonstrates the lattice matched interface to be of higher electronic quality than a disordered interface of the same materials.

New Chemical Dopant and Counterion Mechanism for Organic Electrochemical Transistors and Organic Mixed Ionic-Electronic Conductors

Organic mixed ionic-electronic conductors (OMIECs) have varied performance requirements across a diverse application space. Chemically doping the OMIEC can be a simple, low-cost approach for adapting performance metrics. However, complex challenges, such as identifying new dopant materials and elucidating design rules, inhibit its realization. Here, these challenges are approached by introducing a new n-dopant, tetrabutylammonium hydroxide (TBA-OH), and identifying a new design consideration underpinning its success. TBA-OH behaves as both a chemical n-dopant and morphology additive in donor acceptor co-polymer naphthodithiophene diimide-based polymer, which serves as an electron transporting material in organic electrochemical transistors (OECTs). The combined effects enhance OECT transconductance, charge carrier mobility, and volumetric capacitance, representative of the key metrics underpinning all OMIEC applications. Additionally, when the TBA+ counterion adopts an "edge-on" location relative to the polymer backbone, Coulombic interaction between the counterion and polaron is reduced, and polaron delocalization increases. This is the first time such mechanisms are identified in doped-OECTs and doped-OMIECs. The work herein therefore takes the first steps toward developing the design guidelines needed to realize chemical doping as a generic strategy for tailoring performance metrics in OECTs and OMIECs.

Hydrogen-Bonding Trends in a Bithiophene with 3- and/or 4-Pyridyl Substituents

To improve the charge-carrier transport capabilities of thin-film organic materials, the intermolecular electronic couplings in the material should be maximized. Decreasing intermolecular distance while maintaining proper orbital overlap in highly conjugated aromatic molecules has so far been a successful way to increase electronic coupling. We attempted to decrease the intermolecular distance in this study by synthesizing cocrystals of simple benzoic acid coformers and dipyridyl-2,2'-bithiophene molecules to understand how the coformer identity and pyridine N atom placement affected solid-state properties. We found that with the 5-(3-pyridyl)-5'-(4-pyridyl)-isomer, the 4-pyridyl ring interacted with electrophiles and protons more strongly. Synthesized cocrystal powders were found to have reduced average crystallite size in reference to the parent compounds. The opposite was found for the intermolecular electronic couplings, as determined via density functional theory (DFT) calculations, which were relatively large in some of the cocrystals.

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.

Relating reorganization energies, exciton diffusion length and non-radiative recombination to the room temperature UV-vis absorption spectra of NF-SMA

Understanding excited-state reorganization energies, exciton diffusion lengths and non-radiative (NR) recombination, and the overall optoelectronic responses of nonfullerene small molecule acceptors (NF-SMAs) is important in order to rationally design new materials with controlled properties. While the effects of structural modifications on the optical gaps and electron affinities of NF-SMAs have been studied extensively, analyses of their absorption spectra that carefully characterize electronic and vibrational contributions that allow comparisons of reorganization energies and their implications for exciton diffusion lengths and NR recombination have yet to be reported. Here, we study the room temperature absorption spectra of three structural classes of NF-SMAs in dilute solutions through multiparameter Franck Condon (MFC) analyses and density functional theory (DFT) calculations. We show that the absorption spectra of these NF-SMAs can be categorized based on molecular structure-spectra correlation. The absorption spectra of curved, Y6-like structures can be described using an MFC model with two electronic transitions and two effective vibrational modes. The results of MFC/DFT analyses reveal that Y6 exhibits the smallest intra-molecular reorganization energy among the materials studied. Linear ITIC-like molecular structures reveal larger reorganization energies and reduced conformational uniformity compared to Y6. Meanwhile structures such as IDTBR and IEICO, which have an extra π-conjugated moiety between the donor and acceptor moieties, have large excited-state reorganization energies and low degrees of conformational uniformity. Since the intra-molecular reorganization energy is correlated with exciton diffusion length and nonradiative voltage losses (ΔV_nr), our results highlight the power of RT absorption spectroscopy and DFT calculations as simple tools to designing improved OSCs materials with small reorganization energies, small ΔV_nr, large exciton diffusion length and low energetic disorder (due to a strongly dominant conformation).

Electronic, redox, and optical property prediction of organic π-conjugated molecules through a hierarchy of machine learning approaches

Accelerating the development of π-conjugated molecules for applications such as energy generation and storage, catalysis, sensing, pharmaceuticals, and (semi)conducting technologies requires rapid and accurate evaluation of the electronic, redox, or optical properties. While high-throughput computational screening has proven to be a tremendous aid in this regard, machine learning (ML) and other data-driven methods can further enable orders of magnitude reduction in time while at the same time providing dramatic increases in the chemical space that is explored. However, the lack of benchmark datasets containing the electronic, redox, and optical properties that characterize the diverse, known chemical space of organic π-conjugated molecules limits ML model development. Here, we present a curated dataset containing 25k molecules with density functional theory (DFT) and time-dependent DFT (TDDFT) evaluated properties that include frontier molecular orbitals, ionization energies, relaxation energies, and low-lying optical excitation energies. Using the dataset, we train a hierarchy of ML models, ranging from classical models such as ridge regression to sophisticated graph neural networks, with molecular SMILES representation as input. We observe that graph neural networks augmented with contextual information allow for significantly better predictions across a wide array of properties. Our best-performing models also provide an uncertainty quantification for the predictions. To democratize access to the data and trained models, an interactive web platform has been developed and deployed.

Data storage architectures to accelerate chemical discovery: data accessibility for individual laboratories and the community

As buzzwords like "big data," "machine learning," and "high-throughput" expand through chemistry, chemists need to consider more than ever their data storage, data management, and data accessibility, whether in their own laboratories or with the broader community. While it is commonplace for chemists to use spreadsheets for data storage and analysis, a move towards database architectures ensures that the data can be more readily findable, accessible, interoperable, and reusable (FAIR). However, making this move has several challenges for those with limited-to-no knowledge of computer programming and databases. This Perspective presents basics of data management using databases with a focus on chemical data. We overview database fundamentals by exploring benefits of database use, introducing terminology, and establishing database design principles. We then detail the extract, transform, and load process for database construction, which includes an overview of data parsing and database architectures, spanning Standard Query Language (SQL) and No-SQL structures. We close by cataloging overarching challenges in database design. This Perspective is accompanied by an interactive demonstration available at https://github.com/D3TaLES/databases_demo. We do all of this within the context of chemical data with the aim of equipping chemists with the knowledge and skills to store, manage, and share their data while abiding by FAIR principles.

The Role of Crystal Packing on the Optical Response of Trialkyltetrelethynyl Acenes

The electronic and optical responses of an organic semiconductor (OSC) are dictated by the chemistries of the molecular or polymer building blocks and how these chromophores pack in the solid state. Understanding the physicochemical natures of these responses are not only critical for determining OSC performance for a particular application, but the UV/visible optical response may also be of potential use to determine aspects of the molecular-scale solid-state packing for crystal polymorphs or thin-film morphologies that are difficult to determine otherwise. To probe these relationships, we report the quantum-chemical investigation of a series of trialkyltetrelethynyl acenes (tetrel = silicon or germanium) that adopt the brickwork (BW), slip-stack (SS), or herringbone (HB) packing configurations; the π-conjugated backbones considered here are pentacene (PEN) and anthradithiophene (ADT). For comparison, HB-packed (unsubstituted) pentacene is also included. Density functional theory (DFT) and G0W0 (single-shot GW) electronic band structures, G0W0-BSE (Bethe-Salpeter Equation)-derived optical spectra, polarized ϵ2 spectra, and distributions of both singlet and triplet exciton wave functions are reported. Configurational disorder is also considered. Further, we evaluate the probability of singlet fission in these materials through energy conservation relationships.

The Solution is the Solution: Data-Driven Elucidation of Solution-to-Device Feature Transfer for π-Conjugated Polymer Semiconductors

The advent of data analytics techniques and materials informatics provides opportunities to accelerate the discovery and development of organic semiconductors for electronic devices. However, the development of engineering solutions is limited by the ability to control thin-film morphology in an immense parameter space. The combination of high-throughput experimentation (HTE) laboratory techniques and data analytics offers tremendous avenues to traverse the expansive domains of tunable variables offered by organic semiconductor thin films. This Perspective outlines the steps required to incorporate a comprehensive informatics methodology into the experimental development of polymer-based organic semiconductor technologies. The translation of solution processing and property metrics to thin-film behavior is crucial to inform efficient HTE for data collection and application of data-centric tools to construct new process-structure-property relationships. We argue that detailed investigation of the solution state prior to deposition in conjunction with thin-film characterization will yield a deeper understanding of the physicochemical mechanisms influencing performance in π-conjugated polymer electronics, with data-driven approaches offering predictive capabilities previously unattainable via traditional experimental means.

Biotinylation as a tool to enhance the uptake of small molecules in Gram-negative bacteria

Antibiotic resistance is a major public health concern. The shrinking selection of effective antibiotics and lack of new development is making the situation worse. Gram-negative bacteria more specifically pose serious threat because of their double layered cell envelope and effective efflux systems, which is a challenge for drugs to penetrate. One promising approach to breach this barrier is the “Trojan horse strategy”. In this technique, an antibiotic molecule is conjugated with a nutrient molecule that helps the antibiotic to enter the cell through dedicated transporters for the nutrient. Here, we explored the approach using biotin conjugation with a florescent molecule Atto565 to determine if biotinylation enhances accumulation. Biotin is an essential vitamin for bacteria and is obtained through either synthesis or uptake from the environment. We found that biotinylation enhanced accumulation of Atto565 in E. coli. However, the enhancement did not seem to be due to uptake through biotin transporters since the presence of free biotin had no observable impact on accumulation. Accumulated compound was mostly in the periplasm, as determined by cell fractionation studies. This was further confirmed through the observation that expression of streptavidin in the periplasm specifically enhanced the accumulation of biotinylated Atto565. This enhancement was not observed when streptavidin was expressed in the cytoplasm indicating no significant distribution of the compound inside the cytoplasm. Using gene knockout strains, plasmid complementation and mutagenesis studies we demonstrated that biotinylation made the compound a better passenger through OmpC, an outer membrane porin. Density functional theory (DFT)-based evaluation of the three-dimensional geometries showed that biotinylation did not directly stabilize the conformation of the compound to make it favorable for the entry through a pore. Further studies including molecular dynamics simulations are necessary to determine the possible mechanisms of enhanced accumulation of the biotinylated Atto565.

Lowering Electrocatalytic CO2 Reduction Overpotential Using N-Annulated Perylene Diimide Rhenium Bipyridine Dyads with Variable Tether Length

We report the design, synthesis, and characterization of four N-annulated perylene diimide (NPDI) functionalized rhenium bipyridine [Re(bpy)] supramolecular dyads. The Re(bpy) scaffold was connected to the NPDI chromophore either directly [Re(py-C0-NPDI)] or via an ethyl [Re(bpy-C2-NPDI)], butyl [Re(bpy-C4-NPDI)], or hexyl [Re(bpy-C6-NPDI)] alkyl-chain spacer. Upon electrochemical reduction in the presence of CO₂ and a proton source, Re(bpy-C2/4/6-NPDI) all exhibited significant current enhancement effects, while Re(py-C0-NPDI) did not. During controlled potential electrolysis (CPE) experiments at E_appl = -1.8 V vs Fc^+/0, Re(bpy-C2/4/6-NPDI) all achieved comparable activity (TON_co ∼ 25) and Faradaic efficiency (FE_co ∼ 94%). Under identical CPE conditions, the standard catalyst Re(dmbpy) was inactive for electrocatalytic CO₂ reduction; only at E_appl = -2.1 V vs Fc^+/0 could Re(dmbpy) achieve the same catalytic performance, representing a 300 mV lowering in overpotential for Re(bpy-C2/4/6-NPDI). At higher overpotentials, Re(bpy-C4/6-NPDI) both outperformed Re(bpy-C2-NPDI), indicating the possibility of coinciding electrocatalytic CO₂ reduction mechanisms that are dictated by tether-length and overpotential. Using UV-vis-nearIR spectroelectrochemistry (SEC), FTIR SEC, and chemical reduction experiments, it was shown that the NPDI-moiety served as an electron-reservoir for Re(bpy), thereby allowing catalytic activity at lower overpotentials. Density functional theory studies probing the optimized geometries and frontier molecular orbitals of various catalytic intermediates revealed that the geometric configuration of NPDI relative to the Re(bpy)-moiety plays a critical role in accessing electrons from the electron-reservoir. The improved performance of Re(bpy-C2/4/6-NPDI)dyads at lower overpotentials, relative to Re(dmbpy), highlights the utility of chromophore electron-reservoirs as a method for lowering the overpotential for CO₂ conversion.

Synthesis, structures, and reactivity of isomers of [RuCp*(1,4-(Me2N)2C6H4)]2

[RuCp*(1,3,5-R₃C₆H₃)]₂ {Cp* = η⁵-pentamethylcyclopentadienyl, R = Me, Et} have previously been found to be moderately air stable, yet highly reducing, with estimated D⁺/0.5D₂ (where D₂ and D⁺ represent the dimer and the corresponding monomeric cation, respectively) redox potentials of ca. -2.0 V vs. FeCp₂^+/0. These properties have led to their use as n-dopants for organic semiconductors. Use of arenes substituted with π-electron donors is anticipated to lead to even more strongly reducing dimers. [RuCp*(1-(Me₂N)-3,5-Me₂C₆H₃)]⁺PF₆^- and [RuCp*(1,4-(Me₂N)₂C₆H₄)]⁺PF₆^- have been synthesized and electrochemically and crystallographically characterized; both exhibit D⁺/D potentials slightly more cathodic than [RuCp*(1,3,5-R₃C₆H₃)]⁺. Reduction of [RuCp*(1,4-(Me₂N)₂C₆H₄)]⁺PF₆^- using silica-supported sodium-potassium alloy leads to a mixture of isomers of [RuCp*(1,4-(Me₂N)₂C₆H₄)]₂, two of which have been crystallographically characterized. One of these isomers has a similar molecular structure to [RuCp*(1,3,5-Et₃C₆H₃)]₂; the central C-C bond is exo,exo, i.e., on the opposite face of both six-membered rings from the metals. A D⁺/0.5D₂ potential of -2.4 V is estimated for this exo,exo dimer, more reducing than that of [RuCp*(1,3,5-R₃C₆H₃)]₂ (-2.0 V). This isomer reacts much more rapidly with both air and electron acceptors than [RuCp*(1,3,5-R₃C₆H₃)]₂ due to a much more cathodic D₂˙⁺/D₂ potential. The other isomer to be crystallographically characterized, along with a third isomer, are both dimerized in an exo,endo fashion, representing the first examples of such dimers. Density functional theory calculations and reactivity studies indicate that the central bonds of these two isomers are weaker than those of the exo,exo isomer, or of [RuCp*(1,3,5-R₃C₆H₃)]₂, leading to estimated D⁺/0.5D₂ potentials of -2.5 and -2.6 V vs. FeCp₂^+/0. At the same time the D₂˙⁺/D₂ potentials for the exo,endo dimers are anodically shifted relative to those of [RuCp*(1,3,5-R₃C₆H₃)]₂, resulting in much greater air stability than for the exo,exo isomer.

OCELOT: An infrastructure for data-driven research to discover and design crystalline organic semiconductors

Materials design and discovery are often hampered by the slow pace and materials and human costs associated with Edisonian trial-and-error screening approaches. Recent advances in computational power, theoretical methods, and data science techniques, however, are being manifest in a convergence of these tools to enable in silico materials discovery. Here, we present the development and deployment of computational materials data and data analytic approaches for crystalline organic semiconductors. The OCELOT (Organic Crystals in Electronic and Light-Oriented Technologies) infrastructure, consisting of a Python-based OCELOT application programming interface and OCELOT database, is designed to enable rapid materials exploration. The database contains a descriptor-based schema for high-throughput calculations that have been implemented on more than 56 000 experimental crystal structures derived from 47 000 distinct molecular structures. OCELOT is open-access and accessible via a web-user interface at https://oscar.as.uky.edu.

Modification of the LiFePO4 (010) Surface Due to Exposure to Atmospheric Gases

First-principles thermodynamics enables the description of the surface chemistry of inorganic materials as a function of temperature and partial pressures of atmospheric gases, providing a framework to connect atomistic simulations with macroscopic materials properties. Here we re-examine the surface chemistry of LiFePO₄ (LFP), a widely studied material for use as the cathode in Li-ion batteries. Our results reveal that at room temperature and under standard pressures the LFP (010) surface is covered with water. At elevated temperatures and reduced H₂ partial pressure, one water molecule loses a hydrogen atom and the preferred binding moieties are OH and H₂O; while further reducing the H₂O partial pressure results in the desorption of water leaving only the OH behind. This work also shines new light on the configuration, and resulting electronic properties, of the LFP (010) surface when molecular oxygen (O₂) is adsorbed. The molecular adsorbates are also shown to have an impact on the LFP surface potentials and magnetic properties. These simulations provide an enhanced picture of the LFP surface chemistry and the potential impact of these adsorbates on understanding the characteristics of LFP in different materials applications.

Suppressing bias stress degradation in high performance solution processed organic transistors operating in air

Solution processed organic field effect transistors can become ubiquitous in flexible optoelectronics. While progress in material and device design has been astonishing, low environmental and operational stabilities remain longstanding problems obstructing their immediate deployment in real world applications. Here, we introduce a strategy to identify the most probable and severe degradation pathways in organic transistors and then implement a method to eliminate the main sources of instabilities. Real time monitoring of the energetic distribution and transformation of electronic trap states during device operation, in conjunction with simulations, revealed the nature of traps responsible for performance degradation. With this information, we designed the most efficient encapsulation strategy for each device type, which resulted in fabrication of high performance, environmentally and operationally stable small molecule and polymeric transistors with consistent mobility and unparalleled threshold voltage shifts as low as 0.1 V under the application of high bias stress in air.

Electrical instability of organic field-effect transistors (OFETs) during operation remains a challenge that limits the device’s real-world technological viability. Here, the authors report a method for diagnosing and suppressing bias stress in solution-processed OFETs operated in air.

A molecular interaction-diffusion framework for predicting organic solar cell stability

Rapid increase in the power conversion efficiency of organic solar cells (OSCs) has been achieved with the development of non-fullerene small-molecule acceptors (NF-SMAs). Although the morphological stability of these NF-SMA devices critically affects their intrinsic lifetime, their fundamental intermolecular interactions and how they govern property-function relations and morphological stability of OSCs remain elusive. Here, we discover that the diffusion of an NF-SMA into the donor polymer exhibits Arrhenius behaviour and that the activation energy E_a scales linearly with the enthalpic interaction parameters χ_H between the polymer and the NF-SMA. Consequently, the thermodynamically most unstable, hypo-miscible systems (high χ) are the most kinetically stabilized. We relate the differences in E_a to measured and selectively simulated molecular self-interaction properties of the constituent materials and develop quantitative property-function relations that link thermal and mechanical characteristics of the NF-SMA and polymer to predict relative diffusion properties and thus morphological stability.

n-type charge transport in heavily p-doped polymers

It is commonly assumed that charge-carrier transport in doped π-conjugated polymers is dominated by one type of charge carrier, either holes or electrons, as determined by the chemistry of the dopant. Here, through Seebeck coefficient and Hall effect measurements, we show that mobile electrons contribute substantially to charge-carrier transport in π-conjugated polymers that are heavily p-doped with strong electron acceptors. Specifically, the Seebeck coefficient of several p-doped polymers changes sign from positive to negative as the concentration of the oxidizing agents FeCl₃ or NOBF₄ increase, and Hall effect measurements for the same p-doped polymers reveal that electrons become the dominant delocalized charge carriers. Ultraviolet and inverse photoelectron spectroscopy measurements show that doping with oxidizing agents results in elimination of the transport gap at high doping concentrations. This approach of heavy p-type doping is demonstrated to provide a promising route to high-performance n-type organic thermoelectric materials.

Steric Manipulation as a Mechanism for Tuning the Reduction and Oxidation Potentials of Phenothiazines

Synthetic chemists customarily tune the redox characteristics of π-conjugated molecules by introducing electron-donating or electron-withdrawing substituents onto the molecular core, or by modifying the length of the π-conjugated pathway. Any steric effects of such efforts on molecular geometry typically affect both the neutral and charged (oxidized or reduced) states indiscriminately. However, in electroactive systems that undergo significant conformational changes upon oxidation or reduction, we can leverage the steric and inductive effects of substitution to attain considerable control over individual redox potentials. Here, we make use of density functional theory to elucidate the interplay between electronic and geometric effects of peripheral substitution on the model system of phenothiazine. For instance, we introduce substituents at positions ortho to the nitrogen atom (positions 1 and 9) to induce steric strain in the radical-cation state without significant effect on the neutral molecule, thereby augmenting the overall ionization potential. Notably, this steric effect persists for electron-donating substituents; the resulting ionization potentials therefore deviate from outcomes foretold by Hammett constants. Moreover, the same procedure has limited effect on electron affinities because of differences in phenothiazines' relaxation process upon reduction compared to oxidation. Our results promote molecular design guidelines for manipulating redox potentials in classes of electroactive compounds that experience dramatic changes in geometry upon ionization.

A Genetic Algorithmic Approach to Determine the Structure of Li-Al Layered Double Hydroxides

Layered double hydroxides (LDH) demonstrate significant potential across a range of applications, including as catalysts, delivery vehicles for pharmaceuticals, environmental remediation, and supercapacitors. Explaining the mechanism of LDH action at the atomic scale in these and other applications is challenging, however, due to the difficulty in precisely defining the bulk and surface structure and chemical compositions. Here, we focus on the determination of the structure of lithium-aluminum (Li-Al) LDH, which has shown promise in the catalytic depolymerization of lignin, both directly as the catalyst and as a support for gold nanoparticles. While the relative positions of the Li and Al metals are generally well resolved by X-ray crystallography, it is the structures of the anionic layers, consisting of water and carbonate, that are less well established. Combinatorial analyses of all possible positions and rotations of the water and carbonate in the three-layered Li-AL LDH polytope reveals that the phase space is much too large to examine in any reasonable time frame in a one-by-one structure exploration. To overcome this limitation, we develop and deploy a genetic algorithm (GA) wherein fitness is determined by matching a calculated X-ray diffraction (XRD) pattern for a given structure to the known experimental XRD pattern. The GA approach results in structures of high fitness that portend the bulk Li-Al LDH structure. Importantly, the GA approach offers the potential to determine the structures of other LDH, and more generally layered materials, which are generally difficult to describe given the large chemical and structural space to be explored.

Determination of the Free Energies of Mixing of Organic Solutions through a Combined Molecular Dynamics and Bayesian Statistics Approach

As new generations of thin-film semiconductors are moving toward solution-based processing, the development of printing formulations will require information pertaining to the free energies of mixing of complex mixtures. From the standpoint of in silico material design, this move necessitates the development of methods that can accurately and quickly evaluate these formulations in order to maximize processing speed and reproducibility. Here, we make use of molecular dynamics (MD) simulations, in combination with the two-phase thermodynamic (2PT) model, to explore the free energy of mixing surfaces for a series of halogenated solvents and high-boiling point solvent additives used in the development of thin-film organic semiconductors. Although the combined methods generally show good agreement with available experimental data, the computational cost to traverse the free-energy landscape is considerable. Hence, we demonstrate how a Bayesian optimization scheme, coupled with the MD and 2PT approaches, can drastically reduce the number of simulations required, in turn shrinking both the computational cost and time.

Nitration of benzothioxanthene: towards a new class of dyes with versatile photophysical properties

The selective mono-nitration of benzothioxanthene (BTXI) is demonstrated here, opening doors to a wide range of structures and reactions. In addition, a new series of derivaties was characterized expanding the scope for electronic organic purposes.

Computationally aided design of a high-performance organic semiconductor: the development of a universal crystal engineering core

Herein, we describe the design and synthesis of a suite of molecules based on a benzodithiophene "universal crystal engineering core". After computationally screening derivatives, a trialkylsilylethyne-based crystal engineering strategy was employed to tailor the crystal packing for use as the active material in an organic field-effect transistor. Electronic structure calculations were undertaken to reveal derivatives that exhibit exceptional potential for high-efficiency hole transport. The promising theoretical properties are reflected in the preliminary device results, with the computationally optimized material showing simple solution processing, enhanced stability, and a maximum hole mobility of 1.6 cm2 V-1 s-1.

Exploiting Excited-State Aromaticity To Design Highly Stable Singlet Fission Materials

Singlet fission, the process of forming two triplet excitons from one singlet exciton, is a characteristic reserved for only a handful of organic molecules due to the atypical energetic requirement for low energy excited triplet states. The predominant strategy for achieving such a trait is by increasing ground state diradical character; however, this greatly reduces ambient stability. Herein, we exploit Baird's rule of excited state aromaticity to manipulate the singlet-triplet energy gap and create novel singlet fission candidates. We achieve this through the inclusion of a [4n] 5-membered heterocycle, whose electronic resonance promotes aromaticity in the triplet state, stabilizing its energy relative to the singlet excited state. Using this theory, we design a family of derivatives of indolonaphthyridine thiophene (INDT) with highly tunable excited state energies. Not only do we access novel singlet fission materials, they also exhibit excellent ambient stability, imparted due to the delocalized nature of the triplet excited state. Spin-coated films retained up to 85% activity after several weeks of exposure to oxygen and light, while analogous films of TIPS-pentacene showed full degradation after 4 days, showcasing the excellent stability of this class of singlet fission scaffold. Extension of our theoretical analysis to almost ten thousand candidates reveals an unprecedented degree of tunability and several thousand potential fission-capable candidates, while clearly demonstrating the relationship between triplet aromaticity and singlet-triplet energy gap, confirming this novel strategy for manipulating the exchange energy in organic materials.

Triperyleno[3,3,3]propellane triimides: achieving a new generation of quasi-D 3h symmetric nanostructures in organic electronics

Rigid three-dimensional (3D) polycyclic aromatic hydrocarbons (PAHs), in particular 3D nanographenes, have garnered interest due to their potential use in semiconductor applications and as models to study through-bond and through-space electronic interactions. Herein we report the development of a novel 3D-symmetric rylene imide building block, triperyleno[3,3,3]propellane triimides (6), that possesses three perylene monoimide subunits fused on a propellane. This building block shows several promising characteristics, including high solubility, large π-surfaces, electron-accepting capabilities, and a variety of reactive sites. Further, the building block is compatible with different reactions to readily yield quasi-D 3h symmetric nanostructures (9, 11, and 13) of varied chemistries. For the 3D nanostructures we observed red-shift absorption maxima and amplification of the absorption coefficients when compared to the individual subunits, indicating intramolecular electronic coupling among the subunits. In addition, the microplates of 9 exhibit comparable mobilities in different directions in the range of 10-3 cm2 V-1 s-1, despite the rather limited intermolecular overlap of the π-conjugated moieties. These findings demonstrate that these quasi-D 3h symmetric rylene imides have potential as 3D nanostructures for a range of materials applications, including in organic electronic devices.

Near-Infrared-Absorbing Indolizine-Porphyrin Push-Pull Dye for Dye-Sensitized Solar Cells

Porphyrins are attractive chromophores for application in dye-sensitized solar cells (DSCs), as judicious tuning of donor-acceptor properties can enable excellent near-infrared (NIR) absorption and exceptional device performance. Here, we report a porphyrin-based dye (SM85) conjugated to the planar strong electron donor, indolizine, designed to extend absorption further into the NIR region by inducing π-π interactions such as head-to-tail dye aggregation. The optoelectronic consequences of indolizine incorporation in SM85 include raising the ground-state oxidation potential and broadening and red-shifting ultraviolet-visible-NIR absorptions, along with increased molar absorptivity when compared to the dye SM315. Density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations confirm the push-pull character of SM85, which features an overlap of frontier occupied and unoccupied orbitals. Steady-state spectrophotometric analyses reveal the presence of solution aggregates via absorption and emission spectroscopies. Aggregate modes were probed by DFT and TD-DFT analyses, and plausible models are presented. SM85-based DSC devices demonstrate a 5.7% power conversion efficiency (PCE) at full sun (7.4% PCE at 10% sun) with an exceptional improvement to the incident photon-to-current conversion onset at ∼850 nm. Current dynamics measurements, time-correlated single photon counting, and computational analyses are used to better understand device performances. This study puts forward a novel intramolecular charge-transfer porphyrin system with a dramatic shift into the NIR region, as is needed for nonprecious metal-based sensitizers, and provides an example of controlled, donor-acceptor-mediated aggregation as a complementary strategy to traditional donor-acceptor modifications to single-molecule π-systems in accessing enhancements in long wavelength light harvesting in molecular-based optoelectronic devices.

Oxidation Pathways Involving a Sulfide-Endcapped Donor-Acceptor-Donor π-Conjugated Molecule and Antimony(V) Chloride

The oxidation pathways and products of a discrete, sulfide-endcapped donor-acceptor-donor (D/A/D) molecule, namely, propylenedioxythiophene-benzothiadiazole-propylenedioxythiophene, are investigated. The electrochemical and chemical oxidations proceed by two distinct routes. Specifically, electrochemical oxidation undergoes a sequential two-step, one-electron (1e^-) oxidation route with a 117 mV difference between consecutive half-wave potentials. In contrast, chemical oxidation by antimony(V) chloride (SbCl₅) causes the generation of four different oxidized species: (a) the 1e^- oxidation state, (b) a decomposition product, (c) the 2e^- oxidation state, and (d) a chloride adduct of the 2e^- oxidation state. The decomposition product is generated by the reaction of the 1e^- oxidation state with residual water, resulting in nucleophilic aromatic substitution at the sulfide group terminal positions. This reaction leads to the formation of a 2e^- oxidized, oxygen atom (ketone) terminated decomposed molecule. The chloride adduct is determined to be produced by electrophilic chloronium ion (2e^-) oxidation by the SbCl₄⁺ complex, which is a product of SbCl₅ ligand disproportionation. The formation of the 2e^- oxidized chlorine adduct shows to be linearly dependent on the molarity of SbCl₅ in dichloromethane, giving new insight into the concentration dependent reactivity of SbCl₅ as a 2e^- oxidant. The electronic, optical, and magnetic properties and geometric structures of the 1e^- and 2e^- oxidized hexachloroantimonate salts are fully characterized by a combination of electrochemistry, X-ray crystallography, UV-vis-NIR, electron paramagnetic resonance, NMR spectroscopies, and density functional theory calculations. The aim of this study is to provide a thorough understanding of the redox pathways of a D/A/D π-conjugated organic molecule for potential application in organic electrochromic devices.

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.

Magnetic ordering in a vanadium-organic coordination polymer using a pyrrolo[2,3-d:5,4-d']bis(thiazole)-based ligand

Here we present the synthesis and characterization of a hybrid vanadium-organic coordination polymer with robust magnetic order, a Curie temperature T C of ∼110 K, a coercive field of ∼5 Oe at 5 K, and a maximum mass magnetization of about half that of the benchmark ferrimagnetic vanadium(tetracyanoethylene)∼2 (V·(TCNE)∼2). This material was prepared using a new tetracyano-substituted quinoidal organic small molecule 7 based on a tricyclic heterocycle 4-hexyl-4H-pyrrolo[2,3-d:5,4-d']bis(thiazole) (C6-PBTz). Single crystal X-ray diffraction of the 2,6-diiodo derivative of the parent C6-PBTz, showed a disordered hexyl chain and a nearly linear arrangement of the substituents in positions 2 and 6 of the tricyclic core. Density functional theory (DFT) calculations indicate that C6-PBTz-based ligand 7 is a strong acceptor with an electron affinity larger than that of TCNE and several other ligands previously used in molecular magnets. This effect is due in part to the electron-deficient thiazole rings and extended delocalization of the frontier molecular orbitals. The ligand detailed in this study, a representative example of fused heterocycle aromatic cores with extended π conjugation, introduces new opportunities for structure-magnetic-property correlation studies where the chemistry of the tricyclic heterocycles can modulate the electronic properties and the substituent at the central N-position can vary the spatial characteristics of the magnetic polymer.

Bond Ellipticity Alternation: An Accurate Descriptor of the Nonlinear Optical Properties of π-Conjugated Chromophores

Well-defined structure-property relationships offer a conceptual basis to afford a priori design principles to develop novel π-conjugated molecular and polymer materials for nonlinear optical (NLO) applications. Here, we introduce the bond ellipticity alternation (BEA) as a robust parameter to assess the NLO characteristics of organic chromophores and illustrate its effectiveness in the case of streptocyanines. BEA is based on the symmetry of the electron density, a physical observable that can be determined from experimental X-ray electron densities or from quantum-chemical calculations. Through comparisons to the well-established bond-length alternation and π-bond order alternation parameters, we demonstrate the generality of BEA to foreshadow NLO characteristics and underline that, in the case of large electric fields, BEA is a more reliable descriptor. Hence, this study introduces BEA as a prominent descriptor of organic chromophores of interest for NLO applications.

Correction: An unsymmetrical non-fullerene acceptor: synthesis via direct heteroarylation, self-assembly, and utility as a low energy absorber in organic photovoltaic cells

Correction for 'An unsymmetrical non-fullerene acceptor: synthesis via direct heteroarylation, self-assembly, and utility as a low energy absorber in organic photovoltaic cells' by Abby-Jo Payne et al., Chem. Commun., 2017, 53, 10168-10171.

Unusual Electronic Structure of the Donor-Acceptor Cocrystal Formed by Dithieno[3,2-a:2',3'-c]phenazine and 7,7,8,8-Tetracyanoquinodimethane

Mixed cocrystals derived from electron-rich donor (D) and electron-deficient acceptor (A) molecules showcase electronic, optical, and magnetic properties of interest for a wide range of applications. We explore the structural and electronic properties of a cocrystal synthesized from dithieno[3,2-a:2',3'-c]phenazine (DTPhz) and 7,7,8,8-tetracyanoquinodimethane (TCNQ), which has a mixed-stack packing arrangement of the (π-electronic) face-to-face stacks in a 2:1 D:A stoichiometry. Density functional theory investigations reveal that the primary electronic characteristics of the cocrystal are not determined by electronic interactions along the face-to-face stacks, but rather they are characterized by stronger electronic interactions orthogonal to these stacks that follow the edge-to-edge donor-donor or acceptor-acceptor contacts. These distinctive electronic characteristics portend semiconducting properties that are unusual for semiconducting mixed cocrystals and suggest further potential to design organic semiconductors with orthogonal transport characteristics for different charge carriers.

An unsymmetrical non-fullerene acceptor: synthesis via direct heteroarylation, self-assembly, and utility as a low energy absorber in organic photovoltaic cells

This study reports on the design and synthesis of an unsymmetrical π-conjugated organic molecule composed of perylene diimide, thienyl diketopyrrolopyrrole, and indoloquinoxaline pieced together using direct heteroarylation.

Polarization Energies at Organic-Organic Interfaces: Impact on the Charge Separation Barrier at Donor-Acceptor Interfaces in Organic Solar Cells

We probe the energetic landscape at a model pentacene/fullerene (C60) interface to investigate the interactions between positive and negative charges, which are critical to the processes of charge separation and recombination in organic solar cells. Using a polarizable force field, we find that polarization energy, i.e., the stabilization a charge feels due to its environment, is larger at the interface than in the bulk for both a positive and a negative charge. The combination of the charge being more stabilized at the interface and the Coulomb attraction between the charges results in a barrier to charge separation at the pentacene/C60 interface that can be in excess of 0.7 eV for static configurations of the donor and acceptor locations. However, the impact of molecular motions, i.e., the dynamics, at the interface at room temperature results in a distribution of polarization energies and in charge separation barriers that can be significantly reduced. The dynamic nature of the interface is thus critical, with the polarization energy distributions indicating that sites along the interface shift in time between favorable and unfavorable configurations for charge separation.

Impact of Molecular Orientation and Packing Density on Electronic Polarization in the Bulk and at Surfaces of Organic Semiconductors

The polarizable environment surrounding charge carriers in organic semiconductors impacts the efficiency of the charge transport process. Here, we consider two representative organic semiconductors, tetracene and rubrene, and evaluate their polarization energies in the bulk and at the organic-vacuum interface using a polarizable force field that accounts for induced-dipole and quadrupole interactions. Though both oligoacenes pack in a herringbone motif, the tetraphenyl substituents on the tetracene backbone of rubrene alter greatly the nature of the packing. The resulting change in relative orientations of neighboring molecules is found to reduce the bulk polarization energy of holes in rubrene by some 0.3 eV when compared to tetracene. The consideration of model organic-vacuum interfaces highlights the significant variation in the electrostatic environment for a charge carrier at a surface although the net change in polarization energy is small; interestingly, the environment of a charge even just one layer removed from the surface can be viewed already as representative of the bulk. Overall, it is found that in these herringbone-type layered crystals the polarization energy has a much stronger dependence on the intralayer packing density than interlayer packing density.

Indacenodibenzothiophenes: synthesis, optoelectronic properties and materials applications of molecules with strong antiaromatic character

Indeno[1,2-b]fluorenes (IFs), while containing 4n π-electrons, are best described as two aromatic benzene rings fused to a weakly paratropic s-indacene core. In this study, we find that replacement of the outer benzene rings of an IF with benzothiophenes allows the antiaromaticity of the central s-indacene to strongly reassert itself. Herein we report a combined synthetic, computational, structural, and materials study of anti- and syn-indacenodibenzothiophenes (IDBTs). We have developed an efficient and scalable synthesis for preparation of a series of aryl- and ethynyl-substituted IDBTs. NICS-XY scans and ACID calculations reveal an increasingly antiaromatic core from [1,2-b]IF to anti-IDBT, with syn-IDBT being nearly as antiaromatic as the parent s-indacene. As an initial evaluation, the intermolecular electronic couplings and electronic band structure of a diethynyl anti-IDBT derivative reveal the potential for hole and / or electron transport. OFETs constructed using this molecule show the highest hole mobilities yet achieved for a fully conjugated IF derivative.

Intrinsic Properties of Two Benzodithiophene-Based Donor--Acceptor Copolymers Used in Organic Solar Cells: A Quantum-Chemical Approach

Conjugated donor-acceptor (D-A) copolymers show tremendous promise as active components in thin-film organic bulk heterojunction solar cells and transistors, as appropriate combinations of D-A units enable regulation of the intrinsic electronic and optical properties of the polymer. Here, the structural, electronic, and optical properties of two D-A copolymers that make use of thieno[3,4-c]pyrrole-4,6-dione as the acceptor and differ by their donor unit-benzo[1,2-b:4,5-b']dithiophene (BDT) vs the ladder-type heptacyclic benzodi(cyclopentadithiophene)-are compared using density functional theory methods. Our calculations predict some general similarities, although the differences in the donor structures lead also to clear differences. The extended conjugation of the stiff ladder-type donor destabilizes both the highest occupied and lowest unoccupied molecular orbital energies of the ladder copolymer and results in smaller gap energies compared to its smaller counterpart. However, more significant charge transfer nature is predicted for the smaller BDT-based copolymer by natural transition orbitals than for the ladder copolymer. That is, the influence of the acceptor on the copolymer properties is "diluted" to some extent by the already extended conjugation of the ladder-type donor. Thus, the use of stronger acceptor units with the ladder-type donors would benefit the future design of new D-A copolymers.

Strain effects on the work function of an organic semiconductor

Establishing fundamental relationships between strain and work function (WF) in organic semiconductors is important not only for understanding electrical properties of organic thin films, which are subject to both intrinsic and extrinsic strains, but also for developing flexible electronic devices. Here we investigate tensile and compressive strain effects on the WF of rubrene single crystals. Mechanical strain induced by thermal expansion mismatch between the substrate and rubrene is quantified by X-ray diffraction. The corresponding WF change is measured by scanning Kelvin probe microscopy. The WF of rubrene increases (decreases) significantly with in-plane tensile (compressive) strain, which agrees qualitatively with density functional theory calculations. An elastic-to-plastic transition, characterized by a steep rise of the WF, occurs at ∼0.05% tensile strain along the rubrene π-stacking direction. The results provide the first concrete link between mechanical strain and WF of an organic semiconductor and have important implications for understanding the connection between structural and electronic disorder in soft organic electronic materials.