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.

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.

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.

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 Ea 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 Ea 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 FeCl3 or NOBF4 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.

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 nature of these responses is not only critical for determining the 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, slip-stack, or herringbone (HB) packing configurations; the π-conjugated backbones considered here are pentacene and anthradithiophene. For comparison, HB-packed (unsubstituted) pentacene is also included. Density functional theory and G0W0 (single-shot Green's function G and/or screened Coulomb function W) electronic band structures, G0W0-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. Furthermore, we evaluate the probability of singlet fission in these materials through energy conservation relationships.

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 CO2 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 Eappl = -1.8 V vs Fc+/0, Re(bpy-C2/4/6-NPDI) all achieved comparable activity (TONco ∼ 25) and Faradaic efficiency (FEco ∼ 94%). Under identical CPE conditions, the standard catalyst Re(dmbpy) was inactive for electrocatalytic CO2 reduction; only at Eappl = -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 CO2 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 CO2 conversion.

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.

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.

Rational Functionalization of a C70 Buckybowl To Enable a C70:Buckybowl Cocrystal for Organic Semiconductor Applications.

Fullerene fragments, referred to as buckybowls, are garnering interest due to their distinctive molecular shapes and optoelectronic properties. Here, we report the synthesis and characterization of a novel C70 subunit, diindeno[4,3,2,1-fghi:4',3',2',1'-opqr]perylene, that is substituted with either triethylsilyl(TES)-ethynyl or 2,4,6-triisopropylphenyl groups at the meta-positions. The resulting compounds (1 and 2) display a bowl-to-bowl inversion at room temperature. Notably, the substituent groups on the meta-positions alter both the geometric and the electronic properties as well as the crystal packing of the buckybowls. In contrast to the 2,4,6-triisopropylphenyl groups in 2, the TES-ethynyl groups in 1 lead to enhanced bond length alternation, resulting in weaker aromaticity of the six-membered rings of the buckybowl skeleton. 1 forms one-dimensional (1D) concave-in-convex stacking columns, and when 1 is blended with C70, the buckybowls encapsulate C70 and result in two-dimensional cocrystals. Organic field-effect transistor (OFET) measurements demonstrate that 1 displays a hole mobility of 0.31 cm2 V-1 s-1, and the 1-C70 cocrystal exhibits ambipolar transport characteristics with electron and hole mobilities approaching 0.40 and 0.07 cm2 V-1 s-1, respectively. This work demonstrates the potential of buckybowls for the development of organic semiconductors.

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.

Crossover from band-like to thermally activated charge transport in organic transistors due to strain-induced traps.

The temperature dependence of the charge-carrier mobility provides essential insight into the charge transport mechanisms in organic semiconductors. Such knowledge imparts critical understanding of the electrical properties of these materials, leading to better design of high-performance materials for consumer applications. Here, we present experimental results that suggest that the inhomogeneous strain induced in organic semiconductor layers by the mismatch between the coefficients of thermal expansion (CTE) of the consecutive device layers of field-effect transistors generates trapping states that localize charge carriers. We observe a universal scaling between the activation energy of the transistors and the interfacial thermal expansion mismatch, in which band-like transport is observed for similar CTEs, and activated transport otherwise. Our results provide evidence that a high-quality semiconductor layer is necessary, but not sufficient, to obtain efficient charge-carrier transport in devices, and underline the importance of holistic device design to achieve the intrinsic performance limits of a given organic semiconductor. We go on to show that insertion of an ultrathin CTE buffer layer mitigates this problem and can help achieve band-like transport on a wide range of substrate platforms.

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.

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

For organic semiconductors, the solid-state packing of the π-conjugated molecules or polymers dictate the material electronic, optical, and mechanical characteristics. Combinations of solution and solid-state investigations are often used to establish structure-function relationships, though these connections are often loosely correlated, and experiments in different laboratories can lead to widely variable interpretations. Hence, there remains a need to develop a deeper, more robust understanding of the connections between molecular and polymer chemistry, structure, processing, solid-state order, and materials properties to enable judicious materials design principles. Towards this goal, we employ fully-atomistic molecular dynamics (MD) simulations of poly[4,8-bis(5-(2-ethylhexyl)thiophen-2-yl)benzo[1,2-b;4,5-b']dithiophene-2,6-diyl-alt-(4-(2-ethylhexyl)-3-fluorothieno[3,4-b]thiophene-)-2-carboxylate-2-6-diyl] (PTB7), a donor-acceptor copolymer that has been widely investigated in the organic solar cell literature, to unravel some of these associations. The MD simulations make use of polymer lengths (molecular weights) and solution concentrations that are consistent with those used in experiment, allowing for a detailed picture to arise as to how variations in the polymer environment can direct polymer structure. Comparisons between experiment and theory suggest that processing history can be an important factor in the polymer structures presumed experimentally that are used to interpret optical and electronic responses. The results of these simulations provide specific information into the behavior of PTB7 under different conditions, and showcase how atomistic MD simulations that approach experimentally relevant sizes can be used to develop broader chemical insight that can aid in the design, processing, and characterization of polymer-based organic semiconductors.

Presence of Short Intermolecular Contacts Screens for Kinetic Stability in Packing Polymorphs.

Polymorphism is pervasive in molecular solids. While computational predictions of the molecular polymorphic landscape have improved significantly, identifying which polymorphs are preferentially accessed and experimentally stable remains a challenge. We report a framework that correlates short intermolecular contacts with polymorphic stability. The presence of short contacts between neighboring molecules prevents structural rearrangement and stabilizes the packing arrangement, even when the stabilized polymorph is not enthalpically favored. In the absence of such intermolecular short contacts, the molecules have added degrees of freedom for structural rearrangement, and solid-solid polymorphic transformations occur readily. Starting with a series of core-halogenated naphthalene tetracarboxylic diimides, we establish this framework with the packing polymorphs of more than 20 compounds, ranging from molecular semiconductors to pharmaceutics and biological building blocks. This framework, widely applicable across molecular solids, can help refine computational predictions by identifying the polymorphs that are kinetically stable.

Beyond the Hammett Effect: Using Strain to Alter the Landscape of Electrochemical Potentials.

The substitution of sterically bulky groups at precise locations along the periphery of fused-ring aromatic systems is demonstrated to increase electrochemical oxidation potentials by preventing relaxation events in the oxidized state. Phenothiazines, which undergo significant geometric relaxation upon oxidation, are used as fused-ring models to showcase that electron-donating methyl groups, which would generally be expected to lower oxidation potential, can lead to increased oxidation potentials when used as the steric drivers. Reduction events remain inaccessible through this molecular design route, a critical characteristic for electrochemical systems where high oxidation potentials are required and in which reductive decomposition must be prevented, as in high-voltage lithium-ion batteries. This study reveals a new avenue to alter the redox characteristics of fused-ring systems that find wide use as electroactive elements across a number of developing technologies.

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.

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

Rubrene is one of the most studied molecular semiconductors; its chemical structure consists of a tetracene backbone with four phenyl rings appended to the two central fused rings. Derivatization of these phenyl rings can lead to two very different solid-state molecular conformations and packings: One in which the tetracene core is planar and there exists substantive overlap among neighboring π-conjugated backbones; and another where the tetracene core is twisted and the overlap of neighboring π-conjugated backbones is completely disrupted. State-of-the-art electronic structure calculations show for all isolated rubrene derivatives that the twisted conformation is more favorable (by -1.7 to -4.1 kcal mol(-1)), which is a consequence of energetically unfavorable exchange-repulsion interactions among the phenyl side groups. Calculations based on available crystallographic structures reveal that planar conformations of the tetracene core in the solid state result from intermolecular interactions that can be tuned through well-chosen functionalization of the phenyl side groups and lead to improved intermolecular electronic couplings. Understanding the interplay of these intramolecular and intermolecular interactions provides insight into how to chemically modify rubrene and similar molecular semiconductors to improve the intrinsic materials electronic properties.

Nonlinear Optical Properties of X(C6H5)4 (X = B(-), C, N(+), P(+)): A New Class of Molecules with a Negative Third-Order Polarizability.

Organic π-conjugated materials have been widely used for a variety of nonlinear optical (NLO) applications. Molecules with negative real components Re(γ) of the third-order polarizability, which leads to nonlinear refraction in macroscopic systems, have important benefits for several NLO applications. However, few organic systems studied to date have negative Re(γ) in the long wavelength limit, and all inorganic materials show positive nonlinear refraction in this limit. Here, we introduce a new class of molecules of the form X(C6H5)4, where X = B(-), C, N(+), and P(+), that have negative Re(γ). The molecular mechanism for the NLO properties in these systems is very different from those in typical linear conjugated systems: These systems have a band of excited states involving single-electron excitations within the π-system, several of which have significant coupling to the ground state. Thus, Re(γ) cannot be understood in terms of a simplified essential-state model and must be analyzed in the context of the full sum-over-states expression. Although Re(γ) is significantly smaller than that of other commonly studied NLO chromophores, the introduction of a new molecular architecture offering the potential for a negative Re(γ) introduces new avenues of molecular design for NLO applications.