Universal fluctuation-induced corrections to the Kelvin equation for capillary condensation

Breaking Cassie's Law for Condensation in a Nanopatterned Slit

We study the phase transitions of a fluid confined in a capillary slit made from two adjacent walls, each of which are a periodic composite of stripes of two different materials. For wide slits the capillary condensation occurs at a pressure which is described accurately by a combination of the Kelvin equation and the Cassie law for an averaged contact angle. However, for narrow slits the condensation occurs in two steps involving an intermediate bridging phase, with the corresponding pressures described by two new Kelvin equations. These are characterised by different contact angles due to interfacial pinning, with one larger and one smaller than the Cassie angle. We determine the triple point and predict two types of dispersion force induced Derjaguin-like corrections due to mesoscopic volume reduction and the singular free-energy contribution from nanodroplets and bubbles. We test these predictions using a fully microscopic density functional model which confirms their validity even for molecularly narrow slits. Analogous mesoscopic corrections are also predicted for two-dimensional systems arising from thermally induced interfacial wandering.

Phase diagrams and solvation forces of a uniaxial ferromagnet in a slit--the double-parabola approach

We obtain the phase diagrams of a uniaxial ferromagnet enclosed in a slit and analyse the properties of the solvation forces acting between the slit walls. Our analysis is based on the mean-field theory applied to the double-parabola Hamiltonian for a uniaxial ferromagnet in a slit, which allows us to derive several results in analytic form. In particular, the Kelvin law and the law describing the influence of the slit width on the position of the delocalization transition as compared to the wetting transition are derived and the corresponding critical exponent is evaluated. The influence of the thermodynamic state of the system on the properties of the solvation forces, in particular their discontinuous behaviour at coexistence lines, is discussed.

The order of condensation in capillary grooves

We consider capillary condensation in a deep groove of width L. The transition occurs at a pressure p(co)(L) described, for large widths, by the Kelvin equation p(sat) - p(co)(L) = 2σ cosθ/L, where θ is the contact angle at the side walls and σ is the surface tension. The order of the transition is determined by the contact angle of the capped end θcap; it is continuous if the liquid completely wets the cap, and first-order otherwise. When the transition is first-order, corner menisci at the bottom of the capillary lead to a pronounced metastability, determined by a complementary Kelvin equation Δp(L) = 2σ sinθcap/L. On approaching the wetting temperature of the capillary cap, the corner menisci merge and a single meniscus unbinds from the bottom of the groove. Finite-size scaling shifts, crossover behaviour and critical singularities are determined at mean-field level and beyond. Numerical and experimental results showing the continuous nature of condensation for θcap = 0 and the influence of corner menisci on adsorption isotherms are presented.

Capillary condensation in a square geometry with surface fields

We study the influence of wetting on capillary condensation for a simple fluid in a square geometry with surface fields, where the reference system is an infinitely long slit. The corner transfer matrix renormalization group method has been extended to study a two-dimensional Ising model confined in an L × L geometry with equal surface fields. Our results have confirmed that in both geometries the coexistence line shift is governed by the same scaling powers, but their prefactors are different.

Finite-size effects on the behavior of the susceptibility in van der Waals films bounded by strongly absorbing substrates

We study critical point finite-size effects in the case of the susceptibility of a film in which interactions are characterized by a van der Waals-type power law tail. The geometry is appropriate to a slablike system with two bounding surfaces. Boundary conditions are consistent with surfaces that both prefer the same phase in the low temperature, or broken symmetry, state. We take into account both interactions within the system and interactions between the constituents of the system and the material surrounding it. Specific predictions are made with respect to the behavior of 3He and 4He films in the vicinity of their respective liquid-vapor critical points.

Capillary evaporation in pores

We combine a density functional theory (DFT) treatment of capillary evaporation in a cylindrical pore with the morphometric approach in order to study the formation and breaking of bubbles in a hydrophobically lined part of a cone. The morphometric approach, in which the grand potential of a system is described in four geometrical terms with corresponding thermodynamical coefficients, allows extrapolation or scaling from macroscopic system sizes to nanoscales. Since only a small number of fluid particles are involved in bubble formation, it is a pseudo phase transition, and the system is subjected to fluctuations between states with and without a bubble. Fluctuations are not included in a DFT treatment, which makes it possible to explore both states of the system in great detail, in contrast to computer simulations, in which averages might be obscured by fluctuations.

Critical behavior of a colloid-polymer mixture confined between walls

We investigate the influence of confinement on phase separation in colloid-polymer mixtures. To describe the particle interactions, the colloid-polymer model of Asakura and Oosawa [J. Chem. Phys. 22, 1255 (1954)] is used. Grand canonical Monte Carlo simulations are then applied to this model confined between two parallel hard walls, separated by a distance D = 5 colloid diameters. We focus on the critical regime of the phase separation and look for signs of crossover from three-dimensional (3D) Ising to two-dimensional (2D) Ising universality. To extract the critical behavior, finite size scaling techniques are used, including the recently proposed algorithm of Kim et al [Phys. Rev. Lett. 91, 065701 (2003)]. Our results point to "effective" critical exponents that differ profoundly from 3D Ising values, and that are already very close to 2D Ising values. In particular, we observe that the critical exponent of the order parameter in the confined system is smaller than in 3D bulk, yielding a "flatter" binodal. Our results also show an increase in the critical colloid packing fraction in the confined system with respect to the bulk. The latter seems consistent with theoretical expectations, although subtleties due to singularities in the critical behavior of the coexistence diameter cannot be ruled out.

Adsorption of fluids in slitlike pores containing a small amount of mobile ions

We apply density functional theory to investigate changes in the phase behavior of a fluid caused by the presence of mobile ions inside the pore. The approach has been based on the fundamental measure density functional theory and on the theory of nonuniform electrolytes developed recently by O. Pizio, A. Patrykiejew, S. Sokołowski [J. Chem. Phys. 121 (2005) 11,957]. We have evaluated capillary condensation phase diagrams for pores of different widths and for different concentrations of confined ions. The calculations have demonstrated that the presence of ions leads to lowering the critical temperature and to an increase of the value of the chemical potential at the capillary condensation point.

Capillary condensation of short-chain molecules

A density-functional study of capillary condensation of fluids of short-chain molecules confined to slitlike pores is presented. The molecules are modeled as freely jointed tangent spherical segments with a hard core and with short-range attractive interaction between all the segments. We investigate how the critical parameters of capillary condensation of the fluid change when the pore width decreases and eventually becomes smaller than the nominal linear dimension of the single-chain molecule. We find that the dependence of critical parameters for a fluid of dimers and of tetramers on pore width is similar to that of the monomer fluid. On the other hand, for a fluid of chains consisting of a larger number of segments we observe an inversion effect. Namely, the critical temperature of capillary condensation decreases with increasing pore width for a certain interval of values of the pore width. This anomalous behavior is also influenced by the interaction between molecules and pore walls. We attribute this behavior to the effect of conformational changes of molecules upon confinement.

Corrections to the Kelvin equation for long-range boundary fields

The properties of a simple fluid, or Ising magnet, confined in an L x infinity geometry, are studied by means of numerical density-matrix renormalization-group techniques. Whereas the particle-particle interactions are short ranged, the wall-particle interactions can be long ranged. Using a few complementary criteria the wetting transition line has been found for boundary fields with different ranges (including a marginal case) and its scaling form has been analyzed. A main goal was to analyze the influence of wetting on the scaling of the Kelvin equation, focusing on the leading correction term. We have found that at complete wetting the larger the range of the boundary fields, the smaller the region of anomalous corrections of type 1/L(4/3).

Near-critical confined fluids and Ising films: density-matrix renormalization-group study

Two-dimensional Ising strips subject to identical surface fields h(1)=h(2) > or =0 are studied for temperatures above and below the bulk critical temperature T(c) and a range of bulk fields h by means of the density-matrix renormalization-group method. In the case of nonvanishing surface fields, the near-critical behavior of the solvation force f(solv), total adsorption Gamma, inverse longitudinal correlation length xi(parallel)( -1) and specific heat C(H) is strongly influenced by the (pseudo) capillary condensation that occurs below T(c). We obtain scaling functions of f(solv), Gamma, and xi(parallel)(-1). C(H) exhibits a weakly rounded singularity on crossing the pseudocoexistence line. We contrast these results with those for the case of free boundaries where, for temperatures slightly below T(c), f(solv) and C(H) exhibit a sharp extremum away from h=0. Our results have direct repercussions for the properties of near-critical Ising films in three dimensions and we argue that the long-ranged solvation (Casimir) force in confined fluids should be more attractive in the neighborhood of the capillary critical point than exactly at the bulk critical point.

Influence of capillary condensation on the near-critical solvation force

We argue that in a fluid, or magnet, confined by adsorbing walls which favor liquid, or the (+) phase, the solvation (Casimir) force in the vicinity of the critical point is strongly influenced by capillary condensation which occurs below the bulk critical temperature T(c). At T slightly below and above T(c), a small bulk field h<0, which favors gas, or the (-) phase, leads to residual condensation and a solvation force which is much more attractive (at the same large wall separation) than that found exactly at the critical point. Our predictions are supported by results obtained from density-matrix renormalization-group calculations in a two-dimensional Ising strip subject to identical surface fields.

Coexistence of excited states in confined ising systems

Using the density-matrix renormalization-group method we study the two-dimensional Ising model in an infinite strip geometry with free-boundary conditions. The renormalization scheme enables us to consider systems of width up to 300 (lattice spacings) and study the influence of the bulk magnetic field on correlation function structure for all temperatures. From our numerical results we have determined the crossover line for the correlation length related to the coexistence of the excited states. A detailed scaling study of this line is performed. Our numerical results support and further specify previous conclusions reached by Abraham, Parry, and Upton based on the bubble model of correlations.