Ion exchange selectivity in clay is controlled by nanoscale chemical-mechanical coupling

Uncovering the predictive pathways of lithium and sodium interchange in layered oxides

Ion exchange is a powerful method to access metastable materials with advanced functionalities for energy storage applications. However, high concentrations and unfavourably large excesses of lithium are always used for synthesizing lithium cathodes from parent sodium material, and the reaction pathways remain elusive. Here, using layered oxides as model materials, we demonstrate that vacancy level and its corresponding lithium preference are critical in determining the accessible and inaccessible ion exchange pathways. Taking advantage of the strong lithium preference at the right vacancy level, we establish predictive compositional and structural evolution at extremely dilute and low excess lithium based on the phase equilibrium between Li0.94CoO2 and Na0.48CoO2. Such phase separation behaviour is general in both surface reaction-limited and diffusion-limited exchange conditions and is accomplished with the charge redistribution on transition metals. Guided by this understanding, we demonstrate the synthesis of NayCoO2 from the parent LixCoO2 and the synthesis of Li0.94CoO2 from NayCoO2 at 1-1,000 Li/Na (molar ratio) with an electrochemical assisted ion exchange method by mitigating the kinetic barriers. Our study opens new opportunities for ion exchange in predictive synthesis and separation applications.

Electrostatic coupling and water bridging in adsorption hierarchy of biomolecules at water-clay interfaces

Clay minerals are implicated in the retention of biomolecules within organic matter in many soil environments. Spectroscopic studies have proposed several mechanisms for biomolecule adsorption on clays. Here, we employ molecular dynamics simulations to investigate these mechanisms in hydrated adsorbate conformations of montmorillonite, a smectite-type clay, with ten biomolecules of varying chemistry and structure, including sugars related to cellulose and hemicellulose, lignin-related phenolic acid, and amino acids with different functional groups. Our molecular modeling captures biomolecule-clay and biomolecule-biomolecule interactions that dictate selectivity and competition in adsorption retention and interlayer nanopore trapping, which we determine experimentally by NMR and X-ray diffraction, respectively. Specific adsorbate structures are important in facilitating the electrostatic attraction and Van der Waals energies underlying the hierarchy in biomolecule adsorption. Stabilized by a network of direct and water-bridged hydrogen bonds, favorable electrostatic interactions drive this hierarchy whereby amino acids with positively charged side chains are preferentially adsorbed on the negatively charged clay surface compared to the sugars and carboxylate-rich aromatics and amino acids. With divalent metal cations, our model adsorbate conformations illustrate hydrated metal cation bridging of carboxylate-containing biomolecules to the clay surface, thus explaining divalent cation-promoted adsorption from our experimental data. Adsorption experiments with a mixture of biomolecules reveal selective inhibition in biomolecule adsorption, which our molecular modeling attributes to electrostatic biomolecule-biomolecule pairing that is more energetically favorable than the biomolecule-clay complex. In sum, our findings highlight chemical and structural features that can inform hypotheses for predicting biomolecule adsorption at water-clay interfaces.

Ammonium concentration in stream sediments resulting from decades of discharge from a wastewater treatment plant

A study of ammonium pollution in the sediments of a stream that receives wastewater treatment plant (WWTP) discharge has been carried out. It is urgently necessary to find environmental indicators that can help prevent and detect potential contamination of water, as water is an increasingly scarce resource. To understand the behaviour of ammonium ions introduced by a historical (50-year) contamination process, vertical boreholes were drilled in the stream banks to depths between 30 and 120 cm. Moisture, pH, ammonium (soluble and exchangeable), and clay fraction content were analysed. The variation profile of these parameters was evaluated as a function of depth to determine factors related to the distribution of ammonium in several locations along the stream banks. The ammonium concentration was asymmetrically distributed among samples collected in near-surface locations, with ammonium concentrations between 0.3048 mmol/kg soil and 0.0007 mmol/kg soil. Ammonium was typically concentrated at sediment depths of 30-40 cm, which also exhibited the highest clay fraction content. High positive correlations were detected (r > 0.8; p < 0.0001) among the different ammonium variables (exchanged and dissolved species). No contamination effect was observed below 60-70 cm depth, which was due to ammonium retention in a natural barrier layer of clayey sediment. The clays in our study area (previously identified as smectite, a 2:1 sheet silicate) were able to control the contamination by retaining ammonium in the interlayers, which retarded nitrification. It is suggested that clay could serve as a geo-indicator of ammonium pollution evolution.

Smectite phase separation is driven by hydration-mediated interfacial charge

Smectite clay minerals have an outsize impact on the response of clay-rich media to common stimuli, such as hydration and ion exchange, motivating extensive effort to understand behaviors resulting from these processes such as swelling and exfoliation. Smectites are common and historic systems for investigating colloidal and interfacial phenomena, with two swelling regimes commonly identified across myriad clays: osmotic swelling at high water activity and crystalline swelling at low water activity. However, no current swelling model seamlessly spans the full ranges of water, salt and clay content encountered in natural or engineered settings. Here, we show that structures previously rationalized as either osmotic or crystalline coexist as a rich array of distinct colloidal phases that differ by water content, layer stacking thickness, and curvature. We present an analytical model for intermolecular potentials among water, salt and clay in both mono- and divalent electrolytes that predicts swelling pressures across high and low water activities. Our results indicate that all clay swelling is osmotic swelling, but that the osmotic pressure of charged mineral interfaces becomes attractive and dominates that of the electrolyte at high clay activities. Global energy minima are often not reached on experimental timescales due to many local energy minima that promote long-lived intermediate states with vast differences in clay, ion, and water mobilities, leading to hyperdiffusive layer dynamics driven by variable hydration-mediated interfacial charge. Teaser Distinct colloidal phases of swelling clays emerge via ion (de)hydration at mineral interfaces that drives hyperdiffusive layer dynamics as metastable smectites approach equilibrium.

Interaction between Hydrated Smectite Clay Particles as a Function of Salinity (0-1 M) and Counterion Type (Na, K, Ca)

Swelling clay minerals control the hydrologic and mechanical properties of many soils, sediments, and sedimentary rocks. This important and well-known phenomenon remains challenging to predict because it emerges from complex multiscale couplings between aqueous chemistry and colloidal interaction mechanics in nanoporous clay assemblages, for which predictive models remain elusive. In particular, the predominant theory of colloidal interactions across fluid films, the widely used Derjaguin-Landau-Verwey-Overbeek model, fails to predict the ubiquitous existence of stable swelling states at interparticle distances below 3 nm that are stabilized by specific inter-atomic interactions in overlapping electrical double layers between the charged clay surfaces. Atomistic simulations have the potential to generate detailed insights into the mechanisms of these interactions. Recently, we developed a metadynamics-based molecular dynamics simulation methodology that can predict the free energy of interaction between parallel smectite clay particles in a wide range of interparticle distances (from 0.3 to 3 nm) and salinities (from 0.0 to 1.0 M NaCl). Here, we extend this work by characterizing the sensitivity of interparticle interactions to counterion type (Na, K, Ca). We establish a detailed picture of the free energy of interaction of parallel clay particles across water films as the sum of five interaction mechanisms with different sensitivities to salinity, counterion type, and interparticle distance.

Critical role of water structure around interlayer ions for ion storage in layered double hydroxides

Water-containing layered materials have found various applications such as water purification and energy storage. The highly structured water molecules around ions under the confinement between the layers determine the ion storage ability. Yet, the relationship between the configuration of interlayer ions and water structure in high ion storage layered materials is elusive. Herein, using layered double hydroxides, we demonstrate that the water structure is sensitive to the filling density of ions in the interlayer space and governs the ion storage. For ion storage of dilute nitrate ions, a 24% decrease in the filling density increases the nitrate storage capacity by 300%. Quartz crystal microbalance with dissipation monitoring studies, combined with multimodal ex situ experiments and theoretical calculations, reveal that the decreasing filling density effectively facilitates the 2D hydrogen-bond networking structure in water around interlayer nitrate ions along with minimal change in the layered structure, leading to the high storage capacity.

Cation hydration by confined water and framework-atoms have crucial role on thermodynamics of clay swelling

The swelling capacity and stability of clay play a crucial role in various areas ranging from cosmetics to oil extraction; hence change in their swelling behaviour after cation exchange with the surrounding medium is important for their efficient utilisation. Here we focus on understanding the role of different hydration properties of cation on the thermodynamics of clay swelling by water adsorption. We have used mica as the reference clay, Na\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+, Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+, and H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ ions as the interstitial cations, and performed grand canonical Monte Carlo simulations of water adsorption in mica pores (of widths \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d = 4-40$$\end{document}d=4-40 Å). The disjoining pressure (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Pi$$\end{document}Π), swelling free energy (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \Omega ^{ex}$$\end{document}ΔΩex), and several structural properties of confined water and ions were calculated to perform a thermodynamic analysis of the system. We expected higher water density in H-mica pores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho_{ \hbox{H}}$$\end{document}ρH) due to the smaller size of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {H}^+$$\end{document}H+ ions having higher hydration energy. However, the counter-intuitive trend of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho _{\hbox{Li}}> \rho _{\hbox{Na}} > \rho_b$$\end{document}ρLi>ρNa>ρb (bulk density) \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$> \rho_{\hbox{H}}$$\end{document}>ρH was observed due to adsorption energy, where the interaction of water with mica framework atoms was also found to be significant. All three mica systems exhibited oscillatory behaviour in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Pi$$\end{document}Π and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \Omega ^{ex}$$\end{document}ΔΩex profiles, diminishing to zero for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d \ge 11$$\end{document}d≥11 Å. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \Omega ^{ex}$$\end{document}ΔΩex for Na-mica is characterised by global minima at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=6 {\hbox {\AA}}$$\end{document}d=6Å   corresponding to crystalline swelling with significant and multiple barriers for crystalline swelling to osmotic swelling (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d > 12$$\end{document}d>12 Å). A shift in the location of global minima of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \Omega ^{ex}$$\end{document}ΔΩex towards the higher d values and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \Omega ^{ex}$$\end{document}ΔΩex becoming more repulsive is observed in the increasing order of hydration energy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Na}^+$$\end{document}Na+, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {Li}^+$$\end{document}Li+, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {H}^+$$\end{document}H+ ions. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta \Omega ^{ex} > 0$$\end{document}ΔΩex>0 for all d in the H-mica system thus favours osmotic swelling. We found that the Na\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ ions hydrate more surface oxygens, act as anchors, and hold the mica pore together (at smaller d), by sharing hydrating water with ions of the opposite side, forming an electrostatically connected mica-Na-water-Na-mica bridge. The Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ ions do hydrate surface oxygen atoms, albeit in lesser numbers, and sharing of hydration shell with nearby Li\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ ions is also minimum. Hydration by surface atoms and water sharing, both, are minimum in the H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^+$$\end{document}+ ion case, as they are mostly present in the center of the pore as diffusive ions, thus exerting a consistent osmotic pressure on the mica frameworks, favouring swelling.

Bioaccessibility and Toxicity Assessment of Polycyclic Aromatic Hydrocarbons in Two Contaminated Sites

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous soil contaminants, and their bioaccessibility determines their environmental risks in contaminated land. In the present study, the residual concentrations of PAHs in the soils of two industrial sites were determined, and their bioaccessibility was estimated by the hydroxypropyl-β-cyclodextrin extraction (HPCD) extraction method. The results showed heavy PAH contamination at both site S1 (0.38-3342.5 mg kg^-1) and site S2 (0.2-138.18 mg kg^-1), of which high molecular weight (HMW) PAHs (4-, 5-, and 6-ring compounds) accounted for approximately 80%. The average bioaccessibility of PAHs at sites S1 and S2 was 52.02% and 29.28%, respectively. The bioaccessibility of certain PAH compounds decreased with increasing ring number of the molecule. Lower PAH bioaccessibility was detected in loamy and silty soil textures than in sandy soil. Moreover, among the soil properties, the dissolved organic matter, total organic carbon, total potassium, and total manganese concentrations had significant effects on the bioaccessibility of PAHs. The toxicity analysis showed that the composition and bioaccessibility of PAHs could affect their potential toxicity in soil. We suggest that bioaccessibility should be taken into consideration when assessing the toxicity of PAHs in soil, and more attention should be given to low-ring PAHs with high bioaccessibility.

Surface acidity modulates the peroxidase-like activity of nanoclay

A novel surface acidity modulation strategy allows us to obtain modified nanoclay with specific peroxidase (POD)-like catalytic activity. Fe³⁺ exchange could increase the surface acidity of modified montmorillonite (MMT), resulting in a significant enhancement of its POD-like activity. We proposed that the POD-like catalytic reaction followed the electron transfer pathway and ping-pong mechanism. Correspondingly the constructed colorimetric sensor for H₂O₂ exhibited high sensitivity and specificity.

Structural Alteration, Hydration Stability, Heavy Metal Removal Efficiency, and Montmorillonite Porosity Fate by Coupling the Soil Solution pH and a Thermal Gradient

Clay minerals are considered as a promising material in the context of geological barrier for the confinement of radioactive and industrial waste. Understanding the usefulness of the smectite mineral, in this approach, remains insufficient. The deep investigation about structural response/changes, hydrate stability, cation exchange process, permeability, and heavy metal/radionuclide adsorption/removal efficiency under external constraints is needed. To explore the structural alteration, the hydration stability, and the evolution of montmorillonite porosity under a first order of external applied constraints coupling, a reference Na-rich montmorillonite specimen is used as a starting material, and three exchangeable heavy metal cations (Ba2+, Cu2+, and Co2+) have been selected. The applied constraint coupling is realized at laboratory scale and assured by simultaneously varying of the soil solution pH and the thermal gradient. The evaluation of the mineral fraction response is established by correlation of quantitative XRD analysis results, thermal analysis, and porosity measurements. The quantitative XRD analysis allows rebuilding of the theoretical model describing the interlamellar space (IS) configurations/damages, structural heterogeneity degrees, and hydrous stability. Obtained results show a dominant interstratified hydration character, for all studied complexes, attributed to a new IS organization versus the applied constraint strength. Furthermore, all samples exhibit a heterogeneous hydration behavior traduced by the coexistence of different layers of type population within the crystallite. Additionally, the theoretical XRD profile decomposition allowed us to prove link between the domination of the segregated stacking layers mode against the constraint strength. Thermal analysis allowed us to develop theoretical models describing the decrease of the water molecule amounts as a function of the increase in temperature and soil solution pH. Moreover, a specific hydration footprint response and an interstratification mapping are assigned for each corresponding stress degrees. The evolution of montmorillonite porosity is determined by adsorption measurement, based on Brunauer, Emmett, and Teller (BET) and Barrett-Joyner-Halenda (BJH) pore size distribution analyses which confirms for all samples, the exfoliation process, and the mesopore diameter rise by increasing the constraint intensity.

Toward a mechanistic understanding of cesium adsorption to todorokite: A molecular dynamics simulation study

A mechanistic understanding of cesium (Cs) adsorption to soil mineral phases is essential for effective mitigation of Cs mobility in the subsurface environment. Todorokite, a common tunnel-structured manganese oxide in soil, exhibits sorption capacity for Cs comparable to the capacities of clay minerals. However, the adsorption sites and molecular species of Cs⁺ adsorbed to todorokite remain uncertain in comparison with those of clay minerals. In this study, we explored adsorption of Cs⁺ to hydrated todorokite surfaces via atomistic molecular dynamics (MD) simulations. We performed the first MD simulations based on atomic pair potentials for Mn-oxide edge surfaces interfaced with an aqueous solution. MD simulations predicted that Cs⁺ forms only inner-sphere (IS) complexes within todorokite tunnels; however, Cs⁺ forms both IS and outer-sphere (OS) complexes at the external (010) and (100)/(001) external surfaces. On the (010) surface, the positions between IS and OS complexes of Cs⁺ were interchangeable during MD simulations. Detailed molecular structures of IS and OS Cs⁺ surface complexes are compared to those of Cs⁺ in an aqueous solution. The current MD simulation results can be used as an atomistic structural proxy for spectroscopic analysis of adsorbed metal speciation and surface complexation modeling of metal adsorption to Mn oxides.

Ion complexation waves emerge at the curved interfaces of layered minerals

Visualizing hydrated interfaces is of widespread interest across the physical sciences and is a particularly acute need for layered minerals, whose properties are governed by the structure of the electric double layer (EDL) where mineral and solution meet. Here, we show that cryo electron microscopy and tomography enable direct imaging of the EDL at montmorillonite interfaces in monovalent electrolytes with ångstrom resolution over micron length scales. A learning-based multiple-scattering reconstruction method for cryo electron tomography reveals ions bound asymmetrically on opposite sides of curved, exfoliated layers. We observe conserved ion-density asymmetry across stacks of interacting layers in cryo electron microscopy that is associated with configurations of inner- and outer-sphere ion-water-mineral complexes that we term complexation waves. Coherent X-ray scattering confirms that complexation waves propagate at room-temperature via a competition between ion dehydration and charge interactions that are coupled across opposing sides of a layer, driving dynamic transitions between stacked and aggregated states via layer exfoliation.

The structure of hydrated interfaces is essential for understanding of geochemical processes and behavior of layered minerals. The authors show that waves of hydrated ions emerge at curved aqueous interfaces and couple mineral deformation to the chemistry of the solution.

Control of Structural Hydrophobicity and Cation Solvation on Interlayer Water Transport during Clay Dehydration

Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1-2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. This mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.

Thermodynamics of ion exchange coupled with swelling reactions in hydrated clay minerals

Crystalline hydrates of swelling clay minerals (smectites) exhibit a strong coupling between their ion exchange and hydration/dehydration reactions. The uptake or removal of water from smectite interlayers as a result of a change in the environmental conditions also leads to the partitioning of cations. Three factors, the solid ion composition, the solid basal spacing/water content, and the aqueous solution composition, are all implicated in controlling the thermodynamics of ion exchange. However, conventional approaches to measuring the exchange free energy cannot separate the influence of each of these individual factors. Here, we explore the energetics of the swelling and ion exchange reactions in montmorillonite using a potential of mean force approach and the thermodynamic integration method within molecular simulations. We investigate the influence of solution and clay composition on the spontaneity of the reactions, focusing on the 2 water-layer hydration state. The swelling simulations provide the equilibrium water content, interlayer water structure, and basal spacings, while thermodynamic integration of sodium-potassium exchange in the aqueous solution and solid phase are combined to calculate ion exchange free energies as a function of solution composition. Results confirm the tendency of the clay to collapse to lower hydration states as the concentration of the solution increases. Changes to the equilibrium water content, even at fixed hydration states, and the composition of the mixed electrolyte solution play a critical role in driving ion exchange and the selectivities of the clay to the exchanged cation, while the composition of the solid phase is shown to be insignificant. These findings underscore the extreme sensitivity of clay swelling and ion exchange thermodynamics to small (tenths of an Angstrom) deviations in layer spacing.

Modulation of soft glassy dynamics in aqueous suspensions of an anisotropic charged swelling clay through pH adjustment

HYPOTHESIS

Sodium-montmorillonite (Na-Mt) particles are geometrically anisometric that carry a pH dependent anisotropic surface charge. Therefore, it should be possible to manipulate the particle-particle interaction of colloidal range Na-Mt suspensions through pH changes which in turn should alter the soft glassy dynamics of Na-Mt suspensions.

EXPERIMENTS

Rheological experiments were used to probe the impact of pH mediated colloidal particle-particle interaction on the physical aging, linear viscoelastic response, and yield stress behavior of Na-Mt suspension.

FINDINGS

The temporal evolution of the storage modulus (G') was stronger in the acid regime (pH < 9.5) than the base (pH ≥ 9.5) pH regime. Horizontal shifting of the aging curves in the acid and base regimes led to aging time-H⁺ concentration and aging time-OH^- concentration superposition. An aging time-Na-Mt concentration superposition was also observed in both pH regimes. The critical stress associated with the viscosity bifurcation behavior increased linearly with G' but with different slopes for acid and base regime. We propose that positively charged patches on the Na-Mt particle edge merge with the characteristic surface as a function of H⁺ ions in the system. This leads to a strongly associated microstructure at low pH and a relatively weak but associated microstructure at natural pH, hence confirming the hypothesis.

Molecular Determination of Organic Adsorption Sites on Smectite during Fe Redox Processes Using ToF-SIMS Analysis

Turnover of soil organic carbon (SOC) is strongly affected by a balance between mineral protection and microbial degradation. However, the mechanisms controlling the heterogeneous and preferential adsorption of different types of SOC remain elusive. In this work, the heterogeneous adsorption of humic substances (HSs) and microbial carbon (MC) on a clay mineral (nontronite NAu-2) during microbial-mediated Fe redox cycling was determined using time-of-flight secondary ion mass spectrometry (ToF-SIMS). The results revealed that HSs pre-adsorbed on NAu-2 would partially inhibit structural modification of NAu-2 by microbial Fe(III) reduction, thus retarding the subsequent adsorption of MC. In contrast, NAu-2 without precoated HSs adsorbed a significant amount of MC from microbial polysaccharides as a result of Fe(III) reduction. This was attributed to the deposition of a thin Al-rich layer on the clay surface, which provided active sites for MC adsorption. This study provides direct and detailed molecular evidence for the first time to explain the preferential adsorption of MC over HSs on the surface of clay minerals in iron redox processes, which could be critical for the preservation of MC in soil. The results also indicate that ToF-SIMS is a unique tool for understanding complex organic-mineral-microbe interactions.

Role of cation size on swelling pressure and free energy of mica pores

The ion exchange capacity of clay plays an important role in many industrial applications ranging from radioactive waste disposal to cosmetics. However, swelling or shrinking of clay platelets due to water and ions adsorption in the interstitial zone is also a well-known phenomenon. For their applications, it is crucial to understand the stability of these layered materials, especially after exchange of interstitial ions with surrounding ions having different properties. Here, we probed the role of cation size on swelling pressure and free energy profile. We used molecular simulations to investigate the stability of mica pore, having K⁺, Rb⁺, and Cs⁺ ions. We performed a series of grand canonical Monte Carlo simulations at various pore widths. We probed water adsorption in mica pores from which disjoining pressure, grand potential (swelling free energy), and structural properties of confined water and ions were calculated. While the behavior of these three systems is similar qualitatively because of similar hydration properties of ions, significant differences are observed at the quantitative level due to changes in the hydration structure of cations. The global minimum in swelling free energy is found to be at the smaller pore widths (first minimum) for Rb- and K-mica and at bigger pore widths (second minimum) for Cs-mica pores. We find that ±0.1 Å change in the interstitial cation size leads to a -15 to 5% change in equilibrium loading of adsorbed water and -2 to 35% change in swelling. Our thermodynamic analysis reveals an intricate interplay between enthalpic and entropic contributions caused by the structural change of water in the pores due to the hydration of interstitial cations.

Sorption Mechanisms of Chemicals in Soils

Sorption of chemicals onto soil particle surfaces is an important process controlling their availability for uptake by organisms and loss from soils to ground and surface waters. The mechanisms of chemical sorption are inner- and outer-sphere adsorption and precipitation onto mineral surfaces. Factors that determine the sorption behavior are properties of soil mineral and organic matter surfaces and properties of the sorbing chemicals (including valence, electron configuration, and hydrophobicity). Because soils are complex heterogeneous mixtures, measuring sorption mechanisms is challenging; however, advancements analytical methods have made direct determination of sorption mechanisms possible. In this review, historical and modern research that supports the mechanistic understanding of sorption mechanisms in soils is discussed. Sorption mechanisms covered include cation exchange, outer-sphere adsorption, inner-sphere adsorption, surface precipitation, and ternary adsorption complexes.

Molecular dynamics simulations of the colloidal interaction between smectite clay nanoparticles in liquid water

Colloidal interactions between clay nanoparticles have been studied extensively because of their strong influence on the hydrology and mechanics of many soils and sedimentary media. The predominant theory used to describe these interactions is the Derjaguin-Landau-Verwey-Overbeek (DLVO) model, a framework widely applied in colloidal and interfacial science that accurately predicts the interactions between charged surfaces across water films at distances greater than ~ 3 nm (i.e., ten water monolayers). Unfortunately, the DLVO model is inaccurate at the shorter interparticle distances that predominate in most subsurface environments. For example, it inherently cannot predict the existence of equilibrium states wherein clay particles adopt interparticle distances equal to the thickness of one, two, or three water monolayers. Molecular dynamics (MD) simulations have the potential to provide detailed information on the free energy of interaction between clay nanoparticles; however, they have only been used to examine clay swelling and aggregation at interparticle distances below 1 nm. We present the first MD simulation predictions of the free energy of interaction of smectite clay nanoparticles in the entire range of interparticle distances from the large interparticle distances where the DLVO model is accurate (>3 nm) to the short-range swelling states where non-DLVO interactions predominate (<1 nm). Our simulations examine a range of salinities (0.0 to 1.0 M NaCl) and counterion types (Na, K, Ca) and establish a detailed picture of the breakdown of the DLVO model. In particular, they confirm previous theoretical suggestions of the existence of a strong non-DLVO attraction with a range of ~ 3 nm arising from specific ion-clay Coulomb interactions in the electrical double layer.

Novel light-emitting clays with structural Tb3+ and Eu3+ for chromate anion detection

Tb3+ and Eu3+ ions were encapsulated for the first time in the inorganic layers of a synthetic saponite clay following a one-pot synthetic approach. The co-presence of the two metal ions led to tuneable light-emitting properties, promoted by an efficient Tb3+ → Eu3+ energy transfer and enhanced Stokes shift character. To our knowledge, the so-prepared luminescent material was tested for the first time as an optical sensor for the detection of chromate anions in water.