Herstellung und physikochemische Charakterisierung von planaren gestützten Lipid-Modellmembran-Systemen

Es wurden funktionalisierte polymerunterstützte planare Phospholipid-Modellmembran-Systeme hergestellt und auf jeder Präparationsstufe eingehend charakterisiert. Dünne Polysaccharidfilme wurden in der Form von quellbaren Gelen auf oxidische Oberflächen aufgebracht und bezüglich ihres Quellungsverhaltens und der Oberflächeneigenschaften in Abhängigkeit vom Wassergehalt untersucht. Lipidmonoschichten unterschiedlicher Zusammensetzung wurden mittels Langmuir-Blodgett-Tranfer auf Polymersubstrate übertragen und bezüglich der Stärke der Lipid/Polymer Wechselwirkung, der lateralen Selbstdiffusion in Abhängigkeit von der Wasseraktivität, dem Spreitverhalten der monomolekularen Membran auf dem Substrat in Abhängigkeit von der Wasseraktivität und dem Lateraldruck der Monoschicht, sowie des Ausmaßes der Hydratation im Kopfgruppenbereich der Lipidmembran in Abhängigkeit von der Wasseraktivität mittels Fluoreszensondenmethoden (Fluoreszenzerholung nach Photobleichung (FRAP), Fluoreszenzmikroskopie und Fluoreszenzspektroskopie) untersucht. Diffusions- und Spreitverhalten von amphiphilen Monoschichten auf Polymersubstraten wurden auf der Basis von in dieser Arbeit entwickelten physikalischen Modellen diskutiert. Mittels Langmuir-Schäfer Transfer wurde auf polymerunterstützte Lipidmonoschichten eine zweite Monoschicht übertragen. Die somit erhaltenen Lipid-Doppelschichtmembranen wurden bezüglich ihrer Stabilität, der lateralen Struktur, der lateralen Selbstdiffusion, des Spreitverhaltens auf unbedeckte Bereiche sowie der Stärke der Membran/Substrat Wechselwirkung vermittels Fluoreszenzmikroskopie, FRAP und Interferenz-Kontrast-Mikroskopie (RICM) untersucht. Schließlich wurden substratgestützte Doppelschicht-Lipidmembranen mit als Protonenpumpen fungierenden integralen Membranproteinen versehen. Die laterale Selbstdiffusion der rekonstituierten Proteinmoleküle wurde mittels FRAP, die funktionale Aktivität der Protonenpumpen mit einem Ionen-sensitiven Feldeffekttransistor-Array analysiert.

Untersuchung der Filmbildung aus Polymerdispersionen mit Hilfe der forcierten Rayleighstreuung

Der Filmbildungsprozess wässriger Polymerdispersionen wurde mit forcierter Rayleighstreuung untersucht.Diffusionskoeffizienten D von Farbstoffsonden sind damit zwischen 10^-21 und 10^-9 m^2s^-1 zugänglich. Das Diffusionsverhalten der hydrophoben FarbstoffsondeAberchrome in feuchten und trockenen Filmen aus Poly(n-butylmethacrylat-co-acrylsäure)-Dispersionen sowie bei der Wiederbefeuchtung trockener Filme wurde untersucht.Die Dynamik von Aberchrome weicht in feuchten Filmen von Fickscher Diffusion ab. Dies äußert sich in Abweichungen vomcharakteristischen q^2-Verhalten der Relaxationsrate tau^-1 (tau^-1 = Dq^2; q:Streuvektor) und im Auftreten gestrecktexponentieller Intensitätsrelaxationskurven. Diese Anomalie wurde mit einem apparenten, längenskalenabhängigen Diffusionskoeffizienten Dapp(Lambda) (Lambda= 2Pi/q) beschrieben, der für Lambda -> 0 Werte annimmt, die einem homogen hydroplastifizierten Polymermaterial entsprechen, während Dapp(Lambda) für Lambda -> Unendlich stark anwächst. Diese Anomalien verschwinden bei Wassergehalten entsprechend der Polymerwasserlöslichkeit. Weiteres Trocknen führt zum Absinken des Fickschen-Diffusionskoeffizienten auf einen Grenzwert für trockene Filme. Die Ergebnisse konnten mit einem Zwei-Zustands-Modell beschrieben werden: Die Sonde diffundiert Ficksch in einer hydrophilen Grenzflächenphase und einer langsamen Polymerphase. Austausch zwischendiesen Phasen ist ohne Einschränkung möglich. Das Modell erlaubt die Quantifizierung des Einflusses des Trocknungsprozesses auf Polymer- und Grenzflächenphaseneigenschaften.Dies wurde durch systematische Veränderungen der Grenzflächeneigenschaften demonstriert. Dies geschah durch Acrylsäuregehaltvariationen in den Poly(n-butylmethacrylat-co-acrylsäure)-Dispersionen und Emulgatorbedeckungsgradvariationen. In beiden Fällen reflektieren sich Dispersionsveränderungen in einer Veränderung der Modellparameter in der Grenzflächenphase.

Thermodynamics of aqueous biopolymer solutions and fractionation of Dextran

Flory-Huggins interaction parameters and thermal diffusion coefficients were measured for aqueous biopolymer solutions. Dextran (a water soluble polysaccharide) and bovine serum albumin (BSA, a water soluble protein) were used for this study. The former polymer is representative for chain macromolecules and the latter is for globular macromolecules. The interaction parameters for the systems water/dextran and water/BSA were determined as a function of composition by means of vapor pressure measurements, using a combination of headspace sampling and gas chromatography (HS-GC). A new theoretical approach, accounting for chain connectivity and conformational variability, describes the observed dependencies quantitatively for the system water/dextran and qualitatively for the system water/BSA. The phase diagrams of the ternary systems water/methanol/dextran and water/dextran/BSA were determined via cloud point measurements and modeled by means of the direct minimization of the Gibbs energy using the information on the binary subsystems as input parameters. The thermal diffusion of dextran was studied for aqueous solutions in the temperature range 15 < T < 55 oC. The effects of the addition of urea were also studied. In the absence of urea, the Soret coefficient ST changes its sign as T is varied; it is positive for T > 45.0 oC, but negative for T < 45.0 oC. The positive sign of ST means that the dextran molecules migrate towards the cold side of the fluid; this behavior is typical for polymer solutions. While a negative sign indicates the macromolecules move toward the hot side; this behavior has so far not been observed with any other binary aqueous polymer solutions. The addition of urea to the aqueous solution of dextran increases ST and reduces the inversion temperature. For 2 M urea, the change in the sign of ST is observed at T = 29.7 oC. At higher temperature ST is always positive in the studied temperature range. To rationalize these observations it is assumed that the addition of urea opens hydrogen bonds, similar to that induced by an increase in temperature. For a future extension of the thermodynamic studies to the effects of poly-dispersity, dextran was fractionated by means of a recently developed technique called Continuous Spin Fractionation (CSF). The solvent/precipitant/polymer system used for the thermodynamic studies served as the basis for the fractionation of dextran The starting polymer had a weight average molar mass Mw = 11.1 kg/mol and a molecular non-uniformity U= Mw / Mn -1= 1.0. Seventy grams of dextran were fractionated using water as the solvent and methanol as the precipitant. Five fractionation steps yielded four samples with Mw values between 4.36 and 18.2 kg/mol and U values ranging from 0.28 to 0.48.

Influence of spatial constraints in non-crystalline regions on the polymer dynamics in semi-crystalline polyethylene: a solid state NMR study

A broad variety of solid state NMR techniques were used to investigate the chain dynamics in several polyethylene (PE) samples, including ultrahigh molecular weight PEs (UHMW-PEs) and low molecular weight PEs (LMW-PEs). Via changing the processing history, i.e. melt/solution crystallization and drawing processes, these samples gain different morphologies, leading to different molecular dynamics. Due to the long chain nature, the molecular dynamics of polyethylene can be distinguished in local fluctuation and long range motion. With the help of NMR these different kinds of molecular dynamics can be monitored separately. In this work the local chain dynamics in non-crystalline regions of polyethylene samples was investigated via measuring 1H-13C heteronuclear dipolar coupling and 13C chemical shift anisotropy (CSA). By analyzing the motionally averaged 1H-13C heteronuclear dipolar coupling and 13C CSA, the information about the local anisotropy and geometry of motion was obtained. Taking advantage of the big difference of the 13C T1 relaxation time in crystalline and non-crystalline regions of PEs, the 1D 13C MAS exchange experiment was used to investigate the cooperative chain motion between these regions. The different chain organizations in non-crystalline regions were used to explain the relationship between the local fluctuation and the long range motion of the samples. In a simple manner the cooperative chain motion between crystalline and non-crystalline regions of PE results in the experimentally observed diffusive behavior of PE chain. The morphological influences on the diffusion motion have been discussed. The morphological factors include lamellar thickness, chain organization in non-crystalline regions and chain entanglements. Thermodynamics of the diffusion motion in melt and solution crystallized UHMW-PEs is discussed, revealing entropy-controlled features of the chain diffusion in PE. This thermodynamic consideration explains the counterintuitive relationship between the local fluctuation and the long range motion of the samples. Using the chain diffusion coefficient, the rates of jump motion in crystals of the melt crystallized PE have been calculated. A concept of "effective" jump motion has been proposed to explain the difference between the values derived from the chain diffusion coefficients and those in literatures. The observations of this thesis give a clear demonstration of the strong relationship between the sample morphology and chain dynamics. The sample morphologies governed by the processing history lead to different spatial constraints for the molecular chains, leading to different features of the local and long range chain dynamics. The knowledge of the morphological influence on the microscopic chain motion has many implications in our understanding of the alpha-relaxation process in PE and the related phenomena such as crystal thickening, drawability of PE, the easy creep of PE fiber, etc.

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).

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.

Development and applications of hollow fiber flow field-flow fractionation in the bioanalytical field. Studies of aggregation phenomena in complex protein samples

Recent advances in the fast growing area of therapeutic/diagnostic proteins and antibodies - novel and highly specific drugs - as well as the progress in the field of functional proteomics regarding the correlation between the aggregation of damaged proteins and (immuno) senescence or aging-related pathologies, underline the need for adequate analytical methods for the detection, separation, characterization and quantification of protein aggregates, regardless of the their origin or formation mechanism. Hollow fiber flow field-flow fractionation (HF5), the miniaturized version of FlowFFF and integral part of the Eclipse DUALTEC FFF separation system, was the focus of this research; this flow-based separation technique proved to be uniquely suited for the hydrodynamic size-based separation of proteins and protein aggregates in a very broad size and molecular weight (MW) range, often present at trace levels. HF5 has shown to be (a) highly selective in terms of protein diffusion coefficients, (b) versatile in terms of bio-compatible carrier solution choice, (c) able to preserve the biophysical properties/molecular conformation of the proteins/protein aggregates and (d) able to discriminate between different types of protein aggregates. Thanks to the miniaturization advantages and the online coupling with highly sensitive detection techniques (UV/Vis, intrinsic fluorescence and multi-angle light scattering), HF5 had very low detection/quantification limits for protein aggregates. Compared to size-exclusion chromatography (SEC), HF5 demonstrated superior selectivity and potential as orthogonal analytical method in the extended characterization assays, often required by therapeutic protein formulations. In addition, the developed HF5 methods have proven to be rapid, highly selective, sensitive and repeatable. HF5 was ideally suitable as first dimension of separation of aging-related protein aggregates from whole cell lysates (proteome pre-fractionation method) and, by HF5-(UV)-MALS online coupling, important biophysical information on the fractionated proteins and protein aggregates was gathered: size (rms radius and hydrodynamic radius), absolute MW and conformation.

Radiation-induced defects in calcium fluoride and their influence on material properties under 193 nm laser irradiation

Calcium fluoride (CaF2) is one of the key lens materials in deep-ultraviolet microlithography because of its transparency at 193 nm and its nearly perfect optical isotropy. Its physical and chemical properties make it applicable for lens fabrication. The key feature of CaF2 is its extreme laser stability. rnAfter exposing CaF2 to 193 nm laser irradiation at high fluences, a loss in optical performance is observed, which is related to radiation-induced defect structures in the material. The initial rapid damage process is well understood as the formation of radiation-induced point defects, however, after a long irradiation time of up to 2 months, permanent damage of the crystals is observed. Based on experimental results, these permanent radiation-induced defect structures are identified as metallic Ca colloids.rnThe properties of point defects in CaF2 and their stabilization in the crystal bulk are calculated with density functional theory (DFT). Because the stabilization of the point defects and the formation of metallic Ca colloids are diffusion-driven processes, the diffusion coefficients for the vacancy (F center) and the interstitial (H center) in CaF2 are determined with the nudged elastic band method. The optical properties of Ca colloids in CaF2 are obtained from Mie-theory, and their formation energy is determined.rnBased on experimental observations and the theoretical description of radiation-induced point defects and defect structures, a diffusion-based model for laser-induced material damage in CaF2 is proposed, which also includes a mechanism for annealing of laser damage. rn

An investigation into complex inorganic materials with Mössbauer spectroscopy

Mössbauer Spektroskopie ist ein unverzichtbares Instrument für die Bestimmung von Oxidationszuständen und für die Analyse von lokalen Ordnungsphänomenen von Mössbauer aktiven Atomen. Weil es sich um eine lokale Methode handelt können sowohl kristalline als auch amorphe Materialien untersucht werden. Die Kombination von lokaler Prüfung mit Mössbauer Spektroskopie und globaler Untersuchung z.B. mit Röntgendiffraktometrie ermöglicht die Studie von Ordnungseffekten von statistisch besetzten Positionen in einer geordneten Matrix. Das wurde hier eingesetzt um die lokale Umgebung in zwei Serien von Heuslerverbindungen, Co2-xFe1+xSi and Co2Mn1-xFexAl zu untersuchen. Für die Co2Mn1-xFexAl Serie wurde eine L21 geordnete Phase in einer insgesamt B2 geordneten Probe detektiert. Ein Wechsel von der AlCu2Mn zu der CuHg2Ti Struktur wurde für die Co2-xFe1+xSi Proben gefunden. Die Transformation von einem Glas zu einem keramischen Material wurde mit 119Sn Mössbauer Spektroskopie untersucht. Die höhere Ordnung in der Keramik wurde von einer kleiner werdenden Mössbauerlinienbreite begleitet. Demzufolge geben die Modifikationen der Sn Umgebungen klar die Transformation des gesamten Materials wieder. Ist die lokale Umgebung von unregelmäßig auftretenden Atomen in einer amorphen Matrix von Interesse, sind lokal prüfende Methoden die zuverlässigsten Methoden die zur Verfügung stehen. In dieser Arbeit wurde 119Sn Mössbauer Spektroskopie eingesetzt um die Oxidationszustände, die lokalen Umgebungen und relativen Intensitäten von Zinn Atomen in einer Silikatmatrix zu bestimmen. Modifikationen dieser Parameter als Funktion von Prozess bestimmenden Parametern wie der Sauerstoffpartialdruck, die Temperatur, die Behandlungsdauer und der Abkühlprozess genauso wie der SnO2 Gehalt sind von Interesse, weil durch Reduktions- und Diffusionsprozesse Änderungen des Koordinations- und des Oxidationszustands der Zinnatome auftreten. Da diese Änderungen in der Glasmatrix verursachen, die das fertige Produkt im industriellen Fertigungsprozess ruinieren können sind diese feinen Veränderungen sehr wichtig. Wenigstens zwei Mössbauerlinien korrespondierend mit zwei verschiedenen Umgebungen für Sn2+ und Sn4+ sind für eine Analyse mit ausreichender Qualität notwendig. Durch Vergleich von den bestimmten Hyperfein Parametern mit den Parametern von Modelsubstanzen werden lokale Umgebungen der Zinnatome entworfen. Für Sn2+ werden zwei auf einer trigonalen Pyramide basierende Umgebungen mit variierender Anzahl von bindenden und nicht-bindenden Sauerstoffatomen formuliert. Für Sn4+ wurde eine tetraedrische und eine oktaedrische Umgebung postuliert. Die relativen Intensitäten der vier Mössbauerlinien wurden um ein Diffusions- und Reaktionsmodell zu entwickeln und um einen Satz von Diffusions- und Transferkoeffizienten zu bestimmen eingesetzt. Die bestimmten Diffusionskoeffizienten stimmen mit den Literaturdaten überein. Der Massentransferkoeffizient ist kleiner als der bestimmte Wert, aber immer noch in der gleichen Größenordnung. Im Gegensatz zu den Erwartungen ist der präsentierte Diffusionskoeffizient für Sn4+ bestimmt als der von Sn2+. Das wiederum kann durch Berücksichtigung von Elektronhoppingprozessen erklärt werden.

Assessing the Validity of Brain Glucose Concentration Measurement Using Microdialysis

Brain functions, such as learning, orchestrating locomotion, memory recall, and processing information, all require glucose as a source of energy. During these functions, the glucose concentration decreases as the glucose is being consumed by brain cells. By measuring this drop in concentration, it is possible to determine which parts of the brain are used during specific functions and consequently, how much energy the brain requires to complete the function. One way to measure in vivo brain glucose levels is with a microdialysis probe. The drawback of this analytical procedure, as with many steadystate fluid flow systems, is that the probe fluid will not reach equilibrium with the brain fluid. Therefore, brain concentration is inferred by taking samples at multiple inlet glucose concentrations and finding a point of convergence. The goal of this thesis is to create a three-dimensional, time-dependent, finite element representation of the brainprobe system in COMSOL 4.2 that describes the diffusion and convection of glucose. Once validated with experimental results, this model can then be used to test parameters that experiments cannot access. When simulations were run using published values for physical constants (i.e. diffusivities, density and viscosity), the resulting glucose model concentrations were within the error of the experimental data. This verifies that the model is an accurate representation of the physical system. In addition to accurately describing the experimental brain-probe system, the model I created is able to show the validity of zero-net-flux for a given experiment. A useful discovery is that the slope of the zero-net-flux line is dependent on perfusate flow rate and diffusion coefficients, but it is independent of brain glucose concentrations. The model was simplified with the realization that the perfusate is at thermal equilibrium with the brain throughout the active region of the probe. This allowed for the assumption that all model parameters are temperature independent. The time to steady-state for the probe is approximately one minute. However, the signal degrades in the exit tubing due to Taylor dispersion, on the order of two minutes for two meters of tubing. Given an analytical instrument requiring a five μL aliquot, the smallest brain process measurable for this system is 13 minutes.

Dynamics of supercooled water in highly compacted clays studied by neutron scattering

The freezing behavior of water confined in compacted charged and uncharged clays (montmorillonite in Na-and Ca-forms, illite in Na-and Ca-forms, kaolinite and pyrophyllite) was investigated by neutron scattering. Firstly, the amount of frozen (immobile) water was measured as a function of temperature at the IN16 backscattering spectrometer, Institute Laue-Langevin (ILL). Water in uncharged, partly hydrophobic (kaolinite) and fully hydrophobic (pyrophyllite) clays exhibited a similar freezing and melting behavior to that of bulk water. In contrast, water in charged clays which are hydrophilic could be significantly supercooled. To observe the water dynamics in these clays, further experiments were performed using quasielastic neutron scattering. At temperatures of 250, 260 and 270 K the diffusive motion of water could still be observed, but with a strong reduction in the water mobility as compared with the values obtained above 273 K. The diffusion coefficients followed a non-Arrhenius temperature dependence well described by the Vogel-Fulcher-Tammann and the fractional power relations. The fits revealed that Na-and Ca-montmorillonite and Ca-illite have similar Vogel-Fulcher-Tammann temperatures (T-VFT, often referred to as the glass transition temperature) of similar to 120 K and similar temperatures at which the water undergoes the 'strong-fragile' transition, T-s similar to 210 K. On the other hand, Na-illite had significantly larger values of T-VFT similar to 180 K and T-s similar to 240 K. Surprisingly, Ca-illite has a similar freezing behavior of water to that of montmorillonites, even though it has a rather different structure. We attribute this to the stronger hydration of Ca ions as compared with the Na ions occurring in the illite clays.

Multi-scale micro-structure generation strategy for up-scaling transport in clays

Clays and claystones are used as backfill and barrier materials in the design of waste repositories, because they act as hydraulic barriers and retain contaminants. Transport through such barriers occurs mainly by molecular diffusion. There is thus an interest to relate the diffusion properties of clays to their structural properties. In previous work, we have developed a concept for up-scaling pore-scale molecular diffusion coefficients using a grid-based model for the sample pore structure. Here we present an operational algorithm which can generate such model pore structures of polymineral materials. The obtained pore maps match the rock’s mineralogical components and its macroscopic properties such as porosity, grain and pore size distributions. Representative ensembles of grains in 2D or 3D are created by a lattice Monte Carlo (MC) method, which minimizes the interfacial energy of grains starting from an initial grain distribution. Pores are generated at grain boundaries and/or within grains. The method is general and allows to generate anisotropic structures with grains of approximately predetermined shapes, or with mixtures of different grain types. A specific focus of this study was on the simulation of clay-like materials. The generated clay pore maps were then used to derive upscaled effective diffusion coefficients for non-sorbing tracers using a homogenization technique. The large number of generated maps allowed to check the relations between micro-structural features of clays and their effective transport parameters, as is required to explain and extrapolate experimental diffusion results. As examples, we present a set of 2D and 3D simulations and investigated the effects of nanopores within particles (interlayer pores) and micropores between particles. Archie’s simple power law is followed in systems with only micropores. When nanopores are present, additional parameters are required; the data reveal that effective diffusion coefficients could be described by a sum of two power functions, related to the micro- and nanoporosity. We further used the model to investigate the relationships between particle orientation and effective transport properties of the sample.

Water–montmorillonite systems: Neutron scattering and tracer throughdiffusion studies

Designs for deep geological respositories of nuclear waste include bentonite as a hydraulic and chemisorption buffer material to protect the biosphere from leakage of radionuclides. Bentonite is chosen because it is a cheap, naturally occurring material with the required properties. It consists essentially of montmorillonite, a swelling clay mineral. Upon contact with groundwater such clays can seal the repository by incorporating water in the interlayers of their crystalline structure. The intercalated water exhibits significantly different properties to bulk water in the surrounding interparticle pores, such as lower diffusion coefficients (González Sánchez et. al. 2008). This doctoral thesis presents water distribution and diffusion behavior on various time and space scales in montmorillonite. Experimental results are presented for Na- and Cs-montmorillonite samples with a range of bulk dry densities (0.8 to 1.7 g/cm3). The experimental methods employed were neutron scattering (backscattering, diffraction, time-of-flight), adsorption measurements (water, nitrogen) and tracer-through diffusion. For the tracer experiments the samples were fully saturated via the liquid phase under volume-constrained conditions. In contrast, for the neutron scattering experiments, the samples were hydrated via the vapor phase and subsequently compacted, leaving a significant fraction of interparticle pores unfilled with water. Owing to these differences in saturation, the water contents of the samples for neutron scattering were characterized by gravimetry whereas those for the tracer experiments were obtained from the bulk dry density. The amount of surface water in interlayer pores could be successfully discriminated from the amount of bulk-like water in interparticle pores in Na- and Csmontmorillonite using neutron spectroscopy. For the first time in the literature, the distribution of water between these two pore environments was deciphered as a function of gravimetric water content. The amount was compared to a geometrical estimation of the amount of interlayer and interparticle water determined by neutron diffraction and adsorption measurements. The relative abundances of the 1 to 4 molecular water layers in the interlayer were determined from the area ratios of the (001)-diffraction peaks. Depending on the characterization method, different fractions of surface water and interlayer water were obtained. Only surface and interlayer water exists in amontmorillonite with water contents up to 0.18 g/g according to spectroscopic measurements and up to 0.32 g/g according to geometrical estimations, respectively. At higher water contents, bulk-like and interparticle water also exists. The amounts increase monotonically, but not linearly, from zero to 0.33 g/g for bulk-like water and to 0.43 g/g for interparticle water. It was found that water most likely redistributes between the surface and interlayer sites during the spectroscopic measurements and therefore the reported fraction is relevant only below about -10 ºC (Anderson, 1967). The redistribution effect can explain the discrepancy in fractions between the methods. In a novel approach the fractions of water in different pore environments were treated as a fixed parameter to derive local diffusion coefficients for water from quasielastic neutron scattering data, in particular for samples with high water contents. Local diffusion coefficients were obtained for the 1 to 4 molecular water layers in the interlayer of 0.5·10–9, 0.9·10–9, 1.5·10–9 and 1.4·10–9 m²/s, respectively, taking account of the different water fractions (molecular water layer, bulk-like water). The diffusive transport of 22Na and HTO through Na-montmorillonite was measured on the laboratory experimental scale (i.e. cm, days) by tracer through-diffusion experiments. We confirmed that diffusion of HTO is independent of the ionic strength of the external solution in contact with the clay sample but dependent on the bulk dry density. In contrast, the diffusion of 22Na was found to depend on both the ionic strength of the pore solution and on the bulk dry density. The ratio of the pore and surface diffusion could be experimentally determined for 22Na from the dependence of the diffusion coefficient on the ionic strength. Activation energies were derived from the temperaturedependent diffusion coefficients via the Arrhenius relation. In samples with high bulk dry density the activation energies are slightly higher than those of bulk water whereas in low density samples they are lower. The activation energies as a function of ionic strengths of the pore solutions are similar for 22Na and HTO. The facts that (i) the slope of the logarithmic effective diffusion coefficients as a function of the logarithmic ionic strength is less than unity for low bulk dry densities and (ii) two water populations can be observed for high gravimetric water contents (low bulk dry densities) support the interlayer and interparticle porosity model proposed by Glaus et al. (2007), Bourg et al. (2006, 2007) and Gimmi and Kosakowski (2011).

Chemistry of basalts of ODP Hole 111-504B

The effect of oxygen fugacity (fO2) on the partition relationship of Mg and Fe between Plagioclase and sillicate liquid was investigated at 1 atm for basaltic samples recovered during ODP Leg 111 from Hole 504B. Samples 111-504B-143R-2 (Piece 8) and 111-504B-169R-1 (Piece 1) have Plagioclase as the liquidus phase. The distribution coefficient of Mg between Plagioclase and melt is constant at about 0.04 against the variation of fO2, whereas that of Fe (total Fe) varies from 0.3 at f(O2) = 0.2 atm to 0.03 at f(o2) = 10**-11.5 at 1200°C. The distribution coefficient of Mg is slightly higher than that calculated from the phenocryst and bulk-rock compositions, suggesting a kinetic disequilibrium effect on the distribution of Mg in Plagioclase. Because Mg, Fe, and Fe3+ have similar diffusion coefficients in silicate melt, the disequilibrium effect is greatly reduced for the exchange reaction of Mg and total Fe between Plagioclase and liquid. The exchange partition coefficient is highly dependent on fo2, with log fo2 ranging from -0.7 to - 11.5 at approximately 1200°C. Using this relationship, the f(O2) of crystallization of the magmas is estimated to be near the one defined by the fayalite-quartz-magnetite assemblage.

Liquid lead and Lithium. A Comprehensive molecular dynamics study

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. One of the main issues in these programs is the problem of liquid metals breeder blanket behavior. Structural material of the blanket should meet high requirements because of extreme operating conditions. Therefore the knowledge of eutectic properties like optimal composition, physical and thermodynamic behavior or diffusion coefficients of Tritium are extremely necessary for current designs. In particular, the knowledge of the function linking the tritium concentration dissolved in liquid materials with the tritium partial pressure at a liquid/gas interface in equilibrium, CT=f(PT), is of basic importance because it directly impacts all functional properties of a blanket determining: tritium inventory, tritium permeation rate and tritium extraction efficiency. Nowadays, understanding the structure and behavior of this compound is a real goal in fusion engineering and materials science. Simulations of liquids can provide much information to the community; not only supplementing experimental data, but providing new tests of theories and ideas, making specific predictions that require experimental tests, and ultimately helping to lead to the deeper understanding and better predictive behavior.