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.

Electronic energy transfer processes in pi-conjugated polymers

Currently pi-conjugated polymers are considered as technologically interesting materials to be used as functional building elements for the development of the new generation of optoelectronic devices. More specifically during the last few years, poly-p-phenylene materials have attracted considerable attention for their blue photoluminescence properties. This Thesis deals with the optical properties of the most representative blue light poly-p-phenylene emitters such as poly(fluorene), oligo(fluorene), poly(indenofluorene) and ladder-type penta(phenylene) derivatives. In the present work, laser induced photoluminescence spectroscopy is used as a major tool for the study of the interdependence between the dynamics of the probed photoluminescence, the molecular structures of the prepared polymeric films and the presence of chemical defects. Complementary results obtained by two-dimensional wide-angle X-ray diffraction are reported. These findings show that the different optical properties observed are influenced by the intermolecular solid-state interactions that in turn are controlled by the pendant groups of the polymer backbone. A significant feedback is delivered regarding the positive impact of a new synthetic route for the preparation of a poly(indenofluorene) derivative on the spectral purity of the compound. The energy transfer mechanisms that operate in the studied systems are addressed by doping experiments. After the evaluation of the structure/property interdependence, a new optical excitation pathway is presented. An efficient photon low-energy up-conversion that sensitises the blue emission of poly(fluorene) is demonstrated. The observed phenomenon takes place in poly(fluorene) derivatives hosts doped with metallated octaethyl porphyrins, after quasi-CW photoexcitation of intensities in the order of kW/cm2. The up-conversion process is parameterised in terms of temperature, wavelength excitation and central metal cation in the porphyrin ring. Additionally the observation of the up-conversion is extended in a broad range of poly-p-phenylene blue light emitting hosts. The dependence of the detected up-conversion intensity on the excitation intensity and doping concentration is reported. Furthermore the dynamics of the up-conversion intensity are monitored as a function of the doping concentration. These experimental results strongly suggest the existence of triplet-triplet annihilation events into the porphyrin molecules that are subsequently followed by energy transfer to the host. After confirming the occurrence of the up-conversion in solutions, cyclic voltammetry is used in order to show that the up-conversion efficiency is partially determined from the energetic alignment between the HOMO levels of the host and the dopant.

Electronic transport in the heavy fermion superconductors UPd 2 Al 3 and UNi 2 Al 3

This work addresses the electronical properties of the superconductors UPd2Al3 and UNi2Al3 on the basis of thin film experiments. These isotructural compounds are ideal candiates to study the interplay of magnetism and superconductivity due to the differences of their magnetically ordered states, as well as the experimental evidence for a magnetic pairing mechanism in UPd2Al3. Epitaxial thin film samples of UPd2Al3 and UNi2Al3 were prepared using UHV Molecular Beam Epitaxy (MBE). For UPd2Al3, the change of the growth direction from the intrinsic (001) to epitaxial (100) was predicted and sucessfully demonstrated using LaAlO3 substrates cut in (110) direction. With optimized deposition process parameters for UPd2Al3 (100) on LaAlO3 (110) superconducting samples with critical temperatures up to Tc = 1.75K were obtained. UPd2Al3-AlOx-Ag mesa junctions with superconducting base electrode were prepared and shown to be in the tunneling regime. However, no signatures of a superconducting density of states were observed in the tunneling spectra. The resistive superconducting transition was probed for a possible dependence on the current direction. In contrast to UNi2Al3, the existence of such feature was excluded in UPd2Al3 (100) thin films. The second focus of this work is the dependence of the resisitive transition in UNi2Al3 (100) thin films on the current direction. The experimental fact that the resisitive transition occurs at slightly higher temperatures for I║a than for I║c can be explained within a model of two weakly coupled superconducting bands. Evidence is presented for the key assumption of the two-band model, namely that transport in and out of the ab-plane is generated on different, weakly coupled parts of the Fermi surface. Main indications are the angle dependence of the superconducting transition and the dependence of the upper critical field Bc2 on current and field orientation. Additionally, several possible alternative explanations for the directional splitting of the transition are excluded in this work. An origin due to scattering on crystal defects or impurities is ruled out, likewise a relation to ohmic heating or vortex dynamics. The shift of the transition temperature as function of the current density was found to behave as predicted by the Ginzburg-Landau theory for critical current depairing, which plays a significant role in the two-band model. In conclusion, the directional splitting of the resisitive transition has to be regarded an intrinsic and unique property of UNi2Al3 up to now. Therefore, UNi2Al3 is proposed as a role model for weakly coupled multiband superconductivity. Magnetoresistance in the normalconducting state was measured for UPd2Al3 and UNi2Al3. For UNi2Al3, a negative contribution was observed close to the antiferromagnetic ordering temperature TN only for I║a, which can be associated to reduced spin-disorder scattering. In agreement with previous results it is concluded that the magnetic moments have to be attributed to the same part of the Fermi surface which generates transport in the ab-plane.

Cellular localization and mechanism of amyloid precursor protein (APP) homodimer formation in an oxidizing environment

The amyloid precursor protein (APP) is a type I transmembrane glycoprotein, which resembles a cell surface receptor, comprising a large ectodomain, a single spanning transmembrane part and a short C-terminal, cytoplasmic domain. It belongs to a conserved gene family, with over 17 members, including also the two mammalian APP homologues proteins APLP1 and APLP2 („amyloid precursor like proteins“). APP is encoded by 19 exons, of which exons 7, 8, and 15 can be alternatively spliced to produce three major protein isoforms APP770, APP751 and APP695, reflecting the number of amino acids. The neuronal APP695 is the only isoform that lacks a Kunitz Protease Inhibitor (KPI) domain in its extracellular portion whereas the two larger, peripheral APP isoforms, contain the 57-amino-acid KPI insert. rnRecently, research effort has suggested that APP metabolism and function is thought to be influenced by homodimerization and that the oligomerization state of APP could also play a role in the pathology of Alzheimer's disease (AD), by regulating its processing and amyloid beta production. Several independent studies have shown that APP can form homodimers within the cell, driven by motifs present in the extracellular domain, as well as in the juxtamembrane (JM) and transmembrane (TM) regions of the molecule, whereby the exact molecular mechanism and the origin of dimer formation remains elusive. Therefore, we focused in our study on the actual subcellular origin of APP homodimerization within the cell, an underlying mechanism, and a possible impact on dimerization properties of its homologue APLP1. Furthermore, we analyzed homodimerization of various APP isoforms, in particular APP695, APP751 and APP770, which differ in the presence of a Kunitz-type protease inhibitor domain (KPI) in the extracellular region. In order to assess the cellular origin of dimerization under different cellular conditions, we established a mammalian cell culture model-system in CHO-K1 (chinese hamster ovary) cells, stably overexpressing human APP, harboring dilysine based organelle sorting motifs at the very C-terminus [KKAA-Endoplasmic Reticulum (ER); KKFF-Golgi]. In this study we show that APP exists as disulfide-bound, SDS-stable dimers, when it was retained in the ER, unlike when it progressed further to the cis-Golgi, due to the KKFF ER exit determinant. These stable APP complexes were isolated from cells, and analyzed by SDS–polyacrylamide gel electrophoresis under non-reducing conditions, whereas strong denaturing and reducing conditions completely converted those dimers to monomers. Our findings suggested that APP homodimer formation starts early in the secretory pathway and that the unique oxidizing environment of the ER likely promotes intermolecular disulfide bond formation between APP molecules. We particularly visualized APP dimerization employing a variety of biochemical experiments and investigated the origin of its generation by using a Bimolecular Fluorescence Complementation (BiFC) approach with split GFP-APP chimeras. Moreover, using N-terminal deletion constructs, we demonstrate that intermolecular disulfide linkage between cysteine residues, exclusively located in the extracellular E1 domain, represents another mechanism of how an APP sub-fraction can dimerize within the cell. Additionally, mutational studies revealed that cysteines at positions 98 and 105, embedded in the conserved loop region within the E1 domain, are critical for interchain disulfide bond formation. Using a pharmacological treatment approach, we show that once generated in the oxidative environment of the ER, APP dimers remain stably associated during transport, reaching the plasma membrane. In addition, we demonstrate that APP isoforms, encompassing the KPI domain, exhibit a strongly reduced ability to form cis-directed dimers in the ER, whereas trans-directed cell aggregation of Drosophila Schneider (S2)-cells was isoform independent, mediating cell-cell contacts. Thus, suggesting that steric properties of KPI-APP might be the cause for weaker cis-interaction in the ER, compared to APP695. Finally, we provide evidence that APP/APLP1 heterointeractions are likewise initiated in the ER, suggesting a similar mechanism for heterodimerization. Therefore, dynamic alterations of APP between monomeric, homodimeric, and possibly heterodimeric status could at least partially explain some of the variety in the physiological functions of APP.rn

Variational Monte Carlo study of the hydrogen electronic structure at ultra-high pressure

Die causa finalis der vorliegenden Arbeit ist das Verständnis des Phasendiagramms von Wasserstoff bei ultrahohen Drücken, welche von nichtleitendem H2 bis hin zu metallischem H reichen. Da die Voraussetzungen für ultrahohen Druck im Labor schwer zu schaffen sind, bilden Computersimulationen ein wichtiges alternatives Untersuchungsinstrument. Allerdings sind solche Berechnungen eine große Herausforderung. Eines der größten Probleme ist die genaue Auswertung des Born-Oppenheimer Potentials, welches sowohl für die nichtleitende als auch für die metallische Phase geeignet sein muss. Außerdem muss es die starken Korrelationen berücksichtigen, die durch die kovalenten H2 Bindungen und die eventuellen Phasenübergänge hervorgerufen werden. Auf dieses Problem haben unsere Anstrengungen abgezielt. Im Kontext von Variationellem Monte Carlo (VMC) ist die Shadow Wave Function (SWF) eine sehr vielversprechende Option. Aufgrund ihrer Flexibilität sowohl lokalisierte als auch delokalisierte Systeme zu beschreiben sowie ihrer Fähigkeit Korrelationen hoher Ordnung zu berücksichtigen, ist sie ein idealer Kandidat für unsere Zwecke. Unglücklicherweise bringt ihre Formulierung ein Vorzeichenproblem mit sich, was die Anwendbarkeit limitiert. Nichtsdestotrotz ist es möglich diese Schwierigkeit zu umgehen indem man die Knotenstruktur a priori festlegt. Durch diesen Formalismus waren wir in der Lage die Beschreibung der Elektronenstruktur von Wasserstoff signifikant zu verbessern, was eine sehr vielversprechende Perspektive bietet. Während dieser Forschung haben wir also die Natur des Vorzeichenproblems untersucht, das sich auf die SWF auswirkt, und dabei ein tieferes Verständnis seines Ursprungs erlangt. Die vorliegende Arbeit ist in vier Kapitel unterteilt. Das erste Kapitel führt VMC und die SWF mit besonderer Ausrichtung auf fermionische Systeme ein. Kapitel 2 skizziert die Literatur über das Phasendiagramm von Wasserstoff bei ultrahohem Druck. Das dritte Kapitel präsentiert die Implementierungen unseres VMC Programms und die erhaltenen Ergebnisse. Zum Abschluss fasst Kapitel 4 unsere Bestrebungen zur Lösung des zur SWF zugehörigen Vorzeichenproblems zusammen.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.