Condensation in a capped capillary is a continuous critical phenomenon

A New Statistical Theory for Constructing Sorption Isotherms in Mesoporous Structures Represented by Bethe Lattices

In the present work, a new statistical theory is developed to describe adsorption and desorption in mesoporous materials (pore sizes ranging from 2 to 50 nm) represented by pore networks in the form of Bethe lattices. The new theory is an extension of a previous theory applied for Statistically Disordered Chain Model (SDCM) structures and incorporates the cooperative effects emerging during phase transitions in pore networks. The theory is validated against simulations and algorithmic models that describe sorption of lattice and real fluids in Bethe lattices. It is seen that the pore network coordination number, or pore connectivity, z, has a significant impact on two important processes observed in pore networks: pore assisting condensation during adsorption and evaporation by percolation during desorption. The inclusion of pore connectivity in the earlier developed framework accounting for cooperativity effects is an important step, rendering the existing models to mimic fluid behavior in real materials more accurately. Hence, the new theory inherently contains all essential elements that may offer the extraction of more reliable pore size distributions utilizing both the adsorption and desorption branches of the isotherm.

Critical effects and scaling at meniscus osculation transitions

We propose a simple scaling theory describing critical effects at rounded meniscus osculation transitions which occur when the Laplace radius of a condensed macroscopic drop of liquid coincides with the local radius of curvature R_{w} in a confining parabolic geometry. We argue that the exponent β_{osc} characterizing the scale of the interfacial height ℓ_{0}∝R_{w}^{β_{osc}} at osculation, for large R_{w}, falls into two regimes representing fluctuation-dominated and mean-field-like behavior, respectively. These two regimes are separated by an upper critical dimension, which is determined here explicitly and depends on the range of the intermolecular forces. In the fluctuation-dominated regime, representing the universality class of systems with short-range forces, the exponent is related to the value of the interfacial wandering exponent ζ by β_{osc}=3ζ/(4-ζ). In contrast, in the mean-field regime, which was not previously identified and which occurs for systems with longer-range forces (and higher dimensions), the exponent β_{osc} takes the same value as the exponent β_{s}^{co} for complete wetting, which is determined directly by the intermolecular forces. The prediction β_{osc}=3/7 in d=2 for systems with short-range forces (corresponding to ζ=1/2) is confirmed using an interfacial Hamiltonian model which determines the exact scaling form for the decay of the interfacial height probability distribution function. A numerical study in d=3, based on a microscopic model density-functional theory, determines that β_{osc}≈β_{s}^{co}≈0.326 close to the predicted value of 1/3 appropriate to the mean-field regime for dispersion forces.

Phase behavior of fluids in undulated nanopores

The geometry of walls forming a narrow pore may qualitatively affect the phase behavior of the confined fluid. Specifically, the nature of condensation in nanopores formed of sinusoidally shaped walls (with amplitude A and period P) is governed by the wall mean separation L as follows. For L>L_{t}, where L_{t} increases with A, the pores exhibit standard capillary condensation similar to planar slits. In contrast, for L<L_{t}, the condensation occurs in two steps, such that the fluid first condenses locally via bridging transition connecting adjacent crests of the walls, before it condenses globally. For the marginal value of L=L_{t}, all the three phases (gaslike, bridge, and liquidlike) may coexist. We show that the locations of the phase transitions can be described using geometric arguments leading to modified Kelvin equations. However, for completely wet walls, to which we focus on, the phase boundaries are shifted significantly due to the presence of wetting layers. In order to take this into account, mesoscopic corrections to the macroscopic theory are proposed. The resulting predictions are shown to be in a very good agreement with a density-functional theory even for molecularly narrow pores. The limits of stability of the bridge phase, controlled by the pore geometry, is also discussed in some detail.

Effect of Confinement on Capillary Phase Transition in Granular Aggregates

Using a 3D mean-field lattice-gas model, we analyze the effect of confinement on the nature of capillary phase transition in granular aggregates with varying disorder and their inverse porous structures obtained by interchanging particles and pores. Surprisingly, the confinement effects are found to be much less pronounced in granular aggregates as opposed to porous structures. We show that this discrepancy can be understood in terms of the surface-surface correlation length with a connected path through the fluid domain, suggesting that this length captures the true degree of confinement. We also find that the liquid-gas phase transition in these porous materials is of second order nature near capillary critical temperature, which is shown to represent a true critical temperature, i.e., independent of the degree of disorder and the nature of the solid matrix, discrete or continuous. The critical exponents estimated here from finite-size scaling analysis suggest that this transition belongs to the 3D random field Ising model universality class as hypothesized by F. Brochard and P.G. de Gennes, with the underlying random fields induced by local disorder in fluid-solid interactions.

Filling, depinning, unbinding: Three adsorption regimes for nanocorrugated substrates

We study adsorption at periodically corrugated substrates formed by scoring rectangular grooves into a planar solid wall which interacts with the fluid via long-range (dispersion) forces. The grooves are assumed to be macroscopically long but their depth, width, and separations can all be molecularly small. We show that the entire adsorption process can be divided into three parts consisting of (i) filling the grooves by a capillary liquid; (ii) depinning of the liquid-gas interface from the wall edges; and (iii) unbinding of the interface from the top of the wall, which is accompanied by a rapid but continuous flattening of its shape. Using a nonlocal density functional theory and mesoscopic interfacial models all the regimes are discussed in some detail to reveal the complexity of the entire process and subtle aspects that affect its behavior. In particular, it is shown that the nature of the depinning phenomenon is governed by the width of the wall pillars (separating grooves), while the width of the grooves only controls the location of the depinning first-order transition, if present.

Scaling of wetting and prewetting transitions on nanopatterned walls

We consider a nanopatterned planar wall consisting of a periodic array of stripes of width L, which are completely wet by liquid (contact angle θ=0), separated by regions of width D which are completely dry (contact angle θ=π). Using microscopic density functional theory, we show that, in the presence of long-ranged dispersion forces, the wall-gas interface undergoes a first-order wetting transition, at bulk coexistence as the separation D is reduced to a value D_{w}∝lnL, induced by the bridging between neighboring liquid droplets. Associated with this is a line of prewetting transitions occurring off coexistence. By varying the stripe width L, we show that the prewetting line shows universal scaling behavior and data collapse. This verifies predictions based on mesoscopic models for the scaling properties associated with finite-size effects at complete wetting including the logarithmic singular contribution to the surface free energy.

Bridging of liquid drops at chemically structured walls

Using mesoscopic interfacial models and microscopic density functional theory we study fluid adsorption at a dry wall decorated with three completely wet stripes of width L separated by distances D_{1} and D_{2}. The stripes interact with the fluid with long-range forces inducing a large finite-size contribution to the surface free energy. We show that this nonextensive free-energy contribution scales with lnL and drives different types of bridging transition corresponding to the merging of liquid drops adsorbed at neighboring wetting stripes when the separation between them is molecularly small. We determine the surface phase diagram and show that this exhibits two triple points, where isolated drops, double drops, and triple drops coexist. For the symmetric case, D_{1}=D_{2}≡D, our results also confirm that the equilibrium droplet configuration always has the symmetry of the substrate corresponding to either three isolated drops when D is large or a single triple drop when D is small; however, symmetry-broken configurations do occur in a metastable part of the phase diagram which lies very close to the equilibrium-bridging phase boundary. Implications for phase transitions on other types of patterned surface are considered.

Adsorption on alumina nanopores with conical shape

Adsorption on porous solids depends on the morphology of the pores, the cylindrical one being the most studied in the literature. In this work, we present the first experimental investigation of adsorption and evaporation on conical nanopores produced by anodization of aluminium oxide. The pores are about 50 μm long, with the wide ends having a diameter of ∼79 nm and the narrow ones of ∼30 nm. Three different pores configurations are considered: open at both ends, open only at the narrow end and open only at the wide end. Despite the very small value of the conical angle α, estimated to be ∼0.06°, just barely above α = 0° corresponding to a cylindrical pore, the adsorption isotherms look strikingly different from those measured on cylindrical pores of similar size. First of all, the hysteresis loops of the conical pores with two open ends and with open wide ends practically coincide. Furthermore, they are narrower and the adsorption and evaporation branches are broader than those of the cylindrical pores with similar size. Finally, conical pores with open narrow ends exhibit a large hysteresis indicative of pore blocking. To unravel the mechanisms underlying adsorption and evaporation in such conical pores, we also report complementary results obtained using on-lattice grand canonical Monte Carlo simulations.

Microscopic aspects of wetting using classical density functional theory

Wetting is a rather efficient mechanism for nucleation of a phase (typically liquid) on the interface between two other phases (typically solid and gas). In many experimentally accessible cases of wetting, the interplay between the substrate structure, and the fluid-fluid and fluid-substrate intermolecular interactions brings about an entire 'zoo' of possible fluid configurations, such as liquid films with a thickness of a few nanometers, liquid nanodrops and liquid bridges. These fluid configurations are often associated with phase transitions occurring at the solid-gas interface and at lengths of just several molecular diameters away from the substrate. In this special issue article, we demonstrate how a fully microscopic classical density-functional framework can be applied to the efficient, rational and systematic exploration of the rich phase space of wetting phenomena. We consider a number of model prototype systems such as wetting on a planar wall, a chemically patterned wall and a wedge. Through density-functional computations we demonstrate that for these simply structured substrates the behaviour of the solid-gas interface is already highly complex and non-trivial.

Continuous condensation in nanogrooves

We consider condensation in a capillary groove of width L and depth D, formed by walls that are completely wet (contact angle θ=0), which is in a contact with a gas reservoir of the chemical potential μ. On a mesoscopic level, the condensation process can be described in terms of the midpoint height ℓ of a meniscus formed at the liquid-gas interface. For macroscopically deep grooves (D→∞), and in the presence of long-range (dispersion) forces, the condensation corresponds to a second-order phase transition, such that ℓ∼(μ_{cc}-μ)^{-1/4} as μ→μ_{cc}^{-} where μ_{cc} is the chemical potential pertinent to capillary condensation in a slit pore of width L. For finite values of D, the transition becomes rounded and the groove becomes filled with liquid at a chemical potential higher than μ_{cc} with a difference of the order of D^{-3}. For sufficiently deep grooves, the meniscus growth initially follows the power law ℓ∼(μ_{cc}-μ)^{-1/4}, but this behavior eventually crosses over to ℓ∼D-(μ-μ_{cc})^{-1/3} above μ_{cc}, with a gap between the two regimes shown to be δ[over ¯]μ∼D^{-3}. Right at μ=μ_{cc}, when the groove is only partially filled with liquid, the height of the meniscus scales as ℓ^{*}∼(D^{3}L)^{1/4}. Moreover, the chemical potential (or pressure) at which the groove is half-filled with liquid exhibits a nonmonotonic dependence on D with a maximum at D≈3L/2 and coincides with μ_{cc} when L≈D. Finally, we show that condensation in finite grooves can be mapped on the condensation in capillary slits formed by two asymmetric (competing) walls a distance D apart with potential strengths depending on L. All these predictions, based on mesoscopic arguments, are confirmed by fully microscopic Rosenfeld's density functional theory with a reasonable agreement down to surprisingly small values of both L and D.

Modified Kelvin Equations for Capillary Condensation in Narrow and Wide Grooves

We consider the location and order of capillary condensation transitions occurring in deep grooves of width L and depth D. For walls that are completely wet by liquid (contact angle θ=0) the transition is continuous and its location is not sensitive to the depth of the groove. However, for walls that are partially wet by liquid, where the transition is first order, we show that the pressure at which it occurs is determined by a modified Kelvin equation characterized by an edge contact angle θ_{E} describing the shape of the meniscus formed at the top of the groove. The dependence of θ_{E} on the groove depth D relies, in turn, on whether corner menisci are formed at the bottom of the groove in the low density gaslike phase. While for macroscopically wide grooves these are always present when θ<45° we argue that their formation is inhibited in narrow grooves. This has a number of implications including that the local pinning of the meniscus and location of the condensation transition is different depending on whether the contact angle is greater or less than a universal value θ^{*}≈31°. Our arguments are supported by detailed microscopic density functional theory calculations that show that the modified Kelvin equation remains highly accurate even when L and D are of the order of tens of molecular diameters.

Adsorption on Nanopores of Different Cross Sections Made by Electron Beam Nanolithography

Adsorption on nanoporous matrices is characterized by a pronounced hysteresis loop in the adsorption isotherm, when the substrate is loaded and unloaded with adsorbate, the origin of which is a matter of immense debate in the literature. In this work, we report a study of argon adsorption at 85 K on nonconnecting nanopores with one end closed to the surrounding where the effects of different pore cross sections fabricated by electron beam lithography (EBL) are investigated. A polymethylmethacrylate (PMMA) resist is deposited on the electrodes of a sensitive quartz crystal microbalance without degradation of the resonance quality factor or the long-term and short-term stabilities of the device even at cryogenic temperatures. Four different pores' cross sections: circular, square, rectangular, and triangular, are produced from EBL, and the isotherms for these pore shapes exhibit pronounced hysteresis loops whose adsorption and desorption branches are nearly vertical and have almost the same slopes. No difference is observed in the hysteresis loops of the isotherms for the pores with triangular and square cross sections, whereas the hysteresis loop for the pore with circular cross sections is much narrower, suggesting that they are more regular than the other pores. All of these observations suggest that the hysteresis behavior resulted mainly from microscopic geometric irregularities present in these porous matrices.

Scaling behavior of thin films on chemically heterogeneous walls

We study the adsorption of a fluid in the grand canonical ensemble occurring at a planar heterogeneous wall which is decorated with a chemical stripe of width L. We suppose that the material of the stripe strongly preferentially adsorbs the liquid in contrast to the outer material which is only partially wet. This competition leads to the nucleation of a droplet of liquid on the stripe, the height h_{m} and shape of which (at bulk two-phase coexistence) has been predicted previously using mesoscopic interfacial Hamiltonian theory. We test these predictions using a microscopic Fundamental Measure Density Functional Theory which incorporates short-ranged fluid-fluid and fully long-ranged wall-fluid interactions. Our model functional accurately describes packing effects not captured by the interfacial Hamiltonian but still we show that there is excellent agreement with the predictions h_{m}≈L^{1/2} and for the scaled circular shape of the drop even for L as small as 50 molecular diameters. For smaller stripes the droplet height is considerably lower than that predicted by the mesoscopic interfacial theory. Phase transitions for droplet configurations occurring on substrates with multiple stripes are also discussed.

Edge contact angle and modified Kelvin equation for condensation in open pores

We consider capillary condensation transitions occurring in open slits of width L and finite height H immersed in a reservoir of vapor. In this case the pressure at which condensation occurs is closer to saturation compared to that occurring in an infinite slit (H=∞) due to the presence of two menisci that are pinned near the open ends. Using macroscopic arguments, we derive a modified Kelvin equation for the pressure p_{cc}(L;H) at which condensation occurs and show that the two menisci are characterized by an edge contact angle θ_{e} that is always larger than the equilibrium contact angle θ, only equal to it in the limit of macroscopic H. For walls that are completely wet (θ=0) the edge contact angle depends only on the aspect ratio of the capillary and is well described by θ_{e}≈sqrt[πL/2H] for large H. Similar results apply for condensation in cylindrical pores of finite length. We test these predictions against numerical results obtained using a microscopic density-functional model where the presence of an edge contact angle characterizing the shape of the menisci is clearly visible from the density profiles. Below the wetting temperature T_{w} we find very good agreement for slit pores of widths of just a few tens of molecular diameters, while above T_{w} the modified Kelvin equation only becomes accurate for much larger systems.

Classical density functional study of wetting transitions on nanopatterned surfaces

Even simple fluids on simple substrates can exhibit very rich surface phase behaviour. To illustrate this, we consider fluid adsorption on a planar wall chemically patterned with a deep stripe of a different material. In this system, two phase transitions compete: unbending and pre-wetting. Using microscopic density-functional theory, we show that, for thin stripes, the lines of these two phase transitions may merge, leading to a new two-dimensional-like wetting transition occurring along the walls. The influence of intermolecular forces and interfacial fluctuations on this phase transition and at complete pre-wetting are considered in detail.

Complete prewetting

We study continuous interfacial transitions, analagous to two-dimensional complete wetting, associated with the first-order prewetting line, which can occur on steps, patterned walls, grooves and wedges, and which are sensitive to both the range of the intermolecular forces and interfacial fluctuation effects. These transitions compete with wetting, filling and condensation producing very rich phase diagrams even for relatively simple prototypical geometries. Using microscopic classical density functional theory to model systems with realistic Lennard-Jones fluid-fluid and fluid-substrate intermolecular potentials, we compute mean-field fluid density profiles, adsorption isotherms and phase diagrams for a variety of confining geometries.

Density functional study of condensation in capped capillaries

We study liquid adsorption in narrow rectangular capped capillaries formed by capping two parallel planar walls (a slit pore) with a third wall orthogonal to the two planar walls. The most important transition in confined fluids is arguably condensation, where the pore becomes filled with the liquid phase which is metastable in the bulk. Depending on the temperature T, the condensation in capped capillaries can be first-order (at T≤Tcw) or continuous (at T>Tcw), where Tcw is the capillary wetting temperature. At T>Tcw, the capping wall can adsorb mesoscopic amounts of metastable under-condensed liquid. The onset of condensation is then manifested by the continuous unbinding of the interface between the liquid adsorbed on the capping wall and the gas filling the rest of the capillary volume. In wide capped capillaries there may be a remnant of wedge filling transition, which is manifested by the adsorption of liquid drops in the corners. Our classical statistical mechanical treatment predicts a possibility of three-phase coexistence between gas, corner drops and liquid slabs adsorbed on the capping wall. In sufficiently wide capillaries we find that thick prewetting films of finite length may be nucleated at the capping wall below the boundary of the prewetting transition. Prewetting then proceeds in a continuous manner manifested by the unbinding interface between the thick and thin films adsorbed on the side walls. Our analysis is based on a detailed numerical investigation of the density functional theory for the fluid equilibria for a number of illustrative case studies.

Filling transitions in acute and open wedges

We present numerical studies of first-order and continuous filling transitions in wedges of arbitrary opening angle ψ, using a microscopic fundamental measure density functional model with short-ranged fluid-fluid forces and long-ranged wall-fluid forces. In this system the wetting transition characteristic of the planar wall-fluid interface is always first order regardless of the strength of the wall-fluid potential ɛ(w). In the wedge geometry, however, the order of the filling transition depends not only on ɛ(w) but also on the opening angle ψ. In particular we show that even if the wetting transition is strongly first order the filling transition is continuous for sufficient acute wedges. We show further that the change in the order of the transition occurs via a tricritical point as opposed to a critical end point. These results extend previous effective Hamiltonian predictions which were limited only to shallow wedges.

Adsorption in alumina pores open at one and at both ends

We have studied adsorption in regular, self-ordered alumina pores open at both ends or only at one end. The straight, non-connected pores have diameters ranging from 22 to 83 nm, with a relative dispersion below 1% in the pore size. Adsorption isotherms measured in open pores with a torsional microbalance show pronounced hysteresis loops characterized by nearly vertical and parallel adsorption and desorption branches. Blocking one end of the pores with glue has a strong influence on adsorption, as expected from classical macroscopic arguments. However, the experimental measurements show an unexpectedly rich phenomenology dependent on the pore size. For large pores (Dp ≥ 67 nm), the isotherms for closed end pores present much narrower hysteresis loops whose adsorption and desorption boundaries envelop the desorption branches of the isotherms for the corresponding open pores of the same size. The loop for small closed end pores (Dp = 22 nm) is slightly wider than that for open pores while the adsorption branches coincide. For large pores, in contrast, the desorption branches of pores with the same Dp overlap regardless of the pore opening. These observations are in agreement with our grand canonical Monte Carlo (GCMC) simulations for a cylindrical pore model with constrictions, suggesting that the alumina pores could be modeled using a constricted pore model whose adsorption isotherm depends on the ratio of the constriction size to the pore size (Dc/Dp).

Filling and wetting transitions on sinusoidal substrates: a mean-field study of the Landau-Ginzburg model

We study the interfacial phenomenology of a fluid in contact with a one-dimensional array of infinitely long grooves of sinusoidal section, characterized by the periodicity length L and amplitude A. The system is modelled by the Landau-Ginzburg-Wilson functional, with fluid-substrate couplings which control the wettability of the substrate. We investigate the filling and wetting phenomena within the mean-field approximation, and compare with the predictions of the macroscopic and interfacial Hamiltonian theories. For large values of L and under bulk coexistence conditions, we observe first-order filling transitions between dry (D) and partially filled (F) interfacial states, and wetting transitions between partially filled F and completely wet (W) interfacial states of the same order as for the flat substrate. Depending on the order of the wetting transition, the transition temperature is either shifted towards lower temperatures for first-order wetting or it coincides with the wetting temperature on the flat substrate for continuous wetting. On the other hand, if the groove height is of order of the correlation length, only wetting transitions between D and W states are observed under bulk coexistence conditions. For this case, the transition temperature shift obeys approximately Wenzel's phenomenological law if the substrate favors first-order wetting, but it remains unshifted for continuous wetting. The borderline between the small and large L regimes correspond to a D - F - W triple point if wetting is first-order, and a D - F critical point for continuous wetting. Beyond bulk coexistence conditions, filling and first-order wetting transitions continue into off-coexistence filling and prewetting lines, which end up at critical points. Our findings show that the macroscopic theory only describes accurately the filling transition close to bulk coexistence and large L, while microscopic structure of the fluid is essential to understand wetting and filling away from bulk coexistence.

Capillary contact angle in a completely wet groove

We consider the phase equilibria of a fluid confined in a deep capillary groove of width L with identical side walls and a bottom made of a different material. All walls are completely wet by the liquid. Using density functional theory and interfacial models, we show that the meniscus separating liquid and gas phases at two phase capillary coexistence meets the bottom capped end of the groove at a capillary contact angle θ(cap)(L) which depends on the difference between the Hamaker constants. If the bottom wall has a weaker wall-fluid attraction than the side walls, then θ(cap) > 0 even though all the isolated walls are themselves completely wet. This alters the capillary condensation transition which is now first order; this would be continuous in a capped capillary made wholly of either type of material. We show that the capillary contact angle θ(cap)(L) vanishes in two limits, corresponding to different capillary wetting transitions. These occur as the width (i) becomes macroscopically large, and (ii) is reduced to a microscopic value determined by the difference in Hamaker constants. This second wetting transition is characterized by large scale fluctuations and essential critical singularities arising from marginal interfacial interactions.

Condensation and evaporation transitions in deep capillary grooves

We study the order of capillary condensation and evaporation transitions of a simple fluid adsorbed in a deep capillary groove using a fundamental measure density functional theory (DFT). The walls of the capillary interact with the fluid particles via long-ranged, dispersion, forces while the fluid-fluid interaction is modelled as a truncated Lennard-Jones-like potential. We find that below the wetting temperature Tw condensation is first-order and evaporation is continuous with the metastability of the condensation being well described by the complementary Kelvin equation. In contrast above Tw both phase transitions are continuous and their critical singularities are determined. In addition we show that for the evaporation transition above Tw there is an elegant mapping, or covariance, with the complete wetting transition occurring at a planar wall. Our numerical DFT studies are complemented by analytical slab model calculations which explain how the asymmetry between condensation and evaporation arises out of the combination of long-ranged forces and substrate geometry.

Filling dynamics of closed end nanocapillaries

We have studied the filling dynamics of model capillaries using dynamic mean field theory for a confined lattice gas and Kawasaki dynamics simulations. We have found two different scenarios for filling of capped nanocapillaries from the vapor phase. As compared to channels with macroscopic width, in which the filling process occurs by the detachment of the meniscus from the cap, in mesoscopic channels there is an alternative mechanism associated with the spontaneous condensation of the liquid close to the pore opening and its subsequent growth toward the closed pore end. We show that these two scenarios have totally different filling dynamics, providing an additional mechanism for slow capillary condensation kinetics in nanoscopic objects.

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.

On the irreversibility of the adsorption isotherm in a closed-end pore

We present a simulation study of argon adsorption in a closed-end mesopore of uniform diameter in order to investigate the occurrence of hysteresis and propose two principal reasons for its existence: the variation in the shape and radius of curvature of the meniscus and the change in the packing of adsorbate during adsorption and desorption. This interpretation differs from classical theories that neglect both of these factors, and therefore find that adsorption-desorption in a closed-end pore is reversible. A detailed simulation study of the effects of temperature on the microscopic behavior of the adsorbate supports the interpretation proposed here.

Geometry-induced phase transition in fluids: capillary prewetting

We report a new first-order phase transition preceding capillary condensation and corresponding to the discontinuous formation of a curved liquid meniscus. Using a mean-field microscopic approach based on the density functional theory we compute the complete phase diagram of a prototypical two-dimensional system exhibiting capillary condensation, namely that of a fluid with long-ranged dispersion intermolecular forces which is spatially confined by a substrate forming a semi-infinite rectangular pore exerting long-ranged dispersion forces on the fluid. In the T-μ plane the phase line of the new transition is tangential to the capillary condensation line at the capillary wetting temperature T(cw). The surface phase behavior of the system maps to planar wetting with the phase line of the new transition, termed capillary prewetting, mapping to the planar prewetting line. If capillary condensation is approached isothermally with T>T(cw), the meniscus forms at the capping wall and unbinds continuously, making capillary condensation a second-order phenomenon. We compute the corresponding critical exponent for the divergence of adsorption.

Nematic wetting and filling of crenellated surfaces

We investigate nematic wetting and filling transitions of crenellated surfaces (rectangular gratings) by numerical minimization of the Landau-de Gennes free energy as a function of the anchoring strength, for a wide range of the surface geometrical parameters: depth, width, and separation of the crenels. We have found a rich phase behavior that depends in detail on the combination of the surface parameters. By comparison to simple fluids, which undergo a continuous filling or unbending transition, where the surface changes from a dry to a filled state, followed by a wetting or unbinding transition, where the thickness of the adsorbed fluid becomes macroscopic and the interface unbinds from the surface, nematics at crenellated surfaces reveal an intriguingly rich behavior: in shallow crenels only wetting is observed, while in deep crenels, only filling transitions occur; for intermediate surface geometrical parameters, a new class of filled states is found, characterized by bent isotropic-nematic interfaces, which persist for surfaces structured on large scales, compared to the nematic correlation length. The global phase diagram displays two wet and four filled states, all separated by first-order transitions. For crenels in the intermediate regime re-entrant filling transitions driven by the anchoring strength are observed.

The order of filling transitions in acute wedges

Using a square-gradient density functional model we test the prediction that the filling transition for a fluid in a wedge geometry changes from continuous to first-order as the wedge becomes more acute. Our numerical findings confirm such a change of order, but the value of the tilt angle at which it occurs, α* ≈ 45°, is considerably smaller than the original theoretical prediction. We critically reassess this work, which was based on allowing for the self-interaction of the fluid interface, and argue that the interfacial curvature and effective wavevector dependent surface tension can further lower the predicted value of α*, in keeping with our numerical findings. Interfacial fluctuation effects, occurring beyond mean-field level, are also discussed using effective Hamiltonian theory and are shown to substantially increase the value of α*.

Complete wetting of elastically responsive substrates

We analyze theoretically complete wetting of a substrate supporting an array of parallel, vertical plates which can tilt elastically. The adsorbed liquid tilts the plates, inducing clustering, and thus modifies the substrate geometry. In turn, this change in geometry alters the wetting properties of the substrate and, consequently, the adsorption of liquid. This geometry-wetting feedback loop leads to stepped adsorption isotherms with each step corresponding to an abrupt change in the substrate geometry. We discuss how this can be used for constructing substrates with tunable wetting and adsorption properties.

Scaling properties of fluid adsorption near the base of a cylinder

We consider the adsorption of fluid at the foot of a cylinder that protrudes from a flat substrate made of the same material. Provided the contact angle θ is small enough, a drop of liquid condenses near the base, the size of which can be determined using simple macroscopic arguments. The adsorption in this geometry shows scaling behavior related to a number of different interfacial phase transitions and, for systems with short-ranged forces, shows a remarkable property; for small θ, the height of the drop (measured from the base) and the width (measured from the cylinder axis) are near identical to expressions for the thickness and parallel correlation length for microscopic wetting films (at planar walls). The only difference is that the bulk correlation length is replaced by the radius of the cylinder. By taking into account the correct singular behavior of the line tension we show that this geometrical amplification of the microscopic lengths occurs for second-order, first-order, and complete wetting transitions, and is specific to three dimensions. Similar phenomena occurs for long-ranged forces, and shows crossover scaling behavior.

Capillary condensation of adsorbates in porous materials

Hysteresis in capillary condensation is important for the fundamental study and application of porous materials, and yet experiments on porous materials are sometimes difficult to interpret because of the many interactions and complex solid structures involved in the condensation and evaporation processes. Here we make an overview of the significant progress in understanding capillary condensation and hysteresis phenomena in mesopores that have followed from experiment and simulation applied to highly ordered mesoporous materials such as MCM-41 and SBA-15 over the last few decades.

The self-interaction of a fluid interface, the wavevector dependent surface tension and wedge filling

We argue that whenever an interface, separating bulk fluid phases, adopts a non-planar configuration (induced by a confining geometry or thermal fluctuations, say), the energy cost of it will contain a non-local self-interaction term. For systems with short-ranged forces and Ising symmetry, we determine the self-interaction by integrating out bulk-like degrees of freedom from a more microscopic Landau-Ginzburg-Wilson model. The self-interaction can be written in a simple diagrammatic form involving integrals over effective two-body forces acting at the interface and consistently accounts for a number of known features of the microscopic model, including the wavevector dependence of the surface tension describing the fluctuations of a near planar interface. When applied to wedge filling transitions, the self-interaction describes the attraction between the wetting films on either side of the wedge. We show that, for sufficiently acute wedges, this can alter the order of the filling phase transition.

Nanogeometry matters: unexpected decrease of capillary adhesion forces with increasing relative humidity

The sticking effect between hydrophilic surfaces occurring at increasing relative humidity (RH) is an everyday phenomenon with uncountable implications. Here experimental evidence is presented for a counterintuitive monotonous decrease of the capillary adhesion forces between hydrophilic surfaces with increasing RH for the whole humidity range. It is shown that this unexpected result is related to the actual shape of the asperity at the nanometer scale: a model based on macroscopic thermodynamics predicts this decrease in the adhesion force for a sharp object ending in an almost flat nanometer-sized apex, in full agreement with experiments. This anomalous decrease is due to the fact that a significant growth of the liquid meniscus formed at the contact region with increasing humidity is hindered for this geometry. These results are relevant in the analysis of the dynamical behavior of nanomenisci. They could also have an outstanding value in technological applications, since the undesirable sticking effect between surfaces occurring at increasing RH could be avoided by controlling the shape of the surface asperities at the nanometric scale.

Wetting of nanopatterned grooved surfaces

The wetting by perfluoromethylcyclohexane of a well-defined silicon grating with a channel width of 16 nm has been studied using transmission small angle x-ray scattering. Prefilling, capillary filling, and postfilling wetting regimes have been identified. A detailed comparison of the data with theory reveals the importance of long-ranged substrate-fluid and fluid-fluid interactions for determining the wetting behavior on these length scales, especially at the onset of capillary condensation and in the prefilling regime.

Normal capillary forces

A liquid meniscus between two lyophilic solid surfaces causes an attractive force, the capillary force. The meniscus can form by capillary condensation or by accumulation of adsorbed liquid. Under ambient conditions and between hydrophilic surfaces, capillary forces usually dominate over other surface forces. They are relevant in many processes occurring in nature and technical applications, for example the flow of granular materials and friction between surfaces. Here we review normal capillary forces, focusing on a quantitative description with continuum theory. After introducing the capillary force between spherical surfaces, we extend the discussion to other regular and irregular surfaces. The influence of surface roughness is considered. In addition to capillary forces at equilibrium, we also describe the process of meniscus formation. Assumptions, limits, and perspectives for future work are discussed.

Liquid-vapor coexistence at a mesoporous substrate

The condensation of hexane vapor onto a mesoporous Si substrate with a pore radius of 3.5 nm has been studied by means of volumetry and ellipsometry. The filling fraction of the pores and the coverage of the substrate have been determined. The coverage of the regime after the completion of capillary condensation has been compared to recent theoretical work.

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.

Intrusion of fluids into nanogrooves: how geometry determines the shape of the gas-liquid interface

We study the shape of gas-liquid interfaces forming inside rectangular nanogrooves (i.e., slit-pores capped on one end). On account of purely repulsive fluid-substrate interactions the confining walls are dry (i.e., wet by vapor) and a liquid-vapor interface intrudes into the nanogrooves to a distance determined by the pressure (i.e., chemical potential). By means of Monte Carlo simulations in the grand-canonical ensemble (GCEMC) we obtain the density rho(z) along the midline (x = 0) of the nanogroove for various geometries (i.e., depths D and widths L) of the nanogroove. We analyze the density profiles with the aid of an analytic expression which we obtain through a transfer-matrix treatment of a one-dimensional effective interface Hamiltonian. Besides geometrical parameters such as D and L , the resulting analytic expression depends on temperature T , densities of coexisting gas and liquid phases in the bulk rho g,l(x) and the interfacial tension gamma. The latter three quantities are determined in independent molecular dynamics simulations of planar gas-liquid interfaces. Our results indicate that the analytic formula provides an excellent representation of rho(z) as long as L is sufficiently small. At larger L the meniscus of the intruding liquid flattens. Under these conditions the transfer-matrix analysis is no longer adequate and the agreement between GCEMC data and the analytic treatment is less satisfactory.