Eric Ferrage

Dehydration of Ca-montmorillonite at the crystal scale. Part I: Structure evolution

Abstract The dehydration dynamics of the of the Ca-saturated <1 μm size fraction of SWy-1 (low-charge montmorillonite) were studied at the crystal scale under isothermal conditions using X-ray diffraction with a position-sensitive detector (XRD-PSD) in the 30.170 °C temperature range. A total of 630 XRD patterns were modeled between 30 and 125 °C using a trial-and-error approach based on the direct comparison of experimental and calculated XRD patterns. The proportion of layers with different hydration states (bihydrated, mono-hydrated, and dehydrated) were determined in the temperature-time space as well as small variations of layer thickness within each hydration state. The results showed that dehydration produces complex structures with heterogeneous hydration states, some of which are stable (not transient) and remain at the end of the experiment. The evolution of other structural parameters (interlayer water content, layer thickness fluctuation) was consistent with previous reports of smectite hydration. For bihydrated layers, the amount of water molecules per interlayer cation indicated the presence of water molecules both coordinated and non-coordinated to the interlayer cation. The transition from bi- to mono-hydrated layers produced the maximum structural heterogeneity, with (1) strong interlayer thickness fluctuation (in individual layers), and (2) the presence of several elementary mixed-layer structures. In contrast, the transition from mono-hydrated to dehydrated layers occurs homogeneously within layers. Finally, the decrease in thickness of mono-hydrated layers only implied the removal of some water molecules forming the hydration shell of the interlayer cation.

A reinvestigation of smectite illitization in experimental hydrothermal conditions: Results from X-ray diffraction and transmission electron microscopy

Abstract The hydrothermal reactivity of the <1 μm fraction, K-saturated SWy-2 Wyoming low-charge montmorillonite was studied in the 250-400 °C temperature range with reactions lasting between 5 and 120 days, with a solid/solution mass ratio of 1:10, and in 1 mol/L KCl solution. From X-ray diffraction (XRD) profile modeling results on K-saturated and ethylene-glycol solvated samples the illitization process appears to occur as a progressive replacement of expandable layers by layers with illitic behavior, in a single illite-smectite phase. However this treatment overestimates the amount of illite layers because of the presence of smectitic non-expandable layers. This was revealed by calcium exchange of the products, which causes re-expansion of the apparent illite layers. The illitization model then obtained consists of four phases with increasing illite content: discrete smectite, a randomly interstratified mixed-layered mineral (MLM) and two ordered MLMs that progressively replace the phases with less illite content. This polyphase model is consistent with results from transmission electron microscopy and chemical microanalysis, which show several types of particle morphology with different interlayer K content and Al-for-Si substitution that seem to correspond to the several phases detected by XRD. Thus, the actual illitization mechanism of smectite in our experiments is not a progressive reaction but a dissolution-precipitation process following the Ostwald step rule in which metastable smectite transforms into illite through a series of metastable illite-smectite phases.

Influence of the ionic strength and solid/solution ratio on Ca(II)-for-Na+ exchange on montmorillonite. Part 1: Chemical measurements, thermodynamic modeling and potential implications for trace elements geochemistry

Na/Ca ion-exchange isotherms were performed on a Na-saturated montmorillonite for different constant normalities of the aqueous chloride solution and at two solid/solution ratios. The experimental data suggest that the affinity of Na(+) for the solid increases with total normality and m/V ratio and that a significant proportion of CaCl(+) is also sorbed. A thermodynamic modeling procedure with one sorption site and three sorbed cations (Na(+), Ca(2+) and CaCl(+)) was applied to interpret the data. We show that by accounting for the activities of aqueous species, the experimental data obtained at different total normalities for a given solid/solution ratio can be fitted using a unique set of selectivity coefficient values. However, when the m/V ratio is decreased from 25 to 2.5 g/L, an increase in the log Kc(Na(+)/Ca(2+)) of up to 0.4 ± 0.05 was required to interpret the data with a constant log Kc(Na(+)/CaCl(+)) value. The same behavior concerning the increase in the log Kc(Na(+)/Ca(2+)) with a m/V ratio decrease was observed when using a multi-site model taking into account the specific sorption of H(+). The results clearly indicate that the Na(+)/Ca(2+) selectivity coefficients strongly depend on the solid/solution ratio but are independent of the exchanger composition and the total normality of the solution. Such findings provide possible pitfalls when using selectivity coefficients obtained in diluted suspension to assess the sorption in compacted clayrocks and allow a better prediction of the sorption of trace elements in competition with major cations.

Unraveling complex <2 μm clay mineralogy from soils using X-ray diffraction profile modeling on particle-size sub-fractions: Implications for soil pedogenesis and reactivity

Abstract A specific methodology was developed to refine the complex clay mineralogy commonly encountered in soil environments. The soil examined was a Cambisol developed into a ferralitic paleosol. The sample was split into four sub-fractions of different particle sizes (<0.05, 0.05-0.1, 0.1-0.2, and 0.2-2 μm), and their respective mass contributions to the overall <2 μm clay fraction were determined. For each sub-fraction, X-ray diffraction (XRD) patterns were modeled using a trial-and-error approach based on the direct comparison of experimental and calculated profiles. Quantitative information derived from the fitting procedure for the different sub-fractions allowed for the determination of the complex mineralogy of the <2 μm clay fraction through the identification and quantification of eight clay phases. The results show that the finest and most reactive clay fraction (<0.05 μm) was totally hidden in the XRD pattern of the <2 μm fraction, the fraction commonly considered in soil mineralogical analyses. Similarly, this procedure revealed the presence of illite-smectite-chlorite and kaolinite-illite mixed-layer minerals seldom described in soil literature using classical methods. The use of this methodology improved our understanding of the pedogenesis of this soil through the identification and quantification of clay phases structural properties. The analysis of the evolution of structural parameters with particle size allowed for the detection of local modifications in the interlayer composition of expandable and hydroxy-interlayered vermiculite layers. Following this approach, key information can be derived to determine subtle changes in clay mineralogical composition that are related to microorganism and/or plant activity.

Dehydration of Ca-montmorillonite at the crystal scale. Part 2. Mechanisms and kinetics

Abstract A kinetic study of Ca-montmorillonite dehydration was performed based on information derived from X-ray diffraction (Ferrage et al. 2007, this issue) and, thus, focusing on interlayer water only. The dehydration was quantified following the two processes that were observed in the X-ray pattern modeling: the transitions between the different hydration states and small thickness decrease observed in the bi- and mono-hydrated layers. The thickness decrease of bihydrated layers with dehydration (activation energy Ea = 16 kJ/mol) was found to be controlled by a mechanism of two-dimensional diffusion of water molecules through the interlayer space, whereas for mono-hydrated layers the variation of thickness (Ea = 18 kJ/mol) occurred as a mechanism of slight local layer collapse and collapse propagation, attributed to a rearrangement of the configuration of the interlayer cation hydration shell. For the transition between the bi- and mono-hydrated state (Ea = 84 kJ/mol), the mechanism of reaction was found to evolve gradually with increasing temperature from local layer collapse and collapse propagation to a two-dimensional diffusion mechanism, as the forced diffusion of water molecules produced by the layer collapse transfers the control of the process to diffusion mechanism. This phenomenon causes the coexistence of two hydration states in a given interlayer. Finally, the transition between mono-hydrated and dehydrated layers (Ea = 132 kJ/mol) indicated the concomitance of water diffusion and local layer collapse and propagation mechanisms, although the structures were found to be homogeneous during this transition. The determination of both mechanisms and the activation energy for these processes were used to establish a model of smectite dehydration at the crystal scale. This model can be used to calculate crystal shrinkage and interlayer water content upon dehydration, and to predict the evolution of the system.

Influence of the ionic strength and solid/solution ratio on Ca(II)-for-Na+ exchange on montmorillonite. Part 2: Understanding the effect of the m/V ratio. Implications for pore water composition and element transport in natural media.

The aim of the present paper is to clarify previous results showing that selectivity coefficients determined for the exchange of Na(+) for Ca(2+) in montmorillonite were dependent on the solid/solution ratio. The organization of montmorillonite suspensions upon Na(+)/Ca(II) exchange was analyzed by combining optical microscopy, small-angle X-ray scattering and X-ray diffraction. All samples displayed flocculated characteristics, eliminating the possibility of contrasting accessibility of sorption sites with the solid/solution ratio. Modeling of experimental X-ray diffraction patterns was used to quantify the relative proportions of interlayer Ca(2+) and Na(+) cations along the exchange isotherm. The results further confirmed the influence of the solid/solution ratio on the degree of interlayer Ca(II)-for-Na(+) exchange, and specific selectivity coefficients for interlayer sites were determined. The effect of the solid/solution ratio was finally interpreted by the resulting local changes in the solution chemistry. We demonstrated that by accounting for the Donnan effect, the different data can be interpreted using a single selectivity coefficient. The obtained Kc constant was successfully applied to interpret existing hydrogeochemical data on a natural aquitard. This most likely represents a more constrained and valid approach for the modeling of reactive element transport in natural media than does the poorly defined Kd parameter.

Nature of the sites involved in the process of cesium desorption from vermiculite

Three particle size fractions of sodium-saturated vermiculite (10-20, 1-2 and 0.1-0.2 μm), differing only in their ratios of external-to-total sorption sites, were used to probe the nature of the sites involved in desorption of cesium ions. The sorption was investigated for initial aqueous concentrations of cesium ranging from 5.6×10(-4) to 1.3×10(-2) mol/L, and the cesium desorption was probed by exchange with ammonium ions. The results showed that (1) the amounts of desorbed cesium were strongly dependent on the particle size for a given initial aqueous cesium concentration and (2) the amounts of desorbed cations (Na(+) and Cs(+)) strongly decreased with increasing initial cesium aqueous concentration, irrespective of the particle size investigated. Quantitative analysis of these results suggested that cesium ions sorbed on external (edge+basal) sorption sites can be desorbed by ammonium ions. As a contrast, most of cesium ions sorbed on interlayer sites remain fixed due to the collapse of the structure under aqueous conditions. This study provides important information, such as the nature of the sites involved in the exchange process, when the thermodynamic formalism is considered to describe the ion-exchange process involving cesium and high-charge swelling clay minerals in polluted soil environments.

Interlayer structure model of tri-hydrated low-charge smectite by X-ray diffraction and Monte Carlo modeling in the Grand Canonical ensemble

Abstract The present study aims primarily at refining a structure model for interlayer cations and H2O molecules in tri-hydrated (3W) smectite (d001 = 18-19 Å). The <2 mm fraction of the SWy-2 source clay (low-charge montmorillonite) was saturated by Mg2+, Ca2+, Ba2+, or Na+ cations, before collection of X-ray diffraction (XRD) patterns at 98% relative humidity. Experimental d001 values derived for the essentially homogeneous 3W hydrates provided volume constraints for Grand Canonical Monte Carlo (GCMC) simulations. Computed atomic density distribution of interlayer species were used in turn to calculate XRD intensities of 00l reflections. The agreement between calculated and experimental 00l intensities allowed validating the GCMC results of both interlayer H2O content and distribution of interlayer species (cations and H2O molecules). Computed atomic density profiles do not correspond to the usual model of three discrete planes of H2O molecules but rather exhibit two sharp planes of H2O molecules wetting the clay surfaces (at ~2.7 Å from the clay layer surface). Additional H2O molecules belong to cation hydration shells or define a poorly organized ensemble filling internal voids. This alternative model suggests that the stability of the 3W hydrate results from the dual interaction of some H2O molecules with interlayer cation, through their second hydration shell, and with the 2:1 clay surface. Computed atomic density profiles were approximated to propose an interlayer structure model for 3W smectite. This simplified model includes two sets of two planes (symmetrical relative to the interlayer mid-plane) for H2O molecules and one set for interlayer cations. This model allows reproducing experimental XRD patterns for the different samples investigated and thus represents a valid set of parameters for routine quantitative analysis of XRD profiles in an effort to determine smectite reactivity close to water-saturated conditions. Implications of such studies are crucial to provide experimental constraints on the behavior of the main vector of element transfer under conditions common in surficial environments and prevailing in waste repositories. In addition, the present study provides an experimental validation of structure models derived from the widely used ClayFF model, and thus allows its use to predict the fate of water in clayey systems close to water-saturated conditions.

INVESTIGATION OF THE INTERLAYER ORGANIZATION OF WATER AND IONS IN SMECTITE FROM THE COMBINED USE OF DIFFRACTION EXPERIMENTS AND MOLECULAR SIMULATIONS. A REVIEW OF METHODOLOGY, APPLICATIONS, AND PERSPECTIVES

Investigation of the organization of interlayer water and cations in smectite is a permanent topic in clay science for environmental science, civil engineering, materials science, and industrial applications. Experimental X-ray (or neutron) diffraction methods and molecular simulations are key techniques to probe the organization of the smectite structure at a similar molecular length scale. The combination of both of these experimental and numerical methods represents a complementary approach to reveal the structural heterogeneity of real samples, design and model a wide range of smectite structures, and validate the simulation results through comparison with experimental data.This paper first revisits establishment of the original interlayer model as developed in the 1930s for the organization of water and ions in the smectite structure using X-ray diffraction (XRD) techniques. Then, based on a simplified approach, key theoretical tools are provided to calculate XRD pattern 00l reflections for a periodic smectite structure with a wide range of interlayer compositions and organizations using conventional spreadsheet software. In addition to educational purposes, this theoretical description is used to describe the principal parameters governing the positions and intensities of experimental XRD 00l reflections. This calculation toolbox is also used to determine better the layer-to-layer distances considered in molecular simulations and to validate these simulations through a detailed collation procedure using experimental data.Recent examples of the application of such a procedure to collate experimental diffraction data and molecular simulations are presented for the specific case of deciphering the molecular organization of interlayer water and cations in the different smectite hydrates (mono-, bi-, and tri-hydrated layers). The extension of this approach to the interlayer refinement of organo-clays is also detailed, and perspectives regarding the characterization of other lamellar compounds are discussed.

Evaluation of talc morphology using FTIR and H/D substitution

Abstract Deuteration (H/D substitution at 200ºC) was performed on powders of two ground talcs of different particle shapes (different basal/lateral surface ratios). Results indicate that the deuteration process is only efficient on lateral talc surfaces, and suggest that the hydrogens located in the hexagonal ring of the talc basal surfaces are not exchanged. The FTIR spectra collected from the two talc samples show that it is possible to discriminate between particles with the same specific surface area but with different basal/lateral surface ratios using the deuteration process.