Hostname: page-component-848d4c4894-nmvwc Total loading time: 0 Render date: 2024-06-23T06:02:24.234Z Has data issue: false hasContentIssue false

Comparative Modeling of Ions and Solvent Properties in Ca-Na Montmorillonite by Atomistic Simulations and Fluid Density Functional Theory

Published online by Cambridge University Press:  01 January 2024

Guomin Yang*
Affiliation:
Laboratory for Waste Management, Paul Scherrer Institute, 5232 Villigen, PSI, Switzerland
Nikolaos I. Prasianakis
Affiliation:
Laboratory for Waste Management, Paul Scherrer Institute, 5232 Villigen, PSI, Switzerland
Sergey V. Churakov
Affiliation:
Laboratory for Waste Management, Paul Scherrer Institute, 5232 Villigen, PSI, Switzerland Institute of Geological Sciences, University of Bern, Baltzerstrasse 1+3, Room 121, CH-3012 Bern, Switzerland
*
*E-mail address of corresponding author: guomin.yang@psi.ch

Abstract

Molecular dynamics (MD) simulations provide an accurate description of the mineral–fluid interface from the perspective of the atomistic level taking into account all atom interactions. This simulation approach is computationally expensive if applied to large molecular systems. Classical Fluid Density Functional Theory (f-DFT) delivers structural and thermodynamic information at comparatively small computational costs. Numerous applications of f-DFT for electrolytes neglect an explicit consideration of solvent. In this work, an unrestricted three-component model (3CM) of f-DFT was applied, which incorporates Lennard-Jones (LJ) attractions for the description of the short-range interactions of fluid–fluid and fluid–wall rather than the hard sphere repulsions, named DFT/LJ-3CM. The DFT/LJ-3CM model considers ions as charged LJ particles and treats solvent molecules as neutral LJ particles. To validate the performance of the DFT/LJ-3CM, the f-DFT calculations were compared with atomistic simulations for montmorillonite (Mnt) with various hydrated states in electrolyte solutions. This benchmarking was used to assess critically the advantages and limitations of the f-DFT model. The calibrated DFT/LJ-3CM model for Na and Ca Mnt was applied to calculate cation selectivity for the ion exchange equilibrium with effective ion radius and swelling behavior of Mnt. The predictions of the DFT/LJ-3CM model were found to be in good agreement with the atomistic simulations and experimental data under a wide range of conditions. At the same time, the DFT calculations were 3–4 orders of magnitude faster than conventional MD simulations. Thus, the DFT/LJ-3CM model can be a computationally effective alternative to atomistic simulation in providing structural and thermodynamic properties of fluid–clay mineral interfaces. The DFT/LJ-3CM model provides a robust approach, which can be used for upscaling in reactive transport simulators and modeling ion migration taking place under more complex thermo-chemo-hydro-mechanical conditions.

Type
Original Paper
Copyright
Copyright © Clay Minerals Society 2020

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Alt-Epping, P., Gimmi, T., Wersin, P., & Jenni, A. (2018). Incorporating electrical double layers into reactive-transport simulations of processes in clays by using the Nernst-Planck equation: A benchmark revisited. Applied Geochemistry, 89, 110.CrossRefGoogle Scholar
Appelo, C.A.J., & Postma, D. (2005). Ion exchange. Pp. 241309 in: Geochemistry, Groundwater and Pollution. Taylor & Francis. doi: https://doi.org/10.1201/9781439833544.ch6CrossRefGoogle Scholar
Åqvist, J. (1990). Ion-water interaction potentials derived from free energy perturbation simulations. The Journal of Physical Chemistry, 94, 80218024.CrossRefGoogle Scholar
Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A., & Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. The Journal of Chemical Physics, 81, 36843690.CrossRefGoogle Scholar
Berendsen, H. J. C., Grigera, J. R., & Straatsma, T. P. (1987). The missing term in effective pair potentials. The Journal of Physical Chemistry, 91, 62696271.CrossRefGoogle Scholar
Berendsen, H. J. C., van der Spoel, D., & van Drunen, R. (1995). GROMACS: a message-passing parallel molecular dynamics implementation. Computer Physics Communications, 91, 4356.CrossRefGoogle Scholar
Bloom, P.R. (2000). Soil pH and pH buffering. Pp. B333B352 in: Handbook of Soil Science. Section B Soil Chemistry (Sumner, M.E., edtor). CRC Press: Boca Raton, Florida, USA.Google Scholar
Boda, D., Henderson, D., Plaschko, P., & Fawcett, W. R. (2004). Monte Carlo and density functional theory study of the electrical double layer: the dependence of the charge/voltage relation on the diameter of the ions. Molecular Simulation, 30, 137.CrossRefGoogle Scholar
Boek, E. S., Coveney, P. V., & Skipper, N. T. (1995a). Molecular modeling of clay hydration: A study of hysteresis loops in the swelling curves of sodium montmorillonites. Langmuir, 11, 46294631.CrossRefGoogle Scholar
Boek, E. S., Coveney, P. V., & Skipper, N. T. (1995b). Monte Carlo Molecular modeling studies of hydrated Li-, Na-, and K-Smectites: Understanding the role of potassium as a clay swelling inhibitor. Journal of the American Chemical Society, 117, 1260812617.CrossRefGoogle Scholar
Botan, A., Rotenberg, B., Marry, V., Turq, P., & Noetinger, B. (2011). Hydrodynamics in Clay Nanopores. The Journal of Physical Chemistry C, 115, 1610916115.CrossRefGoogle Scholar
Botan, A., Marry, V., Rotenberg, B., Turq, P., & Noetinger, B. (2013). How electrostatics influences hydrodynamic boundary conditions: Poiseuille and electro-osmostic flows in clay nanopores. The Journal of Physical Chemistry C, 117, 978985.CrossRefGoogle Scholar
Bourg, I. C., & Sposito, G. (2011). Molecular dynamics simulations of the electrical double layer on smectite surfaces contacting concentrated mixed electrolyte (NaCl-CaCl2) solutions. Journal of Colloid and Interface Science, 360, 701715.CrossRefGoogle ScholarPubMed
Bourg, I. C., Bourg, A. C. M., & Sposito, G. (2003). Modeling diffusion and adsorption in compacted bentonite: a critical review. Journal of Contaminant Hydrology, 61, 293302.CrossRefGoogle ScholarPubMed
Chang, F. R. C., Skipper, N. T., & Sposito, G. (1998). Monte Carlo and molecular dynamics simulations of electrical double-layer structure in potassium-montmorillonite hydrates. Langmuir, 14, 12011207.CrossRefGoogle Scholar
Churakov, S. V. (2013). Mobility of Na and Cs on montmorillonite surface under partially saturated conditions. Environmental Science & Technology, 47, 98169823.CrossRefGoogle ScholarPubMed
Churakov, S. V., & Gimmi, T. (2011). Up-scaling of molecular diffusion coefficients in clays: A two-step approach. The Journal of Physical Chemistry C, 115, 67036714.CrossRefGoogle Scholar
Churakov, S. V., Gimmi, T., Unruh, T., Loon, L. R. V., & Juranyi, F. (2014). Resolving diffusion in clay minerals at different time scales: Combination of experimental and modeling approaches. Applied Clay Science, 96, 3644.CrossRefGoogle Scholar
Cygan, R. T., Liang, J. J., & Kalinichev, A. G. (2004). Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field. The Journal of Physical Chemistry, B, 108, 12551266.CrossRefGoogle Scholar
Darden, T., York, D., & Pederson, L. (1993). Particle mesh Ewald: An N-log(N) method for Ewald sums in large systems. The Journal of Chemical Physics, 98, 10089.CrossRefGoogle Scholar
Edil, T. B. (2003). A review of aqueous-phase VOC transport in modern landfill liners. Waste Management, 24, 561571.CrossRefGoogle Scholar
Eriksen, T. E., Jansson, M., & Molera, M. (1999). Sorption effects on cation diffusion in compacted bentonite. Engineering Geology, 54, 231236.CrossRefGoogle Scholar
Fenter, P., Cheng, L., Rihs, S., Machesky, M., Bedzyk, M. J., & Sturchio, N. C. (2000). Electrical double-layer structure at the rutile-water interface as observed in situ with small-period X-ray standing waves. Journal of Colloid and Interface Science, 225, 154165.CrossRefGoogle ScholarPubMed
Frink, L. J. D., & van Swol, F. (1998). Solvation forces between rough surfaces. The Journal of Chemical Physics, 108, 5588.CrossRefGoogle Scholar
Gaines, G. L. Jr., & Thomas, H. C. (1953). Adsorption studies on clay minerals. II. A formulation of the thermodynamics of exchange adsorption. The Journal of Chemical Physics, 21, 714.CrossRefGoogle Scholar
Gimmi, T., & Churakov, S. V. (2019). Water retention and diffusion in unsaturated clays: Connecting atomistic and pore scale simulations. Applied Clay Science, 175, 169183.CrossRefGoogle Scholar
Gimmi, T., & Kosakowski, G. (2011). How mobile are sorbed cations in clays and clay rocks? Environmental Science & Technology, 45, 14431449.CrossRefGoogle Scholar
Hansen, J.-P., & McDonald, I. R. (2006). Theory of Simple Liquids. London: Academic Press.Google Scholar
Henderson, D. (Ed.). (1992). Fundamentals of Inhomogeneous Fluids. New York: Marcel Dekker, Inc.Google Scholar
Hensen, E. J. M., & Smit, B. (2002). Why clays swell. The Journal of Physical Chemistry B, 106, 1266412667.CrossRefGoogle Scholar
Hoang, H., & Galliero, G. (2015). Couplings between swelling and shear in saturated slit nanopores: A molecular simulation study. Physical Review E, 91, 012401.CrossRefGoogle ScholarPubMed
Holmboe, M., & Bourg, I. C. (2014). Molecular dynamics simulations of water and sodium diffusion in smectite interlayer nanopores as a function of pore size and temperature. The Journal of Physical Chemistry C, 118, 10011013.CrossRefGoogle Scholar
Jeanmairet, G., Marry, V., Levesque, M., Rotenberg, B., & Borgis, D. (2014). Hydration of clays at the molecular scale: the promising perspective of classical density functional theory. Molecular Physics, 112, 13201329.CrossRefGoogle Scholar
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. W., & Klein, M. L. (1983). Comparison of simple potential functions for simulating liquid water. The Journal of Chemical Physics, 79, 926935.CrossRefGoogle Scholar
Joung, I. S., & Cheatham III, T. E. (2008). Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations. The Journal of Physical Chemistry B, 112, 90209041.CrossRefGoogle ScholarPubMed
Karaborni, S., Smit, B., Heidug, W., Urai, J., & Oort, E. V. (1996). The swelling of clays: Molecular simulations of the hydration of montmorillonite. Science, 271, 11021104.CrossRefGoogle Scholar
Kawamura, K., Ichikawa, Y., Nakano, M., Kitayama, K., & Kawamura, H. (1999). Swelling properties of smectite up to 90°C In situ X-ray diffraction experiments and molecular dynamic simulations. Engineering Geology, 54, 7579.CrossRefGoogle Scholar
Kosakowski, G., Churakov, S. V., & Thoenen, T. (2008). Diffusion of Na and Cs in montmorillonite. Clays and Clay Minerals, 56, 190206.CrossRefGoogle Scholar
Lamperski, S., & Zydor, A. (2007). Monte Carlo study of the electrode| solvent primitive model electrolyte interface. Electrochimica Acta, 52, 24292436.CrossRefGoogle Scholar
Lee, J. W., Nilson, R. H., Templeton, J. A., Griffiths, S.K., Kung, A., & Wong, B. M. (2012). Comparison of molecular dynamics with classical density functional and poisson-boltzmann theories of the electric double layer in nanochannels. Journal of Chemical Theory and Computation, 8, 20122022.CrossRefGoogle ScholarPubMed
Lee, J. W., Templeton, J. A., Mandadapu, K. K., & Zimmerman, J. A. (2013). Comparison of molecular and primitive solvent models for electrical double layers in nanochannels. Journal of Chemical Theory and Computation, 9, 30513061.CrossRefGoogle ScholarPubMed
Lian, C., Zhan, C., Jiang, D., Liu, H., & Wu, J. (2017). Capacitive energy extraction by few-layer graphene electrodes. The Journal of Physical Chemistry C, 121, 1401014018.CrossRefGoogle Scholar
Lo, I. M. C., Mak, R. K. M., & Lee, S. C. H. (1997). Modified clays for waste containment and pollutant attenuation. Journal of Environmental Engineering, 123, 2532.CrossRefGoogle Scholar
Magda, J. J., Tirrell, M., & Davis, H. T. (1985). Molecular dynamics of narrow, liquid-filled pores. The Journal of Chemical Physics, 83, 18881901.CrossRefGoogle Scholar
Malmberg, C.G., & Maryott, A. A. (1956). Dielectric constant of water from 0° to 100°C. Journal of Research of the National Institute of Standards and Technology, 56, 18.Google Scholar
Marry, V., Rotenberg, B., & Turq, P. (2008). Structure and dynamics of water at a clay surface from molecular dynamics simulation. Physical Chemistry Chemical Physics, 10, 4802.CrossRefGoogle Scholar
Marschall, P., Horseman, S., & Gimmi, T. (2005). Characterisation of gas transport properties of the Opalinus Clay, a potential host rock formation for radioactive waste disposal. i Technology, 60, 121139.Google Scholar
Melkior, T., Yahiaoui, S., Motellier, S., Thoby, D., & Tevissen, E. (2005). Cesium sorption and diffusion in Bure mudrock samples. Applied Clay Science, 29, 172186.CrossRefGoogle Scholar
Melkior, T., Yahiaoui, S., Thoby, D., Motellier, S., & Barthes, V. (2007). Diffusion coefficients of alkaline cations in Bure mudrock. Physics and Chemistry of the Earth, 32, 453462.CrossRefGoogle Scholar
Mooney, R. W., Keenan, A. G., & Wood, L. A. (1952). Adsorption of water vapor by montmorillonite. II. Effect of exchangeable ions and lattice swelling as measured by X-Ray diffraction. Journal of the American Chemical Society, 74, 13711374.CrossRefGoogle Scholar
Moucka, F., Svoboda, M., & Lisal, M. (2017). Modelling aqueous solubility of sodium chloride in clays at thermodynamic conditions of hydraulic fracturing by molecular simulations. Physical Chemistry Chemical Physics, 19, 1658616599.CrossRefGoogle ScholarPubMed
Obliger, A., Duvail, M., Jardat, M., Coelho, D., Bekri, S., & Rotenberg, B. (2013). Numerical homogenization of electrokinetic equations in porous media using lattice-Boltzmann simulations. Physical Review E, 88, 013019.CrossRefGoogle ScholarPubMed
Obliger, A., Jardat, M., Coelho, D., Bekri, S., & Rotenberg, B. (2014). Pore network model of electrokinetic transport through charged porous media. Physical Review E, 89, 043013.CrossRefGoogle ScholarPubMed
Oleksy, A., & Hansen, J. P. (2010). Wetting of a solid substrate by a “civilized” model of ionic solutions. The Journal of Chemical Physics, 132, 204702.CrossRefGoogle ScholarPubMed
Pagonabarraga, I., Rotenberg, B., & Frenkel, D. (2010). Recent advances in the modelling and simulation of electrokinetic effects: bridging the gap between atomistic and macroscopic descriptions. Physical Chemistry Chemical Physics, 12, 95669580.CrossRefGoogle ScholarPubMed
Prasianakis, N., Curti, I. E., Kosakowski, G., Poonoosamy, J., & Churakov, S. V. (2017). Deciphering pore-level precipitation mechanisms. Scientific Reports, 7, 13765.CrossRefGoogle ScholarPubMed
Price, D. J., & Brooks, C. L. (2004). A modified TIP3P water potential for simulation with Ewald summation. The Journal of Chemical Physics, 121, 1009610103.CrossRefGoogle ScholarPubMed
Qiao, R., & Aluru, N. R. (2003). Ion concentrations and velocity profiles in nanochannel electroosmotic flows. The Journal of Physical Chemistry, 118, 46924701.CrossRefGoogle Scholar
Rosenfeld, Y. (1989). Free-energy model for the inhomogeneous hardsphere fluid mixture and density-functional theory of freezing. Physical Review Letters, 63, 980.CrossRefGoogle ScholarPubMed
Rotenberg, B., Marry, V., Dufreche, J. F., Malikova, N., Giffaut, E., & Turq, P. (2007a). Modelling water and ion diffusion in clays: A multiscale approach. Comptes Rendus Chimie, 10, 1108.CrossRefGoogle Scholar
Rotenberg, B., Marry, V., Dufrêche, J.-F., Giffaut, E., & Turq, P. (2007b). A multiscale approach to ion diffusion in clays: Building a two-state diffusion-reaction scheme from microscopic dynamics. Journal of Colloid and Interface Science, 309, 289295.CrossRefGoogle ScholarPubMed
Rowe, R. K. (2005). Long-term performance of contaminant barrier systems. Geotechnique, 55, 631678.CrossRefGoogle Scholar
Schaettle, K., Pestana, L. R., Head-Gordon, T., & Lammers, L. N. (2018). A structural coarse-grained model for clays using simple iterative Boltzmann inversion. The Journal of Chemical Physics, 148, 222809.CrossRefGoogle ScholarPubMed
Schlegel, M. L., Nagy, K. L., Fenter, P., Cheng, L., Sturchio, N. C., & Jacobsen, S. D. (2006). Cation sorption on the muscovite (0 0 1) surface in chloride solutions using high-resolution X-ray reflectivity. Geochimica et Cosmochimica Acta, 70, 35493565.CrossRefGoogle Scholar
Schoonheydt, RA. & Johnston, C. (2006). Surface and interface chemistry of clay minerals. Pp. 139165 in: Handbook of Clay Science, Volume 1 (F. Bergaya, B.K.G. Theng, & G Lagaly editors). 1st edition. Elsevier Science.Google Scholar
Segad, M., Jonsson, B., Akesson, T., & Cabane, B. (2010). Ca/Na montmorillonite: Structure, forces and swelling properties. Langmuir, 26, 57825790.CrossRefGoogle ScholarPubMed
Silvestre-Alcantara, W., Kaja, M., Henderson, D., Lamperski, S., & Bhuiyan, L. B. (2016). Structure and capacitance of an electric double layer formed by fused dimer cations and monomer anions: a Monte Carlo simulation study. Molecular Physics, 114, 5360.CrossRefGoogle Scholar
SKB. Swedish Nuclear Fuel and Waste Management Co. (2006). Long-term safety for KBS-3 repositories at Forsmark and Laxemar – a first evaluat1ion. SKB Technical Report, TR-06-09.Google Scholar
Skipper, N. T., Chang, F. R. C., & Sposito, G. (1995a). Monte-Carlo simulation of interlayer molecular-structure in swelling clay-minerals. 1. Methodology. Clays and Clay Minerals, 43, 285293.CrossRefGoogle Scholar
Skipper, N. T., Chang, F. R. C., & Sposito, G. (1995b). Monte Carlo simulation of interlayer molecular structure in swelling clay minerals. 2. Monolayer hydrates. Clays and Clay Minerals, 43, 294303.CrossRefGoogle Scholar
Smith, D. E., Wang, Y., Chaturvedi, A., & Whitley, H. D. (2006). Molecular simulations of the pressure, temperature, and chemical potential dependencies of clay swelling. The Journal of Physical Chemistry B, 110,20046–20054.CrossRefGoogle Scholar
Song, J., & Zhang, D. (2013). Comprehensive review of caprocksealing mechanisms for geologic carbon sequestration. Environmental Science & Technology, 47, 922.CrossRefGoogle ScholarPubMed
Sposito, G. (1984). The Surface Chemistry of Soils. New York: Oxford University Press.Google Scholar
Sposito, G., Holtzclaw, K. M., Charlet, L., Jouany, C., & Page, A. L. (1983). Sodium-calcium and sodium-magnesium exchange on Wyoming bentonite in perchlorate and chloride background ionic media. Soil Science Society of America Journal, 47, 5156.CrossRefGoogle Scholar
Sun, L. L., Hirvi, J. T., Schatz, T., Kasa, S., & Pakkanen, T. A. (2015). Estimation of montmorillonite swelling pressure: A molecular dynamics approach. The Journal of Physical Chemistry C, 119, 1986319868.CrossRefGoogle Scholar
Tambach, T. J., Hensen, E. J. M., & Smit, B. (2004). Molecular simulations of swelling clay minerals. The Journal of Physical Chemistry B, 108, 75867596.CrossRefGoogle Scholar
Tambach, T. J., Bolhuis, P. G., Hensen, E. J. M., & Smit, B. (2006). Hysteresis in clay swelling induced by hydrogen bonding: Accurate prediction of swelling states. Langmuir, 22, 12231234.CrossRefGoogle ScholarPubMed
Tang, Z., Scriven, L. E., & Davis, H. T. (1991). Density-functional perturbation theory of inhomogeneous simple fluids. The Journal of Chemical Physics, 95, 2659.CrossRefGoogle Scholar
Tang, Z., Scriven, L. E., & Davis, H. T. (1992). A three-component model of the electrical double layer. The Journal of Chemical Physics, 97, 494.CrossRefGoogle Scholar
Torrie, G. M., & Valleau, J. P. (1980). Electrical double layers. I. Monte Carlo study of a uniformly charged surface. The Journal of Chemical Physics, 73, 5807.CrossRefGoogle Scholar
Tournassat, C., Chapron, Y., Leroy, P., Bizi, M., & Boulahya, F. (2009). Comparison of molecular dynamics simulations with triple layer and modified Gouy–Chapman models in a 0.1 M NaCl-montmorillonite system. Journal of Colloid and Interface Science, 339, 533541.CrossRefGoogle Scholar
Tyagi, M., Gimmi, T., & Churakov, S. V. (2013). Multi-scale microstructure generation strategy for up-scaling transport in clays. Advances in Water Resources, 59, 181195.CrossRefGoogle Scholar
van der Spoel, D., Lindahl, E., Hess, B., Groenhof, G., Mark, A. E., & Berendsen, H. J. C. (2005). GROMACS: fast, flexible, and free. Journal of Computational Chemistry, 26, 17011718.CrossRefGoogle ScholarPubMed
van Loon, L. R., & Eikenberg, J. (2005). A high-resolution abrasive method for determining diffusion profiles of sorbing radionuclides in dense argillaceous rocks. Applied Radiation and Isotopes, 63, 1121.CrossRefGoogle ScholarPubMed
van Schaik, J. C., Kemper, W. D., & Olsen, S. R. (1966). Contribution of adsorbed cations to diffusion in clay-water systems. Soil Science Society of America Journal, 30, 1722.CrossRefGoogle Scholar
Waisman, E., & Lebowitz, J. L. (1972). Mean spherical model integral equation for charged hard spheres I. Method of solution. The Journal of Chemical Physics, 56, 30863093.CrossRefGoogle Scholar
Yang, G., & Liu, L. (2015). A systematic comparison of different approaches of density functional theory for the study of electrical double layers. The Journal of Chemical Physics, 142, 194110.CrossRefGoogle Scholar
Yang, G., Neretnieks, I., Moreno, L., & Wold, S. (2016). Density functional theory of electrolyte solutions in slit-like nanopores II. Applications to forces and ion exchange. Applied Clay Science, 132–133, 561570.CrossRefGoogle Scholar
Yang, G., Neretnieks, I., Moreno, L., & Wold, S. (2017a). Density functional theory of electrolyte solutions in slit-like nanopores I. The RFD/WCA approach extended to non-restricted primitive model. Applied Clay Science, 135, 526531.CrossRefGoogle Scholar
Yang, G., Neretnieks, I., & Holmboe, M. (2017b). Atomistic simulations of cation hydration in sodium and calcium montmorillonite nanopores. The Journal of Chemical Physics, 147, 084705.CrossRefGoogle ScholarPubMed
Yang, Y., Patel, R. A., Churakov, S. V., Prasianakis, N. I., Kosakowski, G., & Wang, M. (2019). Multiscale modeling of ion diffusion in cement paste: electrical double layer effects. Cement and Concrete Composites, 96, 5565.CrossRefGoogle Scholar
Young, D. A., & Smith, D. E. (2000). Simulations of clay mineral swelling and hydration: Dependence upon interlayer ion size and charge. The Journal of Physical Chemistry B, 104, 91639170.CrossRefGoogle Scholar