Phase diagram of water between hydrophobic surfaces

Ferroelectric hexagonal and rhombic monolayer ice phases

Two new monolayer ice phases are predicted from molecular dynamics simulations, both proven to be ferroelectric.

Model atomistic protrusions favouring the ordering and retention of water

The ordering of water molecules near model linear atomistic protrusions is studied using classical molecular dynamics simulations.

Direct observation of self-assembled chain-like water structures in a nanoscopic water meniscus

Sawtooth-like oscillatory forces generated by water molecules confined between two oxidized silicon surfaces were observed using a cantilever-based optical interfacial force microscope when the two surfaces approached each other in ambient environments. The humidity-dependent oscillatory amplitude and periodicity were 3-12 nN and 3-4 water diameters, respectively. Half of each period was matched with a freely jointed chain model, possibly suggesting that the confined water behaved like a bundle of water chains. The analysis also indicated that water molecules self-assembled to form chain-like structures in a nanoscopic meniscus between two hydrophilic surfaces in air. From the friction force data measured simultaneously, the viscosity of the chain-like water was estimated to be between 10(8) and 10(10) times greater than that of bulk water. The suggested chain-like structure resolves many unexplained properties of confined water at the nanometer scale, thus dramatically improving the understanding of a variety of water systems in nature.

Electromelting of confined monolayer ice

In sharp contrast to the prevailing view that electric fields promote water freezing, here we show by molecular dynamics simulations that monolayer ice confined between two parallel plates can melt into liquid water under a perpendicularly applied electric field. The melting temperature of the monolayer ice decreases with the increasing strength of the external field due to the field-induced disruption of the water-wall interaction induced well-ordered network of the hydrogen bond. This electromelting process should add an important new ingredient to the physics of water.

Thermodynamic anomalies of a network former in a periodic field: Network former in a periodic field

We study the effect of an external one-dimensional periodic field on thermodynamic anomalies associated with a two-dimensional model liquid with anisotropic interactions. The model system, a 50 : 50 binary mixture of two species of particles interacting with an angle-dependent Lennard-Jones potential, has a rich phase diagram and shows many features of network-forming liquids like water and silica such as a prominent minimum in the pressure-temperature isochore. Confining the system by a commensurate one-dimensional periodic field shifts the temperature of minimum pressure to higher temperatures. A mean-field theory of a lattice-gas in an external field which couples to internal orientational states reproduces these results.

Structuring of interfacial water on silica surface in cyclohexane studied by surface forces measurement and sum frequency generation vibrational spectroscopy

We investigated interfacial water, formed by adsorption or phase separation (prewetting transition), on a silica surface in water-cyclohexane binary liquids using a combination of colloidal probe atomic force microscopy (AFM) and sum frequency generation (SFG) vibrational spectroscopy. At 33 ± 9 ppm water, the long-range attraction extending to 19.4 ± 2.9 nm appeared, which was caused by the contact of water layers formed on silica surfaces. The attraction range increased with increasing water concentration and reached 97 ± 17 nm at the saturation concentration of water in cyclohexane (C*), indicating that the thickness of the water layer formed on silica was ca. 50 nm. The interfacial energy between the water adsorption layer and bulk solution (γ = 79.3 ± 2.0 mN/m) was estimated from the pull-off force, and was significantly larger than the value for the bulk water/cyclohexane interface (γ = 50.1 mN/m). SFG spectroscopy demonstrated that the interfacial water formed an icelike structure at C*. These results indicated that the interfacial water molecules formed an icelike ordered structure induced by the hydrogen bonding with surface silanol groups, resulting in the free OH groups being more exposed to the bulk solution. On the other hand, the water adsorption layer induced by phase separation at water concentrations above C* was found to be less ordered and its structure at the adsorption layer/bulk interface was almost the same as that of bulk water, although its thickness was almost the same as that formed at C*. To our knowledge, this is the first report of the observation of liquid adsorption layers formed by chemical interaction up to saturation and by the wetting transition above saturation, and their differences in the structure and properties at the molecular level.

Water on graphene surfaces

In this paper, we summarize the main results obtained in our group about the behavior of water confined inside or close to different graphene surfaces by means of molecular dynamics simulations. These include the inside and outside of carbon nanotubes, and the confinement inside a slit pore or a single graphene sheet. We paid special attention to some thermodynamical (binding energies), structural (hydrogen-bond distributions) and dynamic (infrared spectra) properties, and their comparison to their bulk counterparts.

Guest-free monolayer clathrate and its coexistence with two-dimensional high-density ice

Three-dimensional (3D) gas clathrates are ice-like but distinguished from bulk ices by containing polyhedral nano-cages to accommodate small gas molecules. Without space filling by gas molecules, standalone 3D clathrates have not been observed to form in the laboratory, and they appear to be unstable except at negative pressure. Thus far, experimental evidence for guest-free clathrates has only been found in germanium and silicon, although guest-free hydrate clathrates have been found, in recent simulations, able to grow from cold stretched water, if first nucleated. Herein, we report simulation evidence of spontaneous formation of monolayer clathrate ice, with or without gas molecules, within hydrophobic nano-slit at low temperatures. The guest-free monolayer clathrate ice is a low-density ice (LDI) whose geometric pattern is identical to Archimedean 4.8(2)-truncated square tiling, i.e. a mosaic of tetragons and octagons. At large positive pressure, a second phase of 2D monolayer ice, i.e. the puckered square high-density ice (HDI) can form. The triple point of the LDI/liquid/HDI three-phase coexistence resembles that of the ice-I(h)/water/ice-III three-phase coexistence. More interestingly, when the LDI is under a strong compression at 200 K, it transforms into the HDI via a liquid intermediate state, the first direct evidence of Ostwald's rule of stages at 2D. The tensile limit of the 2D LDI and water are close to that of bulk ice-I(h) and laboratory water.

Water modeled as an intermediate element between carbon and silicon

Water and silicon are chemically dissimilar substances with common physical properties. Their liquids display a temperature of maximum density, increased diffusivity on compression, and they form tetrahedral crystals and tetrahedral amorphous phases. The common feature to water, silicon, and carbon is the formation of tetrahedrally coordinated units. We exploit these similarities to develop a coarse-grained model of water (mW) that is essentially an atom with tetrahedrality intermediate between carbon and silicon. mW mimics the hydrogen-bonded structure of water through the introduction of a nonbond angular dependent term that encourages tetrahedral configurations. The model departs from the prevailing paradigm in water modeling: the use of long-ranged forces (electrostatics) to produce short-ranged (hydrogen-bonded) structure. mW has only short-range interactions yet it reproduces the energetics, density and structure of liquid water, and its anomalies and phase transitions with comparable or better accuracy than the most popular atomistic models of water, at less than 1% of the computational cost. We conclude that it is not the nature of the interactions but the connectivity of the molecules that determines the structural and thermodynamic behavior of water. The speedup in computing time provided by mW makes it particularly useful for the study of slow processes in deeply supercooled water, the mechanism of ice nucleation, wetting-drying transitions, and as a realistic water model for coarse-grained simulations of biomolecules and complex materials.

Ordering of hard spheres inside hard cylindrical pores

Isothermal-isobaric simulations on the ordering behavior of hard spheres upon confinement are presented. The radii of the confining cylinders go from 1.1 to 2 in units of the diameters of the hard spheres adsorbed. In all the range of pressures considered the spheres were located in concentric layers, as many as the radius of the hard cylinder would permit. When the pressure increases, the hard spheres go from being loosely arranged to the formation of ordered structures. This change is marked in all cases by a distinct break in the density of spheres in a narrow pressure range. When the tube radius is smaller than 1.5, the high-pressure ordering is determined by the number of coplanar spheres you can have within a circle of radius equal to that of the confining tube. For wider tubes, the change upon compression is determined by the formation of defected two-dimensional triangular lattices wrapped to fit inside the particular cylinder we are considering.

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.

Absence of a diffusion anomaly of water in the direction perpendicular to hydrophobic nanoconfining walls

We perform molecular dynamics simulations to investigate the diffusive motion of TIP5P (the transferable intermolecular potential with five points) water in the direction perpendicular to the two hydrophobic confining walls. To calculate the diffusion constant, we use the concept of the characteristic residence time which is calculated from the exponential decay of the residence time probability density function. We find that a diffusion anomaly, increase of diffusion upon compression, is absent in the direction perpendicular to the confining walls down to the lowest temperature we simulate, 220K , whereas there is a diffusion anomaly, similar to that in bulk water, in the direction parallel to the walls. The absence of a diffusion anomaly in the direction perpendicular to the walls may arise mainly due to nanoconfinement, rather than due to the hydrophobic property of the confining walls. In addition, we find that the temperature dependence of the diffusion constant along the constant density path in the perpendicular direction shows a Vogel-Fulcher-Tammann form.

Interfacial anisotropy in the transport of liquid crystals confined between flat, structureless walls: a molecular dynamics simulation approach

Molecular dynamics simulations of uniaxial Gay-Berne ellipsoids as prolate liquid crystal molecules confined between two flat, structureless walls have been carried out in order to investigate anisotropy in their dynamic properties. Several physical quantities are profiled as a function of distance from a wall. The walls stimulate ellipsoids into different behaviors from those of the bulk system. The profiles of self-diffusion coefficients, which are distinguished in each direction of a director-based coordinate system, show that the ellipsoids are more diffusive parallel to the walls and less diffusive perpendicular to the walls with decreasing distance from the walls. According to the self-rotation coefficient and rotational viscosity profiles, ellipsoids are easy to rotate parallel to the walls and hard to rotate in the plane perpendicular to the walls. The analyses of velocity autocorrelation functions, angular velocity autocorrelation functions, director angular velocity autocorrelation functions, and their spectra are useful for the investigation of anisotropy near the walls. We conclude that the flat, structureless wall not only prevents ellipsoids from diffusing and rotating in the plane perpendicular to the walls, but also stimulates them to diffuse and rotate in the plane parallel to the walls.

Phase equilibria and interfacial tension of fluids confined in narrow pores

Correlation between phase behaviors of a Lennard-Jones fluid in and outside a pore is examined over wide thermodynamic conditions by grand canonical Monte Carlo simulations. A pressure tensor component of the confined fluid, a variable controllable in simulation but usually uncontrollable in experiment, is related with the pressure of a bulk homogeneous system in equilibrium with the confined system. Effects of the pore dimensionality, size, and attractive potential on the correlations between thermodynamic properties of the confined and bulk systems are clarified. A fluid-wall interfacial tension defined as an excess grand potential is evaluated as a function of the pore size. It is found that the tension decreases linearly with the inverse of the pore diameter or width.

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.

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.

Effect of water-wall interaction potential on the properties of nanoconfined water

Much of the understanding of bulk liquids has progressed through study of the limiting case in which molecules interact via purely repulsive forces, such as a hard-core or "repulsive ramp" potential. In the same spirit, we report progress on the understanding of confined water by examining the behavior of waterlike molecules interacting with planar walls via purely repulsive forces and compare our results with those obtained for Lennard-Jones (LJ) interactions between the molecules and the walls. Specifically, we perform molecular dynamics simulations of 512 waterlike molecules interacting via the TIP5P potential and confined between two smooth planar walls that are separated by 1.1nm . At this separation, there are either two or three molecular layers of water, depending on density. We study two different forms of repulsive confinement, when the water-wall interaction potential is either (i) 1r;{9} or (ii) a WCA-like repulsive potential. We find that the thermodynamic, dynamic, and structural properties of the liquid in purely repulsive confinements qualitatively match those for a system with a pure LJ attraction to the wall. In previous studies that include attractions, freezing into monolayer or trilayer ice was seen for this wall separation. Using the same separation as these previous studies, we find that the crystal state is not stable with 1r;{9} repulsive walls but is stable with WCA-like repulsive confinement. However, by carefully adjusting the separation of the plates with 1r;{9} repulsive interactions so that the effective space available to the molecules is the same as that for LJ confinement, we find that the same crystal phases are stable. This result emphasizes the importance of comparing systems only using the same effective confinement, which may differ from the geometric separation of the confining surfaces.

Freezing of hard spheres confined in narrow cylindrical pores

Monte Carlo simulations for the equation of state and phase behavior of hard spheres confined inside very narrow hard tubes are presented. For pores whose radii are greater than 1.1 hard sphere diameters, a sudden change in the density and the microscopic structure of the fluid is neatly observed, indicating the onset of freezing. In the high-density structure the particles rearrange in such a way that groups of three particles fit in sections across the pore.

Formation and Decomposition of Ethane, Propane, and Carbon Dioxide Hydrates in Silica Gel Mesopores under High Pressure

The experimental data on decomposition temperatures for the gas hydrates of ethane, propane, and carbon dioxide dispersed in silica gel mesopores are reported. The studies were performed at pressures up to 1 GPa. It is shown that the experimental dependence of hydrate decomposition temperature on the size of pores that limit the size of hydrate particles can be described on the basis of the Gibbs-Thomson equation only if one takes into account changes in the shape coefficient that is present in the equation; in turn, the value of this coefficient depends on a method of mesopore size determination. A mechanism of hydrate formation in mesoporous medium is proposed. Experimental data providing evidence of the possibility of the formation of hydrate compounds in hydrophobic matrixes under high pressure are reported. Decomposition temperature of those hydrate compounds is higher than that for the bulk hydrates of the corresponding gases.

Capillary phase transitions of linear and branched alkanes in carbon nanotubes from molecular simulation

Capillary phase transitions of linear (from C(1) to C(12)) and branched (C(5) isomers) alkanes in single-walled carbon nanotubes have been investigated using the gauge-cell Monte Carlo simulation. The isotherm at a supercritical temperature increases monotonically with chemical potential and coincides with that from the traditional grand canonical Monte Carlo simulation, whereas the isotherm at a subcritical temperature exhibits a sigmoid van der Waals loop including stable, metastable, and unstable regions. Along this loop, the coexisting phases are determined using an Maxwell equal-area construction. A generic confinement effect is found that reduces the saturation chemical potential, lowers the critical temperature, increases the critical density, and shrinks the phase envelope. The effect is greater in a smaller diameter nanotube and is greater in a nanotube than in a nanoslit.

Effect of pressure on the phase behavior and structure of water confined between nanoscale hydrophobic and hydrophilic plates

We perform systematic molecular dynamics simulations of water confined between two nanoscale plates at T = 300K. We investigate the effect of pressure (-0.15 GPa< or = P< or =0.2GPa) and plate separation (0.4 nm < or =d < or =1.6 nm) on the phase behavior of water when the plates are either hydrophobic or hydrophilic. When water is confined between hydrophobic plates, capillary evaporation occurs between the plates at low enough P. The threshold value of d at which this transition occurs decreases with P (e.g., 1.6 nm at P approximately equal to -0.05 GPa, 0.5 nm at P approximately equal to 0.1 GPa), until, at high P, no capillary evaporation occurs. For d approximately equal to 0.6 nm and P > or =0.1 GPa, the system crystallizes into a bilayer ice. A P-d phase diagram showing the vapor, liquid, and bilayer ice phases is proposed. When water is confined by hydrophilic (hydroxylated silica) plates, it remains in the liquid phase at all P and d studied. Interestingly, we observe for this case that even at the P at which bulk water cavitates, the confined water remains in the liquid state. We also study systematically the state of hydration at different P for both kinds of plates. For the range of conditions studied here, we find that in the presence of hydrophobic plates the effect of P is to enhance water structure and to push water molecules toward the plates. The average orientation of water molecules next to the hydrophobic plates does not change upon pressurization. In contrast, in the presence of hydrophilic plates, water structure is insensitive to P. Hence, our results suggest that upon pressurization, hydrophobic plates behave as "soft" surfaces (in the sense of accommodating pressure-dependent changes in water structure) while hydrophilic walls behave as "hard" surfaces.

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.

Thermodynamics, structure, and dynamics of water confined between hydrophobic plates

We perform molecular dynamics simulations of 512 waterlike molecules that interact via the TIP5P potential and are confined between two smooth hydrophobic plates that are separated by 1.10 nm. We find that the anomalous thermodynamic properties of water are shifted to lower temperatures relative to the bulk by approximately 40 K. The dynamics and structure of the confined water resemble bulk water at higher temperatures, consistent with the shift of thermodynamic anomalies to lower temperature. Because of this T shift, our confined water simulations (down to T=220 K) do not reach sufficiently low temperature to observe a liquid-liquid phase transition found for bulk water at T approximately 215 K using the TIP5P potential, but we see inflections in isotherms at lower temperatures presumably due to the presence of a liquid-liquid critical point. We find that the different crystalline structures that can form for two different separations of the plates, 0.7 and 1.10 nm, have no counterparts in the bulk system, and we discuss the relevance to experiments on confined water.