Spectroscopic studies of solvated hydrogen and hydroxide ions at aqueous surfaces

Measuring the molecular properties of the surface of acidic and basic aqueous solutions is essential to understanding a wide range of important biological, chemical, and environmental processes on our planet. In the present studies, vibrational sum-frequency spectroscopy (VSFS) is employed in combination with isotopic dilution experiments at the vapor/water interface to elucidate the interfacial water structure as the pH is varied with HCl and NaOH. In acidic solutions, solvated proton species are seen throughout the interfacial region, and they alter the hydrogen bonding between water molecules in ways that reflect their depth in the interfacial region. At the higher frequencies of the OH stretch region, there is spectral evidence for solvated proton species residing in the topmost layers of the interfacial region. As reported in previous VSF studies, more strongly bound solvated proton species are observed at lower OH stretching frequencies. The solvated proton species that have stronger hydrogen bonding are similar in structure to those found in bulk acid solutions and likely reside somewhat deeper in the interfacial region. There is also evidence of OH stretching from solvated protons and relatively strong hydrogen bonding in the solvation sphere that is similar to other solvated ions. In contrast, water molecules solvating OH(-) ions show relatively weak hydrogen bonding and significantly less interfacial order. VSF spectra are acquired under multiple polarizations to provide crucial information for the interpretation of the spectra and for the determination of interfacial structure.     

Structure and dynamics of water confined in silica nanopores

We report the results of molecular simulation of water in silica nanopores at full hydration and room temperature. The model systems are approximately cylindrical pores in amorphous silica, with diameters ranging from 20 to 40 ?. The filled pores are prepared using grand canonical Monte Carlo simulation and molecular dynamics simulation is used to calculate the water structure and dynamics. We found that water forms two distinct molecular layers at the interface and exhibits uniform, but somewhat lower than bulk liquid, density in the core region. The hydrogen bond density profile follows similar trends, with lower than bulk density in the core and enhancements at the interface, due to hydrogen bonds between water and surface non-bridging oxygens and OH groups. Our studies of water dynamics included translational mean squared displacements, orientational time correlations, survival probabilities in interfacial shells, and hydrogen bond population relaxation. We found that the radial-axial anisotropy in translational motion largely follows the predictions of a model of free diffusion in a cylinder. However, both translational and rotational water mobilities are strongly dependent on the proximity to the interface, with pronounced slowdown in layers near the interface. Within these layers, the effects of interface curvature are relatively modest, with only a small increase in mobility in going from the 20 to 40 ? diameter pore. Hydrogen bond population relaxation is nearly bulk-like in the core, but considerably slower in the interfacial region.     

Molecular Structure and Dynamics of Water at the Water–Air Interface Studied with Surface‐Specific Vibrational Spectroscopy

Water interfaces provide the platform for many important biological, chemical, and physical processes. The water–air interface is the most common and simple aqueous interface and serves as a model system for water at a hydrophobic surface. Unveiling the microscopic (<1 nm) structure and dynamics of interfacial water at the water–vapor interface is essential for understanding the processes occurring on the water surface. At the water interface the network of very strong intermolecular interactions, hydrogen‐bonds, is interrupted and the density of water is reduced. A central question regarding water at interfaces is the extent to which the structure and dynamics of water molecules are influenced by the interruption of the hydrogen‐bonded network and thus differ from those of bulk water. Herein, we discuss recent advances in the study of interfacial water at the water–air interface using laser‐based surface‐specific vibrational spectroscopy.     

Hydrophilic carbon nanotube membrane enhanced interfacial evaporation for desalination

Carbon nanotube-based(CNT-based) interfacial evaporation material is one of the most potential materials for solar desalination. Here, we studied the evaporation rate of the CNT-based membranes with different hydrophilic and hydrophobic chemical modified surfaces using molecular dynamic simulations.We found that the hydrogen bonding density among water molecules at the interface is a key factor in enhancing the evaporation rate. For a hydrophilic CNT-based membrane, the strong interactions betwe...     

Protein resistance of (ethylene oxide)n monolayers at the air/water interface: effects of packing density and chain length

Protein adsorption on poly(ethylene oxide) (PEO) and oligo(ethylene oxide) (OEO) monolayers is studied at different packing densities using the Langmuir technique. In the case of a PEO monolayer, a protein adsorption minimum is revealed at sigma(-1) = 10 nm(2) for both lysozyme and fibrinogen. Manifested are two packing density regimes of steric repulsion and compressive attraction between PEO and a protein on top of the overall attraction of the protein to the air/water interface. The observed protein adsorption minimum coincides with the maximum of the surface segment density at sigma(-1) = 10 nm(2). However, OEO monolayer presents a different scenario, namely that the amount of protein adsorbed decreases monotonically with increasing packing density, indicating that the OEO chains merely act as a steric barrier to protein adsorption onto the air/water interface. Besides, in the adsorption of fibrinogen, three distinct kinetic regimes controlled by diffusion, penetration and rearrangement are recognized, whereas only the latter two were made out in the adsorption of lysozyme.     

Investigating the quartz (1010)/water interface using classical and ab initio molecular dynamics

Two different terminations of the (1010) surface of quartz (α and β) interacting with water are simulated by classical (CMD) (using two different force fields) and ab initio molecular dynamics (AIMD) and compared with previously published X-ray reflectivity (XR) experiments. Radial distribution functions between hydroxyl and water show good agreement between AIMD and CMD using the ClayFF force field for both terminations. The Lopes et al. (Lopes, P. E. M.; Murashov, V.; Tazi, M.; Demchuk, E.; MacKerell, A. D. J. Phys. Chem. B2006, 110, 2782-2792) force field (LFF), however, underestimates the extent of hydroxyl-water hydrogen bonding. The β termination is found to contain hydroxyl-hydroxyl hydrogen bonds; the quartz surface hydroxyl hydrogens and oxygens that hydrogen bond with each other exhibit greatly reduced hydrogen bonding to water. Conversely, the hydroxyl hydrogen and oxygens that are not hydrogen bonded to other surface hydroxyls but are connected to those that are show a considerable amount of hydrogen bonding to water. The electron density distribution of an annealed surface of quartz (1010) obtained by XR is in qualitative agreement with electron densities calculated by CMD and AIMD. In all simulation methods, the interfacial water peak appears farther from the surface than observed by XR. Agreement among AIMD, LFF, and XR is observed for the relaxation of the near-surface atoms; however, ClayFF shows a larger discrepancy. Overall, results show that for both terminations of (1010), LFF treats the near-surface structure more accurately whereas ClayFF treats the interfacial water structure more accurately. It is shown that the number of hydroxyl and water hydrogen bonds to the bridging Si-O-Si oxygens connecting the surface silica groups to the rest of the crystal is much greater for the α than the β termination. It is suggested that this may play a role in the greater resistance to dissolution of the β termination than that of the α termination.     

Dynamics of water confined within reverse micelles

We report structural and dynamical properties of water confined within reverse micelles (RMs) ranging in size from R = 10 A to R = 23 A as determined from molecular dynamics simulations. The low-frequency infrared spectra have been calculated using linear response theory and depend linearly on the fraction of bulklike water within the RMs. Furthermore, these spectra show nearly isosbestic behavior in the region near 660 cm(-1). Both of these characteristics are present in previously measured experimental spectra. The single dipole spectra for interfacial trapped, bound, and bulklike water within the RMs have also been calculated and show region-dependent frequency shifts. Specifically, the bound and trapped water spectra have a peak at lower frequencies than that for the inner core water. We therefore assign the low-frequency band in the IR spectra to bound and trapped interfacial water. Finally, region-dependent hydrogen bonding profiles and spatial distribution functions are also presented.     

Effect of Na+, Mg2+, and Zn2+ chlorides on the structural and thermodynamic properties of water/n-heptane interfaces

The effect on the structural and thermodynamic properties in water/n-heptane interfaces on addition of NaCl, MgCl(2), and ZnCl(2) has been examined through five independent 100-ns molecular dynamics simulations. Results indicate that the interfacial thickness within the framework of the capillary-wave model decreases on addition of electrolytes in the order Na(+) < Mg(2+) < Zn(2+), whereas the interfacial tension increases in the same order. Ionic density profiles and self-diffusion coefficients are strongly influenced by the strength of the first hydration shell, which varies in the order Na(+) < Mg(2+) < Zn(2+). On the other hand, the Cl(-) behavior, that is, diffusion and solvation sphere, is influenced by its counterion. Accordingly, cations are strongly expelled from the interface, which is especially remarkable for the small divalent cations. This fact alters the water geometry near the interface and in a lesser extent n-heptane order and number of hydrogen bonds per water molecule close to the interface.     

Interfacial tension and interfacial profiles: an equation-of-state approach

A quasi-thermodynamic approach of inhomogeneous systems is used for modeling the fluid-fluid interface. It is based on the recently introduced QCHB (quasi-chemical hydrogen bonding) equation-of-state model of fluids and their mixtures, which is used for the estimation of the Helmholtz free energy density difference, Deltapsi(0), between the system with interface and another system of the same constitution but without interface. Consistent expressions for the interfacial tension and interfacial profiles for various properties are presented. The interfacial tension is proportional to the integral of Deltapsi(0) along the full height of the system, the proportionality constant being equal to 1, when no density gradient contributions are taken into consideration, 2, when the Cahn-Hilliard approximation is adopted, and 4, when the full density gradient contributions are taken into consideration. A satisfactory agreement is obtained between experimental and calculated surface tensions. Extension of the approach to mixtures is examined along with the associated problems for the numerical calculations of the interfacial profiles. A new equation is derived for the chemical potentials in the interfacial region, which facilitates very much the calculation of the composition profiles across the interface.     

Hydrogen bonded structure and dynamics of liquid-vapor interface of water-ammonia mixture: an ab initio molecular dynamics study

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations.     

Nucleation in a simple model for protein solutions with anisotropic interactions

A lattice analog of density functional theory is used to explore the structural and thermodynamic properties of critical nuclei in mixtures of particles with attractive anisotropic interactions. Protein molecules are assumed to occupy the sites on a regular cubic lattice, with effective directional interactions that mimic hydrogen bonding and the solvation forces induced by water. Interaction parameters are chosen to qualitatively reproduce the phase behavior of protein solutions. Our model predicts that critical nuclei of the solidlike phase have nonspherical shapes, and that their specific geometry depends on the nature of the anisotropic interactions. Molecules tend to align in distinctive ways in the core and in the interfacial region of these critical clusters, and the width and structure of the interface are highly affected by the presence of a metastable fluid-fluid critical point. Close to the critical region, the height of the barrier to nucleation is strongly reduced; this effect is enhanced by increasing the anisotropy of the intermolecular interactions. Unlike systems with short-range isotropic interactions, nucleation in our model is initiated by highly ordered clusters in which the order-disorder transition is confined to the interfacial region.     

Self-accumulation of aromatics at the oil-water interface through weak hydrogen bonding

It is well-known that the amphiphilic solutes are surface-active and can accumulate at the oil-water interface. Here, we have investigated the water and a light-oil model interface by using molecular dynamic simulations. It was found that aromatics concentrated in the interfacial region, whereas the other hydrocarbons were uniformly distributed throughout the oil phase. Similar to previous studies, such concentrations were not observed at pure aromatics-water interfaces. We show that the self-accumulation of aromatics at the oil-water interface is driven by differences in the interfacial tension, which is lower for aromatics-water than between the others. The weak hydrogen bonding between the aromatic rings and the water protons provides the mechanism for lowering the interfacial tension.     

表面电荷及外电场对液固界面热阻的影响

采用非平衡分子动力学方法模拟了外电场及固体表面电荷对水与固体间界面热阻的影响. 结果表明,外加电场平行于界面时, 其对界面热阻几乎没有影响, 而垂直于界面时, 界面热阻将随着电场强度的增大而减小. 壁面带正电荷或负电荷都将使得界面热阻减小. 界面热阻与表面电荷密度及电场强度均满足二次函数关系. 模拟结果表明施加外电场和表面电荷是控制液固界面热阻的有效方法.     

Investigation of interfacial and structural properties of CTAB at the oil/water interface using dissipative particle dynamics simulations

We have used dissipative particle dynamics (DPD) to simulate the system of cetyltrimethylammonium bromide (CTAB) monolayer at the oil/water interface. The interfacial properties (interfacial density, interfacial thickness, and interfacial tension), structural properties (area compressibility modulus, end to end distance, and order parameter), and their dependence on the oil/water ratio and the surfactant concentration were investigated. Three different microstructures, spherical oil in water (o/w), interfacial phase, and water in oil (w/o), can be clearly observed with the oil/water ratio increasing. Both the snapshots and the density profiles of the simulation show that a well defined interface exists between the oil and water phases. The interface thickens with CTAB concentration and oil/water ratio. The area compressibility modulus decreases with an increase in the oil/water ratio. The CTAB molecules are more highly packed at the interface and more upright with both concentration and oil/water ratio. The root mean square end-to-end distance and order parameter have a very weak dependence on the oil/water ratio. But both of them show an increase with CTAB concentration, indicating that the surfactant molecules at the interface become more stretched and more ordered at high concentration. As CTAB concentration increases further, the order parameter decreases instead because the bending of the interface. At the same time, it is shown that CTAB has a high interfacial efficiency at the oil/water interface.     

Structure and dynamics of water at the interface with phospholipid bilayers

We have performed two molecular-dynamics simulations to study the structural and dynamical properties of water at the interface with phospholipid bilayers. In one of the simulations the bilayer contained neutral phospholipid molecules, dioleoylphosphatidylcholine (DOPC); in the second simulation the bilayer contained charged lipid molecules, dioleoylphosphatidylserine (DOPS). From the density profile of water we observe that water next to the DOPS bilayer is more perturbed as compared to water near the DOPC bilayer. Using an energetic criterion for the determination of hydrogen bonding we find that water molecules create strong hydrogen bonds with the headgroups of the phospholipid molecules. Due to the presence of these bonds and also due to the confinement of water, the translational and orientational dynamics of water at the interface are slowed down. The degree of slowing down of the dynamics depends upon the location of water molecules near a lipid headgroup.     

Orientational motions of vibrational chromophores in molecules at the air/water interface with time-resolved sum frequency generation

The first time-resolved experiments in which interfacial molecules are pumped to excited electronic states and probed by vibrational sum frequency generation (SFG) are reported. This method was used to measure the out-of-plane rotation dynamics, i.e. time dependent changes in the polar angle, of a vibrational chromophore of an interfacial molecule. The chromophore is the carbonyl group, the rotation observed is that of the -C=O bond axis, with respect to the interfacial normal, and the interfacial molecule is coumarin 314 (C314) at the air/water interface. The orientational relaxation time was found to be 220+/-20 ps, which is much faster than the orientational relaxation time of the permanent dipole moment axis of C314 at the same interface, as obtained from pump-second harmonic probe experiments. Possible effects on the rotation of the -C=O bond axis due to the carbonyl group hydrogen bonding with interfacial water are discussed. From the measured equilibrium orientation of the permanent dipole moment axis and the carbonyl axis, and knowledge of their relative orientation in the molecule, the absolute orientation of C314 at the air/water interface is obtained.     

Molecular dynamics of phenol at the liquid-vapor interface of water

Molecular dynamics results are presented for phenol at the water liquid-vapor interface at 300 K. The calculated excess free energy of phenol at the interface is -2.8 +/- 0.4 kcal/mol, in good agreement with the recent experimental results of Eisenthal and co-workers. The most probable orientation of the phenol molecule at the surface is such that the aromatic ring is perpendicular to the interface and the OH group is fully immersed in water. The hydroxyl substituent has a preferred orientation which is similar to the orientation of OH bonds of water at the pure water liquid-vapor interface. The transition between interfacial and bulk-like behavior of phenol is abrupt and occurs when the center of mass of the solute is located about 6 angstroms from the Gibbs surface of water. In this region the para carbon atom of the hydrophobic benzene ring can reach the interface and become partially dehydrated. This result suggests that the width of the interfacial region in which the behavior of a simple amphiphilic solute in water is influenced by the presence of the surface depends primarily on the size of its hydrophobic part. The role of the OH substituent was investigated by comparing phenol at the interface with two model systems: benzene with and without partial charges on carbon and hydrogen atoms. It is shown that in the absence of the hydrophilic substituent the solute is located further away from the liquid phase and is more likely to be oriented parallel to the interface. However, when the center of mass of the solute is moved into the interfacial region where the density of water approaches that of the bulk solvent, all three molecules become oriented perpendicularly to the surface. In this orientation the work of cavity formation needed to accommodate the hydrophobic ring in aqueous solvent is minimized.     

Femtosecond Hydrogen Bond Dynamics of Bulk‐like and Bound Water at Positively and Negatively Charged Lipid Interfaces Revealed by 2D HD‐VSFG Spectroscopy

Interfacial water in the vicinity of lipids plays an important role in many biological processes, such as drug delivery, ion transportation, and lipid fusion. Hence, molecular‐level elucidation of the properties of water at lipid interfaces is of the utmost importance. We report the two‐dimensional heterodyne‐detected vibrational sum frequency generation (2D HD‐VSFG) study of the OH stretch of HOD at charged lipid interfaces, which shows that the hydrogen bond dynamics of interfacial water differ drastically, depending on the lipids. The data indicate that the spectral diffusion of the OH stretch at a positively charged lipid interface is dominated by the ultrafast (<～100 fs) component, followed by the minor sub‐picosecond slow dynamics, while the dynamics at a negatively charged lipid interface exhibit sub‐picosecond dynamics almost exclusively, implying that fast hydrogen bond fluctuation is prohibited. These results reveal that the ultrafast hydrogen bond dynamics at the positively charged lipid–water interface are attributable to the bulk‐like property of interfacial water, whereas the slow dynamics at the negatively charged lipid interface are due to bound water, which is hydrogen‐bonded to the hydrophilic head group.