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

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.

COMPAS-3: a dataset of peri-condensed polybenzenoid hydrocarbons

We introduce the third installment of the COMPAS Project - a COMputational database of Polycyclic Aromatic Systems, focused on peri-condensed polybenzenoid hydrocarbons. In this installment, we develop two datasets containing the optimized ground-state structures and a selection of molecular properties of ∼39k and ∼9k peri-condensed polybenzenoid hydrocarbons (at the GFN2-xTB and CAM-B3LYP-D3BJ/cc-pvdz//CAM-B3LYP-D3BJ/def2-SVP levels, respectively). The manuscript details the enumeration and data generation processes and describes the information available within the datasets. An in-depth comparison between the two types of computation is performed, and it is found that the geometrical disagreement is maximal for slightly-distorted molecules. In addition, a data-driven analysis of the structure-property trends of peri-condensed PBHs is performed, highlighting the effect of the size of peri-condensed islands and linearly annulated rings on the HOMO-LUMO gap. The insights described herein are important for rational design of novel functional aromatic molecules for use in, e.g., organic electronics. The generated datasets provide a basis for additional data-driven machine- and deep-learning studies in chemistry.

COMPAS-2: a dataset of cata-condensed hetero-polycyclic aromatic systems

Polycyclic aromatic systems are highly important to numerous applications, in particular to organic electronics and optoelectronics. High-throughput screening and generative models that can help to identify new molecules to advance these technologies require large amounts of high-quality data, which is expensive to generate. In this report, we present the largest freely available dataset of geometries and properties of cata-condensed poly(hetero)cyclic aromatic molecules calculated to date. Our dataset contains ~500k molecules comprising 11 types of aromatic and antiaromatic building blocks calculated at the GFN1-xTB level and is representative of a highly diverse chemical space. We detail the structure enumeration process and the methods used to provide various electronic properties (including HOMO-LUMO gap, adiabatic ionization potential, and adiabatic electron affinity). Additionally, we benchmark against a ~50k dataset calculated at the CAM-B3LYP-D3BJ/def2-SVP level and develop a fitting scheme to correct the xTB values to higher accuracy. These new datasets represent the second installment in the COMputational database of Polycyclic Aromatic Systems (COMPAS) Project.

Data-Driven Discovery of Organic Electronic Materials Enabled by Hybrid Top-Down/Bottom-Up Design

The high-throughput exploration and screening of molecules for organic electronics involves either a 'top-down' curation and mining of existing repositories, or a 'bottom-up' assembly of user-defined fragments based on known synthetic templates. Both are time-consuming approaches requiring significant resources to compute electronic properties accurately. Here, 'top-down' is combined with 'bottom-up' through automatic assembly and statistical models, thus providing a platform for the fragment-based discovery of organic electronic materials. This study generates a top-down set of 117K synthesized molecules containing structures, electronic and topological properties and chemical composition, and uses them as building blocks for bottom-up design. A tool is developed to automate the coupling of these building blocks at their C(sp2/sp)-H bonds, providing a fundamental link between the two dataset construction philosophies. Statistical models are trained on this dataset and a subset of resulting top-down/bottom-up compounds, enabling on-the-fly prediction of ground and excited state properties with high accuracy across organic compound space. With access to ab initio-quality optical properties, this bottom-up pipeline may be applied to any materials design campaign using existing compounds as building blocks. To illustrate this, over a million molecules are screened for singlet fission. tThe leading candidates provide insight into the features promoting this multiexciton-generating process.

Benchmark Data Set of Crystalline Organic Semiconductors

This work reports a Benchmark Data set of Crystalline Organic Semiconductors to test calculations of the structural and electronic properties of these materials in the solid state. The data set contains 67 crystals consisting of mostly rigid molecules with a single dominant conformer, covering the majority of known structural types. The experimental crystal structure is available for the entire data set, whereas zero-temperature unit cell volume can be reliably estimated for a subset of 28 crystals. Using this subset, we benchmark r2SCAN-D3 and PBE-D3 density functionals. Then, for the entire data set, we benchmark approximate density functional theory (DFT) methods, including GFN1-xTB and DFTB3(3ob-3-1), with various dispersion corrections against r2SCAN-D3. Our results show that r2SCAN-D3 geometries are accurate within a few percent, which is comparable to the statistical uncertainty of experimental data at a fixed temperature, but the unit cell volume is systematically underestimated by 2% on average. The several times faster PBE-D3 provides an unbiased estimate of the volume for all systems except for molecules with highly polar bonds, for which the volume is substantially overestimated in correlation with the underestimation of atomic charges. Considered approximate DFT methods are orders of magnitude faster and provide qualitatively correct but overcompressed crystal structures unless the dispersion corrections are fitted by unit cell volume.

Conjugated Polymer Process Ontology and Experimental Data Repository for Organic Field-Effect Transistors

Polymer-based semiconductors and organic electronics encapsulate a significant research thrust for informatics-driven materials development. However, device measurements are described by a complex array of design and parameter choices, many of which are sparsely reported. For example, the mobility of a polymer-based organic field-effect transistor (OFET) may vary by several orders of magnitude for a given polymer as a plethora of parameters related to solution processing, interface design/surface treatment, thin-film deposition, postprocessing, and measurement settings have a profound effect on the value of the final measurement. Incomplete contextual, experimental details hamper the availability of reusable data applicable for data-driven optimization, modeling (e.g., machine learning), and analysis of new organic devices. To curate organic device databases that contain reproducible and findable, accessible, interoperable, and reusable (FAIR) experimental data records, data ontologies that fully describe sample provenance and process history are required. However, standards for generating such process ontologies are not widely adopted for experimental materials domains. In this work, we design and implement an object-relational database for storing experimental records of OFETs. A data structure is generated by drawing on an international standard for batch process control (ISA-88) to facilitate the design. We then mobilize these representative data records, curated from the literature and laboratory experiments, to enable data-driven learning of process-structure-property relationships. The work presented herein opens the door for the broader adoption of data management practices and design standards for both the organic electronics and the wider materials community.

Machine learning of atomic force microscopy images of organic solar cells

The bulk heterojunction structures of organic photovoltaics (OPVs) have been overlooked in their machine learning (ML) approach despite their presumably significant impact on power conversion efficiency (PCE). In this study, we examined the use of atomic force microscopy (AFM) images to construct an ML model for predicting the PCE of polymer : non-fullerene molecular acceptor OPVs. We manually collected experimentally observed AFM images from the literature, applied data curing and performed image analyses (fast Fourier transform, FFT; gray-level co-occurrence matrix, GLCM; histogram analysis, HA) and ML linear regression. The accuracy of the model did not considerably improve even by including AFM data in addition to the chemical structure fingerprints, material properties and process parameters. However, we found that a specific spatial wavelength of FFT (40-65 nm) significantly affects PCE. The GLCM and HA methods, such as homogeneity, correlation and skewness expand the scope of image analysis and artificial intelligence in materials science research fields.

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.

Identification via Virtual Screening of Emissive Molecules with a Small Exciton-Vibration Coupling for High Color Purity and Potential Large Exciton Delocalization

A sequence of quantum chemical computations of increasing accuracy was used in this work to identify molecules with small exciton reorganization energy (exciton-vibration coupling), of interest for light emitting devices and coherent exciton transport, starting from a set of ∼4500 known molecules. We validated an approximate computational approach based on single-point calculations of the force in the excited state, which was shown to be very efficient in identifying the most promising candidates. We showed that a simple descriptor based on the bond order could be used to find molecules with potentially small exciton reorganization energies without performing excited state calculations. A small set of chemically diverse molecules with a small exciton reorganization energy was analyzed in greater detail to identify common features leading to this property. Many such molecules display an A-B-A structure where the bonding/antibonding patterns in the fragments A are similar in HOMO and LUMO. Another group of molecules with small reorganization energy displays instead HOMO and LUMO with a strong nonbonding character.

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.

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.

Delving into guest-free and He-filled sI and sII clathrate hydrates: a first-principles computational study

The dynamics of the formation of a specific clathrate hydrate as well as its thermodynamic transitions depend on the interactions between the trapped molecules and the host water lattice. The molecular-level understanding of the different underlying processes benefits not only the description of the properties of the system, but also allows the development of multiple technological applications such as gas storage, gas separation, energy transport, etc. In this work we investigate the stability of periodic crystalline structures, such as He@sI and He@sII clathrate hydrates by first-principles computations. We consider such host water networks interacting with a guest He atom using selected density functional theory approaches, in order to explore the effects on the encapsulation of a light atom in the sI/sII crystals, by deriving all energy components (guest-water, water-water, guest-guest). Structural properties and energies were first computed by structural relaxations of the He-filled and empty sI/sII unit cells, yielding lattice and compressibility parameters comparable to experimental and theoretical values available for those hydrates. According to the results obtained, the He enclathration in the sI/sII unit cells is a stabilizing process, and both He@sI and He@sII clathrates, considering single cage occupancy, are predicted to be stable whatever the XDM or D4 dispersion correction applied. Our results further reveal that despite the weak underlying interactions the He encapsulation has a rather notable effect on both lattice parameters and energetics, with the He@sII being the most energetically favorable in accord with recent experimental observations.

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.