Simulation of loss mechanisms in organic solar cells: A description of the mesoscopic Monte Carlo technique and an evaluation of the first reaction method

Tying Together Multiscale Calculations for Charge Transport in P3HT: Structural Descriptors, Morphology, and Tie-Chains

Evaluating new, promising organic molecules to make next-generation organic optoelectronic devices necessitates the evaluation of charge carrier transport performance through the semi-conducting medium. In this work, we utilize quantum chemical calculations (QCC) and kinetic Monte Carlo (KMC) simulations to predict the zero-field hole mobilities of ∼100 morphologies of the benchmark polymer poly(3-hexylthiophene), with varying simulation volume, structural order, and chain-length polydispersity. Morphologies with monodisperse chains were generated previously using an optimized molecular dynamics force-field and represent a spectrum of nanostructured order. We discover that a combined consideration of backbone clustering and system-wide disorder arising from side-chain conformations are correlated with hole mobility. Furthermore, we show that strongly interconnected thiophene backbones are required for efficient charge transport. This definitively shows the role "tie-chains" play in enabling mobile charges in P3HT. By marrying QCC and KMC over multiple length- and time-scales, we demonstrate that it is now possible to routinely probe the relationship between molecular nanostructure and device performance.

Optimization and Validation of Efficient Models for Predicting Polythiophene Self-Assembly

We develop an optimized force-field for poly(3-hexylthiophene) (P3HT) and demonstrate its utility for predicting thermodynamic self-assembly. In particular, we consider short oligomer chains, model electrostatics and solvent implicitly, and coarsely model solvent evaporation. We quantify the performance of our model to determine what the optimal system sizes are for exploring self-assembly at combinations of state variables. We perform molecular dynamics simulations to predict the self-assembly of P3HT at ∼350 combinations of temperature and solvent quality. Our structural calculations predict that the highest degrees of order are obtained with good solvents just below the melting temperature. We find our model produces the most accurate structural predictions to date, as measured by agreement with grazing incident X-ray scattering experiments.

Generalized Kinetic Monte Carlo Framework for Organic Electronics

In this paper, we present our generalized kinetic Monte Carlo (kMC) framework for the simulation of organic semiconductors and electronic devices such as solar cells (OSCs) and light-emitting diodes (OLEDs). Our model generalizes the geometrical representation of the multifaceted properties of the organic material by the use of a non-cubic, generalized Voronoi tessellation and a model that connects sites to polymer chains. Herewith, we obtain a realistic model for both amorphous and crystalline domains of small molecules and polymers. Furthermore, we generalize the excitonic processes and include triplet exciton dynamics, which allows an enhanced investigation of OSCs and OLEDs. We outline the developed methods of our generalized kMC framework and give two exemplary studies of electrical and optical properties inside an organic semiconductor.

Simulating charge transport in organic semiconductors and devices: a review

Charge transport simulation can be a valuable tool to better understand, optimise and design organic transistors (OTFTs), photovoltaics (OPVs), and light-emitting diodes (OLEDs). This review presents an overview of common charge transport and device models; namely drift-diffusion, master equation, mesoscale kinetic Monte Carlo and quantum chemical Monte Carlo, and a discussion of the relative merits of each. This is followed by a review of the application of these models as applied to charge transport in organic semiconductors and devices, highlighting in particular the insights made possible by modelling. The review concludes with an outlook for charge transport modelling in organic electronics.

Influence of permittivity and energetic disorder on the spatial charge carrier distribution and recombination in organic bulk-heterojunctions

Gaining insight into the spatial charge carrier distribution and its relation to geminate and nongeminate recombination in bulk-heterojunction organic solar cells.

Charge carrier thermalization in organic diodes

Charge carrier mobilities of organic semiconductors are often characterized using steady-state measurements of space charge limited diodes. These measurements assume that charge carriers are in a steady-state equilibrium. In reality, however, energetically hot carriers are introduces by photo-excitation and injection into highly energetic sites from the electrodes. These carriers perturb the equilibrium density of occupied states, and therefore change the overall charge transport properties. In this paper, we look into the effect of energetically hot carriers on the charge transport in organic semiconductors using steady state kinetic Monte Carlo simulations. For injected hot carriers in a typical organic semiconductor, rapid energetic relaxation occurs in the order of tens of nanoseconds, which is much faster than the typical transit time of a charge carrier throught the device. Furthermore, we investigate the impact of photo-generated carriers on the steady-state mobility. For a typical organic voltaic material, an increase in mobility of a factor of 1.1 is found. Therefore, we conclude that the impact of energetically hot carriers on normal device operation is limited.

Correlation between gate-dielectric morphology at the nanoscale and charge transport properties in organic field-effect transistors

Initial preparation, thermal treatment, and microscopic structure of gate dielectric polymer materials impact on the transport properties of organic field-effect transistors, as revealed by an integrated computational approach.

An Effective Two-Orbital Quantum Chemical Model for Organic Photovoltaic Materials

We present a coarse-grained quantum chemical model of organic photovoltaic materials, which is based on the classic idea that the main physical processes involve the electrons occupying the frontier orbitals (HOMO and LUMO) of each molecule or "site". This translates into an effective electronic Hamiltonian with two electrons and two orbitals per site. The on-site parameters (one- and two-electron integrals) can be rigorously related to the ionization energy, electron affinity, and singlet and triplet first excitation energies of that site. The intersite Hamiltonian parameters are introduced in a way that is consistent with classical electrostatics, and for the one-electron part, we use a simple approximation that could be refined using information from atomistic quantum chemical calculations. The model has been implemented within the GAMESS-US package. This allows the exploration of the physics of these materials using state-of-the art quantum chemical methods on relatively large systems (hundreds of electron-donor and electron-acceptor sites). To illustrate this point, we present ground- and excited-state calculations on dimers and two-dimensional arrays of sites using the Hartree-Fock, configuration interaction, and coupled-cluster methods. The calculations provide evidence for the possibility of low-energy, long-range electron transfer in donor-acceptor heterojunctions characterized by a moderate degree of disorder.

Percolating bulk heterostructures from neutron reflectometry and small-angle scattering data

We generate percolating fullerene-polymer bulk heterostructures that are consistent with the experimental characterization of a nanostructure, in particular neutron reflectometry and small-angle neutron scattering data from as-cast and annealed poly(3-hexylthiophene) (P3HT) and [6,6]-phenyl C61-butyric acid methyl ester systems. Transport simulations correlate changes in exciton dissociation efficiency and charge collection efficiency with morphological features including characteristic domain size, fullerene concentration profile, degree of fullerene sequestration, and degree of P3HT crystallization.

Anomalous thickness-dependence of photocurrent explained for state-of-the-art planar nano-heterojunction organic solar cells

Due to their simple geometry and design, planar heterojunction (PHJ) solar cells have advantages both as potential photovoltaics with more efficient charge extraction than their bulk heterojunction (BHJ) counterparts, and as idealized interfaces to study basic device operation. The main reason for creating BHJs was the limited exciton diffusion length in the active materials of the PHJ: if an exciton is generated at a distance greater than its diffusion length from the hetero-interface of the PHJ, it would be very unlikely to be able to contribute to the photocurrent. Based on this argument one expects a maximum in the photocurrent of PHJs for a thickness of the active layer equal to the exciton diffusion length (~10 nm). However, in two recently developed PHJs that have appeared in the literature, a maximum photocurrent is observed for 60-65 nm of poly(3-hexylthiophene) (P3HT). In this work, we explore this anomaly by combining both an optical T-matrix and a kinetic Monte Carlo simulation that tracks the exciton behavior in the PHJs. The two systems considered are a P3HT/single walled carbon nanotube (SWNT) device, and a P3HT/phenyl-C61-butyric acid methyl ester (PCBM) device. The model demonstrates how a bulk exciton sink can explain the shifted maximum in the P3HT/SWNT case, whereas in the P3HT/PCBM case the maximum is mainly determined by PCBM molecules interdiffusing in the P3HT upon annealing. Based upon the results of this model it will be possible to more intelligently design nanostructured photovoltaics and optimize them toward higher efficiencies.