Simulations of charge transport in organic compounds

To aid the design of organic semiconductors, we study the charge transport properties of organic liquid crystals, i.e. hexabenzocoronene and carbazole macrocycle, and single crystals, i.e. rubrene, indolocarbazole and benzothiophene derivatives (BTBT, BBBT). The aim is to find structure-property relationships linking the chemical structure as well as the morphology with the bulk charge carrier mobility of the compounds. To this end, molecular dynamics (MD) simulations are performed yielding realistic equilibrated morphologies. Partial charges and molecular orbitals are calculated based on single molecules in vacuum using quantum chemical methods. The molecular orbitals are then mapped onto the molecular positions and orientations, which allows calculation of the transfer integrals between nearest neighbors using the molecular orbital overlap method. Thus we obtain realistic transfer integral distributions and their autocorrelations. In case of organic crystals the differences between two descriptions of charge transport, namely semi-classical dynamics (SCD) in the small polaron limit and kinetic Monte Carlo (KMC) based on Marcus rates, are studied. The liquid crystals are investigated solely in the hopping limit. To simulate the charge dynamics using KMC, the centers of mass of the molecules are mapped onto lattice sites and the transfer integrals are used to compute the hopping rates. In the small polaron limit, where the electronic wave function is spread over a limited number of neighboring molecules, the Schroedinger equation is solved numerically using a semi-classical approach. The results are compared for the different compounds and methods and, where available, with experimental data. The carbazole macrocycles form columnar structures arranged on a hexagonal lattice with side chains facing inwards, so columns can closely approach each other allowing inter-columnar and thus three-dimensional transport. When taking only intra-columnar transport into account, the mobility is orders of magnitude lower than in the three-dimensional case. BTBT is a promising material for solution-processed organic field-effect transistors. We are able to show that, on the time-scales of charge transport, static disorder due to slow side chain motions is the main factor determining the mobility. The resulting broad transfer integral distributions modify the connectivity of the system but sufficiently many fast percolation paths remain for the charges. Rubrene, indolocarbazole and BBBT are examples of crystals without significant static disorder. The high mobility of rubrene is explained by two main features: first, the shifted cofacial alignment of its molecules, and second, the high center of mass vibrational frequency. In comparsion to SCD, only KMC based on Marcus rates is capable of describing neighbors with low coupling and of taking static disorder into account three-dimensionally. Thus it is the method of choice for crystalline systems dominated by static disorder. However, it is inappropriate for the case of strong coupling and underestimates the mobility of well-ordered crystals. SCD, despite its one-dimensionality, is valuable for crystals with strong coupling and little disorder. It also allows correct treatment of dynamical effects, such as intermolecular vibrations of the molecules. Rate equations are incapable of this, because simulations are performed on static snapshots. We have thus shown strengths and weaknesses of two state of the art models used to study charge transport in organic compounds, partially developed a program to compute and visualize transfer integral distributions and other charge transport properties, and found structure-mobility relations for several promising organic semiconductors.

Morphology and charge transport in conjugated polymers

To assist rational compound design of organic semiconductors, two problems need to be addressed. First, the material morphology has to be known at an atomistic level. Second, with the morphology at hand, an appropriate charge transport model needs to be developed in order to link charge carrier mobility to structure.rnrnThe former can be addressed by generating atomistic morphologies using molecular dynamics simulations. However, the accessible range of time- and length-scales is limited. To overcome these limitations, systematic coarse-graining methods can be used. In the first part of the thesis, the Versatile Object-oriented Toolkit for Coarse-graining Applications is introduced, which provides a platform for the implementation of coarse-graining methods. Tools to perform Boltzmann inversion, iterative Boltzmann inversion, inverse Monte Carlo, and force-matching are available and have been tested on a set of model systems (water, methanol, propane and a single hexane chain). Advantages and problems of each specific method are discussed.rnrnIn partially disordered systems, the second issue is closely connected to constructing appropriate diabatic states between which charge transfer occurs. In the second part of the thesis, the description initially used for small conjugated molecules is extended to conjugated polymers. Here, charge transport is modeled by introducing conjugated segments on which charge carriers are localized. Inter-chain transport is then treated within a high temperature non-adiabatic Marcus theory while an adiabatic rate expression is used for intra-chain transport. The charge dynamics is simulated using the kinetic Monte Carlo method.rnrnThe entire framework is finally employed to establish a relation between the morphology and the charge mobility of the neutral and doped states of polypyrrole, a conjugated polymer. It is shown that for short oligomers, charge carrier mobility is insensitive to the orientational molecular ordering and is determined by the threshold transfer integral which connects percolating clusters of molecules that form interconnected networks. The value of this transfer integral can be related to the radial distribution function. Hence, charge mobility is mainly determined by the local molecular packing and is independent of the global morphology, at least in such a non-crystalline state of a polymer.

Charge transport properties of organic conjugated polymers for photovoltaic applications

Charge transport in conjugated polymers as well as in bulk-heterojunction (BHJ) solar cells made of blends between conjugated polymers, as electron-donors (D), and fullerenes, as electron-acceptors (A), has been investigated. It is shown how charge carrier mobility of a series of anthracene-containing poly(p-phenylene-ethynylene)-alt-poly(p-phenylene-vinylene)s (AnE-PVs) is highly dependent on the lateral chain of the polymers, on a moderate variation of the macromolecular parameters (molecular weight and polydispersity), and on the processing conditions of the films. For the first time, the good ambipolar transport properties of this relevant class of conjugated polymers have been demonstrated, consistent with the high delocalization of both the frontier molecular orbitals. Charge transport is one of the key parameters in the operation of BHJ solar cells and depends both on charge carrier mobility in pristine materials and on the nanoscale morphology of the D/A blend, as proved by the results here reported. A straight correlation between hole mobility in pristine AnE-PVs and the fill factor of the related solar cells has been found. The great impact of charge transport for the performance of BHJ solar cells is clearly demonstrated by the results obtained on BHJ solar cells made of neat-C70, instead of the common soluble fullerene derivatives (PCBM or PC70BM). The investigation of neat-C70 solar cells was motivated by the extremely low cost of non-functionalized fullerenes, compared with that of their soluble derivatives (about one-tenth). For these cells, an improper morphology of the blend leads to a deterioration of charge carrier mobility, which, in turn, increases charge carrier recombination. Thanks to the appropriate choice of the donor component, solar cells made of neat-C70 exhibiting an efficiency of 4.22% have been realized, with an efficiency loss of just 12% with respect to the counterpart made with costly PC70BM.

Charge transport in amorphous organic semiconductors

Organic semiconductors with the unique combination of electronic and mechanical properties may offer cost-effective ways of realizing many electronic applications, e.g. large-area flexible displays, printed integrated circuits and plastic solar cells. In order to facilitate the rational compound design of organic semiconductors, it is essential to understand relevant physical properties e.g. charge transport. This, however, is not straightforward, since physical models operating on different time and length scales need to be combined. First, the material morphology has to be known at an atomistic scale. For this atomistic molecular dynamics simulations can be employed, provided that an atomistic force field is available. Otherwise it has to be developed based on the existing force fields and first principle calculations. However, atomistic simulations are typically limited to the nanometer length- and nanosecond time-scales. To overcome these limitations, systematic coarse-graining techniques can be used. In the first part of this thesis, it is demonstrated how a force field can be parameterized for a typical organic molecule. Then different coarse-graining approaches are introduced together with the analysis of their advantages and problems. When atomistic morphology is available, charge transport can be studied by combining the high-temperature Marcus theory with kinetic Monte Carlo simulations. The approach is applied to the hole transport in amorphous films of tris(8-hydroxyquinoline)aluminium (Alq3). First the influence of the force field parameters and the corresponding morphological changes on charge transport is studied. It is shown that the energetic disorder plays an important role for amorphous Alq3, defining charge carrier dynamics. Its spatial correlations govern the Poole-Frenkel behavior of the charge carrier mobility. It is found that hole transport is dispersive for system sizes accessible to simulations, meaning that calculated mobilities depend strongly on the system size. A method for extrapolating calculated mobilities to the infinite system size is proposed, allowing direct comparison of simulation results and time-of-flight experiments. The extracted value of the nondispersive hole mobility and its electric field dependence for amorphous Alq3 agree well with the experimental results.

Charge generation, transport and recombination in organic solar cells

This thesis presents a study of the charge generation, transport, and recombination processes in organic solar cells performed with time-resolved experimental techniques. Organic solar cells based on polymers can be solution-processed on large areas and thus promise to become an inexpensive source of renewable energy. Despite significant improvements of the power conversion efficiency over the last decade, the fundamental working principles of organic solar cells are still not fully understood. It is the aim of this thesis to clarify the role of different performance limiting processes in organic solar cells and to correlate them with the molecular structure of the studied materials, i.e. poly(3-hexylthiophene) (P3HT) and [6,6]-phenyl C61 butyric acid methyl ester (PCBM). By combining time-of-flight charge transport measurements, transient absorption spectroscopy, a newly developed experimental technique called time delayed double pulse experiment and drift-diffusion simulations a comprehensive analysis of the working principles of P3HT:PCBM solar cells could be performed. It was found that the molecular structure of P3HT (i.e. the regioregularity) has a pronounced influence on the morphology of thin films of pristine P3HT and of blends of P3HT with PCBM. This morphology in turn affected the charge transport properties as well as the charge generation and recombination kinetics. Well-ordered regioregular P3HT was found to be characterized by a high charge carrier mobility, efficient charge generation and low but field-dependent (non-geminate) recombination. Importantly, the charge generation yield was found to be independent of temperature and applied electric field as opposed to the expectations of the Onsager-Braun model that is commonly applied to describe the temperature and field dependence of charge generation in organic solar cells. These properties resulted in a reasonably good power conversion efficiency. In contrast to this, amorphous regiorandom P3HT was found to show poor charge generation, transport and recombination properties that combine to a much lower power conversion efficiency.

Thermoelectric performance of p-type TiCoSb half-Heusler compounds - Intrinsic phase separation and charge carrier concentration optimization as key to high efficiency

The world's rising demand of energy turns the development of sustainable and more efficient technologies for energy production and storage into an inevitable task. Thermoelectric generators, composed of pairs of n-type and p-type semiconducting materials, di¬rectly transform waste heat into useful electricity. The efficiency of a thermoelectric mate¬rial depends on its electronic and lattice properties, summarized in its figure of merit ZT. Desirable are high electrical conductivity and Seebeck coefficients, and low thermal con¬ductivity. Half-Heusler materials are very promising candidates for thermoelectric applications in the medium¬ temperature range such as in industrial and automotive waste heat recovery. The advantage of Heusler compounds are excellent electronic properties and high thermal and mechanical stability, as well as their low toxicity and elemental abundance. Thus, the main obstacle to further enhance their thermoelectric performance is their relatively high thermal conductivity.rn rnIn this work, the thermoelectric properties of the p-type material (Ti/Zr/Hf)CoSb1-xSnx were optimized in a multistep process. The concept of an intrinsic phase separation has recently become a focus of research in the compatible n-type (Ti/Zr/Hf)NiSn system to achieve low thermal conductivities and boost the TE performance. This concept is successfully transferred to the TiCoSb system. The phase separation approach can form a significant alternative to the previous nanostructuring approach via ball milling and hot pressing, saving pro¬cessing time, energy consumption and increasing the thermoelectric efficiency. A fundamental concept to tune the performance of thermoelectric materials is charge carrier concentration optimization. The optimum carrier concentration is reached with a substitution level for Sn of x = 0.15, enhancing the ZT about 40% compared to previous state-of-the-art samples with x = 0.2. The TE performance can be enhanced further by a fine-tuning of the Ti-to-Hf ratio. A correlation of the microstructure and the thermoelectric properties is observed and a record figure of merit ZT = 1.2 at 710°C was reached with the composition Ti0.25Hf0.75CoSb0.85Sn0.15.rnTowards application, the long term stability of the material under actual conditions of operation are an important issue. The impact of such a heat treatment on the structural and thermoelectric properties is investigated. Particularly, the best and most reliable performance is achieved in Ti0.5Hf0.5CoSb0.85Sn0.15, which reached a maximum ZT of 1.1 at 700°C. The intrinsic phase separation and resulting microstructure is stable even after 500 heating and cooling cycles.

Giant enhancement of n-type carrier mobility in highly strained germanium nanostructures

First-principles electronic structure methods are used to predict the rate of n-type carrier scattering due to phonons in highly-strained Ge. We show that strains achievable in nanoscale structures, where Ge becomes a direct bandgap semiconductor, cause the phonon-limited mobility to be enhanced by hundreds of times that of unstrained Ge, and over a thousand times that of Si. This makes highly tensile strained Ge a most promising material for the construction of channels in CMOS devices, as well as for Si-based photonic applications. Biaxial (001) strain achieves mobility enhancements of 100 to 1000 with strains over 2%. Low temperature mobility can be increased by even larger factors. Second order terms in the deformation potential of the Gamma valley are found to be important in this mobility enhancement. Although they are modified by shifts in the conduction band valleys, which are caused by carrier quantum confinement, these mobility enhancements persist in strained nanostructures down to sizes of 20 nm.

Lead oxide-modified TiO2 photocatalyst: tuning light absorption and charge carrier separation by lead oxidation state

Modification of TiO2 with metal oxide nanoclusters such as FeOx, NiOx has been shown to be a promising approach to the design of new photocatalysts with visible light absorption and improved electron–hole separation. To study further the factors that determine the photocatalytic properties of structures of this type, we present in this paper a first principles density functional theory (DFT) investigation of TiO2 rutile(110) and anatase(001) modified with PbO and PbO2 nanoclusters, with Pb2+ and Pb4+ oxidation states. This allows us to unravel the effect of the Pb oxidation state on the photocatalytic properties of PbOx-modified TiO2. The nanoclusters adsorb strongly at all TiO2 surfaces, creating new Pb–O and Ti–O interfacial bonds. Modification with PbO and PbO2 nanoclusters introduces new states in the original band gap of rutile and anatase. However the oxidation state of Pb has a dramatic impact on the nature of the modifications of the band edges of TiO2 and on the electron–hole separation mechanism. PbO nanocluster modification leads to an upwards shift of the valence band which reduces the band gap and upon photoexcitation results in hole localisation on the PbO nanocluster and electron localisation on the surface. By contrast, for PbO2 nanocluster modification the hole will be localised on the TiO2 surface and the electron on the nanocluster, thus giving rise to two different band gap reduction and electron–hole separation mechanisms. We find no crystal structure sensitivity, with both rutile and anatase surfaces showing similar properties upon modification with PbOx. In summary the photocatalytic properties of heterostructures of TiO2 with oxide nanoclusters can be tuned by oxidation state of the modifying metal oxide, with the possibility of a reduced band gap causing visible light activation and a reduction in charge carrier recombination.

Logic-gate devices based on printed polymer semiconducting nanostripes

The applications of organic semiconductors in complex circuitry such as printed CMOS-like logic circuits demand miniaturization of the active structures to the submicrometric and nanoscale level while enhancing or at least preserving the charge transport properties upon processing. Here, we addressed this issue by using a wet lithographic technique, which exploits and enhances the molecular order in polymers by spatial confinement, to fabricate ambipolar organic field effect transistors and inverter circuits based on nanostructured single component ambipolar polymeric semiconductor. In our devices, the current flows through a precisely defined array of nanostripes made of a highly ordered diketopyrrolopyrrole-benzothiadiazole copolymer with high charge carrier mobility (1.45 cm2 V-1 s-1 for electrons and 0.70 cm2 V-1 s-1 for holes). Finally, we demonstrated the functionality of the ambipolar nanostripe transistors by assembling them into an inverter circuit that exhibits a gain (105) comparable to inverters based on single crystal semiconductors.

A study of the effects metal residues in poly(9,9-dioctylfluorene) have on field-effect transistor device characteristics

The amount of metal residues from organometallic reagents used in preparation of poly(9,9-dioctylfluorene) by palladium catalysed Suzuki and nickel-induced Yamamoto polycondensations have been determined, and their effect upon the behaviour of the polymer in field-effect transistors (FETs) has been measured. The metal levels from material polymerised by Suzuki method were found to be much higher than from that made by the Yamamoto procedure. Simple treatment of the polymers with suitable metal trapping reagents lowered the metal levels significantly, with EDTA giving best results for nickel and triphenylphosphine for palladium. Comparison of the behaviour of FETs using polyfluorenes with varying levels of metal contamination, showed that the metal residues have little effect upon the mobility values, but often affect the degree of hysteresis, possibly acting as charge traps. Satisfactory device performances were obtained from polymer with palladium levels of 2000 μg/g suggesting that complete removal of metal residues may not be necessary for satisfactory device performance.

Measurement of contact voltage drop and resistance in organic solar cells

Bulk heterojunction organic solar cells based on poly[4,7-bis(3- dodecylthiophene-2-yl) benzothiadiazole-co-benzothiadiazole] and [6,6]-phenyl C71-butyric acid methyl ester are investigated. A prominent kink is observed in the fourth quadrant of the current density-voltage (J-V) response. Annealing the active layer prior to cathode deposition eliminates the kink. The kink is attributed to an extraction barrier. The J-V response in these devices is well described by a power law. This behavior is attributed to an imbalance in charge carrier mobility. An expected photocurrent for the device displaying a kink in the J-V response is determined by fitting to a power law. The difference between the expected and measured photocurrent allows for the determination of a voltage drop within the device. Under simulated 1 sun irradiance, the peak voltage drop and contact resistance at short circuit are 0.14 V and 90 Ω, respectively. © 2012 American Institute of Physics.

Enhancing current density using vertically oriented organic photovoltaics

We report a new organic photovoltaics (OPV) design, a wrapped OPV, which can circumvent both challenges of short exciton diffusion length [1], and low charge carrier mobility [2] of organic semiconductors by orienting the OPV vertically, to capture; manage; guide and use all incident photons and therefore, generate higher current. Resonant light, on being transmitted into a wrapped OPV, makes multiple passes through the photoactive layer and is absorbed completely, thus achieving benefits of thick photoactive layer while maintaining its ultra-thin thickness requirement. The current density generated from a wrapped OPV is twice than that generated by a similar OPV with flat orientation.

Unusual effects of anisotropic bonding in Cu(II) and Ni(II) oxides with K2NiF4 structure

Geometric constraints present in A2BO4 compounds with the tetragonal-T structure of K2NiF4 impose a strong pressure on the B---OII---B bonds and a stretching of the A---OI---A bonds in the basal planes if the tolerance factor is t congruent with RAO/√2 RBO < 1, where RAO and RBO are the sums of the A---O and B---O ionic radii. The tetragonal-T phase of La2NiO4 becomes monoclinic for Pr2NiO4, orthorhombic for La2CuO4, and tetragonal-T′ for Pr2CuO4. The atomic displacements in these distorted phases are discussed and rationalized in terms of the chemistry of the various compounds. The strong pressure on the B---OII---B bonds produces itinerant σ*x2−y2 bands and a relative stabilization of localized dz2 orbitals. Magnetic susceptibility and transport data reveal an intersection of the Fermi energy with the d2z2 levels for half the copper ions in La2CuO4; this intersection is responsible for an intrinsic localized moment associated with a configuration fluctuation; below 200 K the localized moment smoothly vanishes with decreasing temperature as the d2z2 level becomes filled. In La2NiO4, the localized moments for half-filled dz2 orbitals induce strong correlations among the σ*x2−y2 electrons above Td reverse similar, equals 200 K; at lower temperatures the σ*x2−y2 electrons appear to contribute nothing to the magnetic susceptibility, which obeys a Curie-Weiss law giving a μeff corresponding to S = 1/2, but shows no magnetic order to lowest temperatures. These surprising results are verified by comparison with the mixed systems La2Ni1−xCuxO4 and La2−2xSr2xNi1−xTixO4. The onset of a charge-density wave below 200 K is proposed for both La2CuO4 and La2NiO4, but the atomic displacements would be short-range cooperative in mixed systems. The semiconductor-metallic transitions observed in several systems are found in many cases to obey the relation Ea reverse similar, equals kTmin, where varrho = varrho0exp(−Ea/kT) and Tmin is the temperature of minimum resistivity varrho. This relation is interpreted in terms of a diffusive charge-carrier mobility with Ea reverse similar, equals ΔHm reverse similar, equals kT at T = Tmin.

Investigation of phase separation in bulk heterojunction solar cells via supramolecular chemistry

In this work, we have prepared two donor-acceptor-donor (D-A-D) pi-conjugated oligomers to investigate the effect of phase separation on the performance of bulk heterojunction (BHJ) solar cells. These charge transfer low band gap pi-conjugated oligomers (TTB and NMeTTB) were synthesized by Knoevenagel condensation of terthiophenecarbaldehyde and barbiturate appended pyran derivative. The thin film morphology of both the oligomers and along with electron acceptor 6,6]-phenyl-C60-butyric acid methyl ester (PC61BM) was investigated by atomic force microscopy (AFM) and transmission electron microscopy (TEM). The blend of NMeTTB and PC61BM thin film yield highly ordered thin film, whereas there was clear phase separation between TTB and PC61BM in thin film. The BHJ solar cell was fabricated using a blend of NMeTTB and TTB with PC61BM acceptor in 1:1 ratio as active layer, and a power conversion efficiency of 1.8% was obtained. This device characteristic was compared with device having TTB:PC61BM as active layer, and large difference is observed in photocurrents. This poor performance of TTB in BHJ devices was attributed to the difference in the nanoscale morphology of the corresponding derivatives. We rationalize our findings based on the low charge carrier mobility in organic field-effect transistors and miscibility/phase separation parameter of binary components (oligomers and PC61BM) in the active layer of bulk heterojunction solar cells.

Controlling Conformations of Diketopyrrolopyrrole-Based Conjugated Polymers: Role of Torsional Angle

Transport of charge carriers through conjugated polymers is strongly influenced by the presence and distribution of structural disorders. In the present work, structural defects caused by the presence of torsional.. angle were investigated in a diketopyrrolopyrrole (DPP)-based conjugated polymer. Two new copolymers of DPP were synthesized with varying torsional angles to trace the role of structural disorder. The optical properties of these copolymers in solution and thin film reveal the strong influence of torsional angle on their photophysical properties. A strong influence was observed on carrier transport properties of polymers in organic field-effect transistors (OFET) device geometry. The polymers based on phenyl DPP with higher torsional angle (PPTDPP-OD-TEG) resulted in high threshold voltage with less charge carrier mobility as compared to the polymer based on thiophene DPP (2DPP-OD-TEG) bearing a lower torsional angle. Carrier mobility and the molecular orientation of the conjugated polymers were correlated on the basis of grazing incidence X-ray scattering measurements showing the strong role of torsional angle introduced in the form of structural disorder. The results presented in this Article provide a deep insight into the sensitivity of structural disorder and its impact on the device performance of DPP-based conjugated polymers.