Nucleation of liquid bridges and bubbles in nanoscale capillaries

On the role of history-dependent adsorbate distribution and metastable states in switchable mesoporous metal-organic frameworks

A unique feature of metal-organic frameworks (MOFs) in contrast to rigid nanoporous materials is their structural switchabilty offering a wide range of functionality for sustainable energy storage, separation and sensing applications. This has initiated a series of experimental and theoretical studies predominantly aiming at understanding the thermodynamic conditions to transform and release gas, but the nature of sorption-induced switching transitions remains poorly understood. Here we report experimental evidence for fluid metastability and history-dependent states during sorption triggering the structural change of the framework and leading to the counterintuitive phenomenon of negative gas adsorption (NGA) in flexible MOFs. Preparation of two isoreticular MOFs differing by structural flexibility and performing direct in situ diffusion studies aided by in situ X-ray diffraction, scanning electron microscopy and computational modelling, allowed assessment of n-butane molecular dynamics, phase state, and the framework response to obtain a microscopic picture for each step of the sorption process.

Unveiling the Molecular Origin of Vapor-Liquid Phase Transition of Bulk and Confined Fluids

At temperatures below the critical temperature, discontinuities in the isotherms are one critical issue in the design and construction of separation units, affecting the level of confidence for a prediction of vapor–liquid equilibriums and phase transitions. In this work, we study the molecular mechanisms of fluids that involve the vapor–liquid phase transition in bulk and confinement, utilizing grand canonical (GCE) and meso-canonical (MCE) ensembles of the Monte Carlo simulation. Different geometries of the mesopores, including slit, cylindrical, and spherical, were studied. During phase transitions, condensation/evaporation hysteretic isotherms can be detected by GCE simulation, whereas employing MCE simulation allows us to investigate van der Waals (vdW) loop with a vapor spinodal point, intermediate states, and a liquid spinodal point in the isotherms. Depending on the system, the size of the simulation box, and the MCE method, we are able to identify three distinct groups of vdW-type isotherms for the first time: (1) a smooth S-shaped loop, (2) a stepwise S-shaped loop, and (3) a stepwise S-shaped loop with just a vertical segment. The first isotherm type is noticed in the bulk and pores having small box sizes, in which vapor and liquid phases are close and not clearly identified. The second and the third types occurred in the bulk, cylindrical, and slit mesopores with sufficiently large spaces, where vapor and liquid phases are distinctly separated. Results from our studies provide an insight analysis into vapor–liquid phase transitions, elucidating the effect of the confinement of fluid behaviors in a visual manner.

Pore Size Distribution Analysis Using Developing Hysteresis of Nitrogen in the Cylindrical Pores of Silica

A new method of mesopore size analysis was developed for the cylindrical pores of silica using a developing hysteresis of nitrogen measured at liquid nitrogen temperature. The method is based on the semimacroscopic approach of a modified Bonnet-Wolf model that deals with the grand potential of a vapor bubble in the cylindrical pore. It is capable of assessing correctly the pore structures of mesoporous materials with cylindrical pores using the narrow hysteresis loop of the developing hysteresis. When the mesoporous materials possess cylindrical pores of minor imperfections, two pore size distributions (PSDs) from the adsorption and desorption branches overlap. On the other hand, for cylindrical pores with an enhanced amplitude of pore corrugations and/or constrictions, PSD from the desorption branch is shifted to smaller pore sizes compared to the one from the adsorption branch, and at the same time, both the PSDs evaluated are shifted to lower pore sizes compared to the actual ones. The actual PSD can be assessed from the reversible isotherm measured at a hysteresis critical temperature. In principle, the present method may enable the determination of the PSDs from the adsorption hysteresis measured at any given temperature for mesoporous materials with cylindrical pores.

Machine-Learned Free Energy Surfaces for Capillary Condensation and Evaporation in Mesopores

Using molecular simulations, we study the processes of capillary condensation and capillary evaporation in model mesopores. To determine the phase transition pathway, as well as the corresponding free energy profile, we carry out enhanced sampling molecular simulations using entropy as a reaction coordinate to map the onset of order during the condensation process and of disorder during the evaporation process. The structural analysis shows the role played by intermediate states, characterized by the onset of capillary liquid bridges and bubbles. We also analyze the dependence of the free energy barrier on the pore width. Furthermore, we propose a method to build a machine learning model for the prediction of the free energy surfaces underlying capillary phase transition processes in mesopores.

Revisiting the Nature of Adsorption and Desorption Branches: Temperature Dependence of Adsorption Hysteresis in Ordered Mesoporous Silica

To gain a deeper understanding as to the nature of the adsorption hysteresis due to capillary condensation of nitrogen in ordered mesoporous silicas, we calculated the temperature dependences of the activated condensation, equilibrium transition, and activated desorption pressures for nitrogen in spherical and cylindrical silica pores with several different pore sizes on the basis of semimacroscopic continuum models. The results clearly indicate that the models capture the exact nature of capillary condensation and evaporation phenomena of a fluid in cagelike and cylindrical mesopores. The temperature dependences of the adsorption hysteresis of nitrogen measured confirm previous theoretical predictions for cylindrical pores: for the ordered mesoporous silicas with cylindrical mesopores at least greater than ∼7 nm in diameter, the capillary condensation takes place via a nucleation process followed by a growth process of a bridging meniscus at pressures higher than the equilibrium transition, while the capillary evaporation takes place via a receding meniscus from pore ends at the equilibrium. For SBA-15 and MCM-41 with smaller mesopore sizes, on the other hand, the capillary condensation takes place close to the equilibrium transition pressures, while the capillary evaporation takes place at pressures lower than the equilibrium, owing to single pore blocking due to corrugation of the cylindrical pores. We discuss the effect of curvature on surface tension in capillary condensation, as well as the relation between a change in the mechanisms of adsorption and desorption and the pore corrugation in the cylindrical pores.

Direct Observation of Homogeneous Cavitation in Nanopores

We report on the evaporation of hexane from porous alumina and silicon membranes. These membranes contain billions of independent nanopores tailored to an ink-bottle shape, where a cavity several tens of nanometers in diameter is separated from the bulk vapor by a constriction. For alumina membranes with narrow enough constrictions, we demonstrate that cavity evaporation proceeds by cavitation. Measurements of the pressure dependence of the cavitation rate follow the predictions of the bulk, homogeneous, classical nucleation theory, definitively establishing the relevance of homogeneous cavitation as an evaporation mechanism in mesoporous materials. Our results imply that porous alumina membranes are a promising new system to study liquids in a deeply metastable state.

Nucleation of Capillary Bridges and Bubbles in Nanoconfined CO2

Using molecular simulation, we examine the capillary condensation and the capillary evaporation of CO₂ in cylindrical nanopores. More specifically, we employ the recently developed μV T-S method to determine the microscopic mechanism associated with these processes and the corresponding free energy profiles. We calculate the free energy barrier for capillary condensation and identify that the key step consists in the nucleation of a liquid bridge of a critical size. Similarly, the free energy maximum found for the capillary evaporation process is found to correspond to the nucleation of a vapor bubble of a critical size. In addition, we assess the impact of the strength of the wall-fluid on the height of the free energy barrier and on the critical size of liquid bridges (condensation process) and vapor bubbles (evaporation process). We observe that the height of the free energy barrier increases with the strength of the wall-fluid interactions. Finally, we build a theoretical model, based on capillary theory, to rationalize our findings. In particular, the simulation results reveal a linear scaling of the free energy barrier with the critical size, in excellent agreement with the theoretical predictions.

Stochastic analysis of capillary condensation in disordered mesopores

Most mesoporous materials of practical interest are inherently disordered, which has a significant impact on the condensation and evaporation of vapours in their pores. Traditionally, the effect of disorder is theoretically analyzed in a perturbative approach whereby slight elements of disorder (constriction, corrugation) are added to geometrically ideal pores. We propose an alternative approach, which consists of using a stochastic geometrical model to describe both the porous material and the condensate within the pores. This is done through a multiphase generalisation of the standard Gaussian random field model of disordered materials. The model parameters characterising the condensate provide a low-dimensional approximation of its configuration space, and we use a Derjaguin-Broekhoff-de Boer approximation to calculate the free-energy landscape. Our analysis notably questions the existence of vapour-like metastable states in realistically disordered mesoporous materials. Beyond capillary condensation, our general methodology is applicable to a broad array of confined phenomena.

Critical energy barrier for capillary condensation in mesopores: Hysteresis and reversibility

Capillary condensation in the regime of developing hysteresis occurs at a vapor pressure, Pcond, that is less than that of the vapor-like spinodal. This is because the energy barrier for the vapor-liquid transition from a metastable state at Pcond becomes equal to the energy fluctuation of the system; however, a detailed mechanism of the spontaneous transition has not been acquired even through extensive experimental and simulation studies. We therefore construct accurate atomistic silica mesopore models for MCM-41 and perform molecular simulations (gauge cell Monte Carlo and grand canonical Monte Carlo) for argon adsorption on the models at subcritical temperatures. A careful comparison between the simulation and experiment reveals that the energy barrier for the capillary condensation has a critical dimensionless value, Wc (*) = 0.175, which corresponds to the thermal fluctuation of the system and depends neither on the mesopore size nor on the temperature. We show that the critical energy barrier Wc (*) controls the capillary condensation pressure Pcond and also determines a boundary between the reversible condensation/evaporation regime and the developing hysteresis regime.

Formation of Periodically Arranged Nanobubbles in Mesopores: Capillary Bridge Formation and Cavitation during Sorption and Solidification in an Hierarchical Porous SBA-15 Matrix

We report synchrotron-based small-angle X-ray scattering experiments on a template-grown porous silica matrix (Santa Barbara Amorphous-15) upon in situ sorption of fluorinated pentane C5F12 along with volumetric gas sorption isotherm measurements. Within the mean-field model of Saam and Cole for vapor condensation in cylindrical pores, a nitrogen and C5F12 sorption isotherm is well described by a bimodal pore radius distribution dominated by meso- and micropores with 3.4 and 1.6 nm mean radius, respectively. In the scattering experiments, two different periodicities become evident. One of them (d1 = 11.5 nm) reflects the next nearest neighbor distance in a 2D-hexagonal lattice of tubular mesopores. A second periodicity (d2 = 11.4 nm) found during in situ sorption and freezing experiments is traced back to a superstructure along the cylindrical mesopores. It is compatible with periodic pore corrugations found in electron tomograms of empty SBA-15 by Gommes et al. ( Chem. Mater. 2009, 21, 1311 - 1317). A Rayleigh-Plateau instability occurring at the cylindrical blockcopolymer micelles characteristic of the SBA-15 templating process quantitatively accounts for the superstructure and thus the spatial periodicity of the pore wall corrugation. The consequences of this peculiar morphological feature on the spatial arrangement of C5F12, in particular the formation of periodically arranged nanobubbles (or voids) upon adsorption, desorption, and freezing of liquids, are discussed in terms of capillary bridge formation and cavitation in tubular but periodically corrugated pores.

Influence of system size on the properties of a fluid adsorbed in a nanopore: Physical manifestations and methodological consequences

Hysteresis and discontinuities in the isotherms of a fluid adsorbed in a nanopore in general hamper the determination of equilibrium thermodynamic properties, even in computer simulations. A way around this has been to consider both a reservoir of small size and a pore of small extent in order to restrict the fluctuations of density and approach a classical van der Waals loop. We assess this suggestion by thoroughly studying through Monte Carlo simulations and density functional theory the influence of system size on the equilibrium configurations of the adsorbed fluid and on the resulting isotherms. We stress the importance of pore-symmetry-breaking states that even for modest pore sizes lead to discontinuous isotherms and we discuss the physical relevance of these states and the methodological consequences for computing thermodynamic quantities.

Determination of phase equilibria in confined systems by open pore cell Monte Carlo method

We present a modification of the molecular dynamics simulation method with a unit pore cell with imaginary gas phase [M. Miyahara, T. Yoshioka, and M. Okazaki, J. Chem. Phys. 106, 8124 (1997)] designed for determination of phase equilibria in nanopores. This new method is based on a Monte Carlo technique and it combines the pore cell, opened to the imaginary gas phase (open pore cell), with a gas cell to measure the equilibrium chemical potential of the confined system. The most striking feature of our new method is that the confined system is steadily led to a thermodynamically stable state by forming concave menisci in the open pore cell. This feature of the open pore cell makes it possible to obtain the equilibrium chemical potential with only a single simulation run, unlike existing simulation methods, which need a number of additional runs. We apply the method to evaluate the equilibrium chemical potentials of confined nitrogen in carbon slit pores and silica cylindrical pores at 77 K, and show that the results are in good agreement with those obtained by two conventional thermodynamic integration methods. Moreover, we also show that the proposed method can be particularly useful for determining vapor-liquid and vapor-solid coexistence curves and the triple point of the confined system.

Spherical seed mediated vapor condensation of Lennard-Jones fluid: a density functional theory approach

Vapor to liquid condensation in presence of spherical seed particle of any arbitrary radius ranging from zero to infinity has been investigated using density functional theory, by modeling the local Helmholtz free energy density functional as well as the density profile at the vapor-liquid interface. A general theory is, thus, obtained which provides the different modes of nucleation based on the size of the seed ranging from zero (corresponding to the homogeneous mode of nucleation) to infinity (corresponding to the heterogeneous nucleation on flat surface). The theory is applied to the Lennard-Jones fluid and the optimized shape (i.e., contact angle) and formation free energy of droplets of any arbitrary size have been obtained in this work. The change of the shape (optimized) with the variation of the size of the liquid droplet as well as with the size of the solid substrate has been studied, thus predicting the shape-size relationship in the course of vapor to liquid heterogeneous nucleation on a spherical solid substrate of any particular size. The spinodal decomposition of vapor has also been observed at higher strength of the solid-fluid interaction. The results have been compared with the results of the conventional classical nucleation theory.

Translocation dynamics of freely jointed Lennard-Jones chains into adsorbing pores

Polymer translocation into adsorbing nanopores is studied by using the Fokker-Planck equation of chain diffusion along the energy landscape calculated with Monte Carlo simulations using the incremental gauge cell method. The free energy profile of a translocating chain was found by combining two independent sub-chains, one free but tethered to a hard wall, and the other tethered inside an adsorbing pore. Translocation dynamics were revealed by application of the Fokker-Planck equation for normal diffusion. Adsorption of polymer chains into nanopores involves a competition of attractive adsorption and repulsive steric hindrance contributions to the free energy. Translocation times fell into two regimes depending on the strength of the adsorbing pore. In addition, we found a non-monotonic dependence of translocation times with increasing adsorption strength, with sharp peak associated with local free energy minima along the translocation coordinate.

Adsorption, capillary bridge formation, and cavitation in SBA-15 corrugated mesopores: a Derjaguin-Broekhoff-de Boer analysis

A Derjaguin-Broekhoff-de Boer analysis of adsorption and desorption in SBA-15 mesoporous silica is presented, using realistic geometrical models that account for the pore corrugation in these materials. The model parameters are derived from independent electron tomography and small-angle scattering characterization. A geometrical characteristic of the pore that is found to be important for adsorption is the corrugation length, l(C), which describes the longitudinal size of the geometrical defects along a given pore. Capillary bridges are possible only for large values of l(C). The results are explained in terms of two spinodal and two equilibrium pressures, characterizing the wide and the narrow sections of the pores. Simplified analytical expressions are obtained, which provide necessary conditions for bridge formation and for cavitation in terms of the radii of the narrow and wide sections of the pores, as well as of l(C). Quite generally, the results show that the deviation of the pore shape from that of ideal cylinders is key to understanding adsorption and desorption in corrugated mesopores, notably in SBA-15.

Size dependence of the contact angle of a nanodrop in a nanocavity: density functional theory considerations

The dependence of the contact angles of nanodrops of Lennard-Jones type fluids in nanocavities on their sizes are calculated using a nonlocal density functional theory in a canonical ensemble. Cavities of various radii and depths, various temperatures, as well as various values of the energy parameter of the fluid-solid interactions were considered. It is argued that this dependence might affect strongly, for instance, the rate of heterogeneous nucleation on rough surfaces, which is usually calculated under the assumption of constant contact angle.

Microscopic description of a drop on a solid surface

Two approaches recently suggested for the treatment of macro- or nanodrops on smooth or rough, planar or curved, solid surfaces, based on fluid-fluid and fluid-solid interaction potentials are reviewed. The first one employs the minimization of the total potential energy of a drop by assuming that the drop has a well defined profile and a constant liquid density in its entire volume with the exception of the monolayer nearest to the surface where the density has a different value. As a result, a differential equation for the drop profile as well as the necessary boundary conditions are derived which involve the parameters of the interaction potentials and do not contain such macroscopic characteristics as the surface tensions. As a consequence, the macroscopic and microscopic contact angles which the drop profile makes with the surface can be calculated. The macroscopic angle is obtained via the extrapolation of the circular part of the drop profile valid at some distance from the surface up to the solid surface. The microscopic angle is formed at the intersection of the real profile (which is not circular near the surface) with the surface. The theory provides a relation between these two angles. The ranges of the microscopic parameters of the interaction potentials for which (i) the drop can have any height (volume), (ii) the drop can have a restricted height but unrestricted volume, and (iii) a drop cannot be formed on the surface were identified. The theory was also extended to the description of a drop on a rough surface. The second approach is based on a nonlocal density functional theory (DFT), which accounts for the inhomogeneity of the liquid density and temperature effects, features which are missing in the first approach. Although the computational difficulties restrict its application to drops of only several nanometers, the theory can be applied indirectly to macrodrops by calculating the surface tensions and using the Young equation to determine the contact angle. Employing the canonical ensemble version of the DFT, nanodrops on smooth and rough solid surfaces could be investigated and their characteristics, such as the drop profile, contact angle, as well as the fluid density distribution inside the drop can be determined as functions of the parameters of the interaction potentials and temperature. It was found that the contact angle of the drop has a simple (quasi)universal dependence on the energy parameter epsilon(fs) of the fluid-solid interaction potential and temperature. The main feature of this dependence is the existence of a fixed value theta(0) of the contact angle theta which separates the solid substrates (characterized by the energy parameter epsilon(fs) of the fluid-solid interaction potential) into two classes with respect to their temperature dependence. For theta>theta(0) the contact angle monotonously increases and for theta<theta(0) monotonously decreases with increasing temperature. For theta=theta(0) the contact angle is independent of temperature. The results obtained via DFT were also applied to check the validity of the macroscopic phenomenological equations (Cassie-Baxter and Wenzel equations) for drops on rough surfaces, and of the equation for the sticking force of a drop on an inclined surface.

Symmetry breaking in confined fluids

The recent progress in the theoretical investigation of the symmetry breaking (the existence of a stable state of a system, in which the symmetry is lower than the symmetry of the system itself) for classical and quantum fluids is reviewed. The emphasis is on the conditions which cause symmetry breaking in the density distribution for one component fluids and binary mixtures confined in a closed nanoslit between identical solid walls. The existing studies have revealed that two kinds of symmetry breaking can occur in such systems. First, a one-dimensional symmetry breaking occurs only in the direction normal to the walls as a fluid density profile asymmetric with respect of the middle of the slit and uniform in any direction parallel to the walls. Second, a two-dimensional symmetry breaking occurs in the fluid density distribution which is nonuniform in one of the directions parallel to the walls and asymmetrical in the direction normal to the walls. It manifests through liquid bumps and bridges in the fluid density distribution. For one component fluids, conditions of existence of symmetry breaking are provided in terms of the average fluid density, strength of fluid-solid interactions, distance at which the solid wall generates a hard core repulsion, and temperature. In the case of binary mixtures, the occurrence of symmetry breaking also depends on the composition of the confined mixtures.

Tuning colloidal interactions in subcritical solvents by solvophobicity: explicit versus implicit modeling

The distance-resolved effective interaction between two colloidal particles in a subcritical solvent is explored both by an explicit and implicit modeling. An implicit solvent approach based on a simple thermodynamic interface model is tested against grand-canonical Monte Carlo computer simulations using explicit Lennard-Jones solvent molecules. Close to liquid-gas coexistence, a joint gas bubble surrounding the colloidal particle pair yields an effective attraction between the colloidal particles, the strength of which can be vastly tuned by the solvophobicity of the colloids. The implicit model is in good agreement with our explicit computer simulations, thus enabling an efficient modeling and evaluation of colloidal interactions and self-assembly in subcritical solvent environments.

Nucleation and hysteresis of vapor-liquid phase transitions in confined spaces: effects of fluid-wall interaction

In this work, we propose a method to stabilize a nucleus in the framework of lattice density-functional theory (LDFT) by imposing a suitable constraint. Using this method, the shape of critical nucleus and height of the nucleation barrier can be determined without using a predefined nucleus as input. As an application of this method, we study the nucleation behavior of vapor-liquid transition in nanosquare pores with infinite length and relate the observed hysteresis loop on an adsorption isotherm to the nucleation mechanism. According to the dependence of hysteresis and the nucleation mechanism on the fluid-wall interaction, w , in this work, we have classified w into three regions ( w>0.9 , 0.1< or =w< or =0.9 , and w<0.1 ), which are denoted as strongly, moderately, and weakly attractive fluid-wall interaction, respectively. The dependence of hysteresis on the fluid-wall interaction is interpreted by the different nucleation mechanisms. Our constrained LDFT calculations also show that the different transition paths may induce different nucleation behaviors. The transition path dependence should be considered if morphological transition of nuclei exists during a nucleation process.

Nanodrop on a nanorough solid surface: density functional theory considerations

The density distributions and contact angles of liquid nanodrops on nanorough solid surfaces are determined on the basis of a nonlocal density functional theory. Two kinds of roughness, chemical and physical, are examined. The former considers the substrate as a sequence of two kinds of semi-infinite vertical plates of equal thicknesses but of different natures with different strengths for the liquid-solid interactions. The physical roughness involves an ordered set of pillars on a flat homogeneous surface. Both hydrophobic and hydrophilic surfaces were considered. For the chemical roughness, the contact angle which the drop makes with the flat surface increases when the strength of the liquid-solid interaction for one kind of plates decreases with respect to the fixed value of the other kind of plates. Such a behavior is in agreement with the Cassie-Baxter expression derived from macroscopic considerations. For the physical roughness on a hydrophobic surface, the contact angle which a drop makes with the plane containing the tops of the pillars increases with increasing roughness. Such a behavior is consistent with the Wenzel formula developed for macroscopic drops. For hydrophilic surfaces, as the roughness increases the contact angle first increases, in contradiction with the Wenzel formula, which predicts for hydrophilic surfaces a decrease of the contact angle with increasing roughness. However, a further increase in roughness changes nonmonotonously the contact angle, and at some roughness, the drop disappears and only a liquid film is present on the surface. It was also found that the contact angle has a periodic dependence on the volume of the drop.

Macrophase and microphase separations for surfactants adsorbed on solid surfaces: a gauge cell monte carlo study in the lattice model

By combining the gauge cell method and lattice model, we study the surface phase transition and adsorption behaviors of surfactants on a solid surface. Two different cases are considered in this work: macrophase transition and adsorption in a single-phase region. For the case of macrophase transition, where two phases coexist, we investigate the shape and size of the critical nuclei and determine the height of the nucleation barrier. It is found that the nucleation depends on the bulk surfactant concentration. Our simulations show that there exist a critical temperature and critical adsorption energy, below which the transition from low-affinity adsorption to the bilayer structure shows the characteristic of a typical first-order phase transition. Such a surface phase transition in the adsorption isotherm is featured by a hysteresis loop. The hysteresis loop becomes narrower at higher temperature and weaker adsorption energy and finally disappears at the critical value. For the case where no macrophase transition occurs, we study the adsorption isotherm and microphase separation in a single-phase region. The simulation results indicate that the adsorption isotherm in adsorption processes is divided into four regions in a log-log plot, being in agreement with experimental observations. In this work, the four regions are called the low-affinity adsorption region, the hemimicelle region, the morphological transition region, and the plateau region. Simulation results reveal that in the second region the adsorbed monomers aggregate and nucleate hemimicelles, while adsorption in the third region is accompanied by morphological transitions.

Symmetry breaking in binary mixtures in closed nanoslits

The symmetry breaking (SB) of the fluid density distribution (FDD) in closed nanoslits between two identical parallel solid walls described by Berim and Ruckenstein [J. Chem. Phys. 128, 024704 (2008)] for a single component fluid is examined for binary mixtures on the basis of a nonlocal canonical ensemble density functional theory. As in Monte Carlo simulations, the periodicity of the FDD in one of the lateral (parallel to the wall surfaces) directions, denoted as the x direction, was assumed. In the other lateral direction, y direction, the FDD was considered to be uniform. The molecules of the two components have different diameters and their Lennard-Jones interaction potentials have different energy parameters. It was found that depending on the average fluid density in the slit and mixture composition, SB can occur for both or none of the components but never for only one of them. In the direction perpendicular to the walls (h direction), the FDDs of both components can be asymmetrical about the middle plane between walls. In the x direction, the SB occurs as bumps and bridges enriched in one of the components, whereas the composition of the mixture between them is enriched in the other component. The dependence of the SB states on the length Lx of the FDD period at fixed average densities of the two components was examined for Lx in the range from 10 to 120 molecular diameters of the smaller size component. It was shown that for large Lx, the stable state of the system corresponds to a bridge. Because the free energy of that state decreases monotonically with increasing Lx, one can conclude that the real period is very large (infinite) and that a single bridge exists in the slit.

Two-dimensional symmetry breaking of fluid density distribution in closed nanoslits

Stable and metastable fluid density distributions (FDDs) in a closed nanoslit between two identical parallel solid walls have been identified on the basis of a nonlocal canonical ensemble density functional theory. Similar to Monte Carlo simulations, periodicity of the FDD in one of the lateral (parallel to the walls surfaces) directions, denoted as the x direction, was assumed. In the other lateral direction, y direction, the FDD was considered uniform. It was found that depending on the average fluid density in the slit, both uniform as well as nonuniform FDDs in the x direction can occur. The uniform FDDs are either symmetric or asymmetric about the middle plane between walls; the latter FDD being the consequence of a symmetry breaking across the slit. The nonuniform FDDs in the x direction occur either in the form of a bump on a thin liquid film covering the walls or as a liquid bridge between those walls and provide symmetry breaking in the x direction. For small and large average densities, the stable state is uniform in the x direction and is symmetric about the middle plane between walls. In the intermediate range of the average density and depending on the length L(x) of the FDD period, the stable state can be represented either by a FDD, which is uniform in the x direction and asymmetric about the middle of the slit (small values of L(x)), or by a bump- and bridgelike FDD for intermediate and large values of L(x), respectively. These results are in agreement with the Monte Carlo simulations performed earlier by other authors. Because the free energy of the stable state decreases monotonically with increasing L(x), one can conclude that the real period is very large (infinite) and that for the values of the parameters employed, a single bridge of finite length over the entire slit is generated.

Effect of fluid-solid interactions on symmetry breaking in closed nanoslits

The possibility of symmetry breaking of the fluid (argon) density distribution across a long closed slit with identical walls composed of solid carbon dioxide was noted in previous papers by the authors. The main conclusion was that there is a range of average densities in which symmetry breaking occurs and that outside that range the fluid density profile is symmetrical. A critical temperature T(sb) was also identified below which symmetry breaking can occur. In this paper, symmetry breaking is examined for walls made of other materials and it is shown that it occurs only when the energy parameter epsilon(fw) of the fluid-wall interaction in the Lennard-Jones potential satisfies the inequalities epsilon(fw1) < or = epsilon(fw) < or = epsilon(fw2), where epsilon(fw1) and epsilon(fw2) are temperature-dependent critical values of epsilon(fw). The value of epsilon(fw1) increases and that of epsilon(fw2) decreases with increasing temperature. The comparison of the theory with Monte Carlo simulations confirms the existence of symmetry breaking across the slit. The possibility of symmetry breaking along the slit is also noted.

Lower closure point of adsorption hysteresis in ordered mesoporous silicas

To examine the nature of the lower closure point of adsorption hysteresis in ordered mesoporous silicas, we measured the temperature dependence of the adsorption-desorption isotherm of nitrogen for three kinds of ordered silicas with cagelike pores and three kinds of ordered silicas with cylindrical pores. The lower closure point pressure of nitrogen in the cagelike pores with sufficiently small necks, that is, the cavitation pressure of a confined liquid, did not depend appreciably on the cage size in the temperature region far away from a hysteresis critical temperature (Tch) but its cage-size dependence was noticeable in the vicinity of Tch. The lower closure point in the cylindrical pores depended on the pore size, and its thermal behavior was totally different from that in the cagelike pores. Nevertheless, the hysteresis critical points of nitrogen in the ordered mesoporous silicas, which are defined as a threshold of temperatures (Tch) and pressure above which reversible capillary condensation takes place in a given size and shape of pores, fell on a common line in a temperature-pressure diagram regardless of the pore geometries. We consider this finding as evidence that capillary evaporation in the cylindrical pores follows a cavitation process in the vicinity of Tch in the same way as that in the cagelike pores and also that the low limit of the hysteresis loop that has been long recognized since 1965 is due to the occurrence of a vapor bubble in a stretched metastable liquid confined to the pores with decreasing pressure (cavitation).