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

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely, dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data are relevant to another type of non-DLVO interaction, namely, the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. Our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

Atomic-scale structure of interfacial water on gel and liquid phase lipid membranes

Hydration of biological membranes is essential to a wide range of biological processes. In particular, it is intrinsically linked to lipid thermodynamic properties, which in turn influence key cell functions such as ion permeation and protein mobility. Experimental and theoretical studies of the surface of biomembranes have revealed the presence of an interfacial repulsive force, which has been linked to hydration or steric effects. Here, we directly characterise the atomic-scale structure of water near supported lipid membranes of 1,2-dimyristoyl-sn-glycero-3-phosphocholine in their gel and liquid phase through three-dimensional atomic force microscopy (3D AFM). First, we demonstrate the ability to probe the morphology of interfacial water of lipid bilayers in both phases with sub-molecular resolution by using ultrasharp tips. We then visualise the molecular arrangement of water at the lipid surface at different temperatures. Our experiments reveal that water is organised in multiple hydration layers on both the solid-ordered and liquid-disordered lipid phases. Furthermore, we observe a monotonic repulsive force, which becomes relevant only in the liquid phase. These results offer new insights into the water structuring near soft biological surfaces, and demonstrate the importance of investigating it with vertical and lateral sub-molecular resolution.

Hygroscopic Growth of Adsorbed Water Films on Smectite Clay Particles

Hygroscopic growth of adsorbed water films on clay particles underlies a number of environmental science questions, from the air quality and climate impacts of mineral dust aerosols to the hydrology and mechanics of unsaturated soils and sedimentary rocks. Here, we use molecular dynamics (MD) simulations to establish the relation between adsorbed water film thickness (h) and relative humidity (RH) or disjoining pressure (Π), which has long been uncertain due to factors including sensitivity to particle shape, surface roughness, and aqueous chemistry. We present a new MD simulation approach that enables precise quantification of Π in films up to six water monolayers thick. We find that the hygroscopicity of phyllosilicate mineral surfaces increases in the order mica < K-smectite < Na-smectite. The relationship between Π and h on clay surfaces follows a double exponential decay with e-folding lengths of 2.3 and 7.5 Å. The two decay length scales are attributed to hydration repulsion and osmotic phenomena in the electrical double layer (EDL) at the clay-water interface.

Adsorbing DNA to Mica by Cations: Influence of Valency and Ion Type

Ion-mediated attraction between DNA and mica plays a crucial role in biotechnological applications and molecular imaging. Here, we combine molecular dynamics simulations and single-molecule atomic force microscopy experiments to characterize the detachment forces of single-stranded DNA at mica surfaces mediated by the metal cations Li+, Na+, K+, Cs+, Mg2+, and Ca2+. Ion-specific adsorption at the mica/water interface compensates (Li+ and Na+) or overcompensates (K+, Cs+, Mg2+, and Ca2+) the bare negative surface charge of mica. In addition, direct and water-mediated contacts are formed between the ions, the phosphate oxygens of DNA, and mica. The different contact types give rise to low- and high-force pathways and a broad distribution of detachment forces. Weakly hydrated ions, such as Cs+ and water-mediated contacts, lead to low detachment forces and high mobility of the DNA on the surface. Direct ion-DNA or ion-surface contacts lead to significantly higher forces. The comprehensive view gained from our combined approach allows us to highlight the most promising cations for imaging in physiological conditions: K+, which overcompensates the negative mica charge and induces long-ranged attractions. Mg2+ and Ca2+, which form a few specific and long-lived contacts to bind DNA with high affinity.

Ion adsorption and hydration forces: a comparison of crystalline mica vs. amorphous silica surfaces

Hydration forces are ubiquitous in nature and technology. Yet, the characterization of interfacial hydration structures and their dependence on the nature of the substrate and the presence of ions have remained challenging and controversial. We present a systematic study using dynamic Atomic Force Microscopy of hydration forces on mica surfaces and amorphous silica surfaces in aqueous electrolytes containing chloride salts of various alkali and earth alkaline cations of variable concentrations at pH values between 3 and 9. Our measurements with ultra-sharp AFM tips demonstrate the presence of both oscillatory and monotonically decaying hydration forces of very similar strength on both atomically smooth mica and amorphous silica surfaces with a roughness comparable to the size of a water molecule. The characteristic range of the forces is approximately 1 nm, independent of the fluid composition. Force oscillations are consistent with the size of water molecules for all conditions investigated. Weakly hydrated Cs+ ions are the only exception: they disrupt the oscillatory hydration structure and induce attractive monotonic hydration forces. On silica, force oscillations are also smeared out if the size of the AFM tip exceeds the characteristic lateral scale of the surface roughness. The observation of attractive monotonic hydration forces for asymmetric systems suggests opportunities to probe water polarization.

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale. Our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

Adsorption effects and mechanisms of phosphorus by nanosized laponite

Phosphorus (P), an important macroelement for crops, may be lost into water systems by human activities and subsequently cause serious environmental problems such as eutrophication. Thus, the recovery of P from wastewater is essential. P can be adsorbed and recovered from wastewater using many natural, environmentally friendly clay minerals, however the adsorption ability is limited. Here we applied a synthesis nanosized clay mineral, laponite, to evaluate the P adsorption ability and molecular mechanisms of the adsorption process. We apply X-ray Photoelectron Spectroscopy (XPS) to observe the adsorption of inorganic phosphate onto laponite, and then measure the adsorption content of phosphate by laponite via batch experiments in different solution conditions, including pH, ionic species and concentrations. Then the molecular mechanisms of adsorption are analyzed by Transmission Electron Microscopy (TEM) and molecular modeling using Density Functional Theory (DFT). The results show that phosphate adsorbs to the surface and interlayer of laponite via hydrogen bonding, and the adsorption energies of the interlayer are greater than those of the surface. These bulk solution and molecular-scale results in a model system may provide new insights into the recovery of phosphorus by nanosized clay, with possible environmental engineering applications for P-pollution control and sustainable utilization of P sources.

Interaction of Nanoparticles in Electrolyte Solutions

The interaction between nanoparticles includes several components; however, the most frequently used are electrostatic, caused by overlapping double electrical layers, and London-van der Waals interactions, caused by quantum and thermodynamic fluctuations of electromagnetic fields. Only these two kinds of interaction are considered below. The electrostatic interaction is calculated based on the linearized Poisson-Boltzmann equation for particles with constant electrical potential of the surfaces (constant ζ potentials). An exact solution of the problem is obtained for both identical particles and particles of different sizes. For the London-van der Waals interaction, the screening of static fluctuations and the retardation of electromagnetic fields for the dispersive part of the interaction are taken into account. The total interaction energy for two particles was calculated for a range of possible nanoparticle sizes from 1 to 10³ nm and electrolyte concentration from 10^-2 to 10^-6 mol/L. The predominance of the London-van der Waals force over the shielded electrostatic repulsion force was found at high electrolyte concentrations in the range from 10^-2 to 10^-3 mol/L at large interparticle distances.

Long-Range Electrification of an Air/Electrolyte Interface and Probing Potential of Zero Charge by Conductive Amplitude-Modulated Atomic Force Microscopy

The structure of an electrical double layer (EDL) at the interface of electrode/electrolyte or air/electrode/electrolyte is a fundamental aspect, however not fully understood. The potential of zero charge (PZC) is one of the clues to dictate the EDL, where the excess charge on the electrode surface is zero. Here, a nanoscale configuration of immersion method was proposed by integrating an electrochemical system into conductive atomic force spectroscopy under the amplitude modulation (AM) mode and agarose gel as the solid-liquid electrolyte. The PZC of boron-doped diamond was determined to be at 0.2 V (vs Ag/AgCl). By AM spectroscopy, the capacitive force shows remote electrification without direct electrode/electrolyte contact, which is dependent on the population of ions at the air/electrolyte interface. The surface potential by alignment of water is also evaluated. Prospectively, our study could benefit applications such as PZC measurement and non-electrode electrochemical processes at the air/electrolyte interface.

RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations

The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 μs of simulation time. For the mono- and divalent cations Li⁺, Na⁺, K⁺, Cs⁺, Ca²⁺, Sr²⁺, and Ba²⁺, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.

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.

Water Adsorption on Mica Surfaces with Hydrophilicity Tuned by Counterion Types (Na, K, and Cs) and Structural Fluorination

The stability of adsorbed water films on mineral surfaces has far-reaching implications in the Earth, environmental, and materials sciences. Here, we use the basal plane of phlogopite mica, an atomically smooth surface of a natural mineral, to investigate water film structure and stability as a function of two features that modulate surface hydrophilicity: the type of adsorbed counterions (Na, K, and Cs) and the substitution of structural OH groups by F atoms. We use molecular dynamics simulations combined with in situ high-resolution X-ray reflectivity to examine surface hydration over a range of water loadings, from the adsorption of isolated water molecules to the formation of clusters and films. We identify four regimes characterized by distinct adsorption energetics and different sensitivities to cation type and mineral fluorination: from 0 to 0.5 monolayer film thickness, the hydration of adsorbed ions; from 0.5 to 1 monolayer, the hydration of uncharged regions of the siloxane surface; from 1 to 1.5 monolayer, the attachment of isolated water molecules on the surface of the first monolayer; and for >1.5 monolayer, the formation of an incipient electrical double layer at the mineral-water interface.

Hydration Structure at the Calcite-Water (10.4) Interface in the Presence of Rubidium Chloride

Solid-liquid interfaces are of significant importance in a multitude of geochemical and technological fields. More specifically, the solvation structure plays a decisive role in the properties of the interfaces. Atomic force microscopy (AFM) has been used to resolve the interfacial hydration structure in the presence and absence of ions. Despite many studies investigating the calcite-water interface, the impact of ions on the hydration structure at this interface has rarely been studied. Here, we investigate the calcite-water interface at various concentrations (ranging from 0 to 5 M) of rubidium chloride (RbCl) using three-dimensional atomic force microscopy (3D AFM). We present molecularly resolved images of the hydration structure at the interface. Interestingly, the characteristic pattern of the hydration structure appears similar regardless of the RbCl concentration. The presence of the ions is detected in an indirect manner by more frequent contrast changes and slice displacements.

Absence of anomalous underscreening in highly concentrated aqueous electrolytes confined between smooth silica surfaces

Recent surface forces apparatus experiments that measured the forces between two mica surfaces and a series of subsequent theoretical studies suggest the occurrence of universal underscreening in highly concentrated electrolyte solutions. We performed a set of systematic Atomic Force Spectroscopy measurements for aqueous salt solutions in a concentration range from 1 mM to 5 M using chloride salts of various alkali metals as well as mixed concentrated salt solutions (involving both mono- and divalent cations and anions), that mimic concentrated brines typically encountered in geological formations. Experiments were carried out using flat substrates and submicrometer-sized colloidal probes made of smooth oxidized silicon immersed in salt solutions at pH values of 6 and 9 and temperatures of 25 °C and 45 °C. While strong repulsive forces were observed for the smallest tip-sample separations, none of the conditions explored displayed any indication of anomalous long range electrostatic forces as reported for mica surfaces. Instead, forces are universally dominated by attractive van der Waals interactions at tip-sample separations of ≈2 nm and beyond for salt concentrations of 1 M and higher. Complementary calculations based on classical density functional theory for the primitive model support these experimental observations and display a consistent decrease in screening length with increasing ion concentration.

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.

Long-Range Ionic and Short-Range Hydration Effects Govern Strongly Anisotropic Clay Nanoparticle Interactions

The aggregation of clay particles in aqueous solution is a ubiquitous everyday process of broad environmental and technological importance. However, it is poorly understood at the all-important atomistic level since it depends on a complex and dynamic interplay of solvent-mediated electrostatic, hydrogen bonding, and dispersion interactions. With this in mind, we have performed an extensive set of classical molecular dynamics simulations (included enhanced sampling simulations) on the interactions between model kaolinite nanoparticles in pure and salty water. Our simulations reveal highly anisotropic behavior, in which the interaction between the nanoparticles varies from attractive to repulsive depending on the relative orientation of the nanoparticles. Detailed analysis reveals that at large separation (>1.5 nm), this interaction is dominated by electrostatic effects, whereas at smaller separations, the nature of the water hydration structure becomes critical. This study highlights an incredible richness in how clay nanoparticles interact, which should be accounted for in, for example, coarse-grained models of clay nanoparticle aggregation.

Correlation between Electrostatic and Hydration Forces on Silica and Gibbsite Surfaces: An Atomic Force Microscopy Study

The balance between hydration and Derjaguin-Landau-Verwey-Overbeek (DLVO) forces at solid-liquid interfaces controls many processes, such as colloidal stability, wetting, electrochemistry, biomolecular self-assembly, and ion adsorption. Yet, the origin of molecular scale hydration forces and their relation to the surface charge density that controls the continuum scale electrostatic forces is poorly understood. We argue that these two types of forces are largely independent of each other. To support this hypothesis, we performed atomic force microscopy experiments using intermediate-sized tips that enable the simultaneous detection of DLVO and molecular scale oscillatory hydration forces at the interface between composite gibbsite:silica-aqueous electrolyte interfaces. We extract surface charge densities from forces measured at tip-sample separations of 1.5 nm and beyond using DLVO theory in combination with charge regulation boundary conditions for various pH values and salt concentrations. We simultaneously observe both colloidal scale DLVO forces and oscillatory hydration forces for an individual crystalline gibbsite particle and the underlying amorphous silica substrate for all fluid compositions investigated. While the diffuse layer charge varies with pH as expected, the oscillatory hydration forces are found to be largely independent of pH and salt concentration, supporting our hypothesis that both forces indeed have a very different origin. Oscillatory hydration forces are found to be distinctly more pronounced on gibbsite than on silica. We rationalize this observation based on the distribution of hydroxyl groups available for H bonding on the two distinct surfaces.

Tip Charge Dependence of Three-Dimensional AFM Mapping of Concentrated Ionic Solutions

A molecular scale understanding of the organization and structure of a liquid near a solid surface is currently a major challenge in surface science. It has implications across different fields from electrochemistry and energy storage to molecular biology. Three-dimensional AFM generates atomically resolved maps of solid-liquid interfaces. The imaging mechanism behind those maps is under debate, in particular, for concentrated ionic solutions. Theory predicts that the observed contrast should depend on the tip's charged state. Here, by using neutrally, negatively, and positively charged tips, we demonstrate that the 3D maps depend on the tip's polarization. A neutral tip will explore the total particle density distribution (water and ions) while a charged tip will reveal the charge density distribution. The experimental data reproduce the key findings of the theory.

Anomalous interfacial dynamics of single proton charges in binary aqueous solutions

Our understanding of the dynamics of charge transfer between solid surfaces and liquid electrolytes has been hampered by the difficulties in obtaining interface, charge, and solvent-specific information at both high spatial and temporal resolution. Here, we measure at the single charge scale the dynamics of protons at the interface between an hBN crystal and binary mixtures of water and organic amphiphilic solvents (alcohols and acetone), evidencing a marked influence of solvation on interfacial dynamics. Applying single-molecule localization microscopy to emissive crystal defects, we observe correlated activation between adjacent ionizable surface defects, mediated by the transport of single excess protons along the solid/liquid interface. Solvent content has a nontrivial effect on interfacial dynamics, leading at intermediate water fraction to an increased surface diffusivity, as well as an increased affinity of the proton charges to the solid surface. Our measurements evidence the notable role of solvation on interfacial proton charge transport.

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations

Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca²⁺ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg²⁺ improves the structure of two specifically bound Mg²⁺ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.

Influence of surfactants on a pre-adsorbed cationic layer: Removal and modification

Removal of organic species from solid surfaces is a crucial process. The use of oppositely charged surfactants provides a potential method for enhanced removal. Neutron reflectometry has been used to investigate the complex behaviour of a pre-adsorbed and tenacious layer of the cationic surfactant didodecyldimethylammonium bromide (DDAB) on a mica surface, during exposure to different organic species in solution. The anionic surfactant sodium dodecylsulfate (SDS) was shown to be able to remove the cationic layer, but only if anionic micelles were present in solution. To facilitate comparison with the behaviour of a non-ionic surfactant, the direct adsorption of pentaethylene glycol monododecyl ether (C₁₂E₅) to mica was also studied; low surface coverage adsorption was seen at the critical micelle concentration and above. C₁₂E₅ was then found not to remove the cationic layer, but did include into the layer to some degree. The presence of cationic surfactant on the mica was however shown to significantly modify the adsorption behaviour of the non-ionic surfactant.

Atomically resolved interfacial water structures on crystalline hydrophilic and hydrophobic surfaces

Hydration layers are formed on hydrophilic crystalline surfaces immersed in water. Their existence has also been predicted for hydrophobic surfaces, yet the experimental evidence is controversial. Using 3D-AFM imaging, we probed the interfacial water structure of hydrophobic and hydrophilic surfaces with atomic-scale spatial resolution. We demonstrate that the atomic-scale structure of interfacial water on crystalline surfaces presents two antagonistic arrangements. On mica, a common hydrophilic crystalline surface, the interface is characterized by the formation of 2 to 3 hydration layers separated by approximately 0.3 nm. On hydrophobic surfaces such as graphite or hexagonal boron nitride (h-BN), the interface is characterized by the formation of 2 to 4 layers separated by about 0.5 nm. The latter interlayer distance indicates that water molecules are expelled from the vicinity of the surface and replaced by hydrocarbon molecules. This creates a new 1.5-2 nm thick interface between the hydrophobic surface and the bulk water. Molecular dynamics simulations reproduced the experimental data and confirmed the above interfacial water structures.

Long-range hydrodynamic forces in liquid FM-AFM

We study the effects of hydrodynamic forces in frequency-modulation AFM experiments (FM-AFM) in liquid. We first establish the theoretical equations needed to derive the interaction stiffness k _int and the damping β _int due to the hydrodynamic forces from the frequency shift and the excitation amplitude. We develop specific FM-AFM experiments to measure the variation of k _int and β _int over a large range of distance in water up to 200 µm. Comparison between theory and experiments point out that the evolution of k _int at short and long distance arises from unsteady hydrodynamic forces on the cantilever. On the other hand, β _int is small at long distance and diverges at short probe-surface distance, as predicted by the classical Reynolds sphere model.

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.

Modeling hydration-mediated ion-ion interactions in electrolytes through oscillating Yukawa potentials

Classical Poisson-Boltzmann theory represents a mean-field description of the electric double layer in the presence of only Coulomb interactions. However, aqueous solvents hydrate ions, which gives rise to additional hydration-mediated ion-ion interactions. Experimental and computational studies suggest damped oscillations to be a characteristic feature of these hydration-mediated interactions. We have therefore incorporated oscillating Yukawa potentials into the mean-field description of the electric double layer. This is accomplished by allowing the decay length of the Yukawa potential to be complex valued. Ion specificity emerges from assigning individual strengths and phases to the Yukawa potential for anion-anion, anion-cation, and cation-cation pairs as well as for anions and cations interacting with an electrode or macroion. Excluded volume interactions between ions are approximated by replacing the ideal gas entropy by that of a lattice gas. We derive mean-field equations for the Coulomb and Yukawa potentials and use their solutions to compute the differential capacitance for an isolated planar electrode and the pressure that acts between two planar, like-charged macroion surfaces. Attractive interactions appear if the surface charge density of the macroions is sufficiently small.

Monovalent - divalent cation competition at the muscovite mica surface: Experiment and theory

HYPOTHESIS

Ion adsorption on mineral surfaces depends on several factors, such as the mineral surface structure and the valency, size and hydration of the ion. In order to understand competitive adsorption at mineral surfaces, experimental techniques are required that can probe multiple ionic species at the same time. By comparing adsorption of two different cations, it should be possible to derive the factors governing ion adsorption. Divalent cations are expected to bind stronger to the negatively-charged muscovite surface than monovalent cations.

EXPERIMENTS

Here, the competition between the monovalent Cs⁺ and the divalent Ca²⁺ cation for adsorption at the muscovite mica basal plane was investigated using surface X-ray diffraction. Using an extended surface complexation model, we simultaneously fit the measured cation coverages and net surface charges reported in literature.

FINDINGS

In order to reproduce those complementary data sets, both cation adsorption and anion coadsorption were included in the surface complexation model. Moreover, the intrinsic muscovite surface charge and the maximum of available adsorption sites had to be reduced compared to existing literature values. Competition experiments revealed that the affinity of Cs⁺ for the muscovite surface is larger than the affinity of Ca²⁺, showing that hydration forces are more important than electrostatics.

Forces between solid surfaces in aqueous electrolyte solutions

This review addresses experimental findings obtained with direct force measurements between two similar or dissimilar solid surfaces in aqueous electrolyte solutions. Interpretation of these measurements is mainly put forward in terms of the classical theory of Derjaguin, Landau, Verwey, and Overbeek (DLVO). This theory invokes a superposition of attractive van der Waals forces and repulsive double layer forces. DLVO theory is shown to be extremely reliable, even in the case of multivalent ions. However, such a description is only successful, when appropriate surface charge densities, charge regulation characteristics, and ion pairing or complexation equilibria in solution are considered. Deviations from DLVO theory only manifest themselves at distances of typically below few nm. More long-ranged non-DLVO forces can be observed in some situations, particularly, in concentrated electrolyte solutions, in the presence of strongly adsorbed layers, or for hydrophobic surfaces. The latter forces probably originate from patch-charge surface heterogeneities, which can be induced by ion-ion correlation effects, charge fluctuations, or other types of surface heterogeneities.

Molecular-Scale Investigations Reveal Noncovalent Bonding Underlying the Adsorption of Environmental DNA on Mica

Mineral-soil organic matter (SOM including DNA, proteins, and polysaccharides) associations formed through various interactions, play a key role in regulating long-term SOM preservation. The mechanisms underlying DNA-mineral and DNA-protein/polysaccharide interactions at nanometer and molecular scales in environmentally relevant solutions remain uncertain. Here, we present a model mineral-SOM system consisting of mineral (mica)-nucleic acid (environmental DNA, eDNA)/protein (bovine serum albumin)/polysaccharide (alginate), and combine atomic force microscopy (AFM)-based dynamic force spectroscopy and PeakForce quantitative nanomechanical mapping using DNA-decorated tips. Single-molecule binding and adhesion force of eDNA to mineral and to mineral adsorbed by protein/polysaccharide reveal the noncovalent bonds and that systematically changing ion compositions, ionic strength, and pH result in significant differences in organic-organic and organic-mineral binding energies. Consistent with the bond-strength measurements, protein, rather than polysaccharide, promotes mineral-bound DNA molecules by ex situ AFM deposition observations in relatively high concentrations of divalent cation-containing acidic solutions. These molecular-scale determinations and nanoscale observations should substantially improve our understanding of how environmental factors influence the organic-mineral interfacial interactions through the synergy of collective noncovalent and/or covalent bonds in mineral-organic associations.