Solvation forces and liquid–solid phase equilibria for water confined between hydrophobic surfaces

How do water-mediated interactions and osmotic second virial coefficients vary with particle size?

We examine quantitatively the solute-size dependences of the effective interactions between nonpolar solutes in water and in a simple liquid. The potential w(r) of mean force and the osmotic second virial coefficients B are calculated with high accuracy from molecular dynamics simulations. As the solute diameter increases from methane's to C60's with the solute-solute and solute-solvent attractive interaction parameters fixed to those for the methane-methane and methane-water interactions, the first minimum of w(r) lowers from -1.1 to -4.7 in units of the thermal energy kT. Correspondingly, the magnitude of B (<0) increases proportional to σα with some power close to 6 or 7, which reinforces the solute-size dependence of B found earlier for a smaller range of σ [H. Naito, R. Okamoto, T. Sumi and K. Koga, J. Chem. Phys., 2022, 156, 221104]. We also demonstrate that the strength of the attractive interactions between solute and solvent molecules can qualitatively change the characteristics of the effective pair interaction between solute particles, both in water and in a simple liquid. If the solute-solvent attractive force is set to be weaker (stronger) than a threshold, the effective interaction becomes increasingly attractive (repulsive) with increasing solute size.

Molecular Shape and the Hydrophobic Effect

This review focuses on papers published since 2000 on the topic of the properties of solutes in water. More specifically, it evaluates the state of the art of our understanding of the complex relationship between the shape of a hydrophobe and the hydrophobic effect. To highlight this, we present a selection of references covering both empirical and molecular dynamics studies of small (molecular-scale) solutes. These include empirical studies of small molecules, synthetic hosts, crystalline monolayers, and proteins, as well as in silico investigations of entities such as idealized hard and soft spheres, small solutes, hydrophobic plates, artificial concavity, molecular hosts, carbon nanotubes and spheres, and proteins.

Dynamics and energetics of hydrophobically confined water

The effects of water confined in regions between self-assembling entities is relevant to numerous contexts such as macromolecular association, protein folding, protein-ligand association, and nanomaterials self-assembly. Thus assessing the impact of confined water, and the ability of current modeling techniques to capture the salient features of confined water is important and timely. We present molecular dynamics simulation results investigating the effect of confined water on qualitative features of potentials of mean force describing the free energetics of self-assembly of large planar hydrophobic plates. We consider several common explicit water models including the TIP3P, TIP4P, SPC/E, TIP4P-FQ, and SWM4-NDP, the latter two being polarizable models. Examination of the free energies for filling and unfilling the volume confined between the two plates (both in the context of average number of confined water molecules and "depth" of occupancy) suggests TIP4P-FQ water molecules generally occupy the confined volume at separation distances larger than observed for other models under the same conditions. The connection between this tendency of TIP4P-FQ water and the lack of a pronounced barrier in the potential of mean force for plate-plate association in TIP4P-FQ water is explored by artificially, but systematically, populating the confined volume with TIP4P-FQ water at low plate-plate separation distances. When the critical separation distance [denoting the crossover from an unoccupied (dry) confined interior to a filled (wet) interior] for TIP4P-FQ is reduced by 0.5 Å using this approach, a barrier is observed; we rationalize this effect based on increased resistant forces introduced by confined water molecules at these low separations. We also consider the dynamics of water molecules in the confined region between the hydrophobes. We find that the TIP4P-FQ water model exhibits nonbulklike dynamics, with enhanced lateral diffusion relative to bulk. This is consistent with the reduced intermolecular water-water interaction indicated by a decreased molecular dipole moment in the interplate region. Analysis of velocity autocorrelation functions and associated power spectra indicate that the interplate region for TIP4P-FQ at a plate separation of 14.4 Å approaches characteristics of the pure water liquid-vapor interface. This is in stark contrast to the other water models (including the polarizable SWM4-NDP model).

Molecular simulation of water confined in nanoporous silica

This paper reports on a molecular simulation study of the thermodynamics, structure and dynamics of water confined at ambient temperature in hydroxylated silica nanopores of a width H = 10 and 20 Å. The adsorption isotherms for water in these nanopores resemble those observed for experimental samples; the adsorbed amount increases continuously in the multilayer adsorption regime until a jump occurs due to capillary condensation of the fluid within the pore. Strong layering of water in the vicinity of the silica surfaces is observed as marked density oscillations are observed up to 8 Å from the surface in the density profiles for confined water. Our results indicate that water molecules within the first adsorbed layer tend to adopt a H-down orientation with respect to the silica substrate. For all pore sizes and adsorbed amounts, the self-diffusivity of confined water is lower than the bulk, due to the hydrophilic interaction between the water molecules and the hydroxylated silica surface. Our results also suggest that the self-diffusivity of confined water is sensitive to the adsorbed amount.

Effect of pressure on the freezing of pure fluids and mixtures confined in nanopores

Monte Carlo simulations combined with the parallel tempering technique are used to study the freezing of Ar, CH(4), and their mixtures in a slit graphite nanopore. For all systems, the solid/liquid coexistence line is located at higher temperature than that for the bulk phase, as expected for fluids for which the wall/fluid interaction is stronger than the fluid/fluid interaction. In the case of the mixtures, the phase diagram for the confined system is of the same type as that for the bulk (azeotropic). It is also found that the freezing temperatures for the confined fluids and mixture are much more affected by pressure than those for the bulk phase. By calculating the isothermal compressibility of the confined fluids and determining the slope of the solid/liquid coexistence line (T,P) from the Clapeyron equation, we show that such a strong effect of pressure is not related to reduced compressibility within the pores. On the other hand, the pressure dependence of the in-pore freezing temperature is correctly described in the frame of the model proposed by Miyahara et al. [ Miyahara , M. ; Kanda , H. ; Shibao , M. ; Higashitani , K. J. Chem. Phys. 2000 , 112 , 9909. ], which is based on the pressure difference between the bulk and confined phases (capillary effect). In this model, a change in the in-pore freezing temperature with pressure is explained by a drastic change in the in-pore pressure, which varies very sharply with the bulk external pressure. We present an extended version of this model to confined systems for which an increase in the freezing temperature is observed.

Energy effects on the structure and thermodynamic properties of nanoconfined fluids (a density functional theory study)

The structure and properties of fluids confined in nanopores may show a dramatic departure from macroscopic bulk fluids. The main reason for this difference lies in the influence of system walls. In addition to the entropic wall effect, system walls can significantly change the energy of the confined fluid compared to macroscopic bulk fluids. The energy effect of the walls on a nanoconfined fluid appears in two forms. The first effect is the cutting off of the intermolecular interactions by the walls, which appears for example in the integrals for calculation of the thermodynamic properties. The second wall effect involves the wall-molecule interactions. In such confined fluids, the introduction of wall forces and the competition between fluid-wall and fluid-fluid forces could lead to interesting thermodynamic properties, including new kinds of phase transitions not observed in the macroscopic fluid systems. In this article, we use the perturbative fundamental measure density functional theory to study energy effects on the structure and properties of a hard core two-Yukawa fluid confined in a nanoslit. Our results show the changes undergone by the structure and phase transition of the nanoconfined fluids as a result of energy effects.

Phase transitions induced by nanoconfinement in liquid water

We present results from molecular dynamics simulations of water confined by two parallel atomically detailed hydrophobic walls. Simulations are performed at T = 300 K and wall-wall separation d = 0.6-1.6 nm. At 0.7 < or = d < or = 0.9 nm, a first order transition occurs between a bilayer liquid (BL) and a trilayer heterogeneous fluid (THF) as water density increases. The THF is characterized by a liquid (central) layer and two crystal-like layers next to the walls. The BL-THF transition involves freezing of the two surface layers in contact with the walls. At d = 0.6 nm, the THF transforms into a bilayer ice (BI) upon decompression. Both the BL-THF and BI-THF transitions are induced by the surface regular atomic-scale structure.

A computer simulation study of stick-slip transitions in water films confined between model hydrophilic surfaces. 1. Monolayer films

The shear behavior of monolayer water films confined in a slit-like pore between hydrophilic walls is simulated in the quasistatic regime using the grand canonical Monte Carlo technique. Each wall is represented as a hexagonal lattice of force sites that interact with water through an orientation-dependent hydrogen-bonding potential. When the walls are in registry, the water oxygen atoms form either a crystal- or fluid-like structure, depending on the period of the wall's lattice. In both cases, however, the monolayer structure is orientationally disordered. Both the crystal- and fluid-like monolayers prove to be capable of experiencing well-defined stick-slip transitions, with the largest yield stress occurring in the crystal-like case. Beyond the yield point, the crystal-like monolayers "melt", but their structure and molecular motion differ in many respects from those characteristic of normal fluids.

The hydrophobic effect: molecular dynamics simulations of water confined between extended hydrophobic and hydrophilic surfaces

Structural and dynamic properties of water confined between two parallel, extended, either hydrophobic or hydrophilic crystalline surfaces of n-alkane C(36)H(74) or n-alcohol C(35)H(71)OH, are studied by molecular dynamics simulations. Electron density profiles, directly compared with corresponding experimental data from x-ray reflectivity measurements, reveal a uniform weak de-wetting characteristic for the extended hydrophobic surface, while the hydrophilic surface is weakly wetted. These microscopic data are consistent with macroscopic contact angle measurements. Specific water orientation is present at both surfaces. The ordering is characteristically different between the surfaces and of longer range at the hydrophilic surface. Furthermore, the dynamic properties of water are different at the two surfaces and different from the bulk behavior. In particular, at the hydrophobic surface, time-correlation functions reveal that water molecules have characteristic diffusive behavior and orientational ordering due to the lack of hydrogen bonding interactions with the surface. These observations suggest that the altered dynamical properties of water in contact with extended hydrophobic surfaces together with a partial drying of the surfaces are more indicative of the hydrophobic effect than structural ordering, which we suggest to be independent of surface topology.

Interfacial water at hydrophobic and hydrophilic surfaces: depletion versus adsorption

The structure of the water-solid interface for widely varying surface properties is investigated with Monte Carlo simulations using the SPC/E water model. Of particular interest is the relation between the wetting coefficient as a measure of the hydrophobicity of the substrate and the density depletion close to the solid surface. The substrates are modeled as rigid ordered lattices of sites that interact with water molecules through an orientation-independent Lennard-Jones potential of varying strength. Hydrophilic character is obtained by addition of polar hydroxyl groups on the substrate surface, and the influence of density, spatial distribution, and angular orientation of the polar groups on the interfacial water structure is studied.

Temperature and structural changes of water clusters in vacuum due to evaporation

This paper presents a study on evaporation of pure water clusters. Molecular dynamics simulations between 20 ns and 3 micros of clusters ranging from 125 to 4096 molecules in vacuum were performed. Three different models (SPC, TIP4P, and TIP5P) were used to simulate water, starting at temperatures of 250, 275, and 300 K. We monitored the temperature, the number of hydrogen bonds, the tetrahedral order, the evaporation, the radial distribution functions, and the diffusion coefficients. The three models behave very similarly as far as temperature and evaporation are concerned. Clusters starting at a higher temperature show a higher initial evaporation rate and therefore reach the point where evaporation stop (around 240 K) sooner. The radius of the clusters is decreased by 0.16-0.22 nm after 0.5 micros (larger clusters tend to decrease their radius slightly more), which corresponds to around one evaporated molecule per nm(2). The cluster temperature seems to converge towards 215 K independent of cluster size, when starting at 275 K. We observe only small structural changes, but the clusters modeled by TIP5P show a larger percentage of molecules with low diffusion coefficient as t-->infinity, than those using the two other water models. TIP4P seems to be more structured and more hydrogen bonds are formed than in the other models as the temperature falls. The cooling rates are in good agreement with experimental results, and evaporation rates agree well with a phenomenological expression based on experimental observations.

Effects of confinement on freezing and melting

We present a review of experimental, theoretical, and molecular simulation studies of confinement effects on freezing and melting. We consider both simple and more complex adsorbates that are confined in various environments (slit or cylindrical pores and also disordered porous materials). The most commonly used molecular simulation, theoretical and experimental methods are first presented. We also provide a brief description of the most widely used porous materials. The current state of knowledge on the effects of confinement on structure and freezing temperature, and the appearance of new surface-driven and confinement-driven phases are then discussed. We also address how confinement affects the glass transition.

Phase diagram of water between hydrophobic surfaces

Molecular dynamics simulations demonstrate that there are at least two classes of quasi-two-dimensional solid water into which liquid water confined between hydrophobic surfaces freezes spontaneously and whose hydrogen-bond networks are as fully connected as those of bulk ice. One of them is the monolayer ice and the other is the bilayer solid which takes either a crystalline or an amorphous form. Here we present the phase transformations among liquid, bilayer amorphous (or crystalline) ice, and monolayer ice phases at various thermodynamic conditions, then determine curves of melting, freezing, and solid-solid structural change on the isostress planes where temperature and intersurface distance are variable, and finally we propose a phase diagram of the confined water in the temperature-pressure-distance space.

Water at hydrophobic substrates: curvature, pressure, and temperature effects

We studied the water density profile close to spherical and planar hydrophobic objects using molecular dynamics (MD) simulations. For normal pressure and room temperature, the depletion layer thickness of a planar substrate is approximately 2.5 Angstroms. Even for quite large spherical solutes with a radius of R = 18 Angstroms, the depletion layer thickness is reduced by 30%, which shows that substrate curvature and roughness is an experimentally important factor. Rising temperature leads to a substantial increase of the depletion layer thickness. The compressibility of the depletion layer is found to be surprisingly small and only approximately 5 times higher than that of bulk water. A high electrostatic surface potential of 0.5 V is found, which presumably plays an important role in the presence of charged solutes, since it can promote adsorption into the interfacial layer.

Dynamical behavior of one-dimensional water molecule chains in zeolites: Nanosecond time-scale molecular dynamics simulations of bikitaite

Nanosecond scale molecular dynamics simulations of the behavior of the one-dimensional water molecule chains adsorbed in the parallel nanochannels of bikitaite, a rare lithium containing zeolite, were performed at different temperatures and for the fully and partially hydrated material. New empirical potential functions have been developed for representing lithium-water interactions. The structure and the vibrational spectrum of bikitaite were in agreement both with experimental data and Car-Parrinello molecular dynamics results. Classical molecular dynamics simulations were extended to the nanosecond time scale in order to study the flip motion of water molecules around the hydrogen bonds connecting adjacent molecules in the chains, which has been observed by NMR experiments, and the dehydration mechanism at high temperature. Computed relaxation times of the flip motion follow the Arrhenius behavior found experimentally, but the activation energy of the simulated system is slightly underestimated. Based on the results of the simulations, it may be suggested that the dehydration proceeds by a defect-driven stepwise diffusion. The diffusive mechanism appears as a single-file motion: the molecules never pass one another, even at temperatures as high as about 1000 K, nor can they switch between different channels. However, the mean square displacement (MSD) of the molecules, computed with respect to the center of mass of the simulated system, shows an irregular trend from which the single-file diffusion cannot be clearly evidenced. If the MSDs are evaluated with respect to the center of mass of the molecules hosted in each channel, the expected dependence on the square root of time finally appears.