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.     

The effect of the counterion on water mobility in reverse micelles studied by molecular dynamics simulations

In this study, mobility and structure of water molecules in Aerosol OT (bis(2-ethylhexyl) sulfosuccinate, AOT) reverse micelles with water content w0 = 5 and Na+, K+, Cs+ counterions have been explored with molecular dynamics (MD) simulations. Using the Faeder/Ladanyi model (J. Phys. Chem. B, 2000, 104, 1033) of the reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function, FS(Q,t), for water hydrogen atoms that could be measured in a quasielastic neutron scattering experiment. Separate intermediate scattering functions FRS(Q,t) and FCMS(Q,t) were determined for rotational and translational motion. We find that the decay of FCMS(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior arises from decreased water mobility for molecules close to the interface and from confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay, which is consistent with relatively rapid restricted rotation and slower rotational relaxation over the full angular range. Rotational relaxation is anisotropic, with the O-H bond short-time rotational mobility considerably higher than that of the molecular dipole. This behavior is related to the decreased density of water-water hydrogen bonds in the vicinity of the interface compared to core or bulk water. We find that the interfacial mobility of water molecules is quite different for the three counterion types, but that the core mobility exhibits weak counterion dependence. Differences in interfacial mobility are strongly correlated with structural features, especially ion-water coordination, and the extent of disruption by the counterions of the water hydrogen bond network.     

A first principles simulation study of fluctuations of hydrogen bonds and vibrational frequencies of water at liquid-vapor interface

We present a theoretical study of the structure and dynamics of water-vapor interface by means of ab initio molecular dynamics simulations. The inhomogeneous density, hydrogen bond and orientational profiles, voids and vibrational frequency distributions are investigated. We have also studied various dynamical properties of the interface such as diffusion, orientational relaxation, hydrogen bond dynamics and vibrational frequency fluctuations. The diffusion and orientational relaxation of water molecules are found to be faster at the interface which can be correlated with the voids present in the system. The hydrogen bond dynamics, however, is found to be slightly slower at the interface than that in bulk water. The correlations of hydrogen bond relaxation with the dynamics of vibrational frequency fluctuations are also discussed.     

Hydrogen bond dynamics and microscopic structure of confined water inside carbon nanotubes

We have investigated the density and temperature dependences of microscopic structure and hydrogen bond dynamics of water inside carbon nanotubes (CNTs) using molecular dynamics simulation. The CNTs are treated as rigid, and smoothly truncated extended simple point charge water model is adopted. The results show that as the overall density increases, the atomic density profiles of water inside CNTs become sharper, the peaks shift closer to the wall, and a new peak of hydrogen atomic density appears between the first (outermost) and second layer. The intermittent hydrogen bond correlation function C(HB)(t) of water inside CNTs decays slower than that of bulk water, and the rate of decay decreases as the tube diameter decreases. C(HB)(t) clearly decays more slowly for the first layer of water than for other regions inside CNTs. The C(HB)(t) of the interlayer hydrogen bonds decays faster than those of the other regions and even faster than that of the bulk water. On the other hand, the hydrogen bond lifetimes of the first layer are shorter than those of the inner layer(s). Interlayer hydrogen bond lifetimes are clearly shorter than those of the constituent layers. As a whole, the hydrogen bond lifetimes of water inside CNTs are shorter than those of bulk water, while the relaxation of C(HB)(t) is slower for the confined water than for bulk water. In other words, hydrogen bonds of water inside CNTs break more easily than those of bulk water, but the water molecules remain in each other's vicinity and can easily reform the bonds.     

Reorientational Dynamics of Water Confined in Zeolites

We present a detailed molecular‐dynamics study of water reorientation and hydrogen‐bond dynamics in a strong confinement situation, within the narrow pores of an all‐silica Linde type A (LTA) zeolite. Two water loadings of the zeolite are compared with the bulk case. Water dynamics are retarded in this extreme hydrophobic confinement and the slowdown is more pronounced at higher water loading. We show that water reorientation proceeds mainly by large‐amplitude angular jumps, whose mechanism is similar to that determined in the bulk. The slowdown upon hydrophobic confinement arises predominantly from an excluded‐volume effect on the large fraction of water molecules lying at the interface with the zeolite matrix, with an additional minor contribution coming from a structuring effect induced by the confinement.     

Aqueous Heterogeneity at the Air/Water Interface Revealed by 2D‐HD‐SFG Spectroscopy

Water molecules interact strongly with each other through hydrogen bonds. This efficient intermolecular coupling causes strong delocalization of molecular vibrations in bulk water. We study intermolecular coupling at the air/water interface and find intermolecular coupling 1) to be significantly reduced and 2) to vary strongly for different water molecules at the interface—whereas in bulk water the coupling is homogeneous. For strongly hydrogen‐bonded OH groups, coupling is roughly half of that of bulk water, due to the lower density in the near‐surface region. For weakly hydrogen‐bonded OH groups that absorb around 3500 cm^?1, which are assigned to the outermost, yet hydrogen‐bonded OH groups pointing towards the liquid, coupling is further reduced by an additional factor of 2. Remarkably, despite the reduced structural constraints imposed by the interfacial hydrogen‐bond environment, the structural relaxation is slow and the intermolecular coupling of these water molecules is weak.     

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.     

Distribution, structure, and dynamics of cesium and iodide ions at the H2O-CCl4 and H2O-vapor interfaces

Molecular dynamics simulations utilizing many-body potentials of H2O-CCl4 and H2O-vapor interfaces were carried out at different cesium and iodide ion concentrations to compare ion distribution, interfacial orientational and structural properties, and dynamics. It was found that cesium was repelled by both interfaces, and iodide was active at both interfaces, but to a much greater degree at the H2O-vapor interface. At the interface, the iodide dipole was strongly induced, orienting perpendicular to the interface for both systems, leading to stronger hydrogen bonds with water. For the H2O-CCl4 interface, though, there was a compensation between these strong hydrogen bonds and short to moderate ranged repulsion between iodide and CCl4. Hydrogen bond distance and angular distributions showed weaker water-water hydrogen bonds at both interfaces, but generally stronger water-iodide hydrogen bonds. Both translational and rotational dynamics of water were faster at the interface, while for CCl4, its translational dynamics was slower, but rotational dynamics faster at the interface. For many of the studied systems and species, translational diffusion was found to be anisotropic at both interfacial and bulk regions.     

Hydrogen-bond dynamics in the air-water interface

Hydrogen-bond (H-bond) dynamics in the air-water interface is studied by molecular dynamics simulations. The analysis reveals that the dynamics of breaking and forming hydrogen bonds in the air-water interface is faster than that in bulk water for the polarizable water models. This is in contrast to the results found on a protein surface. We show that the difference stems from more rapid translational diffusion in the interface. When the effect of pair diffusion is eliminated, the hydrogen-bond dynamics in the interface is observed to be slower than that in the bulk. This occurs because the number of water molecules adjacent to a hydrogen-bonded pair and available to accept or donate a hydrogen bond is smaller in the interface than in the bulk. The comparison between polarizable water models and fixed-charge models highlights the potential importance of the polarization effect in the water-vapor interface.     

Molecular dynamics study of the liquid-vapor interface of acetonitrile: equilibrium and dynamical properties

The equilibrium and dynamical properties of the liquid-vapor interface of pure acetonitrile are studied by means of molecular dynamics simulations. Both nonpolarizable and polarizable models are employed in the present study. For the nonpolarizable model, the simulations are carried out for two different system sizes and at two different temperatures whereas the simulation with the polarizable model is done for a single system. The inhomogeneous density, anisotropic orientational profile, the width of the interface, and also the surface tension are calculated at room temperature and also at a lower temperature of 273 K. The dynamical aspects of the interface are investigated in terms of the single-particle dynamical properties such as the relaxation of velocity autocorrelation and the translational diffusion coefficients along the perpendicular and parallel directions and the dipole orientational relaxation of the interfacial acetonitrile molecules. The results of the interfacial dynamics are compared with those of the corresponding bulk phases at both temperatures. The convergence of the calculated results with respect to the length of simulation runs and the system size are also discussed.     

Anisotropic dynamics of dipolar liquids in narrow slit pores

We report molecular dynamics simulation results for Stockmayer fluids confined to narrow slitlike pores with structureless, nonconducting walls. The translational and rotational dynamics of the dipolar particles have been investigated by calculating autocorrelation functions, diffusion coefficients, and relaxation times for various pore widths (five or less particle diameters) and directions parallel and perpendicular to the walls. The dynamic properties of the confined systems are compared to bulk properties, where corresponding bulk and pore states at the same temperature and chemical potential are determined in parallel grand canonical Monte Carlo simulations. We find that the dynamic behavior inside the pore depends on the distance from the walls and can be strongly anisotropic even in globally isotropic systems. This concerns especially the particles in the surface layers close to the walls, where the single particle and collective dipolar relaxation resemble that of true two-dimensional dipolar fluids with different in-plane and out-of-plane relaxations. On the other hand, bulklike relaxation is observed in the pore center of sufficiently wide pores.     

Hydrogen bond lifetime dynamics at the interface of a surfactant monolayer

The dynamics of water near the polar headgroups of surfactants in a monolayer adsorbed at the air/water interface is likely to play a decisive role in determining the physical behavior of such organized assemblies. We have carried out an atomistic molecular dynamics (MD) simulation of a monolayer of the anionic surfactant sodium bis(2-ethyl-1-hexyl) sulfosuccinate (aerosol-OT or AOT) adsorbed at the air/water interface. The simulation is performed at room temperature with a surface coverage of that at the critical micelle concentration (78 Angstrom(2)/molecule). Detailed analyses of the lifetime dynamics of surfactant-water (SW) and water-water (WW) hydrogen bonds at the interface have been carried out. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the bound states at the interface and quantification of the dynamic equilibrium between bound and quasi-free water molecules, in terms of time-dependent relaxation rates. It is observed that the water molecules present in the first hydration layer form strong hydrogen bonds with the surfactant headgroups and hence have longer lifetimes. Importantly, it is found that the overall relaxation of the SW hydrogen bonds is faster for those water molecules which form two hydrogen bonds with the surfactant headgroups than those forming one such hydrogen bond. Equally interestingly, it is further noticed that water molecules beyond the first hydration layer form weaker hydrogen bonds than pure bulk water.     

Structure and dynamics of the aqueous liquid-vapor interface: a comprehensive particle-based simulation study

This research addresses a comprehensive particle-based simulation study of the structural, dynamic, and electronic properties of the liquid-vapor interface of water utilizing both ab initio (based on density functional theory) and empirical (fixed charge and polarizable) models. Numerous properties such as interfacial width, hydrogen bond populations, dipole moments, and correlation times will be characterized with identical schemes to draw useful conclusions on the strengths and weakness of the proposed models for interfacial water. Our findings indicate that all models considered in this study yield similar results for the radial distribution functions, hydrogen bond populations, and orientational relaxation times. Significant differences in the models appear when examining both the dipole moments and surface relaxation near the aqueous liquid-vapor interface. Here, the ab initio interaction potential predicts a significant decrease in the molecular dipole moment and expansion in the oxygen-oxygen distance as one approaches the interface in accordance with recent experiments. All classical polarizable interaction potentials show a less dramatic drop in the molecular dipole moment, and all empirical interaction potentials studied yield an oxygen-oxygen contraction as the interface is approached.     

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.     

Water structure and dynamics at a silica surface: Pake doublets in 1H NMR spectra

Detailed knowledge about the dynamics and structure of liquids in the vicinity of a solid surface is important in several fields of research. In this study a homogeneous model system of colloidal and nonporous silica particles with a narrow particle size distribution was used to examine such properties of adsorbed water and 1-heptanol. Doublet (1)H water resonances ("Pake doublets") indicate a preferred spatial orientation for the water molecules, as well as a lower molecular density in the surface-induced water structures compared to bulk water. These surface-induced structures are found to extend at least 8 nm from the silica surface. T(1) relaxation measurements at several temperatures indicate weaker H-bonding in the adsorbed water compared to bulk water. T(2) relaxation measurements at several temperatures reveal the presence of two water phases and give quantitative information on the mobility of water molecules and proton exchange processes. The presence of 1-heptanol changes the water characteristics, primarily in the water phase closer to the surface, where water molecules experience decreased translational and increased rotational freedom. In the absence of water, adsorbed 1-heptanol forms surface aggregates encompassing several molecular layers, where the first adsorbed layer shows severe restrictions in mobility and subsequent layers are more mobile.     

Hydrogen bond dynamics and water structure in glucose-water solutions by depolarized Rayleigh scattering and low-frequency Raman spectroscopy

The effect of glucose on the relaxation process of water at picosecond time scales has been investigated by depolarized Rayleigh scattering (DRS) experiments. The process is assigned to the fast hydrogen bonding dynamics of the water network. In DRS spectra this contribution can be safely separated from the slower relaxation process due to the sugar. The detected relaxation time is studied at different glucose concentrations and modeled considering bulk and hydrating water contributions. As a result, it is found that in diluted conditions the hydrogen bond lifetime of proximal water molecules becomes about three times slower than that of the bulk. The effect of the sugar on the hydrogen bond water structure is investigated by analyzing the low-frequency Raman (LFR) spectrum sensitive to intermolecular modes. The addition of glucose strongly reduces the intensity of the band at 170 cm(-1) assigned to a collective stretching mode of water molecules arranged in cooperative tetrahedral domains. These findings indicate that proximal water molecules partially lose the tetrahedral ordering typical of the bulk leading to the formation of high density environments around the sugar. Thus the glucose imposes a new local order among water molecules localized in its hydration shell in which the hydrogen bond breaking dynamics is sensitively retarded. This work provides new experimental evidences that support recent molecular dynamics simulation and thermodynamics results.     

Relaxation and jump dynamics of water at the mica interface

The orientational relaxation dynamics of water confined between mica surfaces is investigated using molecular dynamics simulations. The study illustrates the wide heterogeneity that exists in the dynamics of water adjacent to a strongly hydrophilic surface such as mica. Analysis of the survival probabilities in different layers is carried out by normalizing the corresponding relaxation times with bulk water layers of similar thickness. A 10-fold increase in the survival times is observed for water directly in contact with the mica surface and a non-monotonic variation in the survival times is observed moving away from the mica surface to the bulk-like interior. The orientational relaxation time is highest for water in the contact layer, decreasing monotonically away from the surface. In all cases the ratio of the relaxation times of the 1st and 2nd rank Legendre polynomials of the HH bond vector is found to lie between 1.5 and 1.9 indicating that the reorientational relaxation in the different water layers is governed by jump dynamics. The orientational dynamics of water in the contact layer is particularly novel and is found to undergo distinct two-dimensional hydrogen bond jump reorientational dynamics with an average waiting time of 4.97 ps. The waiting time distribution is found to possess a long tail extending beyond 15 ps. Unlike previously observed jump dynamics in bulk water and other surfaces, jump events in the mica contact layer occur between hydrogen bonds formed by the water molecule and acceptor oxygens on the mica surface. Despite slowing down of the water orientational relaxation near the surface, life-times of water in the hydration shell of the K(+) ion are comparable to that observed in bulk salt solutions.     

Reorientation dynamics of nanoconfined water: power-law decay, hydrogen-bond jumps, and test of a two-state model

The reorientation dynamics of water confined within nanoscale, hydrophilic silica pores are investigated using molecular dynamics simulations. The effect of surface hydrogen-bonding and electrostatic interactions are examined by comparing with both a silica pore with no charges (representing hydrophobic confinement) and bulk water. The OH reorientation in water is found to slow significantly in hydrophilic confinement compared to bulk water, and is well-described by a power-law decay extending beyond one nanosecond. In contrast, the dynamics of water in the hydrophobic pore are more modestly affected. A two-state model, commonly used to interpret confined liquid properties, is tested by analysis of the position-dependence of the water dynamics. While the two-state model provides a good fit of the orientational decay, our molecular-level analysis evidences that it relies on an over-simplified picture of water dynamics. In contrast with the two-state model assumptions, the interface dynamics is markedly heterogeneous, especially in the hydrophilic pore and there is no single interfacial state with a common dynamics.     

Ultrafast energy transfer in water-AOT reverse micelles

A spectroscopic investigation of the vibrational dynamics of water in a geometrically confined environment is presented. Reverse micelles of the ternary microemulsion H2O/AOT/n-octane (AOT = bis-2-ethylhexyl sulfosuccinate or aerosol-OT) with diameters ranging from 1 to 10 nm are used as a model system for nanoscopic water droplets surrounded by a soft-matter boundary. Femtosecond nonlinear infrared spectroscopy in the OH-stretching region of H2O fully confirms the core/shell model, in which the entrapped water molecules partition onto two molecular subensembles: a bulk-like water core and a hydration layer near the ionic surfactant headgroups. These two distinct water species display different relaxation kinetics, as they do not exchange vibrational energy. The observed spectrotemporal ultrafast response exhibits a local character, indicating that the spatial confinement influences approximately one molecular layer located near the water-amphiphile boundary. The core of the encapsulated water droplet is similar in its spectroscopic properties to the bulk phase of liquid water, i.e., it does not display any true confinement effects such as droplet-size-dependent vibrational lifetimes or rotational correlation times. Unlike in bulk water, no intermolecular transfer of OH-stretching quanta occurs among the interfacial water molecules or from the hydration shell to the bulk-like core, indicating that the hydrogen bond network near the H2O/AOT interface is strongly disrupted.     

Water dynamics--the effects of ions and nanoconfinement

Hydrogen bond dynamics of water in highly concentrated NaBr salt solutions and reverse micelles are studied using ultrafast 2D-IR vibrational echo spectroscopy and polarization-selective IR pump-probe experiments performed on the OD hydroxyl stretch of dilute HOD in H(2)O. The vibrational echo experiments measure spectral diffusion, and the pump-probe experiments measure orientational relaxation. Both experimental observables are directly related to the structural dynamics of water's hydrogen bond network. The measurements performed on NaBr solutions as a function of concentration show that the hydrogen bond dynamics slow as the NaBr concentration increases. The most pronounced change is in the longest time scale dynamics which are related to the global rearrangement of the hydrogen bond structure. Complete hydrogen bond network randomization slows by a factor of approximately 3 in approximately 6 M NaBr solution compared to that in bulk water. The hydrogen bond dynamics of water in nanoscopically confined environments are studied by encapsulating water molecules in ionic head group (AOT) and nonionic head group (Igepal CO 520) reverse micelles. Water dynamics in the nanopools of AOT reverse micelles are studied as a function of size by observing orientational relaxation. Orientational relaxation dynamics deviate significantly from bulk water when the size of the reverse micelles is smaller than several nm and become nonexponential and slower as the size of the reverse micelles decreases. In the smallest reverse micelles, orientational relaxation (hydrogen bond structural randomization) is almost 20 times slower than that in bulk water. To determine if the changes in dynamics from bulk water are caused by the influence of the ionic head groups of AOT or the nanoconfinement, the water dynamics in 4 nm nanopools in AOT reverse micelles (ionic) and Igepal reverse micelles (nonionic) are compared. It is found that the water orientational relaxation in the 4 nm diameter nanopools of the two types of reverse micelles is almost identical, which indicates that confinement by an interface to form a nanoscopic water pool is a primary factor governing the dynamics of nanoscopic water rather than the presence of charged groups at the interface.