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

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

Austauschprozesse an Tropopausenfalten extratropischer Zyklonen

Troposphärisches Ozon ist bekannt als wichtiges Oxidationsmittel und als Vorläufergas hoch reaktiver Radikale. Es zählt zu den wichtigsten Treibhausgasen und wirkt bei hohen Konzentrationen an der Erdoberfläche giftig für alle Lebewesen. Zwar wird der Großteil des troposphärischen Ozons photochemisch produziert, ein erheblicher Anteil hat aber stratosphärischen Ursprung und wird entlang von Tropopausenfalten in Zyklonen in die Troposphäre transportiert. Dieser Transport von Luftmassen aus der Stratosphäre in diernTroposphäre (STT) kann zu einem kurzzeitigen, starken Ozonanstieg am Boden führen und langfristig die Chemie der Troposphäre beeinflussen. Die Quantifizierung des Ozoneintrages und die Identifizierung der dafür verantwortlichen Prozesse ist mit großen Unsicherheiten belastet und ein aktuelles Forschungsthema.rnAufgrund ihrer groben Auflösung ist es mit globalen Modellen nicht möglich, die Details dieser STT-Prozesse zu erfassen. Deshalb wird in dieser Arbeit das Modellsystem MECO(n) genutzt, welches das regionale Atmosphärenchemie- und Klimamodell COSMO/MESSy mit dem globalen Chemie-Klimamodell ECHAM5/MESSy (EMAC) koppelt. Eine einheitliche Prozessparametrisierung ermöglicht konsistente, simultane Simulationen in verschiedenen Modellauflösungen. Ein als Teil dieser Arbeit neu entwickeltes Submodell erlaubt die Initialisierung künstlicher, passiver Tracer in Abhängigkeit verschiedener Variablen. Mit einem auf diese Weise freigesetzten, stratosphärischen Tracer lässt sich Ozon mit stratosphärischer Herkunft von solchem, das photochemisch produziert wurde, unterscheiden.rnIm Rahmen einer Fallstudie werden die Austauschprozesse an einer Tropopausenfalte sowohl aus der Eulerischen, als auch aus der Lagrangeschen Perspektive betrachtet. Die Analyse der STT-Prozesse zeigt, dass Luftmassen aus der Stratosphäre durch turbulente und diabatische Prozesse am Rand der Tropopausenfalte in die Troposphäre gelangen und anschließend bis zum Boden transportiert werden. Diese absinkenden, stratosphärischen Luftmassen führen in den Simulationen zu Ozonanstiegen am Boden, die mit Beobachtungsdaten evaluiert werden können. Es wird gezeigt, dass die Ergebnisse der feiner auflösendenrnModellinstanz gut mit den Messungen übereinstimmen.rnIn einer Lagrangeschen Analyse lassen sich Mischungszeitskalen für STT-Prozesse bestimmen. Es wird gezeigt, dass Luftpakete, die sich länger als zehn Stunden in der Troposphäre aufhalten, diese durch den Eintrag ihrer stratosphärischen Tracereigenschaften beeinflussen und daher nicht vernachlässigbar sind. Eine weitere Studie gibt Aufschluss über die Effektivität der Mischung an Tropopausenfalten: Fast die gesamte Luftmasse, die sich zu einem bestimmten Zeitpunkt in der Tropopausenfalte befindet, gelangt innerhalb von zwei Tagen in die Troposphäre.

Nucleation studies on graphics processing units

Ein System in einem metastabilen Zustand muss eine bestimmte Barriere in derrnfreien Energie überwinden um einen Tropfen der stabilen Phase zu formen.rnHerkömmliche Untersuchungen nehmen hierbei kugelförmige Tropfen an. Inrnanisotropen Systemen (wie z.B. Kristallen) ist diese Annahme aber nicht ange-rnbracht. Bei tiefen Temperaturen wirkt sich die Anisotropie des Systems starkrnauf die freie Energie ihrer Oberfläche aus. Diese Wirkung wird oberhalb derrnAufrauungstemperatur T R schwächer. Das Ising-Modell ist ein einfaches Mo-rndell, welches eine solche Anisotropie aufweist. Wir führen großangelegte Sim-rnulationen durch, um die Effekte, die mit einer endlichen Simulationsbox ein-rnhergehen, sowie statistische Ungenauigkeiten möglichst klein zu halten. DasrnAusmaß der Simulationen die benötigt werden um sinnvolle Ergebnisse zu pro-rnduzieren, erfordert die Entwicklung eines skalierbaren Simulationsprogrammsrnfür das Ising-Modell, welcher auf verschiedenen parallelen Architekturen (z.B.rnGrafikkarten) verwendet werden kann. Plattformunabhängigkeit wird durch ab-rnstrakte Schnittstellen erreicht, welche plattformspezifische Implementierungs-rndetails verstecken. Wir benutzen eine Systemgeometrie die es erlaubt eine Ober-rnfläche mit einem variablen Winkel zur Kristallebene zu untersuchen. Die Ober-rnfläche ist in Kontakt mit einer harten Wand, wobei der Kontaktwinkel Θ durchrnein Oberflächenfeld eingestellt werden kann. Wir leiten eine Differenzialglei-rnchung ab, welche das Verhalten der freien Energie der Oberfläche in einemrnanisotropen System beschreibt. Kombiniert mit thermodynamischer Integrationrnkann die Gleichung benutzt werden, um die anisotrope Oberflächenspannungrnüber einen großen Winkelbereich zu integrieren. Vergleiche mit früheren Mes-rnsungen in anderen Geometrien und anderen Methoden zeigen hohe Überein-rnstimung und Genauigkeit, welche vor allem durch die im Vergleich zu früherenrnMessungen wesentlich größeren Simulationsdomänen erreicht wird. Die Temper-rnaturabhängigkeit der Oberflächensteifheit κ wird oberhalb von T R durch diernKrümmung der freien Energie der Oberfläche für kleine Winkel gemessen. DiesernMessung lässt sich mit Simulationsergebnissen in der Literatur vergleichen undrnhat bessere Übereinstimmung mit theoretischen Voraussagen über das Skalen-rnverhalten von κ. Darüber hinaus entwickeln wir ein Tieftemperatur-Modell fürrndas Verhalten um Θ = 90 Grad weit unterhalb von T R. Der Winkel bleibt bis zu einemrnkritischen Feld H C quasi null; oberhalb des kritischen Feldes steigt der Winkelrnrapide an. H C wird mit der freien Energie einer Stufe in Verbindung gebracht,rnwas es ermöglicht, das kritische Verhalten dieser Größe zu analysieren. Die harternWand muss in die Analyse einbezogen werden. Durch den Vergleich freier En-rnergien bei geschickt gewählten Systemgrößen ist es möglich, den Beitrag derrnKontaktlinie zur freien Energie in Abhängigkeit von Θ zu messen. Diese Anal-rnyse wird bei verschiedenen Temperaturen durchgeführt. Im letzten Kapitel wirdrneine 2D Fluiddynamik Simulation für Grafikkarten parallelisiert, welche u. a.rnbenutzt werden kann um die Dynamik der Atmosphäre zu simulieren. Wir im-rnplementieren einen parallelen Evolution Galerkin Operator und erreichen

Evolutionäre Veränderungen in der Nutzung redox-aktiver Aminosäuren und ihre Ursachen

Gewebe, Zellen und speziell Zellkompartimente unterscheiden sich in ihrer Sauerstoffkonzentration, Stoffwechselrate und in der Konzentration an gebildeten reaktiven Sauerstoffspezies. Um eine mögliche Änderung in der Aminosäurennutzung durch den Einfluss von Sauerstoff und seinen reaktiven Spezies untersuchen zu können wurden, Bereiche bzw. Kompartimente der menschlichen Zelle definiert, die einen Referenzrahmen bildeten und bekannt dafür sind, einen relativ hohen Grad an reaktiven Sauerstoffspezies aufzuweisen. Aus dem Vergleich wurde deutlich, dass vor allem die beiden redox-aktiven und schwefeltragenden Aminosäuren Cystein und Methionin durch eine besondere Verteilung und Nutzung charakterisiert sind. Cystein ist hierbei diejenige Aminosäure mit den deutlichsten Änderungen in den fünf untersuchten Modellen der oxidativen Belastung. In all diesen Modellen war die Nutzung von Cystein deutlich reduziert, wohingegen Methionin in Proteinen des Mitochondriums und der Elektronentransportkette angereichert war. Dieser auf den ersten Blick paradoxe Unterschied zwischen Cystein und Methionin wurde näher untersucht, indem die differenzierte Methioninnutzung in verschiedenen Zellkompartimenten von Homo sapiens charakterisiert wurde.rnDie sehr leicht zu oxidierende Aminosäure Methionin zeigt ein ungewöhnliches Verteilungsmuster in ihrer Nutzungshäufigkeit. Entgegen mancher Erwartung wird Methionin in zellulären Bereichen hoher oxidativer Belastung und starker Radikalproduktion intensiv verwendet. Dieses Verteilungsmuster findet man sowohl im intrazellulären Vergleich, als auch im Vergleich verschiedener Spezies untereinander, was daraufhin deutet, dass es einen lokalen Bedarf an redox-aktiven Aminosäuren gibt, der einen sehr starken Effekt auf die Nutzungshäufigkeit von Methionin ausübt. Eine hohe Stoffwechselrate, die im Allgemeinen mit einer erhöhten Produktion von Oxidantien assoziiert wird, scheint ein maßgeblicher Faktor der Akkumulation von Methionin in Proteinen der Atmungskette zu sein. Die Notwendigkeit, oxidiertes Antioxidans wieder zu reduzieren, findet auch bei Methionin Anwendung, denn zu Methioninsulfoxid oxidiertes Methionin wird durch die Methioninsulfoxidreduktase wieder zu Methionin reduziert. Daher kann die spezifische Akkumulation von Methionin in Proteinen, die verstärkt reaktiven Sauerstoffspezies ausgesetzt sind, als eine systematische Strategie angesehen werden, um andere labile Strukturen vor ungewollter Oxidation zu schützen. rnDa Cystein in allen untersuchten Modellen der oxidativen Belastung und im Besonderen in Membranproteinen der inneren Mitochondrienmembran lebensspannenabhängig depletiert war, wurde dieses Merkmal näher untersucht. Deshalb wurde die Hypothese getestet, ob ein besonderer Redox-Mechanismus der Thiolfunktion für diese selektive Depletion einer im Allgemeinen als harmlos oder antioxidativ geltenden Aminosäure verantwortlich ist. Um den Effekt von Cysteinresten in Membranen nachzustellen, wurden primäre humane Lungenfibroblasten (IMR90) mit diversen Modellsubstanzen behandelt. Geringe Konzentrationen der lipophilen Substanz Dodecanthiol verursachten eine signifikante Toxizität in IMR90-Zellen, die von einer schnellen Zunahme an polyubiquitinierten Proteinen und anderen Indikatoren des proteotoxischen Stresses, wie Sequestosom 1 (P62), HSP70 und HSP90 begleitet wurde. Dieser Effekt konnte spezifisch der Chemie der Thiolfunktion in Membranen zugeordnet werden, da Dodecanol (DOH), Dodecylmethylsulfid (DMS), Butanthiol oder wasserlösliche Thiole weder eine cytotoxische Wirkung noch eine Polyubiquitinierung von Proteinen verursachten. Die Ergebnisse stimmen mit der Hypothese überein, dass Thiole innerhalb von biologischen Membranen als radikalische Kettentransferagentien wirken. Diese Eigenschaft wird in der Polymerchemie durch Nutzung von lipophilen Thiolen in hydrophoben Milieus technisch für die Produktion von Polymeren benutzt. Da die Thiylradikal-spezifische Reaktion von cis-Fettsäuren zu trans-Fettsäuren in 12SH behandelten Zellen verstärkt ablief, kann gefolgert werden, dass 12SH zellulär radikalisiert wurde. In lebenden Organismen kann demnach die Oxidation von Cystein die Schädigung von Membranen beschleunigen und damit Einfallstore für die laterale Radikalisierung von integralen Membranproteinen schaffen, welche möglicherweise der Langlebigkeit abträglich ist, zumindest, wenn sie in der inneren Mitochondrienmembran auftritt.

Mesozoic fold structure of the Otago and Alpine Schist Belt and its implications on the tectonic evolution of South Island, New Zealand

Erneute Untersuchungen der mesozoischen Faltenstruktur des Otago Schiefergürtels, Südinsel, Neuseeland, zeigen, dass diese aus zwei aufeinander folgenden, ähnlichen, asymmetrischen, offenen bis mäßig engen Großfaltengenerationen im km- Größenbereich besteht anstatt aus den vorher angenommenen Decken- oder Halbfalten. Hauptproblem der Großfaltenstruktur sind Zonen von durchgreifender Boudinage, die in der Nähe der Großfaltenscharniere entstanden sind. Vorherige Bearbeiter deuteten diese Zonen als 'starke Verformungszonen' oder Überschiebungszonen. Diese Arbeit zeigt, dass in diesen Zonen nur durch die asymmetrische Faltung die unteren liegenden Schenkel der Großfalten boudiniert und somit häufig die ansonsten typischen Faltenstrukturen des liegenden Schenkels einer symmetrischen Faltung überprägt wurden. Ein weiteres Problem dieser mesozoischen Großfaltenstruktur ist die Überprägung einer Faltengeneration auf eine frühere. Weil die Verkürzungsrichtung der überprägenden Faltengeneration nicht subparallel zur älteren Faltenachse ist, sondern einen Winkel von rund 30 Grad einschließt, ist ein Wechsel von orthogonalen zu koaxialen Interferenzmustern der Kleinfalten beobachtbar. Folglich ist die Orientierung der Scheitellinie einer überprägenden und überprägten Kleinfalte nicht unbedingt subparallel zur Orientierung der Faltenachse der Großfalte trotz zylindrischer Faltung. Im letzten Teil dieser Arbeit wird die Überprägung der mesozoischen Großfaltenstruktur durch das känozoisch entstandene, transpressionale Alpine Störungssystem, das einen zweiseitigen Falten- und Überschiebungsgürtel im Otago und im Nordwesten anschließenden Alpinen Schiefergürtel bildet, beschrieben.