Monte Carlo simulations of Potts glasses

A complete understanding of the glass transition isstill a challenging problem. Some researchers attributeit to the (hypothetical) occurrence of a static phasetransition, others emphasize the dynamical transitionof mode coupling-theory from an ergodic to a non ergodicstate. A class of disordered spin models has been foundwhich unifies both scenarios. One of these models isthe p-state infinite range Potts glass with p>4, whichexhibits in the thermodynamic limit both a dynamicalphase transition at a temperature T_D, and a static oneat T_0 < T_D. In this model every spins interacts withall the others, irrespective of distance. Interactionsare taken from a Gaussian distribution.In order to understand better its behavior forfinite number N of spins and the approach to thethermodynamic limit, we have performed extensive MonteCarlo simulations of the p=10 Potts glass up to N=2560.The time-dependent spin-autocorrelation function C(t)shows strong finite size effects and it does not showa plateau even for temperatures around the dynamicalcritical temperature T_D. We show that the N-andT-dependence of the relaxation time for T > T_D can beunderstood by means of a dynamical finite size scalingAnsatz.The behavior in the spin glass phase down to atemperature T=0.7 (about 60% of the transitiontemperature) is studied. Well equilibratedconfigurations are obtained with the paralleltempering method, which is also useful for properlyestablishing static properties, such as the orderparameter distribution function P(q). Evidence is givenfor the compatibility with a one step replica symmetrybreaking scenario. The study of the cumulants of theorder parameter does not permit a reliable estimation ofthe static transition temperature. The autocorrelationfunction at low T exhibits a two-step decay, and ascaling behavior typical of supercooled liquids, thetime-temperature superposition principle, is observed. Inthis region the dynamics is governed by Arrheniusrelaxations, with barriers growing like N^{1/2}.We analyzed the single spin dynamics down to temperaturesmuch lower than the dynamical transition temperature. We found strong dynamical heterogeneities, which explainthe non-exponential character of the spin autocorrelationfunction. The spins seem to relax according to dynamicalclusters. The model in three dimensions tends to acquireferromagnetic order for equal concentration of ferro-and antiferromagnetic bonds. The ordering has differentcharacteristics from the pure ferromagnet. The spinglass susceptibility behaves like chi_{SG} proportionalto 1/T in the region where a spin glass is predicted toexist in mean-field. Also the analysis of the cumulantsis consistent with the absence of spin glass orderingat finite temperature. The dynamics shows multi-scalerelaxations if a bimodal distribution of bonds isused. We propose to understand it with a model based onthe local spin configuration. This is consistent with theabsence of plateaus if Gaussian interactions are used.

Partial melting within the lower crust – Constraints from bulk rock geochemical data and trace element mineral analyses of granulites and migmatites from Central Finland

This PhD thesis concerns geochemical constraints on recycling and partial melting of Archean continental crust. A natural example of such processes was found in the Iisalmi area of Central Finland. The rocks from this area are Middle to Late Archean in age and experienced metamorphism and partial melting between 2.7-2.63 Ga. The work is based on extensive field work. It is furthermore founded on bulk rock geochemical data as well as in-situ analyses of minerals. All geochemical data were obtained at the Institute of Geosciences, University of Mainz using X-ray fluorescence, solution ICP-MS and laser ablation-ICP-MS for bulk rock geochemical analyses. Mineral analyses were accomplished by electron microprobe and laser ablation ICP-MS. Fluid inclusions were studied by microscope on a heating-freezing-stage at the Geoscience Center, University Göttingen. Part I focuses on the development of a new analytical method for bulk rock trace element determination by laser ablation-ICP-MS using homogeneous glasses fused from rock powder on an Iridium strip heater. This method is applicable for mafic rock samples whose melts have low viscosities and homogenize quickly at temperatures of ~1200°C. Highly viscous melts of felsic samples prevent melting and homogenization at comparable temperatures. Fusion of felsic samples can be enabled by addition of MgO to the rock powder and adjustment of melting temperature and melting duration to the rock composition. Advantages of the fusion method are low detection limits compared to XRF analyses and avoidance of wet-chemical processing and use of strong acids as in solution ICP-MS as well as smaller sample volumes compared to the other methods. Part II of the thesis uses bulk rock geochemical data and results from fluid inclusion studies for discrimination of melting processes observed in different rock types. Fluid inclusion studies demonstrate a major change in fluid composition from CO2-dominated fluids in granulites to aqueous fluids in TTG gneisses and amphibolites. Partial melts were generated in the dry, CO2-rich environment by dehydration melting reactions of amphibole which in addition to tonalitic melts produced the anhydrous mineral assemblages of granulites (grt + cpx + pl ± amph or opx + cpx + pl + amph). Trace element modeling showed that mafic granulites are residues of 10-30 % melt extraction from amphibolitic precursor rocks. The maximum degree of melting in intermediate granulites was ~10 % as inferred from modal abundances of amphibole, clinopyroxene and orthopyroxene. Carbonic inclusions are absent in upper-amphibolite facies migmatites whereas aqueous inclusion with up to 20 wt% NaCl are abundant. This suggests that melting within TTG gneisses and amphibolites took place in the presence of an aqueous fluid phase that enabled melting at the wet solidus at temperatures of 700-750°C. The strong disruption of pre-metamorphic structures in some outcrops suggests that the maximum amount of melt in TTG gneisses was ~25 vol%. The presence of leucosomes in all rock types is taken as the principle evidence for melt formation. However, mineralogical appearance as well as major and trace element composition of many leucosomes imply that leucosomes seldom represent frozen in-situ melts. They are better considered as remnants of the melt channel network, e.g. ways on which melts escaped from the system. Part III of the thesis describes how analyses of minerals from a specific rock type (granulite) can be used to determine partition coefficients between different minerals and between minerals and melt suitable for lower crustal conditions. The trace element analyses by laser ablation-ICP-MS show coherent distribution among the principal mineral phases independent of rock composition. REE contents in amphibole are about 3 times higher than REE contents in clinopyroxene from the same sample. This consistency has to be taken into consideration in models of lower crustal melting where amphibole is replaced by clinopyroxene in the course of melting. A lack of equilibrium is observed between matrix clinopyroxene / amphibole and garnet porphyroblasts which suggests a late stage growth of garnet and slow diffusion and equilibration of the REE during metamorphism. The data provide a first set of distribution coefficients of the transition metals (Sc, V, Cr, Ni) in the lower crust. In addition, analyses of ilmenite and apatite demonstrate the strong influence of accessory phases on trace element distribution. Apatite contains high amounts of REE and Sr while ilmenite incorporates about 20-30 times higher amounts of Nb and Ta than amphibole. Furthermore, trace element mineral analyses provide evidence for magmatic processes such as melt depletion, melt segregation, accumulation and fractionation as well as metasomatism having operated in this high-grade anatectic area.

Molecular dynamics simulations of silicate and borate glasses and melts : structure, diffusion dynamics and vibrational properties

Molecular dynamics simulations of silicate and borate glasses and melts: Structure, diffusion dynamics and vibrational properties. In this work computer simulations of the model glass formers SiO2 and B2O3 are presented, using the techniques of classical molecular dynamics (MD) simulations and quantum mechanical calculations, based on density functional theory (DFT). The latter limits the system size to about 100−200 atoms. SiO2 and B2O3 are the two most important network formers for industrial applications of oxide glasses. Glass samples are generated by means of a quench from the melt with classical MD simulations and a subsequent structural relaxation with DFT forces. In addition, full ab initio quenches are carried out with a significantly faster cooling rate. In principle, the structural properties are in good agreement with experimental results from neutron and X-ray scattering, in all cases. A special focus is on the study of vibrational properties, as they give access to low-temperature thermodynamic properties. The vibrational spectra are calculated by the so-called ”frozen phonon” method. In all cases, the DFT curves show an acceptable agreement with experimental results of inelastic neutron scattering. In case of the model glass former B2O3, a new classical interaction potential is parametrized, based on the liquid trajectory of an ab initio MD simulation at 2300 K. In this course, a structural fitting routine is used. The inclusion of 3-body angular interactions leads to a significantly improved agreement of the liquid properties of the classical MD and ab initio MD simulations. However, the generated glass structures, in all cases, show a significantly lower fraction of 3-membered planar boroxol rings as predicted by experimental results (f=60%-80%). The largest boroxol ring fraction of f=15±5% is observed in the full ab initio quenches from 2300 K. In case of SiO2, the glass structures after the quantum mechanical relaxation are the basis for calculations of the linear thermal expansion coefficient αL(T), employing the quasi-harmonic approximation. The striking observation is a change change of sign of αL(T) going along with a temperature range of negative αL(T) at low temperatures, which is in good agreement with experimental results.

Constraints on mantle melting from major and trace element systematics in residual abyssal peridotites

ZusammenfassungDie Bildung von mittelozeanischen Rückenbasalten (MORB) ist einer der wichtigsten Stoffflüsse der Erde. Jährlich wird entlang der 75.000 km langen mittelozeanischen Rücken mehr als 20 km3 neue magmatische Kruste gebildet, das sind etwa 90 Prozent der globalen Magmenproduktion. Obwohl ozeanische Rücken und MORB zu den am meisten untersuchten geologischen Themenbereichen gehören, existieren weiterhin einige Streit-fragen. Zu den wichtigsten zählt die Rolle von geodynamischen Rahmenbedingungen, wie etwa Divergenzrate oder die Nähe zu Hotspots oder Transformstörungen, sowie der absolute Aufschmelzgrad, oder die Tiefe, in der die Aufschmelzung unter den Rücken beginnt. Diese Dissertation widmet sich diesen Themen auf der Basis von Haupt- und Spurenelementzusammensetzungen in Mineralen ozeanischer Mantelgesteine.Geochemische Charakteristika von MORB deuten darauf hin, dass der ozeanische Mantel im Stabilitätsfeld von Granatperidotit zu schmelzen beginnt. Neuere Experimente zeigen jedoch, dass die schweren Seltenerdelemente (SEE) kompatibel im Klinopyroxen (Cpx) sind. Aufgrund dieser granatähnlichen Eigenschaft von Cpx wird Granat nicht mehr zur Erklärung der MORB Daten benötigt, wodurch sich der Beginn der Aufschmelzung zu geringeren Drucken verschiebt. Aus diesem Grund ist es wichtig zu überprüfen, ob diese Hypothese mit Daten von abyssalen Peridotiten in Einklang zu bringen ist. Diese am Ozeanboden aufgeschlossenen Mantelfragmente stellen die Residuen des Aufschmelz-prozesses dar, und ihr Mineralchemismus enthält Information über die Bildungs-bedingungen der Magmen. Haupt- und Spurenelementzusammensetzungen von Peridotit-proben des Zentralindischen Rückens (CIR) wurden mit Mikrosonde und Ionensonde bestimmt, und mit veröffentlichten Daten verglichen. Cpx der CIR Peridotite weisen niedrige Verhältnisse von mittleren zu schweren SEE und hohe absolute Konzentrationen der schweren SEE auf. Aufschmelzmodelle eines Spinellperidotits unter Anwendung von üblichen, inkompatiblen Verteilungskoeffizienten (Kd's) können die gemessenen Fraktionierungen von mittleren zu schweren SEE nicht reproduzieren. Die Anwendung der neuen Kd's, die kompatibles Verhalten der schweren SEE im Cpx vorhersagen, ergibt zwar bessere Resultate, kann jedoch nicht die am stärksten fraktionierten Proben erklären. Darüber hinaus werden sehr hohe Aufschmelzgrade benötigt, was nicht mit Hauptelementdaten in Einklang zu bringen ist. Niedrige (~3-5%) Aufschmelzgrade im Stabilitätsfeld von Granatperidotit, gefolgt von weiterer Aufschmelzung von Spinellperidotit kann jedoch die Beobachtungen weitgehend erklären. Aus diesem Grund muss Granat weiterhin als wichtige Phase bei der Genese von MORB betrachtet werden (Kapitel 1).Eine weitere Hürde zum quantitativen Verständnis von Aufschmelzprozessen unter mittelozeanischen Rücken ist die fehlende Korrelation zwischen Haupt- und Spuren-elementen in residuellen abyssalen Peridotiten. Das Cr/(Cr+Al) Verhältnis (Cr#) in Spinell wird im Allgemeinen als guter qualitativer Indikator für den Aufschmelzgrad betrachtet. Die Mineralchemie der CIR Peridotite und publizierte Daten von anderen abyssalen Peridotiten zeigen, dass die schweren SEE sehr gut (r2 ~ 0.9) mit Cr# der koexistierenden Spinelle korreliert. Die Auswertung dieser Korrelation ergibt einen quantitativen Aufschmelz-indikator für Residuen, welcher auf dem Spinellchemismus basiert. Damit kann der Schmelzgrad als Funktion von Cr# in Spinell ausgedrückt werden: F = 0.10×ln(Cr#) + 0.24 (Hellebrand et al., Nature, in review; Kapitel 2). Die Anwendung dieses Indikators auf Mantelproben, für die keine Ionensondendaten verfügbar sind, ermöglicht es, geochemische und geophysikalischen Daten zu verbinden. Aus geodynamischer Perspektive ist der Gakkel Rücken im Arktischen Ozean von großer Bedeutung für das Verständnis von Aufschmelzprozessen, da er weltweit die niedrigste Divergenzrate aufweist und große Transformstörungen fehlen. Publizierte Basaltdaten deuten auf einen extrem niedrigen Aufschmelzgrad hin, was mit globalen Korrelationen im Einklang steht. Stark alterierte Mantelperidotite einer Lokalität entlang des kaum beprobten Gakkel Rückens wurden deshalb auf Primärminerale untersucht. Nur in einer Probe sind oxidierte Spinellpseudomorphosen mit Spuren primärer Spinelle erhalten geblieben. Ihre Cr# ist signifikant höher als die einiger Peridotite von schneller divergierenden Rücken und ihr Schmelzgrad ist damit höher als aufgrund der Basaltzusammensetzungen vermutet. Der unter Anwendung des oben erwähnten Indikators ermittelte Schmelzgrad ermöglicht die Berechnung der Krustenmächtigkeit am Gakkel Rücken. Diese ist wesentlich größer als die aus Schweredaten ermittelte Mächtigkeit, oder die aus der globalen Korrelation zwischen Divergenzrate und mittels Seismik erhaltene Krustendicke. Dieses unerwartete Ergebnis kann möglicherweise auf kompositionelle Heterogenitäten bei niedrigen Schmelzgraden, oder auf eine insgesamt größere Verarmung des Mantels unter dem Gakkel Rücken zurückgeführt werden (Hellebrand et al., Chem.Geol., in review; Kapitel 3).Zusätzliche Informationen zur Modellierung und Analytik sind im Anhang A-C aufgeführt

Geochronology, geochemistry and isotopic compositions of the volcanic rocks from oceanic (Hawaii) and continental (Eifel) intra-plate environments

The Eifel volcanism is part of the Central European Volcanic Province (CEVP) and is located in the Rhenish Massif, close to the Rhine and Leine Grabens. The Quaternary Eifel volcanism appears to be related to a mantle plume activity. However, the causes of the Tertiary Hocheifel volcanism remain debated. We present geochronological, geochemical and isotope data to assess the geotectonic settings of the Tertiary Eifel volcanism. Based on 40Ar/39Ar dating, we were able to identify two periods in the Hocheifel activity: from 43.6 to 39.0 Ma and from 37.5 to 35.0 Ma. We also show that the pre-rifting volcanism in the northernmost Upper Rhine Graben (59 to 47 Ma) closely precede the Hocheifel volcanic activity. In addition, the volcanism propagates from south to north within the older phase of the Hocheifel activity. At the time of Hocheifel volcanism, the tectonic activity in the Hocheifel was controlled by stress field conditions identical to those of the Upper Rhine Graben. Therefore, magma generation in the Hocheifel appears to be caused by decompression due to Middle to Late Eocene extension. Our geochemical data indicate that the Hocheifel magmas were produced by partial melting of a garnet peridotite at 75-90 km depth. We also show that crustal contamination is minor although the magmas erupted through a relatively thick continental lithosphere. Sr, Nd and Pb isotopic compositions suggest that the source of the Hocheifel magmas is a mixing between depleted FOZO or HIMU-like material and enriched EM2-like material. The Tertiary Hocheifel and the Quaternary Eifel lavas appear to have a common enriched end-member. However, the other sources are likely to be distinct. In addition, the Hocheifel lavas share a depleted component with the other Tertiary CEVP lavas. Although the Tertiary Hocheifel and the Quaternary Eifel lavas appear to originate from different sources, the potential involvement of a FOZO-like component would indicate the contribution of deep mantle material. Thus, on the basis of the geochemical and isotope data, we cannot rule out the involvement of plume-type material in the Hocheifel magmas. The Ko’olau Scientific Drilling Project (KSDP) has been initiated in order to evaluate the long-term evolution of Ko’olau volcano and obtain information about the Hawaiian mantle plume. High precision Pb triple spike data, as well as Sr and Nd isotope data on KSDP lavas and Honolulu Volcanics (HVS) reveal compositional source variations during Ko’olau growth. Pb isotopic compositions indicate that, at least, three Pb end-members are present in Ko’olau lavas. Changes in the contributions of each component are recorded in the Pb, Sr and Nd isotopes stratigraphy. The radiogenic component is present, at variable proportion, in all three stages of Ko’olau growth. It shows affinities with the least radiogenic “Kea-lo8” lavas present in Mauna Kea. The first unradiogenic component was present in the main-shield stage of Ko’olau growth but its contribution decreased with time. It has EM1 type characteristics and corresponds to the “Ko’olau” component of Hawaiian mantle plume. The second unradiogenic end-member, so far only sampled by Honololu lavas, has isotopic characteristics similar to those of a depleted mantle. However, they are different from those of the recent Pacific lithosphere (EPR MORB) indicating that the HVS are not derived from MORB-related source. We suggest, instead, that the HVS result from melting of a plume material. Thus the evolution of a single Hawaiian volcano records the geochemical and isotopic changes within the Hawaiian plume.

Nature and origin of ultramafic lamprophyres and carbonatites from the borders of the Labrador Sea

Die Ränder des Labrador Meeres wurden während des späten Neoproterozoikums intensiv von karbonatreichen silikatischen Schmelzen durchsetzt. Diese Schmelzen bildeted sich bei Drucken zwischen ca. 4-6 GPa (ca. 120-180 km Tiefe) an der Basis der kontinentalen Mantel-Lithosphäre. Diese Magmengenerierung steht in zeitlichem und räumlichem Zusammenhang mit kontinentalen Extensionsprozessen, welche zu beiden Seiten des sich öffnenden Iapetus-Ozeans auftraten.

Petrogenesis of dunite, wehrlite and websterite xenoliths from Kimberley, South Africa: Origin as mantle peridotites of cumulates

Dunite, wehrlite and websterite xenoliths occur amongst a large abundance of mantle xenoliths in kimberlites of the Kimberley cluster in South Africa. Up to know they have mostly been neglected. On the basis of texture, major and trace elements, oxygen isotopes as well as Re-Os isotope characteristics, they can be subdivided into two groups. A coarse-grained mantle peridotite group, comprising dunite, wehrlite and websterite xenoliths, that are similar to fertile peridotites and represent upper mantle assemblages that are differently influenced by mantle metasomatism. And a cumulate group, containing fine-grained Fe-rich dunite xenoliths that represent cumulates of flood basalt magmatism related to ~183 Ma Karoo and ~2.7 Ga Ventersdorp events in southern Africa. Dunite, wehrlite and websterite xenoliths have preserved a complex history of melt depletion and metasomatic re-enrichment events, which gives information about the different re-enrichment stages of the subcratonic lithospheric mantle and the spatial differences within the Kaapvaal craton upper mantle. Websterite xenoliths comprise orthopyroxene (40-85 Vol. %), clinopyroxene (5-42 Vol. %), garnet (4-10 Vol. %) and subordinately olivine, while dunite and wehrlite xenoliths contain predominantly olivine (65-100 Vol %) and subordinately orthopyroxene, clinopyroxene and garnet. High melt depletion and a dunitic to harzburgitic protolith composition are reflected by high forsterite (Fo90-92) and high olivine NiO contents (2800-5000 ppm) and high orthopyroxene Mg# (Mg/(Mg+Fe)) of 0.91-0.93. Re-depletion ages of predominantly 2.9 Ga reflect a minimum age of melt depletion. Melt depletion ceased in conjunction with collision of the Kimberley block with the Witwatersrand block ~2.9 Ga ago. Subduction related re-fertilisation of the previously depleted mantle xenoliths is documented by i) amoeboid textured orthopyroxene, clinopyroxene and garnet, which crystallized in schlieren along olivine grain boundaries, ii) high whole-rock SiO2, Al2O3, CaO, TiO2, FeO contents, iii) low oxygen isotope ratios in clinopyroxene and garnet of 4.8-5.4 ‰ and 4.7-5.3 ‰, respectively and iv) trace element compositions of wehrlitic clinopyroxene and garnet in equilibrium with high-pressure partial melts of eclogite. Trace element disequilibrium of orthopyroxene with clinopyroxene and garnet indicates a separate origin for orthopyroxene, on one side as primary mantle orthopyroxene in dunite and wehrlite xenoliths and on the other side as reaction product with Si-rich melts produced by partial melting of eclogite. This reaction triggered replacement of olivine by orthopyroxene in the surrounding mantle and produced the typical Si-rich composition of Kaapvaal mantle peridotites. Partial melting of eclogite at higher temperatures produced a second metasomatic melt with lower SiO2, but higher Al2O3, CaO, FeO, Ti, Zr, Hf and a low oxygen isotope ratio. This melt triggered clinopyroxene and locally garnet and rutile crystallization in percolation veins, replacing olivine and orthopyroxene in the Kaapvaal upper mantle. Additionally, websterite xenoliths have experienced late stage cryptic metasomatism by the host kimberlite melt, changing the trace element composition of clinopyroxene, orthopyroxene and garnet to different extent. Hence websterite and most fertile lherzolite xenoliths have experienced three metasomatic events: i) reaction with high-Si melt, ii) percolation of subduction related silica melt with lower SiO2 content and iii) cryptic metasomatism by kimberlite. In contrast, dunite and wehrlite xenoliths have only experienced the second metasomatic event. They represent mantle lithologies further away from metasomatising agents. The Fe-rich dunites comprise olivine neoblasts with subordinate olivine porphyroclasts and parallel-orientated needles of ilmenite, which may enclose spinel. The lower forsterite and NiO contents of olivine in Fe-rich dunites compared to mantle peridotite xenoliths (Fo87-89 vs. Fo93-95 and 1300-2800ppm vs. 2200-3900 ppm, respectively), rules out a restitic origin. Cr-rich spinels are remnants of the original cumulate mineralogy that survived a late stage metasomatic overprint related to the production of the host kimberlite, producing ilmenite and phlogopite in some samples. Olivine porphyroclasts and neoblasts have different trace element compositions, the latter having high Ti, V, Cr and Ni and low Zn, Zr and Nb contents, indicating contrasting origins for neoblasts and porphyroclasts. The dunites have high 187Os/188Os ratios (0.11-0.15) indicating young (Phanerozoic) model ages for most samples, whereas three samples show isotopic mixtures between Phanerozoic neoblasts and ancient porphyroclastic material. Most Fe-rich dunite xenoliths can be interpreted as cumulates of fractional crystallization of Karoo magmatism, whereas the porphyroclasts are interpreted to be remnants from the much earlier Archaean Ventersdorp magmatism.

UHP metamorphic rocks of the Eastern Rhodope Massif, NE Greece: new constraints from petrology, geochemistry and zircon ages

In this PhD thesis, a multidisciplinary study has been carried out on metagranitoids and paragneisses from the Eastern Rhodope Massif, northern Greece, to decipher the pre-Alpine magmatic and geodynamic evolution of the Rhodope Massif and to correlate the eastern part with the western/central parts of the orogen. The Rhodope Massif, which occupies the major part of NE Greece and S Bulgaria, represents the easternmost part of the Internal Hellenides. It is regarded as a nappe stack of high-grade units, which is classically subdivided into an upper unit and a lower unit, separated by a SSE-NNW trending thrust plane, the Nestos thrust. Recent research in the central Greek Rhodope Massif revealed that the two units correspond to two distinct terranes of different age, the Permo-Carboniferous Thracia Terrane, which was overthrusted by the Late Jurassic/Early Cretaceous Rhodope Terrane. These terranes are separated by the Nestos suture, a composite zone comprising metapelites, metabasites, metagranitoids and marbles, which record high-pressure and even ultrahigh-pressure metamorphism in places. Similar characteristic rock associations were investigated during this study along several well-constrained cross sections in vincity to the Ada, Sidiro and Kimi villages in the Greek Eastern Rhodope Massif. Field evidence revealed that the contact zone of the two terranes in the Eastern Rhodope Massif is characterized by a mélange of metapelites, migmatitic amphibolites/eclogites, strongly sheared orthogneisses and marbles. The systematical occurrence of this characteristic rock association between the terranes implies that the Nestos suture is a continuous belt throughout the Greek Rhodope Massif. In this study, a new UHP locality could be established and for the first time in the Greek Rhodope, metamorphic microdiamonds were identified in situ in their host zircons using Laser-Raman spectroscopy. The presence of the diamonds as well as element distribution patterns of the zircons, obtained by TOF-SIMS, indicate metamorphic conditions of T > 1000 °C and P > 4 GPa. The high-pressure and ultrahigh-pressure rocks of the mélange zone are considered to have formed during the subduction of the Nestos Ocean in Jurassic times at ~150 Ma. Melting of metapelitic rocks at UHP conditions facilitated the exhumation to lower crustal levels. To identify major crust forming events, basement granitoids were dated by LA-SF-ICPMS and SHRIMP-II U-Pb analyses of zircons. The geochronological results revealed that the Eastern Rhodope Massif consists of two crustal units, a structurally lower Permo-Carboniferous unit corresponding to the Thracia Terrane and a structurally upper Late Jurassic/Early Cretaceous unit corresponding to the Rhodope Terrane, like it was documented for the Central Rhodope Massif. Inherited zircons in the orthogneisses from the Thracia Terrane of the Eastern Rhodope Massif indicate the presence of a pre-existing Neoproterozoic and Ordovician-Silurian basement in this region. Triassic magmatism is witnessed by the zircons of few orthogneisses from the easternmost Rhodope Massif and is interpreted to be related to rifting processes. Whole-rock major and trace element analyses indicate that the metagranitoids from both terranes originated in a subduction-related magmatic-arc environment. The Sr-Nd isotope data for both terranes of the Eastern and Central Rhodope Massif suggest a mixed crust-mantle source with variable contributions of older crustal material as already indicated by the presence of inherited zircons. Geochemical and isotopic similarity of the basement of the Thracia Terrane and the Pelagonian Zone implies that the Thracia Terrane is a fragment of a formerly unique Permo-Carboniferous basement, separated by rifting and opening of the Meliata-Maliac ocean system in Triassic times. A branch of the Meliata-Maliac ocean system, the Nestos Ocean, subducted northwards in Late Jurassic times leading to the formation of the Late Jurassic/Early Cretaceous Rhodope magmatic arc on remnants of the Thracia Terrane as suggested by inherited Permo-Carboniferous zircons. The ~150 Ma zircon ages of the orthogneisses from the Rhodope Terrane indicate that subduction-related magmatism and HP/UHP metamorphism occurred during the same subduction phase. Subduction ceased due to the closure of the Nestos Ocean in the Late Jurassic/Early Cretaceous. The post-Jurassic evolution of the Rhodope Massif is characterized by the exhumation of the Rhodope core complex in the course of extensional tectonics associated with late granite intrusions in Eocene to Miocene times.

Granulite facies metamorphism and partial melting processes in the Ivrea Zone, Northern Italy

The Ivrea Zone in northern Italy has been the focus of numerous petrological, geochemical and structural studies. It is commonly inferred to represent an almost complete section through the mid to lower continental crust, in which metamorphism and partial melting of the abundant metapelites was the result of magmatic underplating by a large volume of mantle-derived magma. This study concerns amphibolite and granulite facies metamorphism in the Ivrea Zone with focus on metapelites and metapsammites/metagreywackes from Val Strona di Omegna and metapelites from Val Sesia and Val Strona di Postua, with the aim to better constrain their metamorphic evolution as well as their pressure and temperature conditions via phase equilibria modelling.rnrnIn Val Strona di Omegna, the metapelites show a structural and mineralogical change from mica-schists with the common assemblage bi-mu-sill-pl-q-ilm ± liq at the lowest grades, through metatexitic migmatites (g-sill-bi-ksp-pl-q-ilm-liq) at intermediate grades, to complex diatexitic migmatites (g-sill-ru-bi-ksp-pl-q-ilm-liq) at the highest grades. Within this section several mappable isograds occur, including the first appearance of K-feldspar in the metapelites, the first appearance of orthopyroxene in the metabasites and the disappearance of prograde biotite from the metapelites. The inferred onset of partial melting in the metapelites occurs around Massiola. The prograde suprasolidus evolution of the metapelites is consistent with melting via the breakdown of first muscovite then biotite. Maximum modelled melt fractions of 30–40 % are predicted at the highest grade. The regional metamorphic field gradient in Val Strona di Omegna is constrained to range from conditions of 3.5–6.5 kbar at T = 650–730 °C to P &gt; 9 kbar at T &gt; 900 °C. The peak P–T estimates, particularly for granulite facies conditions, are significantly higher (around 100 °C) than those of most previous studies. In Val Sesia and Val Strona di Postua to the south the exposure is more restricted. P–T estimates for the metapelites are 750–850 °C and 5–6.5 kbar in Val Sesia and approximately 800–900 °C and 5.5–7 kbar in Val Strona di Postua. These results show similar temperatures but lower pressure than metapelites in Val Strona di Omegna. Metapelites in Val Sesia in contact with the Mafic Complex exhibit a metatexitic structure, while in Val Strona di Postua diatexitic structures occur. Further, metapelites at the contact with the Mafic Complex contain cordierite (± spinel) that overprint the regional metamorphic assemblages and are interpreted to have formed during contact metamorphism related to intrusion of the Mafic Complex. The lower pressures in the high-grade rocks in Val Sesia and Val Strona di Postua are consistent with some decompression from the regional metamorphic peak prior to the intrusion of the Mafic Complex, suggesting the rocks followed a clockwise P–T path. In contrast, the metapelites in Val Strona di Omegna, especially in the granulite facies, do not contain any cordierite or any evidence for a contact metamorphic overprint. The extrapolated granulite facies mineral isograds are cut by the rocks of the Mafic Complex to the south. Therefore, the Mafic Complex cannot have caused the regional metamorphism and it is unlikely that the Mafic Complex occurs in Val Strona di Omegna.

Forward modelling of petrological crust-forming processes on the early Earth

Tonalite-trondhjemite-granodiorite (TTG) gneisses form up to two-thirds of the preserved Archean continental crust and there is considerable debate regarding the primary magmatic processes of the generation of these rocks. The popular theories indicate that these rocks were formed by partial melting of basaltic oceanic crust which was previously metamorphosed to garnet-amphibolite and/or eclogite facies conditions either at the base of thick oceanic crust or by subduction processes.rnThis study investigates a new aspect regarding the source rock for Archean continental crust which is inferred to have had a bulk compostion richer in magnesium (picrite) than present-day basaltic oceanic crust. This difference is supposed to originate from a higher geothermal gradient in the early Archean which may have induced higher degrees of partial melting in the mantle, which resulted in a thicker and more magnesian oceanic crust. rnThe methods used to investigate the role of a more MgO-rich source rock in the formation of TTG-like melts in the context of this new approach are mineral equilibria calculations with the software THERMOCALC and high-pressure experiments conducted from 10–20 kbar and 900–1100 °C, both combined in a forward modelling approach. Initially, P–T pseudosections for natural rock compositions with increasing MgO contents were calculated in the system NCFMASHTO (Na2O–CaO–FeO–MgO–Al2O3–SiO2–H2O–TiO2) to ascertain the metamorphic products from rocks with increasing MgO contents from a MORB up to a komatiite. A small number of previous experiments on komatiites showed the development of pyroxenite instead of eclogite and garnet-amphibolite during metamorphism and established that melts of these pyroxenites are of basaltic composition, thus again building oceanic crust instead of continental crust.rnThe P–T pseudosections calculated represent a continuous development of their metamorphic products from amphibolites and eclogites towards pyroxenites. On the basis of these calculations and the changes within the range of compositions, three picritic Models of Archean Oceanic Crust (MAOC) were established with different MgO contents (11, 13 and 15 wt%) ranging between basalt and komatiite. The thermodynamic modelling for MAOC 11, 13 and 15 at supersolidus conditions is imprecise since no appropriate melt model for metabasic rocks is currently available and the melt model for metapelitic rocks resulted in unsatisfactory calculations. The partially molten region is therfore covered by high-pressure experiments. The results of the experiments show a transition from predominantly tonalitic melts in MAOC 11 to basaltic melts in MAOC 15 and a solidus moving towards higher temperatures with increasing magnesium in the bulk composition. Tonalitic melts were generated in MAOC 11 and 13 at pressures up to 12.5 kbar in the presence of garnet, clinopyroxene, plagioclase plus/minus quartz (plus/minus orthopyroxene in the presence of quartz and at lower pressures) in the absence of amphibole but it could not be explicitly indicated whether the tonalitic melts coexisting with an eclogitic residue and rutile at 20 kbar do belong to the Archean TTG suite. Basaltic melts were generated predominantly in the presence of granulite facies residues such as amphibole plus/minus garnet, plagioclase, orthopyroxene that lack quartz in all MAOC compositions at pressures up to 15 kbar. rnThe tonalitic melts generated in MAOC 11 and 13 indicate that thicker oceanic crust with more magnesium than that of a modern basalt is also a viable source for the generation of TTG-like melts and therefore continental crust in the Archean. The experimental results are related to different geologic settings as a function of pressure. The favoured setting for the generation of early TTG-like melts at 15 kbar is the base of an oceanic crust thicker than existing today or by melting of slabs in shallow subduction zones, both without interaction of tonalic melts with the mantle. Tonalitic melts at 20 kbar may have been generated below the plagioclase stability by slab melting in deeper subduction zones that have developed with time during the progressive cooling of the Earth, but it is unlikely that those melts reached lower pressure levels without further mantle interaction.rn

Quantensimulationen physisorbierter Molekülschichten auf Graphit

Mit Hilfe der Pfadintegral-Monte Carlo-Methode werdenPhasendiagramme von physisorbierten Molekülschichten aufGraphit untersucht. Die Verwendung von realistischen Potenzialen sowie dieBehandlung aller translatorischen und rotatorischenFreiheitsgrade erlaubt einen quantitativen Vergleich mit denExperimenten.Krypton-Atome bilden in der Monolage ein kommensurablesGitter mit den Atomen über der Mitte jeder drittenGraphitwabe.Die Vorgänge am Schmelzübergang werden von der Desorptioneiniger Atome dominiert. Die Argon-Schicht auf Graphit ist dagegen inkommensurabel.Zweiatomigen Stickstoff-Moleküle bilden eineorientierungsgeordnete Tieftemperaturphase(Fischgrät-Struktur). Quantenfluktuationen führen zu einer Erniedrigung der mitklassischen Methoden berechneten Phasenübergangstemperaturum 12%.Damit wird der experimentelle Wert von 28 K erreicht.Die Anisotropie und das Dipolmoment von Kohlenmonoxid führenzu einer dipolar geordneten Tieftemperaturphase.Die experimentell nicht geklärte Struktur kann in derQuantensimulation als antiferroelektrischeFischgrät-Struktur identifiziert werden.Der Phasenübergang liegt mit 6 K sehr nahe am Experiment(5.2 K).Für die Argon-Stickstoff-Mischsysteme wird dasPhasendiagramm in der Konzentrations(x)-Temperatur(T)-Ebeneerstellt. Die Übergangstemperaturen decken sich mit denen desExperiments.In Konfigurationen mit zufälliger Teilchenbesetzung weisen die linearen Moleküle ab Argon-Konzentrationen von10% ein Orientierungsglas-Verhalten auf.Durch einen zusätzlichen Teilchenaustausch wird in denMischsystemen die Bildung einer Windrad-Phase ermöglicht, inder die Argon-Atome eine Überstruktur annehmen.Diese Phase wird experimentell imArgon-Kohlenmonoxid-Mischsystem vorgefunden, dessenx-T-Phasendiagramm in guter Übereinstimmung mit denSimulationsergebnissen steht.Die explizite Berücksichtigung der Quantenmechanik in denComputersimulationen liefert wesentliche Beiträge zurKlärung des Phasenverhaltens und der Bestimmung vonÜbergangstemperaturen der Tieftemperaturstrukturen.

Molekulardynamik-Untersuchungen zur atomistischen Struktur & Dynamik von binären Mischgläsern

Mit Hilfe von Molekulardynamik-Computersimulationen werdenin dieser Arbeit die Struktur und Dynamik von Gläsern undSchmelzen der Systeme 'NSx'=(Na2O)(xSiO2), mit x=2,3,5, und'AS2'=(Al2O3)(2SiO2) untersucht. Zur Beschreibung dermikroskopischen Wechselwirkungen dient ein Modellpotenzial,das auf einem effektiven Paarpotenzial aus der Literaturaufbaut. Simuliert wurden Teilchentrajektorien über mehrereNanosekunden im Bereich 6100 K >= T >= 2100 K sowie dieGlasstruktur bei 100 K (NSx) bzw. 300 K (AS2). Das Aufbrechen der tetraedrischen Netzwerkstruktur durch denNetzwerkmodifikator Na2O führt zur Ausbildung einerzusätzlichen Struktur auf intermediären Längenskalen, diebei allen Systemen NSx etwa dem Abstand übernächster Na-bzw. Si-Nachbarn eines Na-Atoms entspricht. Die diffusiveDynamik ist in allen Systeme NSx bis zu drei Größenordnungenschneller als in SiO2 . Sie nimmt mit wachsenderNa-Konzentration zu. Die Na-Diffusion zeigtArrhenius-Verhalten; hierbei vollführen die Na-Atome einaktiviertes Hüpfen durch eine erstarrte Si-O-Matrix. DieZustandsdichten werden bis ca. 20 THz durch dominanteNa-Moden bestimmt. Bei hohen Frequenzen weichen die für SiO2 typischen intratetraedrischen Si-O-Schwingungsmodenauf.Im Gegensatz zu Na2O wird Al2O3 in die tetraedrischeNetzwerkstruktur eingebunden. AS2 zeigt eine überwiegend ausAlO4- und SiO4-Tetraedern verknüpfte Polyederstruktur, beider sich die AlO4-Tetraeder lokal anders anordnen als dieSiO4-Tetraeder, um Ladungsausgleich zu gewährleisten. Esbilden sich typische 3(Si,Al)O4-Bausteine ('3-Cluster'). Sie führen auf intermediären Längenskalen zur Ausbildungeines Al-reichen perkolierenden Netzwerks, das dieSiO4-Struktur durchdringt. Im Vergleich zu SiO2 erhöht sichdie diffusive Dynamik aller Komponenten in AS2 um ca. zweiGrößenordnungen. Die intratetraedrischenAl-O-Steckschwingungsmoden sind wesentlich weicher als die typischen intratetraedrischen Si-O-Moden des SiO2.

Experimentelle Untersuchung von Phasenübergängen im kolloidalen Modellsystem harter Kugeln

Die vorliegende Arbeit behandelt den fluid-kristallinen Phasenübergang sowie den Glasübergang anhand von kolloidalen Hart-Kugel(HK)-Modellsystemen. Die Untersuchungen erfolgen dabei im Wesentlichen mit unterschiedlichen Lichtstreumethoden und daher im reziproken Raum. rnDie Analyse der Kristallisationskinetik zeigt, dass es bei der Kristallisation zu signifikanten Abweichungen vom Bild der klassischen Nukleationstheorie (CNT) kommt. Diese geht von einem einstufigen Nukleationsprozess aus, wohingegen bei den hier durchgeführten Experimenten ein mehrstufiger Prozess beobachtet wird. Vor der eigentlichen Kristallisation kommt es zunächst zur Nukleation einer metastabilen Zwischenphase, sogenannter Precursor. In einer zweiten Stufe erfolgt innerhalb der Precursor die eigentliche Nukleation der Kristallite. rnDurch weitere Analyse und den Vergleich des Kristallisations- und Verglasungsszenarios konnte das Konzept der Precursornukleation auf den Vorgang der Verglasung erweitert werden. Während die Kristallnukleation oberhalb des Glasübergangspunktes zum Erliegen kommt, bleibt der Prozess der Precursornukleation auch bei verglasenden Proben erhalten. Ein Glas erstarrt somit in einem amorphen Zustand mit lokalen Precursorstrukturen. Die Korrelation der gemessenen zeitlichen Entwicklung der strukturellen sowie der dynamischen Eigenschaften zeigt darüber hinaus, dass das bisher unverstandene Ageing-Phänomen von HK-Gläsern mit der Nukleation von Precursorn zusammenhängt.rnEin solches mehrstufiges Szenario wurde bereits in früheren Veröffentlichungen beobachtet. Die im Rahmen dieser Arbeit durchgeführten Messungen ermöglichten erstmals die Bestimmung von Kristallnukleationsratendichten (Kristall-NRD) und Ratendichten für die Precursornukleation bis über den Glasübergangspunkt hinaus. Die Kristall-NRD bestätigen die Resultate aus anderen experimentellen Arbeiten. Die weiteren Analysen der Kristall-NRD belegen, dass die fluid-kristalline Grenzflächenspannung bei der Nukleation entgegen den Annahmen der CNT nicht konstant ist, sondern mit ansteigendem Volumenbruch linear zunimmt. Die Erweiterung der CNT um eine linear zunehmende Grenzflächenspannung ermöglichte eine quantitative Beschreibung der gemessenen Kristall- sowie der Precursor-NRD, was den Schluss zulässt, dass es sich in beiden Fällen um einen Boltzmann-aktivierten Prozess handelt. rnUm die beobachteten Abweichungen des Nukleationsprozesses vom Bild der CNT näher zu untersuchen, wurden die kollektiven Partikeldynamiken in stabilen Fluiden und metastabilen Schmelzen analysiert. Im klassischen Bild wird angenommen, dass die kollektive Partikeldynamik beim Vorgang der Nukleation keine Rolle spielt. Anhand der Resultate zeigen sich Abweichungen in der Dynamik stabiler Fluide und metastabiler Schmelzen. Während die kollektive Partikeldynamik in der stabilen Schmelze von der Struktur entkoppelt ist, tritt oberhalb des Phasenübergangspunktes eine Kopplung von Struktur und Dynamik auf. Dabei treten die Abweichungen zunächst in der Umgebung des ersten Strukturfaktormaximums und somit bei den am stärksten besetzten Moden auf. Mit steigender Unterkühlung nehmen die Anzahl der abweichenden Moden sowie die Stärke der Abweichungen zu. Dieses Phänomen könnte einen signifikanten Einfluss auf den Nukleationsprozess und somit auf die Kristallisationskinetik haben. Die Analyse der Dynamik im stabilen Fluid zeigt darüber hinaus Hinweise auf eine Singularität bei Annäherung an den fluid-kristallinen Phasenübergangspunkt.rnDes Weiteren wurden im Rahmen der vorliegenden Arbeit erstmals Ratendichten für die heterogene Nukleation eines HK-Systems an einer flachen Wand mittels statischer Lichtstreuung (SLS) bestimmt. Die Ergebnisse der Messung zeigen, dass die Nukleationsbarriere der heterogenen Nukleation annähernd Null ist und folglich eine vollständige Benetzung der Wand mit einer kristallinen Monolage vorliegt. Die Erweiterung der Untersuchungen auf gekrümmte Oberflächen in Form von sphärischen Partikeln (Seeds) stellt die erste experimentelle Arbeit dar, die den Einfluss eines Ensembles von Seeds auf die Kristallisationskinetik in HK-Systemen untersucht. Die Kristallisationskinetik und die Mikrostruktur werden abhängig von Größe und Anzahldichte der Seed-Partikel signifikant beeinflusst. In Übereinstimmung mit konfokalmikroskopischen Experimenten und Simulationen spielt dabei das Radienverhältnis der Majoritäts- zur Minoritätskomponente eine entscheidende Rolle.