The Surface Potential of the Water–Vapor Interface from Classical Simulations

Interfacial Effects of Catalysis in Pickering Emulsions

Liquid-liquid or gas-liquid interfaces are ubiquitous in nature and in industrial production. Understanding the unique effects arising from the asymmetric interfaces and controlling the catalytic reactions are frontiers of physical chemistry. However, our understanding of the reactivity and selectivity at the interfaces remains scant. Pickering emulsions are emerging as a stable biphasic reaction system, which provides a new opportunity for clarifying the inherent features responsible for prominent interfacial reactivity or selectivity. This Perspective tentatively discusses the unique effects of interfacial adsorption, hydrogen bonding of water molecules, and strong electric field at the interfaces. Additionally, it highlights key insights into the fundamental mechanisms of reaction kinetic and thermodynamic alterations, molecular orientations, and the spontaneous generation of reactive species at the interfaces through representative examples. Finally, we delineate the current challenges and propose future research directions. The perspectives advanced herein may serve as valuable guidance for the design of efficient interfacial catalytic systems.

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}.

Identifying the acidic or basic behavior of surface water: a QM/MM-MD study

Controversies on the water surface were theoretically addressed with the help of large scale quantum mechanical molecular dynamics (QMMD) simulations on water surface model systems with and without excess hydroniums and hydroxides. It was revealed that the thermodynamic surface structures of these ions strongly depend on their location and dipole orientation. Fast hydronium diffusion by proton transfer establishes a wider kinetic depth distribution (∼6 Å) than that predicted by its thermodynamic affinity for the water surface, while slow hydroxide is shallowly trapped below the outermost molecular layer (3-4 Å). In addition, the anisotropic orientation of surface water dipole can generate a substantial magnitude of surface potential, which extends to a depth of a few molecular layers. With these distinctively different surface properties of two ions and water molecules, the seemingly contradictory observations of acidic and negatively charged water surfaces may be successfully explained. That is, the negative surface charge of neutral water mostly stems from intrinsic water properties such as water dipole orientation and electron density spillage at the surface, rather than surface OH- ions. The enhanced acidity of the water surface can be attributed in large part to the kinetic depth profile of ion density in addition to static thermodynamic origin. Furthermore, the different depth profiles of the two ions may differently affect the surface-sensitive spectroscopic observations.

Simultaneous and Spontaneous Oxidation and Reduction in Microdroplets by the Water Radical Cation/Anion Pair

Microdroplets show unique chemistry, especially in their intrinsic redox properties, and to this we here add a case of simultaneous and spontaneous oxidation and reduction. We report the concurrent conversions of several phosphonates to phosphonic acids by reduction (R-P → H-P) and to pentavalent phosphoric acids by oxidation. The experimental results suggest that the active reagent is the water radical cation/anion pair. The water radical cation is observed directly as the ionized water dimer while the water radical anion is only seen indirectly though the spontaneous reduction of carbon dioxide to formate. The coexistence of oxidative and reductive species in turn supports the proposal of a double-layer structure at the microdroplet surface, where the water radical cation and radical anion are separated and accumulated.

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.

Spontaneous Oxidation of Aromatic Sulfones to Sulfonic Acids in Microdroplets

Reactions in microdroplets can be accelerated and can present unique chemistry compared to reactions in bulk solution. Here, we report the accelerated oxidation of aromatic sulfones to sulfonic acids in microdroplets under ambient conditions without the addition of acid, base, or catalyst. The experimental data suggest that the water radical cation, (H₂O)^+•, derived from traces of water in the solvent, is the oxidant. The substrate scope of the reaction indicates the need for a strong electron-donating group (e.g., p-hydroxyl) in the aromatic ring. An analogous oxidation is observed in an aromatic ketone with benzoic acid production. The shared mechanism is suggested to involve field-assisted ionization of water at the droplet/air interface, its reaction with the sulfone (M) to form the radical cation adduct, (M + H₂O)^+•, followed by 1,2-aryl migration and C-O cleavage. A remarkably high reaction rate acceleration (∼10³) and regioselectivity (∼100-fold) characterize the reaction.

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.

Can electric fields drive chemistry for an aqueous microdroplet?

Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O-H bonds at the surface are ~16 MV/cm larger on average than that found for O-H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.

Divalent Ion Specific Outcomes on Stern Layer Structure and Total Surface Potential at the Silica:Water Interface

The second-order nonlinear susceptibility, χ⁽²⁾, in the Stern layer and the total interfacial potential drop, Φ(0)_tot, across the oxide:water interface are estimated from SHG amplitude and phase measurements for divalent cations (Mg²⁺, Ca²⁺, Sr²⁺, and Ba²⁺) at the silica:water interface at pH 5.8 and various ionic strengths. We find that interfacial structure and total potential depend strongly on ion valency. We observe statistically significant differences between the experimentally determined χ⁽²⁾ value for NaCl and that of the alkali earth series but smaller differences between ions of the same valency in that series. These differences are particularly pronounced at intermediate salt concentrations, which we attribute to the influence of hydration structure in the Stern layer. Furthermore, we corroborate the differences by examining the effects of anion substitution (SO₄^2- for Cl^-). Finally, we identify that hysteresis in measuring the reversibility of ion adsorption and desorption at fused silica in forward and reverse titrations manifests itself both in Stern layer structure and in total interfacial potential for some of the salts, most notably for CaCl₂ and MgSO₄ but less so for BaCl₂ and NaCl.

Accurate vertical ionization energy and work function determinations of liquid water and aqueous solutions

The absolute-scale electronic energetics of liquid water and aqueous solutions, both in the bulk and at associated interfaces, are the central determiners of water-based chemistry. However, such information is generally experimentally inaccessible. Here we demonstrate that a refined implementation of the liquid microjet photoelectron spectroscopy (PES) technique can be adopted to address this. Implementing concepts from condensed matter physics, we establish novel all-liquid-phase vacuum and equilibrated solution–metal-electrode Fermi level referencing procedures. This enables the precise and accurate determination of previously elusive water solvent and solute vertical ionization energies, VIEs. Notably, this includes quantification of solute-induced perturbations of water's electronic energetics and VIE definition on an absolute and universal chemical potential scale. Defining and applying these procedures over a broad range of ionization energies, we accurately and respectively determine the VIE and oxidative stability of liquid water as 11.33 ± 0.03 eV and 6.60 ± 0.08 eV with respect to its liquid-vacuum-interface potential and Fermi level. Combining our referencing schemes, we accurately determine the work function of liquid water as 4.73 ± 0.09 eV. Further, applying our novel approach to a pair of exemplary aqueous solutions, we extract absolute VIEs of aqueous iodide anions, reaffirm the robustness of liquid water's electronic structure to high bulk salt concentrations (2 M sodium iodide), and quantify reference-level dependent reductions of water's VIE and a 0.48 ± 0.13 eV contraction of the solution's work function upon partial hydration of a known surfactant (25 mM tetrabutylammonium iodide). Our combined experimental accomplishments mark a major advance in our ability to quantify electronic–structure interactions and chemical reactivity in liquid water, which now explicitly extends to the measurement of absolute-scale bulk and interfacial solution energetics, including those of relevance to aqueous electrochemical processes.

A generalised liquid-phase photoelectron spectroscopy approach is reported, allowing accurate, absolute energy scale ionisation energies of liquid water and aqueous solutions, as well as liquid water's work function to be reported.

A New Imaginary Term in the Second-Order Nonlinear Susceptibility from Charged Interfaces

Nonresonant second harmonic generation (SHG) phase and amplitude measurements obtained from the silica-water interface at varying pH values and an ionic strength of 0.5 M point to the existence of a nonlinear susceptibility term, which we call χ_X⁽³⁾, that is associated with a 90° phase shift. Including this contribution in a model for the total effective second-order nonlinear susceptibility produces reasonable point estimates for interfacial potentials and second-order nonlinear susceptibilities when χ_X⁽³⁾ ≈ 1.5χ_water⁽³⁾. A model without this term and containing only traditional χ⁽²⁾ and χ⁽³⁾ terms cannot recapitulate the experimental data. The new model also provides a demonstrated utility for distinguishing apparent differences in the second-order nonlinear susceptibility when the electrolyte is NaCl versus MgSO₄, pointing to the possibility of using heterodyne-detected SHG to investigate ion specificity in interfacial processes.

Charge State Dependence of Amino Acid Propensity at Water Surface: Mechanisms Elucidated by Molecular Dynamics Simulations

Atmospheric aerosols contain a variety of compounds, among them free amino acids and salt ions. The pH of the aerosol droplets depends on their origin and environment. Consequently, compounds like free amino acids found in the droplets will be at different charge states, since these states to a great extent depend on the surrounding pH condition. In droplets of marine origin, amino acids are believed to drive salt ions to the water surface and a pH-dependent amino acid surface propensity will, therefore, indirectly affect many processes in atmospheric chemistry and physics such as for instance cloud condensation. To understand the surface propensity of glycine, valine, and phenylalanine at acidic, neutral, and basic pH, we used molecular dynamics (MD) simulations to investigate them at three different charge states in water. Their respective surface propensities were obtained by the means of a potential of mean force (PMF) in an umbrella sampling approach. Glycine was found to have no preference for the surface, while both valine and phenylalanine showed high propensities. Among the charge states of the surface-enriched ones, the cation, representing the amino acids at low pH, was found to have the highest affinity. Free energy decomposition revealed that the driving forces depend strongly on the nature of the amino acid and its charge state. In phenylalanine, the main factor was found to be a substantial entropy gain, likely related to the side chain, whereas in valine, hydrogen bonding to the functional groups leads to favorable energies and, in turn, affects the surface propensity. A significant gain in water-water enthalpy was seen for both valine and phenylalanine.

Kosmotropic Electrolyte (Na2CO3, NaF) Perturbs the Air/Water Interface through Anion Hydration Shell without Forming a Well-Defined Electric Double Layer

The ion-driven electric double layer (EDL) and the structural transformation of interfacial water are implicated in unusual reaction kinetics at the air/water interface. By combining heterodyne-detected vibrational sum frequency generation (HD-VSFG) with differential spectroscopy involving simultaneous curve fitting (DS-SCF) analysis, we retrieve electrolyte (Na₂CO₃ and NaF)-correlated OH-stretch spectra of water at the air/water interface. Vibrational mapping of the perturbed interfacial water with the hydration shell spectra (obtained by DS-SCF analysis of Raman spectra) of the corresponding anion discloses that the kosmotropic electrolytes do not form well-defined EDL at the air/water interface. Instead, the interfacial water forms a stronger hydrogen-bond with the surface-expelled anions (CO₃^2- and F^-) and becomes more inhomogeneous than the pristine air/water interface. Together, the results reveal that the perturbation of interfacial water by the kosmotropic electrolyte is a "local phenomenon" confined within the hydration shell of the surface-expelled anion.

"Breaking" and "Making" of Water Structure at the Air/Water-Electrolyte (NaXO3; X = Cl, Br, I) Interface

The prevalence of ions at the aqueous interface has been widely recognized, but their effect on the structure of interfacial water (e.g., hydrogen (H)-bonding) remains enigmatic. Using heterodyne-detected vibrational sum frequency generation (HD-VSFG) and Raman difference spectroscopy with simultaneous curve fitting (DS-SCF) analysis, we show that the ion-induced perturbations of H-bonding at the air/water interface and in the bulk water are strongly correlated. Specifically, the structure-breaking anions such as ClO₃^- decrease the average H-bonding of water at the air/water interface, as it does to the water in its hydration shell in the bulk. The structure-making anion of the same series (IO₃^-) does exactly the reverse. None of the electrolytes (NaXO₃; X = Cl, Br, I) form well-defined electric double layers that significantly increase or reverse the hydrogen-down (H-down) orientation of water at the air/water interface. These results provide a unified picture of specific anion effect at the air/water interface and in the bulk water.

Absolute ion hydration free energy scale and the surface potential of water via quantum simulation

With a goal of determining an absolute free energy scale for ion hydration, quasi-chemical theory and ab initio quantum mechanical simulations are employed to obtain an accurate value for the bulk hydration free energy of the Na⁺ ion. The free energy is partitioned into three parts: 1) the inner-shell or chemical contribution that includes direct interactions of the ion with nearby waters, 2) the packing free energy that is the work to produce a cavity of size λ in water, and 3) the long-range contribution that involves all interactions outside the inner shell. The interfacial potential contribution to the free energy resides in the long-range term. By averaging cation and anion data for that contribution, cumulant terms of all odd orders in the electrostatic potential are removed. The computed total is then the bulk hydration free energy. Comparison with the experimentally derived real hydration free energy produces an effective surface potential of water in the range -0.4 to -0.5 V. The result is consistent with a variety of experiments concerning acid-base chemistry, ion distributions near hydrophobic interfaces, and electric fields near the surface of water droplets.

Restricted intramolecular rotation of fluorescent molecular rotors at the periphery of aqueous microdroplets in oil

Fluorescent molecular rotor dyes, including Cy3, Cy5, and Alexa Fluor 555, dissolved in micron-sized aqueous droplets (microdroplets) in oil were excited, and the fluorescence intensity was recorded as function of time. We observed lengthening of the fluorescence lifetime of these dyes at the water-oil periphery, which extended several microns inward. This behavior shows that intramolecular rotation is restricted at and near the microdroplet interface. Lengthened lifetimes were observed in water microdroplets but not in microdroplets composed of organic solvents. This lifetime change was relatively insensitive to added glycerol up to 60%, suggesting that solution viscosity is not the dominant mechanism. These restricted intramolecular rotations at and near the microdroplet periphery are consistent with the reduced entropy observed in chemical reactions in microdroplets compared to the same reaction conditions in bulk solution and helps us further understand why microdroplet chemistry differs so markedly from bulk-phase chemistry.

Strong Electric Field Observed at the Interface of Aqueous Microdroplets

Chemical reactions in aqueous microdroplets often exhibit unusual kinetic and thermodynamic properties not observed in bulk solution. While an electric field has been implicated at the water interface, there has been no direct measurement in aqueous microdroplets, largely due to the lack of proper measurement tools. Herein, we employ newly developed stimulated Raman excited fluorescence microscopy to measure the electric field at the water-oil interface of microdroplets. As determined by the vibrational Stark effect of a nitrile-bearing fluorescent probe, the strength of the electric field is found to be on the order of 10⁷ V/cm. This strong electric field aligns probe dipoles with respect to the interface. The formation of the electric field likely arises from charge separation caused by the adsorption of negative ions at the water-oil interface of microdroplets. We suggest that this strong electric field might account in part for the unique properties of chemical reactions reported in microdroplets.

Massive generation of metastable bulk nanobubbles in water by external electric fields

Nanobubbles (NBs) are nanoscopic gaseous domains than can exist on solid surfaces or in bulk liquids. They have attracted substantial attention due to their long-time (meta)stability and a high potential for real-world applications. Using an approach not previously investigated, we exploit surface-electrostatic NB formation and stabilization via application of external electric fields in gas-liquid systems, with the marked result of massively increased gas uptake into the liquid in NB form. The de facto gas solubility enhancement (over many months) ranges from 2.5-fold for oxygen to 30-fold for methane vis-à-vis respective Henry's law values for gas solubility; the more hydrophobic the gas, the more spectacular the increase. Molecular dynamics simulations reveal that the origin of NBs' movement lies in dielectrophoresis, while substantial NB stabilization arises from a surface-polarization interaction.

Electromechanics of the liquid water vapour interface

Two collective properties distinguishing the thin liquid water vapour interface from the bulk liquid are the anisotropy of the pressure tensor giving rise to surface tension and the orientational alignment of the molecules leading to a finite dipolar surface potential. Both properties can be regarded as capillary phenomena and are likely to be coupled. We have investigated this coupling by determining the response of the tangential component of the surface tension to the application of an electric field normal to the surface using finite field molecular dynamics simulations. We find an upside down parabola with a maximum shifted away from zero field. Comparing the molecular dynamics results to a phenomenological electromechanical model we relate the zero field derivative of the tangential part of the surface tension to the electrostatic potential generated by the spontaneous interface polarization. When interpreted with this model our simulations also indicate that Kelvin forces due to electric field gradients at a polarized interface play an important role in the effective dielectric response.

Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution

We estimate both single ion hydration Gibbs free energies in water droplets, comprising from 50 to 1000 molecules, and water/vacuum surface potentials in pure water droplets comprising up to 10 000 molecules. We consider four ions, namely, Li⁺, NH₄ ⁺, F^-, and Cl^-, and we model their hydration process and water/water interactions using polarizable force fields based on an induced point dipole approach. We show both ion hydration Gibbs free energies and water surface potentials to obey linear functions of the droplet radius as soon as droplets comprising a few hundred water molecules. Moreover, we also show that the differences in anion/cation hydration Gibbs free energies in droplets obey a different regime in large droplets than in small clusters comprising no more than six water molecules, in line with the earlier results computed from standard additive point charge force fields. Hence, both point charge and more sophisticated induced point dipole molecular modeling approaches suggest that methods considering only the thermodynamical properties of small ion/water clusters to estimate the absolute proton hydration Gibbs free energy in solution are questionable. In particular, taking into account the data of large ion/water droplets may yield a proton hydration Gibbs free energy in solution value to be shifted by several k_BT units compared to small clusters-based approaches.

Second-Order Vibrational Lineshapes from the Air/Water Interface

We explore by means of modeling how absorptive-dispersive mixing between the second- and third-order terms modifies the imaginary χ_total⁽²⁾ responses from air/water interfaces under conditions of varying charge densities and ionic strength. To do so, we use published Im(χ⁽²⁾) and χ⁽³⁾ spectra of the neat air/water interface that were obtained either from computations or experiments. We find that the χ_total⁽²⁾ spectral lineshapes corresponding to experimentally measured spectra contain significant contributions from both interfacial χ⁽²⁾ and bulk χ⁽³⁾ terms at interfacial charge densities equivalent to less than 0.005% of a monolayer of water molecules, especially in the 3100 to 3300 cm^-1 frequency region. Additionally, the role of short-range static dipole potentials is examined under conditions mimicking brine. Our results indicate that surface potentials, if indeed present at the air/water interface, manifest themselves spectroscopically in the tightly bonded H-bond network observable in the 3200 cm^-1 frequency range.

Energy level alignment at semiconductor–water interfaces from atomistic and continuum solvation models

Efficient electronic energy level alignment at solid–liquid interfaces with continuum solvation models.