Virginie Marry

Structure and dynamics of water at a clay surface from molecular dynamics simulation

We report a molecular dynamics study of the structure and dynamics of water at a clay surface. The negative charge of the surface and the presence of surface oxygen atoms perturbs water over two to three molecular layers, while the nature of the counterions (Na(+)or Cs(+)) has only a small effect. In the first molecular layer, approximately half of the water molecules are H-bonded to the surface. We also analyze the H-bond network between surface water molecules. The diffusion of water molecules along the surface is slowed down compared to the bulk case. As far as the orientational order and dynamics of the water dipole are concerned, only the component normal to the clay surface is perturbed. We investigate the surface H-bond formation and dissociation dynamics and their coupling to the release of molecules from the first molecular layer. We introduce a simple kinetic model in the spirit of Luzar and Chandler [Nature, 1996, 379, 55] to allow for a comparison with bulk water dynamics. This model semi-quantitatively reproduces the molecular simulation results and suggests that H-bond formation is faster with the surface than in the bulk, while H-bond dissociation is slower.

Microscopic simulation of structure and dynamics of water and counterions in a monohydrated montmorillonite

Models of swelling clays are studied by computer simulations (Monte Carlo and molecular dynamics). We focus on the comparison of structural and dynamic properties of two montmorillonites with different kinds of counterions Na+ and Cs+. The calculated values are compared with available experimental quantities such as interlayer spacing as a function of water content and diffusion coefficients of both water molecules and counterions in the monohydrated state. The results are consistent with experimental values and previous simulations. For the dynamics, the short time behavior of water as observed with quasielastic neutron scattering is in agreement with simulated one. For the ions, the experimental values are related to macroscopic long time motions and are much smaller than the short time values calculated from MD. Thus, the present study provides a detailed insight in the microscopic dynamics of ions related to the structure of the clay: it is shown that Cs+ diffuse faster than Na+ and that the arrangeme...

Water dynamics in hectorite clays: influence of temperature studied by coupling neutron spin echo and molecular dynamics.

Within the wider context of water behavior in soils, and with a particular emphasis on clays surrounding underground radioactive waste packages, we present here the translational dynamics of water in clays in low hydrated states as studied by coupling molecular dynamics (MD) simulations and quasielastic neutron scattering experiments by neutron spin echo (NSE). A natural montmorillonite clay of interest is modeled by a synthetic clay which allows us to understand the determining parameters from MD simulations by comparison with the experimental values. We focus on temperatures between 300 and 350 K, i.e., the range relevant to the highlighted application. The activation energy Ea experimentally determined is 6.6 kJ/mol higher than that for bulk water. Simulations are in good agreement with experiments for the relevant set of conditions, and they give more insight into the origin of the observed dynamics.

Diffusion coefficient and shear viscosity of rigid water models

We report the diffusion coefficient and viscosity of popular rigid water models: two non-polarizable ones (SPC/E with three sites, and TIP4P/2005 with four sites) and a polarizable one (Dang-Chang, four sites). We exploit the dependence of the diffusion coefficient on the system size (Yeh and Hummer 2004 J. Phys. Chem. B 108 15873) to obtain the size-independent value. This also provides an estimate of the viscosity of all water models, which we compare to the Green-Kubo result. In all cases, a good agreement is found. The TIP4P/2005 model is in better agreement with the experimental data for both diffusion and viscosity. The SPC/E and Dang-Chang models overestimate the diffusion coefficient and underestimate the viscosity.

Temperature effect in a montmorillonite clay at low hydration—microscopic simulation

The effect of temperature in the range 0–150°C was studied for homo-ionic montmorillonite clays with Na+ and Cs+ compensating ions in low hydration states. Monte Carlo and molecular dynamics simulations were employed to provide both static and dynamic information concerning the interlayer ions and water molecules, and emphasis was laid on the temperature activation of the diffusion coefficients. Principal structural changes were limited to the interlayer water phase. In the monohydrated systems, neither of the cations was seen to enter into the hexagonal cavities of the clay. Cs+ exhibited clear site-to-site diffusion between sites allowing coordination to six oxygen atoms of the clay sheets, this behaviour persisting to high temperatures. Preferential sites for the Na+ counterion were much less well-defined, even at low temperatures. The behaviour of the water phase in the monohydrated states was similar for the two ions. A rapid approach to bulk dynamics was seen in the transition from monohydrated to bihydrated Na-montmorillonite. A detailed quantitative comparison of the temperature activation of diffusion for a two-dimensional water phase and three-dimensional bulk water is presented for the first time.

Molecular simulation of aqueous solutions at clay surfaces

We report a molecular simulation study of aqueous solutions at montmorillonite clay surfaces. Unlike most previous studies, ours does not focus on the interlayer nanopores, but looks at both kinds of external surfaces of clay particles: basal surfaces along the clay layers, and lateral surfaces through which interlayer and larger interparticle pores are linked. We present results on structural, dynamic and thermodynamic properties and phenomena, including hydration complexes of ions, H bonding networks, modification of the water dynamics with respect to the bulk, and the role of water in the cation exchange between interlayer and interparticle pores.

Salt exclusion in charged porous media: a coarse-graining strategy in the case of montmorillonite clays

We study the exclusion of salt from charged porous media (Donnan effect), by using a coarse-grained approach. The porous medium is a lamellar system, namely a Montmorillonite clay, in contact with a reservoir, which contains an electrolyte solution. We develop a specific coarse-graining strategy to investigate the structural properties of this system. Molecular simulations are used to calibrate a mesoscopic model of the clay micropore in equilibrium with a reservoir. Brownian Dynamics simulations are then used to predict the structure of ions in the pore and the amount of NaCl salt entering the pore as a function of the pore size (the distance L between clay surfaces) and of the electrolyte concentration in the reservoir. These results are also compared to the predictions of a Density Functional Theory, which takes into account the excluded volumes of ions. We show that the calibration of the mesoscopic model is a key point and has a strong influence on the result. We observe that the salt exclusion increases when kappaL decreases (where kappa is the inverse of the Debye length) and that this effect is modulated by the correlations between ions. Two different regimes are revealed. At low concentrations in the reservoir, we observe a regime controlled by electrostatics: the Coulomb attraction between ions increases the amount of salt in the interlayer space. On the opposite, at high concentrations in the reservoir, the excluded volume effect dominates.

Equilibrium and electrokinetic phenomena in charged porous media from microscopic and mesoscopic models: electro-osmosis in montmorillonite

Ion concentration profiles and electro-osmosis have been theoretically investigated by complementary methods. The counterion profiles of these charged porous media are obtained from microscopic (molecular dynamics) and mesoscopic (Brownian dynamics and Poisson-Boltzmann equation) models. Thus, the limit of validity of the more macroscopic approaches has been quantified. The electro-osmotic flow is evaluated from a similar multi-scale description. A proper evaluation of the Kubo relations adapted to this dynamical process allowed the value of this transport quantity to be obtained from the microscopic model. The comparison with the classical Smoluchowski theory of this electrokinetic phenomenon shows that the latter cannot be applied unless slip hydrodynamic boundary conditions are taken into account. The corresponding slipping length is roughly equal to 6 Å.

Dynamics in Clays - Combining Neutron Scattering and Microscopic Simulation

Abstract Mobility of ions and water in clays is at the heart of their remarkable properties of water retention and ion-exchange. It has been addressed here using two microscopic techniques: neutron scattering and molecular dynamics simulations. Neutron scattering gives access exclusively to water dynamics in clays, due to the exceptional sensitivity of neutrons to H atoms. The data interpretation can be challenging, especially for natural clays such as montmorillonite, with inhomogeneous swelling characteristics. A great improvement is achieved with the use of synthetic materials, as demonstrated here on the case of synthetic (fluoro)hectorite. The standard analytical models for long-range diffusive motion, isotropic translation and its derivative, powder averaged two dimensional translation, have been used to interpret the neutron scattering data. They both agree on the order of magnitude for the diffusion coefficient of water in monohydrated and bihydrated clays, 10−10 m2s−1 and 10−9 m2s−1 respectively. While the two-dimensional nature of water diffusion in clays is seen clearly from molecular dynamics simulations, its signature in neutron scattering data is obscured by the powder-averaging of the signal. A novel method, based on a multi-resolution analysis of scattering functions from powder samples, allows never-the-less a clear determination of the dimensionality of water motion in the system. Extracting information on local water motion is difficult on the basis of neutron scattering data only. Various models for localised motion, rotation on a sphere or jump diffusion, have been proposed and used to interpret the observed neutron data, however their applicability is questionable in light of information from molecular dynamics simulations. Aside from aiding the interpretation of neutron scattering data, MD simulations are most valuable in providing information on the behaviour of ions in clays. MD estimates the interlayer ion coefficients as of the some order of magnitude as water, even if the details of ionic motion are strikingly different between the two ions considered here, Na+ and Cs+. Further, MD has also allowed to address the topic of ion exchange between clay interlayers and bulk aqueous solution. The microscopic picture of water and ion motion in clays, emerging from neutron scattering and MD simulations, should be treated as a building block of the overall modelling of macroscopic transport in clays, the ultimate property of interest for many clay applications.

Water diffusion in a synthetic hectorite by neutron scattering—beyond the isotropic translational model

A quasi-elastic neutron scattering study of water dynamics confined in a model clay system (a synthetic hectorite with Na+ compensating counterions) is presented. The neutron spin echo (NSE) data obtained at ambient temperature for the monohydrated and bihydrated hectorite powder are analysed with the help of models taking into account the anisotropy of the system. A powder-averaged two-dimensional (2D) diffusion model applied to the monohydrated state yields a 2D diffusion coefficient in the plane of the clay layers () of 2.8 × 10−10 m2 s−1, where an isotropic analysis of the data would have underestimated the diffusion coefficient by approximately 25%. In the case of the bihydrated state, the analysis by a model including a possible diffusion perpendicular to the clay layers () leads to a range of values for between 3.5 and 8.7 × 10−10 m2 s−1, depending on the ratio . We conclude that simple isotropic analysis can only give a rough estimate (order of magnitude) of the water dynamics in the medium.