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.

Solid-state NMR on columnar supramolecular structures

The present work deals with the characterisation of three columnar self-assembled systems, that is, benzene-1,3,5-tricarboxamides, a peripherally thioalkyl-substituted phthalocyanine, and several oligo-(p-phenylenevinylene)s. In order to probe the supramolecular organisation solid-state NMR has been used as the main technique, supported by X-ray measurements, theoretical methods, and thermal analysis. rnrnBenzene-1,3,5-tricarboxamides (BTAs) turned out to be well suited model compounds to study various fundamental supramolecular interactions, such as π-π-interactions, hydrogen bonding, as well as dynamic and steric effects of attached side chains. Six BTAs have been investigated in total, five with a CO-centred amide group bearing different side chains and one with an inverted N-centred amide group. The physical properties of these BTAs have been investigated as a function of temperature. The results indicated that in case of the CO-centred BTAs the stability of the columnar mesophase depends strongly on the nature of the side chains. Further experiments revealed a coplanar orientation of adjacent BTA molecules in the columnar assembly of CO-centred BTAs, whereas the N-centred BTA, showed a deviating not fully coplanar arrangement. These differences were ascribed to distinct hydrogen bonding schemes, involving a parallel alignment of hydrogen bonds in case of CO-centred BTAs and an antiparallel alignment in case of the N-centred counterpart.rnrn The fundamental insights of the supramolecular organisation of BTAs could be partially adapted to an octa-substituted phthalocyanine with thiododecyl moieties. Solid-state NMR in combination with chemical shift calculations determined a tilted herringbone arrangement of phthalocyanine rings in the crystalline phase as well as in the mesophase. Moreover, 1H NMR measurements in the mesophase of this compound suggested an axial rotation of molecules, which is inhibited in the crystalline phase.rnrnAs a third task, the supramolecular assembly of oligo-(p-phenylenevinylene)s of varying length and with different polar head groups have been investigated by a combined X-ray and solid-state NMR study. The results revealed a columnar structure formation of these compounds, being promoted by phase separation of alkyl side chains and aromatic rigid rods. In this system solid-state NMR yielded meaningful insight into the isotropisation process of butoxy and 2-S-methylbutoxy substituted oligo-(p-phenylenevinylene) rods.rn

Polymer-nanoparticle composites: from orientation phenomena to optoelectronic applications

Polymer-nanoparticle hybrids show synergistic effects, demonstrating both, the unique properties of nanosized structures and the good processability and functionalities of polymeric materials. This work shows the synthesis and application of block copolymers containing a soluble, functional block and a short anchor block, which efficiently binds to the surface of nanocrystals. We functionalized anisotropic, semiconducting nanoparticles, which can be dissolved in organic and polymeric matrices upon modification. The modified nanorods have the ability to form liquid crystalline phases, which behave similar to low molecular liquid crystals with a reversible clearing behaviour. These liquid crystalline phases could also be obtained in hole conducting matrices. For a macroscopic orientation of the nanorods, electric fields were applied and a switching (in analogy to known liquid crystals) to a homeotropic orientation was observed.rnBy introduction of dye molecules in the anchor block of a hole conducting block copolymer, all essential components of a solar cell can be combined in a single particle. Light absorption of the dye induces the injection of electrons into the particles, followed by a charging, that was monitored by a special AFM technique.rnLight emitting nanocrystals were functionalized analogously with a hole transporting polymer. The stability of the particles could be enhanced by the sterically stabilizing polymer corona and the particles showed improved properties in terms of processing. We applied these hybrid materials in light emitting devices, which showed better characteristics due to an improved hole injection and well dispersed emitting particles in the active device layer.rnThe work shows the broad spectrum of properties and applications based on the synergistic effects in hybrid and composite materials.

Integration von Heizern in thermotrope flüssigkristalline Elastomer-Aktoren und deren Anwendung beim technischen Nachbau eines menschlichen Auges

Flüssigkristalline Elastomere (LCE) zeigen eine reversible Kontraktion und werden in der Literatur auch als „künstliche Muskeln“ bezeichnet. In dieser Arbeit werden sie mit einem integrierten Heizer versehen, um eine schnelle und präzise Ansteuerung zu ermöglichen. Anschließend werden diese als Aktoren zur Realisierung eines technischen Nachbaus des menschlichen Auges verwendet. rnDas einzigartige Verhalten der flüssigkristallinen Elastomere beruht auf der Kombination der Entropie Elastizität des Elastomers mit der Selbstorganisation der flüssigkristallinen Einheiten (Mesogene). Diese beiden Eigenschaften ermöglichen eine reversible, makroskopische Verformung beim Phasenübergang des Flüssigkristalls in die isotrope Phase. Hierbei ist es wichtig eine homogene Orientierung der Mesogene zu erzeugen, was in dieser Arbeit durch ein Magnetfeld erreicht wird. Da es sich um ein thermotropes flüssigkristallines Elastomer handelt, werden in dieser Arbeit zwei Ansätze vorgestellt, um den LCE intern zu heizen. Zum einen werden Kohlenstoffnanoröhren integriert, um diese über Strahlung oder Strom zu heizen und zum anderen wird ein flexibler Heizdraht integriert, welcher ebenfalls über Strom geheizt wird. rnUm den technischen Nachbau des menschlichen Auges zu realisieren, ist die Herstellung einer flüssigkristallinen Iris gezeigt. Hierzu wird ein radiales Magnetfeld aufgebaut, welches eine radiale Orientierung des Mesogene ermöglicht, wodurch wiederum eine radiale Kontraktion ermöglicht wird. Außerdem sind zwei Konzepte vorgestellt, um eine Elastomer Linse zu verformen. Zum einen wird diese mit einem ringförmigen LCE auseinandergezogen und somit abgeflacht. Zum anderen sind acht Aktoren über Anker an einer Linse angebracht, welche ebenfalls eine Vergrößerung der Linse bewirken. In beiden Fällen werden LCE mit dem zuvor präsentierten integrierten Heizdraht verwendet. Abschließend ist das Zusammensetzen des technische Nachbaus des menschlichen Auges dargestellt, sowie Aufnahmen, welche mit diesem erzeugt wurden.

The impact of ice crystals on radiative forcing and remote sensing of arctic boundary-layer mixed-phase clouds

This PhD thesis is embedded into the Arctic Study of Tropospheric Aerosol, Clouds and Radiation (ASTAR) and investigates the radiative transfer through Arctic boundary-layer mixed-phase (ABM) clouds. For this purpose airborne spectral solar radiation measurements and simulations of the solar and thermal infrared radiative transfer have been performed. This work reports on measurements with the Spectral Modular Airborne Radiation measurement sysTem (SMART-Albedometer) conducted in the framework of ASTAR in April 2007 close to Svalbard. For ASTAR the SMART-Albedometer was extended to measure spectral radiance. The development and calibration of the radiance measurements are described in this work. In combination with in situ measurements of cloud particle properties provided by the Laboratoire de M¶et¶eorologie Physique (LaMP) and simultaneous airborne lidar measurements by the Alfred Wegener Institute for Polar and Marine Research (AWI) ABM clouds were sampled. The SMART-Albedometer measurements were used to retrieve the cloud thermodynamic phase by three different approaches. A comparison of these results with the in situ and lidar measurements is presented in two case studies. Beside the dominating mixed-phase clouds pure ice clouds were found in cloud gaps and at the edge of a large cloud field. Furthermore the vertical distribution of ice crystals within ABM clouds was investigated. It was found that ice crystals at cloud top are necessary to describe the observed SMART-Albedometer measurements. The impact of ice crystals on the radiative forcing of ABM clouds is in vestigated by extensive radiative transfer simulations. The solar and net radiative forcing was found to depend on the ice crystal size, shape and the mixing ratio of ice crystals and liquid water droplets.

Computer simulation methods to study interfacial tensions : from the Ising model to colloidal crystals

In condensed matter systems, the interfacial tension plays a central role for a multitude of phenomena. It is the driving force for nucleation processes, determines the shape and structure of crystalline structures and is important for industrial applications. Despite its importance, the interfacial tension is hard to determine in experiments and also in computer simulations. While for liquid-vapor interfacial tensions there exist sophisticated simulation methods to compute the interfacial tension, current methods for solid-liquid interfaces produce unsatisfactory results.rnrnAs a first approach to this topic, the influence of the interfacial tension on nuclei is studied within the three-dimensional Ising model. This model is well suited because despite its simplicity, one can learn much about nucleation of crystalline nuclei. Below the so-called roughening temperature, nuclei in the Ising model are not spherical anymore but become cubic because of the anisotropy of the interfacial tension. This is similar to crystalline nuclei, which are in general not spherical but more like a convex polyhedron with flat facets on the surface. In this context, the problem of distinguishing between the two bulk phases in the vicinity of the diffuse droplet surface is addressed. A new definition is found which correctly determines the volume of a droplet in a given configuration if compared to the volume predicted by simple macroscopic assumptions.rnrnTo compute the interfacial tension of solid-liquid interfaces, a new Monte Carlo method called ensemble switch method'' is presented which allows to compute the interfacial tension of liquid-vapor interfaces as well as solid-liquid interfaces with great accuracy. In the past, the dependence of the interfacial tension on the finite size and shape of the simulation box has often been neglected although there is a nontrivial dependence on the box dimensions. As a consequence, one needs to systematically increase the box size and extrapolate to infinite volume in order to accurately predict the interfacial tension. Therefore, a thorough finite-size scaling analysis is established in this thesis. Logarithmic corrections to the finite-size scaling are motivated and identified, which are of leading order and therefore must not be neglected. The astounding feature of these logarithmic corrections is that they do not depend at all on the model under consideration. Using the ensemble switch method, the validity of a finite-size scaling ansatz containing the aforementioned logarithmic corrections is carefully tested and confirmed. Combining the finite-size scaling theory with the ensemble switch method, the interfacial tension of several model systems, ranging from the Ising model to colloidal systems, is computed with great accuracy.

Monte Carlo simulations of nucleation of colloidal crystals

A crystal nucleus in a finite volume may exhibit phase coexistence with a surrounding fluid. The thermodynamic properties of the coexisting fluid (pressure and chemical potential) are enhanced relative to their coexistence values. This enhancement is uniquely related to the surface excess free energy. rnA model for weakly attractive soft colloidal particles is investigated, the so called Asakura-Oosawa model. In simulations, this model allows for the calculation of the pressure in the liquid using the virial formula directly. The phase coexistence pressure in the thermodynamic limit is obtained from the interface velocity method. We introduce a method by which the chemical potential in dense liquids can be measured. There is neither a need to locate the interface nor to compute the anisotropic interfacial tension to obtain nucleation barriers. Therefore, our analysis is appropriate for nuclei of arbitrary shape. Monte Carlo simulations over a wide range of nucleus volumes yield to nucleation barriers independent from the total system volume. The interfacial tension is determined via the ensemble-switch method, hence a detailed test of classical nucleation theory is possible. The anisotropy of the interfacial tension and the resulting non-spherical shape has only a minor effect on the barrier for the Asakura-Oosawa model.