Read between the Molecules: Computational Insights into Organic Semiconductors

Data-driven tailoring of molecular dipole polarizability and frontier orbital energies in chemical compound space

Understanding correlations - or lack thereof - between molecular properties is crucial for enabling fast and accurate molecular design strategies. In this contribution, we explore the relation between two key quantities describing the electronic structure and chemical properties of molecular systems: the energy gap between the frontier orbitals and the dipole polarizability. Based on the recently introduced QM7-X dataset, augmented with accurate molecular polarizability calculations as well as analysis of functional group compositions, we show that polarizability and HOMO-LUMO gap are uncorrelated when considering sufficiently extended subsets of the chemical compound space. The relation between these two properties is further analyzed on specific examples of molecules with similar composition as well as homooligomers. Remarkably, the freedom brought by the lack of correlation between molecular polarizability and HOMO-LUMO gap enables the design of novel materials, as we demonstrate on the example of organic photodetector candidates.

Four Slip-Stacked Arrangements, Three Types of Photophysics: Crystal Structure and Solid-State Fluorescence of 3,6-Diaryl Substituted Furo[3,4-c]furanone Polymorphs and Regioisomers

Six symmetrical 3,6-diaryl (aryl=phenyl, 2-, 3- and 4-tolyl, 2,4- and 3,5-xylyl) substituted furo[3,4-c]furanones (DFF) were synthesized. The computational analysis, based on density functional theory, found eight possible centrosymmetrical slipped π-stack arrangements, formed according to electron repulsion minimization principle, as for previously reported for π-isoelectronic diketopyrrolopyrroles (DPP). One of these slipped stack arrangements was found to form infinite columns in the crystals of a new polymorph of parent phenyl derivative (with centre-to-centre distance CC=6.975 Å), other three types of stacks were found for 3-tolyl (CC=6.153 Å), 4-tolyl (CC=3.849 Å) and 2,4-xylyl (CC=4.856 Å) derivatives by single crystal X-ray diffractometry. All six derivatives show intense solution fluorescence in blue/green region, with a maximum driven entirely by a number and position of methyl substituents on phenyl rings. On the other hand, the solid-state fluorescence from yellow over orange to red is observed only for four derivatives and its presence/absence, spectral position and vibronic structure is driven exclusively by the slips in π-stacks (with interplanar distance always less than 3.5 Å) of almost planar DFF molecules, resulting in J-type emission, H-type excimer-like emission and H-type quenching.

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.

Recent Research Progress of Organic Small-Molecule Semiconductors with High Electron Mobilities

Organic electronics has made great progress in the past decades, which is inseparable from the innovative development of organic electronic devices and the diversity of organic semiconductor materials. It is worth mentioning that both of these great advances are inextricably linked to the development of organic high-performance semiconductor materials, especially the representative n-type organic small-molecule semiconductor materials with high electron mobilities. The n-type organic small molecules have the advantages of simple synthesis process, strong intermolecular stacking, tunable molecular structure, and easy to functionalize structures. Furthermore, the n-type semiconductor is a remarkable and important component for constructing complementary logic circuits and p-n heterojunction structures. Therefore, n-type organic semiconductors play an extremely important role in the field of organic electronic materials and are the basis for the industrialization of organic electronic functional devices. This review focuses on the modification strategies of organic small molecules with high electron mobility at molecular level, and discusses in detail the applications of n-type small-molecule semiconductor materials with high mobility in organic field-effect transistors, organic light-emitting transistors, organic photodetectors, and gas sensors.

High-throughput property-driven generative design of functional organic molecules

The design of molecules and materials with tailored properties is challenging, as candidate molecules must satisfy multiple competing requirements that are often difficult to measure or compute. While molecular structures produced through generative deep learning will satisfy these patterns, they often only possess specific target properties by chance and not by design, which makes molecular discovery via this route inefficient. In this work, we predict molecules with (Pareto-)optimal properties by combining a generative deep learning model that predicts three-dimensional conformations of molecules with a supervised deep learning model that takes these as inputs and predicts their electronic structure. Optimization of (multiple) molecular properties is achieved by screening newly generated molecules for desirable electronic properties and reusing hit molecules to retrain the generative model with a bias. The approach is demonstrated to find optimal molecules for organic electronics applications. Our method is generally applicable and eliminates the need for quantum chemical calculations during predictions, making it suitable for high-throughput screening in materials and catalyst design.

Polymer Photoelectrodes for Solar Fuel Production: Progress and Challenges

Converting solar energy to fuels has attracted substantial interest over the past decades because it has the potential to sustainably meet the increasing global energy demand. However, achieving this potential requires significant technological advances. Polymer photoelectrodes are composed of earth-abundant elements, e.g. carbon, nitrogen, oxygen, hydrogen, which promise to be more economically sustainable than their inorganic counterparts. Furthermore, the electronic structure of polymer photoelectrodes can be more easily tuned to fit the solar spectrum than inorganic counterparts, promising a feasible practical application. As a fast-moving area, in particular, over the past ten years, we have witnessed an explosion of reports on polymer materials, including photoelectrodes, cocatalysts, device architectures, and fundamental understanding experimentally and theoretically, all of which have been detailed in this review. Furthermore, the prospects of this field are discussed to highlight the future development of polymer photoelectrodes.

Complete Set of Diketopyrrolopyrrole Centrosymmetrical Cofacial Stacked Pairs

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

Synergistic effects of side-functionalization and aza-substitution on the charge transport and optical properties of perylene-based organic materials: a DFT study

The physicochemical properties of organic materials are subject to the chemical structure of the molecular unit and the arrangement of molecules in a crystal.

Is a Single Molecule Sufficient to Determine the Internal Charge Trapping Energy in Crystalline Organic Semiconductors?

In silico diagnoses of charge-transfer efficiencies in organic semiconductors require the accurate computations of transport parameters. We show here that ignoring the molecular packing effects when computing the internal charge trapping energy may cause a severe deviation to the result deduced from the periodic crystal structure. This deviation can reach up to 100 meV in common organic materials. According to the semiclassical Marcus theory, this energy difference can lead to orders of magnitude change in the charge transfer rate. Studying from a total of 45 organic crystals, we find that single-molecule approximation though is adequate for rigid planar molecules yet it fails to describe the trapping energy for molecules that have inter-ring single bond(s) or are subjected to planarity changes during the transition from the isolated state to the embedded state. These results and conclusions may shed light on the removal of theory-experiment discrepancies on charge mobilities and lay the basis for the future fast and precise screening of high-performance organic semiconductors.

Ultrafast estimation of electronic couplings for electron transfer between pi-conjugated organic molecules. II

The development of highly efficient methods for the calculation of electronic coupling matrix elements between the electron donor and acceptor is an important goal in theoretical organic semiconductor research. In Paper I [F. Gajdos, S. Valner, F. Hoffmann, J. Spencer, M. Breuer, A. Kubas, M. Dupuis, and J. Blumberger, J. Chem. Theory Comput. 10, 4653 (2014)], we introduced the analytic overlap method (AOM) for this purpose, which is an ultrafast electronic coupling estimator parameterized to and orders of magnitude faster than density functional theory (DFT) calculations at a reasonably small loss in accuracy. In this work, we reparameterize and extend the AOM to molecules containing nitrogen, oxygen, fluorine, and sulfur heteroatoms using 921 dimer configurations from the recently introduced HAB79 dataset. We find again a very good linear correlation between the frontier orbital overlap, calculated ultrafast in an optimized minimum Slater basis, and DFT reference electronic couplings. The new parameterization scheme is shown to be transferable to sulfur-containing polyaromatic hydrocarbons in experimentally resolved dimeric configurations. Our extension of the AOM enables high-throughput screening of very large databases of chemically diverse organic crystal structures and the application of computationally intense non-adiabatic molecular dynamics methods to charge transport in state-of-the-art organic semiconductors, e.g., non-fullerene acceptors.

HAB79: A new molecular dataset for benchmarking DFT and DFTB electronic couplings against high-level ab initio calculations

A new molecular dataset called HAB79 is introduced to provide ab initio reference values for electronic couplings (transfer integrals) and to benchmark density functional theory (DFT) and density functional tight-binding (DFTB) calculations. The HAB79 dataset is composed of 79 planar heterocyclic polyaromatic hydrocarbon molecules frequently encountered in organic (opto)electronics, arranged to 921 structurally diverse dimer configurations. We show that CASSCF/NEVPT2 with a minimal active space provides a robust reference method that can be applied to the relatively large molecules of the dataset. Electronic couplings are largest for cofacial dimers, in particular, sulfur-containing polyaromatic hydrocarbons, with values in excess of 0.5 eV, followed by parallel displaced cofacial dimers. V-shaped dimer motifs, often encountered in the herringbone layers of organic crystals, exhibit medium-sized couplings, whereas T-shaped dimers have the lowest couplings. DFT values obtained from the projector operator-based diabatization (POD) method are initially benchmarked against the smaller databases HAB11 (HAB7-) and found to systematically improve when climbing Jacob's ladder, giving mean relative unsigned errors (MRUEs) of 27.7% (26.3%) for the generalized gradient approximation (GGA) functional BLYP, 20.7% (15.8%) for hybrid functional B3LYP, and 5.2% (7.5%) for the long-range corrected hybrid functional omega-B97X. Cost-effective POD in combination with a GGA functional and very efficient DFTB calculations on the dimers of the HAB79 database give a good linear correlation with the CASSCF/NEVPT2 reference data, which, after scaling with a multiplicative constant, gives reasonably small MRUEs of 17.9% and 40.1%, respectively, bearing in mind that couplings in HAB79 vary over 4 orders of magnitude. The ab initio reference data reported here are expected to be useful for benchmarking other DFT or semi-empirical approaches for electronic coupling calculations.

Ion-Driven Self-Assembly of Lanthanide Bis-phthalocyaninates into Conductive Quasi-MOF Nanowires: an Approach toward Easily Recyclable Organic Electronics

Controlled self-assembly and rapid disintegration of supramolecular nanowires is potentially useful for ecology-friendly organic electronics. Herein, a novel method exploiting the binding between crown-substituted double-decker lanthanide phthalocyaninates (ML₂, M = Lu, Ce, Tb) and K⁺ ions is applied for the one-step fabrication of macroscopically long conductive one-dimensional quasi-metal-organic frameworks. Their properties are controlled by the size of the lanthanide ion guiding the assembly through either intra- or intermolecular interactions. A LuL₂ linker with a small interdeck distance yields fully conjugated intermolecular-bonded K⁺-LuL₂ nanowires with a thickness of 10-50 nm, a length of up to 50 μm, and a conductivity of up to 11.4 S cm^-1, the highest among them being reported for phthalocyanine assemblies. The large size of CeL₂ and TbL₂ leads to the formation of mixed intra- and intermolecular K⁺-ML₂ phases with poor electric properties. A field-assisted method is developed to deposit aligned conductive K⁺-LuL₂ assemblies on solids. The solid-supported nanowires can be disintegrated into starting components in a good aprotic solvent for further recycling.

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

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

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

Physics-Inspired Structural Representations for Molecules and Materials

The first step in the construction of a regression model or a data-driven analysis, aiming to predict or elucidate the relationship between the atomic-scale structure of matter and its properties, involves transforming the Cartesian coordinates of the atoms into a suitable representation. The development of atomic-scale representations has played, and continues to play, a central role in the success of machine-learning methods for chemistry and materials science. This review summarizes the current understanding of the nature and characteristics of the most commonly used structural and chemical descriptions of atomistic structures, highlighting the deep underlying connections between different frameworks and the ideas that lead to computationally efficient and universally applicable models. It emphasizes the link between properties, structures, their physical chemistry, and their mathematical description, provides examples of recent applications to a diverse set of chemical and materials science problems, and outlines the open questions and the most promising research directions in the field.

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.

Multi-level aggregation of conjugated small molecules and polymers: from morphology control to physical insights

Aggregation of molecules is a multi-molecular phenomenon occurring when two or more molecules behave differently from discrete molecules due to their intermolecular interactions. Moving beyond single molecules, aggregation usually demonstrates evolutive or wholly emerging new functionalities relative to the molecular components. Conjugated small molecules and polymers interact with each other, resulting in complex solution-state aggregates and solid-state microstructures. Optoelectronic properties of conjugated small molecules and polymers are sensitively determined by their aggregation states across a broad range of spatial scales. This review focused on the aggregation ranging from molecular structure, intermolecular interactions, solution-state assemblies, and solid-state microstructures of conjugated small molecules and polymers. We addressed the importance of such aggregation in filling the gaps from the molecular level to device functions and highlighted the multi-scale structures and properties at different scales. From the view of multi-level aggregation behaviors, we divided the whole process from the molecule to devices into several parts: molecular design, solvation, solution-state aggregation, crystal engineering, and solid-state microstructures. We summarized the progress and challenges of relationships between optoelectronic properties and multi-level aggregation. We believe aggregation science will become an interdisciplinary research field and serves as a general platform to develop future materials with the desired functions.

Active discovery of organic semiconductors

The versatility of organic molecules generates a rich design space for organic semiconductors (OSCs) considered for electronics applications. Offering unparalleled promise for materials discovery, the vastness of this design space also dictates efficient search strategies. Here, we present an active machine learning (AML) approach that explores an unlimited search space through consecutive application of molecular morphing operations. Evaluating the suitability of OSC candidates on the basis of charge injection and mobility descriptors, the approach successively queries predictive-quality first-principles calculations to build a refining surrogate model. The AML approach is optimized in a truncated test space, providing deep methodological insight by visualizing it as a chemical space network. Significantly outperforming a conventional computational funnel, the optimized AML approach rapidly identifies well-known and hitherto unknown molecular OSC candidates with superior charge conduction properties. Most importantly, it constantly finds further candidates with highest efficiency while continuing its exploration of the endless design space.

Roughness evolution in strongly interacting donor:acceptor mixtures of molecular semiconductors. An in situ, real-time growth study using x-ray reflectivity

The evolution of surface roughness in binary mixtures of the two molecular organic semiconductors (OSCs) diindenoperylene (DIP) as electron-donor and 1, 3, 4, 5, 7, 8-hexafluoro-tetracyano naphthoquinodimethane (F6TCNNQ) as electron-acceptor is studied. We co-deposit DIP and F6TCNNQ in vacuum with varying relative molar content while keeping a molar excess of DIP in order to produce phase-heterogeneous mixtures. The excess DIP phase segregates in pristine crystallites, whereas the remaining mixed phase is constituted by DIP:F6TCNNQ co-crystallites. We calculate the surface roughness as function of film thickness by modelling x-ray reflectivity data acquired in situ and in real-time during film growth. To model the experimental data, two distinct approaches, namely the kinematic approximation and the Parratt formalism, are applied. A comparative study of surface roughness evolution as function of DIP:F6TCNNQ mixing ratio is carried out implementing the Trofimov growth model within the kinematic approximation. Depending on the thickness regime, mixing ratio-specific trends are identified and discussed. To explain them, a growth mechanism for binary heterogeneous mixtures of strongly interacting OSCs is proposed.

Structure-Property Relationships in Bithiophenes with Hydrogen-Bonded Substituents

The use of crystal engineering to control the supramolecular arrangement of π-conjugated molecules in the solid-state is of considerable interest for the development of novel organic electronic materials. In this study, we investigated the effect of combining of two types of supramolecular interaction with different geometric requirements, amide hydrogen bonding and π-interactions, on the π-overlap between calamitic π-conjugated cores. To this end, we prepared two series of bithiophene diesters and diamides with methylene, ethylene, or propylene spacers between the bithiophene core and the functional groups in their terminal substituents. The hydrogen-bonded bithiophene diamides showed significantly denser packing of the bithiophene cores than the diesters and other known α,ω-disubstituted bithiophenes. The bithiophene packing density reach a maximum in the bithiophene diamide with an ethylene spacer, which had the smallest longitudinal bithiophene displacement and infinite 1D arrays of electronically conjugated, parallel, and almost linear N-H⋅⋅⋅O=C hydrogen bonds. The synergistic hydrogen bonding and π-interactions were attributed to the favorable conformation mechanics of the ethylene spacer and resulted in H-type spectroscopic aggregates in solid-state absorption spectroscopy. These results demonstrate that the optoelectronic properties of π-conjugated materials in the solid-state may be tailored systematically by side-chain engineering, and hence that this approach has significant potential for the design of organic and polymer semiconductors.

Quantifying Planarity in the Design of Organic Electronic Materials

Planarity is essential for many organic electronic materials as it maximizes the intramolecular π-orbital overlap and enables efficient intermolecular interactions through π-stacking. We propose a statistical way of quantifying the planarity of a wide range of conjugated systems. The quantification takes into account all torsional conformations and their relative contribution to the overall structural disorder, through a planarity index ⟨cos² ϕ⟩. The propensity for planarization and the effect of rotational disorder were examined for a series of commonly used building blocks. The application of the analysis to extended conjugated systems and the correlations between the gas-phase ⟨cos² ϕ⟩ and crystallographically observed planarity in the solid state were explored. Our calculations also reveal a previously unrecognized effect of increasing band gap upon planarization for conjugated systems coupling strong electron donor and acceptor units.

Radical cations and dications of bis[1]benzothieno[1,4]thiazine isomers

Radical cations and dications of three isomeric bis[1]benzothieno[1,4]thiazines are (electro)chemically generated, spectroscopically characterized and assigned by (TD)DFT calculations.

Mutually exclusive hole and electron transfer coupling in cross stacked acenes

Acenes in the Greek cross (+) stack orientation exhibit selective hole and electron transfer coupling based on gerade symmetry in frontier molecular orbitals.

Unilaterally Fluorinated Acenes: Synthesis and Solid-State Properties

The rapid development of organic electronics is closely related to the availability of molecular materials with specific electronic properties. Here, we introduce a novel synthetic route enabling a unilateral functionalization of acenes along their long side, which is demonstrated by the synthesis of 1,2,10,11,12,14-hexafluoropentacene (1) and the related 1,2,9,10,11-pentafluorotetracene (2). Quantum chemical DFT calculations in combination with optical and X-ray absorption spectroscopy data indicate that the single-molecule properties of 1 are a connecting link between the organic semiconductor model systems pentacene (PEN) and perfluoropentacene (PFP). In contrast, the crystal structure analysis reveals a different packing motif than for the parent molecules. This can be related to distinct F⋅⋅⋅H interactions identified in the corresponding Hirshfeld surface analysis and also affects solid-state properties such as the exciton binding energy and the sublimation enthalpy.

Designing Organic Semiconductors with Ultrasmall Reorganization Energies: Insights from Molecular Symmetry, Aromaticity and Energy Gap

An intuitive design strategy for organic semiconductors with ultrasmall reorganization energy (λ) is proposed. Learning from a total of 98 molecules condensed by benzene and/or thiophene rings, we find that linear compounds in D_2h symmetry have the smallest λ in each of the three molecular categories (PAHs, thienothiophenes, benzothiophenes). 2D expanded analogues that contain these D_2h building blocks also give unusually small λ (<100 meV). λ of 1D elongated polycyclics show an approximate linear correlation with the ring-averaged HOMA indices and the HOMO-LUMO gaps. Compared to the symmetry principle, the HOMA and energy gap, though much less intuitive to design a priori, provide additional quantitative guidelines to further optimize λ through substitutions, for example, when molecules have the same symmetries. Our results indicate that ring-fused π-conjugates that have narrower HOMO-LUMO gaps and are less aromatic are better candidates to achieve ultrasmall λ.

FB-REDA: fragment-based decomposition analysis of the reorganization energy for organic semiconductors

We present a fragment-based decomposition analysis tool (FB-REDA) for the reorganisation energy (λ). This tool delivers insights on how to rationally design low-λ organic semiconductors. The contribution of the fragment vibrational modes to the reorganization energy is exploited to identity the individual contributions of the molecular building blocks. The usefulness of the approach is demonstrated by offering three strategies to reduce the reorganization energy of a promising dopant-free hole transport material (TPA1PM, λ = 213 meV). A reduction of nearly 50% (TPD3PM, λ = 108 meV) is achieved. The proposed design principles are likely transferable to other organic semiconductors exploiting common molecular building blocks.

Bent-Shaped p-Type Small-Molecule Organic Semiconductors: A Molecular Design Strategy for Next-Generation Practical Applications

Significant progress has been made in both molecular design and fundamental scientific understanding of organic semiconductors (OSCs) in recent years. Suitable charge-carrier mobilities (μ) have been obtained by many high-performance OSCs (μ > 10 cm² V^-1 s^-1), but drawbacks remain, including low solution processability and poor thermal durability. In addition, since aggregation of OSCs involves weak intermolecular interactions, the molecules are perpetually in thermal motion, even in the solid state, which disrupts charge-carrier transport. These issues limit potential applications of OSCs. The present work examines a molecular design for hole-transporting (p-type) OSCs based on the "bent-shaped" geometry with specific molecular orbital configurations, which aims to enhance effective intermolecular orbital overlaps, stabilize crystal phases, suppress detrimental molecular motions in the solid state, and improve solution processability. The results indicated that such OSCs have high μ and suitable solution processability, and are resistant to ambient and thermal conditions, making them suitable for practical applications.

Fluorinated Thiophene-Phenylene Co-Oligomers for Optoelectronic Devices

Organic optoelectronics requires materials combining bright luminescence and efficient ambipolar charge transport. Thiophene-phenylene co-oligomers (TPCOs) are promising highly emissive materials with decent charge-carrier mobility; however, they typically show poor electron injection in devices, which is usually assigned to high energies of their lowest unoccupied molecular orbitals (LUMOs). A widely used approach to lower the frontier orbitals energy levels of a conjugated molecule is its fluorination. In this study, we synthesized three new fluorinated derivatives of one of the most popular TPCOs, 2,2'-(1,4-phenylene)bis[5-phenylthiophene] (PTPTP) and studied them by cyclic voltammetry, absorption, photoluminescence, and Raman spectroscopies. The obtained data reveal a positive effect of fluorination on the optoelectronic properties of PTPTP: LUMO levels are finely tuned, and photoluminescence quantum yield and absorbance are increased. We then grew crystals from fluorinated PTPTPs, resolved their structures, and showed that fluorination dramatically affects the packing motif and facilitates π-stacking. Finally, we fabricated thin-film organic field-effect transistors (OFETs) and demonstrated a strong impact of fluorination on charge injection/transport for both types of charge carriers, namely, electrons and holes. Specifically, balanced ambipolar charge transport and electroluminescence were observed only in the OFET active channel based on the partially fluorinated PTPTP. The obtained results can be extended to other families of conjugated oligomers and highlight the efficiency of fluorination for rational design of organic semiconductors for optoelectronic devices.

Charge mobility calculation of organic semiconductors without use of experimental single-crystal data

Prediction of material properties of newly designed molecules is a long-term goal in organic electronics. In general, it is a difficult problem, because the material properties are dominated by the unknown packing structure. We present a practical method to obtain charge transport properties of organic single crystals, without use of experimental single-crystal data. As a demonstration, we employ the promising molecule C10-DNBDT. We succeeded in quantitative evaluation of charge mobility of the single crystal using our quantum wave-packet dynamical simulation method. Here, the single-crystal data is computationally obtained by searching possible packing structures from structural formula of the molecule. We increase accuracy in identifying the actual crystal structure from suggested ones by using not only crystal energy but also similarity between calculated and experimental powder X-ray diffraction patterns. The proposed methodology can be a theoretical design technique for efficiently developing new high-performance organic semiconductors, since it can estimate the charge transport properties at early stage in the process of material development.

Stereochemistry, Stereodynamics, and Redox and Complexation Behaviors of 2,2'-Diaryl-1,1'-Biazulenes

2,2'-Diaryl-1,1'-biazulenes were synthesized and electronic communication between the azulene subunits was suggested based on redox measurements. The linkage of azulene at the 1-position also appeared to increase the HOMO levels. In addition, cyclic voltammetry measurements of 2-arylazulenes showed a return peak associated with the oxidation, which was not observed for azulene. The stabilization of the single-electron oxidant may be due to the SOMO-HOMO energy inversion phenomenon. X-ray crystallography of the azulene dimers revealed that this species possessed a syn-type structure in which both aryl groups in the 2-positions formed π-stacks. The twisted structure was indicated to be in the (R)- or (S)-configuration for all molecules in the unit cell. Spontaneous resolution was also shown. Furthermore, from the solid circular dichroism (CD) spectral measurements, the relationship between the absolute configuration of the molecules and the CD spectra was determined. A racemization rotational barrier of ca. 27 kcal mol^-1 was calculated. Moreover, the pyridylazulene dimer cyclized upon reaction with PdCl₂ to form a 3 : 3 complex, in which the biazulene units cyclized to give ratios between the (R)- and (S)-forms of either 2 : 1 or 1 : 2.

Towards materials discovery: assays for screening and study of chemical interactions of novel corrosion inhibitors in solution and coatings

Using chemical design, computational modelling, structure–activity and structure–property relationships, with high-throughput solution and coating assays we can rapidly identify prime corrosion inhibitor candidates for further detailed evaluation.

Impact of N-substitution on structural, electronic, optical, and vibrational properties of a thiophene–phenylene co-oligomer

Properties of the organic semiconductors can be finely tuned via changes in their molecular structure.

Structural order and charge transfer in highly strained carbon nanobelts

We present a computational study of the atomic morphology, structural order and charge transfer properties of radially π-conjugated, closed-loop, and highly strained chiral carbon nanobelts (CNBs).

High performance single-crystalline organic field-effect transistors based on molecular-modified dibenzo[a,e]pentalenes derivatives

Modulating the charge transport properties realized by a controllable molecular structure resulted in different packing arrangements.

Machine Learning for Predicting Electron Transfer Coupling

Electron transfer coupling is a critical factor in determining electron transfer rates. This coupling strength can be sensitive to details in molecular geometries, especially intermolecular configurations. Thus, studying charge transporting behavior with a full first-principle approach demands a large amount of computation resources in quantum chemistry (QC) calculation. To address this issue, we developed a machine learning (ML) approach to evaluate electronic coupling. A prototypical ML model for an ethylene system was built by kernel ridge regression with Coulomb matrix representation. Since the performance of the ML models highly dependent on their building strategies, we systematically investigated the generality of the ML models, the choice of features and target labels. The best ML model trained with 40 000 samples achieved a mean absolute error of 3.5 meV and greater than 98% accuracy in predicting phases. The distance and orientation dependence of electronic coupling was successfully captured. Bypassing QC calculation, the ML model saved 10-10⁴ times the computation cost. With the help of ML, reliable charge transport models and mechanisms can be further developed.

Trisulfide-Bond Acenes for Organic Batteries

The molecular design of organic battery electrodes is a big challenge. Here, we synthesize two metal-free organosulfur acenes and shed insight into battery properties using first-principles calculations. A new zone-melting chemical-vapor-transport (ZM-CVT) apparatus was fabricated to provide a simple, solvent-free, and continuous synthetic protocol, and produce single crystals of tetrathiotetracene (TTT) and hexathiapentacene (HTP) at a large scale. Single crystals of HTP showed better Li-ion battery performance and higher cycling stability than those of TTT. A two-step, three-electron lithiation mechanism instead of the commonly depicted two-electron mechanism is proposed for the HTP Li-ion battery. The superior performance of HTP is linked to unique trisulfide bonding scenarios, which are also responsible for the formation of empty channels along the stacking direction. In-depth theoretical analysis suggests that organosulfur acenes are potential prototypes for organic battery materials with tunable properties, and that the tuning of sulfur bonds is critical in designing these new materials.

The interplay among molecular structures, crystal symmetries and lattice energy landscapes revealed using unsupervised machine learning: a closer look at pyrrole azaphenacenes

Using unsupervised machine learning and CSPs to help crystallographers better understand how crystallizations are affected by molecular structures.