Quantenkorrekturen an magnetischen Domänenwänden auf Gitterstrukturen

Wir betrachten die eindimensionale Heisenberg-Spinkette aus einem neuen und aktuelleren Blickwinkel. Experimentelle Techniken der Herstellung und selbstverständlich auch experimentelle Meßmethoden erlauben nicht nur die Herstellung von Nanopartikeln und Nanodrähten, sondern gestatten es auch, Domänenwände in diesen Strukturen auszumessen. Die meisten heute verwendeten Theorien und Simulationsmethoden haben ihre Grundlage im mikromagnetischen Kontinuumsmodell, daß schon über Jahrzehnte hinweg erforscht und erprobt ist. Wir stellen uns jedoch die Frage, ob die innere diskrete Struktur der Substrate und die quantenmechanischen Effekte bei der Genauigkeit heutiger Messungen in Betracht gezogen werden müssen. Dazu wählen wir einen anderen Ansatz. Wir werden zunächst den wohlbekannten klassischen Fall erweitern, indem wir die diskrete Struktur der Materie in unseren Berechnungen berücksichtigen. Man findet in diesem Formalismus einen strukturellen Phasenübergang zwischen einer Ising-artigen und einer ausgedehnten Wand. Das führt zu bestimmten Korrekturen im Vergleich zum Kontinuumsfall. Der Hauptteil dieser Arbeit wird sich dann mit dem quantenmechanischen Fall beschäftigen. Wir rotieren das System zunächst mit einer Reihe lokaler Transformationen derart, daß alle Spins in die z-Richtung ausgerichtet sind. Im Rahmen einer 1/S-Entwicklung läßt sich der erhaltene neue Hamilton-Operator diagonalisieren. Setzt man hier die klassische Lösung ein, so erhält man Anregungsmoden in diesem Grenzfall. Unsere Resultate erweitern und bestätigen frühere Berechnungen. Mit Hilfe der Numerik wird schließlich der Erwartungswert der Energie minimiert und somit die Form der Domänenwand im quantenmechanischen Fall berechnet. Hieraus ergeben sich auch bestimmte Korrekturen zum kritischen Verhalten des Systems. Diese Ergebnisse sind vollkommen neu.

Contribution of tectonic processes to the exhumation of the Cycladic blueschist unit, Greece and Turkey

In dieser Studie werden strukturgeologische, metamorphe und geochronologische Daten benutzt, um eine Quantifizierung tektonischer Prozesse vorzunehmen, die für die Exhumierung der Kykladischen Blauschiefereinheit in der Ägäis und der Westtürkei verantwortlich waren. Bei den beiden tektonischen Prozessen handelt es sich um: (1) Abschiebungstektonik und (2) vertikale duktile Ausdünnung. Eine finite Verformungsanalyse an Proben der Kykladischen Blauschiefereinheit ermöglicht eine Abschätzung des Beitrags von vertikaler duktiler Ausdünnung an der gesamten Exhumierung. Kalkulationen mit einem eindimensionalen, numerischen Model zeigt, daß vertikale duktile Ausdünnung nur ca. 10% an der gesamten Exhumierung ausmacht. Kinematische, metamorphe und geochronologische Daten erklären die tektonische Natur und die Evolution eines extensionalen Störungssystems auf der Insel Ikaria in der östlichen Ägäis. Thermobarometrische Daten lassen erkennen, daß das Liegende des Störungssystems aus ca. 15 km Tiefe exhumiert wurde. Sowohl Apatit- und Zirkonspaltspurenalter als auch Apatit (U-Th)/He-Alter zeigen, daß sich das extensionale Störungssystem zwischen 11-3 Ma mit einer Geschwindigkeit von ca. 7-8 km/Ma bewegte. Spät-Miozäne Abschiebungen trugen zur Exhumierung der letzten ~5-15 km der Hochdruckgesteine bei. Ein Großteil der Exhumierung der Kykladischen Blauschiefereinheit muß vor dem Miozän stattgefunden haben. Dies wird durch einen Extrusionskeil erklärt, der ca. 30-35 km der Kykladischen Blauschiefereinheit in der Westtürkei exhumierte. 40Ar/39Ar und 87Rb/86Sr Datierungen an Myloniten des oberen Abschiebungskontakts zwischen der Selçuk Decke und der darunterliegenden Ampelos/Dilek Decke der Kykladischen Blauschiefereinheit als auch des unteren Überschiebungskontakts zwischen der Ampelos/Dilek Decke und den darunterliegenden Menderes Decken zeigt, daß sich beide mylonitische Zonen um ca. ~35 Ma formten, was die Existenz eines Spät-Eozänen/Früh-Oligozänen Extrusionskeils beweist.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Untersuchung magnetischer Nanostrukturen mittels Rastertunnelspektroskopie und Kerr-Magnetometrie am Beispiel von Fe, Co, Co-Fe und Fe-Mn Nanostrukturen

Nach einer kurzen Einführung in die Entwicklung der magnetischen Anwendungen, werden in Kapitel 2 und 3 die physikalischen Grundlagen der Messmethoden, insbesondere die Rastertunnelspektroskopie und Kerr-Magnetometrie, sowie der experimentelle Aufbau erläutert. Kapitel 4 beschäftigt sich mit den magnetischen Eigenschaften von quasi ein-dimensionalen ferromagnetischen Nanostreifen und Monolagen, die durch Selbstorganisation auf einem Wolfram(110)-Einkristall mit vizinaler und glatter Oberfläche präpariert werden. Hierbei wird die Temperaturabhängigkeit der magnetischen Größen, wie Remanenz, Sättigungsmagnetisierung und Suszeptibilität, sowie die Auswirkung einer Abdeckung des Systems auf die Domänenwandenergie und Anisotropie untersucht. Zusätzlich wird die Kopplung von parallelen Nanostreifen in Abhängigkeit des Streifenabstandes betrachtet. In Kapitel 5 werden das Wachstum und die Morphologie von Co-Monolagen auf W(110) untersucht. Der Übergang von pseudomorphem zu dicht gepacktem Wachstum in der Monolage wird mit Hilfe der Rastertunnelspektroskopie sichtbar gemacht, ebenso wie unterschiedliche Stapelfolgen in Tripellagen Co-Systemen. Atomar aufgelöste Rastertunnelmikroskopie erlaubt die genauen Atompositionen der Oberfläche zu bestimmen und mit theoretischen Wachstumsmodellen zu vergleichen. Auf die Untersuchung zwei-dimensionaler binärer Co-Fe und Fe-Mn Legierungen auf W(110) wird in Kapitel 6 eingegangen. Mit einer Präparationstemperatur von T=520 K ist es möglich, atomar geordnete Co-Fe Legierungsmonolagen wachsen zu lassen. Ein direkter Zusammenhang zwischen der Magnetisierung und der lokalen Zustandsdichte in Abhängigkeit der Legierungszusammensetzung wird gezeigt.

New discotic liquid crystals based on large polycyclic aromatic hydrocarbons as materials for molecular electronics

A series of new columnar discotic liquid crystalline materials based on the superphenalene (C96) core has been synthesized by oxidative cyclodehydrogenation with iron(III) chloride of suitable three-dimensional oligophenylene precursors. These compounds were investigated by means of differential scanning calorimetry (DSC), polarized optical microscopy (POM) and wide angle X-ray scattering (WAXS), and showed highly ordered supramolecular arrays and mesophase behavior over a broad temperature range. Good solubility, through the introduction of long alkyl chains, and the fact that these new superphenalene derivatives were found to be liquid crystalline at room temperature enabled the formation of highly ordered films (using the zone-casting technique), a requirement for application in organic electronic devices. The one-dimensional, intracolumnar charge carrier mobilities of superphenalene derivatives were determined using the pulse-radiolysis time-resolved microwave conductivity technique (PR-TRMC). Electrical properties of different C96-C12 architectures on mica surfaces were examined by using Electrostatic Force Microscopy (EFM) and Kelvin Probe Force Microscopy (KPFM). Hexa-peri-hexabenzocoronene (C42) derivatives substituted at the periphery with six branched alkyl ether chains were also synthesized. It was found that the introduction of ether groups within the side chains enhances the affinity of the discotic molecules towards polar surfaces, resulting in homeotropic self-assembly (as shown by POM and 2D-WAXS) when the compounds are processed from the isotropic state between two surfaces. A new, insoluble, superphenalene building block bearing six reactive sites was prepared, and was further used for the preparation of dendronized superphenalenes with bulky dendritic substituents around the core. UV/Vis and fluorescence experiments suggest reduced π-π stacking of the superphenalene cores as a result of steric hindrance between the peripheral dendritic units. A new family of graphitic molecules with partial ”zig-zag” periphery has been established. The incorporation of ”zig-zag” edges was shown to have a strong influence on the electronic properties of the new molecules (as studied by solution and solid-state UV/Vis, and fluorescence spectroscopy), leading to a significant bathochromic shift with respect to the parent PAHs (C42 and C96). The reactivity of the additional double bonds was examined. The attachment of long alkyl chains to a ”zig-zag” superphenalene core afforded a new, processable, liquid crystalline material.

Detailed study of the seismotectonic setting in the Lucanian Apennines and surrounding areas (Southern Italy)

In this research work I analyzed the instrumental seismicity of Southern Italy in the area including the Lucanian Apennines and Bradano foredeep, making use of the most recent seismological database available so far. I examined the seismicity occurred during the period between 2001 and 2006, considering 514 events with magnitudes M ≥ 2.0. In the first part of the work, P- and S-wave arrival times, recorded by the Italian National Seismic Network (RSNC) operated by the Istituto Nazionale di Geofisica e Vulcanologia (INGV), were re-picked along with those of the SAPTEX temporary array (2001–2004). For some events located in the Upper Val d'Agri, I also used data from the Eni-Agip oil company seismic network. I computed the VP/VS ratio obtaining a value of 1.83 and I carried out an analysis for the one-dimensional (1D) velocity model that approximates the seismic structure of the study area. After this preliminary analysis, making use of the records obtained in the SeSCAL experiment, I incremented the database by handpicking new arrival times. My final dataset consists of 15,666 P- and 9228 S-arrival times associated to 1047 earthquakes with magnitude ML ≥ 1.5. I computed 162 fault-plane solutions and composite focal mechanisms for closely located events. I investigated stress field orientation inverting focal mechanism belonging to the Lucanian Apennine and the Pollino Range, both areas characterized by more concentrated background seismicity. Moreover, I applied the double difference technique (DD) to improve the earthquake locations. Considering these results and different datasets available in the literature, I carried out a detailed analysis of single sub-areas and of a swarm (November 2008) recorded by SeSCAL array. The relocated seismicity appears more concentrated within the upper crust and it is mostly clustered along the Lucanian Apennine chain. In particular, two well-defined clusters were located in the Potentino and in the Abriola-Pietrapertosa sector (central Lucanian region). Their hypocentral depths are slightly deeper than those observed beneath the chain. I suggest that these two seismic features are representative of the transition from the inner portion of the chain with NE-SW extension to the external margin characterized by dextral strike-slip kinematics. In the easternmost part of the study area, below the Bradano foredeep and the Apulia foreland, the seismicity is generally deeper and more scattered and is associated to the Murge uplift and to the small structures present in the area. I also observed a small structure NE-SW oriented in the Abriola-Pietrapertosa area (activated with a swarm in November 2008) that could be considered to act as a barrier to the propagation of a potential rupture of an active NW-SE striking faults system. Focal mechanisms computed in this study are in large part normal and strike-slip solutions and their tensional axes (T-axes) have a generalized NE-SW orientation. Thanks to denser coverage of seismic stations and the detailed analysis, this study is a further contribution to the comprehension of the seismogenesis and state of stress of the Southern Apennines region, giving important contributions to seismotectonic zoning and seismic hazard assessment.

Partial reconstruction of the trajectories of a discretely observed branching diffusion with immigration and an application to inference

In this treatise we consider finite systems of branching particles where the particles move independently of each other according to d-dimensional diffusions. Particles are killed at a position dependent rate, leaving at their death position a random number of descendants according to a position dependent reproduction law. In addition particles immigrate at constant rate (one immigrant per immigration time). A process with above properties is called a branching diffusion withimmigration (BDI). In the first part we present the model in detail and discuss the properties of the BDI under our basic assumptions. In the second part we consider the problem of reconstruction of the trajectory of a BDI from discrete observations. We observe positions of the particles at discrete times; in particular we assume that we have no information about the pedigree of the particles. A natural question arises if we want to apply statistical procedures on the discrete observations: How can we find couples of particle positions which belong to the same particle? We give an easy to implement 'reconstruction scheme' which allows us to redraw or 'reconstruct' parts of the trajectory of the BDI with high accuracy. Moreover asymptotically the whole path can be reconstructed. Further we present simulations which show that our partial reconstruction rule is tractable in practice. In the third part we study how the partial reconstruction rule fits into statistical applications. As an extensive example we present a nonparametric estimator for the diffusion coefficient of a BDI where the particles move according to one-dimensional diffusions. This estimator is based on the Nadaraya-Watson estimator for the diffusion coefficient of one-dimensional diffusions and it uses the partial reconstruction rule developed in the second part above. We are able to prove a rate of convergence of this estimator and finally we present simulations which show that the estimator works well even if we leave our set of assumptions.

Control of the long-range self-organization of polycyclic aromatic hydrocarbons for device applications

Electronic devices based on organic semiconductors have gained increased attention in nanotechnology, especially applicable to the field of field-effect transistors and photovoltaic. A promising class of materials in this reseach field are polycyclic aromatic hydrocarbons (PAHs). Alkyl substitution of these graphenes results in the selforganization into one-dimensional columnar superstructures and provides solubility and processibility. The nano-phase separation between the π-stacking aromatic cores and the disordered peripheral alkyl chains leads to the formation of thermotropic mesophases. Hexa-peri-hexabenzocoronenes (HBC), as an example for a PAH, exhibits some of the highest values for the charge carrier mobility for mesogens, which makes them promising candidates for electronic devices. Prerequisites for efficient charge carrier transport between electrodes are a high purity of the material to reduce possible trapping sites for charge carriers and a pronounced and defect-free, long-range order. Appropriate processing techniques are required to induce a high degree of aligned structures in the discotic material over macroscopic dimensions. Highly-ordered supramolecular structures of different discotics, in particular, of HBC derivatives have been obtained by solution processing using the zone-casting technique, zone-melting or simple extrusion. Simplicity and fabrication of highly oriented columnar structures over long-range are the most essential advantages of these zone-processing methods. A close relation between the molecular design, self-aggregation and the processing conditions has been revealed. The long-range order achieved by the zone-casting proved to be suitable for field effect transistors (FET).

Charakterisierung von plasmafunktionalisierten Oberflächen mittels ToF-SIMS und multivariaten Analysemethoden

Die grundlegenden Prinzipien und Möglichkeiten der Oberflächencharakterisierung mittels ToF-SIMS (Flugzeit-Sekundärionen Massenspektrometrie) werden an ausgewählten Beispielen aus einem aktuell laufenden und vom BMBF geförderten Verbundforschungsprojekt (Fkz: 03N8022A) zum Thema Nanofunktionalisierung von Grenzflächen vorgestellt. Ein Schwerpunkt innerhalb des Projekts stellen die nichtgeschlossenen Schichtsysteme dar, die entweder über Domänenstrukturen oder einer definierten Einzelfunktionalisierung neuartige funktionelle Oberflächen bereitstellen. Mithilfe der sehr oberflächensensitiven ToF-SIMS Methode sowie der Möglichkeit einer graphischen Darstellung lateraler Molekülionenverteilungen auf funktionalisierten Oberflächen können Informationen über Struktur und Belegungsdichte der Funktionsschicht gewonnen werden. Die Kombination des ToF-SIMS Experimentes und eines multivariaten Algorithmus (partial least squares, PLS) liefert eine interessante Möglichkeit zur quantitativen und simultanen Bestimmung von Oberflächeneigenschaften (Element- und molekulare Konzentrationen sowie Kontaktwinkelwerte). Da das ToF-SIMS Spektrum einer plasmafunktionalisierten Oberfläche im Allgemeinen eine Vielzahl unterschiedlicher Fragmentsignale enthält, lässt eine einfache eindimensionale Korrelation (z.B. CF3 - Fragmentintensität ßà CF3-Konzentration) den größten Teil der im Spektrum prinzipiell enthaltenen Information unberücksichtigt. Aufgrund der großen Anzahl von atomaren und molekularen Signalen, die repräsentativ für die chemische Struktur der analysierten Oberflächen sind, ist es sinnvoll, diese Fülle von Informationen zur Quantifizierung der Oberflächeneigenschaften (Elementkonzentrationen, Kontaktwinkel etc.) zu verwenden. Zusätzlich ermöglicht diese Methode eine quantitative Bestimmung der Oberflächeneigenschaften auf nur µm-großen Bereichen. Das ist vorteilhaft für Untersuchungen chemisch strukturierter Oberflächen, da die Größe der Strukturierung für viele Anwendungen in einem Bereich von mehreren µm liegt. Anhand eines Beispieles aus dem biologisch-medizinischen Fachgebiet, soll der erfolgreiche Einsatz multivariater Modelle aufgezeigt werden. In diesem Experiment wurden menschlichen Bindegewebs- (Fibroblasten) und Pankreaszellen auf plasmafunktionalisiserten Oberflächen kultiviert, um die Beeinflussung der Funktionalisierung auf das Zellwachstum zu untersuchen. Die plasmabehandelten Oberflächen wurden durch die Verwendung von TEM-Gittern mit µm-großen Öffnungen chemisch strukturiert und das Wachstumsverhalten der Zellen beobachtet. Jedem dieser µm-großen Bereiche können mithilfe der multivariaten Modelle quantitative Größen zugeordnet werden (Konzentrationen und Kontaktwinkelwerte), die zur Interpretation des Wachstumsverhaltens der Zellen beitragen.

Hypersonic elastic excitations in soft mesoscopic structures

This dissertation is devoted to the experimental exploration of the propagation of elastic waves in soft mesoscopic structures with submicrometer dimensions. A strong motivation of this work is the large technological relevance and the fundamental importance of the subject. Elastic waves are accompanied by time-dependent fluctuations of local stress and strain fields in the medium. As such, the propagation phase velocities are intimately related to the elastic moduli. Knowledge of the elastic wave propagation directly provides information about the mechanical properties of the probed mesoscopic structures, which are not readily accessible experimentally. On the other hand, elastic waves, when propagating in an inhomogeneous medium with spatial inhomogeneities comparable to their wavelength, exhibit rather rich behavior, including the appearance of novel physical phenomena, such as phononic bandgap formation. So far, the experimental work has been restricted to macroscopic structures, which limit wave propagation below the KHz range. It was anticipated that an experimental approach capable of probing the interplay of the wave propagation with the controlled mesoscopic structures would contribute to deeper insights into the fundamental problem of elastic wave propagation in inhomogeneous systems. The mesoscopic nature of the structures to be studied precludes the use of traditional methods, such as sound transmission, for the study of elastic wave propagation. In this work, an optical method utilizing the inelastic scattering of photons by GHz frequency thermally excited elastic waves, known as Brillouin light scattering spectroscopy (BLS), was employed. Two important classes of soft structures were investigated: thin films and colloidal crystals. For the former, the main interest was the effect of the one-dimensional (1D) confinement on the wave propagation due to the presence of the free-surface or interface of the layer and the utilization of these waves to extract relevant material parameters. For the second system, the primary interest was the interaction of the elastic wave and the strong scattering medium with local resonance units in a three-dimensional (3D) periodic arrangement.

Quantum hydrodynamic models for semiconductors with and without collisions

In dieser Arbeit werden Quantum-Hydrodynamische (QHD) Modelle betrachtet, die ihren Einsatz besonders in der Modellierung von Halbleiterbauteilen finden. Das QHD Modell besteht aus den Erhaltungsgleichungen für die Teilchendichte, das Momentum und die Energiedichte, inklusive der Quanten-Korrekturen durch das Bohmsche Potential. Zu Beginn wird eine Übersicht über die bekannten Ergebnisse der QHD Modelle unter Vernachlässigung von Kollisionseffekten gegeben, die aus ein­em Schrödinger-System für den gemischten-Zustand oder aus der Wigner-Glei­chung hergeleitet werden können. Nach der Reformulierung der eindimensionalen QHD Gleichungen mit linearem Potential als stationäre Schrö­din­ger-Gleichung werden die semianalytischen Fassungen der QHD Gleichungen für die Gleichspannungs-Kurve betrachtet. Weiterhin werden die viskosen Stabilisierungen des QHD Modells be­rück­sich­tigt, sowie die von Gardner vorgeschlagene numerische Viskosität für das {sf upwind} Finite-Differenzen Schema berechnet. Im Weiteren wird das viskose QHD Modell aus der Wigner-Glei­chung mit Fokker-Planck Kollisions-Ope­ra­tor hergeleitet. Dieses Modell enthält die physikalische Viskosität, die durch den Kollision-Operator eingeführt wird. Die Existenz der Lösungen (mit strikt positiver Teilchendichte) für das isotherme, stationäre, eindimensionale, viskose Modell für allgemeine Daten und nichthomogene Randbedingungen wird gezeigt. Die dafür notwendigen Abschätzungen hängen von der Viskosität ab und erlauben daher den Grenzübergang zum nicht-viskosen Fall nicht. Numerische Simulationen der Resonanz-Tunneldiode modelliert mit dem nichtisothermen, stationären, eindimensionalen, viskosen QHD Modell zeigen den Einfluss der Viskosität auf die Lösung. Unter Verwendung des von Degond und Ringhofer entwickelten Quanten-Entropie-Minimierungs-Verfahren werden die allgemeinen QHD-Gleichungen aus der Wigner-Boltzmann-Gleichung mit dem BGK-Kollisions-Operator hergeleitet. Die Herleitung basiert auf der vorsichtige Entwicklung des Quanten-Max­well­ians in Potenzen der skalierten Plankschen Konstante. Das so erhaltene Modell enthält auch vertex-Terme und dispersive Terme für die Ge­schwin­dig­keit. Dadurch bleibt die Gleichspannungs-Kurve für die Re­so­nanz-Tunnel­diode unter Verwendung des allgemeinen QHD Modells in einer Dimension numerisch erhalten. Die Ergebnisse zeigen, dass der dispersive Ge­schwin­dig­keits-Term die Lösung des Systems stabilisiert.

Supramolecular order and dynamics of functional materials studied by solid state NMR

The goal of this thesis was the investigation of the structure, conformation, supramolecular order and molecular dynamics of different classes of functional materials (phthalocyanine, perylene and hexa-peri-hexabenzocoronene derivatives and mixtures of those), all having planar aromatic cores modified with various types of alkyl chains. The planar aromatic systems are known to stack in the solid and the liquid-crystalline state due to p-p interactions forming columnar superstructures with high one-dimensional charge carrier mobility and potential application in photovoltaic devices. The different functionalities attached to the aromatic cores significantly influence the behavior of these systems allowing the experimentalists to modify the structures to fine-tune the desired thermotropic properties or charge carrier mobility. The aim of the presented studies was to understand the interplay between the driving forces causing self-assembly by relating the structural and dynamic information about the investigated systems. The supramolecular organization is investigated by applying 1H solid state NMR recoupling techniques. The results are related with DSC and X-ray scattering data. Detailed information about the site-specific molecular dynamics is gained by recording spinning sideband patterns using 1H-1H and 13C-1H solid state NMR recoupling techniques. The determined dipole-dipole coupling constants are then related with the coupling constants of the respective rigid pairs, thus providing local dynamic order parameters for the respective moieties. The investigations presented reveal that in the crystalline state the preferred arrangement in the columnar stack of discotic molecules modified with alkyl chains is tilted. This leads to characteristic differences in the 1H chemical shifts of otherwise chemically equivalent protons. Introducing branches and increasing the length of the alkyl chains results in lower mesophase transitions and disordered columnar stacks. In the liquid-crystalline state some of the discs lose the tilted orientation, others do not, but all start a rapid rotation about the columnar axis.

Interactive simulation of flexible parts

Computer simulations play an ever growing role for the development of automotive products. Assembly simulation, as well as many other processes, are used systematically even before the first physical prototype of a vehicle is built in order to check whether particular components can be assembled easily or whether another part is in the way. Usually, this kind of simulation is limited to rigid bodies. However, a vehicle contains a multitude of flexible parts of various types: cables, hoses, carpets, seat surfaces, insulations, weatherstrips... Since most of the problems using these simulations concern one-dimensional components and since an intuitive tool for cable routing is still needed, we have chosen to concentrate on this category, which includes cables, hoses and wiring harnesses. In this thesis, we present a system for simulating one dimensional flexible parts such as cables or hoses. The modeling of bending and torsion follows the Cosserat model. For this purpose we use a generalized spring-mass system and describe its configuration by a carefully chosen set of coordinates. Gravity and contact forces as well as the forces responsible for length conservation are expressed in Cartesian coordinates. But bending and torsion effects can be dealt with more effectively by using quaternions to represent the orientation of the segments joining two neighboring mass points. This augmented system allows an easy formulation of all interactions with the best appropriate coordinate type and yields a strongly banded Hessian matrix. An energy minimizing process accounts for a solution exempt from the oscillations that are typical of spring-mass systems. The use of integral forces, similar to an integral controller, allows to enforce exactly the constraints. The whole system is numerically stable and can be solved at interactive frame rates. It is integrated in the DaimlerChrysler in-house Virtual Reality Software veo for use in applications such as cable routing and assembly simulation and has been well received by users. Parts of this work have been published at the ACM Solid and Physical Modeling Conference 2006 and have been selected for the special issue of the Computer-Aided-Design Journal to the conference.

Numerical simulation of interdigitated back contact hetero-junction solar cells

The goal of this thesis is the application of an opto-electronic numerical simulation to heterojunction silicon solar cells featuring an all back contact architecture (Interdigitated Back Contact Hetero-Junction IBC-HJ). The studied structure exhibits both metal contacts, emitter and base, at the back surface of the cell with the objective to reduce the optical losses due to the shadowing by front contact of conventional photovoltaic devices. Overall, IBC-HJ are promising low-cost alternatives to monocrystalline wafer-based solar cells featuring front and back contact schemes, in fact, for IBC-HJ the high concentration doping diffusions are replaced by low-temperature deposition processes of thin amorphous silicon layers. Furthermore, another advantage of IBC solar cells with reference to conventional architectures is the possibility to enable a low-cost assembling of photovoltaic modules, being all contacts on the same side. A preliminary extensive literature survey has been helpful to highlight the specific critical aspects of IBC-HJ solar cells as well as the state-of-the-art of their modeling, processing and performance of practical devices. In order to perform the analysis of IBC-HJ devices, a two-dimensional (2-D) numerical simulation flow has been set up. A commercial device simulator based on finite-difference method to solve numerically the whole set of equations governing the electrical transport in semiconductor materials (Sentuarus Device by Synopsys) has been adopted. The first activity carried out during this work has been the definition of a 2-D geometry corresponding to the simulation domain and the specification of the electrical and optical properties of materials. In order to calculate the main figures of merit of the investigated solar cells, the spatially resolved photon absorption rate map has been calculated by means of an optical simulator. Optical simulations have been performed by using two different methods depending upon the geometrical features of the front interface of the solar cell: the transfer matrix method (TMM) and the raytracing (RT). The first method allows to model light prop-agation by plane waves within one-dimensional spatial domains under the assumption of devices exhibiting stacks of parallel layers with planar interfaces. In addition, TMM is suitable for the simulation of thin multi-layer anti reflection coating layers for the reduction of the amount of reflected light at the front interface. Raytracing is required for three-dimensional optical simulations of upright pyramidal textured surfaces which are widely adopted to significantly reduce the reflection at the front surface. The optical generation profiles are interpolated onto the electrical grid adopted by the device simulator which solves the carriers transport equations coupled with Poisson and continuity equations in a self-consistent way. The main figures of merit are calculated by means of a postprocessing of the output data from device simulation. After the validation of the simulation methodology by means of comparison of the simulation result with literature data, the ultimate efficiency of the IBC-HJ architecture has been calculated. By accounting for all optical losses, IBC-HJ solar cells result in a theoretical maximum efficiency above 23.5% (without texturing at front interface) higher than that of both standard homojunction crystalline silicon (Homogeneous Emitter HE) and front contact heterojuction (Heterojunction with Intrinsic Thin layer HIT) solar cells. However it is clear that the criticalities of this structure are mainly due to the defects density and to the poor carriers transport mobility in the amorphous silicon layers. Lastly, the influence of the most critical geometrical and physical parameters on the main figures of merit have been investigated by applying the numerical simulation tool set-up during the first part of the present thesis. Simulations have highlighted that carrier mobility and defects level in amorphous silicon may lead to a potentially significant reduction of the conversion efficiency.

Numerical coupling of thermal-electric network models and energy-transport equations including optoelectronic semiconductor devices

In this work the numerical coupling of thermal and electric network models with model equations for optoelectronic semiconductor devices is presented. Modified nodal analysis (MNA) is applied to model electric networks. Thermal effects are modeled by an accompanying thermal network. Semiconductor devices are modeled by the energy-transport model, that allows for thermal effects. The energy-transport model is expandend to a model for optoelectronic semiconductor devices. The temperature of the crystal lattice of the semiconductor devices is modeled by the heat flow eqaution. The corresponding heat source term is derived under thermodynamical and phenomenological considerations of energy fluxes. The energy-transport model is coupled directly into the network equations and the heat flow equation for the lattice temperature is coupled directly into the accompanying thermal network. The coupled thermal-electric network-device model results in a system of partial differential-algebraic equations (PDAE). Numerical examples are presented for the coupling of network- and one-dimensional semiconductor equations. Hybridized mixed finite elements are applied for the space discretization of the semiconductor equations. Backward difference formluas are applied for time discretization. Thus, positivity of charge carrier densities and continuity of the current density is guaranteed even for the coupled model.