Hydration structure of water confined between mica surfaces

We report further molecular dynamics simulations on the structure of bound hydration layers under extreme confinement between mica surfaces. We find that the liquid phase of water is maintained down to 2 monolayer (ML) thick, whereas the structure of the K(+) ion hydration shell is close to the bulk structure even under D = 0.92 nm confinement. Unexpectedly, the density of confined water remains approximately the bulk value or less, whereas the diffusion of water molecules decreases dramatically. Further increase in confinement leads to a transition to a bilayer ice, whose density is much less than that of ice Ih due to the formation of a specific hydrogen-bonding network.     

Structure of methanol confined in MCM-41 investigated by large-angle X-ray scattering technique

Large-angle X-ray scattering (LAXS) measurements over a temperature range from 223 to 298 K have been made on methanol confined in mesoporous silica MCM-41 with two different pore diameters, 28 A (C14) and 21 A (C10), under both monolayer and capillary-condensed adsorption conditions. To compare the structure of methanol in the MCM-41 pores with that of bulk methanol, X-ray scattering intensities for bulk methanol in the same temperature range have also been measured. The radial distribution functions (RDFs) for the monolayer methanol samples showed that methanol molecules are strongly hydrogen bonded to the silanol groups on the MCM-41 surface, resulting in no significant change in the structure of adsorbed methanol with respect to the pore size and temperature. On the other hand, the RDFs for the capillary-condensed methanol samples showed that hydrogen-bonded chains of methanol molecules are formed in both pores. However, the distance and number of hydrogen bonds estimated from the RDFs suggested that hydrogen bonds between methanol molecules in the pores are significantly distorted or partly disrupted. It has been found that the hydrogen bonds are more distorted in the smaller pores of MCM-41. With decreasing temperature, however, the hydrogen-bonded chains of methanol in the pores were gradually ordered. A comparison of the present results on methanol in MCM-41 pores with those on water in the same pores revealed that the structural change with temperature is less significant for confined methanol than for confined water.     

Gas solubility in hydrophobic confinement

Measured forces between apolar surfaces in water have often been found to be sensitive to exposure to atmospheric gases despite low gas solubilities in bulk water. This raises questions as to how significant gas adsorption is in hydrophobic confinement, whether it is conducive to water depletion at such surfaces, and ultimately if it can facilitate the liquid-to-gas phase transition in the confinement. Open Ensemble molecular simulations have been used here to determine saturated concentrations of atmospheric gases in water-filled apolar confinements as a function of pore width at varied gas fugacities. For paraffin-like confinements of widths barely exceeding the mechanical instability threshold (spinodal) of the liquid-to-vapor transition of confined water (aqueous film thickness between three and four molecular diameters), mean gas concentrations in the pore were found to exceed the bulk values by a factor of approximately 30 or approximately 15 in cases of N2 and CO2, respectively. At ambient conditions, this does not result in visible changes in the water density profile next to the surfaces. Whereas the barrier to capillary evaporation has been found to decrease in the presence of dissolved gas (Leung, K.; Luzar, A.; and Bratko, D. Phys. Rev. Lett. 2003, 90, 065502), gas concentrations much higher than those observed at normal atmospheric conditions would be needed to produce noticeable changes in the kinetics of capillary evaporation. In simulations, dissolved gas concentrations corresponding to fugacities above approximately 40 bar for N2, or approximately 2 bar for CO2, were required to trigger expulsion of water from a hydrocarbon slit as narrow as 1.4 nm. For nanosized pore widths corresponding to the mechanical instability threshold or above, no significant coupling between adsorption layers at opposing confinement walls was observed. This finding explains the approximately linear increase in gas solubility with inverse confinement width and the apparent validity of Henry's law in the pores over a broad fugacity range.     

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.     

Melting behavior of water in cylindrical pores: carbon nanotubes and silica glasses

We report a study of the effects of confinement in multi-walled carbon nanotubes and mesoporous silica glasses (SBA-15) on the solid structure and melting of both H(2)O and D(2)O ice, using differential scanning calorimetry, dielectric relaxation spectroscopy, and neutron diffraction. Multi-walled nanotubes of 2.4, 3.9 and 10 nm are studied, and the SBA-15 studied has pores of mean diameter 3.9 nm; temperatures ranging from approximately 110 to 290 K were studied. We find that the melting point is depressed relative to the bulk water for all systems studied, with the depression being greater in the case of the silica mesopores. These results are shown to be consistent with molecular simulation studies of freezing in silica and carbon materials. The neutron diffraction data show that the cubic phase of ice is stabilized by the confinement in carbon nanotubes, as well as in silica mesopores, and persists up to temperatures of about 240 K, above which there is a transition to the hexagonal ice structure.     

Water confined in nanotubes and between graphene sheets: a first principle study

Water confined at the nanoscale has been the focus of numerous experimental and theoretical investigations in recent years, yet there is no consensus on such basic properties as diffusion and the nature of hydrogen bonding (HB) under confinement. Unraveling these properties is important to understand fluid flow and transport at the nanoscale, and to shed light on the solvation of biomolecules. Here we report on a first principle, computational study focusing on water confined between prototypical nonpolar substrates, i.e., single-wall carbon nanotubes and graphene sheets, 1-2.5 nm apart. The results of our molecular dynamics simulations show the presence of a thin, interfacial liquid layer (approximately 5 A) whose microscopic structure and thickness are independent of the distance between confining layers. The properties of the HB network are very similar to those of the bulk outside the interfacial region, even in the case of strong confinement. Our findings indicate that the perturbation induced by the presence of confining media is extremely local in liquid water, and we propose that many of the effects attributed to novel phases under confinement are determined by subtle electronic structure rearrangements occurring at the interface with the confining medium.     

Effects of confinement on small water clusters structure and proton transport

Analyses of the structure of two to four water molecule clusters confined between two benzene and between two naphthalene molecules have been performed using ab initio methods. The water clusters tend to maximize the number of hydrogen bonds via formation of a cyclic network. The oxygen atoms locate approximately in the middle of the confined geometry, and the dipole vectors arrange either parallel or pointing to the surfaces. Energy barriers for proton transfer calculated for H3O+-(H2O) complexes in the same confined geometries suggest that there is a specific range of confinement that helps to lower the energy barriers of the proton transfer. When the walls are too close to each other, at a separation of 4 A, the energy barriers are extremely high. Confinement does not lower the barrier energies of proton transfer when the H3O+-(H2O) complexes are located further from each of the surfaces by more than approximately 8 A.     

Effect of salts on water viscosity in narrow membrane pores

The effect of various salts on the viscosity, and by implication structure, of water in polymeric membrane pores of radius approximately 1.69 nm and low charge density has been studied. Permeation of pure water and various electrolyte solutions was analyzed using the Hagen-Poiseuille equation expressed in a ratio form to exclude membrane-specific quantities such as pore radius and length. The analysis produced viscosity ratios of electrolyte to pure water inside the membrane pores. Comparing the viscosity ratios inside the pores with their bulk counterparts showed that confinement significantly increased the sensitivity of water structure to the presence of ions. It has been found that, in relative terms in the pores, Cl- was a strong structure breaker, K+ was a moderately strong structure breaker, Na+ was a weak structure breaker, SO4(2-) was a weak structure maker, and Mg2+ was a strong structure maker. Predictive modeling of membrane separation performance would benefit from such effects being taken into account in cases where the pore ion concentrations may be high.     

Ab Initio H2O in Realistic Hydrophilic Confinement

A protocol for the ab initio construction of a realistic cylindrical pore in amorphous silica, serving as a geometric nanoscale confinement for liquids and solutions, is presented. Upon filling the pore with liquid water at different densities, the structure and dynamics of the liquid inside the confinement can be characterized. At high density, the pore introduces long‐range oscillations into the water density profile, which makes the water structure unlike that of the bulk across the entire pore. The tetrahedral structure of water is also affected up to the second solvation shell of the pore wall. Furthermore, the effects of the confinement on hydrogen bonding and diffusion, resulting in a weakening and distortion of the water structure at the pore walls and a slowdown in diffusion, are characterized.     

Adsorption of water molecules inside a Au nanotube: a molecular dynamics study

A molecular dynamics simulation of water molecules through a Au nanotube with a diameter of 20 A at bulk densities 0.8, 1, and 1.2 gcm(3) has been carried out. The water molecules inside a nanoscale tube, unlike those inside a bulk tube, have a confined effect. The interaction energy of the Au nanotube wall has a direct influence on the distribution of water molecules inside the Au tube in that the adsorption of the water molecules creates shell-like formations of water. Moreover, the high number of adsorbed molecules has already achieved saturation at the wall of the Au nanotube at three bulk densities. This work compares the distribution percentage profiles of hydrogen bonds for different regions inside the tube. The structural characteristics of water molecules inside the tube have also been studied. The results reveal that the numbers of hydrogen bonds per water molecule influence the orientational order parameter q. In addition, the phenomenon of a group of molecules bonded inside the tube can be observed as the number of hydrogen bonds increase.     

Radiowave dielectric investigation of water confined in channels of carbon nanotubes

Structure and dynamics of water confined in channels of diameter of few nanometer in size strongly differ from the ones of water in the bulk phase. Here, we present radiowave dielectric relaxation measurements on water-filled single-walled carbon nanotubes, with the aim of highlighting some aspects on the molecular electric dipole organization of water responding to high spatial confinement in a hydrophobic environment. The observed dielectric spectra, resulting into two contiguous relaxation processes, allow us to separate the confined water in the interior of the nanotubes from external water, providing support for the existence in the confinement region of water domains held together by hydrogen bonds. Our results, based on the deconvolution of the dielectric spectra due to the presence of a bulk and a confined water phase, furnish a significantly higher Kirkwood correlation factor, larger than the one of water in bulk phase, indicating a strong correlation between water molecules inside nanotubes, not seen in bulk water.     

Thermodynamic functions of water and ice confined to 2 nm radius pores

The heat capacity C(p) of the liquid state of water confined to 2 nm radius pores in Vycor glass was measured by temperature modulation calorimetry in the temperature range of 253-360 K, with an accuracy of 0.5%. On nanoconfinement, C(p) of water increases, and the broad minimum in the C(p) against T plot shifts to higher temperature. The increase in the C(p) of water is attributed to an increase in the phonon and configurational contributions. The apparent heat capacity of the liquid and partially frozen state of confined water was measured by temperature scanning calorimetry in the range of 240-280 K with an accuracy of 2%, both on cooling or heating at 6 K h(-1) rate. The enthalpy, entropy, and free energy of nanoconfined liquid water have been determined. The apparent heat capacity remains higher than that of bulk ice at 240 K and it is concluded that freezing is incomplete at 240 K. This is attributed to the intergranular-water-ice equilibrium in the pores. The nanoconfined sample melts over a 240-268 K range. For 9.6 wt % nanoconfined water concentration ( approximately 50% of the maximum filling) at 280 K, the enthalpy of water is 81.6% of the bulk water value and the entropy is 88.5%. For 21.1 wt % (100% filling) the corresponding values are 90.7% and 95.0%. The enthalpy decrease on nanoconfinement is a reflection of the change in the H-bonded structure of water. The use of the Gibbs-Thomson equation for analyzing the data has been discussed and it is found that a distribution of pore size does not entirely explain our results.     

Thermal properties and vibrational density of states of a nanoconfined discotic liquid crystal

Neutron scattering is employed to investigate the vibrational density of states (VDOS) of the discotic liquid crystal 2,3,6,7,10,11-hexakis[hexyloxy] triphenylene (HAT6) confined to the pores of alumina oxide membranes with different pore sizes. Additionally, the phase transitions were studied by differential scanning calorimetry. The transitions were observed down to the smallest pore size. The decrease of the transition enthalpies versus inverse pore size for both transitions implies an increase of the amount of disordered amorphous material. By extrapolation of its pore size dependence, a critical pore diameter for structure formation of 17 nm is estimated. Similar to the bulk, excess contributions to the VDOS (Boson peak) are also observed for confined HAT6. The Boson peak gains in intensity and shifts to lower frequencies with decreasing pore diameter. This is discussed in the framework of a softening of HAT6 induced by the confinement due to a less-developed plastic crystalline state inside the pores compared to the bulk.     

Effect of confinement on the liquid-liquid phase transition of supercooled water

We report on an observation of the phase transition between two liquid phases of supercooled confined water in simulations. The temperature of the liquid-liquid transition of water at zero pressure slightly decreases due to confinement in the hydrophobic pore. The hydrophilic confinement affects this temperature in the opposite direction and shifts the critical point of the liquid-liquid transition to a higher pressure. As a result, in a strongly hydrophilic pore the liquid-liquid phase transition becomes continuous at zero pressure, indicating the shift of its critical point from negative to a positive pressure. These findings indicate that experimental studies of water confined in the pores of various hydrophobicity/hydrophilicity may clarify the location of the liquid-liquid critical point of bulk water.     

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.     

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.     

Does water condense in carbon pores?

Using grand canonical Monte Carlo (GCMC) simulations of molecular models, we investigate the nature of water adsorption and desorption in slit pores with graphitelike surfaces. Special emphasis is placed on the question of whether water exhibits capillary condensation (i.e., condensation when the external pressure is below the bulk vapor pressure). Three models of water have been considered. These are the SPC and SPC/E models and a model where the hydrogen bonding is described by tetrahedrally coordinated square-well association sites. The water-carbon interaction was described by the Steele 10-4-3 potential. In addition to determining adsorption/desorption isotherms, we also locate the states where vapor-liquid equilibrium occurs for both the bulk and confined states of the models. We find that for wider pores (widths >1 nm), condensation does not occur in the GCMC simulations until the pressure is higher than the bulk vapor pressure, P0. This is consistent with a physical picture where a lack of hydrogen bonding with the graphite surface destabilizes dense water phases relative to the bulk. For narrow pores where the slit width is comparable to the molecular diameter, strong dispersion interactions with both carbon surfaces can stabilize dense water phases relative to the bulk so that pore condensation can occur for P < P0 in some cases. For the narrowest pores studied--a pore width of 0.6 nm--pore condensation is again shifted to P > P0. The phase-equilibrium calculations indicate vapor-liquid coexistence in the slit pores for P < P0 for all but the narrowest pores. We discuss the implications of our results for interpreting water adsorption/desorption isotherms in porous carbons.     

Density minimum of confined water at low temperatures: a combined study by small-angle scattering of X-rays and neutrons

A simple explanation is given for the low-temperature density minimum of water confined within cylindrical pores of ordered nanoporous materials of different pore size. The experimental evidence is based on combined data from in-situ small-angle scattering of X-rays (SAXS) and neutrons (SANS), corroborated by additional wide-angle X-ray scattering (WAXS). The combined scattering data cannot be described by a homogeneous density distribution of water within the pores, as was originally suggested from SANS data alone. A two-step density model reveals a wall layer covering approximately two layers of water molecules with higher density than the residual core water in the central part of the pores. The temperature-induced changes of the scattering signal from both X-rays and neutrons are consistent with a minimum of the average water density. We show that the temperature at which this minimum occurs depends monotonically on the pore size. Therefore we attribute this minimum to a liquid-solid transition of water influenced by confinement. For water confined in the smallest pores of only 2 nm in diameter, the density minimum is explained in terms of a structural transition of the surface water layer closest to the hydrophilic pore walls.     

High-density hydration layer of lysozymes: molecular dynamics decomposition of solution scattering data

A characterization of the physical properties of protein hydration water is critical for understanding protein structure and function. Recent small-angle X-ray and neutron scattering data indicate that the density of water on the surface of lysozyme is significantly higher than in bulk water. Here, we provide an interpretation of the scattering results using a molecular dynamics simulation, which allows us to make quantitative predictions about density variations in the first hydration shell. The perturbation relative to bulk water involves statistically significant changes in the average water structure in the first hydration layer. The water density in the first hydration shell is increased by 5% with respect to the bulk. In regions of higher water density, the water dipoles align more parallel to each other and the number of hydrogen bonds per water molecule is higher. Increased water density is found for water molecules interacting with hydrogen and carbon atoms in the backbone or with nonpolar or negatively charged side-chain groups.     

Neutron scattering study on dynamics of water molecules in MCM-41. 2. Determination of translational diffusion coefficient

Quasielastic neutron scattering (QENS) spectra of water-filled MCM-41 samples (pore diameters: 21.4 and 28.4 Angstrom) were measured over the temperature range 238-298 K and the momentum transfer range 0.31-0.99 A(-1) to investigate the dynamics of confined water molecules. The spectra, which consist mainly of contributions from the translational diffusion of water molecules, were analyzed by using the Lorentzian and the stretched exponential functions. Comparison of the fits indicated that the latter analysis is more reliable than the former one. The fraction of immobile water molecules located in the vicinity of the pore walls, which give an elastic component, was found to be 0.044-0.061 in both pores. The stretch exponent beta was determined as 0.66-0.80. It was shown that the translational diffusion of water molecules in the pores is decelerated by confinement and that the deceleration becomes marked with a decrease in pore size. The ratios of the translational diffusion coefficient D(T) of confined water to that of bulk water at room temperature were within a range of 0.47-0.63.