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.

Strukturen und Phasen des beta-Eukryptits sowie die Sammlung von Beugungsdaten mit axialen q-Scans

Der 'gestopfte Hochquarz' ß-Eukryptit (LiAlSiO4) ist bekannt für seine außergewöhnliche anisotrope Li-Ionenleitfähigkeit und die nahe Null liegende thermische Ausdehnung.Untersucht wurde die temperaturabhängige ß-Eukryptit-Phasenabfolge, insbesondere die modulierte Phase. Deren Satellitenreflexe sind gegenüber den normalen Reflexen erheblich verbreitert, überlappen miteinander sowie mit den dazwischen liegenden 'a-Reflexen' zu Tripletts. Für die Separation der Triplett-Intensitäten waren bisherige Standardverfahren zur Beugungsdatensammlung ungeeignet. Mit 'axialen q-Scans' wurde ein neuartiges Verfahren entwickelt. Intensitäten wurden mit dem neu-entwickelten least squares-Programm GKLS aus 2000 Profilen seriell und automatisch gewonnen und erfolgreich auf Standarddaten skaliert. Die Verwendung verbreiterter Reflexprofile erwies sich als zulässig.Die Verbreiterung wurde auf eine verminderte Fernordnung der Modulation von 11 bis 16 Perioden zurückgeführt (Analyse mit der Gitterfunktion), womit ein ungewöhnliches beugungswinkelabhängiges Verhalten der Reflexbreiten korrespondiert und mit typischen Antiphasendomänendurchmessern (andere Autoren) korreliert.Eine verminderte Si-/ Al-Ordnung wird als ursächlich für geringe Domänengrößen und Fernordnung angesehen, sowie für Eigenschaften wie z.B. a/c-Verhältnisse, Ausdehnungskoeffizienten, Ionenleitfähigkeit, Strukturtyp und Umwandlungstemperaturen. Änderungen des SiO2-Gehaltes, der Temperatur oder der Si- /Al-Ordnung zeitigen für einige Eigenschaften ähnliche Wirkungen.Die gemittelte Struktur der modulierten Phase wurde erstmals zuverlässig bestimmt, die Rolle des Li charakterisiert, Zweifel an der hexagonalen Symmetrie des ß-Eukryptits wurden ausgeräumt und die Bestimmung der modulierten Struktur wurde weitgehend vorbereitet.

Quantification of deformation processes in the Torlesse accretionary wedge, New Zealand

Zusammenfassung:In dieser Studie werden Deformationsprozesse im mesozoischen Torlesse Akkretionskeil (Neuseeland) quantifiziert, um Aufschluß über die Dynamik in Akkretionskeilen zu erhalten. Absolute und relative Verformungsmessungen zeigen sowohl im lokalen als auch regionalen Maßstab eine stark heterogene Deformation des Torlesse Keils. Die regionale Deformation wurde mit Hilfe einer Tensordurchschnittsberechnung, unter Benutzung einzelner lokaler Verformungsdaten, als uniaxiale Verkürzung entlang einer subvertikalen, maximalen Verkürzungsachse charakterisiert. Absolute Verformungsmessungen an niedriggradigen Metasandsteinen belegen darüber hinaus durchschnittliche Volumenverluste von ca. 20% SiO2. Volumenveränderungen in tieferkrustalen Aufschlüssen wurden mittels einer geochemischen Massenbilanzanalyse abgeschätzt. Chemische Zusammensetzungen höhergradiger Zonen weichen je nach Grad der Volumenverformung von der Protolitzusammensetzung ab und zeigen somit Verluste von 15% SiO2 an. Da Speicherorte für das gelöste Material nicht bekannt sind, muss angenommen werden, dass das Material aus dem Keil abtransportiert wurde. Die Verformungsergebnisse geben weiterhin Aufschluß über den Grad der Kopplung zwischen Akkretionskeil und subduzierter Platte. Die ermittelten Scherwerte in den Gesteinen liegen deutlich unter den zu erwartenden Scherwerten, die mittels eines einfachen Modells berechnet wurden, das sowohl verschiedene Konvergenzgeschwindigkeiten als auch Exhumierungsraten berücksichtigt. Dies belegt, dass der Torlesse Keil stark von der subduzierten pazifischen Platte entkoppelt war und die Deformation hauptsächlich durch den Fluß der Sedimente in und aus dem Keil bestimmt wurde.

Comparison of model potentials for molecular dynamics simulation of crystalline silica

Simulationen von SiO2 mit dem von van Beest, Kramer und vanSanten (BKS) entwickelten Paarpotenzial erzeugen vielezufriedenstellende Ergebnisse, aber auch charakteristischeSchwachstellen. In dieser Arbeit wird das BKS-Potenzial mitzwei kürzlich vorgeschlagenen Potenzialen verglichen, dieeffektiv Mehrteilchen-Wechselwirkungen beinhalten. Der ersteAnsatz erlaubt dazu fluktuierende Ladungen, der zweiteinduzierbare Polarisierungen auf den Sauerstoffatomen. Die untersuchten Schwachstellen des BKS Potenzialsbeinhalten das Verhältnis der zwei Gitterkonstanten a und cim Quarzübergang, das von BKS falsch beschrieben wird.Cristobalit und Tridymit erscheinen instabil mit BKS.Weiterhin zeigt die BKS-Zustandsdichte charakteristischeAbweichungen von der wahren Zustandsdichte. DerÜbergangsdruck für den Stishovit I-II Übergang wird deutlichüberschätzt. Das Fluktuierende-Ladungs-Modell verbesserteinige der genannten Punkte, reproduziert aber viele andereEigenschaften schlechter als BKS. DasFluktierende-Dipol-Modell dagegen behebt alle genanntenArtefakte. Zusätzlich wird der druckinduzierte Phasenübergang imalpha-Quarz untersucht. Alle Potentiale finden die selbeStruktur für Quarz II. Bei anschliessender Dekompressionerzeugt BKS eine weitere Phase, während die beiden anderenPotentiale wieder zum alpha-Quarz zurückkehren. Weiterhinwerden zwei Methoden entwickelt, um die piezoelektrischenKonstanten bei konstantem Druck zu bestimmen. Die Ergebnissegeben Hinweise auf eine möglicherweisenicht-elektrostatische Natur der Polarisierungen imFluktuierende-Dipole-Modell. Mit dieser Interpretation scheint das Fluktuierende-DipolPotential alle verfügbaren experimentellen Daten am bestenvon allen drei untersuchten Ansätzen zu reproduzieren.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Molekulardynamik-Computersimulation einer amorph-kristallinen SiO 2 Grenzschicht

In dieser Arbeit werden Molekulardynamik-Computersimulationen zur Untersuchung der statischen und dynamischen Eigenschaften einer amorph/kristallinen Siliziumdioxid(SiO2)-Grenzschicht durchgefuehrt.Die Grenzflaeche wird von der [100]-Ebene des beta-Kristobalit-Kristalls und der fluessigen SiO2-Phase gebildet und in einem Temperaturbereich zwischen 2900K und 3100K im Zustand eines metastabilen Gleichgewichts untersucht. Als Modellpotential zur Beschreibung der mikroskopischen Wechselwirkungen zwischen den Teilchen wird ein einfaches Paarpotential aus der Literatur verwendet, das sowohl die Struktur der kristallinen Phase als auch die der fluessigen Phase gut reproduziert. Bezogen auf die Dichte und die potentielle Energie der Teilchen erstreckt sich der Uebergang von der fluessigen in die kristalline Phase ueber 3-5 Atomlagen. Ein Layering-Effekt der Dichte in der fluessigen Phase in der Naehe der Grenzschicht wird nicht beobachtet. Der Einfluss der Grenzschicht auf statische Groessen, welche das System auf einer mittelreichweitigen Laengenskala beschreiben (z. B. Koordinationszahlverteilung und Ringverteilung) reicht im Vergleich dazu weiter in die fluessige Phase hinein und manifestiert sich in Defektstrukturen, wie z. B. der Erhoehung der Wahrscheinlichkeit fuer das Auftreten von 5-fach koordiniertem Silizium und der vermehrten Bildung von 2er-Ringen in der Fluessigkeit. Dies beguenstigt das Aufbrechen und Umklappen von Si-O-Bindungen und fuehrt zu einer Beschleunigung der Dynamik und einer Erhoehung der Diffusionsgeschwindigkeit in der Fluessigkeit. Im weiteren wird die Hochfrequenzdynamik der reinen SiO2-Fluessigkeit untersucht. Dazu berechnen wir die vibratorische Zustandsdichte in harmonischer Naeherung aus der inhaerenten Struktur. Wir finden einen stark ausgepraegten Peak bei einer Frequenz von 0.6 THz. Dieser Peak kann der niederenergetischsten transversalen akustischen Mode zugeordnet werden, die auch als Scherschwingung des Systems direkt sichtbar ist.

Petrology and geochemistry of Lipari Island (Aeolian archipelago): constraints on magma genesis and evolution

A full set of geochemical and Sr, Nd and Pb isotope data both on bulk-rock and mineral samples is provided for volcanic rocks representative of the whole stratigraphic succession of Lipari Island in the Aeolian archipelago. These data, together with petrographic observations and melt/fluid inclusion investigations from the literature, give outlines on the petrogenesis and evolution of magmas through the magmatic and eruptive history of Lipari. This is the result of nine successive Eruptive Epochs developing between 271 ka and historical times, as derived from recentmost volcanological and stratigraphic studies, combined with available radiometric ages and correlation of tephra layers and marine terrace deposits. These Eruptive Epochs are characterized by distinctive vents partly overlapping in space and time, mostly under control of the main regional tectonic trends (NNW-SSE, N-S and minor E-W). A large variety of lava flows, scoriaceous deposits, lava domes, coulees and pyroclastics are emplaced, ranging in composition through time from calcalkaline (CA) and high-K (HKCA) basaltic andesites to rhyolites. CA and HKCA basaltic andesitic to dacitic magmas were erupted between 271 and 81 ka (Eruptive Epochs 1-6) from volcanic edifices located along the western coast of the island (and subordinately the eastern Monterosa) and the M.Chirica and M.S.Angelo stratocones. These mafic to intermediate magmas mainly evolved through AFC and RAFC processes, involving fractionation of mafic phases, assimilation of wall rocks and mixing with newly injected mafic magmas. Following a 40 ka-long period of volcanic quiescence, the rhyolitic magmas were lately erupted from eruptive vents located in the southern and north-eastern sectors of Lipari between 40 ka and historical times (Eruptive Epochs 7-9). They are suggested to derive from the previous mafic to intermediate melts through AFC processes. During the early phases of rhyolitic magmatism (Eruptive Epochs 7-8), enclaves-rich rocks and banded pumices, ranging in composition from HKCA dacites to low-SiO2 rhyolites were erupted, representing the products of magma mixing between fresh mafic magmas and the fractionated rhyolitic melts. The interaction of mantle-derived magmas with the crust represents an essential process during the whole magmatic hystory of Lipari, and is responsible for the wide range of observed geochemical and isotopic variations. The crustal contribution was particularly important during the intermediate phases of activity of Lipari when the cordierite-bearing lavas were erupted from the M. S.Angelo volcano (Eruptive Epoch 5, 105 ka). These lavas are interpreted as the result of mixing and subsequent hybridization of mantle-derived magmas, akin to the ones characterizing the older phases of activity of Lipari (Eruptive Epochs 1-4), and crustal anatectic melts derived from dehydration-melting reactions of metapelites in the lower crust. A comparison between the adjacent islands of Lipari and Vulcano outlines that their mafic to intermediate magmas seem to be genetically connected and derive from a similar mantle source affected by different degrees of partial melting (and variable extent of crustal assimilation) producing either the CA magmas of Lipari (higher degrees) or the HKCA to SHO magmas of Vulcano (lower degrees). On a regional scale, the most primitive rocks (SiO2<56%, MgO>3.5%) of Lipari, Vulcano, Salina and Filicudi are suggested to derive from a similar MORB-like source, variably metasomatized by aqueous fluids coming from the slab and subordinately by the additions of sediments.

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.

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.

The influence of heat treatment on the phase relations in mineral growth systems

ABSTRACT Corundum is one of the most famous gems materials. Different heat treatment methods for enhancement purposes are commonly applied and accepted in the gem market. With this reason, the identification of the natural, unheated corundum is intensively investigated. In this study, aluminium hydroxide minerals and zircon are focused to observe the crystallization and phase change of these minerals during heat treatment procedures. Aluminium hydroxide minerals can be transformed to alumina with the corundum structure by heating. The reaction history of aluminium hydroxide minerals containing corundum was investigated comparing it with diaspore, boehmite, gibbsite and bayerite by TG and DTA methods. These hydroxide minerals were entirely transformed to corundum after heating at 600°C. Zircon inclusions in corundums from Ilakaka, Madagascar, were investigated for the influence of different heat-treatment temperatures on the recovery of their crystalline structure and on possible reactions within and with the host crystals. The host corundum was heated at 500, 800, 1000, 1200, 1400, 1600 and 1800°C. The crystallinity, the trapped pressure, and the decomposition of the zircon inclusions within the host corundum have been investigated by Raman spectroscopy. Radiation-damaged zircon inclusions may be used as an indicator for unheated Ilakaka corundum crystals. They are fully recrystallized after heating at 1000°C influencing the lowering of the 3 Raman band shift, the decreasing of FWHM of the 3 Raman band and the decreasing of the trapped pressure between the inclusion and the host corundum. Under microscopic observation, surface alterations of the inclusions can be firstly seen from transparent into frosted-like appearance at 1400°C. Then, between 1600°C and 1800 °C, the inclusion becomes partly or even completely molten. The decomposition of the zircon inclusion to m-ZrO2 and SiO2-glass phases begins at the rim of the inclusion after heating from 1200°C to 1600°C which can be detected by the surface change, the increase of the 3 Raman band position and the trapped pressure. At 1800°C, the zircon inclusions entirely melt transforming to solid phases during cooling like m-ZrO2 and SiO2-glass accompanied by an increase of pressure between the transformed inclusion and its host.

Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ≃ 147◦, while the CP yields to aSiOSi angle centered around 135◦. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of α-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

Effetto dell'addizione di fibre corte di carbonio su microstruttura e proprietà di materiali a base di ZrB2

Il lavoro di tesi, svolto presso l’Istituto di Scienza e Tecnologia dei Materiali Ceramici (ISTEC-CNR, Faenza, RA), ha affrontato la produzione e la caratterizzazione di ceramici a base di boruro di zirconio (ZrB2) con lo scopo di valutare l’efficacia delle fibre corte di carbonio come potenziale rinforzo. Il boruro di zirconio appartiene a una famiglia di materiali noti come UHTC (Ultra-High Temperature Ceramics) caratterizzati da elevato punto di fusione e in grado di mantenere la resistenza meccanica e operare con limitata ossidazione a temperature superiori ai 2000°C. Il principale ostacolo nella produzione dei materiali a base di ZrB2 è il processo di sintesi, infatti, a causa della loro elevata temperatura di fusione, per ottenere un materiale completamente denso è necessario utilizzare processi a temperatura e pressione elevati (T > 2000°C e P > 30 MPa), condizioni che vanno ad influenzare la microstruttura della matrice e delle fibre e di conseguenza le proprietà meccaniche del materiale. L’aggiunta di additivi di sinterizzazione idonei permette di ottenere materiali perfettamente densi anche a temperature e pressioni inferiori. Tuttavia lo ZrB2 non viene ampiamente utilizzato per applicazioni strutturali a causa della sua fragilità, per far fronte alla sua bassa tenacità il materiale viene spesso rinforzato con una fase allungata (whiskers o fibre). È già oggetto di studi l’utilizzo di fibre corte e whiskers di SiC per tenacizzare lo ZrB2, tuttavia la forte interfaccia che viene a crearsi tra fibra e matrice, che non permette il pull-out delle fibre, ci porta a credere che una fibra che non tenda a reagire con la matrice, presentando un’interfaccia più debole, possa portare ad una tenacizzazione più efficace. Per questo scopo sono stati realizzati mediante pressatura a caldo due materiali rinforzati con fibre corte di carbonio: ZrB2 + 5% vol MoSi2 + 8% vol fibre di carbonio e [ZrB2 + 2 % peso C] + 8% vol fibre di carbonio, indicati rispettivamente con Z5M_Cf e Z2C_Cf. Sono stati analizzati e discussi diversi aspetti del materiale rinforzato tra cui: il comportamento di densificazione durante la pressatura a caldo, l’evoluzione della microstruttura della matrice, la distribuzione e la morfologia delle fibre, l’influenza del rinforzo sulle proprietà meccaniche di durezza e tenacità e sulla resistenza all’ossidazione. L’elaborato è strutturato come segue: inizialmente sono state introdotte le caratteristiche generali dei ceramici avanzati tra cui le proprietà, la produzione e le applicazioni; successivamente è stata approfondita la descrizione dei materiali a base di boruro di zirconio, in particolare i processi produttivi e l’influenza degli additivi di sinterizzazione sulla densificazione e sulle proprietà; ci si è poi concentrati sull’effetto di una seconda fase allungata per il rinforzo del composito. Per quanto riguarda la parte sperimentale vengono descritte le principali fasi della preparazione e caratterizzazione dei materiali: le materie prime, disperse in un solvente, sono state miscelate mediante ball-milling, successivamente è stato evaporato il solvente e la polvere ottenuta è stata formata mediante pressatura uniassiale. I campioni, dopo essere stati sinterizzati mediante pressatura uniassiale a caldo, sono stati tagliati e lucidati a specchio per poter osservare la microstruttura. Quest’ultima è stata analizzata al SEM per studiare l’effetto dell’additivo di sinterizzazione (MoSi2 e carbonio) e l’interfaccia tra matrice e fase rinforzante. Per approfondire l’effetto del rinforzo sulle proprietà meccaniche sono state misurate la durezza e la tenacità del composito; infine è stata valutata la resistenza all’ossidazione mediante prove in aria a 1200°C e 1500°C. L’addizione di MoSi2 ha favorito la densificazione a 1800°C mediante formazione di una fase liquida transiente, tuttavia il materiale è caratterizzato da una porosità residua di ~ 7% vol. L’addizione del carbonio ha favorito la densificazione completa a 1900°C grazie alla reazione dall’additivo con gli ossidi superficiali dello ZrB2. La microstruttura delle matrici è piuttosto fine, con una dimensione media dei grani di ~ 2 μm per entrambi i materiali. Nel caso del materiale con Z5M_Cf sono presenti nella matrice particelle di SiC e fasi MoB derivanti dalla reazione dell’additivo con le fibre e con la matrice; invece nel materiale Z2C_Cf sono presenti grani di carbonio allungati tra i bordi grano, residui delle reazioni di densificazione. In entrambi i materiali le fibre sono distribuite omogeneamente e la loro interfaccia con la matrice è fortemente reattiva. Nel caso del materiale Z5M_Cf si è formata una struttura core-shell con lo strato più esterno formato da SiC, formato dalla reazione tra il siliciuro e la fibra di C. Nel caso del materiale Z2C_Cf non si forma una vera e propria interfaccia, ma la fibra risulta fortemente consumata per via dell’alta temperatura di sinterizzazione. I valori di durezza Vickers dei materiali Z5M_Cf e Z2C_Cf sono rispettivamente 11 GPa e 14 GPa, valori inferiori rispetto al valore di riferimento di 23 GPa dello ZrB2, ma giustificati dalla presenza di una fase meno dura: le fibre di carbonio e, nel caso di Z5M_Cf, anche della porosità residua. I valori di tenacità dei materiali Z5M_Cf e Z2C_Cf, misurati con il metodo dell’indentazione, sono rispettivamente 3.06 MPa·m0.5 e 3.19 MPa·m0.5. L’osservazione, per entrambi i materiali, del fenomeno di pull-out della fibra, sulla superficie di frattura, e della deviazione del percorso della cricca, all’interno della fibra di carbonio, lasciano supporre che siano attivi questi meccanismi tenacizzanti a contributo positivo, unitamente al contributo negativo legato allo stress residuo. La resistenza all’ossidazione dei due materiali è confrontabile a 1200°C, mentre dopo esposizione a 1500°C il materiale Z5M_Cf risulta più resistente rispetto al materiale Z2C_Cf grazie alla formazione di uno strato di SiO2 protettivo, che inibisce la diffusione dell’ossigeno all’interno della matrice. Successivamente, sono stati considerati metodi per migliorare la densità finale del materiale e abbassare ulteriormente la temperatura di sinterizzazione in modo da minimizzare la degenerazione della fibra. Da ricerca bibliografica è stato identificato il siliciuro di tantalio (TaSi2) come potenziale candidato. Pertanto è stato prodotto un terzo materiale a base di ZrB2 + Cf contenente una maggiore quantità di siliciuro (10% vol TaSi2) che ha portato ad una densità relativa del 96% a 1750°C. Questo studio ha permesso di approcciare per la prima volta le problematiche legate all’introduzione delle fibre di carbonio nella matrice di ZrB2. Investigazioni future saranno mirate alla termodinamica delle reazioni che hanno luogo in sinterizzazione per poter analizzare in maniera più sistematica la reattività delle fibre nei confronti della matrice e degli additivi. Inoltre riuscendo ad ottenere un materiale completamente denso e con fibre di carbonio poco reagite si potrà valutare la reale efficacia delle fibre di carbonio come possibili fasi tenacizzanti.

In-situ DRIFT/MS study of the conversion of ethanol on mixed Mg/Si oxides

Due to the high price of natural oil and harmful effects of its usage, as the increase in emission of greenhouse gases, the industry focused in searching of sustainable types of the raw materials for production of chemicals. Ethanol, produced by fermentation of sugars, is one of the more interesting renewable materials for chemical manufacturing. There are numerous applications for the conversion of ethanol into commodity chemicals. In particular, the production of 1,3-butadiene whose primary source is ethanol using multifunctional catalysts is attractive. With the 25% of world rubber manufacturers utilizing 1,3-butadiene, there is an exigent need for its sustainable production. In this research, the conversion of ethanol in one-step process to 1,3-butadiene was studied. According to the literature, the mechanisms which were proposed to explain the way ethanol transforms into butadiene require to have both acid and basic sites. But still, there are a lot of debate on this topic. Thus, the aim of this research work is a better understanding of the reaction pathways with all the possible intermediates and products which lead to the formation of butadiene from ethanol. The particular interests represent the catalysts, based on different ratio Mg/Si in comparison to bare magnesia and silica oxides, in order to identify a good combination of acid/basic sites for the adsorption and conversion of ethanol. Usage of spectroscopictechniques are important to extract information that could be helpful for understanding the processes on the molecular level. The diffuse reflectance infrared spectroscopy coupled to mass spectrometry (DRIFT-MS) was used to study the surface composition of the catalysts during the adsorption of ethanol and its transformation during the temperature program. Whereas, mass spectrometry was used to monitor the desorbed products. The set of studied materials include MgO, Mg/Si=0.1, Mg/Si=2, Mg/Si=3, Mg/Si=9 and SiO2 which were also characterized by means of surface area measurements.

Proton-conducting copolymers, blends and composites with phosphonic acid as protogenic group

ABSTRACT: In this work, proton conducting copolymers, polymer blends and composites containing phosphonic acid groups have been prepared. Proton conduction mechanisms in these materials are discussed respectively in both, the anhydrous and humidified state. Atom transfer radical copolymerization (ATRCP) of diisopropyl-p-vinylbenzyl phosphonate (DIPVBP) and 4-vinyl pyridine (4VP) is studied for the first time in this work. The kinetic parameters are obtained by using the 1H-NMR online technique. Proton conduction in poly(vinylbenzyl phosphonic acid) (PVBPA) homopolymer and its statistical copolymers with 4-vinyl pyridine (poly(VBPA-stat-4VP)s) are comprehensively studied in both, the “dry” and “wet” state. Effects of temperature, water content and polymer composition on proton conductivities are studied and proton transport mechanisms under various conditions are discussed. The proton conductivity of the polymers is in the range of 10-6-10-8 S/cm in nominally dry state at 150 oC. However, proton conductivity of the polymers increases rapidly with water content in the polymers which can reach 10-2 S/cm at the water uptake of 25% in the polymers. The highest proton conductivity obtained from the polymers can even reach 0.3 S/cm which was measured at 85oC with 80% relative humidity in the measuring atmosphere. Poly(4-vinyl pyridine) was grafted from the surface of SiO2 nanoparticles using ATRP in this work for the first time. Following this approach, silica nanoparticles with a shell of polymeric layer are used as basic particles in a polymeric acidic matrix. The proton conductivities of the composites are studied under both, humidified and dry conditions. In dry state, the conductivity of the composites is in the range of 10-10~10-4 S/cm at 150 oC. While in humid state, the composites show much higher proton conductivity. The highest proton conductivity obtained with the composites is 0.5 S/cm measured at 85oC with 80% relative humidity in the measuring atmosphere. The miscibility of poly (vinyl phosphonic acid) and PEO is studied for the first time in this work and a phase diagram is plotted based on a DSC study and optical microscopy. With this knowledge, homogeneous PVPA/PEO mixtures are prepared as proton-conducting polymer blends. The mobility of phosphonic acid groups and PEO in the blends is determined by 1H-MAS-NMR in temperature dependent measurements. The effect of composition and the role of PEO on proton conduction are discussed.

Nicht geträgerte Radioarsenisotope: ihre Herstellung, Abtrennung und Markierung von Proteinen

Das Element Arsen besitzt eine Reihe von Isotopen, die in nahezu trägerfreier Form (nca) produziert werden können und deshalb in der Radiopharmazie für die Diagnose oder Endoradiotherapie Verwendung finden können. Bei der Positronenemissionstomographie (PET) gibt es eine gewisse Lücke bei der Versorgung mit langlebigen Positronenemittern, die zur Untersuchung von langsamen physiologischen Prozessen wie z.B. der Biodistribution und Anreicherung von Antikörpern in Tumorgewebe eingesetzt werden können. Die beiden Arsenisotope 72As (T1/2 = 26 h, 88 % beta+) und 74As (T1/2 = 17,8 d, 29 % beta+) vereinen eine lange physikalische Halbwertszeit mit einer hohen Positronenemissionsrate und sind daher geeignete Kandidaten. Da das Verhalten von radioaktivem Arsen und seine Verwendung in der molekularen Bildgebung international relativ wenig bearbeitet sind, wurde die Radiochemie des Arsens von der Isotopenproduktion an Kernreaktor und Zyklotron, über die Entwicklung von Abtrennungsmethoden für Germanium und Arsen, bis hin zur Entwicklung einer soliden Markierungschemie für Antikörper weiterentwickelt. Die in dieser Arbeit bearbeiteten Felder sind: 1. Die Isotopenproduktion der relevanten Arsenisotope (72/74/77As) wurde an Kernreaktor und Zyklotron durch Bestrahlung von GeO2- und Germaniummetalltargets durchgeführt. Pro 6 h Bestrahlung von 100 mg Germanium konnten ca. 2 MBq 77As am TRIGA Reaktor in Mainz hergestellt werden. Am Zyklotron des DKFZ in Heidelberg konnten unter optimierten Bedingungen bei der Bestrahlung von Germaniummetall (EP = 15 Mev, 20 µA, 200 µAh) ca. 4 GBq 72As und ca. 400 MBq 74As produziert werden. 2. Die Entwicklung neuer Abtrennungsmethoden für nca 72/74/77As von makroskopischen Mengen Germanium wurde vorangetrieben. Für die Aufarbeitung von GeO2- und Germaniummetalltargets kamen insgesamt 8 verschiedene Methoden wie Festphasenextraktion, Flüssig-Flüssig-Extraktion, Destillation, Anionenaustauschchromatographie zum Einsatz. Die erzielten Ausbeuten lagen dabei zwischen 31 und 56 %. Es wurden Abtrennungsfaktoren des Germaniums zwischen 1000 und 1•10E6 erreicht. Alle erfolgreichen Abtrennungsmethoden lieferten *As(III) in 500 µl PBS-Puffer bei pH 7. Diese Form des Radioarsens ist für die Markierung von SH-modifizierten Molekülen, wie z.B. Antikörpern geeignet. 3. Die Entwicklung von Methoden zur Bestimmung des Oxidationszustandes von nca *As in organischem, neutralem wässrigen, oder stark sauren Medium mittels Radio-DC und Anionenaustauschchromatographie wurde durchgeführt und führte zu einem besseren Verständnis der Redoxchemie des nca *As. 4. SH-modifizierte Antikörper wurden mit 72/74/77As(III) markiert. Dabei wurden zwei Methoden (Modifizierung mit SATA und TCEP) miteinander verglichen. Während das *As(III) bei Verwendung von TCEP in Ausbeuten > 90 % mit dem Antikörper reagierte, wurde für SATA-modifizierte Antikörper in Abhängigkeit von der verwendeten Abtrennungsmethode eine breite Spanne von 0 % bis > 90 % beobachtet. 5. Es wurden Phantommessungen mit 18F, 72As und 74As am µ-PET-Scanner durchgeführt, um erste Aussagen über die zu erwartende Auflösung der Arsenisotope zu erhalten. Die Auflösung von 74As ist mit 18F vergleichbar, während die von 72As erkennbar schlechter ist.