Profile of the static permittivity tensor of water at interfaces: consequences for capacitance, hydration interaction and ion adsorption

How Molecular Orientation Affects the Static Permittivity Profile of the Polar and Nonpolar Liquid-Liquid Interface

The dielectric permittivity across the liquid-liquid interface presents an intrinsic response, with respect to the instantaneous interface reference. We hypothesize that dielectric responses across the nonpolar and polar liquid-liquid interfaces have different behaviors and underlying mechanisms. Molecular dynamics simulations were used to compare and contrast the dielectric response of a nonpolar (1,2-dichloroethane/water) and a polar (1-octanol/water) liquid-liquid interface system. We found that the enhanced dielectric permittivity at the nonpolar interface is attributed to the increased water dipole orientation and polarization density. In the case of the polar interface, strong association of the immiscible solvents inhibits the molecular dipole orientation, counteracting the effect from the enhanced surface water polarization density and resulting in a standard dielectric response. Detailed knowledge of the hydrogen bond networks and molecular dipole orientation with respect to the specific instantaneous interfacial and bulk regions reveals the effect of molecular proximity and the interaction with the opposing interfacial molecules on the mechanism of the dielectric permittivity response across the liquid-liquid interface phase boundary.

Water, Not Salt, Causes Most of the Seebeck Effect of Nonisothermal Aqueous Electrolytes

A temperature difference between two electrolyte-immersed electrodes often yields a voltage Δψ between them. This electrolyte Seebeck effect is usually explained by cations and anions flowing differently in thermal gradients. However, using molecular simulations, we found almost the same Δψ for cells filled with pure water as with aqueous alkali halides. Water layering and orientation near polarizable electrodes cause a large temperature-dependent potential drop χ there. The difference in χ of hot and cold electrodes captures most of the thermovoltage, Δψ≈χ_{hot}-χ_{cold}.

Wetting Behavior of Kerogen Surfaces: Insights from Molecular Dynamics

In this study, the wettability of a kerogen surface, a key component of shale reservoirs, is investigated by using molecular dynamics simulations. Specifically, we examined the impact of droplet size and morphology as well as surface roughness on the water contact angles. The findings highlighted that the contact angle dependency on the droplet size intensifies with increased rigidity of the surface. Conversely, as the surface becomes more flexible and rougher, it gains hydrophilicity. The higher hydrophilicity stems from the ability of water molecules to penetrate the kerogen corrugations and form more hydrogen bonds with heteroatoms, particularly oxygen. Notably, the contact angle of kerogen hovers between 65 and 75°, thereby crossing the transition from an underoil hydrophilic to an underoil hydrophobic state. Consequently, minor alterations in the kerogen nanostructure can dramatically alter the wetting preference between water and oil. This insight is of paramount significance for refining strategies in managing fluid interactions in shale reservoirs such as geological carbon storage or oil extraction.

Collective modes and quantum effects in two-dimensional nanofluidic channels

Nanoscale fluid transport is typically pictured in terms of atomic-scale dynamics, as is natural in the real-space framework of molecular simulations. An alternative Fourier-space picture, that involves the collective charge fluctuation modes of both the liquid and the confining wall, has recently been successful at predicting new nanofluidic phenomena such as quantum friction and near-field heat transfer, that rely on the coupling of those fluctuations. Here, we study the charge fluctuation modes of a two-dimensional (planar) nanofluidic channel. Introducing confined response functions that generalize the notion of surface response function, we show that the channel walls exhibit coupled plasmon modes as soon as the confinement is comparable to the plasmon wavelength. Conversely, the water fluctuations remain remarkably bulk-like, with significant confinement effects arising only when the wall spacing is reduced to 7 Å. We apply the confined response formalism to predict the dependence of the solid-water quantum friction and thermal boundary conductance on channel width for model channel wall materials. Our results provide a general framework for Coulomb interactions of fluctuating matter under nanoscale confinement.

Multiscale Modeling of Aqueous Electric Double Layers

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Origin of dielectric polarization suppression in confined water from first principles

It has long been known that the dielectric constant of confined water should be different from that in bulk. Recent experiments have shown that it is vanishingly small, however the origin of the phenomenon remains unclear. Here we used ab initio molecular dynamics simulations (AIMD) and AIMD-trained machine-learning potentials to understand water's structure and electronic properties underpinning this effect. For the graphene and hexagonal boron-nitride substrates considered, we find that it originates in the spontaneous anti-parallel alignment of the water dipoles in the first two water layers near the solid interface. The interfacial layers exhibit net ferroelectric ordering, resulting in an overall anti-ferroelectric arrangement of confined water. Together with constrained hydrogen-bonding orientations, this leads to much reduced out-of-plane polarization. Furthermore, we directly contrast AIMD and simple classical force-field simulations, revealing important differences. This work offers insight into a property of water that is critical in modulating surface forces, the electric-double-layer formation and molecular solvation, and shows a way to compute it.

Transparent Colloids of Detonation Nanodiamond: Physical, Chemical and Biological Properties

Aqueous suspensions (colloids) containing detonation nano-diamond (DND) feature in most applications of DND and are an indispensable stage of its production; therefore, the interaction of DND with water is actively studied. However, insufficient attention has been paid to the unique physico-chemical and biological properties of transparent colloids with low DND content (≤0.1%), which are the subject of this review. Thus, such colloids possess giant dielectric permittivity which shows peculiar temperature dependence, as well as quasi-periodic fluctuations during slow evaporation or dilution. In these colloids, DND interacts with water and air to form cottonwool-like fibers comprising living micro-organisms (fungi and bacteria) and DND particles, with elevated nitrogen content due to fixation of atmospheric N2. Prolonged contact between these solutions and air lead to the formation of ammonium nitrate, sometimes forming macroscopic crystals. The latter was also formed during prolonged oxidation of fungi in aqueous DND colloids. The possible mechanism of N2 fixation is discussed, which can be attributable to the high reactivity of DND.

Orientational Ordering in Nano-confined Polar Liquids

Water and other polar liquids exhibit nanoscale structuring near charged interfaces. When a polar liquid is confined between two charged surfaces, the interfacial solvent layers begin to overlap, resulting in solvation forces. Here, we perform molecular dynamics simulations of polar liquids with different dielectric constants and molecular shapes and sizes confined between charged surfaces, demonstrating strong orientational ordering in the nanoconfined liquids. To rationalize the observed structures, we apply a coarse-grained continuum theory that captures the orientational ordering and solvation forces of those liquids. Our findings reveal the subtle behavior of different nanoconfined polar liquids and establish a simple law for the decay length of the interfacial orientations of the solvents, which depends on their molecular size and polarity. These insights shed light on the nature of solvation forces, which are important in colloid and membrane science, scanning probe microscopy, and nano-electrochemistry.

Distance-dependent dielectric constant at the calcite/electrolyte interface: Implication for surface complexation modeling

HYPOTHESIS

The electrical double layer formed at the mineral/electrolyte interface is often modeled using mean-field approaches based on a continuum description of the solvent whose dielectric constant is assumed to decrease monotonically with decreasing distance to the surface. In contrast, molecular simulations show that the solvent polarizability oscillates near the surface similar to the water density profile - as shown previously, for example, by Bonthuis et al. (D.J. Bonthuis, S. Gekle, R.R. Netz, Dielectric Profile of Interfacial Water and its Effect on Double-Layer Capacitance, Phys Rev Lett 107(16) (2011) 166102). We showed that molecular and mesoscale pictures agree by spatially averaging the dielectric constant obtained from molecular dynamics simulations over the distances relevant to the mean-field representation. In addition, the values of capacitances used to describe the electrical double layer in Surface Complexation Models (SCMs) of the mineral/electrolyte interface can be estimated using molecularly informed spatially averaged dielectric constants and positions of hydration layers.

EXPERIMENTS

First, we used molecular dynamics simulations to model the calcite 101¯4/electrolyte interface. Next, by using atomistic trajectories, we calculated the distance-dependent static dielectric constant and water density in the direction normal to the. Finally, we applied spatial compartmentalization consistent with the model of parallel-plate capacitors connected in series to estimate SCM capacitances.

FINDINGS

Computationally expensive simulations are required to determine the dielectric constant profile of interfacial water near the mineral surface. On the other hand, water density profiles are readily assessable from much shorter simulation trajectories. Our simulations confirmed that dielectric and water density oscillations at the interface are correlated. Here, we parametrized linear regression models to estimate the dielectric constant directly from the local water density. This is a significant computational shortcut compared to slowly converging calculations relying on total dipole moment fluctuations. The amplitude of the interfacial dielectric constant oscillation can exceed the dielectric constant of the bulk water, suggesting an ice-like frozen state, but only if there are no electrolyte ions. The interfacial accumulation of electrolyte ions causes a decrease in the dielectric constant due to the reduction of water density and re-orientation of water dipoles in ion hydration shells. Finally, we show how to use the computed dielectric properties to estimate SCM's capacitances.

Fluids and Electrolytes under Confinement in Single-Digit Nanopores

Confined fluids and electrolyte solutions in nanopores exhibit rich and surprising physics and chemistry that impact the mass transport and energy efficiency in many important natural systems and industrial applications. Existing theories often fail to predict the exotic effects observed in the narrowest of such pores, called single-digit nanopores (SDNs), which have diameters or conduit widths of less than 10 nm, and have only recently become accessible for experimental measurements. What SDNs reveal has been surprising, including a rapidly increasing number of examples such as extraordinarily fast water transport, distorted fluid-phase boundaries, strong ion-correlation and quantum effects, and dielectric anomalies that are not observed in larger pores. Exploiting these effects presents myriad opportunities in both basic and applied research that stand to impact a host of new technologies at the water-energy nexus, from new membranes for precise separations and water purification to new gas permeable materials for water electrolyzers and energy-storage devices. SDNs also present unique opportunities to achieve ultrasensitive and selective chemical sensing at the single-ion and single-molecule limit. In this review article, we summarize the progress on nanofluidics of SDNs, with a focus on the confinement effects that arise in these extremely narrow nanopores. The recent development of precision model systems, transformative experimental tools, and multiscale theories that have played enabling roles in advancing this frontier are reviewed. We also identify new knowledge gaps in our understanding of nanofluidic transport and provide an outlook for the future challenges and opportunities at this rapidly advancing frontier.

Quantifying the Molecular Polarization Response of Liquid Water Interfaces at Heterogeneously Charged Surfaces

The hydration shells of proteins mediate interactions, such as small molecule binding, that are vital to their biological function or in some cases their dysfunction. However, even when the structure of a protein is known, the properties of its hydration environment cannot be easily predicted due to the complex interplay between protein surface heterogeneity and the collective structure of water's hydrogen bonding network. This manuscript presents a theoretical study of the influence of surface charge heterogeneity on the polarization response of the liquid water interface. We focus our attention on classical point charge models of water, where the polarization response is limited to molecular reorientation. We introduce a new computational method for analyzing simulation data that is capable of quantifying water's collective polarization response and determining the effective surface charge distribution of hydrated surfaces over atomistic length scales. To illustrate the utility of this method, we present the results of molecular dynamics simulations of liquid water in contact with a heterogeneous model surface and the CheY protein.

Frequency-Dependent Impedance of Nanocapacitors from Electrode Charge Fluctuations as a Probe of Electrolyte Dynamics

The frequency-dependent impedance is a fundamental property of electrical components. We show that it can be determined from the equilibrium dynamical fluctuations of the electrode charge in constant-potential molecular simulations, extending in particular a fluctuation-dissipation relation for the capacitance recovered in the low-frequency limit and provide an illustration on water-gold nanocapacitors. This Letter opens the way to the interpretation of electrochemical impedance measurements in terms of microscopic mechanisms, directly from the dynamics of the electrolyte, or indirectly via equivalent circuit models as in experiments.

Electrostatics of soft (bio)interfaces: Corrections of mean-field Poisson-Boltzmann theory for ion size, dielectric decrement and ion-ion correlation

HYPOTHESIS

Electrostatics of soft (ion-permeable) (bio)particles (e.g. microorganisms, core/shell colloids) in aqueous electrolytes is commonly formulated by the mean-field Poisson-Boltzmann theory and integration of the charge contributions from electrolyte ions and soft material. However, the effects connected to the size of the electrolyte ions and that of the structural charges carried by the particle, to dielectric decrement and ion-ion correlations on soft interface electrostatics have been so far considered at the margin, despite the limits of the Gouy theory for condensed and/or multivalent electrolytes.

EXPERIMENTS

Accordingly, we modify herein the Poisson-Boltzmann theory for core/shell (bio)interfaces to include the aforementioned molecular effects considered separately or concomitantly. The formalism is applicable for poorly to highly charged particles in the thin electric double layer regime and to unsymmetrical multivalent electrolytes.

FINDINGS

Computational examples of practical interests are discussed with emphasis on how each considered molecular effect or combination thereof affects the interfacial potential distribution depending on size and valence of cations and anions, size of particle charges, length scale of ionic correlations and shell-to-Debye layer thickness ratio. The origins of here-evidenced pseudo-harmonic potential profile and ion size-dependent screening of core/shell particle charges are detailed. In addition, the existence and magnitude of the Donnan potential when reached in the shell layer are shown to depend on the excluded volumes of the electrolyte ions.

Coupled Interactions at the Ionic Graphene-Water Interface

We compute ionic free energy adsorption profiles at an aqueous graphene interface by developing a self-consistent approach. To do so, we design a microscopic model for water and put the liquid on an equal footing with the graphene described by its electronic band structure. By evaluating progressively the electronic and dipolar coupled electrostatic interactions, we show that the coupling level including mutual graphene and water screening permits one to recover remarkably the precision of extensive quantum simulations. We further derive the potential of mean force evolution of several alkali cations.

Ionizable lipids in bio-inspired nanocarriers

In applications of bio-inspired nanoparticles (NPs), their composition is often optimised by including ionizable lipids. I use a generic statistical model to describe the charge and potential distributions in lipid nanoparticles (LNPs) containing such lipids. The LNP structure is considered to contain the biophase regions separated by narrow interphase boundaries with water. Ionizable lipids are uniformly distributed at the biophase-water boundaries. The potential is there described at the mean-filed level combining the Langmuir-Stern equation for ionizable lipids and the Poisson-Boltzmann equation for other charges in water. The latter equation is used outside a LNP as well. With physiologically reasonable parameters, the model predicts the scale of the potential in a LNP to be rather low, smaller or about [Formula: see text], and to change primarily near the LNP-solution interface or, more precisely, inside an NP near this interface because the charge of ionizable lipids becomes rapidly neutralized along the coordinate towards the center of a LNP. The extent of dissociation-mediated neutralization of ionizable lipids along this coordinate increases but only slightly. Thus, the neutralization is primarily due to the negative and positive ions related to the ionic strength in solution and located inside a LNP.

Does twist angle affect the properties of water confined inside twisted bilayer graphene?

Graphene nanoslit pores are used for nanofluidic devices, such as, in water desalination, ion-selective channels, ionic transistors, sensing, molecular sieving, blue energy harvesting, and protein sequencing. It is a strenuous task to prepare nanofluidic devices, because a small misalignment leads to a significant alteration in various properties of the devices. Here, we focus on the rotational misalignment between two parallel graphene sheets. Using molecular dynamics simulation, we probe the structure and dynamics of monolayer water confined inside graphene nanochannels for a range of commensurate twist angles. With SPC/E and TIP4P/2005 water models, our simulations reveal the independence of the equilibrium number density- n ∼ 13 nm^-2 for SPC/E and n ∼ 11.5 nm^-2 for TIP4P/2005- across twists. Based on the respective densities of the water models, the structure and dielectric constant are invariant of twist angles. The confined water structure at this density shows square ice ordering for SPC/E water only. TIP4P/2005 shows ordering at the vicinity of a critical density (n ∼ 12.5 nm^-2). The average perpendicular dielectric constant of the confined water remains anomalously low (∼2 for SPC/E and ∼6 for TIP4P/2005) for the studied twist angles. We find that the friction coefficient of confined water molecules varies for small twist angles, while becoming independent for twists greater than 5.1°. Our results indicate that a small, angular misalignment will not impair the dielectric properties of monolayer water within a graphene slit-pore, but can significantly influence its dynamics.

Electrokinetic, electrochemical, and electrostatic surface potentials of the pristine water liquid-vapor interface

Although conceptually simple, the air-water interface displays rich behavior and is subject to intense experimental and theoretical investigations. Different definitions of the electrostatic surface potential as well as different calculation methods, each relevant for distinct experimental scenarios, lead to widely varying potential magnitudes and sometimes even different signs. Based on quantum-chemical density-functional-theory molecular dynamics (DFT-MD) simulations, different surface potentials are evaluated and compared to force-field (FF) MD simulations. As well explained in the literature, the laterally averaged electrostatic surface potential, accessible to electron holography, is dominated by the trace of the water molecular quadrupole moment, and using DFT-MD amounts to +4.35 V inside the water phase, very different from results obtained with FF water models which yield negative values of the order of -0.4 to -0.6 V. Thus, when predicting potentials within water molecules, as relevant for photoelectron spectroscopy and non-linear interface-specific spectroscopy, DFT simulations should be used. The electrochemical surface potential, relevant for ion transfer reactions and ion surface adsorption, is much smaller, less than 200 mV in magnitude, and depends specifically on the ion radius. Charge transfer between interfacial water molecules leads to a sizable surface potential as well. However, when probing electrokinetics by explicitly applying a lateral electric field in DFT-MD simulations, the electrokinetic ζ-potential turns out to be negligible, in agreement with predictions using continuous hydrodynamic models. Thus, interfacial polarization charges from intermolecular charge transfer do not lead to significant electrokinetic mobility at the pristine vapor-liquid water interface, even assuming these transfer charges are mobile in an external electric field.

Effects of surface rigidity and metallicity on dielectric properties and ion interactions at aqueous hydrophobic interfaces

Using classical molecular dynamics simulations we investigate the dielectric properties at interfaces of water with graphene, graphite, hexane and water vapor. For graphite we compare metallic and non-metallic versions. At the vapor-liquid water and hexane-water interfaces the laterally averaged dielectric profiles are significantly broadened due to interfacial roughness and only slightly anisotropic. In contrast, at the rigid graphene surface the dielectric profiles are strongly anisotropic and the perpendicular dielectric profile exhibits pronounced oscillations and sign changes. The interfacial dielectric excess, characterized by the shift of the dielectric-dividing-surface with respect to the Gibbs-dividing-surface, is positive for all surfaces, showing that water has an enhanced dielectric response at hydrophobic surfaces. The dielectric-dividing-surface positions vary significantly among the different surfaces, which points to pronounced surface-specific dielectric behavior. The interfacial repulsion of a chloride ion is shown to be dominated by electrostatic interactions for the soft fluid-fluid interfaces and by non-electrostatic Lennard-Jones interactions for the rigid graphene-water interface. A linear tensorial dielectric model for the ion-interface interaction with sharp dielectric interfaces located on the dielectric-dividing-surface positions works well for graphene but fails for vapor and hexane, because these interfaces are smeared out. The repulsion of chloride from the metallic and non-metallic graphite versions differs very little, which reflects the almost identical interfacial water structure and can be understood based on linear continuum dielectric theory. Interface flexibility shows up mostly in the non-linear Coulomb part of the ion-interface interaction, which changes significantly close to the interfaces and signals the breakdown of linear dielectric continuum theory.

The Role of Surface Chemistry in the Orientational Behavior of Water at an Interface

Molecular dynamics studies have demonstrated that molecular water at an interface, with either a gas or a solid, displays anisotropic orientational behavior in contrast to its bulk counterpart. This effect has been recently implicated in the like-charge attraction problem for colloidal particles in solution. Here, negatively charged particles in solution display a long-ranged attraction where continuum electrostatic theory predicts monotonically repulsive interactions, particularly in solutions with monovalent salt ions at low ionic strength. Anisotropic orientational behavior of solvent molecules at an interface gives rise to an excess interfacial electrical potential which we suggest generates an additional solvation contribution to the total free energy that is traditionally overlooked in continuum descriptions of interparticle interactions in solution. In the present investigation we perform molecular dynamics simulation based calculations of the interfacial potential using realistic surface models representing various chemistries as well as different solvents. Similar to previous work that focused on simple model surfaces constructed by using oxygen atoms, we find that solvents at more realistic model surfaces exhibit substantial anisotropic orientational behavior. We explore the dependence of the interfacial solvation potential on surface properties such as surface group chemistry and group density at silica and carboxylated polystyrene interfaces. For water, we note surprisingly good agreement between results obtained for a simple O-atom wall and more complex surface models, suggesting a general qualitative consistency of the interfacial solvation effect for surfaces in contact with water. In contrast, for an aprotic solvent such as DMSO, surface chemistry appears to exert a stronger influence on the sign and magnitude of the interfacial solvation potential. The study carries broad implications for molecular-scale interactions and may find relevance in explaining a range of phenomena in soft-matter physics and cell biology.

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

Electrically Suppressed Outflow of Confined Liquid in Hydrophobic Nanopores

Confining liquid in a hydrophobic nanoenvironment has enabled a broad spectrum of applications in biomedical sensors, mechanical actuators, and energy storage and converters, where the outflow of confined liquid is spontaneous and fast due to the intrinsic hydrophobic nature of nanopores with extremely low interfacial friction, challenging design capacity and control tolerance of structures and devices. Here, we present a facile approach of suppressing the outflow of water confined in hydrophobic nanopores with an electric field. Extensive molecular dynamics simulations show that the presence of an electric field could significantly strengthen hydrogen bonds and retard degradations of the associated networks during the outflow. The outflow deformation and strength are extracted to quantitatively characterize the electrical suppression to outflow and agree well with simulations. This study proposes a practical means of impeding the fast liquid outflow in hydrophobic nanopores, potentially useful for devising nanofluidics-based functional structures and devices with controllable performance.

Microscopic Simulations of Electrochemical Double-Layer Capacitors

Electrochemical double-layer capacitors (EDLCs) are devices allowing the storage or production of electricity. They function through the adsorption of ions from an electrolyte on high-surface-area electrodes and are characterized by short charging/discharging times and long cycle-life compared to batteries. Microscopic simulations are now widely used to characterize the structural, dynamical, and adsorption properties of these devices, complementing electrochemical experiments and in situ spectroscopic analyses. In this review, we discuss the main families of simulation methods that have been developed and their application to the main family of EDLCs, which include nanoporous carbon electrodes. We focus on the adsorption of organic ions for electricity storage applications as well as aqueous systems in the context of blue energy harvesting and desalination. We finally provide perspectives for further improvement of the predictive power of simulations, in particular for future devices with complex electrode compositions.

Electrokinetic behavior of artificial and natural calcites: A review of experimental measurements and surface complexation models

The surface charge of calcite in aqueous environments is essential to many industrial and environmental applications. Electrokinetic measurements are usually used to assess the calcite charging behavior and characterize its electrical double layer (EDL). Numerous surface complexation models (SCMs) have been proposed to interpret the effect of different surface interactions on the zeta potential. Because of their versatility, SCMs have also become important tools in reactive transport modeling. The research on enhanced oil recovery within the last decade has led to an increased number of publications reporting both zeta potential measurements and SCMs for calcite. Nonetheless, the measurements are often inconsistent and the reasons for choosing one model over another are unclear. In this work, we review the models proposed for calcite and address their main differences. We first collect a large number of published zeta potential measurements and then we fit a Diffuse Layer, Basic Stern, and Charge-Distribution Multi-Site Complexation models to a selected reliable dataset. For each model, we maintain a similar number of adjustable parameters. After optimizing the parameters of the models, we systematically compare their prediction capabilities against data obtained in monovalent and divalent electrolyte systems containing calcium, magnesium, sulfate, or carbonate. We show that, often, the discrepancies between the models and the experimental data can be explained by different levels of disequilibrium. Nonetheless, assumptions used in the development of the models may significantly reduce their extrapolability to variable chemical conditions. The poor agreement between the models tuned to electrokinetic data with surface charge measurements and dynamic retention from single-phase flowthrough tests show that zeta potential may not be the best type of data to characterize ion binding at the calcite surface. Including the effect of mineral impurities and temperature on the calcite surface speciation and electrokinetic behavior prevail as main challenges for reactive transport modeling.

Nonmonotonic Pair Potentials in the Interaction of Like-Charged Objects in Solution

We consider the long-standing like-charge attraction problem, wherein under certain conditions, similarly charged spheres suspended in aqueous electrolyte have been observed to display a minimum in their interaction potential, contrary to the intuitively expected monotonically varying repulsion. Recently, we described an interfacial mechanism invoking the molecular nature of the solvent that explains this anomalous experimental observation. In our model for the interaction of negatively charged particles in water, the minimum in the pair potential results from the superposition of competing contributions to the total free energy. One of these contributions is the canonical repulsive electrostatic term, whereas the other is a solvation-induced attractive contribution. We find that whereas both contributions grow approximately exponentially with decreasing interparticle separation, the occurrence of a stable, long-ranged minimum in the pair potential arises from differences in the precise interparticle separation dependence of the two terms. Specifically, the interfacial solvation term exhibits a more gradual decay with distance than the electrostatic repulsion, permitting the attractive contribution to dominate the interaction at large distances. Importantly, these disparities become evident in quantities calculated from exact numerical solutions to the governing nonlinear Poisson-Boltzmann (PB) equation for the spatial electrical potential distribution in the system. In marked contrast, we find that the linearized PB equation, applicable in the regime of low surface electrical potentials, does not support nonmonotonic trends in the total interaction free energy within the present model. Our results point to the importance of exact descriptions of electrostatic interactions in real systems that most often do not subscribe to particular mathematical limits where analytical approximations may provide a sufficiently accurate description of the problem.

Dielectric response of thin water films: a thermodynamic perspective

The surface of a polar liquid presents a special environment for the solvation and organization of charged solutes, which differ from bulk behaviors in important ways. These differences have motivated many attempts to understand electrostatic response at aqueous interfaces in terms of a spatially varying dielectric permittivity, typically concluding that the dielectric constant of interfacial water is significantly lower than in the bulk liquid. Such analyses, however, are complicated by the potentially nonlocal nature of dielectric response over the short length scales of interfacial heterogeneity. Here we circumvent this problem for thin water films by adopting a thermodynamic approach. Using molecular simulations, we calculate the solvent's contribution to the reversible work of charging a parallel plate capacitor. We find good agreement with a simple dielectric continuum model that assumes bulk dielectric permittivity all the way up to the liquid's boundary, even for very thin (∼1 nm) films. This comparison requires careful attention to the placement of dielectric boundaries between liquid and vapor, which also resolves apparent discrepancies with dielectric imaging experiments.

Free energy calculations from molecular simulations reveal that water's interfacial dielectric response is well-described by bulk properties.

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Implicit solvation is an effective, highly coarse-grained approach in atomic-scale simulations to account for a surrounding liquid electrolyte on the level of a continuous polarizable medium. Originating in molecular chemistry with finite solutes, implicit solvation techniques are now increasingly used in the context of first-principles modeling of electrochemistry and electrocatalysis at extended (often metallic) electrodes. The prevalent ansatz to model the latter electrodes and the reactive surface chemistry at them through slabs in periodic boundary condition supercells brings its specific challenges. Foremost this concerns the difficulty of describing the entire double layer forming at the electrified solid–liquid interface (SLI) within supercell sizes tractable by commonly employed density functional theory (DFT). We review liquid solvation methodology from this specific application angle, highlighting in particular its use in the widespread ab initio thermodynamics approach to surface catalysis. Notably, implicit solvation can be employed to mimic a polarization of the electrode’s electronic density under the applied potential and the concomitant capacitive charging of the entire double layer beyond the limitations of the employed DFT supercell. Most critical for continuing advances of this effective methodology for the SLI context is the lack of pertinent (experimental or high-level theoretical) reference data needed for parametrization.

Hydrophobicity of Self-Assembled Monolayers of Alkanes: Fluorination, Density, Roughness, and Lennard-Jones Cutoffs

The interplay of fluorination and structure of alkane self-assembled monolayers and how these affect hydrophobicity are explored via molecular dynamics simulations, contact angle goniometry, and surface-enhanced infrared absorption spectroscopy. Wetting coefficients are found to grow linearly in the monolayer density for both alkane and perfluoroalkane monolayers. The larger contact angles of monolayers of perfluorinated alkanes are shown to be primarily caused by their larger molecular volume, which leads to a larger nearest-neighbor grafting distance and smaller tilt angle. Increasing the Lennard-Jones force cutoff in simulations is found to increase hydrophilicity. Specifically, wetting coefficients scale like the inverse square of the cutoff, and when extrapolated to the infinite cutoff limit, they yield contact angles that compare favorably to experimental values. Nanoscale roughness is also found to reliably increase monolayer hydrophobicity, mostly via the reduction of the entropic part of the work of adhesion. Analysis of depletion lengths shows that droplets on nanorough surfaces partially penetrate the surface, intermediate between Wenzel and Cassie-Baxter states.

A molecular perspective on induced charges on a metallic surface

Understanding the response of the surface of metallic solids to external electric field sources is crucial to characterize electrode-electrolyte interfaces. Continuum electrostatics offer a simple description of the induced charge density at the electrode surface. However, such a simple description does not take into account features related to the atomic structure of the solid and to the molecular nature of the solvent and of the dissolved ions. In order to illustrate such effects and assess the ability of continuum electrostatics to describe the induced charge distribution, we investigate the behavior of a gold electrode interacting with sodium or chloride ions fixed at various positions, in a vacuum or in water, using all-atom constant-potential classical molecular dynamics simulations. Our analysis highlights important similarities between the two approaches, especially under vacuum conditions and when the ion is sufficiently far from the surface, as well as some limitations of the continuum description, namely, neglecting the charges induced by the adsorbed solvent molecules and the screening effect of the solvent when the ion is close to the surface. While the detailed features of the charge distribution are system-specific, we expect some of our generic conclusions on the induced charge density to hold for other ions, solvents, and electrode surfaces. Beyond this particular case, the present study also illustrates the relevance of such molecular simulations to serve as a reference for the design of improved implicit solvent models of electrode-electrolyte interfaces.

Probing the size-dependent polarizability of mesoscopic ionic clusters and their induced-dipole interactions

Mesoscopic clusters composed of oppositely charged particles are ubiquitous in synthetic and biological soft materials. The effective interaction between these clusters is influenced by their polarizability, that is, the ability of their constituent charges to re-arrange in response to an external electrical field. Here, using coarse-grained simulations, we show that the polarizability of electrically neutral ionic clusters decreases as the number of constituent charges increases and/or their Coulombic interaction strength increases for various ion valencies, ion densities, and degrees of cluster boundary hardness. For clusters of random ionomers and their counterions, their polarizability is shown to depend on the number of polymer chains. The variation of the cluster polarizability with the cluster size indicates that throughout the assembly, the induced-dipole interactions between the clusters may be reduced substantially as they acquire more charges while maintaining zero net charge. Under certain conditions, the induced-dipole interactions may become repulsive, as inferred from our simulations with a polarizable solvent. As a result, the dipole-induced related interactions can serve as a counterbalancing force that contributes to the self-limiting aggregation of charge-containing assemblies.

Competition between Born solvation, dielectric exclusion, and Coulomb attraction in spherical nanopores

The recent measurement of a very low dielectric constant, ε, of water confined in nanometric slit pores leads us to reconsider the physical basis of ion partitioning into nanopores. For confined ions in chemical equilibrium with a bulk of dielectric constant ε_{b}>ε, three physical mechanisms, at the origin of ion exclusion in nanopores, are expected to be modified due to this dielectric mismatch: dielectric exclusion at the water-pore interface (with membrane dielectric constant, ε_{m}<ε), the solvation energy related to the difference in Debye-Hückel screening parameters in the pore, κ, and in the bulk κ_{b}, and the classical Born solvation self-energy proportional to ε^{-1}-ε_{b}^{-1}. Our goal is to clarify the interplay between these three mechanisms and investigate the role played by the Born contribution in ionic liquid-vapor (LV) phase separation in confined geometries. We first compute analytically the potential of mean force (PMF) of an ion of radius R_{i} located at the center of a nanometric spherical pore of radius R. Computing the variational grand potential for a solution of confined ions, we then deduce the partition coefficients of ions in the pore versus R and the bulk electrolyte concentration ρ_{b}. We show how the ionic LV transition, directly induced by the abrupt change of the dielectric contribution of the PMF with κ, is favored by the Born self-energy and explore the decrease of the concentration in the pore with ε both in the vapor and liquid states. Phase diagrams are established for various parameter values and we show that a signature of this phase transition can be detected by monitoring the total osmotic pressure as a function of R. For charged nanopores, these exclusion effects compete with the electrostatic attraction that imposes the entry of counterions into the pore to enforce electroneutrality. This study will therefore help in deciphering the respective roles of the Born self-energy and dielectric mismatch in experiments and simulations of ionic transport through nanopores.

Weakly Anisotropic Dielectric Properties of Water Droplets at the Nanoscale

The slab-confined water at the nanoscale exhibits anomalous dielectric properties compared to bulk water, for example, significantly low dielectric constant. In this work, we study the dielectric properties of nanoscale water droplets at room temperature using molecular dynamics simulations. We find that the nanoscale water droplets feature weakly anisotropic dielectric constant: the radial component of dielectric constants is distinctly smaller than the tangential component although they both decrease with reducing droplet size in a similar way. Such dielectric behavior is closely related to the orientational preference of water molecules near the convex surface. The molecular dipole prefers to slightly orientate toward the interior of droplets in contrast to the out-of-plane preference for free-standing water films and slab-confined water, which suppresses the fluctuation of dipole moments in the radial direction. Meanwhile, it facilitates the formation of the open hydrogen-bond network in the surface layer and ultimately leads to the relatively weak suppression of tangential fluctuations. The differential suppression is responsible for the anisotropic dielectric constant of water droplets. This anisotropic characteristic is also found in dielectric relaxation: both the radial and the tangential relaxation are consistently slowed down upon approaching surface but the latter is universally slower.

Effects of electrical image potentials and solvation energy on salt ions near a metallic or dielectric wall

Electrical image potentials near a metallic or a dielectric wall of higher dielectric constant than that of the solution are attractive, and therefore, could concentrate salt ions near the wall. In fact, ions in room temperature ionic liquids have been observed to precipitate near a metallic surface (but not near a nonmetallic surface). It will be argued that a likely reason for why precipitation of ions in salt water due to electrical image forces has not as yet been observed is that the solvation of the ions is reduced near the wall. This results in an energy barrier. This reduction occurs because of the large decrease near the wall of the dielectric constant of water normal to the wall. The conditions under which ions are able to get past the resulting energy barrier and concentrate at a solid wall, either as a result of a reduction in this barrier due to screening at high ion concentration or as a result of thermal activation over the barrier will be explored.

Electrophoresis of Liquid-Layer Coated Particles: Impact of Ion Partitioning and Ion Steric Effects

The biomimetic core-shell nanoparticles coated with membranes of various biological cells have attracted significant research interest, because of their extensive applications in targeted drug delivery systems. The cell membrane consists of a lipid bilayer, which can be regarded as a two-dimensional oriented viscous liquid with low dielectric permittivity, compared to a bulk aqueous medium. Such a liquid layer comprised of cell membrane may bear additional mobile charges, because of the presence of free lipid molecules or charged surfactant molecules, which further results in nonzero charge along the surface of the peripheral layer. In this article, we present an analytical theory for electrophoresis of such cell membrane coated functionalized nanoparticles in the extent of electrolyte solution, considering the combined effects of finite ion size and of ion partitioning. Going beyond the Debye-Huckel approximations, we propose an analytical theory for Donnan potential and electrophoretic mobility. The derived expressions are applicable for moderate to highly charged undertaken core-shell particles when the thickness of the peripheral liquid layer greatly exceeds the electric double layer thickness. The impact of pertinent parameters on the electrophoretic response of such a particle is further discussed.

Dielectric Susceptibility of Water in the Interface

It has long been anticipated that dielectric constants of polar liquids are reduced in the interfacial layer. Recent experiments and computer simulations support these expectations. A strong reduction of the dielectric constant is found in the direction perpendicular to a planar substrate, while the parallel response is bulk-like. This Perspective highlights recent theoretical calculations and simulations with an eye on relating them to properties observable in the laboratory. The average interface dielectric constant computed from simulations connects to thin films experiments, but this cannot be extended to screening of charges. In contrast to dielectric theories where a single dielectric constant gauges both the polarization energy and screening, these two signatures of dielectric polarization diverge on the molecular scale. The reduction of the dielectric constant of water in thin films is currently viewed as a combined effect of geometric confinement imposed by the substrate and the reconstruction of water hydrogen bonds in the surface layer.

Dipolar Poisson models in a dual view

In this work, we study the continuum theories of dipolar-Poisson models. Both the standard dipolar-Poisson model and the dipolar-Poisson-Langevin model, which keeps the dipolar density fixed, are non-convex functionals of the scalar electrostatic potential ϕ. Applying the Legendre transform approach introduced by Maggs [Europhys. Lett. 98, 16012 (2012)], the dual functionals of these models are derived and are given by convex vector-field functionals of the dielectric displacement D and the polarization field P. We compare the convex functionals in P-space to the non-convex functionals in electric field E-space and apply them to the classic problem of the solvation of point-like ions. Since the dipolar-Poisson model does not properly describe polarization saturation, we argue that only the dipolar-Poisson-Langevin functional can be used to provide a nonlinear generalization of the harmonic polarization functional used in the theory of Marcus for the electron transfer rate to nonlinear regimes. We show that the model can be quantitatively parameterized by molecular dynamics simulations.

Dielectric Profile and Electromelting of a Monolayer of Water Confined in Graphene Slit Pore

A monolayer of water confined between two parallel graphene sheets exists in many different phases and exhibits fascinating dielectric properties that have been studied in experiments. In this work, we use molecular dynamics simulations to study how the dielectric properties of a confined monolayer of water is affected by its structure. We consider six of the popular nonpolarizable water models-SPC/E, SPC/Fw, TIP3P, TIP3P_M (modified), TIP4P-2005, and TIP4P-2005f-and find that the in-plane structure of the water molecules at ambient temperature and pressure is strongly dependent on the water model: all the 3-point water models considered here show square ice formation, whereas no such structural ordering is observed for the 4-point water models. This allows us to investigate the role of the in-plane structure of the water monolayer on its dielectric profile. Our simulations show an anomalous perpendicular dielectric constant compared to the bulk, and the models that do not exhibit ice formation show very different dielectric response along the channel width compared to models that exhibit square ice formation. We also demonstrate the occurrence of electromelting of the in-plane ordered water under the application of a perpendicular electric field and find that the critical field for electromelting strongly depends on the water model. Together, we have shown the dependence of confined water properties on the different water structures that it may take when sandwiched between bilayer graphene. These remarkable properties of confined water can be exploited in various nanofluidic devices, artificial ion channels, and molecular sieving.

Confinement-Controlled Aqueous Chemistry within Nanometric Slit Pores

In this Focus Review, we put the spotlight on very recent insights into the fascinating world of wet chemistry in the realm offered by nanoconfinement of water in mechanically rather rigid and chemically inert planar slit pores wherein only monolayer and bilayer water lamellae can be hosted. We review the effect of confinement on different aspects such as hydrogen bonding, ion diffusion, and charge defect migration of H⁺(aq) and OH^-(aq) in nanoconfined water depending on slit pore width. A particular focus is put on the strongly modulated local dielectric properties as quantified in terms of anisotropic polarization fluctuations across such extremely confined water films and their putative effects on chemical reactions therein. The stunning findings disclosed only recently extend wet chemistry in particular and solvation science in general toward extreme molecular confinement conditions.

Nonlocal Dielectric Response of Water in Nanoconfinement

Recent experiments reporting a very low dielectric permittivity for nanoconfined water have renewed the interest in the structure and dielectric properties of water in narrow gaps. Here, we describe such systems with a minimal Landau-Ginzburg field theory composed of a nonlocal bulk-determined term and a local water-surface interaction term. We show how the interplay between the boundary conditions and intrinsic bulk correlations encodes the dielectric properties of confined water. Our theoretical analysis is supported by molecular dynamics simulations and comparison with the experimental data.

Confined Water's Dielectric Constant Reduction Is Due to the Surrounding Low Dielectric Media and Not to Interfacial Molecular Ordering

Liquid water confined within nanometer-sized channels exhibits a surprisingly low dielectric constant along the direction orthogonal to the channel walls. This is typically assumed to result from a pronounced heterogeneity across the sample: the dielectric constant would be bulk-like everywhere except at the interface, where it would be dramatically reduced by strong restrictions on interfacial molecules. Here we study the dielectric properties of water confined within graphene slit channels via classical molecular dynamics simulations. We show that the permittivity reduction is not due to any important alignment of interfacial water molecules, but instead to the long-ranged anisotropic dipole correlations combined with an excluded-volume effect of the low-dielectric confining material. The bulk permittivity is gradually recovered only over several nanometers due to the impact of long-range electrostatics, rather than structural features. This has important consequences for the control of, e.g., ion transport and chemical reactivity in nanoscopic channels and droplets.

Interfacial, Electroviscous, and Nonlinear Dielectric Effects on Electrokinetics at Highly Charged Surfaces

The dielectric constant and the viscosity of water at the interface of hydrophilic surfaces differ from their bulk values, and it has been proposed that the deviation is caused by the strong electric field and the high ion concentration in the interfacial layer. We calculate the dependence of the dielectric constant and the viscosity of bulk electrolytes on the electric field and the salt concentration. Incorporating the concentration and field-dependent dielectric constant and viscosity in the extended Poisson–Boltzmann and Stokes equations, we calculate the electro-osmotic mobility. We compare the results to literature experimental data and explicit molecular dynamics simulations of OH-terminated surfaces and show that it is necessary to additionally include the presence of a subnanometer wide interfacial water layer, the properties of which are drastically transformed by the sheer presence of the interface. We conclude that the origin of the anomalous behavior of aqueous interfacial layers cannot be found in electrostriction or electroviscous effects caused by the interfacial electric field and ion concentration. Instead, it is primarily caused by the intrinsic ordering and orientation of the interfacial water layer.

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS₂ nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS₂ channel can be significantly altered. We found that larger cations with weaker hydration energy (K⁺, Rb⁺, and Cs⁺) exhibited a distinctive CN decrease under confinement, while smaller cations (Li⁺ and Na⁺) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K⁺, Rb⁺, and Cs⁺) in MoS₂ showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS₂ channel.

Molecular Simulation of Electrode-Solution Interfaces

Many key industrial processes, from electricity production, conversion, and storage to electrocatalysis or electrochemistry in general, rely on physical mechanisms occurring at the interface between a metallic electrode and an electrolyte solution, summarized by the concept of an electric double layer, with the accumulation/depletion of electrons on the metal side and of ions on the liquid side. While electrostatic interactions play an essential role in the structure, thermodynamics, dynamics, and reactivity of electrode-electrolyte interfaces, these properties also crucially depend on the nature of the ions and solvent, as well as that of the metal itself. Such interfaces pose many challenges for modeling because they are a place where quantum chemistry meets statistical physics. In the present review, we explore the recent advances in the description and understanding of electrode-electrolyte interfaces with classical molecular simulations, with a focus on planar interfaces and solvent-based liquids, from pure solvent to water-in-salt electrolytes.

Abnormal in-plane permittivity and ferroelectricity of confined water: From sub-nanometer channels to bulk

Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (ε_∥) reaching values of about 120 for channels with height 8 Å < h < 10 Å. With the increase in the width of the channel, we predict that ε_∥ decreases nonlinearly and reaches the bulk value for h > 70 Å. A stratified continuum model is proposed that reproduces the h > 10 Å dependence of ε_∥. For sub-nanometer height channels, abnormal behavior of ε_∥ is found with two orders of magnitude reduction of ε_∥ around h ∼ 7.5 Å, which is attributed to the formation of a particular ice phase that exhibits long-time (∼μs) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.

Static dielectric permittivity of ionic liquids ultraconfined in carbon nanotubes

In this work the parallel component of the static dielectric permittivity, ϵ II of ionic liquids ultraconfined into flexible carbon nanotubes of radius of 1.2 nm and 2.4 nm is evaluated from molecular dynamics simulations. We show an enhancement of ϵ II with respect to bulk value and a counter-intuitive temperature dependence. Indeed an increase of ϵ II as a function of the temperature opposed to a bulk behavior is evidenced. Increase in ϵ II is the result of the strong orientation of ionic liquid close to the pore wall. The temperature dependence is the consequence of the thermal fluctuations increasing the dipolar fluctuations such that the strong orientation is conserved. Eventually, we show a molecular stacking between [C4mim+][Tf2N−] and CNT decreasing dipolar fluctuations close to the CNT surface reducing ϵ II .

Dielectric Behavior near a Spherical Ion

The dielectric response of a polar solvent to an ion is analyzed in terms of the bound charge, the net charge that accumulates near the ion as a consequence of the inhomogeneous polarization of the surrounding solvent. We demonstrate that the total bound charge arising in a full molecular treatment is identical to the total bound charge from standard continuum theory. In continuum theory, the bound charge resides in an infinitely thin layer, while in a molecular description the bound charge is spread over a region of finite width. Near simple atomic ions, the width of the bound charge distribution is roughly 1.3 nm. By simulating a sequence of ion charges from 0.1 to 2 e, where e is the magnitude of the electron charge, we analyze the applicability of linear response theory, which has been used by several authors. With increasing charge, the nonlinear response extends to an increasing distance from the ion. However, outside the region containing bound charge, the response is linear and in accord with continuum theory. Previous attempts to assign a dielectric constant to a solvent in the interfacial region are analyzed.

Anomalously Low Dielectric Constant of Ordered Interfacial Water

Despite considerable effort, the dielectric constant of interfacial water at solid surfaces is still not fully understood, thus hindering our understanding of the ubiquitous physical interactions in many materials and biological surfaces. In this study, we used molecular dynamics simulations to show that the parallel dielectric constant at the solid/water interface depends on solid-water interactions as well as the interfacial water structure on various solid crystal faces. In particular, ordered water structures can lead to a significant reduction (∼44%) in the parallel dielectric constant at the solid/water interface compared with that of bulk water. This sharp decrease in the parallel dielectric constant can be attributed to the specific antiferroelectric ordered structure of interfacial water molecules, which significantly suppresses the amplitude of the dipolar fluctuation associated with both the number of hydrogen bonds and the degree of order of interfacial water.

Uncovering a Universal Molecular Mechanism of Salt Ion Adsorption at Solid/Water Interfaces

Solid/water interfaces, in which salt ions come in close proximity to solids, are ubiquitous in nature. Because water is a polar solvent and salt ions are charged, a long-standing puzzle involving solid/water interfaces is how do the electric fields exerted by the salt ions and the interfacial water molecules polarize the charge distribution in the solid and how does this polarization, in turn, influence ion adsorption at any solid/water interface. Here, using state-of-the-art polarizable force fields derived from quantum chemical simulations, we perform all-atomistic molecular dynamics simulations to investigate the adsorption of various ions comprising the well-known Hofmeister series at the graphene/water interface, including comparing with available experimental data. Our findings reveal that, in vacuum, the ionic electric field-induced polarization of graphene results in a significantly large graphene-ion polarization energy, which drives all salt ions to adsorb to graphene. On the contrary, in the presence of water molecules, we show that the ions and the water molecules exert waves of molecular electric fields on graphene which destructively interfere with each other. This remarkable phenomenon is shown to cause a water-mediated screening of more than 85% of the graphene-ion polarization energy. Finally, by investigating superhydrophilic and superhydrophobic model surfaces, we demonstrate that this phenomenon occurs universally at all solid/water interfaces and results in a significant weakening of the ion-solid interactions, such that ion specific effects are governed primarily by a competition between the ion-water and water-water interactions, irrespective of the nature of the solid/water interface.

How the hydroxylation state of the (110)-rutile TiO2 surface governs its electric double layer properties

The hydroxylation state of an oxide surface is a central property of its solid/liquid interface and its corresponding electrical double layer. This study integrated both a reactive force field (ReaxFF) and a non-reactive potential into a hierarchical framework within molecular dynamics (MD) simulations to reveal how the hydroxylation state of the (110)-rutile TiO2 surface affects the electrical double layer properties. The simulation results obtained in the ReaxFF framework have shown that, while water dissociation occurs only at the under-coordinated Ti5c sites on the pristine TiO2 surface, the presence of point defects on the surface facilitates water dissociation at the oxygen vacancy sites, leading to two protonated oxygen bridge atoms for each vacancy site. As a consequence of enhanced water dissociation at the vacancy sites, water dissociation is quenched at the under-coordinated Ti5c sites resulting in two competitive hydroxylation mechanisms on the (110)-TiO2 surface. Using non-reactive MD simulations with hydroxylation states derived from the ReaxFF analysis, we demonstrate that water dissociation at the vacancy sites is a central mechanism governing the structuring of water near the interface. While the structuring of water near the interface is the main contribution to the electric field, water dissociation at the vacancy site enhances the adsorption of the electrolyte ions at the interface. The adsorbed ions lead to an increase of the effective surface charge as well as surface (zeta) potentials which are in the range of experimental observations. Our work provides a hierarchical multiscale simulation approach, covering a series of results with in-depth discussion for atomic/molecular level understanding of water dissociation and its effect on electric double layer properties of TiO2 to advance water splitting.