Skip to Main content Skip to Navigation
Poster communications

Effects of surface cations on the structure and dynamics of the hydrogen-bonding network at the illite-water interface: A molecular dynamics simulation study

L. Narasimhan 1, * Andrey G. Kalinichev 1, *
* Corresponding author
1 SUBATECH - Radiochimie
SUBATECH - Laboratoire SUBATECH Nantes
Abstract : Safe and sustainable management of nuclear energy poses major scientific and engineering challenges, one of which is the necessity to make the environmental impacts of the long-term nuclear waste storage as small as possible. This requires detailed understanding and prediction of the behaviour of radionuclides and their migration and retention properties in the geological formations of nuclear waste repositories. The Callovo-Oxfordien rock formation of the French nuclear repository site is mainly composed of clay minerals (illite, smectite and interstratified illite/smectite), quartz, calcite, with some non-negligible amount of organic matter. The adsorption of water can change the properties of mineral surfaces, including protonation state, surface charge, structure, and reactivity [1]. Similarly, the properties of interfacial water are strongly affected by the mineral substrate structure and composition. Recent advances in experimental techniques such as FTIR [2], ellipsometry [3], synchrotron X-ray scattering [4], sum-frequency vibrational spectroscopy [5] are capable of probing the properties of mineral-water interfaces at different levels of hydration. However, the surface-specific results of these experiments are often difficult to quantitatively interpret without having a reliable molecular scale picture of the underlying physical and chemical processes. Molecular computer simulations have become one of the most important tools in the study of such interfacial systems and phenomena by providing invaluable atomistic information on the underlying chemical and physical processes. The present study is aimed at investigating the structural and dynamics effects of three different cations (K+, NH4+, and H3O+) exchanged at the hydrated surface of muscovite mica, which is taken here as a model illite. Molecular dynamics computer simulations were performed using the CLAYFF force field [6] to investigate the important differences of the H-bonding configurations formed by the sorbed species, including H2O, H3O+, and NH4+, in contrast to the behavior of spherical metal ions, such as K+. At the muscovite (001) surface, H2O can donate 2 H-bonds (to other H2O and/or to the surface O atoms) and accept 2 H-bonds (from other H2O), but it can also partially replace surface K+, because the hydrogens of H2O bear some positive charge. This behavior was observed in previous MD simulations of the mica surface [7]. Such surface-adsorbed H2O molecules have their negatively charged oxygen atoms exposed to the fluid phase and accessible for either H-bond acceptance from other H2Os or hydration of metal cations in outer-sphere coordination. For surface H3O+, the charges on the hydrogens are slightly higher than those of H2O, but the oxygen atom of hydronium is now almost hydrophobic, and cannot participate in a H-bond network (e.g., [8]). In contrast, NH4+ can equally well donate H-bonds to the surface O atoms and to the neighboring H2O molecules, but it cannot participate in the hydration shell of a displaced metal cation. Thus, three similar species (H2O - two HB donors and 2 HB acceptors; H3O+ - 3 HB donors and no acceptors; NH4+ - 4 HB donors, no acceptors) can provide for three greatly different structural, energetic, and dynamical situations at the muscovite-water interface. Since the hydrogen-bonding network in any aqueous media provides a natural mechanism of forming low-barrier reaction paths for proton transfer in such systems, it is also an important phenomenon controlling the surface reactivity under various pH conditions. In addition, a detailed study of the structural characteristics of surface-adsorbed NH4+ provides a way for better understanding of the mechanisms of adsorption for organic molecules having amino-groups in their structure, which is quite common for natural organic matter (e.g., [9]). Each of the three systems was simulated at 7 different hydration states providing information on the structure and dynamics of the adsorbed water film in a wide range of relative humidity conditions. The atomic density profiles of water show significant layering at all hydration levels and the layering strongly depends upon the nature of the ionic species present on the surface. Our studies support the fact that the H3O+ ion is less strongly bound when compared to K+ on the muscovite surface as observed in earlier studies [7]. At muscovite surfaces, both NH4+ and H3O+ cations establish strong hydrogen bonds with the surface bridging oxygen atoms and also with the neighbouring H2O molecules at all hydration levels. However, we observed that the interactions are different for both species at low hydration levels (<< molecular monolayer). At low hydration levels, H3O+ prefers to strongly bind as 3-cordinated species to the surface than with the neighbouring waters molecules. However, as the hydration levels increase, H3O+ binds as 2-cordinated species with the surface as is indicated by hydrogen bonding analysis (Figure 1). In contrast, irrespective of the hydration levels, NH4+ ion strongly interacts with the surface as 3-cordinated species because of its tetrahedral geometry. At the same time, we observe from hydrogen bond analysis that the hydrogen bonding network of water has been strongly influenced by the nature of the surface cations present at the mineral-water interface. The dynamics of water molecules were examined by self-diffusion coefficients from the mean square displacement of water oxygen. The diffusion mechanism is similar for K+ and NH4+ but was different for H3O+, in particular at the low hydration states. Furthermore, the spatial and orientation distributions of H2O and ions at the muscovite-water surface are analyzed in quantitative detail. All the simulation results are compared with available experimental data and the results of previous molecular simulations to provide reliable molecular view of the ions and water at the muscovite surface. References [1] Henderson, M. A. Surf. Sci. Rep. 2002, 46, 5-308. [2] Cantrell and G. E. Ewing. J. Phys. Chem. B, 2001, 105, 5434-5439. [3] Beaglehole, D and Christenson, H. K. J. Phys. Chem, 1992, 96, 3395-3403. [4] Fenter, P. and Sturchio, N. C. Progress in Surface Science, 2004, 77, 171-258. [5] Shen, Y. R. and Ostroverkhov, V. Chem. Rev., 2006, 106, 1140-1154. [6] Cygan, R.T., Liang, J.J., and Kalinichev, A.G. J. Phys. Chem. B, 2004, 108, 1255-1266. [7] Wang, J., Kalinichev, A., Kirkpatrick, R., Cygan, R. J. Phys. Chem B., 2005, 109, 15893-15905. [8] Petersen, P.B. and Saykally, R.J. J. Phys. Chem. B, 2005, 109, 7976-7980. [9] Leenheer, J.A. Annals of Environmental Science, 2009, 3, 1-130.
Complete list of metadatas
Contributor : Andrey Kalinichev <>
Submitted on : Friday, December 28, 2012 - 8:09:54 PM
Last modification on : Friday, June 7, 2019 - 6:10:03 PM


  • HAL Id : in2p3-00769153, version 1


L. Narasimhan, Andrey G. Kalinichev. Effects of surface cations on the structure and dynamics of the hydrogen-bonding network at the illite-water interface: A molecular dynamics simulation study. 5th International meeting "Clays in Natural and Engineered Barriers for Radioactive Waste Confinement", Oct 2012, Montpellier, France. ⟨in2p3-00769153⟩



Record views