Phase behavior and local structure of a binary mixture in pores: mean-field lattice model calculations for analyzing neutron scattering data

We investigate the phase behavior of an asymmetric binary liquid A-W mixture confined between two planar homogenous substrates (slit pore). Molecules of species W interact preferentially with the solid walls via a long-range potential. Assuming nearest-neighbor attractions between the liquid molecules, we employ a lattice-gas model and a mean-field approximation for the grand potential. Minimization of this potential yields the density profiles of thermodynamically stable phases for fixed temperature, chemical potentials of both species, pore width and strengths of attraction. This model is used to analyze experimental small-angle neutron-scattering (SANS) data on the microscopic structure of the binary system isobutyric acid (iBA)+heavy water (D2O) inside a mesoscopic porous matrix (controlled-pore glass of about 10 nm mean pore width). Confinement-independent model parameters are adjusted so that the theoretical liquid-liquid coexistence curve in the bulk matches its experimental counterpart. By choosing appropriate values of the pore width and the attraction strength between substrates and water we analyze the effect of confinement on the phase diagram. In addition to a depression of the liquid-liquid critical point we observe surface induced phase transitions as well as water-film adsorption near the walls. The temperature dependence of the structure of water-rich and iBA-rich phases of constant composition are discussed in detail. The theoretical predictions are consistent with results of the SANS study and assist their interpretation.     

Density measurement of 1-d confined water by small angle neutron scattering method: pore size and hydration level dependences

Small angle neutron scattering (SANS) is used to measure the absolute density of water contained in 1-D cylindrical pores of a silica material MCM-41-S with pore diameters of 19 and 15 A. By being able to suppress the homogeneous nucleation process inside the narrow pore, one can keep water in the liquid state down to at least 160 K. From a combined analysis of SANS data from both H(2)O and D(2)O hydrated samples, we determined the absolute value of the density of 1-D confined water. We found that the average density of water inside the fully hydrated 19 A pore is 8% higher than that of the bulk water at room temperature. The temperature derivative of the density shows a pronounced peak at T(L) = 235 K signaling the crossing of the Widom line at ambient pressure and confirming the existence of a liquid-liquid phase transition at an elevated pressure. Pore size and hydration level dependences of the density are also studied.     

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.     

Water in nanopores. I. Coexistence curves from Gibbs ensemble Monte Carlo simulations

Coexistence curves of water in cylindrical and slitlike nanopores of different size and water-substrate interaction strength were simulated in the Gibbs ensemble. The two-phase coexistence regions cover a wide range of pore filling level and temperature, including ambient temperature. Five different kinds of two-phase coexistence are observed. A single liquid-vapor coexistence is observed in hydrophobic and moderately hydrophilic pores. Surface transitions split from the main liquid-vapor coexistence region, when the water-substrate interaction becomes comparable or stronger than the water-water pair interaction. In this case prewetting, one and two layering transitions were observed. The critical temperature of the first layering transition decreases with strengthening water-substrate interaction towards the critical temperature expected for two-dimensional systems and is not sensitive to the variation of pore size and shape. Liquid-vapor phase transition in a pore with a wall which is already covered with two water layers is most typical for hydrophilic pores. The critical temperature of this transition is very sensitive to the pore size, in contrast to the liquid-vapor critical temperature in hydrophobic pores. The observed rich phase behavior of water in pores evidences that the knowledge of coexistence curves is of crucial importance for the analysis of experimental results and a prerequiste of meaningful simulations.     

Thermodynamic and dielectric studies concerning the influence of cylindrical submicrometer confinement on heptyloxycyanobiphenyl

Measurements of the specific heat and the static dielectric permittivity of heptyloxycyanobiphenyl (7OCB) confined to the 0.2 microm diameter parallel cylindrical pores of Anopore membranes in the isotropic phase and nematic mesophase, are presented. A comparison between the bulk and the confined 7OCB in treated and untreated pore wall surfaces using a chemical surfactant (HTBA) is performed. Both the treated and untreated membrane confinements seem to affect the nematic-to-isotropic phase transition by a downshift in transition temperature and some rounding at the specific-heat maximum, in a way similar to that which was earlier published for other liquid crystals confined in the same geometry. The static dielectric measurements clearly point out that untreated membrane confinement is axial, with the nematic director aligned parallel to the pore axis being homeotropic bulklike, i.e., with the nematic director aligned perpendicular to the electrode cell surfaces. After chemical surfactant treatment, the nematic director is constrained in a radial alignment being perpendicular to the pore walls. The dielectric measurements are revealed to be specially sensible to analyze the surface-induced nematic order due to the pore wall. The tricritical nature of the nematic-to-isotropic phase transition in bulk 7OCB as well as in treated and untreated Anopore confined geometries is discussed through both the specific heat and the static dielectric data.     

Transitions between disordered phases in supercooled liquid silicon

We have investigated the transitions between disordered phases in supercooled liquid silicon using computer simulations. The thermodynamic properties were directly obtained from the free energy, which was computed using the recently proposed reversible scaling method. The calculated free energies of the crystalline and liquid phases of silicon at zero pressure, obtained using the environment dependent interatomic potential, are in excellent agreement with the available experimental data. The results show that, at zero pressure, a weak first-order liquid-liquid transition occurs at 1135 K and a continuous liquid-amorphous transition takes place at 843 K. These results are consistent with the existence of a second critical point for the liquid-liquid transition at a negative pressure.     

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.     

Liquid-liquid phase transitions in supercooled water studied by computer simulations of various water models

Liquid-liquid and liquid-vapor coexistence regions of various water models were determined by Monte Carlo (MC) simulations of isotherms of density fluctuation-restricted systems and by Gibbs ensemble MC simulations. All studied water models show multiple liquid-liquid phase transitions in the supercooled region: we observe two transitions of the TIP4P, TIP5P, and SPCE models and three transitions of the ST2 model. The location of these phase transitions with respect to the liquid-vapor coexistence curve and the glass temperature is highly sensitive to the water model and its implementation. We suggest that the apparent thermodynamic singularity of real liquid water in the supercooled region at about 228 K is caused by an approach to the spinodal of the first (lowest density) liquid-liquid phase transition. The well-known density maximum of liquid water at 277 K is related to the second liquid-liquid phase transition, which is located at positive pressures with a critical point close to the maximum. A possible order parameter and the universality class of liquid-liquid phase transitions in one-component fluids are discussed.     

Effect of temperature on the adsorption of water in porous carbons

We report experimental and simulation studies to investigate the effect of temperature on the adsorption isotherms for water in carbons. Adsorption isotherms are measured by a gravimetric technique in carbon-fiber monoliths at 378 and 423 K and studied by molecular simulation in ideal carbon pores in the temperature range 298-600 K. Experimental adsorption isotherms show a gradual water uptake, as the pressure increases, and narrow adsorption-desorption hysteresis loops. In contrast, simulated adsorption isotherms at room temperature are characterized by negligible uptake at low pressures, sudden and complete pore filling once a threshold pressure is reached, and wide adsorption-desorption hysteresis loops. As the temperature increases, the relative pressure at which pore filling occurs increases and the size of the hysteresis loop decreases. Experimental adsorption-desorption hysteresis loops are narrower than those from simulation. Discrepancies between simulation and experimental results are attributed to heterogeneities in chemical composition, pore connectivity, and nonuniform pore-size distribution, which are not accounted for in the simulation model. The hysteresis phase diagram for confined water is obtained by recording the pressure-density conditions that bound the simulated hysteresis loop at each temperature. We find that the hysteresis critical temperature, i.e., the lowest temperature at which no hysteresis is detected, can be hundreds of degrees lower than the vapor-liquid critical temperature for bulk model water. The properties of confined water are discussed with the aid of simulation snapshots and by analyzing the structure of the confined fluid.     

Pressure effects on the transitions between disordered phases in supercooled liquid silicon

We investigate the pressure effects on the transitions between the disordered phases in supercooled liquid silicon through Monte Carlo simulations and efficient methods to compute free energies. Our calculations, using an environment dependent interatomic potential for Si, indicate that at zero pressure the liquid-liquid phase transition, between the high density liquid and the low density liquid, occurs at a temperature 325K below melting. We found that the liquid-liquid transition temperature decreases with increasing pressure, following the liquid-solid coexistence curve. As pressure increases, the liquid-liquid coexistence curve approaches the region where the glass transition between the low density liquid and the low density amorphous takes place. Above 5 GPa, our calculations show that the liquid-liquid transition is suppressed by the glassy dynamics of the system. We also found that above 5 GPa, the glass transition temperature is lower than that at lower pressures, suggesting that under these conditions the glass transition occurs between the high density liquid and the high density amorphous.     

Local structure of a phase-separating binary mixture in a mesoporous glass matrix studied by small-angle neutron scattering

The mesoscopic structure of the binary system isobutyric acid + heavy water (D(2)O) confined in a porous glass (controlled-pore silica glass, mean pore width ca. 10 nm) was studied by small-angle neutron scattering at off-critical compositions in a temperature range above and below the upper critical solution point. The scattering data were analyzed in terms of a structure factor model similar to that proposed by Formisano and Teixeira [Eur. Phys. J. E 1, 1 (2000)], but allowing for both Ornstein-Zernike-type composition fluctuations and domainlike structures in the microphase-separated state of the pore liquid. The results indicate that the phase separation in the pores is shifted by ca. 10 K and spread out in temperature. Microphase separation is pictured as a transition from partial segregation at high temperature, due to the strong preferential adsorption of water at the pore wall, to a tube or capsule configuration of the two phases at low temperatures, depending on the overall composition of the pore liquid. Results for samples in which the composition of the pore liquid can vary with temperature due to equilibration with extra-pore liquid are consistent with this picture.     

Confinement Effects on the Magnetic Ionic Liquid 1-Ethyl-3-methylimidazolium Tetrachloroferrate(III)

Confinement effects for the magnetoresponsive ionic liquid 1-ethyl-3-methylimidazolium tetrachloroferrate(III), [C₂mim]FeCl₄, are explored from thermal, spectroscopic, and magnetic points of view. Placing the ionic liquid inside SBA-15 mesoporous silica produces a significant impact on the material’s response to temperature, pressure, and magnetic fields. Isobaric thermal experiments show melting point reductions that depend on the pore diameter of the mesopores. The confinement-induced reductions in phase transition temperature follow the Gibbs–Thomson equation if a 1.60 nm non-freezable interfacial layer is postulated to exist along the pore wall. Isothermal pressure-dependent infrared spectroscopy reveals a similar modification to phase transition pressures, with the confined ionic liquid requiring higher pressures to trigger phase transformation than the unconfined system. Confinement also impedes ion transport as activation energies are elevated when the ionic liquid is placed inside the mesopores. Finally, the antiferromagnetic ordering that characterizes unconfined [C₂mim]FeCl₄ is suppressed when the ionic liquid is confined in 5.39-nm pores. Thus, confinement provides another avenue for manipulating the magnetic properties of this compound.     

Molecular modeling of freezing of simple fluids confined within carbon nanotubes

We report Monte Carlo simulation results for freezing of Lennard-Jones carbon tetrachloride confined within model multiwalled carbon nanotubes of different diameters. The structure and thermodynamic stability of the confined phases, as well as the transition temperatures, were determined from parallel tempering grand canonical Monte Carlo simulations and free-energy calculations. The simulations show that the adsorbate forms concentric molecular layers that solidify into defective quasi-two-dimensional hexagonal crystals. Freezing in such concentric layers occurs via intermediate phases that show remnants of hexatic behavior, similar to the freezing mechanism observed for slit pores in previous works. The adsorbate molecules in the inner regions of the pore also exhibit changes in their properties upon reduction of temperature. The structural changes in the different regions of adsorbate occur at temperatures above or below the bulk freezing point, depending on pore diameter and distance of the adsorbate molecules from the pore wall. The simulations show evidence of a rich phase behavior in confinement; a number of phases, some of them inhomogeneous, were observed for the pore sizes considered. The multiple transition temperatures obtained from the simulations were found to be in good agreement with recent dielectric relaxation spectroscopy experiments for CCl(4) confined within multiwalled carbon nanotubes.     

Effect of Nanoscale Confinement on Freezing of Modified Water at Room Temperature and Ambient Pressure

Understanding the phase behavior of confined water is central to fields as diverse as heterogeneous catalysis, corrosion, nanofluidics, and to emerging energy technologies. Altering the state points (temperature, pressure, etc.) or introduction of a foreign surface can result in the phase transformation of water. At room temperature, ice nucleation is a very rare event and extremely high pressures in the GPa–TPa range are required to freeze water. Here, we perform computer experiments to artificially alter the balance between electrostatic and dispersion interactions between water molecules, and demonstrate nucleation and growth of ice at room temperature in a nanoconfined environment. Local perturbations in dispersive and electrostatic interactions near the surface are shown to provide the seed for nucleation (nucleation sites), which lead to room temperature liquid–solid phase transition of confined water. Crystallization of water occurs over several tens of nanometers and is shown to be independent of the nature of the substrate (hydrophilic oxide vs. hydrophobic graphene and crystalline oxide vs. amorphous diamond‐like carbon). Our results lead us to hypothesize that the freezing transition of confined water can be controlled by tuning the relative dispersive and electrostatic interaction.     

An experimental study of melting of CCl4 in carbon nanotubes

We report dielectric relaxation spectroscopy measurements of the melting point of carbon tetrachloride confined within open-tip multi-walled carbon nanotubes with two different pore diameters, 4.0 and 2.8 nm. In both cases, a single transition temperature well above the bulk melting point was obtained for confined CCl4. These results contrast with what was obtained in our previous measurements using carbon nanotubes with a pore diameter of 5.0 nm, where multiple transition temperatures both above and below the bulk melting point of CCl4 were observed. Our experimental measurements are consistent with our recent molecular simulation results (F. R. Hung, B. Coasne, E. E. Santiso, K. E. Gubbins, F. R. Siperstein and M. Sliwinska-Bartkowiak, J. Chem. Phys., 2005, 122, 144706). Although the simulations overestimate the temperatures in which melting upon confinement occurs, both simulations and experiments suggest that all regions of adsorbate freeze at the same temperature, and that freezing occurs at higher temperatures upon reduction of the pore diameter.     

Direct measurements of pore fluid density by vibrating tube densimetry

The densities of pore-confined fluids were measured for the first time by means of vibrating tube densimetry (VTD). A custom-built high-pressure, high-temperature vibrating tube densimeter was used to measure the densities of propane at subcritical and supercritical temperatures (between 35 and 97 °C) and carbon dioxide at supercritical temperatures (between 32 and 50 °C) saturating hydrophobic silica aerogel (0.2 g/cm(3), 90% porosity) synthesized inside Hastelloy U-tubes. Additionally, supercritical isotherms of excess adsorption for CO(2) and the same porous solid were measured gravimetrically using a precise magnetically coupled microbalance. Pore fluid densities and total adsorption isotherms increased monotonically with increasing density of the bulk fluid, in contrast to excess adsorption isotherms, which reached a maximum and then decreased toward zero or negative values above the critical density of the bulk fluid. The isotherms of confined fluid density and excess adsorption obtained by VTD contain additional information. For instance, the maxima of excess adsorption occur below the critical density of the bulk fluid at the beginning of the plateau region in the total adsorption, marking the end of the transition of pore fluid to a denser, liquidlike pore phase. Compression of the confined fluid significantly beyond the density of the bulk fluid at the same temperature was observed even at subcritical temperatures. The effect of pore confinement on the liquid-vapor critical temperature of propane was less than ~1.7 K. The results for propane and carbon dioxide showed similarity in the sense of the principle of corresponding states. Good quantitative agreement was obtained between excess adsorption isotherms determined from VTD total adsorption results and those measured gravimetrically at the same temperature, confirming the validity of the vibrating tube measurements. Thus, it is demonstrated that vibrating tube densimetry is a novel experimental approach capable of providing directly the average density of pore-confined fluids, and hence complementary to the conventional gravimetric or volumetric/piezometric adsorption techniques, which yield the excess adsorption (the Gibbsian surface excess).     

2H-solid state NMR and DSC study of isobutyric acid in mesoporous silica materials

Solid state deuterium NMR has been used to study the molecular motion of d(6)-isobutyric acid (d(6)-iBA) in the pure (unconfined) state and confined in the cylindrical pores of two periodic mesoporous silica materials (MCM-41, pore size 3.3 nm and SBA-15, pore size 8 nm), and in a controlled pore glass (CPG-10-75, pore size ca. 10 nm). The line shape analysis of the spectra at different temperatures revealed three rotational states of the iBA molecules: liquid (fast anisotropic reorientation of the molecule), solid I (rotation of the methyl group) and solid II (no rotational motion on the time scale of the experiment). Transition temperatures between these states were determined from the temperature dependence of the fraction of molecules in these states. Whereas the solid I-solid II transition temperature is not affected by confinement, a significant lowering of the liquid-solid I transition temperature in the pores relative to the bulk acid was found for the three matrix materials, exhibiting an unusual dependence on pore size and pore morphology. Complementary DSC measurements on the same systems show that the rotational melting (solid I-liquid) of d(6)-iBA in the pores occurs at a temperature 20-45 K below the thermodynamic melting point. This finding indicated that the decoupling of rotational and translational degrees of freedom in phase transitions in confined systems previously found for benzene is not restricted to molecules with non-specific interactions, but represents a more general phenomenon.     

Phase diagram of supercooled water confined to hydrophilic nanopores

We present a phase diagram for water confined to cylindrical silica nanopores in terms of pressure, temperature, and pore radius. The confining cylindrical wall is hydrophilic and disordered, which has a destabilizing effect on ordered water structure. The phase diagram for this class of systems is derived from general arguments, with parameters taken from experimental observations and computer simulations and with assumptions tested by computer simulation. Phase space divides into three regions: a single liquid, a crystal-like solid, and glass. For large pores, radii exceeding 1 nm, water exhibits liquid and crystal-like behaviors, with abrupt crossovers between these regimes. For small pore radii, crystal-like behavior is unstable and water remains amorphous for all non-zero temperatures. At low enough temperatures, these states are glasses. Several experimental results for supercooled water can be understood in terms of the phase diagram we present.     

Confinement induced critical micelle concentration shift

In this work, extensive lattice Monte Carlo simulations were performed to investigate the influence of confinement on critical micelle concentration (CMC). It is found that the CMC of surfactants in a confined space is shifted from its bulk value, and the shift is affected by the presence of the confining boundaries, which induces both the finite size effect and the wall-surfactant interaction. In general, for strongly confined system (the system with narrow pore size), the finite size effect dominates the CMC shift because the confined space cannot accommodate fully developed micelles, and the rapid increase of the entropic loss due to the decrease of the pore size results in the rapid increase of CMC. In contrast, for a weakly confined space, the CMC shift depends on the interaction between the walls and surfactants. For the systems with two weakly hydrophilic surfaces, the local density depletion of the surfactants near the walls results in lower CMCs than the bulk value, and the CMC shifts to a higher value as the pore size increases. For the systems with moderately hydrophilic surfaces, the shifts of CMCs show a similar behavior as those for weakly hydrophilic surfaces, but the CMCs are near their bulk values in the range of weak confinement. For the systems with strongly attractive wall-surfactant interactions, the strong adsorption also results in lower CMCs than their bulk value, but the CMCs decrease with the increase of pore size.     

Confinement effect on phase transitions of a discotic liquid crystal: a calorimetric study

Differential scanning calorimetry was used to investigate the confinement effects on the phase transition behaviour of a discotic liquid crystal. The liquid crystal studied is the hexa-n-octanoate of rufigallol (RHO); Millipore membranes of various pore sizes were the confining materials. The polymorphism of RHO is affected by confinement. The transition from an enantiotropic columnar phase (D1) to a monotropic columnar phase (D2) is supressed in membranes with pore sizes 500 A. The transformation from D1 to the crystalline phase is also perturbed, particularly in the membrane having an average pore size of 250 A. In the first case the crystal formed displays a double-melting endotherm, with a distinct structure melting at lower temperatures; in the other, the induction period of isothermal crystallization becomes longer and the global rate of crystallization is slowed. However, confinement shows no effect on the overall crystallization mechanism; a similar Avrami constant of n ~ 3 was obtained for both confined and bulk RHO. An analysis of the results is presented.