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.

Conditional and transgenic mouse models as tools to study the role of TGFbeta, BMP signaling in skeletal development

Die TGFbeta/BMP Signaltransduktionskaskade ist wichtig für viele Entwicklungsprozesse fast aller embryonaler sowie extraembryonaler Gewebe und sie ist ebenso essentiell bei der Aufrechterhaltung der Homöostase im adulten Organismus. In vielen Mausmodellen und Zellkulturversuchen wurde gezeigt, dass Liganden dieses Signalweges in verschiedene Stadien der Knorpel- und Knochenentwicklung involviert sind. BMPs sind beispielsweise maßgeblich an der frühen Kondensation und Bildung des Knorpels und später an Proliferation und Hypertrophie der Chondrozyten beteiligt. BMPs können ektopisch Knochenbildung auslösen und das Expressionsmuster der Liganden und spezifischen Rezeptoren in der Wachstumsfuge lässt auf eine wichtige Rolle der BMPs in der Wachstumsfuge schließen. Der gezielte knock out der BMP-Rezeptoren Bmpr1a und Bmpr1b in proliferierenden Chondrozyten führt zur Ausbildung einer generellen Chondrodysplasie. Smad1, Smad5 und Smad8 sind die Mediatoren der BMP-Signalkaskade. Im Rahmen der vorliegenden Arbeit sollte die Rolle und Funktion der Smad1- und Smad5-Proteine in der Wachstumsfuge untersucht werden. Hierzu wurden konditionale Smad1-knock out-Mäuse mit einer transgenen Mauslinie gekreuzt, die die Cre-Rekombinase spezifisch in proliferierenden Chondrozyten exprimiert. Diese Mäuse wurden mit und ohne heterozygotem Smad5-Hintergrund charakterisiert. Bei einem knock out von Smad1 allein konnte ein leichte Verkürzung der Wachstumsfuge beobachtet werden, wobei prähypertrophe und hypertrophe Zone gleichermaßen betroffen waren. Dieser Phänotyp war verstärkt in Mäusen mit zusätzlichem heterozygotem Smad5-Hintergrund. Eine Verringerung der Proliferationsrate konnte zusammen mit einer verminderten Ihh-Expression nachgewiesen werden. Zusätzlich konnte anhand von Röntgenaufnahmen eine Dysorganisation der nasalen Region und ein fehlendes nasales Septum beobachtet werden. Produktion und Mineralisation der extrazellulären Matrix waren nicht beeinträchtigt. Um die Rolle der BMP- und TGFbeta-Signalkaskaden während der endochondralen Ossifikation zu vergleichen, wurden transgene Mäuse generiert, in denen die TGFbeta-Signalkaskade spezifisch in proliferierenden Chondrozyten gestört war. Zwei Mauslinien, die ähnliche Phänotypen zeigten, wurden untersucht. Esl1 ist ein TGFbeta-bindendes Protein, von dem man annimmt, dass es die TGFbeta-Signalkaskade inhibieren kann. Esl1-knock out-Mäuse sind kleiner als Wildtypmäuse und die Überexpression von Esl1 in proliferierenden Chondrozyten führt zu einer Verlängerung der Wachstumsfuge und einer verstärkten Proliferationsrate. Knorpelmarker, wie Col2a1 und Sox9 sind in diesen Mäusen herunterreguliert, während Col10a1 und Ihh als Marker für die hypertrophe und prähypertrophe Zone herunterreguliert waren. Dies führt zu der Annahme, dass mehr Zellen in die terminale Differenzierung eintreten. Bei transgenen Mäusen, in denen ein dominant-negativer (dn) TGFbeta-Rezeptor in proliferierenden Chondrozyten überexprimiert wurde, konnte eine verlängerte prähypertrophe Zone, eine erhöhte Ihh-Expression, sowie eine verstärkte Proliferationsrate beobachtet werden. Zusätzlich konnte in homozygoten Tieren ein craniofacialer Phänotyp beschrieben werden, der zu Problemen bei der Nahrungsaufnahme und damit zu einer starken Wachstumsbeeinträchtigung führte. Die BMP- und TGFbeta-Signalkaskaden haben möglicherweise antagonistische Effekte in der Wachstumsfuge. Während der Ausfall von BMP in proliferierenden Chondrozyten aufgrund einer gesunkenen Proliferationsrate zu einer Verkürzung der Wachstumsfuge führte, kann man in Mäusen mit einer Störung der TGFbeta-Signalkaskade eine verstärkte Proliferation in einer daher verlängerten Wachstumsfuge beobachten. Ein weiteres Ziel dieser Arbeit war die Generation einer transgenen Mauslinie, die die Cre-Rekombinase spezifisch in hypertrophen Chondrozyten exprimiert. Promoterstudien mit transgenen Mäusen weisen darauf hin, dass ein putatives AP1-Element, etwa 4 kb vor dem ersten Exon des Col10a1 gelegen, wichtig für die spezifische Expression in hypertrophen Chondrozyten ist. Ein Konstrukt, dass vier Kopien dieses Elements und den basalen Promoter enthält, wurde benutzt, um die Cre-Rekombinase spezifisch zu exprimieren. Diese Mauslinie befindet sich in der Testphase und erste Daten deuten auf eine spezifische Expression der Cre-Rekombinase in hypertrophen Chondrozyten hin.

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.

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.

A simple and novel method to monitor breathing and heart rate in awake and urethane-anesthetized newborn rodents

Rodents are most useful models to study physiological and pathophysiological processes in early development, because they are born in a relatively immature state. However, only few techniques are available to monitor non-invasively heart frequency and respiratory rate in neonatal rodents without restraining or hindering access to the animal. Here we describe experimental procedures that allow monitoring of heart frequency by electrocardiography (ECG) and breathing rate with a piezoelectric transducer (PZT) element without hindering access to the animal. These techniques can be easily installed and are used in the present study in unrestrained awake and anesthetized neonatal C57/Bl6 mice and Wistar rats between postnatal day 0 and 7. In line with previous reports from awake rodents we demonstrate that heart rate in rats and mice increases during the first postnatal week. Respiratory frequency did not differ between both species, but heart rate was significantly higher in mice than in rats. Further our data indicate that urethane, an agent that is widely used for anesthesia, induces a hypoventilation in neonates whilst heart rate remains unaffected at a dose of 1 g per kg body weight. Of note, hypoventilation induced by urethane was not detected in rats at postnatal 0/1. To verify the detected hypoventilation we performed blood gas analyses. We detected a respiratory acidosis reflected by a lower pH and elevated level in CO2 tension (pCO2) in both species upon urethane treatment. Furthermore we found that metabolism of urethane is different in P0/1 mice and rats and between P0/1 and P6/7 in both species. Our findings underline the usefulness of monitoring basic cardio-respiratory parameters in neonates during anesthesia. In addition our study gives information on developmental changes in heart and breathing frequency in newborn mice and rats and the effects of urethane in both species during the first postnatal week.

Analytische Ableitungen der Energie für den "Equation-of-Motion-Coupled-Cluster"-Ansatz zur Berechnung von Ionisationspotentialen

In der vorliegenden Arbeit wird die Theorie der analytischen zweiten Ableitungen für die EOMIP-CCSD-Methode formuliert sowie die durchgeführte Implementierung im Quantenchemieprogramm CFOUR beschrieben. Diese Ableitungen sind von Bedeutung bei der Bestimmung statischer Polarisierbarkeiten und harmonischer Schwingungsfrequenzen und in dieser Arbeit wird die Genauigkeit des EOMIP-CCSD-Ansatzes bei der Berechnung dieser Eigenschaften für verschiedene radikalische Systeme untersucht. Des Weiteren können mit Hilfe der ersten und zweiten Ableitungen vibronische Kopplungsparameter berechnet werden, welche zur Simulation von Molekülspektren in Kombination mit dem Köppel-Domcke-Cederbaum (KDC)-Modell - in der Arbeit am Beispiel des Formyloxyl (HCO2)-Radikals demonstriert - benötigt werden.rnrnDer konzeptionell einfache EOMIP-CC-Ansatz wurde gewählt, da hier die Wellenfunktion eines Radikalsystems ausgehend von einem stabilen geschlossenschaligen Zustand durch die Entfernung eines Elektrons gebildet wird und somit die Problematik der Symmetriebrechung umgangen werden kann. Im Rahmen der Implementierung wurden neue Programmteile zur Lösung der erforderlichen Gleichungen für die gestörten EOMIP-CC-Amplituden und die gestörten Lagrange-Multiplikatoren zeta zum Quantenchemieprogramm CFOUR hinzugefügt. Die unter Verwendung des Programms bestimmten Eigenschaften werden hinsichtlich ihrer Leistungsfähigkeit im Vergleich zu etablierten Methoden wie z.B. CCSD(T) untersucht. Bei der Berechnung von Polarisierbarkeiten und harmonischen Schwingungsfrequenzen liefert die EOMIP-CCSD-Theorie meist gute Resultate, welche nur wenig von den CCSD(T)-Ergebnissen abweichen. Einzig bei der Betrachtung von Radikalen, für die die entsprechenden Anionen nicht stabil sind (z.B. NH2⁻ und CH3⁻), liefert der EOMIP-CCSD-Ansatz aufgrund methodischer Nachteile keine aussagekräftige Beschreibung. rnrnDie Ableitungen der EOMIP-CCSD-Energie lassen sich auch zur Simulation vibronischer Kopplungen innerhalb des KDC-Modells einsetzen.rnZur Kopplung verschiedener radikalischer Zustände in einem solchen Modellpotential spielen vor allem die Ableitungen von Übergangsmatrixelementen eine wichtige Rolle. Diese sogenannten Kopplungskonstanten können in der EOMIP-CC-Theorie besonders leicht definiert und berechnet werden. Bei der Betrachtung des Photoelektronenspektrums von HCO2⁻ werden zwei Alternativen untersucht: Die vertikale Bestimmung an der Gleichgewichtsgeometrie des HCO2⁻-Anions und die Ermittlung adiabatischer Kraftkonstanten an den Gleichgewichtsgeometrien des Radikals. Lediglich das adiabatische Modell liefert bei Beschränkung auf harmonische Kraftkonstanten eine qualitativ sinnvolle Beschreibung des Spektrums. Erweitert man beide Modelle um kubische und quartische Kraftkonstanten, so nähern sich diese einander an und ermöglichen eine vollständige Zuordnung des gemessenen Spektrums innerhalb der ersten 1500 cm⁻¹. Die adiabatische Darstellung erreicht dabei nahezu quantitative Genauigkeit.