Real-time observation of cation exchange kinetics and dynamics at the muscovite-water interface

Water and ions in electrified silica nano-pores: a molecular dynamics study

Solid-liquid interfaces (SLIs) are ubiquitous in science and technology from the development of energy storage devices to the chemical reactions occurring in the biological milieu. In systems involving aqueous saline solutions as the liquid, both the water and the ions are routinely exposed to an electric field, whether the field is externally applied, or originating from the natural surface charges of the solid. In the current study a molecular dynamics (MD) framework is developed to study the effect of an applied voltage on the behaviour of ionic solutions located in a ∼7 nm pore between two uncharged hydrophilic silica slabs. We systematically investigate the dielectric properties of the solution and the organisation of the water and ions as a function of salt concentration. In pure water, the interplay between interfacial hydrogen bonds and the applied field can induce a significant reorganisation of the water orientation and densification at the interface. In saline solutions, at low concentrations and voltages the interface dominates the whole system due to the extended Debye length resulting in a dielectric constant lower than that for the bulk solution. An increase in salt concentration or voltage brings about more localized interfacial effects resulting in dielectric properties closer to that of the bulk solution. This suggests the possibility of tailoring the system to achieve the desired dielectric properties. For example, at a specific salt concentration, interfacial effects can locally increase the dielectric constant, something that could be exploited for energy storage.

Insights into the structure and dynamics of K+ ions at the muscovite-water interface from machine learning potential simulations

The surfaces of many minerals are covered by naturally occurring cations that become partially hydrated and can be replaced by hydronium or other cations when the surface is exposed to water or an aqueous solution. These ion exchange processes are relevant to various chemical and transport phenomena, yet elucidating their microscopic details is challenging for both experiments and simulations. In this work, we make a first step in this direction by investigating the behavior of the native K+ ions at the interface between neat water and the muscovite mica (001) surface with ab-initio-based machine learning molecular dynamics and enhanced sampling simulations. Our results show that the desorption of the surface K+ ions in pure ion-free water has a significant free energy barrier irrespective of their local surface arrangement. In contrast, facile K+ diffusion between mica's ditrigonal cavities characterized by different Al/Si orderings is observed. This behavior suggests that the K+ ions may favor a dynamic disordered surface arrangement rather than complete desorption when exposed to deionized water.

Two-dimensional ice-like water adlayers on a mica surface with and without a graphene coating under ambient conditions

Water tends to wet all hydrophilic surfaces under ambient conditions, and the first water adlayers on solids are important for a broad range of physicochemical phenomena and technological processes, including corrosion, wetting, lubrication, anti-icing, catalysis, and electrochemistry. Unfortunately, challenges in characterizing the first water adlayer in the laboratory have hampered molecular-level understanding of the contact water structure. Herein, we present the first ab initio molecular dynamics simulation evidence of a previously unreported ice-like adlayer structure (named as Ice-AL-II) on a prototype mica surface under ambient conditions. Calculation showed that the newly identified Ice-AL-II structure is more stable than the widely recognized ice-adlayer structure on mica surfaces (named as Ice-AL-I). Ice-AL-II exhibited a face-centered corner-cut tetragon (or a face-centered irregular pentagon) pattern of a hydrogen-bonded network. The center of the corner-cut tetragon was occupied by either a K+ cation or a water molecule with two H atoms pinned by the mica (100) via double hydrogen bonds. Our simulation also suggested that bilayer Ice-AL-II favors AA stacking rather than AB stacking. Interestingly, when a graphene sheet was coated on top of the ice-like adlayer, the stability of Ice-AL-II was further enhanced. In contrast, due to its strongly puckered structure, the Ice-AL-I structure could be crushed into a near-Ice-AL-II structure by the graphene coating. Ice-AL-II is thus proposed as a promising candidate for the ice-like structure on a mica surface detected by scanning polarization force microscopy and by atomic force microscopy between a graphene coating and a mica surface.

Superconducting Two-Dimensional FeSe Grown on the Fe-Enriched Interface

Two-dimensional (2D) tetragonal FeSe has sparked extensive research interest owing to its tunable superconductivity, providing valuable insights into the design of high-temperature superconductors. Currently, the intricate Fe-Se phase diagram poses a challenge to the controlled synthesis of superconducting 2D FeSe in a pure tetragonal phase. Here, we exploit the ion-exchange property of fluorophlogopite mica to devise a straightforward approach for the phase-controlled synthesis of tetragonal FeSe on an Fe-enriched mica surface within a molten salt environment. This method successfully produces highly crystalline FeSe in a pure tetragonal phase with adjustable thickness. We investigated the surface composition of the postgrowth mica substrate using various microscopic and spectroscopic characterizations to highlight the importance of the Fe-enriched growth interface in the phase-selective synthesis of 2D tetragonal FeSe. The obtained 2D FeSe exhibited 2D superconductivity, comparable to that of FeSe mechanically exfoliated from bulk crystals, confirming the high quality of our samples. Beyond tetragonal FeSe, 2D antiferromagnetic FeTe and superconducting FeSxSeyTe1-x-y have been phase-selectively synthesized via this approach. Our study elucidates the significance of the growth interface on the phase-selective synthesis of 2D materials and presents potential opportunities for the phase-controlled synthesis of 2D multiphase materials via the rational design of the growth interface.

NaCl, MgCl2, and AlCl3 Surface Coverages on Fused Silica and Adsorption Free Energies at pH 4 from Nonlinear Optics

We employ amplitude- and phase-resolved second harmonic generation experiments to probe interactions of fused silica:aqueous interfaces with Al3+, Mg2+, and Na+ cations at pH 4 and as a function of metal cation concentration. We quantify the second-order nonlinear susceptibility and the total interfacial potential in the presence and absence of a 10 mM screening electrolyte to understand the influence of charge screening on cation adsorption. Strong cation:surface interactions are observed in the absence of the screening electrolyte. The total potential is then employed to estimate the total number of absorbed cations cm-2. The contributions to the total potential from the bound and mobile charges were separated using Gouy-Chapman-Stern model estimates. All three cations bind fully reversibly, indicating physisorption as the mode of interaction. Of the isotherm models tested, the Kd adsorption model fits the data with binding constants of 3-30 and ∼300 mol-1 for the low (<0.1 mM) and high (0.1-3 mM) concentration regimes, corresponding to adsorption free energies of -13 to -18 and -24 kJ mol-1 at room temperature, respectively. The maximum surface coverages are around 1013 cations cm-2, matching the number of deprotonated silanol groups on silica at pH 4. Clear signs of decoupled Stern and diffuse layer nonlinear optical responses are observed and found to be cation-specific.

Energy harvesting from water streaming at charged surface

Water flowing at a charged surface may produce electricity, known as streaming current/potentials, which may be traced back to the 19th century. However, due to the low gained power and efficiencies, the energy conversion from streaming current was far from usable. The emergence of micro/nanofluidic technology and nanomaterials significantly increases the power (density) and energy conversion efficiency. In this review, we conclude the fundamentals and recent progress in electrical double layers at the charged surface. We estimate the generated power by hydrodynamic energy dissipation in multi-scaling flows considering the viscous systems with slipping boundary and inertia systems. Then, we review the coupling of volume flow and current flow by the Onsager relation, as well as the figure of merits and efficiency. We summarize the state-of-the-art of electrokinetic energy conversions, including critical performance metrics such as efficiencies, power densities, and generated voltages in various systems. We discuss the advantages and possible constraints by the figure of merits, including single-phase flow and flying droplets.

Direct Experimental Observations of Ion Distributions during Overcharging at the Muscovite-Water Interface by Adsorption of Rb+ and Halides (Cl- , Br- , I- ) at High Salinity

Classical electric double layer (EDL) models have been widely used to describe ion distributions at charged solid-water interfaces in dilute electrolytes. However, the chemistry of EDLs remains poorly constrained at high ionic strength where ion-ion correlations control non-classical behavior such as overcharging, i. e., the accumulation of counter-ions in amounts exceeding the substrate's surface charge. Here, we provide direct experimental observations of correlated cation and anion distributions adsorbed at the muscovite (001)-aqueous electrolyte interface as a function of dissolved RbBr concentration ([RbBr]=0.01-5.8 M) using resonant anomalous X-ray reflectivity. Our results show alternating cation-anion layers in the EDL when [RbBr]≳100 mM, whose spatial extension (i. e., ~20 Å from the surface) far exceeds the dimension of the classical Stern layer. Comparison to RbCl and RbI electrolytes indicates that these behaviors are sensitive to the choice of co-ion. This new in-depth molecular-scale understanding of the EDL structure during transition from classical to non-classical regimes supports the development of realistic EDL models for technologies operating at high salinity such as water purification applications or modern electrochemical storage.

Nanoscale probing of local dielectric changes at the interface between solids and aqueous saline solutions

The mobility of dissolved ions and charged molecules at interfaces underpins countless processes in science and technology. Experimentally, this is typically measured from the averaged response of the charges to an electrical potential. High-resolution Atomic Force Microscopy (AFM) can image single adsorbed ions and molecules at solid-liquid interfaces, but probing the associated dynamics remains highly challenging. One possible strategy is to investigate the response of the species of interest to a highly localized AC electric field in an approach analogous to dielectric spectroscopy. The dielectric force experienced by the AFM tip apex is modulated by the dielectric properties of the sample probed, itself sensitive to the mobilities of solvated charges and dipoles. Previous work successfully used this approach to quantify the dielectric constant of thin samples, but with limited spatial resolution. Here we propose a strategy to simultaneously map the nanoscale topography and local dielectric variations across a range of interfaces by conducting high-resolution AFM imaging concomitantly with electrical AC measurements in a multifrequency approach. The strategy is tested over a 500 MHz bandwidth in pure liquids with different dielectric constants and in saline aqueous solutions. In liquids with higher dielectric constants, the system behaves as inductive-resistive-capacitive but the adjunction of ions removes the inductive resonances and precludes measurements at higher frequencies. High-resolution imaging is demonstrated over single graphene oxide (GrO) flakes with simultaneous but decoupled dielectric measurements. The dielectric constant is consistent and reproducible across liquids, except at higher salt concentrations where frequency-dependent effects occur. The results suggest the strategy is suitable for nanometre-level mapping of the dielectric properties of solid-liquid interfaces, but more work is needed to fully understand the different physical effects underpinning the measurements.

Quantification of Stern Layer Water Molecules, Total Potentials, and Energy Densities at Fused Silica:Water Interfaces for Adsorbed Alkali Chlorides, CTAB, PFOA, and PFAS

We have employed amplitude- and phase-resolved second-harmonic generation spectroscopy to investigate ion-specific effects of monovalent cations at the fused silica:water interface maintained under acidic, neutral, and alkaline conditions. We find a negligible dependence of the total potential (as negative as -400 mV at pH 14), the second-order nonlinear susceptibility (as large as 1.5 × 10-21 m2 V-1 at pH 14), the number of Stern layer water molecules (1 × 1015 cm-2 at pH 5.8), and the energy associated with water alignment upon going from neutral to high pH (ca. -24 kJ mol-1 to -48 kJ mol-1 at pH 13 and 14, close to the cohesive energy of liquid water but smaller than that of ice) on chlorides of the alkali series (M+ = Li+, Na+, K+, Rb+, and Cs+). Attempts are presented to provide estimates for the molecular hyperpolarizability of the cations and anions in the Stern layer at high pH, which arrive at ca. 20-fold larger values for αtotal ions(2) = αM+(2) + αOH-(2) + αCl-(2) when compared to water's molecular hyperpolarizability estimate from theory and point to a sizable contribution of deprotonated silanol groups at high pH. In contrast to the alkali series, a pronounced dependence of the total potential and the second-order nonlinear susceptibility on monovalent cationic (cetrimonium bromide, CTAB) and anionic (perfluorooctanoic and perfluorooctanesulfonic acid, PFOA and PFOS) surfactants was quantifiable. Our findings are consistent with a low surface coverage of the alkali cations and a high surface coverage of the surfactants. Moreover, they underscore the important contribution of Stern layer water molecules to the total potential and second-order nonlinear susceptibility. Finally, they demonstrate the applicability of heterodyne-detected second-harmonic generation spectroscopy for identifying perfluorinated acids at mineral:water interfaces.

Robotic pendant drop: containerless liquid for μs-resolved, AI-executable XPCS

The dynamics and structure of mixed phases in a complex fluid can significantly impact its material properties, such as viscoelasticity. Small-angle X-ray Photon Correlation Spectroscopy (SA-XPCS) can probe the spontaneous spatial fluctuations of the mixed phases under various in situ environments over wide spatiotemporal ranges (10-6-103 s /10-10-10-6 m). Tailored material design, however, requires searching through a massive number of sample compositions and experimental parameters, which is beyond the bandwidth of the current coherent X-ray beamline. Using 3.7-μs-resolved XPCS synchronized with the clock frequency at the Advanced Photon Source, we demonstrated the consistency between the Brownian dynamics of ~100 nm diameter colloidal silica nanoparticles measured from an enclosed pendant drop and a sealed capillary. The electronic pipette can also be mounted on a robotic arm to access different stock solutions and create complex fluids with highly-repeatable and precisely controlled composition profiles. This closed-loop, AI-executable protocol is applicable to light scattering techniques regardless of the light wavelength and optical coherence, and is a first step towards high-throughput, autonomous material discovery.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Resolving the intrinsic short-range ordering of K+ ions on cleaved muscovite mica

Muscovite mica, KAl2(Si3Al)O10(OH)2, is a common layered phyllosilicate with perfect cleavage planes. The atomically flat surfaces obtained through cleaving lend themselves to scanning probe techniques with atomic resolution and are ideal to model minerals and clays. Despite the importance of the cleaved mica surfaces, several questions remain unresolved. It is established that K+ ions decorate the cleaved surface, but their intrinsic ordering - unaffected by the interaction with the environment - is not known. This work presents clear images of the K+ distribution of cleaved mica obtained with low-temperature non-contact atomic force microscopy (AFM) under ultra-high vacuum (UHV) conditions. The data unveil the presence of short-range ordering, contrasting previous assumptions of random or fully ordered distributions. Density functional theory (DFT) calculations and Monte Carlo simulations show that the substitutional subsurface Al3+ ions have an important role for the surface K+ ion arrangement.

Selenite and Selenate Sequestration during Coprecipitation with Barite: Insights from Mineralization Processes of Adsorption, Nucleation, and Growth

Coprecipitation of selenium oxyanions with barite is a facile way to sequester Se in the environments. However, the chemical composition of Se-barite coprecipitates usually deviates from that predicted from thermodynamic calculations. This discrepancy was resolved by considering variations in nucleation and growth rates controlled by ion-mineral interactions, solubility, and interfacial energy. For homogeneous precipitation, ∼10% of sulfate, higher than thermodynamic predictions (<0.3%), was substituted by Se(IV) or Se(VI) oxyanion, which was attributed to adsorption-induced entrapment during crystal growth. For heterogeneous precipitation, thiol- and carboxylic-based organic films, utilized as model interfaces to mimic the natural organic-abundant environments, further enhanced the sequestration of Se(VI) oxyanions (up to 41-92%) with barite. Such enhancement was kinetically driven by increased nucleation rates of selenate-rich barite having a lower interfacial energy than pure barite. In contrast, only small amounts of Se(IV) oxyanions (∼1%) were detected in heterogeneous coprecipitates mainly due to a lower saturation index of BaSeO₃ and deprotonation degree of Se(IV) oxyanion at pH 5.6. These roles of nanoscale mineralization mechanisms observed during composition selection of Se-barite could mark important steps toward the remediation of contaminants through coprecipitation.

Fukushima and Chernobyl: Similarities and Differences of Radiocesium Behavior in the Soil-Water Environment

In the wake of Chernobyl and Fukushima accidents, radiocesium has become a radionuclide of most environmental concern. The ease with which this radionuclide moves through the environment and is taken up by plants and animals is governed by its chemical forms and site-specific environmental characteristics. Distinctions in climate and geomorphology, as well as ¹³⁷Cs speciation in the fallout, result in differences in the migration rates of ¹³⁷Cs in the environment and rates of its natural attenuation. In Fukushima areas, ¹³⁷Cs was strongly bound to soil and sediment particles, with its bioavailability being reduced as a result. Up to 80% of the deposited ¹³⁷Cs on the soil was reported to be incorporated in hot glassy particles (CsMPs) insoluble in water. Disintegration of these particles in the environment is much slower than that of Chernobyl-derived fuel particles. The higher annual precipitation and steep slopes in Fukushima-contaminated areas are conducive to higher erosion and higher total radiocesium wash-off. Among the common features in the ¹³⁷Cs behavior in Chernobyl and Fukushima are a slow decrease in the ¹³⁷Cs activity concentration in small, closed, and semi-closed lakes and its particular seasonal variations: increase in the summer and decrease in the winter.

Water Structure in the Electrical Double Layer and the Contributions to the Total Interfacial Potential at Different Surface Charge Densities

The electric double layer governs the processes of all charged surfaces in aqueous solutions; however, elucidating the structure of the water molecules is challenging for even the most advanced spectroscopic techniques. Here, we present the individual Stern layer and diffuse layer OH stretching spectra at the silica/water interface in the presence of NaCl over a wide pH range using a combination of vibrational sum frequency generation spectroscopy, heterodyned second harmonic generation, and streaming potential measurements. We find that the Stern layer water molecules and diffuse layer water molecules respond differently to pH changes: unlike the diffuse layer, whose water molecules remain net-oriented in one direction, water molecules in the Stern layer flip their net orientation as the solution pH is reduced from basic to acidic. We obtain an experimental estimate of the non-Gouy-Chapman (Stern) potential contribution to the total potential drop across the insulator/electrolyte interface and discuss it in the context of dipolar, quadrupolar, and higher order potential contributions that vary with the observed changes in the net orientation of water in the Stern layer. Our findings show that a purely Gouy-Chapman (Stern) view is insufficient to accurately describe the electrical double layer of aqueous interfaces.

Beyond Thermodynamics: Assessing the Dynamical Softness of Hydrated Ions from First Principles

Ion (de)hydration is a key rate-determining step in interfacial processes from corrosion to electrochemical energy storage. However, predicting the kinetics of ion (de)hydration remains challenging, prompting the use of static proxies such as hydration energy and valence. While useful for assessing thermodynamic preferences, such descriptors cannot fully capture the dynamical softness of the hydration shell that dictates kinetics. Accordingly, we use first-principles molecular dynamics to analyze hydration shell softness for a diverse set of metal cations. Three dynamic metrics are introduced to intuitively describe the bond rigidity, shape deformability, and exchange fluidity of the solvation shell. Together, these metrics capture the relevant physics in the static descriptors, while offering a far more complete and efficient representation for the overall propensity for (de)hydration. Application to the hydrated ion set demonstrates a weak connection between dynamical softness and hydration energy, confirming that dynamical descriptors of hydration are key for correctly describing ion transfer processes.

Ion exchange in atomically thin clays and micas

The physical properties of clays and micas can be controlled by exchanging ions in the crystal lattice. Atomically thin materials can have superior properties in a range of membrane applications, yet the ion-exchange process itself remains largely unexplored in few-layer crystals. Here we use atomic-resolution scanning transmission electron microscopy to study the dynamics of ion exchange and reveal individual ion binding sites in atomically thin and artificially restacked clays and micas. We find that the ion diffusion coefficient for the interlayer space of atomically thin samples is up to 10⁴ times larger than in bulk crystals and approaches its value in free water. Samples where no bulk exchange is expected display fast exchange at restacked interfaces, where the exchanged ions arrange in islands with dimensions controlled by the moiré superlattice dimensions. We attribute the fast ion diffusion to enhanced interlayer expandability resulting from weaker interlayer binding forces in both atomically thin and restacked materials. This work provides atomic scale insights into ion diffusion in highly confined spaces and suggests strategies to design exfoliated clay membranes with enhanced performance.

Real-time tracking of ionic nano-domains under shear flow

The behaviour of ions at solid-liquid interfaces underpins countless phenomena, from the conduction of nervous impulses to charge transfer in solar cells. In most cases, ions do not operate as isolated entities, but in conjunction with neighbouring ions and the surrounding solution. In aqueous solutions, recent studies suggest the existence of group dynamics through water-mediated clusters but results allowing direct tracking of ionic domains with atomic precision are scarce. Here, we use high-speed atomic force microscopy to track the evolution of Rb⁺, K⁺, Na⁺ and Ca²⁺ nano-domains containing 20 to 120 ions adsorbed at the surface of mica in aqueous solution. The interface is exposed to a shear flow able to influence the lateral motion of single ions and clusters. The results show that, when in groups, metal ions tend to move with a relatively slow dynamics, as can be expected from a correlated group motion, with an average residence timescale of ~ 1-2 s for individual ions at a given atomic site. The average group velocity of the clusters depends on the ions' charge density and can be explained by the ion's hydration state. The lateral shear flow of the fluid is insufficient to desorb ions, but indirectly influences the diffusion dynamics by acting on ions in close vicinity to the surface. The results provide insights into the dynamics of ion clusters when adsorbed onto an immersed solid under shear flow.

Ion correlations drive charge overscreening and heterogeneous nucleation at solid-aqueous electrolyte interfaces

Ion distributions at charged solid–water interfaces, referred to as the electrical double layer (EDL), are poorly understood at high ion concentrations, in part due to the lack of molecular-scale descriptions of the interactions between adsorbed hydrated ions. Here, direct visualization of the salinity-dependent evolution of EDL structure reveals molecular origins of nonclassical transformation of the EDL, in which charge overscreening and heterogeneous nucleation are driven by ion–ion correlations at the interfaces. This manifestation of the atomistic basis of nonclassical behaviors provides a much-needed understanding of the impact of ion cooperativity at charged interfaces for the development of predictive models for element transport in natural environments and advanced technologies for material growth and synthesis in saline environments.

Classical electrical double layer (EDL) models are foundational to the representation of atomistic structure and reactivity at charged interfaces. An important limitation to these models is their dependence on a mean-field approximation that is strictly valid for dilute aqueous solutions. Theoretical efforts to overcome this limitation are severely impeded by the lack of visualization of the structure over a wide range of ion concentration. Here, we report the salinity-dependent evolution of EDL structure at negatively charged mica–water interfaces, revealing transition from the Langmuir-type charge compensation in dilute salt solutions to nonclassical charge overscreening in highly concentrated solutions. The EDL structure in this overcharging regime is characterized by the development of both lateral positional correlation between adsorbed ions and vertical layering of alternating cations and anions reminiscent of the structures of strongly correlated ionic liquids. These EDL ions can spontaneously grow into nanocrystalline nuclei of ionic compounds at threshold ion concentrations that are significantly lower than the bulk solubility limit. These results shed light on the impact of ion cooperativity that drives heterogeneous nonclassical behaviors of the EDL in high-salinity conditions.

Stability of Two-Dimensional Ionic Clusters at Solid-Liquid Interfaces

The stability of two-dimensional clusters (2DCs) at the interface between ionic crystals and their solutions was investigated by molecular dynamics simulations. We found that 2DCs show a remarkable feature of odd-even alternation in stability. In NaCl and NaBr systems, the clusters containing an odd number of ions are more stable than those with an even number of ions, while in KCl systems, it is the other way round. Accordingly, the stability of water molecules in the first hydration shell of 2DCs also shows an odd-even alternation, which is consistent with the associated 2DCs. The odd-even alternation is discussed based on a competition mechanism between two factors: the Coulomb repulsion in charged 2DCs and the interaction between charges and water dipoles. Our discussion indicates that this odd-even alternation should be a universal feature in similar systems and would be important for understanding the nucleation and crystallization of solutions on ionic crystal surfaces.

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.

Effects of Hydration on the Adsorption of Benzohydroxamic Acid on the Lead-Ion-Activated Cassiterite Surface: A DFT Study

The strategy of enhancing the surface activity by preadsorption of metal ions (surface activation) is an effective way to promote the adsorption of surfactant on surfaces, which is very important in surface process engineering. However, the adsorption mechanism of surfactant (collector) on the surface preadsorbed by metal ions in the explicit solution phase is still poorly understood. Herein, the effects of hydration on the adsorption of benzohydroxamic acid (BHA) onto the oxide mineral surface before and after lead-ion activation are investigated by first-principles calculations, owing to its importance in the field of flotation. The results show that the direct adsorption of BHA on the hydrated surface is not thermodynamically allowed in the absence of metal ions. However, the adsorption of BHA onto the lead-ion-activated surface possesses a very low barrier and a very negative reaction energy difference, indicating that the adsorption of BHA on hydrated Pb²⁺ at cassiterite surface is very favorable in both thermodynamics and kinetics. In addition, the adsorption of BHA results in the dehydration of hydrated Pb²⁺. More interestingly, the surface hydroxyl groups could participate in and may promote the coordination adsorption through proton transfer. This work sheds some new lights on understanding the roles of interfacial water and the mechanisms of metal-ion surface activation.

Radiocesium distribution and mid-term dynamics in the ponds of the Fukushima Dai-ichi nuclear power plant exclusion zone in 2015-2019

This study analyzes the ¹³⁷Cs behavior in the ponds of Okuma Town from 2015 to 2019 in the Fukushima Dai-ichi nuclear power plant (FDNPP) exclusion zone. A decline in both particulate and dissolved ¹³⁷Cs activity concentrations was revealed. The decline rate constants for the particulate ¹³⁷Cs activity concentration were found to be higher than for the dissolved ¹³⁷Cs activity concentration. In terms of seasonality the dissolved ¹³⁷Cs concentrations were higher from June to October, depending on the specific pond and year, most likely due to temperature dependence of ¹³⁷Cs desorption from frayed edge sites of micaceous clay minerals. The apparent K_d(¹³⁷Cs) in the studied ponds, in absolute value, appeared to be much higher than that for closed and semi-closed lakes of the Chernobyl contaminated area; however, these were comparable to the values characteristic of the rivers and reservoirs of the FDNPP contaminated area. The apparent K_d(¹³⁷Cs) in the suspended sediment-water system was observed to decrease over time. It was hypothesized that this trend was associated with the decomposition of glassy hot particles. Relying on the theory of selective sorption and fixation, the exchangeable radiocesium interception potential, RIP^ex(K) was estimated using data on ¹³⁷Cs speciation in the surface bottom-sediment layer and its distribution in the sediment-water system. For the studied ponds, RIP^ex(K) was on the average 2050 mEq/kg, which is within the range of values measured in laboratory studies reported in the literature.

Nonclassical Behavior in Competitive Ion Adsorption at a Charged Solid-Water Interface

Ion adsorption at solid-water interfaces is commonly described by interactions between specific surface sites and adsorbed ions in classical models. However, energetic contributions from non-site-specific ion-ion interactions have been less well understood. Here, we report nonclassical behaviors observed during competitive adsorption between Sr²⁺ and Na⁺/Rb⁺ at the negatively charged muscovite mica (001)-water interface, revealing apparent controls of adsorbed ion speciation over the interfacial reactivity. In the absence of competing cations, Sr²⁺ adsorbs in approximately equivalent proportions of inner-sphere and outer-sphere complexes, whereas it adsorbs predominantly as an outer-sphere complex in the presence of Na⁺/Rb⁺. This transformation of adsorbed Sr²⁺ speciation significantly decreases its adsorption strength, as indicated by the ∼15-fold shift in the Sr²⁺ adsorption edge concentration, compared to that calculated from a classical Langmuir isotherm model developed on the basis of site-specific interactions. These observations highlight the importance of non-site-specific interactions in controlling the energetics of chemical reactions at the charged interface.

Radiocesium interaction with clay minerals: Theory and simulation advances Post-Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

Ion Distribution and Hydration Structure at Solid-Liquid Interface between NaCl Crystal and Its Solution

The interface structure between NaCl crystal and its solution has been investigated at the saturated concentration of 298 K by molecular dynamics simulations. We have found that there are many fine structures at this complex interface. Near the surface of crystal, most of Na⁺ only coordinate with water molecules, while almost all Cl^- coordinate with Na⁺ in addition to water molecules. An ion coordinating with more water molecules is farther away from the epitaxial position of lattice. As approaching to the interface, the first hydration shell of ions has the tendency of being ordered, while the orientation of dipole of water molecules in the first hydration shell becomes more disordered than that in the solution. Generally, the first hydration shell of Na⁺ is less affected by nearest Cl^-, whereas the first hydration shell of Cl^- is significantly affected by nearest Na⁺.

Directed nanoscale metal deposition by the local perturbation of charge screening at the solid-liquid interface

Understanding and directing electrochemical reactions below the micrometer scale is a long-standing challenge in electrochemistry. Confining reactions to nanoscale areas paradoxically requires both isolation from and communication with the bulk electrolyte in terms of electrochemical potential and access of ions, respectively. Here, we demonstrate the directed electrochemical deposition of copper nanostructures by using an oscillating nanoelectrode operated with an atomic force microscope (AFM). Strikingly, the writing is only possible in highly dilute electrolytes and for a particular combination of AFM and electrochemical parameters. We propose a mechanism based on cyclic charging and discharging of the electrical double layer (EDL). The extended screening length and slower charge dynamics in dilute electrolytes allow the nanoelectrode to operate inside, and disturb, the EDL even for large oscillation amplitudes (∼100 nm). Our unique approach can not only be used for controlled additive nano-fabrication but also provides insights into ion behavior and EDL dynamics at the solid-liquid interface.

Atomically thin micas as proton-conducting membranes

Monolayers of graphene and hexagonal boron nitride (hBN) are highly permeable to thermal protons^1,2. For thicker two-dimensional (2D) materials, proton conductivity diminishes exponentially, so that, for example, monolayer MoS₂ that is just three atoms thick is completely impermeable to protons¹. This seemed to suggest that only one-atom-thick crystals could be used as proton-conducting membranes. Here, we show that few-layer micas that are rather thick on the atomic scale become excellent proton conductors if native cations are ion-exchanged for protons. Their areal conductivity exceeds that of graphene and hBN by one to two orders of magnitude. Importantly, ion-exchanged 2D micas exhibit this high conductivity inside the infamous gap for proton-conducting materials³, which extends from ∼100 °C to 500 °C. Areal conductivity of proton-exchanged monolayer micas can reach above 100 S cm^-2 at 500 °C, well above the current requirements for the industry roadmap⁴. We attribute the fast proton permeation to ~5-Å-wide tubular channels that perforate micas' crystal structure, which, after ion exchange, contain only hydroxyl groups inside. Our work indicates that there could be other 2D crystals⁵ with similar nanometre-scale channels, which could help close the materials gap in proton-conducting applications.

Ion-Specific and pH-Dependent Hydration of Mica-Electrolyte Interfaces

Hydration forces play a crucial role in a wide range of phenomena in physics, chemistry, and biology. Here, we study the hydration of mica surfaces in contact with various alkali chloride solutions over a wide range of concentrations and pH values. Using atomic force microscopy and molecular dynamics simulations, we demonstrate that hydration forces consist of a superposition of a monotonically decaying and an oscillatory part, each with a unique dependence on the specific type of cation. The monotonic hydration force gradually decreases in strength with decreasing bulk hydration energy, leading to a transition from an overall repulsive (Li+, Na+) to an attractive (Rb+, Cs+) force. The oscillatory part, in contrast, displays a binary character, being hardly affected by the presence of strongly hydrated cations (Li+, Na+), but it becomes completely suppressed in the presence of weakly hydrated cations (Rb+, Cs+), in agreement with a less pronounced water structure in simulations. For both aspects, K+ plays an intermediate role, and decreasing pH follows the trend of increasing Rb+ and Cs+ concentrations.

Long-lived ionic nano-domains can modulate the stiffness of soft interfaces

Metal ions underpin countless processes at bio-interfaces, including maintaining electroneutrality, modifying mechanical properties and driving bioenergetic activity. These processes are typically described by ions behaving as independently diffusing point charges. Here we show that Na+ and K+ ions instead spontaneously form correlated nanoscale networks that evolve over seconds at the interface with an anionic bilayer in solution. Combining single-ion level atomic force microscopy and molecular dynamic simulations we investigate the configuration and dynamics of Na+, K+, and Rb+ at the lipid surface. We identify two distinct ionic states: the well-known direct electrostatic interaction with lipid headgroups and a water-mediated interaction that can drive the formation of remarkably long-lived ionic networks which evolve over many seconds. We show that this second state induces ionic network formation via correlative ion-ion interactions that generate an effective energy well of -0.4kBT/ion. These networks locally reduce the stiffness of the membrane, providing a spontaneous mechanism for tuning its mechanical properties with nanoscale precision. The ubiquity of water-mediated interactions suggest that our results have far-reaching implications for controlling the properties of soft interfaces.

Bottom-up synthesis of 2D Co-based metal–organic framework nanosheets by an ammonia-assisted strategy for tuning the crystal morphology

Two 2D Co₂(bim)₄ and Co(bim)(OAc) nanosheets were directly synthesized by an ammonia-modulated approach.

Atomic- and Molecular-Resolution Mapping of Solid-Liquid Interfaces by 3D Atomic Force Microscopy

Hydration layers are ubiquitous in life and technology. Hence, interfacial aqueous layers have a central role in a wide range of phenomena from materials science to molecular and cell biology. A complete understanding of those processes requires, among other things, the development of very-sensitive and high-resolution instruments. Three-dimensional atomic force microscopy (3D-AFM) represents the latest and most successful attempt to generate atomically resolved three-dimensional images of solid-liquid interfaces. This review provides an overview of the 3D-AFM operating principles and its underlying physics. We illustrate and explain the capability of the instrument to resolve atomic defects on crystalline surfaces immersed in liquid. We also illustrate some of its applications to imaging the hydration structures on DNA or proteins. In the last section, we discuss some perspectives on emerging applications in materials science and molecular biology.

Radiocesium interaction with clay minerals: Theory and simulation advances Post-Fukushima

Insights at the microscopic level of the process of radiocesium adsorption and interaction with clay mineral particles have improved substantially over the past several years, triggered by pressing social issues such as management of huge amounts of waste soil accumulated after the Fukushima Dai-ichi nuclear power plant accident. In particular, computer-based molecular modeling supported by advanced hardware and algorithms has proven to be a powerful approach. Its application can now generally encompass the full complexity of clay particle adsorption sites from basal surfaces to interlayers with inserted water molecules, to edges including fresh and weathered frayed ones. On the other hand, its methodological schemes are now varied from traditional force-field molecular dynamics on large-scale realizations composed of many thousands of atoms including water molecules to first-principles methods on smaller models in rather exacting fashion. In this article, we overview new understanding enabled by simulations across methodological variations, focusing on recent insights that connect with experimental observations, namely: 1) the energy scale for cesium adsorption on the basal surface, 2) progress in understanding the structure of clay edges, which is difficult to probe experimentally, 3) cesium adsorption properties at hydrated interlayer sites, 4) the importance of the size relationship between the ionic radius of cesium and the interlayer distance at frayed edge sites, 5) the migration of cesium into deep interlayer sites, and 6) the effects of nuclear decay of radiocesium. Key experimental observations that motivate these simulation advances are also summarized. Furthermore, some directions toward future solutions of waste soil management are discussed based on the obtained microscopic insights.

Lubricating properties of single metal ions at interfaces

The behaviour of ionic solutions confined in nanoscale gaps is central to countless processes, from biomolecular function to electrochemistry, energy storage and lubrication. However, no clear link exists between the molecular-level behaviour of the liquid and macroscopic observations. The problem mainly comes from the difficulty to interrogate a small number of liquid molecules. Here, we use atomic force microscopy to investigate the viscoelastic behaviour of pure water and ionic solutions down to the single ion level. The results show a glassy-like behaviour for pure water, with single metal ions acting as lubricants by reducing the elasticity of the nano-confined solution and the magnitude of the hydrodynamic friction. At small ionic concentration (<20 mM) the results can be quantitatively explained by the ions moving via a thermally-activated process resisted by the ion's hydration water (Prandtl-Tomlinson model). The model breaks down at higher salt concentrations due to ion-ion interaction effects that can no longer be neglected. The correlations are confirmed by direct sub-nanometre imaging of the interface at equilibrium. The results provide a molecular-level basis for explaining the tribological properties of aqueous solutions and suggest that ion-ion interactions create mesoscale effects that prevent a direct link between nanoscale and macroscopic measurements.

Electroviscous Dissipation in Aqueous Electrolyte Films with Overlapping Electric Double Layers

We use dynamic atomic force microscopy (AFM) to investigate the forces involved in squeezing out thin films of aqueous electrolyte between an AFM tip and silica substrates at variable pH and salt concentration. From amplitude and phase of the AFM signal we determine both conservative and dissipative components of the tip sample interaction forces. The measured dissipation is enhanced by up to a factor of 5 at tip-sample separations of ≈ one Debye length compared to the expectations based on classical hydrodynamic Reynolds damping with bulk viscosity. Calculating the surface charge density from the conservative forces using Derjaguin-Landau-Verwey-Overbeek (DLVO) theory in combination with a charge regulation boundary condition we find that the viscosity enhancement correlates with increasing surface charge density. We compare the observed viscosity enhancement with two competing continuum theory models: (i) electroviscous dissipation due to the electrophoretic flow driven by the streaming current that is generated upon squeezing out the counterions in the diffuse part of the electric double layer, and (ii) visco-electric enhancement of the local water viscosity caused by the strong electric fields within the electric double layer. While the visco-electric model correctly captures the qualitative trends observed in the experiments, a quantitative description of the data presumably requires more sophisticated simulations that include microscopic aspects of the distribution and mobility of ions in the Stern layer.

Effects of Ionic Strength on Arsenate Adsorption at Aluminum Hydroxide–Water Interfaces

Adsorption processes at mineral–water interfaces control the fate and transport of arsenic in soils and aquatic systems. Mechanistic and thermodynamic models to describe this phenomenon only consider inner-sphere complexes but recent observation of the simultaneous adsorption of inner- and outer-sphere arsenate on single crystal surfaces complicates this picture. In this study, we investigate the ionic strength-dependence of the macroscopic adsorption behavior and molecular-scale surface speciation of arsenate bound to gibbsite and bayerite. Arsenate adsorption decreases with increasing ionic strength on both minerals, with a larger effect at pH 4 than pH 7. The observed pH-dependence corresponds with a substantial decrease in surface charge at pH 7, as indicated by ζ-potential measurements. Extended X-ray absorption fine structure (EXAFS) spectroscopy finds that the number of second shell Al neighbors around arsenate is lower than that required for arsenate to occur solely as an inner-sphere surface complex. Together, these observations demonstrate that arsenate displays macroscopic and molecular-scale behavior consistent with the co-occurrence of inner- and outer-sphere surface complexes. This demonstrated that outer-sphere species can be responsible for strong adsorption of ions and suggests that environments experiencing an increase in salt content may induce arsenic release to water, especially under weakly acidic conditions.

Identification of Single Adsorbed Cations on Mica-Liquid Interfaces by 3D Force Microscopy

Force microscope provides atomically resolved images of surfaces immersed in a liquid. The presence of different chemical species in the interface (cations, anions, water, neutral atoms) complicates the adscription of the observed features to a given species. We develop a 3D atomic force microscopy method to identify the cations adsorbed on a mica surface from a potassium chloride solution. The method is based on measuring the peak value of the attractive force within the Stern layer. The maximum of the attractive force shows site-specific variations. The positions with the highest attractive force values are associated with the presence of adsorbed potassium ions, while the other positions are associated with a local depletion of the hydration layer. This criterion provides a surface coverage of K cations that is consistent with the one reported by other techniques.