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.

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.

High frequency acoustics in colloid-based meso- and nanostructures by spontaneous Brillouin light scattering

Materials that can mold the flow of elastic waves of certain energy in certain directions are called phononic materials. The present thesis deals essentially with such phononic systems, which are structured in the mesoscale (<1 µm), and with their individual components. Such systems show interesting phononic properties in the hypersonic region, i.e., at frequencies in the GHz range. It is shown that colloidal systems are excellent model systems for the realization of such phononic materials. Therefore, different structures and particle architectures are investigated by Brillouin light scattering, the inelastic scattering of light by phonons.rnThe experimental part of this work is divided into three chapters: Chapter 4 is concerned with the localized mechanical waves in the individual spherical colloidal particles, i.e., with their resonance- or eigenvibrations. The investigation of these vibrations with regard to the environment of the particles, their chemical composition, and the influence of temperature on nanoscopically structured colloids allows novel insights into the physical properties of colloids at small length scales. Furthermore, some general questions concerning light scattering on such systems, in dispute so far, are convincingly addressed.rnChapter 5 is a study of the traveling of mechanical waves in colloidal systems, consisting of ordered and disordered colloids in liquid or elastic matrix. Such systems show acoustic band gaps, which can be explained geometrically (Bragg gap) or by the interaction of the acoustic band with the eigenvibrations of the individual spheres (hybridization gap).rnWhile the latter has no analogue in photonics, the presence of strong phonon scatterers, when a large elastic mismatch between the composite components exists, can largely impact phonon propagation in analogy to strong multiple light scattering systems. The former is exemplified in silica based phononic structures that opens the door to new ways of sound propagation manipulation.rnChapter 6 describes the first measurement of the elastic moduli in newly fabricated by physical vapor deposition so-called ‘stable organic glasses’. rnIn brief, this thesis explores novel phenomena in colloid-based hypersonic phononic structures, utilizing a versatile microfabrication technique along with different colloid architectures provided by material science, and applying a non-destructive optical experimental tool to record dispersion diagrams.rn

Solid-state NMR investigation of phosphonic acid based proton conducting materials

In this work, new promising proton conducting fuel cell membrane materials were characterized in terms of their structure and dynamic properties using solid-state nuclear magnetic resonance (NMR) spectroscopy and X-ray diffraction. Structurally different, phosphonic acid (PA) containing materials were systematically evaluated for possible high-temperature operation (e.g. at T>100°C). Notably, 1H, 2H and 31P magic angle spinning (MAS) NMR provided insight into local connectivities and dynamics of the hydrogen bonded network, while packing arrangements were identified by means of heteronuclear dipolar recoupling techniques.rnThe first part of this work introduced rather crystalline, low molecular weight ionomers for proton conducting membranes, where six different geometries such as line, triangle, screw, tetrahedron, square and hexagon, were investigated. The hexagon was identified as the most promising geometry with high-temperature bulk proton conductivities in the range of 10-3 Scm-1 at a relative humidity of 50%. However, 2H NMR and TGA-MS data suggest that the bulk proton transport is mainly due to the presence of crystal water. Single crystal X-ray data revealed that in the tetrahedron phosphonic acids form tetrameric clusters isolating the mobile protons while the phosphonic acids in the hexagon form zigzag-type pathways through the sample.rnThe second part of this work demonstrates how acid-base pairing and the choice of appropriate spacers may influence proton conduction. Different ratios of statistical copolymers of poly (vinylphosphonic acid) and poly (4-vinylpyridine) were measured to derive information about the local structure and chemical changes. Though anhydrous proton conductivities of all statistical copolymers are rather poor, the conductivity increases to 10-2 S cm-1 when exposing the sample to relative humidity of 80%. In contrast to PVPA, anhydride formation of phosphonic acids in the copolymer is not reversible even when exposing the sample to a relative humidity of 100%.rnIn addition, the influence of both spacers and degree of backbone crystallinity on bulk proton conductivity was investigated. Unlike in systems such as poly benzimidazole (PBI), spacers were inserted between the protogenic groups along the backbone. It was found that dilution of the protogenic groups decreases the conductivity, but compared to PVPA, similar apparent activation energies for local motions were obtained from both variable temperature 1H NMR and impedance spectroscopy data. These observations suggest the formation of phosphonic acid clusters with high degrees of local proton motion, where only a fraction of motions contribute to the observable bulk proton conductivity. Additionally, it was shown that gradual changes of the spacer length lead to different morphologies.rnIn summary, applying advanced solid-state NMR and X-ray analysis, structural and dynamic phenomena in proton conducting materials were identified on a molecular level. The results were discussed with respect to different proton conduction mechanisms and may contribute to a more rational design or improvement of proton conducting membranes.rn

Electron spectroscopy of novel charge transfer systems based on polycyclic aromatic hydrocarbons

Organische Ladungstransfersysteme weisen eine Vielfalt von konkurrierenden Wechselwirkungen zwischen Ladungs-, Spin- und Gitterfreiheitsgraden auf. Dies führt zu interessanten physikalischen Eigenschaften, wie metallische Leitfähigkeit, Supraleitung und Magnetismus. Diese Dissertation beschäftigt sich mit der elektronischen Struktur von organischen Ladungstransfersalzen aus drei Material-Familien. Dabei kamen unterschiedliche Photoemissions- und Röntgenspektroskopietechniken zum Einsatz. Die untersuchten Moleküle wurden z.T. im MPI für Polymerforschung synthetisiert. Sie stammen aus der Familie der Coronene (Donor Hexamethoxycoronen HMC und Akzeptor Coronen-hexaon COHON) und Pyrene (Donor Tetra- und Hexamethoxypyren TMP und HMP) im Komplex mit dem klassischen starken Akzeptor Tetracyanoquinodimethan (TCNQ). Als dritte Familie wurden Ladungstransfersalze der k-(BEDT-TTF)2X Familie (X ist ein monovalentes Anion) untersucht. Diese Materialien liegen nahe bei einem Bandbreite-kontrollierten Mottübergang im Phasendiagramm.rnFür Untersuchungen mittels Ultraviolett-Photoelektronenspektroskopie (UPS) wurden UHV-deponierte dünne Filme erzeugt. Dabei kam ein neuer Doppelverdampfer zum Einsatz, welcher speziell für Milligramm-Materialmengen entwickelt wurde. Diese Methode wies im Ladungstransferkomplex im Vergleich mit der reinen Donor- und Akzeptorspezies energetische Verschiebungen von Valenzzuständen im Bereich weniger 100meV nach. Ein wichtiger Aspekt der UPS-Messungen lag im direkten Vergleich mit ab-initio Rechnungen.rnDas Problem der unvermeidbaren Oberflächenverunreinigungen von lösungsgezüchteten 3D-Kristallen wurde durch die Methode Hard-X-ray Photoelectron Spectroscopy (HAXPES) bei Photonenenergien um 6 keV (am Elektronenspeicherring PETRA III in Hamburg) überwunden. Die große mittlere freie Weglänge der Photoelektronen im Bereich von 15 nm resultiert in echter Volumensensitivität. Die ersten HAXPES Experimente an Ladungstransferkomplexen weltweit zeigten große chemische Verschiebungen (mehrere eV). In der Verbindung HMPx-TCNQy ist die N1s-Linie ein Fingerabdruck der Cyanogruppe im TCNQ und zeigt eine Aufspaltung und einen Shift zu höheren Bindungsenergien von bis zu 6 eV mit zunehmendem HMP-Gehalt. Umgekehrt ist die O1s-Linie ein Fingerabdruck der Methoxygruppe in HMP und zeigt eine markante Aufspaltung und eine Verschiebung zu geringeren Bindungsenergien (bis zu etwa 2,5eV chemischer Verschiebung), d.h. eine Größenordnung größer als die im Valenzbereich.rnAls weitere synchrotronstrahlungsbasierte Technik wurde Near-Edge-X-ray-Absorption Fine Structure (NEXAFS) Spektroskopie am Speicherring ANKA Karlsruhe intensiv genutzt. Die mittlere freie Weglänge der niederenergetischen Sekundärelektronen (um 5 nm). Starke Intensitätsvariationen von bestimmten Vorkanten-Resonanzen (als Signatur der unbesetzte Zustandsdichte) zeigen unmittelbar die Änderung der Besetzungszahlen der beteiligten Orbitale in der unmittelbaren Umgebung des angeregten Atoms. Damit war es möglich, präzise die Beteiligung spezifischer Orbitale im Ladungstransfermechanismus nachzuweisen. Im genannten Komplex wird Ladung von den Methoxy-Orbitalen 2e(Pi*) und 6a1(σ*) zu den Cyano-Orbitalen b3g und au(Pi*) und – in geringerem Maße – zum b1g und b2u(σ*) der Cyanogruppe transferiert. Zusätzlich treten kleine energetische Shifts mit unterschiedlichem Vorzeichen für die Donor- und Akzeptor-Resonanzen auf, vergleichbar mit den in UPS beobachteten Shifts.rn

Femtosecond laser-based investigations of magnetic circular dichroism in near-threshold photoemission

During the last decades magnetic circular dichroism (MCD) has attracted much interest and evolved into various experimental methods for the investigation of magnetic thin films. For example, synchrotron-based X-ray magnetic circular dichroism (XMCD) displays the absolute values of spin and orbital magnetic moments. It thereby benefits from large asymmetry values of more than 30% due to the excitation of atomic core-levels. Similarly large values are also expected for threshold photoemission magnetic circular dichroism (TPMCD). Using lasers with photon energies in the range of the sample work function this method gives access to the occupied electronic structure close to the Fermi level. However, except for the case of Ni(001) there exist only few studies on TPMCD moreover revealing much smaller asymmetries than XMCD-measurements. Also the basic physical mechanisms of TPMCD are not satisfactorily understood. In this work we therefore investigate TPMCD in one- and two-photon photoemission (1PPE and 2PPE) for ferromagnetic Heusler alloys and ultrathin Co films using ultrashort pulsed laser light. The observed dichroism is explained by a non-conventional photoemission model using spin-resolved band-structure calculations and linear response theory. For the two Heusler alloys Ni2MnGa and Co2FeSi we give first evidence of TPMCD in the regime of two-photon photoemission. Systematic investigations concerning general properties of TPMCD in 1PPE and 2PPE are carried out at ultrathin Co films grown on Pt(111). Here, photon-energy dependent measurements reveal asymmetries of 1.9% in 1PPE and 11.7% in 2PPE. TPMCD measurements at decreased work function even yield larger asymmetries of 6.2% (1PPE) and 17% (2PPE), respectively. This demonstrates that enlarged asymmetries are also attainable for the TPMCD effect on Co(111). Furthermore, we find that the TPMCD asymmetry is bulk-sensitive for 1PPE and 2PPE. This means that the basic mechanism leading to the observed dichroism must be connected to Co bulk properties; surface effects do not play a crucial role. Finally, the enhanced TPMCD asymmetries in 2PPE compared to the 1PPE case are traced back to the dominant influence of the first excitation step and the existence of a real intermediate state. The observed TPMCD asymmetries cannot be interpreted by conventional photoemission theory which only considers direct interband transitions in the direction of observation (Γ-L). For Co(111), these transitions lead to evanescent final states. The excitation to such states, however, is incompatible with the measured bulk-sensitivity of the asymmetry. Therefore, we generalize this model by proposing the TPMCD signal to arise mostly from direct interband transitions in crystallographic directions other than (Γ-L). The necessary additional momentum transfer to the excited electrons is most probably provided by electron-phonon or -magnon scattering processes. Corresponding calculations on the basis of this model are in reasonable agreement with the experimental results so that this approach represents a promising tool for a quantitative description of the TPMCD effect. The present findings encourage an implementation of our experimental technique to time- and spatially-resolved photoemission electron microscopy, thereby enabling a real time imaging of magnetization dynamics of single excited states in a ferromagnetic material on a femtosecond timescale.

Hard x-ray photoelectron spectroscopy of bulk and thin films of Heusler compounds

X-ray photoemission spectroscopy (XPS) is one of the most universal and powerful tools for investigation of chemical states and electronic structures of materials. The application of hard x-rays increases the inelastic mean free path of the emitted electrons within the solid and thus makes hard x-ray photoelectron spectroscopy (HAXPES) a bulk sensitive probe for solid state research and especially a very effective nondestructive technique to study buried layers.rnThis thesis focuses on the investigation of multilayer structures, used in magnetic tunnel junctions (MTJs), by a number of techniques applying HAXPES. MTJs are the most important components of novel nanoscale devices employed in spintronics. rnThe investigation and deep understanding of the mechanisms responsible for the high performance of such devices and properties of employed magnetic materials that are, in turn, defined by their electronic structure becomes feasible applying HAXPES. Thus the process of B diffusion in CoFeB-based MTJs was investigated with respect to the annealing temperature and its influence on the changes in the electronic structure of CoFeB electrodes that clarify the behaviour and huge TMR ratio values obtained in such devices. These results are presented in chapter 6. The results of investigation of the changes in the valence states of buried off-stoichiometric Co2MnSi electrodes were investigated with respect to the Mn content α and its influence on the observed TMR ratio are described in chapter 7.rnrnMagnetoelectronic properties such as exchange splitting in ferromagnetic materials as well as the macroscopic magnetic ordering can be studied by magnetic circular dichroism in photoemission (MCDAD). It is characterized by the appearance of an asymmetry in the photoemission spectra taken either from the magnetized sample with the reversal of the photon helicity or by reversal of magnetization direction of the sample when the photon helicity direction is fixed. Though recently it has been widely applied for the characterization of surfaces using low energy photons, the bulk properties have stayed inaccessible. Therefore in this work this method was integrated to HAXPES to provide an access to exploration of magnetic phenomena in the buried layers of the complex multilayer structures. Chapter 8 contains the results of the MCDAD measurements employing hard x-rays for exploration of magnetic properties of the common CoFe-based band-ferromagnets as well as half-metallic ferromagnet Co2FeAl-based MTJs.rnrnInasmuch as the magnetoresistive characteristics in spintronic devices are fully defined by the electron spins of ferromagnetic materials their direct measurements always attracted much attention but up to date have been limited by the surface sensitivity of the developed techniques. Chapter 9 presents the results on the successfully performed spin-resolved HAXPES experiment using a spin polarimeter of the SPLEED-type on a buried Co2FeAl0.5Si0.5 magnetic layer. The measurements prove that a spin polarization of about 50 % is retained during the transmission of the photoelectrons emitted from the Fe 2p3/2 state through a 3-nm-thick oxide capping layer.rn

Spectroscopic investigations of functional bulk and thin film Heusler alloys

Functional and smart materials have gained large scientific and practical interest in current research and development. The Heusler alloys form an important class of functional materials used in spintronics, thermoelectrics, and for shape memory alloy applications. An important aspect of functional materials is the adaptability of their physical properties. In this work functional polycrystalline bulk and epitaxial thin film Heusler alloys are characterized by means of spectroscopic investigation methods, X-ray magnetic circular dichroism (XMCD) and energy dispersive X-ray analysis (EDX). With EDX the homogeneity of the samples is studied extensively. For some cases of quaternary compounds, for example Co2(MnxTi1−x)Sn and Co2(Mn0.5Dy0.5)Sn, an interesting phase separation in two nearly pure ternary Heusler phases occurs. For these samples the phase separation leads to an improvement of thermoelectric properties. XMCD as the main investigation method was used to study Co2TiZ (Z = Si, Sn, and Sb), Co2(MnxTi1−x)Si, Co2(MnxTi1−x)Ge, Co2Mn(Ga1−xGex), Co2FeAl, Mn2VAl, and Ni2MnGa Heusler compounds. The element-specific magnetic moments are calculated. Also, the spin-resolved unoccupied density of states is determined, for example giving hints for half-metallic ferromagnetism for some Co-based compounds. The systematic change of the magnetic moments and the shift of the Fermi energy is a proof that Heusler alloys are suitable for a controlled tailoring of physical properties. The comparison of the experimental results with theoretical predictions improves the understanding of complex materials needed to optimize functional Heusler alloys.

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.