Isolation and identification of molecular partners of the proteins encoded by the Drosophila tumor suppressor gene lethal(2)tumorous imaginal discs

Ziel dieser Arbeit war es, die funktionelle Bedeutung des Drosophila melanogaster tumor suppressor Gens lethal(2)tumorous imaginal discs (l(2)tid) durch die Identifikation von molekularen Partnern der vom Gen kodierten Proteine zu etablieren. Mit dem Screening einer Expressionsbibliothek mittels des Hefe-Di-Hybrid-Systems wurde das Protein Patched (Ptc) als ein neues Tid-bindendes Protein identifiziert. Ptc ist ein Zentralregulator der Hedhehog-Signalkette. Diese ist in der Entwicklung konserviert und in manchen humanen Krebsarten verwickelt. Die Tid/Ptc-Interaktion wurde mittels unabhängigen biochemischen Methoden wie dem GST-pulldown-Test oder der Immunopräzipitation überprüft. Außerdem ergaben funktionelle Studien in tumorosen Imaginalscheiben einen möglichen inhibitorischen Effekt von Tid über die Hh Signaltransduktion.Im letzten Teil dieser Arbeit wurde die Interaktion zwischen Tid und dem E-APC-Protein (Adenomatous polyposis coli) bewiesen. Polakis und seine Gruppe zeigten durch Studien mit dem Hefe-Di-Hybrid-System und in vitro, dass das hTid mit dem APC-Protein interagiert. Um dies auch auf Drosophila-Ebene zu überprüfen, wurden Immunopräzipitation-Studien mit den Drosophila-Gegenstücken durchgeführt. Diese Studien zeigen zum ersten Mal eine direkte Interaktion beider Proteine in vivo.

Characterization of igneous terranes by zircon dating: implications for the UHP relicts occurrences and suture identification in the Central Rhodope, Northern Greece

One of the key for the understanding of an orogenic belt is the characterization of the terranes involved and the identification of the suture(s) separating crustal blocks: these are essential information for large-scale paleo-reconstructions. In addition, the structural relationships between the terranes involved in the collisional processes and the eventual UHP relicts may provide first order inputs to exhumation models of subducted rocks. The structure of the Rhodope Massif (northern Greece and southern Bulgaria) results from the stacking of high-grade nappes during a continental collision, which age is comprised between Latest-Jurassic and Early-Cenozoic. UHP and HP relicts, associated with oceanic and ultramafic material, suggest the presence of a dismembered suture zone within the massif. The location of this suture remains unclear; furthermore, up to now, the UHP and eclogitic localities represent isolated spots and no synthesis on their structural position within the massif has been proposed. The first aim of this work is to define the relationships between HP-UHP relicts, crustal blocks, shear zones and amphibolitic material. To achieve this objective, we characterized the accreted blocks in terms of protoliths ages of the orthogneisses mainly along two cross sections on the Greek part of the belt. Geochemical affinities of meta-igneous rocks served as a complementary tool for terrane characterization and geodynamic interpretation. Single-zircon Pb-Pb evaporation and zircon U-Pb SHRIMP dating of orthogneiss protoliths define two groups of intrusion-ages: Permo-Carboniferous and Late Jurassic-Early Cretaceous. Structurally, these two groups correspond to distinct units: the Late Jurassic gneissic complex overthrusts the one bearing the Permo-Carboniferous orthogneisses. Mylonites, eclogites, amphibolites of oceanic affinities, and UHP micaschists, mark a “melange” zone, intensively sheared towards the SW, which separates the two units. Thus, we interpret them as two distinct terranes, the Rhodope and Thracia terranes, separated by the Nestos suture. The correlation of our findings in northern Greece to the Bulgarian part of the Massif suggests a northern rooting of the Nestos Suture. This configuration results of the closure of a marginal oceanic basin of the Tethys system by a north-directed subduction. This interpretation is supported by the geochemical affinities of the orthogneisses: the Late-Jurassic igneous rocks formed by subduction-related magmatism, pprobably the same north-directed subduction that gave rise to the UHP metamorphism of the metasediments of the “melange” zone. It is noteworthy that the UHP-HP relicts seem to be restricted to the contact between the two terranes suggesting that the UHP relicts are exhumed only within the suture zone. Furthermore, the singularity of the suture suggests that the Late-Jurassic subduction explains the occurrence of UHP and eclogite relicts in the Central Rhodope despite the large age range previously attributed the UHP and/or HP stage.

Struktur-Dynamik-Beziehungen in amorphen festen Lösungen

Feste Lösungen homogen dispergierter Wirkstoffmoleküle in amorphen Polymermatrizen sind wichtige Materialien in vielen pharmazeutischen Anwendungen, bei denen eine kontrollierte Abgabe wasserunlöslicher Wirkstoffe in wässrige Systeme eine Rolle spielt. Die intermolekulare Bindungs-stärke zwischen Polymer- und Wirkstoffmolekülgruppen bestimmt die Stabilität der festen Lösung und steuert somit die biologische Aktivität der Wirkstoffmoleküle. In festen Lösungen, die aus acryl-säurehaltigen Copolymeren (Protonendonoren) und basischen Wirkstoffmolekülen (Protonenakzepto-ren) hergestellt werden, sind intermolekulare Wasserstoffbrücken zwischen den Systemkomponenten Triebkraft für die Bildung einer stabilen homogenen Dispersion und für die Entstehung struktureller Merkmale zwischen den Molekülgruppen der Systemkomponenten. Zudem ist die Bindungsstärke der Wasserstoffbrücken im Hinblick auf die kontrollierte Abgabe der Wirkstoffe von Bedeutung. Da dynamische chemische Gleichgewichte bei der Bildung der Wasserstoffbrücken eine wichtige Rolle spielen müssen neben strukturellen Parametern auch dynamische Faktoren beleuchtet werden. Ziel dieser Arbeit ist neben der Ermittlung von intermolekularen Bindungsstärken vor allem die Identifika-tion struktureller Verhältnisse zwischen den Systemkomponenten auf molekularer Ebene. Die Be-stimmung der Abhängigkeit dieser Parameter von der Struktur der verwendeten Polymere und einer Vielzahl weiterer Einflüsse wie z.B. Feuchtigkeit, Lagerdauer oder Wirkstoffkonzentration soll ein kontrolliertes Design fester Lösungen mit definierten anwendungsspezifischen Eigenschaften ermögli-chen. Temperaturabhängige 1H-Festkörper-MAS-NMR (Magic Angle Spinning Nuclear Magnetic Resonance) Experimente an festen Lösungen mit unterschiedlichen Copolymer-Zusammensetzungen weisen die Existenz dynamischer chemischer Gleichgewichte in den komplexen Wasserstoffbrücken-netzwerken nach. Veränderungen in der chemischen Verschiebung und in der Linienform der Reso-nanzlinien acider Protonen erlauben einen tiefen Einblick in die Architektur dieser Netzwerke und legen die Bindungsverhältnisse unter Berücksichtigung der Polymerchemie und der Mobilität der Systemkomponenten dar, wobei die Befunde mithilfe quantenchemischer Rechnungen untermauert werden können. Die Gegenwart acider Protonen ermöglicht einen einfachen 1H-2H-Austausch, wor-aufhin mithilfe rotorsynchronisierter temperaturabhängiger 2H-MAS-NMR Experimente die Wasser-stoffbrückenbindungsstärke bestimmt werden kann. Mit 1H-1H-Korrelationsexperimenten (Doppelquantenspektroskopie) stehen Methoden für die Bestimmung homonuklearer dipolarer 1H-1H-Kopplungen zur Verfügung, die strukturelle Aussagen aufgrund von bevorzugten räumlichen Kontak-ten bestimmter Molekülgruppen ermöglichen. Weiterhin können diese Experimente verwendet werden, um Wasserstoffbrücken zwischen Polymergruppen von Polymer-Wirkstoff-Wasserstoffbrücken zu unterscheiden, wodurch eine quantitative Beschreibung des Bindungsnetzwerks und der Konkurrenz-prozesse zwischen den einzelnen wasserstoffverbrückten Spezies ermöglicht wird. Eine Kristallisation der Wirkstoffmoleküle ist in vielen Anwendungen unerwünscht, da sie die biologische Verfügbarkeit des Wirkstoffs reduzieren. Mit 1H-Festkörper-MAS-NMR Experimenten können kristalline von amorph dispergierten Wirkstoffmolekülen unterschieden werden, wodurch eine Quantifizierung der Destabilisierungsprozesse ermöglicht wird, die durch Exposition der festen Lösungen mit Wasserdampf ausgelöst werden können. Die Zeit- und Konzentrationsabhängigkeit der Wasseraufnahme kann mit NMR-Experimenten verfolgt werden, wobei unterschiedlich mobile Was-serspezies an unterschiedlichen Bindungsorten identifiziert werden können, was zum molekularen Verständnis der Destabilisierungsprozesse beiträgt. Zusätzlich wird die Mobilität der Wirkstoffmole-küle bestimmt, die sich – wie auch die Wirkstoffkonzentration - als wichtige Größe in der Beschrei-bung der Destabilisierung erweist. Aufbauend auf den Beobachtungen wird ein Zusammenhang zwischen der Copolymerzusammensetzung und einer kritischen Wirkstoffkonzentration hergestellt, der für die Anwendungen amorpher fester Lösungen in biologischen Systemen von großer Bedeutung ist.

Transgenic mice as a tool for depression research: examples from the endocannabinoid and the corticotropin-releasing hormone system

Major depression belongs to the most serious and widespread psychiatric disorders in today’s society. There is a great need for the delineation of the underlying molecular mechanisms as well as for the identification of novel targets for its treatment. In this thesis, transgenic mice of the endocannabinoid and the corticotropin-releasing hormone (CRH) system were investigated to determine the putative role of these systems for depression-like phenotypes in mice. In the first part of the thesis, we found that the endocannabinoid system was prominently involved in a brain region-specific and temporally controlled manner in acute as well as in chronic stress processing. Genetic deletion in combination with pharmacological intervention revealed the importance of a fully functional endocannabinoid system for efficient neuroendocrine and behavioral stress coping. Accordingly, cannabinoid type 1 (CB1) receptor-deficient mice displayed several depression-like symptoms and molecular alterations, including “behavioral despair”, stress hormone hypersecretion and decreased glucocorticoid receptor and brain-derived neurotrophic factor expression in the hippocampus. However, the endocannabinoid system was dispensable for the efficacy of currently used antidepressant drugs. To facilitate future endocannabinoid research, a transgenic mouse was generated, which overexpressed the CB1 receptor protein fused to a fluorescent protein. In the second part of the thesis, conditional brain region-specific CRH overexpressing mice were evaluated as a model for pathological chronic CRH hyperactivation. Mutant mice showed aberrant neuroendocrine and behavioral stress coping and hyperarousal due to CRH-induced activation of the noradrenergic system in the brain. Mutant mice appeared to share similarities with naturally occurring endogenous CRH activation in wild-type mice and were sensitive to acute pharmacological blockade of CRH receptor type 1 (CRH-R1). Thus, CRH overexpressing mice serve as an ideal in vivo tool to evaluate the efficacy of novel CRH-R1 antagonists. Together, these findings highlight the potential of transgenic mice for the understanding of certain endo-phenotypes (isolated symptoms) of depression and their molecular correlates.

Hox genes and the regulation of programmed cell death in the embryonic central nervous system of Drosophila melanogaster

This study deals with the function and regulation of programmed cell death, or apoptosis, in the development of the embryonic central nervous system of Drosophila melanogaster. The first part provides a description of apoptosis-deficient embryos, which showed that preventing apoptosis does not cause gross morphological defects in the CNS, as it appears well organized despite the presence of too many cells. An analysis of the incidence and pattern of apoptosis over the course of development discloses a partly very orderly pattern suggesting tight spatio-temporal control, but also reveals random apoptotic cells, which suggests a certain amount of plasticity in the embryo. This analysis also allowed precise identification of some of the dying neural cells in the embryo, and establishment of single cell models for studying regulation of segment-specific apoptosis in the embryonic CNS. In the second part of the work, further investigations into mechanisms controlling segment-specific apoptosis revealed the involvement of two Hox genes, Antennapedia (Antp) and Ultrabithorax (Ubx), in this process. Hox genes control the formation of segment-specific structures in their domains of expression, but also regulate organ and tissue morphogenesis. The study presented here shows that Antp and Ubx play antagonistic roles in motoneuron survival in the embryo. Ubx expression in the CNS is strongly upregulated at a late point in development, when most cells have begun to differentiate. This upregulation shortly precedes Ubx-dependent, segment-specific apoptosis of two differentiated motoneurons. It could further be demonstrated that Antp is required for proper development of the NB7-3 lineage and for survival of the NB7-3 motoneuron in the anterior thoracic segments. In segments where Antp and Ubx expression overlaps, Ubx counteracts the anti-apoptotic function of Antp, resulting in cell death. Thus, these two Hox genes play opposing roles in the survival of differentiated neurons in the late developing nervous system. They thereby contribute to establishment of correct connections between outward-projecting neurons and their targets, which is crucial for the assembly of functional neural circuits, as these have to fulfill region-specific locomotion and sensory requirements along the antero-posterior body axis.

A biomimetic model membrane system on oxidic surfaces designed for the investigation of cytochrome c oxidase and other membrane proteins by fluorescence and waveguide spectroscopy

Die transmembrane Potenzialdifferenz Δφm ist direkt mit der katalytischen Aktivität der Cytochrom c Oxidase (CcO) verknüpft. Die CcO ist das terminale Enzym (Komplex IV) in der Atmungskette der Mitochondrien. Das Enzym katalysiert die Reduktion von O2 zu 2 H2O. Dabei werden Elektronen vom natürlichen Substrat Cytochrom c zur CcO übertragen. Der Eleltronentransfer innerhalb der CcO ist an die Protonentranslokation über die Membran gekoppelt. Folglich bildet sich über der inneren Membrane der Mitochondrien eine Differenz in der Protonenkonzentration. Zusätzlich wird eine Potenzialdifferenz Δφm generiert.rnrnDas Transmembranpotenzial Δφm kann mit Hilfe der Fluoreszenzspektroskopie unter Einsatz eines potenzialemfindlichen Farbstoffs gemessen werden. Um quantitative Aussagen aus solchen Untersuchungen ableiten zu können, müssen zuvor Kalibrierungsmessungen am Membransystem durchgeführt werden.rnrnIn dieser Arbeit werden Kalibrierungsmessungen von Δφm in einer Modellmembrane mit inkorporiertem CcO vorgestellt. Dazu wurde ein biomimetisches Membransystem, die Proteinverankerte Doppelschicht (protein-tethered Bilayer Lipid Membrane, ptBLM), auf einem transparenten, leitfähigem Substrat (Indiumzinnoxid, ITO) entwickelt. ITO ermöglicht den simultanen Einsatz von elektrochemischen und Fluoreszenz- oder optischen wellenleiterspektroskopischen Methoden. Das Δφm in der ptBLM wurde durch extern angelegte, definierte elektrische Spannungen induziert. rnrnEine dünne Hydrogelschicht wurde als "soft cushion" für die ptBLM auf ITO eingesetzt. Das Polymernetzwerk enthält die NTA Funktionsgruppen zur orientierten Immobilisierung der CcO auf der Oberfläche der Hydrogels mit Hilfe der Ni-NTA Technik. Die ptBLM wurde nach der Immobilisierung der CcO mittels in-situ Dialyse gebildet. Elektrochemische Impedanzmessungen zeigten einen hohen elektrischen Widerstand (≈ 1 MΩ) der ptBLM. Optische Wellenleiterspektren (SPR / OWS) zeigten eine erhöhte Anisotropie des Systems nach der Bildung der Doppellipidschicht. Cyklovoltammetriemessungen von reduziertem Cytochrom c bestätigten die Aktivität der CcO in der Hydrogel-gestützten ptBLM. Das Membranpotenzial in der Hydrogel-gestützten ptBLM, induziert durch definierte elektrische Spannungen, wurde mit Hilfe der ratiometrischen Fluoreszenzspektroskopie gemessen. Referenzmessungen mit einer einfach verankerten Dopplellipidschicht (tBLM) lieferten einen Umrechnungsfaktor zwischen dem ratiometrischen Parameter Rn und dem Membranpotenzial (0,05 / 100 mV). Die Nachweisgrenze für das Membranpotenzial in einer Hydrogel-gestützten ptBLM lag bei ≈ 80 mV. Diese Daten dienen als gute Grundlage für künftige Untersuchungen des selbstgenerierten Δφm der CcO in einer ptBLM.

A soft-tetramer model for diblock copolymer melts: structure formation and geometric characterization

The ability of block copolymers to spontaneously self-assemble into a variety of ordered nano-structures not only makes them a scientifically interesting system for the investigation of order-disorder phase transitions, but also offers a wide range of nano-technological applications. The architecture of a diblock is the most simple among the block copolymer systems, hence it is often used as a model system in both experiment and theory. We introduce a new soft-tetramer model for efficient computer simulations of diblock copolymer melts. The instantaneous non-spherical shape of polymer chains in molten state is incorporated by modeling each of the two blocks as two soft spheres. The interactions between the spheres are modeled in a way that the diblock melt tends to microphase separate with decreasing temperature. Using Monte Carlo simulations, we determine the equilibrium structures at variable values of the two relevant control parameters, the diblock composition and the incompatibility of unlike components. The simplicity of the model allows us to scan the control parameter space in a completeness that has not been reached in previous molecular simulations.The resulting phase diagram shows clear similarities with the phase diagram found in experiments. Moreover, we show that structural details of block copolymer chains can be reproduced by our simple model.We develop a novel method for the identification of the observed diblock copolymer mesophases that formalizes the usual approach of direct visual observation,using the characteristic geometry of the structures. A cluster analysis algorithm is used to determine clusters of each component of the diblock, and the number and shape of the clusters can be used to determine the mesophase.We also employ methods from integral geometry for the identification of mesophases and compare their usefulness to the cluster analysis approach.To probe the properties of our model in confinement, we perform molecular dynamics simulations of atomistic polyethylene melts confined between graphite surfaces. The results from these simulations are used as an input for an iterative coarse-graining procedure that yields a surface interaction potential for the soft-tetramer model. Using the interaction potential derived in that way, we perform an initial study on the behavior of the soft-tetramer model in confinement. Comparing with experimental studies, we find that our model can reflect basic features of confined diblock copolymer melts.

Transport of nanoparticles into polymersomes : a minimal model system of particles passage through biological membranes

This thesis focuses on the design and characterization of a novel, artificial minimal model membrane system with chosen physical parameters to mimic a nanoparticle uptake process driven exclusively by adhesion and softness of the bilayer. The realization is based on polymersomes composed of poly(dimethylsiloxane)-b-poly(2-methyloxazoline) (PMDS-b-PMOXA) and nanoscopic colloidal particles (polystyrene, silica), and the utilization of powerful characterization techniques. rnPDMS-b-PMOXA polymersomes with a radius, Rh ~100 nm, a size polydispersity, PD = 1.1 and a membrane thickness, h = 16 nm, were prepared using the film rehydratation method. Due to the suitable mechanical properties (Young’s modulus of ~17 MPa and a bending modulus of ~7⋅10-8 J) along with the long-term stability and the modifiability, these kind of polymersomes can be used as model membranes to study physical and physicochemical aspects of transmembrane transport of nanoparticles. A combination of photon (PCS) and fluorescence (FCS) correlation spectroscopies optimizes species selectivity, necessary for a unique internalization study encompassing two main efforts. rnFor the proof of concepts, the first effort focused on the interaction of nanoparticles (Rh NP SiO2 = 14 nm, Rh NP PS = 16 nm; cNP = 0.1 gL-1) and polymersomes (Rh P = 112 nm; cP = 0.045 gL-1) with fixed size and concentration. Identification of a modified form factor of the polymersome entities, selectively seen in the PCS experiment, enabled a precise monitor and quantitative description of the incorporation process. Combining PCS and FCS led to the estimation of the incorporated particles per polymersome (about 8 in the examined system) and the development of an appropriate methodology for the kinetics and dynamics of the internalization process. rnThe second effort aimed at the establishment of the necessary phenomenology to facilitate comparison with theories. The size and concentration of the nanoparticles were chosen as the most important system variables (Rh NP = 14 - 57 nm; cNP = 0.05 - 0.2 gL-1). It was revealed that the incorporation process could be controlled to a significant extent by changing the nanoparticles size and concentration. Average number of 7 up to 11 NPs with Rh NP = 14 nm and 3 up to 6 NPs with Rh NP = 25 nm can be internalized into the present polymersomes by changing initial nanoparticles concentration in the range 0.1- 0.2 gL-1. Rapid internalization of the particles by polymersomes is observed only above a critical threshold particles concentration, dependent on the nanoparticle size. rnWith regard possible pathways for the particle uptake, cryogenic transmission electron microscopy (cryo-TEM) has revealed two different incorporation mechanisms depending on the size of the involved nanoparticles: cooperative incorporation of nanoparticles groups or single nanoparticles incorporation. Conditions for nanoparticle uptake and controlled filling of polymersomes were presented. rnIn the framework of this thesis, the experimental observation of transmembrane transport of spherical PS and SiO2 NPs into polymersomes via an internalization process was reported and examined quantitatively for the first time. rnIn a summary the work performed in frames of this thesis might have significant impact on cell model systems’ development and thus improved understanding of transmembrane transport processes. The present experimental findings help create the missing phenomenology necessary for a detailed understanding of a phenomenon with great relevance in transmembrane transport. The fact that transmembrane transport of nanoparticles can be performed by artificial model system without any additional stimuli has a fundamental impact on the understanding, not only of the nanoparticle invagination process but also of the interaction of nanoparticles with biological as well as polymeric membranes. rn

Identification of target genes of the T-box proteins in Drosophila melanogaster

In a prior bioinformatic analysis by Hüyseyin Binbas, potential Tbx targets sequences in wing-related genes have been identified. Guided by this information, enhancer trap/reporter lacZ insertions were characterized by X-gal staining first in wildtype and then in l(1)omb imaginal discs.rnIn several lines I observed an increase in reporter expression in a l(1)omb mutant background. Since Omb is assumed to function predominantly as a transcriptional repressor, this may indicate direct regulation. Repression by Omb was observed e.g. for brk and tkv. These genes are negatively regulated by Dpp, while omb is induced by Dpp. Omb which mediates the effects of Dpp on proliferation could, thus, also mediate the Dpp effect on patterning of the wing disc. However, brk and tkv were not completely derepressed in l(1)omb indicating that Dpp represses these genes also by an Omb-independent mechanism.rnMore frequently I observed loss of reporter expression in an l(1)omb mutant background. In these cases, regulation by Omb presumably is indirect. For example, STAT92E-lacZ expression in the wildtype eye was symmetrically expressed at the dorsal and ventral margins. In l(1)omb, ventral expression was selectively lost. Loss of omb is known to cause ventral overproliferation of the eye by activation of the Jak/STAT pathway. STAT92E expression is negatively regulated by Jak/STAT signaling suggesting that loss of omb activates Jak/STAT further upstream in the pathway.rnRegional overproliferation of eye and wing in the l(1)omb mutant background proved a complicating issue in the search for Omb targets. This effect made it difficult to decide whether an expanded reporter expression pattern was due to tissue expansion or reporter gene derepression. For instance hth-lacZ appeared to expand along the ventral eye disc margin in l(1)omb. Without addtional experiments it cannot be concluded whether this is due to de-repression or to activation in association with the proliferative state. Parallel to my experiments, evidence accumulated in our laboratory that loss of omb may attenuate Wg and Hegehog signaling. Since these diffusible proteins are the main patterning molecules in the wing imaginal disc, with dpp being downstream of Hh, many of the observed effects could be secondary to reduced Wg and Hh activity. Examples are ab-lacZ, Dll-lacZ and vgBE-lacZ (reduced expression on the dorso-ventral boundary) and inv-lacZ (late larval expression in the anterior wing disc compartment is lost) or sal-lacZ. Epistasis experiment will be required to clarifiy these issues.rnFurthermore, loss of omb appeared to induce cell fate changes. It was reported previously that in an omb null mutant, the dorsal determinant apterous (ap) is ectopically expressed in the ventral compartment (an effect I did not observe with the strongly hypomorphic l(1)omb15, indicating strong dose dependence). Ventral repression of ap is maintained by epigenetic mechanisms. The patchy and variable nature of ectopic expression of ap or grn-1.1-lacZ points to an effect of omb on epigenetic stability.rnIn the second part of my thesis, an analysis of Omb expression in the Drosophila embryonic ventral nervous system was performed. Omb was found co-expressed with Eve in the medial aCC and RP2 motorneurons as well as the fpCC interneuron and the mediolateral CQ neurons. Additionally, Omb was detected in the Eg positive NB7-3 GW serotonergic motoneuron and the N2-4 neurons. Omb was not found in Repo positive glial cells. During embryonic stage 14, Omb showed some coepression with Dpn or Pros. At the embryonic stage 16, Omb was expressed in minor subset of Mid and Wg positive cells.

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.