Modeling charge transport in high-mobility molecular semiconductors: Balancing electronic structure and quantum dynamics methods with the help of experiments.

Computing the charge mobility of molecular semiconductors requires a balanced set of approximations covering both the electronic structure of the Hamiltonian parameters and the modeling of the charge dynamics. For problems of such complexity, it is hard to make progress without independently validating each layer of approximation. In this perspective, we survey how all terms of the model Hamiltonian can be computed and validated by independent experiments and discuss whether some common approximations made to build the model Hamiltonian are valid. We then consider the range of quantum dynamics approaches used to model the charge carrier dynamics stressing the strong and weak points of each method on the basis of the available computational results. Finally, we discuss non-trivial aspects and novel opportunities related to the comparison of theoretical predictions with recent experimental data.

Impact of electron-phonon coupling on the quantum yield of photovoltaic devices.

In describing the charge carriers' separation mechanism in the organic solar cell, providing a method, which considers the impact of all parameters of interest on the same footing within an inexpensive numerical effort, could play an essential role. We use here a simple tight-binding model to describe the dissociation of the charge carriers and investigate their dependence on the physical parameters of the system. We demonstrate that the quantum yield of the cell is subtly controlled by the collective action of the Coulomb interaction of the electron-hole pair, electron-phonon coupling, and the geminate recombination of the charge carriers. This approach should help us understand the performance of organic solar cells and optimize their efficiency.

The Conductance Isotope Effect in Oligophenylene Imine Molecular Wires Depends on the Number and Spacing of 13C-Labeled Phenylene Rings.

We report a strong and structurally sensitive 13C intramolecular conductance isotope effect (CIE) for oligophenyleneimine (OPI) molecular wires connected to Au electrodes. Wires were built from Au surfaces beginning with the formation of 4-aminothiophenol self-assembled monolayers (SAMs) followed by subsequent condensation reactions with 13C-labeled terephthalaldehyde and phenylenediamine; in these monomers the phenylene rings were either completely 13C-labeled or the naturally abundant 12C isotopologues. Alternatively, perdeuterated versions of terephthalaldehyde and phenylenediamine were employed to make 2H(D)-labeled OPI wires. For 13C-isotopologues of short OPI wires (<4 nm) in length where the charge transport mechanism is tunneling, there was no measurable effect, i.e., 13C CIE ≈ 1, where CIE is defined as the ratio of labeled and unlabeled wire resistances, i.e., CIE = Rheavy/Rlight. However, for long OPI wires >4 nm, in which the transport mechanism is polaron hopping, a strong 13C CIE = 4-5 was observed. A much weaker inverse CIE < 1 was evident for the longest D-labeled wires. Importantly, the magnitude of the 13C CIE was sensitive to the number and spacing of 13C-labeled rings, i.e., the CIE was structurally sensitive. The structural sensitivity is intriguing because it may be employed to understand polaron hopping mechanisms and charge localization/delocalization in molecular wires. A preliminary theoretical analysis explored several possible explanations for the CIE, but so far a fully satisfactory explanation has not been identified. Nevertheless, the latest results unambiguously demonstrate structural sensitivity of the heavy atom CIE, offering directions for further utilization of this interesting effect.

Feasibility of p-Doped Molecular Crystals as Transparent Conductive Electrodes via Virtual Screening.

Transparent conducting materials are an essential component of optoelectronic devices. It is proven difficult, however, to develop high-performance materials that combine the often-incompatible properties of transparency and conductivity, especially for p-type-doped materials. In this work, we have employed a large set of molecular semiconductors extracted from the Cambridge Structural Database to evaluate the likelihood of transparent conducting material technology based on p-type-doped molecular crystals. Candidates are identified imposing the condition of high highest occupied molecular orbital (HOMO) energy level (for the material to be easily dopable), high charge carrier mobility (for the material to display large conductivity when doped), and a high threshold for energy absorption (for the material to absorb radiation only in the ultraviolet). The latest condition is found to be the most stringent criterion in a virtual screening protocol on a database composed of structures with sufficiently wide two-dimensional (2D) electronic bands. Calculation of excited-state energy is shown to be essential as the HOMO-lowest unoccupied molecular orbital (LUMO) gap cannot be reliably used to predict the transparency of this material class. Molecular semiconductors with desirable mobility are transparent because they display either forbidden electronic transition(s) to the lower excited states or small exchange energy between the frontier orbitals. Both features are difficult to design but can be found in a good number of compounds through virtual screening.

Organic materials repurposing, a data set for theoretical predictions of new applications for existing compounds.

We present a data set of 48182 organic semiconductors, constituted of molecules that were prepared with a documented synthetic pathway and are stable in solid state. We based our search on the Cambridge Structural Database, from which we selected semiconductors with a computational funnel procedure. For each entry we provide a set of electronic properties relevant for organic materials research, and the electronic wavefunction for further calculations and/or analyses. This data set has low bias because it was not built from a set of materials designed for organic electronics, and thus it provides an excellent starting point in the search of new applications for known materials, with a great potential for novel physical insight. The data set contains molecules used as benchmarks in many fields of organic materials research, allowing to test the reliability of computational screenings for the desired application, "rediscovering" well-known molecules. This is demonstrated by a series of different applications in the field of organic materials, confirming the potential for the repurposing of known organic molecules.

Quantitative Hole Mobility Simulation and Validation in Substituted Acenes.

Knowledge of the full phonon spectrum is essential to accurately calculate the dynamic disorder (σ) and hole mobility (µh) in organic semiconductors (OSCs). However, most vibrational spectroscopy techniques under-measure the phonons, thus limiting the phonon validation. Here, we measure and model the full phonon spectrum using multiple spectroscopic techniques and predict µh using σ from only the Γ-point and the full Brillouin zone (FBZ). We find that only inelastic neutron scattering (INS) provides validation of all phonon modes, and that σ in a set of small molecule semiconductors can be miscalculated by up to 28% when comparing Γ-point against FBZ calculations. A subsequent mode analysis shows that many modes contribute to σ and that no single mode dominates. Our results demonstrate the importance of a thoroughly validated phonon calculation, and a need to develop design rules considering the full spectrum of phonon modes.

Bright Frenkel Excitons in Molecular Crystals: A Survey.

We computed the optical properties of a large set of molecular crystals (∼2200 structures) composed of molecules whose lowest excited states are strongly coupled and generate wide excitonic bands. Such bands are classified in terms of their dimensionality (1-, 2-, and 3-dimensional), the position of the optically allowed state in relation with the excitonic density of states, and the presence of Davydov splitting. The survey confirms that one-dimensional aggregates are rare in molecular crystals highlighting the need to go beyond the simple low-dimensional models. Furthermore, this large set of data is used to search for technologically interesting and less common properties. For instance, we considered the largest excitonic bandwidth that is achievable within known molecular crystals and identified materials with strong super-radiant states. Finally, we explored the possibility that strong excitonic coupling can be used to generate emissive states in the near-infrared region in materials formed by molecules with bright visible absorption and we could identify the maximum allowable red shift in this material class. These insights with the associated searchable database provide practical guidelines for designing materials with interesting optical properties.

High-throughput virtual screening for organic electronics: a comparative study of alternative strategies.

We present a review of the field of high-throughput virtual screening for organic electronics materials focusing on the sequence of methodological choices that determine each virtual screening protocol. These choices are present in all high-throughput virtual screenings and addressing them systematically will lead to optimised workflows and improve their applicability. We consider the range of properties that can be computed and illustrate how their accuracy can be determined depending on the quality and size of the experimental datasets. The approaches to generate candidates for virtual screening are also extremely varied and their relative strengths and weaknesses are discussed. The analysis of high-throughput virtual screening is almost never limited to the identification of top candidates and often new patterns and structure-property relations are the most interesting findings of such searches. The review reveals a very dynamic field constantly adapting to match an evolving landscape of applications, methodologies and datasets.

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials.

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. Here we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.