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.

Dispersive approach to hadronic light-by-light scattering

Based on dispersion theory, we present a formalism for a model-independent evaluation of the hadronic light-by-light contribution to the anomalous magnetic moment of the muon. In particular, we comment on the definition of the pion pole in this framework and provide a master formula that relates the effect from ππ intermediate states to the partial waves for the process γ * γ * → ππ. All contributions are expressed in terms of on-shell form factors and scattering amplitudes, and as such amenable to an experimental determination.

New Constraints on Dark Matter Effective Theories from Standard Model Loops

We consider an effective field theory for a gauge singlet Dirac dark matter particle interacting with the standard model fields via effective operators suppressed by the scale Λ≳1  TeV. We perform a systematic analysis of the leading loop contributions to spin-independent Dirac dark matter–nucleon scattering using renormalization group evolution between Λ and the low-energy scale probed by direct detection experiments. We find that electroweak interactions induce operator mixings such that operators that are naively velocity suppressed and spin dependent can actually contribute to spin-independent scattering. This allows us to put novel constraints on Wilson coefficients that were so far poorly bounded by direct detection. Constraints from current searches are already significantly stronger than LHC bounds, and will improve in the near future. Interestingly, the loop contribution we find is isospin violating even if the underlying theory is isospin conserving.

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

Dispersion relation for hadronic light-by-light scattering: theoretical foundations

In this paper we make a further step towards a dispersive description of the hadronic light-by-light (HLbL) tensor, which should ultimately lead to a data-driven evaluation of its contribution to (g − 2) μ . We first provide a Lorentz decomposition of the HLbL tensor performed according to the general recipe by Bardeen, Tung, and Tarrach, generalizing and extending our previous approach, which was constructed in terms of a basis of helicity amplitudes. Such a tensor decomposition has several advantages: the role of gauge invariance and crossing symmetry becomes fully transparent; the scalar coefficient functions are free of kinematic singularities and zeros, and thus fulfill a Mandelstam double-dispersive representation; and the explicit relation for the HLbL contribution to (g − 2) μ in terms of the coefficient functions simplifies substantially. We demonstrate explicitly that the dispersive approach defines both the pion-pole and the pion-loop contribution unambiguously and in a model-independent way. The pion loop, dispersively defined as pion-box topology, is proven to coincide exactly with the one-loop scalar QED amplitude, multiplied by the appropriate pion vector form factors.

MEMLS3&a: Microwave Emission Model of Layered Snowpacks adapted to include backscattering

The Microwave Emission Model of Layered Snowpacks (MEMLS) was originally developed for microwave emissions of snowpacks in the frequency range 5–100 GHz. It is based on six-flux theory to describe radiative transfer in snow including absorption, multiple volume scattering, radiation trapping due to internal reflection and a combination of coherent and incoherent superposition of reflections between horizontal layer interfaces. Here we introduce MEMLS3&a, an extension of MEMLS, which includes a backscatter model for active microwave remote sensing of snow. The reflectivity is decomposed into diffuse and specular components. Slight undulations of the snow surface are taken into account. The treatment of like- and cross-polarization is accomplished by an empirical splitting parameter q. MEMLS3&a (as well as MEMLS) is set up in a way that snow input parameters can be derived by objective measurement methods which avoid fitting procedures of the scattering efficiency of snow, required by several other models. For the validation of the model we have used a combination of active and passive measurements from the NoSREx (Nordic Snow Radar Experiment) campaign in Sodankylä, Finland. We find a reasonable agreement between the measurements and simulations, subject to uncertainties in hitherto unmeasured input parameters of the backscatter model. The model is written in Matlab and the code is publicly available for download through the following website: http://www.iapmw.unibe.ch/research/projects/snowtools/memls.html.

From quantum to classical dynamics: The relativistic O ( N ) model in the framework of the real-time functional renormalization group

We investigate the transition from unitary to dissipative dynamics in the relativistic O(N) vector model with the λ(φ2)2 interaction using the nonperturbative functional renormalization group in the real-time formalism. In thermal equilibrium, the theory is characterized by two scales, the interaction range for coherent scattering of particles and the mean free path determined by the rate of incoherent collisions with excitations in the thermal medium. Their competition determines the renormalization group flow and the effective dynamics of the model. Here we quantify the dynamic properties of the model in terms of the scale-dependent dynamic critical exponent z in the limit of large temperatures and in 2≤d≤4 spatial dimensions. We contrast our results to the behavior expected at vanishing temperature and address the question of the appropriate dynamic universality class for the given microscopic theory.

Hadronic light-by-light scattering and the muon g−2

The largest uncertainties in the Standard Model calculation of the anomalous magnetic moment of the muon (g − 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, which is based on unitarity, analyticity, crossing symmetry, and gauge invariance. Such a model-independent Approach opens up an avenue towards a data-driven determination of the HLbL contribution to the (g − 2)μ.

Circuital model for the spherical geodesic waveguide perfect drain

The perfect drain for the Maxwell fish eye (MFE) is a non-magnetic dissipative region placed in the focal point to absorb all the incident radiation without reflection or scattering.

Comparison of dynamic effects of high-speed traffic load on ballasted track using a simplified two-dimensional and full three-dimensional model.

The vertical dynamic actions transmitted by railway vehicles to the ballasted track infrastructure is evaluated taking into account models with different degree of detail. In particular, we have studied this matter from a two-dimensional (2D) finite element model to a fully coupled three-dimensional (3D) multi-body finite element model. The vehicle and track are coupled via a non-linear Hertz contact mechanism. The method of Lagrange multipliers is used for the contact constraint enforcement between wheel and rail. Distributed elevation irregularities are generated based on power spectral density (PSD) distributions which are taken into account for the interaction. The numerical simulations are performed in the time domain, using a direct integration method for solving the transient problem due to the contact nonlinearities. The results obtained include contact forces, forces transmitted to the infrastructure (sleeper) by railpads and envelopes of relevant results for several track irregularities and speed ranges. The main contribution of this work is to identify and discuss coincidences and differences between discrete 2D models and continuum 3D models, as wheel as assessing the validity of evaluating the dynamic loading on the track with simplified 2D models

Experimental and theoretical cross sections for positron scattering from the pentane isomers.

Isomerism is ubiquitous in chemistry, physics, and biology. In atomic and molecular physics, in particular, isomer effects are well known in electron-impact phenomena; however, very little is known for positron collisions. Here we report on a set of experimental and theoretical cross sections for low-energy positron scattering from the three structural isomers of pentane: normal-pentane, isopentane, and neopentane. Total cross sections for positron scattering from normal-pentane and isopentane were measured at the University of Trento at incident energies between 0.1 and 50 eV. Calculations of the total cross sections, integral cross sections for elastic scattering, positronium formation, and electronic excitations plus direct ionization, as well as elastic differential cross sections were computed for all three isomers between 1 and 1000 eV using the independent atom model with screening corrected additivity rule. No definitive evidence of a significant isomer effect in positron scattering from the pentane isomers appears to be present. (C) 2016 AIP Publishing LLC.

Effective nonlinear model of resonant tunneling nanostructures

We introduce a model of a nonlinear double-barrier structure to describe in a simple way the effects of electron-electron scattering while remaining analytically tractable. The model is based on a generalized effective-mass equation where a nonlinear local field interaction is introduced to account for those inelastic scattering phenomena. Resonance peaks seen in the transmission coefficient spectra for the linear case appear shifted to higher energies depending on the magnitude of the nonlinear coupling. Our results are in good agreement with self-consistent solutions of the Schrodinger and Poisson equations. The calculation procedure is seen to be very fast, which makes our technique a good candidate for a rapid approximate analysis of these structures.

Integral elastic, electronic-state, ionization, and total cross sections for electron scattering with furfural

We report absolute experimental integral cross sections (ICSs) for electron impact excitation of bands of electronic-states in furfural, for incident electron energies in the range 20-250 eV. Wherever possible, those results are compared to corresponding excitation cross sections in the structurally similar species furan, as previously reported by da Costa et al. [Phys. Rev. A 85, 062706 (2012)] and Regeta and Allan [Phys. Rev. A 91, 012707 (2015)]. Generally, very good agreement is found. In addition, ICSs calculated with our independent atom model (IAM) with screening corrected additivity rule (SCAR) formalism, extended to account for interference (I) terms that arise due to the multi-centre nature of the scattering problem, are also reported. The sum of those ICSs gives the IAM-SCAR+I total cross section for electron-furfural scattering. Where possible, those calculated IAM-SCAR+I ICS results are compared against corresponding results from the present measurements with an acceptable level of accord being obtained. Similarly, but only for the band I and band II excited electronic states, we also present results from our Schwinger multichannel method with pseudopotentials calculations. Those results are found to be in good qualitative accord with the present experimental ICSs. Finally, with a view to assembling a complete cross section data base for furfural, some binary-encounter-Bethe-level total ionization cross sections for this collision system are presented. (C) 2016 AIP Publishing LLC.

Polarimetric SAR model for soil moisture estimation over vineyards at C-band

In this paper we propose a two-component polarimetric model for soil moisture estimation on vineyards suited for C-band radar data. According to a polarimetric analysis carried out here, this scenario is made up of one dominant direct return from the soil and a multiple scattering component accounting for disturbing and nonmodeled signal fluctuations from soil and short vegetation. We propose a combined X-Bragg/Fresnel approach to characterize the polarized direct response from soil. A validation of this polarimetric model has been performed in terms of its consistency with respect to the available data both from RADARSAT-2 and from indoor measurements. High inversion rates are reported for different phenological stages of vines, and the model gives a consistent interpretation of the data as long as the volume component power remains about or below 50% of the surface contribution power. However, the scarcity of soil moisture measurements in this study prevents the validation of the algorithm in terms of the accuracy of soil moisture retrieval and an extensive campaign is required to fully demonstrate the validity of the model. Different sources of mismatches between the model and the data have been also discussed and analyzed.