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.

DNA methylation, histone modification and the transcription factor dE2F1 in Drosophila

Drosophila melanogaster enthält eine geringe Menge an 5-methyl-Cytosin. Die von mir untersuchte männliche Keimbahn von Drosophila weist jedoch keine nachweisbaren Mengen an DNA-Methylierung auf. Eine künstliche Expression der murinen de novo Methyltransferasen, DNMT3A und DNMT3B1, in den Fliegenhoden, führte nicht zu der erwarteten Methylierungszunahme und hatte keinen Effekt auf die Fruchtbarkeit der Männchen. Auch die gewebespezifische Expression unter der Verwendung des UAS/GAL4-Systems zeigte keine phenotypischen Veränderungen. Hingegen fanden wir auf Protein-Ebene des Chromatins von D. melanogaster und D. hydei spezifische Modifikationsmuster der Histone H3 und H4 in der Keimbahn, wie auch in den somatischen Zellen des Hodenschlauches. Die Modifikationsmuster der beiden Zelltypen unterscheiden sich grundlegend und weichen zudem von dem für Eu- und Heterochromatin erwarteten ab, was auf eine größere Komplexität des „Histon-Codes“ als angenommen hindeutet. Folglich liegt die epigenetische Information in Drosophila wahrscheinlich anstatt auf DNA- auf Protein-Ebene, wodurch Genexpression über die Chromatinstruktur reguliert wird. Es wurde gezeigt, dass der Transkriptionsfaktor E2F, der eine Schlüsselfunktion im Zellzyklus hat, durch unterschiedliche Transkripte offenbar quantitativ reguliert wird. Unsere Nachforschungen ergaben, dass die drei E2F1 Genprodukte in Drosophila neben ihrer Zellspezifität auch in unterschiedlichen Expressionsniveaus auftreten, was die Annahme einer quantitativen Expression unterstützt. Die verschiedenen Funktionen der multiplen Gene in Säugern, könnten so funktionell kompensiert werden. Die durch die Expression dreier dE2F1-Transkripte vermutete Synthese verschiedener Proteine konnte nicht bewiesen werden.

Unconventional correlated systems

This work comprises three different types of unconventional correlated systems.rnChapters 3-5 of this work are about the open shell compounds Rb4O6 and Cs4O6. These mixed valent compounds contain oxygen in two different modifications: the closed-shell peroxide anion is nonmagnetic, whereas the hyperoxide anion contains an unpaired electrons in an antibonding pi*-orbital. Due to this electron magnetic ordering is rendered possible. In contrast to theoretical predictions, which suggested half-metallic ferromagnetism for Rb4O6,rndominating antiferromagnetic interactions were found in the experiment. Besidesrna symmetry reduction due to the mixed valency, strong electronic correlations of this highly molecular system determine its properties; it is a magnetically frustrated insulator. The corresponding Cs4O6 was found to show similar properties.rnChapters 6-9 of this work are about intermetallic Heusler superconductors. rnAll of these superconductors were rationally designed using the van Hove scenario as a working recipe. A saddle point in the energy dispersion curve of a solid leads to a van Hove singularity in the density of states. In the Ni-based and Pd-based Heusler superconductors presented in this work this sort of a valence instability occurs at the high-symmetry L point and coincides or nearly coincides with the Fermi level. The compounds escape the high density of states at the Fermi energy through a transition into the correlated superconducting state.rnChapter 10 of this work is about the tetragonally distorted ferrimagnetic DO22 phase of Mn3Ga. This hard-magnetic modification is technologically useful for spin torque transfer applications. The phase exhibits two different crystallographic sites that are occupied by Mn atoms and can thus be written as Mn2MnGa. The competition between the mainly itinerant moments of the Mn atoms at the Wyckoff position 4d and the localized moments of the Mn atoms at the Wyckoff position 2b leads to magnetic correlations. The antiferromagnetic orientation of these moments determines the compound to exhibit a resulting magnetic moment of approximately 1 µB per formula unit in a partially compensated ferrimagnetic configuration.

Structure property relations in complex copolymer systems

A thorough investigation was made of the structure-property relation of well-defined statistical, gradient and block copolymers of various compositions. Among the copolymers studied were those which were synthesized using isobornyl acrylate (IBA) and n-butyl acrylate (nBA) monomer units. The copolymers exhibited several unique properties that make them suitable materials for a range of applications. The thermomechanical properties of these new materials were compared to acrylate homopolymers. By the proper choice of the IBA/nBA monomer ratio, it was possible to tune the glass transition temperature of the statistical P(IBA-co-nBA) copolymers. The measured Tg’s of the copolymers with different IBA/nBA monomer ratios followed a trend that fitted well with the Fox equation prediction. While statistical copolymers showed a single glass transition (Tg between -50 and 90 ºC depending on composition), DSC block copolymers showed two Tg’s and the gradient copolymer showed a single, but very broad, glass transition. PMBL-PBA-PMBL triblock copolymers of different composition ratios were also studied and revealed a microphase separated morphology of mostly cylindrical PMBL domains hexagonally arranged in the PBA matrix. DMA studies confirmed the phase separated morphology of the copolymers. Tensile studies showed the linear PMBL-PBA-PMBL triblock copolymers having a relatively low elongation at break that was increased by replacing the PMBL hard blocks with the less brittle random PMBL-r-PMMA blocks. The 10- and 20-arm PBA-PMBL copolymers which were studied revealed even more unique properties. SAXS results showed a mixture of cylindrical PMBL domains hexagonally arranged in the PBA matrix, as well as lamellar. Despite PMBL’s brittleness, the triblock and multi-arm PBA-PMBL copolymers could become suitable materials for high temperature applications due to PMBL’s high glass transition temperature and high thermal stability. The structure-property relation of multi-arm star PBA-PMMA block copolymers was also investigated. Small-angle X-ray scattering revealed a phase separated morphology of cylindrical PMMA domains hexagonally arranged in the PBA matrix. DMA studies found that these materials possess typical elastomeric behavior in a broad range of service temperatures up to at least 250°C. The ultimate tensile strength and the elastic modulus of the 10- and 20-arm star PBA-PMMA block copolymers are significantly higher than those of their 3-arm or linear ABA type counterparts with similar composition, indicating a strong effect of the number of arms on the tensile properties. Siloxane-based copolymers were also studied and one of the main objectives here was to examine the possibility to synthesize trifluoropropyl-containing siloxane copolymers of gradient distribution of trifluoropropyl groups along the chain. DMA results of the PDMS-PMTFPS siloxane copolymers synthesized via simultaneous copolymerization showed that due to the large difference in reactivity rates of 2,4,6-tris(3,3,3-trifluoropropyl)-2,4,6-trimethylcyclotrisiloxane (F) and hexamethylcyclotrisiloxane (D), a copolymer of almost block structure containing only a narrow intermediate fragment with gradient distribution of the component units was obtained. A more dispersed distribution of the trifluoropropyl groups was obtained by the semi-batch copolymerization process, as the DMA results revealed more ‘‘pure gradient type’’ features for the siloxane copolymers which were synthesized by adding F at a controlled rate to the polymerization of the less reactive D. As with trifluoropropyl-containing siloxane copolymers, vinyl-containing polysiloxanes may be converted to a variety of useful polysiloxane materials by chemical modification. But much like the trifluoropropyl-containing siloxane copolymers, as a result of so much difference in the reactivities between the component units 2,4,6-trivinyl-2,4,6-trimethylcyclotrisiloxane (V) and hexamethylcyclotrisiloxane (D), thermal and mechanical properties of the PDMS-PMVS copolymers obtained by simultaneous copolymerization was similar to those of block copolymers. Only the copolymers obtained by semi-batch method showed properties typical for gradient copolymers.

Mechanical properties and structure of gel systems

Understanding the origins of the mechanical properties and its correlation withrnthe microstructure of gel systems is of great scientific and industrial interest. Inrngeneral, colloidal gels can be classified into chemical and physical gels, accordingrnto the life time of the network bonds. The characteristic di↵erences in gelationrndynamics can be observed with rheological measurements.rnAs a model system, a mixture of sodium silicate and low concentration sulfuric acidrnwas used. Nano-sized silica particles grow and aggregate to a system-spanning gelrnnetwork. The influence of the finite solubility of silica at high pH on the gelationrnwas studied with classical and piezo rheometer. The storage modulus of therngel grew logarithmically with time with two distinct growth laws. A relaxationrnat low frequency was observed in the frequency dependent measurements. I attributernthese two behaviors as a sign of structural rearrangements due to the finiternsolubility of silica at high pH. The reaction equilibrium between formation andrndissolution of bonds leads to a finite life time of the bonds and behavior similar tornphysical gel. The frequency dependence was more pronounced for lower water concentrations,rnhigher temperatures and shorter reaction times. With two relaxationrnmodels, I deduced characteristic relaxation times from the experimental data. Besidesrnrheology, the evolution of silica gels at high pH on di↵erent length scales wasrnstudied by NMR and dynamic light scattering. The results revealed that the primaryrnparticles existed already in sodium silicate and aggregated after the mixingrnof reactants due to a chemical reaction. Throughout the aggregation process thernsystem was in its chemical reaction equilibrium. Applying large oscillatory shearrnstrain to the gel allowed for modifying the gel modulus. The e↵ect of shear andrnshear history on the rheological properties of the gel were investigated. The storagernmodulus of the final gel increased with increasing strain. This behavior can be explained with (i) shear-induced aggregate compaction and (ii) combination ofrnbreakage and new formation of bonds.rnIn comparison with the physical gel-like behavior of the silica gel at high pH, typicalrnchemical gel features were exhibited by other gels formed from various chemicalrnreactions. Influences of the chemical structure modification on the gelation wererninvestigated with the piezo-rheometer. The external stimuli can be applied to tunernthe mechanical properties of the gel systems.

Exact quantum Monte Carlo methods for correlated electron systems

Thema dieser Arbeit ist die Entwicklung und Kombination verschiedener numerischer Methoden, sowie deren Anwendung auf Probleme stark korrelierter Elektronensysteme. Solche Materialien zeigen viele interessante physikalische Eigenschaften, wie z.B. Supraleitung und magnetische Ordnung und spielen eine bedeutende Rolle in technischen Anwendungen. Es werden zwei verschiedene Modelle behandelt: das Hubbard-Modell und das Kondo-Gitter-Modell (KLM). In den letzten Jahrzehnten konnten bereits viele Erkenntnisse durch die numerische Lösung dieser Modelle gewonnen werden. Dennoch bleibt der physikalische Ursprung vieler Effekte verborgen. Grund dafür ist die Beschränkung aktueller Methoden auf bestimmte Parameterbereiche. Eine der stärksten Einschränkungen ist das Fehlen effizienter Algorithmen für tiefe Temperaturen.rnrnBasierend auf dem Blankenbecler-Scalapino-Sugar Quanten-Monte-Carlo (BSS-QMC) Algorithmus präsentieren wir eine numerisch exakte Methode, die das Hubbard-Modell und das KLM effizient bei sehr tiefen Temperaturen löst. Diese Methode wird auf den Mott-Übergang im zweidimensionalen Hubbard-Modell angewendet. Im Gegensatz zu früheren Studien können wir einen Mott-Übergang bei endlichen Temperaturen und endlichen Wechselwirkungen klar ausschließen.rnrnAuf der Basis dieses exakten BSS-QMC Algorithmus, haben wir einen Störstellenlöser für die dynamische Molekularfeld Theorie (DMFT) sowie ihre Cluster Erweiterungen (CDMFT) entwickelt. Die DMFT ist die vorherrschende Theorie stark korrelierter Systeme, bei denen übliche Bandstrukturrechnungen versagen. Eine Hauptlimitation ist dabei die Verfügbarkeit effizienter Störstellenlöser für das intrinsische Quantenproblem. Der in dieser Arbeit entwickelte Algorithmus hat das gleiche überlegene Skalierungsverhalten mit der inversen Temperatur wie BSS-QMC. Wir untersuchen den Mott-Übergang im Rahmen der DMFT und analysieren den Einfluss von systematischen Fehlern auf diesen Übergang.rnrnEin weiteres prominentes Thema ist die Vernachlässigung von nicht-lokalen Wechselwirkungen in der DMFT. Hierzu kombinieren wir direkte BSS-QMC Gitterrechnungen mit CDMFT für das halb gefüllte zweidimensionale anisotrope Hubbard Modell, das dotierte Hubbard Modell und das KLM. Die Ergebnisse für die verschiedenen Modelle unterscheiden sich stark: während nicht-lokale Korrelationen eine wichtige Rolle im zweidimensionalen (anisotropen) Modell spielen, ist in der paramagnetischen Phase die Impulsabhängigkeit der Selbstenergie für stark dotierte Systeme und für das KLM deutlich schwächer. Eine bemerkenswerte Erkenntnis ist, dass die Selbstenergie sich durch die nicht-wechselwirkende Dispersion parametrisieren lässt. Die spezielle Struktur der Selbstenergie im Impulsraum kann sehr nützlich für die Klassifizierung von elektronischen Korrelationseffekten sein und öffnet den Weg für die Entwicklung neuer Schemata über die Grenzen der DMFT hinaus.