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

Analyzing diffuse scattering with supercomputers

Two new approaches to quantitatively analyze diffuse diffraction intensities from faulted layer stacking are reported. The parameters of a probability-based growth model are determined with two iterative global optimization methods: a genetic algorithm (GA) and particle swarm optimization (PSO). The results are compared with those from a third global optimization method, a differential evolution (DE) algorithm [Storn & Price (1997). J. Global Optim. 11, 341–359]. The algorithm efficiencies in the early and late stages of iteration are compared. The accuracy of the optimized parameters improves with increasing size of the simulated crystal volume. The wall clock time for computing quite large crystal volumes can be kept within reasonable limits by the parallel calculation of many crystals (clones) generated for each model parameter set on a super- or grid computer. The faulted layer stacking in single crystals of trigonal three-pointedstar- shaped tris(bicylco[2.1.1]hexeno)benzene molecules serves as an example for the numerical computations. Based on numerical values of seven model parameters (reference parameters), nearly noise-free reference intensities of 14 diffuse streaks were simulated from 1280 clones, each consisting of 96 000 layers (reference crystal). The parameters derived from the reference intensities with GA, PSO and DE were compared with the original reference parameters as a function of the simulated total crystal volume. The statistical distribution of structural motifs in the simulated crystals is in good agreement with that in the reference crystal. The results found with the growth model for layer stacking disorder are applicable to other disorder types and modeling techniques, Monte Carlo in particular.

Accurate evaluation of hadronic uncertainties in spin-independent WIMP-nucleon scattering: Disentangling two- and three-flavor effects

We show how to avoid unnecessary and uncontrolled assumptions usually made in the literature about soft SU(3) flavor symmetry breaking in determining the two-flavor nucleon matrix elements relevant for direct detection of weakly interacting massive particles (WIMPs). Based on SU(2) chiral perturbation theory, we provide expressions for the proton and neutron scalar couplings fp,nu and fp,nd with the pion-nucleon σ term as the only free parameter, which should be used in the analysis of direct detection experiments. This approach for the first time allows for an accurate assessment of hadronic uncertainties in spin-independent WIMP-nucleon scattering and for a reliable calculation of isospin-violating effects. We find that the traditional determinations of Vfpu−fnu and fpd−fnd are off by a factor of 2.

The multiple-scattering series in few-nucleon systems

We discuss under which circumstances the resummation of the multiple-scattering series is justified from an EFT point of view. The application to πd and K̅d scattering is briefly discussed.

Roy-Steiner equations for piN scattering

In this talk, we present a coupled system of integral equations for the πN → πN (s-channel) and ππ → N̅N (t-channel) lowest partial waves, derived from Roy–Steiner equations for pion–nucleon scattering. After giving a brief overview of this system of equations, we present the solution of the t-channel sub-problem by means of Muskhelishvili–Omnès techniques, and solve the s-channel sub-problem after finding a set of phase shifts and subthreshold parameters which satisfy the Roy–Steiner equations.

Dispersive approach to hadronic light-by-light scattering and the muon ɡ − 2

The largest uncertainties in the Standard Model calculation of the anomalous magnetic moment of the muon (ɡ − 2)μ come from hadronic contributions. In particular, it can be expected that in a few years the subleading hadronic light-by-light (HLbL) contribution will dominate the theory uncertainty. We present a dispersive description of the HLbL tensor. This new, model-independent approach opens up an avenue towards a data-driven determination of the HLbL contribution to the (ɡ − 2)μ.

Virtual photon-photon scattering

Based on analyticity, unitarity, and Lorentz invariance the contribution from hadronic vacuum polarization to the anomalous magnetic moment of the muon is directly related to the cross section of e+e− → hadrons. We review the main difficulties that impede such an approach for light-by-light scattering and identify the required ingredients from experiment. Amongst those, the most critical one is the scattering of two virtual photons into meson pairs. We analyze the analytic structure of the process γ*γ* → ππ and show that the usual Muskhelishvili–Omnès representation can be amended in such a way as to remain valid even in the presence of anomalous thresholds.

First πK atom lifetime and πK scattering length measurements

The results of a search for hydrogen-like atoms consisting of π∓K±π∓K± mesons are presented. Evidence for πK atom production by 24 GeV/c protons from CERN PS interacting with a nickel target has been seen in terms of characteristic πK pairs from their breakup in the same target (178±49178±49) as well as in terms of produced πK atoms (653±42653±42). Using these results, the analysis yields a first value for the πK atom lifetime of View the MathML sourceτ=(2.5−1.8+3.0) fs and a first measurement of the S-wave isospin-odd πK scattering length View the MathML source|a0−|=13|a1/2−a3/2|=(0.11−0.04+0.09)Mπ−1 (aIaI for isospin I).

Measurement of the total cross section from elastic scattering in pp collisions at √s=7TeV with the ATLAS detector

A measurement of the total pp cross section at the LHC at √s = 7 TeV is presented. In a special run with high-β* beam optics, an integrated luminosity of 80 μb−1 was accumulated in order to measure the differential elastic cross section as a function of the Mandelstam momentum transfer variable t . The measurement is performed with the ALFA sub-detector of ATLAS. Using a fit to the differential elastic cross section in the |t | range from 0.01 GeV2 to 0.1 GeV2 to extrapolate to |t | →0, the total cross section, σtot(pp→X), is measured via the optical theorem to be: σtot(pp→X) = 95.35± 0.38 (stat.)± 1.25 (exp.)± 0.37 (extr.) mb, where the first error is statistical, the second accounts for all experimental systematic uncertainties and the last is related to uncertainties in the extrapolation to |t | → 0. In addition, the slope of the elastic cross section at small |t | is determined to be B = 19.73 ±0.14 (stat.) ±0.26 (syst.) GeV−2.

Measurement of the parity-violating asymmetry parameter αb and the helicity amplitudes for the decay Λ 0 b →J/ψ⁰ with the ATLAS detector

A measurement of the parity-violating decay asymmetry parameter, αb , and the helicity amplitudes for the decay Λ 0 b →J/ψ(μ + μ − )Λ 0 (pπ − ) is reported. The analysis is based on 1400 Λ 0 b and Λ ¯ 0 b baryons selected in 4.6  fb −1 of proton–proton collision data with a center-of-mass energy of 7 TeV recorded by the ATLAS experiment at the LHC. By combining the Λ 0 b and Λ ¯ 0 b samples under the assumption of CP conservation, the value of α b is measured to be 0.30±0.16(stat)±0.06(syst) . This measurement provides a test of theoretical models based on perturbative QCD or heavy-quark effective theory.

Probing the moduli dependence of refined topological amplitudes

With the aim of providing a worldsheet description of the refined topological string, we continue the study of a particular class of higher derivative couplings Fg,n in the type II string effective action compactified on a Calabi–Yau threefold. We analyse first order differential equations in the anti-holomorphic moduli of the theory, which relate the Fg,n to other component couplings. From the point of view of the topological theory, these equations describe the contribution of non-physical states to twisted correlation functions and encode an obstruction for interpreting the Fg,n as the free energy of the refined topological string theory. We investigate possibilities of lifting this obstruction by formulating conditions on the moduli dependence under which the differential equations simplify and take the form of generalised holomorphic anomaly equations. We further test this approach against explicit calculations in the dual heterotic theory.