Water and ions at interfaces

Solution and Surface Solvation of Nitrate Anions with Iron(III) and Aluminum(III) in Aqueous Environments: A Raman and Vibrational Sum Frequency Generation Study

Hydrated trivalent metal nitrate salts, Fe(NO3)3·9H2O and Al(NO3)3·9H2O, in both solid and aqueous phases are investigated. Raman and surface-selective vibrational sum frequency generation (SFG) spectroscopy, are used to shed light on ion-ion interactions and hydration in several spectral regions spanning low frequency (440-550 cm-1) to higher frequency modes of nitrate and water (720, 1050, 1250-1450, and 2800-3750 cm-1). These frequencies span the metal-water mode, nitrate in-plane deformation, nitrate symmetric and asymmetric modes, and the OH stretch of condensed phase water molecules. Comparison to NaNO3, and in some cases KNO3, is also shown, providing insight. Splitting and frequency shifts are observed and discussed for both the solid state and solution phase. The Lewis acidity of Fe3+ and Al3+ ions plays a significant role in the observed spectra, in particular for the nitrate asymmetric band splitting and frequency shift. The spectral response from water solvation for iron and aluminum nitrates is nonlinear as compared to linear for sodium nitrate, suggesting significantly different solvation environments that are limited by water hydration capacity at higher concentrations. Moreover, a non-hydrogen bonded OH, dangling OH, from hydrating water molecules is observed spectroscopically for Al and Fe nitrate solutions. Furthermore, aluminum nitrate perturbs the surface water structure more than iron nitrate despite aluminum being a weaker Lewis acid. The surface water structure is thus found to be unique for the Al(NO3)3 solutions as compared to both Fe(NO3)3 and NaNO3, such that surface solvation is more pronounced. This observation exemplifies the nature of the Fe(III) and Al(III) ions and their substantial influence on the surface water structure.

Hydrovoltaic Effects from Mechanical-Electric Coupling at the Water-Solid Interface

The natural water cycle on the Earth carries an enormous amount of energy as thirty-five percent of solar energy reaching the Earth's surface goes into water. However, only a very marginal part of the contained energy, mostly kinetic energy of large volume bulk water, is harvested by hydroelectric power plants. Natural processes in the water cycle, such as rainfall, water evaporation, and moisture adsorption, are widespread but have remained underexploited in the past due to the lack of appropriate technologies. In the past decade, the emergence of hydrovoltaic technology has provided ever-increasing opportunities to extend the technical capability for energy harvesting from the water cycle. Featuring electricity generation from mechanical-electric coupling at the water-solid interface, hydrovoltaic technology embraces almost all dynamic processes associated with water, including raining, waving, flowing, evaporating, and moisture adsorbing. This versatility in dealing with various forms of water and associated energy renders hydrovoltaic technology a solution for fossil fuel-caused environmental problems. Here, we review the current progress of hydrovoltaic energy harvesting from water motion, evaporation, and ambient moisture. Device configuration, energy conversion mechanism mediated by mechanical-electric coupling at various water-solid interfaces, as well as materials selection and functionalization are discussed. Useful strategies guided by established mechanisms for device optimization are then covered. Finally, we provide an outlook on this emerging field and outline the challenges of improving output performance toward potential practical applications.

Do soft anions promote protein denaturation through binding interactions? A case study using ribonuclease A

It has long been known that large soft anions like bromide, iodide and thiocyanate are protein denaturing agents, but their mechanism of action is still unclear. In this work we have investigated the protein denaturing properties of these anions using Ribonuclease A (RNase A) as a model protein system. Salt-induced perturbations to the protein folding free energy were determined using differential scanning calorimetry and the results demonstrate that the addition of sodium iodide and sodium thiocyanate significantly decreases the melting temperature of the protein. In order to account for this reduction in protein stability, we show that the introduction of salts that contain soft anions to the aqueous solvent perturbs the protein unfolding free energy through three mechanisms: (a) screening Coulomb interactions that exist between charged protein residues, (b) Hofmeister effects, and (c) specific anion binding to CH and CH2 moieties in the protein polypeptide backbone. Using the micellization of 1,2-hexanediol as a ruler for hydrophobicity, we have devised a practical methodology that separates the Coulomb and Hofmeister contributions of salts to the protein unfolding free energy. This allowing us to isolate the contribution of soft anion binding interactions to the unfolding process. The analysis shows that binding contributions have the largest magnitude, confirming that it is the binding of soft anions to the polypeptide backbone that is the main promoter of protein unfolding.

Effective Polarization in Pairwise Potentials at the Graphene-Electrolyte Interface

At the graphene-electrolyte interface, the polarizability of both the surface and the solution plays a major role in defining the interfacial structure and dynamics of the ions. Current molecular models predict different ion adsorption behavior at the interface depending on whether surface or solution polarization is included in the model. Here, we propose a simple method to parametrize the ion-carbon interaction from density functional theory, implicitly modeling the solution using the conductor-like polarizable continuum model. The new model simultaneously takes into account the polarizability of both the graphene sheet and the solution without the need to use time-consuming polarizable potentials and can predict the ion adsorption trend so far only achievable using first-principles simulations. Simulations performed with 1 M electrolyte solutions of different ions show that cations are strongly adsorbed onto the graphene surface with a trend (Li⁺ < Na⁺ < K⁺) opposite to that predicted by the gas-phase calculations and different from that obtained from the single-ion simulations.

Oppositely Charged Ions at Water-Air and Water-Oil Interfaces: Contrasting the Molecular Picture with Thermodynamics

The surface-active ions tetraphenylarsonium (Ph4As(+)) and tetraphenylboron (Ph4B(-)) have a similar structure but opposite charge. At the solution-air interface, the two ions affect the surface tension in an identical manner, yet sum-frequency generation (SFG) spectra reveal an enhanced surface propensity for Ph4As(+) compared with Ph4B(-), in addition to opposite alignment of interfacial water molecules. At the water-oil interface, the interfacial tension is 7 mN/m lower for Ph4As(+) than for Ph4B(-) salts, but this can be fully accounted for by the different bulk solubility of these ions in the hydrophobic phase, rather than inherently different surface activities. The different solubility can be accounted for by differences in electronic structure, as evidenced by quantum chemical calculations and NMR studies. Our results show that the surface propensity concluded from SFG spectroscopy does not necessarily correlate with interfacial adsorption concluded from thermodynamic measurements.

Liquid–liquid extraction of alkali cations by 18-crown-6: complexation and interface crossing studied by MD and PMF simulations

The 18C6/M⁺Pic⁻complexes form and adsorb “right at the nano-interface” where 18C6 prefers the K⁺guest.

Tuning the isoelectric point of graphene by electrochemical functionalization

The ability to control the charge-potential landscape at solid-liquid interfaces is pivotal to engineer novel devices for applications in sensing, catalysis and energy conversion. The isoelectric point (pI)/point of zero charge (pzc) of graphene plays a key role in a number of physico-chemical phenomena occurring at the graphene-liquid interface. Supported by theory, we present here a methodology to identify the pI/pzc of (functionalized) graphene, which also allows for estimating the nature and extent of ion adsorption. The pI of bare graphene (as-prepared, chemical vapor deposition (CVD)-grown) is found to be less than 3.3, which we can continuously modify up to 7.5 by non-covalent electrochemical attachment of aromatic amino groups, preserving the favorable electronic properties of graphene throughout. Modelling all the observed results with detailed theory, we also show that specific adsorption of ions and the substrate play only an ancillary role in our capability to tune the pI of graphene.

Simple, low-cost styrene-ethylene/butylene-styrene microdevices for electrokinetic applications

Styrene-ethylene/butylene-styrene (SEBS) copolymers combine thermoplastic and elastomeric properties to provide microdevices with the advantageous properties of hard thermoplastics and ease of fabrication similar to PDMS. This work describes the electrical surface properties of SEBS block copolymers using current monitoring experiments to determine zeta potential. We show that SEBS exhibits a stable and relatively high zeta potential magnitude compared to similar polymers. The zeta potential of SEBS is stable when stored in air over time, and no significant differences are observed between different batches and devices, demonstrating reproducibility of results. We show zeta potential trends for varying pH and counterion concentration and demonstrate that SEBS has a repeatable surface potential comparable to glass. Oxygen plasma treatment greatly increases the zeta potential magnitude immediately following treatment before undergoing a moderate hydrophobic recovery to a stable zeta potential. SEBS copolymers also offer simple rapid prototyping fabrication and mass production potential. The presented electrokinetic properties combined with simple, low-cost fabrication of microdevices make SEBS a quality material for electrokinetic research and application development.

The Influence of Dielectric Decrement on Electrokinetics

We treat the dielectric decrement induced by excess ion polarization as a source of ion specificity and explore its impact on electrokinetics. We employ a modified Poisson-Nernst-Planck (PNP) equations accounting for the dielectric decrement. The dielectric decrement is determined by the excess ion polarization parameter α and when α = 0 the standard PNP model is recovered. Our model shows that ions saturate at large zeta potentials (ζ). Because of ion saturation, a condensed counterion layer forms adjacent to the charged surface, introducing a new length scale, the thickness of the condensed layer (l_c ). For the electro-osmotic mobility, the dielectric decrement weakens the electro-osmotic flow owing to the decrease of the dielectric permittivity. At large ζ, when α ≠ 0, the electro-osmotic mobility is found to be proportional to ζ/2, in contrast to ζ predicted by the standard PNP model. This is attributed to ion saturation at large ζ. In terms of the electrophoretic mobility M_e , we carry out both an asymptotic analysis in the thin-double-layer limit and solve the full modified PNP model to compute M_e . Our analysis reveals that the impact of the dielectric decrement is intriguing. At small and moderate ζ, the dielectric decrement decreases M_e with an increasing α. At large ζ, it is well known that the surface conduction becomes significant and plays an important role in determining M_e . It is observed that the dielectric decrement effectively reduces the surface conduction. Hence in stark contrast, M_e increases as α increases. Our predictions of the contrast dependence of the mobility on α at different zeta potentials qualitatively agree with experimental results on the dependence of the mobility among ions and provide a possible explanation for such ion specificity. Finally, the comparisons between the thin-double-layer asymptotic analysis and the full simulations of the modified PNP model suggest that at large ζ the validity of the thin-double-layer approximation is determined by l_c rather than the traditional Debye length.

Self-consistent treatment of the local dielectric permittivity and electrostatic potential in solution for polarizable macromolecular force fields

A self-consistent method is presented for the calculation of the local dielectric permittivity and electrostatic potential generated by a solute of arbitrary shape and charge distribution in a polar and polarizable liquid. The structure and dynamics behavior of the liquid at the solute/liquid interface determine the spatial variations of the density and the dielectric response. Emphasis here is on the treatment of the interface. The method is an extension of conventional methods used in continuum protein electrostatics, and can be used to estimate changes in the static dielectric response of the liquid as it adapts to charge redistribution within the solute. This is most relevant in the context of polarizable force fields, during electron structure optimization in quantum chemical calculations, or upon charge transfer. The method is computationally efficient and well suited for code parallelization, and can be used for on-the-fly calculations of the local permittivity in dynamics simulations of systems with large and heterogeneous charge distributions, such as proteins, nucleic acids, and polyelectrolytes. Numerical calculation of the system free energy is discussed for the general case of a liquid with field-dependent dielectric response.

Structural and technical details of the Kirkwood-Buff integrals from the optimization of ionic force fields: focus on fluorides

Results on the structural details of Kirkwood-Buff integrals obtained from the optimization of ionic force fields are presented. We have proposed and make use of an optimization scheme for ionic force fields, which is based on the modification of the cation-anion mixing rules, the calculation of the thermodynamics properties of various monovalent salt solutions according to the Kirkwood-Buff theory of solutions and the comparison to relevant experimental findings. Here, we complete and extend our calculations and analysis as we focus on the technical details of this optimization procedure and the case of fluorides, which have been proven difficult to handle. Important insight is given on the dependence of the radial distribution functions, the short-ranged potentials of mean force, and the Kirkwood-Buff integrals of the salt solutions on the different scaling factors in the mixing rules. Specifically, the way the structural details and inherent characteristics of the above properties are affected by the quantitative and qualitative differences in the mixing rules for a variety of common biologically relevant monovalent salts is mainly addressed. We conclude on the efficiency of this scheme, again with a focus on the fluorides. In the end, we provide a variation of the ion-pair mixing rules scaling factors with salt concentration to identify regimes for which different mixing rules prefactors lead to well-optimized force fields. All results are obtained through Molecular Dynamics simulations using previously optimized force fields for the monovalent ions.

Polarizable Poisson-Boltzmann equation: the study of polarizability effects on the structure of a double layer

We incorporate ion polarizabilities into the Poisson-Boltzmann equation by modifying the effective dielectric constant and the Boltzmann distribution of ions. The extent of the polarizability effects is controlled by two parameters, γ(1) and γ(2); γ(1) determines the polarization effects in a dilute system and γ(2) regulates the dependence of the polarizability effects on the concentration of ions. For a polarizable ion in an aqueous solution γ(1) ≈ 0.01 and the polarizability effects are negligible. The conditions where γ(1) and/or γ(2) are large and the polarizability is relevant involve the low dielectric constant media, high surface charge, and/or large ionic concentrations.

Effects of surface charge and Gibbs surface energy on the settlement behaviour of barnacle cyprids (Balanus amphitrite)

Gibbs surface energy has long been considered to be an important parameter in the design of fouling-resistant surfaces for marine applications. Rigorous testing of the hypothesis that settlement is related to Gibbs surface energy however has never been accomplished, due mainly to practical limitations imposed by the necessary combination of surface engineering and biological evaluation methods. In this article, the effects of surface charge and Gibbs surface energy on the settlement of cyprids of an important fouling barnacle, Balanus amphitrite, were evaluated. Settlement assays were conducted on a range of self-assembled monolayers (SAMs) (CH(3)-, OH-, COOH-, N(CH(3))(3) (+)-, NH(2)-terminated), presented in gold-coated polystyrene well plates, varying in terms of their surface charge and Gibbs surface energy. Contrary to contemporary theory, settlement was not increased by high-energy surfaces, rather the opposite was found to be the case with cyprids settling in greater numbers on a low-energy CH(3)- SAM compared to a high-energy OH- SAM. Settlement was also greater on negatively-charged SAMs, compared to neutral and positively-charged SAMs. These findings are discussed in the context of data drawn from surfaces that varied in multiple characteristics simultaneously, as have been used previously for such experiments. The finding that surface charge, rather than total surface energy, may be responsible for surface selection by cyprids, will have significant implications for the design of future fouling-resistant materials.

Possible origin of the inverse and direct Hofmeister series for lysozyme at low and high salt concentrations

Protein solubility studies below the isoelectric point exhibit a direct Hofmeister series at high salt concentrations and an inverse Hofmeister series at low salt concentrations. The efficiencies of different anions measured by salt concentrations needed to effect precipitation at fixed cations are the usual Hofmeister series (Cl(-) > NO(3)(-) > Br(-) > ClO(4)(-) > I(-) > SCN(-)). The sequence is reversed at low concentrations. This has been known for over a century. Reversal of the Hofmeister series is not peculiar to proteins. Its origin poses a key test for any theoretical model. Such specific ion effects in the cloud points of lysozyme suspensions have recently been revisited. Here, a model for lysozymes is considered that takes into account forces acting on ions that are missing from classical theory. It is shown that both direct and reverse Hofmeister effects can be predicted quantitatively. The attractive/repulsive force between two protein molecules was calculated. To do this, a modification of Poisson-Boltzmann theory is used that accounts for the effects of ion polarizabilities and ion sizes obtained from ab initio calculations. At low salt concentrations, the adsorption of the more polarizable anions is enhanced by ion-surface dispersion interactions. The increased adsorption screens the protein surface charge, thus reducing the surface forces to give an inverse Hofmeister series. At high concentrations, enhanced adsorption of the more polarizable counterions (anions) leads to an effective reversal in surface charge. Consequently, an increase in co-ion (cations) adsorption occurs, resulting in an increase in surface forces. It will be demonstrated that among the different contributions determining the predicted specific ion effect the entropic term due to anions is the main responsible for the Hofmeister sequence at low salt concentrations. Conversely, the entropic term due to cations determines the Hofmeister sequence at high salt concentrations. This behavior is a remarkable example of the charge-reversal phenomenon.

Development of molecular simulation methods to accurately represent protein-surface interactions: The effect of pressure and its determination for a system with constrained atoms

When performing molecular dynamics simulations for a system with constrained (fixed) atoms, traditional isobaric algorithms (e.g., NPT simulation) often cannot be used. In addition, the calculation of the internal pressure of a system with fixed atoms may be highly inaccurate due to the nonphysical nature of the atomic constraints and difficulties in accurately defining the volume occupied by the unconstrained atoms in the system. The inability to properly set and control pressure can result in substantial problems for the accurate simulation of condensed-phase systems if the behavior of the system (e.g., peptide/protein adsorption) is sensitive to pressure. To address this issue, the authors have developed an approach to accurately determine the internal pressure for a system with constrained atoms. As the first step in this method, a periodically extendable portion of the mobile phase of the constrained system (e.g., the solvent atoms) is used to create a separate unconstrained system for which the pressure can be accurately calculated. This model system is then used to create a pressure calibration plot for an intensive local effective virial parameter for a small volume cross section or "slab" of the system. Using this calibration plot, the pressure of the constrained system can then be determined by calculating the virial parameter for a similarly sized slab of mobile atoms. In this article, the authors present the development of this method and demonstrate its application using the CHARMM molecular simulation program to characterize the adsorption behavior of a peptide in explicit water on a hydrophobic surface whose lattice spacing is maintained with atomic constraints. The free energy of adsorption for this system is shown to be dramatically influenced by pressure, thus emphasizing the importance of properly maintaining the pressure of the system for the accurate simulation of protein-surface interactions.

Charge reversal of surfaces in divalent electrolytes: the role of ionic dispersion interactions

Surface potentials of alkali earth nitrates at a mica surface are calculated using a modified Poisson-Boltzmann approach that includes nonelectrostatic ion-surface dispersion interactions. New ab initio dynamic polarizabilities are used to determine dispersion interactions. A hydration model describing the hydration shell of cations is presented. Excellent agreement with experiment is achieved, including charge reversal at high electrolyte concentration without the need for site binding models. This suggests that specific ionic dispersion forces provide the mechanism for ion surface binding. An asymptotic surface potential is found in the limit of very high concentration. A Hofmeister series is predicted according to the strength of charge reversal, with Mg > Ca > Sr > Ba. The ion-surface dispersion adsorption energies of hydrated ions appear to explain the apparent repulsive secondary hydration forces observed experimentally between mica surfaces when taken with a surface hydration layer.

Simulations of water at the interface with hydrophilic self-assembled monolayers

Simulations of water at hydrophilic self-assembled monolayer (SAM) surfaces are especially relevant for biological interfaces. Well-defined, atomically smooth surfaces that can be continuously varied are possible with SAMs. These characteristics enable more accurate measurements than many other surfaces with the added advantage of tailoring the surface to treat specific chemical groups. A fundamental question is how solid surfaces affect the structure and dynamics of water. Measurements of the structure and dynamics of water at solid surfaces have improved significantly, but there remain differences among the experiments. In this article, the authors review simulations of water at the interface with hydrophilic SAMs. These simulations find that while the interfacial water molecules are slower than the bulk water molecules, the interfacial dynamics remains that of a liquid. A major biological application of SAMs is for making coatings resistant to protein adsorption. SAMs terminated with ethylene glycol monomers have proven to be excellent at resisting protein adsorption. Understanding the mechanisms behind this resistance remains an unresolved issue. Recent simulations suggest a new perspective of the role of interfacial water and the inseparable interplay between the SAM and the water.

Why direct or reversed Hofmeister series? Interplay of hydration, non-electrostatic potentials, and ion size

A modified Poisson-Boltzmann analysis is made of the double layer interaction between two silica surfaces and two alumina surfaces in chloride electrolyte. The analysis incorporates nonelectrostatic ion-surface dispersion interactions based on ab initio ionic excess polarizabilities with finite ion sizes. A hydration model for the tightly held hydration shell of kosmotropic ions is introduced. A direct Hofmeister series (K > Na > Li) is found at the silica surface while the reversed series (Li > Na > K) is found at alumina, bringing theory in line with experiment for the first time. Calculations with unhydrated ions also suggest that surface-induced dehydration may be occurring at the alumina surface.

Ion specific surface charge density of SBA-15 mesoporous silica

Potentiometric titrations were used to estimate the surface charge density of SBA-15 mesoporous silica in different salt solutions. It was found that surface charge depends both on cation type, following a Hofmeister series (Cs(+) < Guanidinium(+) < K(+) < Na(+) < Li(+)), and on salt concentration (in the range 0.05-1 M). The surface charge series is reproduced by theoretical calculations performed using a modified Poisson-Boltzmann equation that includes ionic dispersion forces with ab initio ion polarizabilities and hydrated ions. The hydration model assigns an explicit hydration shell to kosmotropic (strong hydrated) ions only. The Hofmeister series appears to be due to the combination of ion-surface dispersion interactions and ion hydration.

Importance of accurate dynamic polarizabilities for the ionic dispersion interactions of alkali halides

Ab initio quantum mechanical calculations of the dynamic polarizability of alkali metal and halide ions are performed as a function of imaginary frequency. Electron correlation is shown to provide a significant correction to ionic polarizabilities. Ab initio ion-surface dispersion coefficients are compared with single- and multimode London approximations. The commonly employed single-mode model with the characteristic frequency taken from the ionization potential of the ion is shown to be inadequate, underestimating dispersion forces with an average error around 40% or as high as 80% for halide ions. Decomposition of the polarizability data into five modes covers the major modes of each ion adequately (four modes for Li(+)). Illustrative calculations of surface potentials at the mica surface in aqueous alkali halide electrolytes are made. Charge reversal is obtained with the more polarizable cations, K(+) and Rb(+). The error in the single-mode ionization potential models is seen as a strong shift in the surface potential from negative toward positive values.

Analysis of the ion distribution at a charged solid-liquid interface using X-ray standing waves

Functionalized solid-liquid interfaces were analyzed by X-ray standing waves (XSW) combined with streaming current measurements to study surface charges, interfacial potential, and ion distributions. Thin films of aqueous solution containing Br(-) anions and Fe(3+) cations at a concentration of 10 mg/L were prepared on functionalized silicon wafers. Functionalization of Si surfaces was accomplished by aminosilane groups shifting the interfacial potential toward positive values. The ion distribution was measured with nanometer resolution, which allows distinguishing between absorbed and mobile ions at the surface and in the diffusive layer, respectively. For Br(-), different degrees of ion attraction were measured for the pH values 5.7 and 2.8. The ion Debye length values of the diffuse layer were 4 and 2 nm, respectively.

Interfacial thermodynamics of confined water near molecularly rough surfaces

We study the effects of nanoscopic roughness on the interfacial free energy of water confined between solid surfaces. SPC/E water is simulated in confinement between two infinite planar surfaces that differ in their physical topology: one is smooth and the other one is physically rough on a sub-nanometre length scale. The two thermodynamic ensembles considered, with constant pressure either normal or parallel to the walls, correspond to different experimental conditions. We find that molecular-scale surface roughness significantly increases the solid-liquid interfacial free energy compared to the smooth surface. For our surfaces with a water-wall interaction energy minimum of -1.2 kcal mol(-1), we observe a transition from a hydrophilic surface to a hydrophobic surface at a roughness amplitude of about 3 angstroms and a wavelength of 11.6 angstroms, with the interfacial free energy changing sign from negative to positive. In agreement with previous studies of water near hydrophobic surfaces, we find an increase in the isothermal compressibility of water with increasing surface roughness. Interestingly, average measures of the water density and hydrogen-bond number do not contain distinct signatures of increased hydrophobicity. In contrast, a local analysis indicates transient dewetting of water in the valleys of the rough surface, together with a significant loss of hydrogen bonds, and a change in the dipole orientation toward the surface. These microscopic changes in the density, hydrogen bonding, and water orientation contribute to the large increase in the interfacial free energy, and the change from a hydrophilic to a hydrophobic character of the surface.

Zeta potential and electroosmotic mobility in microfluidic devices fabricated from hydrophobic polymers: 2. Slip and interfacial water structure

We discuss the structure of water at hydrophobic interfaces from the standpoint of its impact on electrokinetic phenomena in microfluidic devices fabricated from hydrophobic polymers such as Teflon or Zeonor. Water structuring at hydrophobic interfaces has been described as a source of interfacial charge (see Part 1, this issue), and dewetting phenomena, whether via depletion layers or nanobubbles, contribute to slip and enhanced apparent electrokinetic potentials. Issues concerning the impact of hydrodynamic slip and the role of diffuse interfacial structures are discussed. These issues are coupled with each other and with interfacial charge concerns, providing challenges for measurements of individual parameters.

Bubbles, gating, and anesthetics in ion channels

We suggest that bubbles are the bistable hydrophobic gates responsible for the on-off transitions of single channel currents. In this view, many types of channels gate by the same physical mechanism-dewetting by capillary evaporation-but different types of channels use different sensors to modulate hydrophobic properties of the channel wall and thereby trigger and control bubbles and gating. Spontaneous emptying of channels has been seen in many simulations. Because of the physics involved, such phase transitions are inherently sensitive, unstable threshold phenomena that are difficult to simulate reproducibly and thus convincingly. We present a thermodynamic analysis of a bubble gate using morphometric density functional theory of classical (not quantum) mechanics. Thermodynamic analysis of phase transitions is generally more reproducible and less sensitive to details than simulations. Anesthetic actions of inert gases-and their interactions with hydrostatic pressure (e.g., nitrogen narcosis)-can be easily understood by actions on bubbles. A general theory of gas anesthesia may involve bubbles in channels. Only experiments can show whether, or when, or which channels actually use bubbles as hydrophobic gates: direct observation of bubbles in channels is needed. Existing experiments show thin gas layers on hydrophobic surfaces in water and suggest that bubbles nearly exist in bulk water.

Structure and dynamics of water near the interface with oligo(ethylene oxide) self-assembled monolayers

We performed molecular dynamics simulations of the oligo(ethylene oxide) (OEO) self-assembled monolayers in water to determine the nature of the systems' interfacial structure and dynamics. The density profiles, hydrogen bonding, and water dynamics are calculated as a function of the area per molecule A of OEO. At the highest coverages, the interface is hydrophobic, and a density drop is found at the interface. The interfacial region becomes more like bulk water as A increases. The OEO and water become progressively more mixed, and hydrogen bonding increases within the interfacial region. Water mobility is slower within the interfacial region, but not substantially. The implications of our results on the resistance of OEO SAMs to protein adsorption are discussed. Our principal result is that as A increases the increasingly waterlike interfacial region provides a more protein-resistant surface. This finding supports recent experimental measurements that protein resistance is maximal for less than full coverage on Au.

The isoelectric point/point-of zero-charge of interfaces formed by aqueous solutions and nonpolar solids, liquids, and gases

From our previous work on the role of the electrostatic field strength in controlling the pH of the iso-electric point (iep)/point-of-zero-charge (pzc) of polar solids we have extended the analysis to predict that the pH of the iep/pzc of a nonpolar solid, liquid or gas-aqueous interface should occur at pH 1.0-3.0, dependent on the value assigned to water molecules or clusters at the interface. Consideration of a wide range of experimental results covering nonpolar solids such as molybdenite, stibnite, paraffin, etc. as well as hydrocarbon liquids such as xylene, decalin, and long chain (>C8) alkane oils, as well as nitrogen and hydrogen gases, all in various simple 1:1 electrolyte solutions confirm the general validity of the result. We further consider various models of the origin of the charge on nonpolar material-water interfaces.

Ion-specific anomalous electrokinetic effects in hydrophobic nanochannels

We show with computer simulations that anomalous electrokinetic effects, such as ion specificity and nonzero zeta potentials for uncharged surfaces, are generic features of electro-osmotic flow in hydrophobic channels. This behavior is due to the stronger attraction of larger ions to the "vapor-liquidlike" interface induced by a hydrophobic surface. We propose an analytical model involving a modified Poisson-Boltzmann description for the ion density distributions that quantitatively predicts the anomalous flow profiles. This description includes as a crucial component an ion-size-dependent hydrophobic solvation energy. These results provide an effective framework for predicting ion-specific effects, with potentially important implications for biological modeling.