Investigations on maar-diatreme volcanoes by inversion of magnetic and gravity data from the Eifel area, Germany

A study of maar-diatreme volcanoes has been perfomed by inversion of gravity and magnetic data. The geophysical inverse problem has been solved by means of the damped nonlinear least-squares method. To ensure stability and convergence of the solution of the inverse problem, a mathematical tool, consisting in data weighting and model scaling, has been worked out. Theoretical gravity and magnetic modeling of maar-diatreme volcanoes has been conducted in order to get information, which is used for a simple rough qualitative and/or quantitative interpretation. The information also serves as a priori information to design models for the inversion and/or to assist the interpretation of inversion results. The results of theoretical modeling have been used to roughly estimate the heights and the dip angles of the walls of eight Eifel maar-diatremes — each taken as a whole. Inversemodeling has been conducted for the Schönfeld Maar (magnetics) and the Hausten-Morswiesen Maar (gravity and magnetics). The geometrical parameters of these maars, as well as the density and magnetic properties of the rocks filling them, have been estimated. For a reliable interpretation of the inversion results, beside the knowledge from theoretical modeling, it was resorted to other tools such like field transformations and spectral analysis for complementary information. Geologic models, based on thesynthesis of the respective interpretation results, are presented for the two maars mentioned above. The results gave more insight into the genesis, physics and posteruptive development of the maar-diatreme volcanoes. A classification of the maar-diatreme volcanoes into three main types has been elaborated. Relatively high magnetic anomalies are indicative of scoria cones embeded within maar-diatremes if they are not caused by a strong remanent component of the magnetization. Smaller (weaker) secondary gravity and magnetic anomalies on the background of the main anomaly of a maar-diatreme — especially in the boundary areas — are indicative for subsidence processes, which probably occurred in the late sedimentation phase of the posteruptive development. Contrary to postulates referring to kimberlite pipes, there exists no generalized systematics between diameter and height nor between geophysical anomaly and the dimensions of the maar-diatreme volcanoes. Although both maar-diatreme volcanoes and kimberlite pipes are products of phreatomagmatism, they probably formed in different thermodynamic and hydrogeological environments. In the case of kimberlite pipes, large amounts of magma and groundwater, certainly supplied by deep and large reservoirs, interacted under high pressure and temperature conditions. This led to a long period phreatomagmatic process and hence to the formation of large structures. Concerning the maar-diatreme and tuff-ring-diatreme volcanoes, the phreatomagmatic process takes place due to an interaction between magma from small and shallow magma chambers (probably segregated magmas) and small amounts of near-surface groundwater under low pressure and temperature conditions. This leads to shorter time eruptions and consequently to structures of smaller size in comparison with kimberlite pipes. Nevertheless, the results show that the diameter to height ratio for 50% of the studied maar-diatremes is around 1, whereby the dip angle of the diatreme walls is similar to that of the kimberlite pipes and lies between 70 and 85°. Note that these numerical characteristics, especially the dip angle, hold for the maars the diatremes of which — estimated by modeling — have the shape of a truncated cone. This indicates that the diatreme can not be completely resolved by inversion.

Grenzflächenspannung und Morphologie fließender zweiphasiger Polymermischungen

In der vorliegenden Arbeit wurde die Morphologie von zweiphasigen Polymermischungen unter Scherung in situ mit Hilfe einer Kombination aus optischer Scherzelle, Durchlichtmikroskop und computergestützten CCD-Kamera untersucht. Als Modellblends dienten die unverträglichen, bei Raumtemperatur flüssigen Polymersysteme Polyisobutylen (PIB)/Polydimethylsiloxan (PDMS) (I) und Poly(dimethyl-co-methylphenyl)siloxan/PDMS (II). Alle Komponenten verhalten sich bei den verwendeten Scherraten newtonisch.Eine der wichtigsten Einflussgrößen für die Blendmorphologie ist die Grenzflächenspannung gamma 12. Sie wurde für I und II mit Hilfe der Methode der Tropfenrelaxation (dynamisch) als Funktion der Zeit bestimmt. Diese Methode erlaubt die Messung von gamma 12 für Tropfen der Phase A in B sowie von Tropfen B in A. Bei der Methode des hängenden Tropfens (statisch) muss der Tropfen aus der Phase mit der höheren Dichte bestehen. Wo der Vergleich der beiden Methoden möglich ist, stimmen die Ergebnisse für beide Systeme sehr gut überein. Bei II sind die aus der Tropfenrelaxation erhaltenen gamma 12-Werte der beiden komplementären Zusammensetzungen im Rahmen des Fehlers gleich, bei I zeigt ein PIB-Tropfen in PDMS einen um 40 % niedrigeren Wert als ein PDMS-Tropfen in PIB, dies wird auf die Diffusion von kurzkettigen Anteilen des PDMS in die Grenzschicht zurückgeführt. Die Grenzflächenspannung hängt also unter Umständen auch bei binären Systemen deutlich von der Zusammensetzung ab.Für II wurde die Blendmorphologie über den gesamten Zusammensetzungsbereich untersucht. Die häufig beobachteten cokontinuierlichen Strukturen treten bei keiner Zusammensetzung auf. Die Phaseninversion erkennt man in einer sprunghaften Änderung der Tropfengröße zwischen phiPDMS <= 0,400 und 0,500; zudem lässt sich die Zeitabhängigkeit der Radien durch Auftragung gegen das Produkt aus der Deformation und dem Quadrat des Volumenbruchs der Tropfenphase für 0 <= phiPDMS <= 0,400 sowie 0,500 <= phiPDMS <= 1 normieren. Für I und II wurde die Morphologieentwicklung bei 25 °C nach Vorscherung bei 100 bzw. 50 s-1 und anschließendem Sprung der Scherrate auf deutlich niedrigere Werte als Funktion der Zeit verfolgt. Hierbei erhält man bei genügend langer Messdauer (mindestens 200 000-300 000 Schereinheiten) konstante Tropfengrößen. Zum einen handelt es sich dabei um pseudo-stationäre Werte, die nur durch Koaleszenz bestimmt sind, zum anderen um echte stationäre Radien, die durch gleichzeitig ablaufende Koaleszenz und Zerteilung entstehen. Für I liegen die stationären Mittelwerte auf der Zerteilungskurve, für II hingegen auf der Koaleszenzkurve.Der Einfluss eines grenzflächenwirksamen Additivs wurde anhand von I durch Zugabe des Blockcopolymer PIB-b-PDMS zu PIB untersucht. Der Vergleich des zeitlichen Verlaufs von gamma 12 mit der Morphologieentwicklung zeigt, dass das Additiv eine Stabilisierung der feinen Tropfen/Matrix-Struktur des Blends durch Hinderung der Koaleszenz und nicht durch Reduktion der Grenzflächenspannung bewirkt.

Einfluß von Topologie und Verscherung auf Dynamik und mechanische Parameter bei Homopolymeren

Die FT-Rheologie wird zur Unterscheidung verschiedener Kamm-Topologien in Polystyrollösungen und –schmelzen angewendet. Die Polystyrole werden in Abhängigkeit der Deborahzahl De unter LAOS-Bedingungen vermessen. Die Meßergebnisse zeigen, daß der Schritt von wohldefinierten Systemen (lineare Ketten, Sterne) zu solchen mit statistischer Verteilung wie in Kämmen zu großen Veränderungen sowohl im linearen als auch im nichtlinearen Bereich der rheologischen Messungen führt. Sowohl die Masterkurven als auch die Intenstiäten I3/1 und Phasen Phi3 der Nichtlinearitäten der einzelnen Proben weisen jeweils deutliche Unterschiede untereinander auf. Diese sind durch die bisherigen Ergebnisse noch nicht vollständig mit topologischen Merkmalen in Verbindung zu bringen. Die Messungen wurden mit dem von McLeish eingeführten Pom-pom Modell und daraus weiterentwickelten double convected-Pom-pom Modell (DCPP) simuliert und lieferten gute Übereinstimmung sowie auch Vorhersagen über den experimentell nicht mehr zugänglichen Bereich. Zur Untersuchung des Einflusses von mechanischer Scherung auf die lokale, molekulare Dynamik wird das LAOS-Experiment in situ mit dielektrischer Spektroskopie kombiniert. Dazu wurde eine Apparatur entwickelt, die das hochsensitive ARES-Rheometer mit dem hochauflösenden dielektrischen ALPHA-Analyzer verbindet. Mit dieser Apparatur wurde das Typ-A Polymer 1,4-cis-Polyisopren, mit einem Dipolmoment entlang des Rückgrats, bei oszillatorischer Scherung unter gleichzeitiger Aufnahme eines dielektrischen Spektrums vermessen. Es konnte gezeigt werden, daß die oszillatorische Verscherung weder die charakteristische Relaxationszeit noch die Form des Normal Mode Peaks beeinflußt, wohl aber die dielektrische Stärke Delta epsilon. Diese entspricht der Fläche unter dem e“-Peak und kann mit einer Debye- und einer Cole/Davidson-Funktion angepasst werden. Die Abnahme der dielektrischen Stärke mit zunehmender Scheramplitude kann mit der Orientierungsverteilung der End-zu-End-Vektoren in der Probe erklärt werden.

Untersuchung zur Bestimmung quantitativer Parameter der Lungenventilation mittels CT und MRT

Die zuverlässige Berechnung von quantitativen Parametern der Lungenventilation ist für ein Verständnis des Verhaltens der Lunge und insbesondere für die Diagnostik von Lungenerkrankungen von großer Bedeutung. Nur durch quantitative Parameter sind verlässliche und reproduzierbare diagnostische Aussagen über den Gesundheitszustand der Lunge möglich. Im Rahmen dieser Arbeit wurden neue quantitative Verfahren zur Erfassung der Lungenventilation basierend auf der dynamischen Computer- (CT) und Magnetresonanztomographie (MRT) entwickelt. Im ersten Teil dieser Arbeit wurde die Frage untersucht, ob das Aufblähen der Lunge in gesunden Schweinelungen und Lungen mit Akutem Lungenversagen (ARDS) durch einzelne, diskrete Zeitkonstanten beschrieben werden kann, oder ob kontinuierliche Verteilungen von Zeitkonstanten die Realität besser beschreiben. Hierzu wurden Serien dynamischer CT-Aufnahmen während definierter Beatmungsmanöver (Drucksprünge) aufgenommen und anschließend aus den Messdaten mittels inverser Laplace-Transformation die zugehörigen Verteilungen der Zeitkonstanten berechnet. Um die Qualität der Ergebnisse zu analysieren, wurde der Algorithmus im Rahmen von Simulationsrechnungen systematisch untersucht und anschließend in-vivo an gesunden und ARDS-Schweinelungen eingesetzt. Während in den gesunden Lungen mono- und biexponentielle Verteilungen bestimmt wurden, waren in den ARDS-Lungen Verteilungen um zwei dominante Zeitkonstanten notwendig, um die gemessenen Daten auf der Basis des verwendeten Modells verlässlich zu beschreiben. Es wurden sowohl diskrete als auch kontinuierliche Verteilungen gefunden. Die CT liefert Informationen über das solide Lungengewebe, während die MRT von hyperpolarisiertem 3He in der Lage ist, direkt das eingeatmete Gas abzubilden. Im zweiten Teil der Arbeit wurde zeitlich hochaufgelöst das Einströmen eines 3He-Bolus in die Lunge erfasst. Über eine Entfaltungsanalyse wurde anschließend das Einströmverhalten unter Idealbedingungen (unendlich kurzer 3He-Bolus), also die Gewebeantwortfunktion, berechnet und so eine Messtechnik-unabhängige Erfassung des Einströmens von 3He in die Lunge ermöglicht. Zentrale Fragestellung war hier, wie schnell das Gas in die Lunge einströmt. Im Rahmen von Simulationsrechnungen wurde das Verhalten eines Entfaltungsalgorithmus (basierend auf B-Spline Repräsentationen) systematisch analysiert. Zusätzlich wurde ein iteratives Entfaltungsverfahren eingesetzt. Aus zeitlich hochaufgelösten Messungen (7ms) an einer gesunden und einer ARDS-Schweinelunge konnte erstmals nachgewiesen werden, dass das Einströmen in-vivo in weniger als 0,1s geschieht. Die Ergebnisse zeigen Zeitkonstanten im Bereich von 4ms–50ms, wobei zwischen der gesunden Lungen und der ARDS-Lunge deutliche Unterschiede beobachtet wurden. Zusammenfassend ermöglichen daher die in dieser Arbeit vorgestellten Algorithmen eine objektivere Bestimmung quantitativer Parameter der Lungenventilation. Dies ist für die eindeutige Beschreibung ventilatorischer Vorgänge in der Lunge und somit für die Lungendiagnostik unerlässlich. Damit stehen quantitative Methoden für die Lungenfunktionsdiagnostik zur Verfügung, deren diagnostische Relevanz im Rahmen wissenschaftlicher und klinischer Studien untersucht werden kann.

Realisierung der IceCube MonteCarlo-Produktion im WLCG und Untersuchung von Auswirkungen meteorologischer Parameter auf die Erzeugung atmosphärischer Myonen

Das am Südpol gelegene Neutrinoteleskop IceCube detektiert hochenergetische Neutrinos über die schwache Wechselwirkung geladener und neutraler Ströme. Die Analyse basiert auf einem Vergleich mit Monte-Carlo-Simulationen, deren Produktion global koordiniert wird. In Mainz ist es erstmalig gelungen, Simulationen innerhalb der Architektur des Worldwide LHC Computing Grid (WLCG) zu realisieren, was die Möglichkeit eröffnet, Monte-Carlo-Berechnungen auch auf andere deutsche Rechnerfarmen (CEs) mit IceCube-Berechtigung zu verteilen. Atmosphärische Myonen werden mit einer Rate von über 1000 Ereignissen pro Sekunde aufgezeichnet. Eine korrekte Interpretation dieses dominanten Signals, welches um einen Faktor von 10^6 reduziert werden muss um das eigentliche Neutrinosignal zu extrahieren, ist deswegen von großer Bedeutung. Eigene Simulationen mit der Software-Umgebung CORSIKA wurden durchgeführt um die von Energie und Einfallswinkel abhängige Entstehungshöhe atmosphärischer Myonen zu bestimmen. IceCube Myonraten wurden mit Wetterdaten des European Centre for Medium-Range Weather Forcasts (ECMWF) verglichen und Korrelationen zwischen jahreszeitlichen sowie kurzzeitigen Schwankungen der Atmosphärentemperatur und Myonraten konnten nachgewiesen werden. Zudem wurde eine Suche nach periodischen Effekten in der Atmosphäre, verursacht durch z.B. meteorologische Schwerewellen, mit Hilfe einer Fourieranalyse anhand der IceCube-Daten durchgeführt. Bislang konnte kein signifikanter Nachweis zur Existenz von Schwerewellen am Südpol erbracht werden.

Rheology of arrested colloids : a parameter study using novel experimental methods

The steadily increasing diversity of colloidal systems demands for new theoretical approaches and a cautious experimental characterization. Here we present a combined rheological and microscopical study of colloids in their arrested state whereas we did not aim for a generalized treatise but rather focused on a few model colloids, liquid crystal based colloidal suspensions and sedimented colloidal films. We laid special emphasis on the understanding of the mutual influence of dominant interaction mechanisms, structural characteristics and the particle properties on the mechanical behavior of the colloid. The application of novel combinations of experimental techniques played an important role in these studies. Beside of piezo-rheometry we employed nanoindentation experiments and associated standardized analysis procedures. These rheometric methods were complemented by real space images using confocal microscopy. The flexibility of the home-made setup allowed for a combination of both techniques and thereby for a simultaneous rheological and three-dimensional structural analysis on a single particle level. Though, the limits of confocal microscopy are not reached by now. We show how hollow and optically anisotropic particles can be utilized to quantify contact forces and rotational motions for individual particles. In future such data can contribute to a better understanding of particle reorganization processes, such as the liquidation of colloidal gels and glasses under shear.

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.