The orientational dependence of charge carrier mobilities in organic semiconductor crystals and the correlation with the crystal structure are investigated by means of quantum chemical first principles calculations combined with a model using hopping rates from Marcus theory. A master equation approach is presented which is numerically more efficient than the Monte Carlo method frequently applied in this context. Furthermore, it is shown that the widely used approach to calculate the mobility via the diffusion constant along with rate equations is not appropriate in many important cases. The calculations are compared with experimental data, showing good qualitative agreement for pentacene and rubrene. In addition, charge transport properties of core-fluorinated perylene bisimides are investigated.

Future prospects of the organic light emitting diode (OLED) technology rely on the development of new organic semiconductors with optical and electronic properties outperforming those of presently available materials. Computational materials design is becoming a widely used tool to complement and accelerate experimental efforts. Computational tools were also shown to contribute to the understanding of experimentally observed phenomena. Impurities and charge traps are omnipresent in most currently available organic semiconductors and limit the charge transport and thus the efficiency of the devices. The microscopic cause as well as the chemical nature of these traps is presently not well understood. Using a multiscale model we characterize the influence of impurities on the density of states and charge transport in small-molecule amorphous organic semiconductors. We use the model to quantitatively describe the influence of water molecules and water-oxygen complexes on the electron and hole mobility by influencing the shape of the density of states and at the same time acting as explicit charge traps within the energy gap. Our results show that deep trap states introduced by molecular oxygen mainly determine the electron mobility in widely used materials such as α-NPD. TOC

Exciton diffusion is a critical step for energy conversion in optoelectronic devices. This spawns the desire for theoretical approaches that allow for fast but reliable determinations of the material-dependent exciton transport parameters. For this purpose, the Marcus theory, which is widely used in the context of charge transport, is adapted to exciton diffusion. In contrast to the common approach of calculating the exciton hopping rate via the coupling and the spectral overlap, this alternative approach is less costly, because, instead of the spectral overlap, only the reorganization energy is needed. To demonstrate the capability of the approach, the diffusion constants for naphthalene, anthracene, and diindenoperylene crystals are calculated and compared with both calculations conducted with the well-established exciton hopping rate, including coupling and spectral overlap, and with experimental data. These test calculations show that Marcus-based exciton diffusion properties tend to be too small but are qualitatively correct (i.e., they seem to be useful to predict trends). Nevertheless, for reliable results, high-level quantum chemical approaches are necessary for the computation of the reorganization energies. However, they have to be calculated only once. Coupling constants, which are needed for all pairs of monomers, have a considerably smaller influence, i.e., they can be computed by a lower level approach, which makes the method even less costly.

The exciton diffusion length (LD) is a key parameter for the efficiency of organic optoelectronic devices. Its limitation to the nm length scale causes the need of complex bulk-heterojunction solar cells incorporating difficulties in long-term stability and reproducibility. A comprehensive model providing an atomistic understanding of processes that limit exciton trasport is therefore highly desirable and will be proposed here for perylene-based materials. Our model is based on simulations with a hybrid approach which combines high-level ab initio computations for the part of the system directly involved in the described processes with a force field to include environmental effects. The adequacy of the model is shown by detailed comparison with available experimental results. The model indicates that the short exciton diffusion lengths of α-perylene tetracarboxylicdianhydride (PTCDA) are due to ultrafast relaxation processes of the optical excitation via intermolecular motions leading to a state from which further exciton diffusion is hampered. As the efficiency of this mechanism depends strongly on molecular arrangement and environment, the model explains the strong dependence of LD on the morphology of the materials, for example, the differences between α-PTCDA and diindenoperylene. Our findings indicate how relaxation processes can be diminished in perylene-based materials. This model can be generalized to other organic compounds.

Due to its importance for the function of organic optoelectronic devices, accurate simulations of the singlet exciton diffusion are crucial to predict the performance of new materials. We present a protocol which allows for the efficient directional analysis of exciton transport with high-level ab initio methods. It is based on an alternative to the frequently employed rate equation since the latter was found to be erroneous in some cases. The new approach can be used in combination with the master equation which is considerably faster than the corresponding Monte Carlo approach. The long-range character of the singlet exciton coupling is taken into account by an extrapolation scheme. The approach is applied to singlet exciton diffusion in those substances where these quantities are experimentally best established: naphthalene and anthracene. The high quality of the crystals, furthermore, diminish uncertainties arising from the geometrical structures used in the computations. For those systems, our new approach provides exciton diffusion lengths L for naphthalene and anthracene crystals which show an excellent agreement with their experimental counterparts. For anthracene, for example, the computed L value in a direction is computed to 58 nm while the experimental value is 60 ± 10 nm.