Fact-checked by Grok 2 weeks ago

Capillary condensation

Capillary condensation is the process by which a vapor condenses into within the narrow pores or crevices of a porous at partial pressures below the saturation of the bulk , driven by the of the liquid-vapor in confined geometries. This phenomenon occurs preferentially in wetting systems where the is less than 90°, enabling the formation of stable liquid bridges or menisci that lower the equilibrium required for . The underlying mechanism is thermodynamically described by the , which relates the shift in to the pore radius: \ln(p/p_s) = -\frac{2\gamma V_m \cos \theta}{rRT}, where p is the equilibrium over the curved interface, p_s is the saturation , \gamma is the , V_m is the of the , \theta is the , r is the pore radius, R is the , and T is the . Building upon foundational contributions from Thomas Young (1805) on and Pierre-Simon (1806) on capillary forces, William Thomson () proposed it in linear form in 1871 to explain vapor equilibrium at curved surfaces, which later developed into the exponential form used today, with practical applications to porous media emerging in the early , such as Richard Zsigmondy's 1911 studies on silica gels. A hallmark of capillary condensation is the adsorption-desorption observed in isotherms, arising from the metastable states during filling and emptying of pores, which persists across diverse mesoporous materials like metal-organic frameworks and silicas. Capillary condensation plays a pivotal role in numerous scientific and engineering fields, influencing processes from gas adsorption in and to water management in soils and textiles. In energy applications, it affects recovery by altering permeability in nanopores and enhances in supercapacitors. Emerging uses include atmospheric harvesting via tailored porous structures and the controlled growth of , such as nanowires, underscoring its versatility in addressing challenges in and .

Fundamentals

Definition and Mechanism

Capillary condensation is the phenomenon in which a transitions to a within the narrow confines of pores or capillaries in a porous material, occurring at relative pressures (P_v/P_sat) less than unity, where P_v is the and P_sat is the saturation pressure of the bulk . This process is fundamentally driven by at the liquid-vapor interface and the associated Laplace pressure, which arises from the of the interface in confined geometries. The mechanism initiates with the adsorption of vapor molecules onto the pore walls, forming an initial thin liquid due to attractive intermolecular forces, primarily van der Waals interactions. As the relative rises, the thickens until it reaches an instability point, where of discrete liquid clusters occurs, often described as a morphological transition from a uniform layer to localized liquid phases. These clusters develop into curved , whose shape in pores generates a that reduces the equilibrium needed for , thereby promoting further liquid formation. Key physical drivers include (γ), which dictates the meniscus curvature and pressure differential, and confinement effects that enhance intermolecular forces; in hydrophilic pores, strong favors film stability and rapid filling, whereas hydrophobic pores exhibit delayed or altered dynamics due to reduced adhesion. The process unfolds step-by-step: vapor adsorption first establishes menisci at walls or entrances, creating localized regions of high that draw additional vapor inward via . This leads to the growth of liquid bridges connecting opposite surfaces, stabilized by the balance of and disjoining pressures from intermolecular forces. As pressure increases, these bridges expand and coalesce, ultimately filling the volume with liquid condensate. The quantitative foundation for this pressure reduction is captured by the , relating vapor pressure to interface . This effect was first theoretically described in the 19th century by William Thomson (later Lord Kelvin) in the context of atmospheric vapor equilibrium at curved liquid surfaces.

Thermodynamic Principles

Capillary condensation is governed by thermodynamic principles that minimize the Gibbs free energy of the system, where the phase transition from vapor to liquid in confined geometries reduces the overall free energy through the formation of curved menisci. This process balances the bulk vapor-liquid free energy gain against surface energy contributions from the liquid-vapor interface and the solid-liquid/solid-vapor interactions. In porous media, the curved interface lowers the system's total Gibbs free energy compared to the non-condensed state when the chemical potential difference favors condensation, particularly for wetting surfaces where the spreading coefficient is positive. The Laplace pressure plays a central role, creating a difference across the curved given by \Delta P = \frac{2\gamma \cos\theta}{r}, where \gamma is the liquid-vapor , \theta is the , and r is the pore radius. This difference induces within the pores, as the liquid phase experiences a lower than the surrounding vapor, stabilizing the condensed phase at vapor pressures below the bulk saturation value. For hydrophilic pores with \theta < 90^\circ, this effect promotes condensation by enhancing the mechanical stability of the liquid bridge. At equilibrium, the chemical potential of the vapor must equal that of the condensed liquid phase within the confinement, ensuring no net driving force for phase change. This equality shifts the coexistence condition due to the curvature-induced modifications to the interfacial free energy, distinguishing confined systems from bulk behavior. In small pores, nucleation faces significant energy barriers, as the activation energy for forming a critical liquid nucleus scales with the pore size and is modulated by wetting properties; for partial wetting (\theta > 0), higher barriers delay , leading to metastable states. Unlike bulk condensation, which occurs at the saturation vapor pressure, confinement alters the by the Laplace pressure and surface terms, enabling condensation at lower relative pressures and introducing due to these energy barriers. This thermodynamic shift is particularly pronounced in nanopores, where surface effects dominate over bulk contributions.

Theoretical Framework

Kelvin Equation

The is derived from the thermodynamic requirement that the chemical potentials of the vapor and liquid phases are equal at equilibrium during capillary condensation. This builds on the principle of balance between coexisting phases. For the vapor phase, assumed to behave as an , the is expressed as
\mu_v = \mu^0(T) + RT \ln\left(\frac{P_v}{P_\mathrm{sat}}\right),
where \mu^0(T) is the standard at temperature T, R is the universal , P_v is the equilibrium in the , and P_\mathrm{sat} is the over a flat - .
For the phase in the bulk at , \mu_l^\mathrm{bulk} = \mu^0(T). In a capillary , the experiences a pressure reduction due to the Laplace pressure across the curved at the - , given by \Delta P = \frac{2\gamma \cos\theta}{r} for a cylindrical , where \gamma is the - , \theta is the , and r is the . The in the is thus P_l = P_v - \Delta P. Assuming the is incompressible with constant V_m, the of the confined becomes
\mu_l = \mu_l^\mathrm{bulk} + V_m (P_l - P_\mathrm{sat}) = \mu^0(T) + V_m (P_v - \Delta P - P_\mathrm{sat}).
Setting \mu_v = \mu_l at yields
RT \ln\left(\frac{P_v}{P_\mathrm{sat}}\right) = V_m (P_v - P_\mathrm{sat} - \Delta P).
Given that P_v \approx P_\mathrm{sat} for modest curvatures (small \Delta P), the term V_m (P_v - P_\mathrm{sat}) is second-order small and can be neglected, resulting in the :
\ln\left(\frac{P_v}{P_\mathrm{sat}}\right) = -\frac{2 \gamma V_m \cos\theta}{r RT}.
This form indicates that capillary condensation occurs at a relative P_v / P_\mathrm{sat} < 1, lower than the bulk saturation pressure, due to the concave curvature stabilizing the liquid phase in the pore. The equation originates from the work of William Thomson (Lord Kelvin) in 1871, who first related vapor to curved surfaces.
The derivation relies on several key assumptions: ideal cylindrical pores with smooth walls; partial to complete wetting (\theta \leq 90^\circ); incompressible liquid with bulk thermodynamic properties; ideal gas behavior for the vapor; and negligible gas adsorption or multilayer film effects on the pore walls. For complete wetting (\theta = 0), the equation simplifies to
\ln\left(\frac{P_v}{P_\mathrm{sat}}\right) = -\frac{2 \gamma V_m}{r RT}.
In spherical pore geometry, the meniscus approximates a hemisphere, leading to a mean curvature of $2/r and \Delta P = 2\gamma / r (assuming \theta = 0), which yields the same simplified form as the cylindrical case with complete wetting:
\ln\left(\frac{P_v}{P_\mathrm{sat}}\right) = -\frac{2 \gamma V_m}{r RT}.
The Kelvin equation holds reliably for mesopores, defined as having widths between 2 nm and 50 nm, where continuum thermodynamics applies and surface curvature dominates without significant molecular-scale deviations. It breaks down in micropores (widths < 2 nm), where enhanced wall-fluid interactions, molecular discreteness, and non-bulk liquid densities invalidate the assumptions, necessitating approaches like density functional theory.

Relative Pressure Dependence

The Kelvin equation predicts that the relative vapor pressure P_v / P_{\text{sat}} at which capillary condensation occurs in a pore of radius r follows an exponential dependence given by \ln(P_v / P_{\text{sat}}) \propto -1/r, such that smaller pores exhibit condensation at significantly lower relative pressures than the bulk saturation value of P_v / P_{\text{sat}} = 1. This shift arises from the enhanced curvature of the liquid-vapor meniscus in confined spaces, which lowers the equilibrium vapor pressure required for phase transition. Graphically, this dependence is evident in type IV adsorption isotherms typical of mesoporous materials, where the capillary condensation manifests as a steep uptake step at P_v / P_{\text{sat}} < 1, with the position of the step inversely related to the pore radius; for instance, narrower pores show the transition at lower relative pressures, while bulk-like behavior approaches P_v / P_{\text{sat}} = 1 for very large pores. These isotherms allow visual comparison of confined versus bulk condensation, highlighting how pore confinement stabilizes the liquid phase at sub-saturation conditions. Temperature influences the relative pressure shift through the RT term in the Kelvin equation, where higher temperatures decrease the magnitude of \ln(P_v / P_{\text{sat}}), bringing the condensation closer to bulk saturation due to the dominance of thermal energy over interfacial effects. For quantitative examples, in a 10 nm radius cylindrical pore with water at 25°C (assuming complete wetting), P_v / P_{\text{sat}} \approx 0.90; similarly, for nitrogen at 77 K in the same pore size, P_v / P_{\text{sat}} \approx 0.90, and for argon at 87 K, values are comparably close to 0.90 under standard conditions. This pressure-radius relationship underpins the predictive modeling of adsorption isotherms in porous materials, forming the basis for techniques like the Barrett-Joyner-Halenda (BJH) method to derive pore size distributions from experimental data.

Curvature and Contact Angle Effects

In capillary condensation within cylindrical pores, the curvature of the liquid-vapor meniscus plays a critical role in determining the condensation pressure, with the effective pore radius accounting for the pre-existing adsorbate layer thickness t. Specifically, the effective radius is given by r_{\text{eff}} = r - t, where r is the nominal pore radius, allowing the Kelvin equation to predict the onset of condensation more accurately by considering the reduced space available for the meniscus formation. During the filling process, the meniscus transitions from a cylindrical shape, which dominates at lower relative pressures and features a single principal radius of curvature, to a spherical shape as the pore fills, reflecting the evolving geometry and increasing the mean curvature, which shifts the equilibrium vapor pressure. The contact angle \theta, which quantifies surface wettability, modulates the pressure required for through its influence on the meniscus geometry. For hydrophilic surfaces where \theta < 90^\circ, the cosine term \cos \theta > 0 enhances the capillary forces, promoting at lower relative pressures compared to flat interfaces. In contrast, hydrophobic surfaces with \theta > 90^\circ yield \cos \theta < 0, which inhibits by requiring higher pressures or preventing it altogether in some cases. The standard Kelvin equation is modified to incorporate non-zero contact angles, yielding the form for a spherical meniscus: \ln\left(\frac{P_v}{P_{\text{sat}}}\right) = -\frac{2 \gamma V_m \cos \theta}{r R T} where \gamma is the surface tension, V_m is the molar volume of the liquid, R is the gas constant, and T is the temperature; for cylindrical menisci during initial filling, the factor of 2 is replaced by 1 to reflect the geometry: \ln\left(\frac{P_v}{P_{\text{sat}}}\right) = -\frac{ \gamma V_m \cos \theta}{r R T}. In real materials, \theta varies significantly: hydrophilic silica typically exhibits \theta \approx 0^\circ, facilitating condensation in pores as small as 2-5 nm at relative pressures near 0.4, while hydrophobic carbon surfaces show \theta \approx 80^\circ - 90^\circ, shifting thresholds to higher pressures and altering pore-filling dynamics. Experimentally, the contact angle is often determined by fitting adsorption isotherms to modified Kelvin models, where deviations in pore-filling thresholds from ideal predictions reveal \theta values; for instance, nitrogen adsorption on silica gels at 77 K allows extraction of \theta by comparing observed condensation steps to calculated pressures, confirming hydrophilic behavior that lowers filling thresholds by up to 20% relative to assumed \theta = 0^\circ. This approach validates how \theta directly impacts the relative humidity or pressure at which pores fill, with higher \theta delaying condensation in controlled mesoporous samples. In boundary cases, complete wetting (\theta = 0^\circ) leads to unbounded multilayer adsorption prior to bulk condensation, forming thick precursor films that effectively reduce the pore radius and enhance capillary forces, as seen in systems like water on clean silica where film thicknesses exceed 1 nm at sub-saturation pressures. Partial wetting ($0^\circ < \theta < 90^\circ) limits multilayer growth to finite thicknesses, resulting in sharper condensation transitions without extensive pre-wetting layers, which has implications for adsorption capacity in porous materials by altering the balance between surface coverage and capillary filling.

Real-System Variations

Pore Geometry Influences

In real porous materials, deviations from ideal cylindrical pore geometries significantly alter the dynamics of capillary condensation by influencing the shape and stability of the liquid-vapor meniscus. For instance, ink-bottle pores, characterized by a wide body connected to a narrow neck, lead to condensation initiating in the neck at a higher relative pressure corresponding to the smaller neck radius, followed by rapid filling of the larger body once the meniscus advances, resulting in distinct steps in adsorption isotherms. Slit-like pores, in contrast to cylindrical ones, promote condensation at lower relative pressures due to the attraction between opposing walls forming a flatter meniscus, whereas spherical or more closed geometries require higher pressures for meniscus formation because of increased curvature constraints. These shape-dependent meniscus configurations shift the equilibrium condensation pressure away from predictions based on uniform cylinders, as demonstrated in simulations of linear mesopores where slit-shaped pores exhibit earlier filling than cylindrical or spherical counterparts. Cylindrical pores with conical variations or tapered geometries further modify the filling sequence, with condensation progressing along the pore axis in a pressure-dependent manner that creates multiple plateaus in isotherms, reflecting sequential meniscus advancements from wider to narrower sections. In such non-uniform structures, the local radius variation along the pore axis, denoted as r(z), dictates the progression of the liquid bridge, leading to delayed or accelerated filling compared to straight cylindrical pores of equivalent average size. Experimental isotherms from materials like (cylindrical) versus disordered silicas (conical or irregular) reveal these effects through inflection points indicating geometry-driven transitions. Pre-adsorbed multilayer films on pore walls reduce the effective pore size, shifting the condensation threshold to higher relative pressures by effectively decreasing the radius available for meniscus formation; this multilayer thickness, often derived from t-curves, can be 0.3–0.5 nm for nitrogen at 77 K, altering the apparent geometry. Modeling these influences typically involves modifications to the to account for varying r(z), such as integrating local curvature terms along the pore length, though classical approaches like assume cylindrical symmetry and thus underestimate sizes in irregular pores by up to 20–30%. Numerical simulations using non-local density functional theory () or grand canonical Monte Carlo () provide more accurate predictions for complex shapes, capturing multilayer effects and meniscus evolution in slit-like or ink-bottle geometries without relying on idealized assumptions. Experimental characterization via methods like BET surface area analysis or BJH pore size distribution extracts geometry insights from isotherm steps; for example, sharp plateaus in nitrogen adsorption at 77 K on ordered mesoporous silicas indicate uniform cylindrical pores, while broader or multi-step features signal irregular shapes like slits or bottles, validated against TEM imaging. These techniques, refined in seminal works, highlight how pore geometry controls condensation accessibility and thus material uptake capacities.

Hysteresis Behavior

Hysteresis in capillary condensation manifests as a discrepancy between the adsorption (condensation) and desorption (evaporation) branches of the isotherm, where liquid condenses in pores at a higher relative pressure than the pressure at which it evaporates. This behavior arises from the advancing meniscus during adsorption, which forms at higher pressures due to favorable wetting conditions, versus the receding meniscus during desorption, which persists at lower pressures owing to metastable configurations. The primary causes of this hysteresis include pore blocking in constricted geometries, such as ink-bottle pores with narrow necks, where adsorbed liquid in larger cavities blocks evaporation from the wide body until the meniscus in the neck reaches a critical curvature. Metastable menisci also contribute, as the liquid-vapor interface remains stable in non-equilibrium states during desorption, delaying evaporation. Additionally, cavitation in wider pores during desorption triggers sudden vapor bubble formation when the tensile stress in the metastable liquid exceeds its mechanical stability, leading to abrupt emptying independent of pore connectivity. Hysteresis loops exhibit distinct characteristics depending on the underlying mechanism: Type A loops, associated with pore blocking in neck-dominated structures, show a steep desorption branch at lower pressures with gradual adsorption, often observed in materials with uniform narrow openings. In contrast, Type B loops arise from cavitation in expansive pore bodies, featuring parallel adsorption and desorption branches separated by a wide pressure gap, particularly in disconnected or wide-pore systems. The width of these loops is influenced by the pore size distribution, with broader distributions yielding wider hysteresis due to sequential filling and emptying across varied pore dimensions. Theoretical modeling of hysteresis often employs the independent domain theory, originally developed by Everett for parallel, non-interacting pores, which treats the porous material as an assembly of independent domains that fill or empty at specific threshold pressures, enabling prediction of loop shapes and reversibility. This framework extends to scanning curves in partial adsorption-desorption cycles, where intermediate paths trace reversible filling within domains between primary loop boundaries, providing insights into metastable state distributions. Experimentally, hysteresis is measured through nitrogen physisorption isotherms at 77 K, where the adsorption branch reflects equilibrium capillary condensation via the , while the desorption branch reveals irreversible processes like blocking or cavitation, allowing differentiation between reversible multilayer adsorption and irreversible pore filling. These observations have implications for assessing pore network accessibility, as narrow hysteresis indicates reversible filling suitable for reversible storage applications, whereas wide loops signal irreversible trapping that affects material performance in catalysis or filtration.

Nanoscale Corrections

At nanoscale dimensions, particularly for pore radii below 2 nm, the classical for capillary condensation breaks down due to the overlapping molecular potentials from opposing pore walls, which enhance adsorption beyond continuum predictions. This overlap creates a strong potential field gradient that favors fluid accumulation even at very low relative pressures, shifting the process from meniscus-driven condensation to cooperative filling. Additionally, disjoining pressure arises from wall-fluid interactions, such as , which stabilize thin adsorbed films and modify the effective pressure difference across the interface, leading to deviations from macroscopic capillarity assumptions. To address these limitations, corrections incorporate the curvature dependence of surface tension, parameterized by the Tolman length, which quantifies the shift between the equimolar surface and the surface of tension. For water, the Tolman length is approximately -0.07 nm (-0.71 Å), enabling the Kelvin equation to remain valid down to curvatures of about 0.5 nm when this dependence is included, as demonstrated in molecular simulations of condensation in carbon nanotubes. Furthermore, the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory integrates van der Waals attractions and electrostatic repulsions to model disjoining pressures in thin films, providing a framework for predicting enhanced adsorption in small pores where these long-range forces dominate over Laplace pressure. Advanced models like density functional theory (DFT) offer a molecular-level description for micropores, accounting for inhomogeneous density profiles and fluid-solid interactions that the Kelvin equation overlooks. Non-local DFT kernels accurately predict adsorption isotherms and phase transitions in pores as small as 1-2 nm, bridging simulations and macroscopic thermodynamics by incorporating disjoining effects and potential overlaps. Potential theory complements this by quantifying enhanced adsorption through wall potential gradients, revealing filling mechanisms in confined geometries. In materials like zeolites and metal-organic frameworks (MOFs) with sub-2 nm pores, these corrections manifest as condensation—or more precisely, volume filling—at relative pressures (P_v/P_sat) approaching zero, driven by overlapping field gradients that lower the free energy barrier for fluid ingress. For instance, in hydrophilic zeolites, initial adsorption occurs at low humidities due to strong site interactions, transitioning to pore filling without a distinct capillary step. Similarly, certain MOFs exhibit steep uptake at P_v/P_sat < 0.1, enabling applications in water harvesting. Recent advancements as of 2025 include refined models that describe capillary condensation features in sub-nanometer pores, accounting for extreme confinement and bridging classical thermodynamics with molecular simulations. A key limitation of these nanoscale corrections is the transition from the capillary condensation regime to a regime in ultrasmall pores (D < 4σ, where σ is the molecular diameter), where interfaces blur, and the process becomes a continuous, second-order density increase rather than a first-order phase jump. This crossover, governed by a pore-specific critical temperature, eliminates hysteresis and aligns with enhanced potential effects, as validated by simulations and zeolite experiments.

Associated Phenomena

Capillary Adhesion

Capillary adhesion refers to the attractive force between two solid surfaces mediated by liquid bridges formed via . This phenomenon occurs when vapor condenses in the narrow gap between surfaces, creating a meniscus with a curved interface that generates a pressure difference according to the . The resulting Laplace pressure deficit within the liquid bridge pulls the surfaces together, while the surface tension acts along the three-phase contact line to enhance the attraction. For small liquid bridges where the bridge height is comparable to the neck radius, the adhesion force can be approximated as F \approx 2\pi r \gamma \cos\theta, with r the radius of the bridge neck, \gamma the liquid-vapor surface tension, and \theta the equilibrium contact angle. Liquid bridges form through capillary condensation in gaps between particles or flat surfaces, particularly for hydrophilic materials exposed to ambient vapors like water. This condensation initiates at relative humidities around 50-70% RH for rough hydrophilic surfaces with small gaps, depending on gap size and surface wettability, as the vapor pressure required is reduced in confined geometries per the . Once formed, the bridge volume and shape evolve with humidity, leading to stronger adhesion as the meniscus curvature adjusts. The magnitude of the capillary adhesion force typically ranges from 10 to 100 nN for micron-scale bridges, scaling with bridge volume, separation distance (decreasing sharply beyond a few micrometers), and relative humidity. At fixed separation, higher humidity increases the effective bridge radius and reduces the mean curvature, amplifying the force until saturation near 100% RH. This force dominates van der Waals interactions in humid conditions and can exceed 10 times the dry adhesion for rough or deformable surfaces. In static conditions, the adhesion force represents the equilibrium pull-off force required to break the bridge at constant humidity, often exhibiting hysteresis due to pinning at the contact line. Dynamically, during surface separation, the force evolves through stages of neck elongation, instability, and rupture, with the maximum force occurring just before catastrophic breakup; rupture dynamics can generate transient peaks up to twice the static value. Even in ostensibly dry environments with trace moisture (e.g., <20% RH), capillary adhesion can cause unexpected sticking between closely apposed surfaces, as minute amounts of condensed vapor form nanoscale bridges that persist due to kinetic barriers to evaporation. This effect is particularly relevant in microscale devices and powder handling, where it contributes to stiction or agglomeration despite low ambient humidity.

Bridging Mechanisms

Capillary condensation facilitates the formation of liquid bridges between closely spaced particles or surface features, where vapor condenses into menisci that connect adjacent elements, often resulting in necklace-like chains of these menisci along particle alignments. These bridges arise when the relative humidity reaches levels sufficient to promote condensation in the narrow gaps, typically on the order of micrometers for micron-sized particles, leading to attractive forces that promote aggregation. The process is driven by the minimization of surface free energy, with the liquid phase wetting the particle surfaces and forming concave menisci that exert both surface tension and components. The strength and stability of bridges also depend on the contact angle, with better wetting (lower θ) promoting stronger forces. The stability of such capillary bridges is highly sensitive to the separation distance between particles, with stable configurations generally limited to critical separations of approximately 1-5 times the particle radius r, beyond which the bridge ruptures due to insufficient capillary pressure to sustain the meniscus shape. For separations closer than this threshold, the bridges exhibit robust connectivity, enabling the linking of multiple particles into extended structures. In multi-particle systems, such as loosely packed powders, these bridges induce coalescence by forming networks that enhance overall material cohesion, where the effective strength at each contact point scales with the number of individual bridges formed per particle pair, often increasing nonlinearly with liquid content. Environmental factors, particularly relative humidity, dictate the onset of bridging, with thresholds typically at high relative humidities of 70-95% RH or higher for micron-sized particles with smooth surfaces, and lower (down to ~50%) for rough particles or smaller sizes due to nanoscale crevices, depending on particle size, surface chemistry, and geometry; below these thresholds, insufficient vapor pressure prevents condensation. Upon decreasing humidity, these liquid bridges evaporate; however, in the presence of dissolved impurities, drying can leave solid residues that form permanent adhesive bonds persisting even at low humidity. Analytical models approximate the capillary force F between two spherical particles as F = 2\pi \gamma r^2 / d, where \gamma is the liquid-vapor surface tension and d is the center-to-center separation, capturing the dominant hydrostatic contribution for small bridges; more detailed derivations incorporate contact angle effects for precision. Numerical simulations of random particle packings reveal that bridge distribution follows the local void geometry, with forces varying spatially and leading to anisotropic stress fields that stabilize loose assemblies against shear. In granular materials, these mechanisms underpin cohesion under ambient humidity, enabling phenomena like self-supporting arches in sands or enhanced flow resistance in powders without full saturation.

Practical Applications

Pore Size Analysis

Capillary condensation serves as a foundational principle in pore size analysis, enabling the determination of pore size distributions in porous materials through the analysis of gas adsorption-desorption isotherms. The process relies on the Kelvin equation, which predicts the relative pressure at which vapor condenses in pores of varying radii due to enhanced capillary forces. This approach traces its origins to William Thomson (Lord Kelvin)'s 1871 derivation of the relationship between vapor pressure and the curvature of the liquid-vapor interface in capillaries. Over time, advancements have refined these concepts into standardized methods for characterizing porous structures, culminating in the International Union of Pure and Applied Chemistry (IUPAC) recommendations of 1985, which classify pores as micropores (width < 2 nm), mesopores (2–50 nm), and macropores (> 50 nm) based on adsorption behavior. The seminal Barrett-Joyner-Halenda (BJH) , developed in , applies the to the desorption branch of nitrogen adsorption isotherms to calculate pore size distributions, assuming cylindrical pores and stepwise multilayer adsorption followed by capillary evaporation. In this model, the volume of gas desorbed at each relative pressure corresponds to the emptying of pores with a specific radius, derived from the Kelvin radius adjusted for the thickness of the adsorbed multilayer film. The procedure begins with obtaining the adsorption-desorption isotherm, typically using at 77 K, followed by correction for multilayer thickness using statistical thickness curves such as the de Boer t-curve or those derived from Brunauer-Emmett-Teller (, which models multilayer adsorption on non-porous surfaces. The t-plot further refines this by plotting adsorbed against statistical thickness to isolate micropore filling from multilayer adsorption on external surfaces and mesopore walls. The BJH method excels in analyzing mesopores, providing reliable distributions for materials with pore sizes between 2 and 50 nm, but it has notable limitations for micropores, where enhanced adsorbate-adsorbent interactions and pore filling mechanisms deviate from classical capillary condensation assumptions, often leading to underestimation of pore volumes. For microporous materials, more advanced techniques like (DFT) are preferred, as they account for fluid-wall potentials and non-uniform density profiles without relying on macroscopic thermodynamics. observed in isotherms, arising from differences in adsorption and desorption mechanisms, is typically addressed by focusing on the desorption branch in BJH analysis to approximate equilibrium conditions. In practical applications, pore size analysis via capillary condensation is essential for characterizing catalysts and adsorbents, where mesopore distributions influence reactant diffusion and selectivity. For instance, BJH of reveals broad mesopore networks that enhance gas storage capacity for applications like CO2 capture, while in silica gels, it identifies uniform mesopore sizes around 2–10 nm that optimize moisture adsorption in drying processes. These insights guide , ensuring tailored for industrial performance.

Sintering and Materials Processing

In the sintering process, plays a crucial role by enabling the formation of liquid bridges from ambient vapor at the contact points or grain boundaries of powder particles, particularly under low-temperature conditions. This phenomenon occurs when the relative humidity exceeds the threshold dictated by the , causing water vapor to condense into menisci between particles, thereby generating attractive capillary forces that draw grains together and initiate neck formation. These forces promote and material transport, facilitating initial bonding without requiring high temperatures, as demonstrated in studies of pendular liquid bridges between spherical particles. In ceramics such as (YSZ) and (BaTiO₃), this mechanism enhances particle rearrangement and densification at temperatures below 200°C, often in conjunction with applied pressure. The sintering process unfolds in distinct stages influenced by these liquid bridges. During the initial stage, capillary forces from the condensed bridges act to compact particles and smooth asperities, increasing green body strength and enabling handling prior to full densification; this is particularly evident in humid environments where bridge volumes can reach up to 10% of particle separation distances for submicron grains. In the intermediate stage, the —governed by the Young-Laplace equation, ΔP = 2σ/r (where σ is and r is radius)—drives dissolution-precipitation mechanisms that fill interparticle voids, reducing and promoting for neck growth. Models of bridge geometry, such as the toroidal approximation, link bridge volume V to shrinkage rates, showing that smaller particle radii (e.g., 25 nm) yield higher forces (up to several nN) at 85% relative humidity and 85°C, accelerating densification by 10-20% compared to dry conditions. Examples include cold of , where controlled vapor exposure achieves 92-95% . This approach is applied under conditions of ambient or controlled (e.g., 50-85% relative humidity) and low pressures (350-520 ), making it suitable for energy-efficient processing of ceramics and metals like silver nanoparticles, where vapor-assisted bridges lower activation energies for mass transport. Benefits include enhanced green strength for improved handling of fragile compacts and reduced thermal stresses, enabling compatibility with temperature-sensitive additives or substrates in materials like porous ceramics. However, issues arise from uneven of the condensed liquid, potentially leading to localized stresses and cracking during subsequent or heating phases, which can be mitigated by gradual reduction. Quantitative assessments indicate that contributions can boost shrinkage rates by linking bridge volumes directly to reduction, with models predicting up to 15% volume shrinkage in early stages for humid-assisted processes.

Nanotechnology Devices

Capillary condensation poses significant challenges in the fabrication and operation of microelectromechanical systems () and similar nanoscale devices, primarily through the formation of liquid bridges that cause , where suspended structures adhere irreversibly to underlying surfaces. During wet etching processes, residual rinse liquids evaporate and form menisci in narrow gaps, generating attractive forces that exceed the restoring forces of microstructures, leading to collapse and losses. In humid environments, vapor-induced capillary bridges during operation can similarly pin , compromising functionality. Early surface-micromachined devices often experienced failure rates exceeding 50% due to this in-process stiction, highlighting its role as a dominant yield limiter in initial polysilicon-based designs. To mitigate these issues, dry release techniques such as supercritical CO2 drying have been widely adopted, as the eliminates during the transition from liquid to gas, preventing formation and enabling high-yield release of delicate structures like beams and gears. Additionally, applying hydrophobic coatings with contact angles greater than 90°, such as self-assembled monolayers of perfluorinated silanes on surfaces, reduces wettability and minimizes bridge strength, thereby enhancing device reliability in ambient . These strategies have improved fabrication yields to over 90% in modern processes. Beyond challenges, capillary condensation enables beneficial applications in devices, particularly for humidity-responsive actuation in sensors. In these designs, controlled formation and rupture of bridges between electrodes or compliant elements alter or mechanical deflection in response to relative changes, providing sensitive detection without external power. For instance, carbon nanotube-enhanced capacitive sensors exploit enhanced capillary condensation in nanopores to achieve rapid response times under 1 second for variations from 30% to 90% RH. In , capillary condensation in nanoporous electrodes enhances ion transport in supercapacitors by facilitating distribution and reducing resistance. In resource extraction, it influences recovery by altering permeability in nanopores through phase transitions that affect fluid flow. Emerging applications as of 2025 include atmospheric harvesting using amphiphilic nanoporous films that condense undersaturated vapor for sustainable collection, and controlled growth of arrays for photodetectors via capillary-driven on flexible substrates. Representative examples include microstructures in accelerometers and RF switches, where sub-micron gaps amplify capillary forces relative to restoring forces, following inverse scaling laws that make dominant at nanoscale separations—often orders of magnitude stronger than van der Waals interactions alone. These forces, akin to those in capillary adhesion phenomena, can reach several micro-Newtons in 100 nm gaps at 50% , underscoring the need for precise gap . In the 2020s, advances in integrating vapor barriers, such as atomic layer deposition (ALD) of alumina or parylene thin films, have further bolstered reliability by achieving water vapor transmission rates below 10^{-5} g/m²/day, effectively shielding internal structures from ambient moisture and reducing humidity-induced stiction in implantable and harsh-environment MEMS devices like pressure sensors.

Surface Characterization Techniques

Atomic force microscopy (AFM) leverages capillary condensation to probe nanoscale surface properties by measuring tip-sample interactions under controlled humidity conditions. In AFM setups, capillary forces arise from water menisci forming between the sharp tip (typically with a radius of 10-50 nm) and the sample surface due to ambient vapor condensation, dominating adhesion at relative humidities above 20-30%. These forces are quantified through cantilever deflection, where the tip approaches the surface, and the onset of condensation is detected as a sudden increase in attractive force when the tip-sample separation reaches a critical distance of approximately 9-10 nm. Force-distance curves, obtained by ramping the tip toward and away from the surface at speeds of 25-325 nm/s, reveal this condensation event as an abrupt jump in force, allowing precise mapping of local surface energetics without physical contact instability. This approach provides insights into surface chemistry and topography, as the meniscus force scales with the tip-sample geometry and vapor pressure. The resolution of AFM in detecting capillary bridges enables characterization down to sub-nanometer scales, with vertical sensitivity below 0.1 nm and lateral approaching 1 nm for bridge formation. Studies have visualized and quantified bridges as small as 0.5-1 nm in radius, particularly in confined geometries, by analyzing force jumps and retraction in the curves. Wettability mapping is achieved by correlating force variations with θ, where higher forces indicate more hydrophilic regions (θ < 30°) and lower forces hydrophobic ones (θ > 90°); for instance, on natural fibers, treated surfaces show up to 2-3 times greater due to enhanced . This technique maps heterogeneous wettability across rough or fibrous substrates, resolving θ variations over areas of 10-100 nm², which is crucial for understanding or interfaces. Beyond AFM, environmental scanning electron (ESEM) facilitates in-situ observation of capillary bridge dynamics by maintaining vapor pressures up to 100% relative within the chamber, allowing real-time of meniscus formation and evaporation on surfaces. In ESEM, bridges between particles or features are visualized at magnifications up to 100,000×, revealing growth and stability on superhydrophobic substrates where bridges persist at separations of 1-10 μm. (QCM) complements these by monitoring adsorption of water vapor leading to capillary condensation, measuring mass uptake on thin films or layers with picogram sensitivity. For SiO₂ , QCM detects rapid sorption-desorption cycles, with frequency shifts indicating bridge formation times on the order of seconds under steps from 10% to 80%. Combined QCM-ESEM setups have confirmed that capillary bridges enhance particle , with following Langmuir-type isotherms modified for multilayer adsorption. Data analysis in these techniques often involves fitting experimental force or mass data to the Kelvin equation, which relates the condensation relative humidity to the local meniscus radius r via \ln(RH) = -\frac{2\gamma V_m}{r RT}, where \gamma is surface tension, V_m molar volume, R gas constant, and T temperature; this yields local r values of 0.5-5 nm and \gamma variations up to 20% on heterogeneous surfaces. Humidity control, achieved through sealed chambers with saturated salt solutions or gas flow regulators maintaining RH to ±1%, isolates capillary contributions by minimizing van der Waals forces, which dominate below 10% RH and contribute <10% of total adhesion at higher levels. For example, at 50% RH, capillary forces account for 80-90% of pull-off adhesion on silica, separable via RH-dependent extrapolation. Recent advances in the have pushed AFM toward single-molecule sensitivity for studies on heterogeneous surfaces, incorporating qPlus sensors and low-temperature to achieve sub-Ångstrom force resolution and atomic-scale imaging of interfacial layers. In 3D nanomenisci (80-250 nm³), these setups map structures with signal-to-noise ratios exceeding 100, detecting individual molecule rearrangements during condensation on insulators like . Such capabilities reveal curvature-dependent \gamma down to molecular scales, enabling characterization of biomolecular adsorption sites where forces modulate binding energies by 10-50 .

References

  1. [1]
    Capillary Condensation - an overview | ScienceDirect Topics
    Capillary condensation is defined as the condensation of vapor into small pores or crevices within a porous solid at partial pressures above the vapor ...
  2. [2]
    [PDF] Capillary Condensation in Confined Media - HAL
    Oct 19, 2009 · From the above equation it is clear that capillary condensation can occur only if the liquid wets, at least partially, the solid surfaces.
  3. [3]
  4. [4]
    Historical aspects of capillarity and capillary condensation
    However, this exponential relationship was not the form of the original equation proposed in 1871 by William Thomson (later Lord Kelvin). Zsigmondy (1911) ...
  5. [5]
    On the Universality of Capillary Condensation and Adsorption ...
    Apr 30, 2020 · The underlying mechanism of capillary condensation and hysteresis in a mesoporous metal-organic framework (MOF), MIL-101(Cr) is explored.
  6. [6]
  7. [7]
    [PDF] Particle-Based Porous Material Development for Capillary ...
    Jan 19, 2023 · One of the significant applications of capillary condensation is atmospheric water capture. In the face of rapidly increasing consumption ...
  8. [8]
    Capillary condensation-driven growth of perovskite nanowire arrays ...
    Jan 24, 2025 · We develop a facile method utilizing capillary condensation to grow high-quality centimeter-scale perovskite nanowire arrays.
  9. [9]
    [PDF] Capillary condensation as a morphological transition
    The process of capillary condensationrevaporation in cylindrical pores is considered within the idea of symmetry breaking.
  10. [10]
    LX. On the equilibrium of vapour at a curved surface of liquid
    LX. On the equilibrium of vapour at a curved surface of liquid · Sir William Thomson F.R.S. · Published 1 December 1871 · Physics · Philosophical Magazine Series 1.
  11. [11]
    [PDF] Capillary Condensation in Confined Media - arXiv
    Oct 24, 2009 · We review here the physics of capillary condensation of liquids in confined media, with a special regard to the application in ...
  12. [12]
  13. [13]
    Revisiting Kelvin equation and Peng–Robinson equation of state for ...
    Mar 22, 2021 · Per the Kelvin equation, the alteration of vapour pressure in confined pores depends on capillary size, surface tension, and contact angle ...
  14. [14]
    LX. On the equilibrium of vapour at a curved surface of liquid
    On the equilibrium of vapour at a curved surface of liquid. Sir William Thomson FRS. Pages 448-452 | Published online: 13 May 2009.
  15. [15]
    [PDF] Capillary Condensation in Confined Media
    We see from the Laplace–Kelvin equation that the pres- sure in capillary condensates is usually very low: taking the example of water in ambient conditions ...
  16. [16]
    [PDF] 2015-IUPAC-Technical-Report.pdf - 3P Instruments
    (ii) pores of widths between 2 nm and 50 nm are called mesopores;. (iii) pores with widths not exceeding about 2 nm are called micropores. These limits ...
  17. [17]
    Accurate Relations between Pore Size and the Pressure of Capillary ...
    We study the size and shape dependence on capillary condensation within nanopores. ... We propose an alternative expression of the Kelvin equation to describe the ...<|control11|><|separator|>
  18. [18]
    On the capillary condensation of water in mesoporous silica films ...
    Nov 21, 2005 · The pressure of condensation is directly related to the liquid-vapor interfacial tension and to the pore geometry, its size, and shape. In ...
  19. [19]
    [PDF] A Theory of Capillary Condensation To Model Water Sorption within ...
    Sep 26, 2024 · (5) Can the capillary condensation phase transition be nucleated closer to the equilibrium chemical potential by perturbation of an external ...
  20. [20]
    Capillary Condensation - an overview | ScienceDirect Topics
    This is the so-called Kelvin equation, which was derived originally by Lord Kelvin for vapor pressure of small droplet without the negative sign. View chapter ...<|control11|><|separator|>
  21. [21]
    The role of contact angle and pore width on pore condensation and ...
    Aug 12, 2020 · Our results reveal that the water contact angle plays an important role in predicting the humidity required for pore filling.
  22. [22]
    Review of Kelvin's Equation and Its Modification in Characterization ...
    Apr 27, 2006 · Kelvin's equation is a link between the relative adsorbate pressure, the mean pore radius, and pore capillarity on the basis of macroscopic ...<|separator|>
  23. [23]
    Effect of Capillary Condensation on Friction Force and Adhesion
    The effect of the capillary condensation will be to partially solubilize and hydrate the polymeric material at the fiber joints (for example, strength additives ...Missing: importance | Show results with:importance
  24. [24]
    A review of wetting versus adsorption, complexions, and related ...
    Jun 27, 2013 · To avoid confusion, we will use the terminology “complete” or “perfect” wetting when the contact angle is zero, and “non-wetting” when the ...
  25. [25]
  26. [26]
  27. [27]
    Domain theory for capillary condensation hysteresis | Phys. Rev. B
    Jul 25, 2005 · These predictions over the scanning behavior within capillary condensation hysteresis can be used to check whether real materials are made up of ...
  28. [28]
    [PDF] Equilibrium, Pore Blocking, and Cavitation
    A qualitative explanation of the capillary hysteresis in the ink-bottle pore is as follows. The pore filling in the process of condensation follows the.
  29. [29]
    Understanding adsorption and desorption processes in mesoporous ...
    Sep 22, 2009 · For the case of the type A channel, where the ... In the type B channels with wider segments, the evaporation from the smallest pore ...
  30. [30]
    Cavitation in Metastable Fluids Confined to Linear Mesopores
    In the Si/B/A configuration, the capillary condensation is shifted toward higher pressure and cavitation occurs in layer B at a pressure slightly lower than ...
  31. [31]
    Capillary condensation of adsorbates in porous materials
    The shapes of the adsorption hysteresis loop classified by IUPAC (see Fig. 2 [1]) have been attributed to particular pore shapes. Type H1 loop, which is ...
  32. [32]
    Reminiscent capillarity in subnanopores | Nature Communications
    Oct 11, 2019 · For such blurred interfaces, pore filling becomes reversible and continuous (second order transition) in stark contrast with the capillary- ...
  33. [33]
    None
    Nothing is retrieved...<|separator|>
  34. [34]
    Adsorption and capillary condensation in porous media as a ...
    Feb 24, 2014 · The chemical potential of water may play an important role in adsorption and capillary condensation of water under multiphase conditions at ...
  35. [35]
    [PDF] Density Functional Theory of Capillary Condensation in Nanopores
    Non-local density functional theory (NLDFT) bridges the scale gap from molecular simulations to macroscopic thermodynamics, using inhomogeneous phases.
  36. [36]
    Water Adsorption in Porous Metal–Organic Frameworks and ...
    First, the pore filling or condensation of water into the pores of the solid must occur at low relative pressure (relative humidity) and exhibit a steep uptake ...
  37. [37]
    Article Water Harvesting with Zeolite, MOF, COF, and HOF Adsorbents
    Oct 23, 2025 · The strong hydrophilic sites of zeolite materials promote initial adsorption at low relative pressures, accompanied by a high adsorption ...
  38. [38]
    An analytical theory for the capillary bridge force between spheres
    Dec 22, 2016 · An analytical theory has been developed for properties of a steady, axisymmetric liquid–gas capillary bridge that is present between two ...
  39. [39]
    Capillary Forces in Nanoparticle Adhesion: A Review of AFM Methods
    Aug 24, 2015 · Capillary Forces in Nanoparticle Adhesion: A Review of AFM Methods. Aaron J. HarrisonSchool of Chemical Engineering, Purdue University, West ...
  40. [40]
    Formation and rupture of capillary bridges in atomic scale friction
    Oct 31, 2012 · It has been found that condensation of capillary bridges may be responsible for frictional ageing that is manifested as logarithmic increase of ...
  41. [41]
    Simulation of Capillary Bridges between Particles - ScienceDirect.com
    The dependence of the capillary force on various parameters such as particle sizes, contact angle and humidity was investigated and it is demonstrated that ...
  42. [42]
    Capillary suspensions: Particle networks formed through the ...
    The addition of small amounts of a secondary fluid to a suspension can, through the attractive capillary force, lead to particle bridging and network formation.
  43. [43]
    Effects of relative humidity and particle and surface properties on ...
    In particular, high RH can lead to capillary condensation between a particle and a surface strongly enhancing adhesion and suppressing particle resuspension ( ...
  44. [44]
    The strength of liquid bridges in random granular materials
    The strength of capillary bridges in randomly packed granular media was analyzed by means of computer simulations. A novel simulation method, based on the ...
  45. [45]
    Historical aspects of capillarity and capillary condensation - ADS
    The related capillary condensation equation is now universally known as the Kelvin ... original equation proposed in 1871 by William Thomson (later Lord Kelvin).
  46. [46]
    [PDF] REPORTING PHYSISORPTION DATA FOR with Special Reference ...
    pore size and shape. The upper limit of 2.0 nm for the micropore width was put forward as part of the IUPAC classification of pore size (see ...
  47. [47]
    The Determination of Pore Volume and Area Distributions in Porous ...
    The determination of pore volume and area distributions in porous substances. I. Computations from nitrogen isotherms.
  48. [48]
    [PDF] Brunauer-Emmett-Teller (BET) surface area analysis
    In multilayer adsorption the adsorp- tion space accommodates more than one layer of molecules so that not all the adsorbed molecules are in direct contact with ...
  49. [49]
    Surface Analysis by T-plot - Materials Characterization Lab
    In the t-plot method, a t-curve is plotted which graphs the average thickness of the adsorbed layer on the surface. into the average thickness of the film ...Missing: multilayer | Show results with:multilayer
  50. [50]
    Density functional theory methods for characterization of porous ...
    This review presents the state-of-the-art of adsorption characterization of mesoporous and microporous materials by using the density functional theory (DFT) ...Density Functional Theory... · Introduction · Applications Of Dft Kernels...
  51. [51]
    CO2 Capture by Low-Cost Date Pits-Based Activated Carbon ... - PMC
    Jul 12, 2021 · CO2 capture uses date pits-based activated carbon and silica gel. The activated carbon achieved 73.1 mg/g uptake, while silica gel achieved 8. ...
  52. [52]
    Evaluation of Silica gel adsorbent potential for carbon dioxide capture
    Aug 30, 2022 · In this research, silica gel as a low-cost adsorbent for carbon dioxide uptake was investigated experimentally.
  53. [53]
    A study of pendular liquid bridge between two equal solid spheres
    ### Summary of Liquid Bridge Formation in Sintering (arXiv:2501.08344)
  54. [54]
    Cold sintering of microwave dielectric ceramics and devices
    Jan 13, 2021 · Cold sintering is able to densify ceramics at < 200 °C via a combination of external pressure and a transient liquid phase, reducing the energy consumed.
  55. [55]
    Comparison of different sintering aids in cold sinter‐assisted ...
    May 24, 2021 · Floyd et al. concluded that increasing the relative humidity from 24% to 54% at room temperature leads to a 10% rise in the relative density ...
  56. [56]
    Direct measurement of the capillary condensation time of a water ...
    Sep 14, 2018 · We directly measured the capillary condensation time of a water meniscus, by approaching two surfaces at different speeds and monitoring the relative position ...
  57. [57]
    Adhesion Forces Measured by Atomic Force Microscopy in Humid Air
    The results indicate that even under the condition where capillary condensation occurs there is chemical specificity in the measured pull-off force. Issues such ...
  58. [58]
    In Situ Atomic-Scale Imaging of Interfacial Water under 3D ...
    May 13, 2021 · We develop an AFM-based method to provide in situ atomic-scale resolution maps of the solid–water interface of a nanomeniscus (80–250 nm 3 ).
  59. [59]
    Atomic force microscopy characterization of the surface wettability of ...
    When AFM is conducted under atmospheric conditions, a liquid bridge may form between the tip and sample due to capillary condensation, resulting in a meniscus ...
  60. [60]
    In Situ Observation of Liquid Capillary Bridges on Superhydrophobic ...
    We report a technique for observing the dynamic behavior of a liquid contact line at high magnification using environmental scanning electron microscopy. We ...
  61. [61]
    Study of water adsorption and capillary bridge formation for SiO 2 ...
    Feb 24, 2014 · Water adsorption and capillary bridge formation within a layer of SiO2-nanoparticles were studied in situ by means of a combination of quartz ...Experimental · Results · Afm, Ft-Ir, Xps And Fe-Sem...
  62. [62]
    Direct Evidence for Curvature-Dependent Surface Tension in ...
    Dec 12, 2018 · Capillary condensation is the first-order vapor-to-liquid phase transition taking place in confined geometries.
  63. [63]
    Effect of relative humidity on onset of capillary forces for rough ...
    Dec 1, 2013 · Atomic force microscopy (AFM) was used to investigate the effect of relative humidity (RH) on the adhesion forces between silicon nitride AFM probes, ...Missing: isolate | Show results with:isolate
  64. [64]
    [PDF] Capillary condensation under atomic-scale confinement - arXiv
    Their contact angles are known to be in the range of 0–20° and 55–. 85°, respectively16,28,29. For surfaces exposed to air, θ is close to the above upper ...