Capillary condensation
Capillary condensation is the process by which a vapor condenses into liquid within the narrow pores or crevices of a porous solid at partial pressures below the saturation vapor pressure of the bulk liquid, driven by the curvature of the liquid-vapor meniscus in confined geometries.[1] This phenomenon occurs preferentially in wetting systems where the contact angle is less than 90°, enabling the formation of stable liquid bridges or menisci that lower the equilibrium vapor pressure required for phase transition.[2] The underlying mechanism is thermodynamically described by the Kelvin equation, which relates the shift in vapor pressure to the pore radius: \ln(p/p_s) = -\frac{2\gamma V_m \cos \theta}{rRT}, where p is the equilibrium vapor pressure over the curved interface, p_s is the saturation vapor pressure, \gamma is the surface tension, V_m is the molar volume of the liquid, \theta is the contact angle, r is the pore radius, R is the gas constant, and T is the temperature.[3] Building upon foundational contributions from Thomas Young (1805) on surface tension and Pierre-Simon Laplace (1806) on capillary forces, William Thomson (Lord Kelvin) 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 20th century, such as Richard Zsigmondy's 1911 studies on silica gels.[4] A hallmark of capillary condensation is the adsorption-desorption hysteresis 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.[5] Capillary condensation plays a pivotal role in numerous scientific and engineering fields, influencing processes from gas adsorption in catalysis and chromatography to water management in soils and textiles.[1] In energy applications, it affects shale gas recovery by altering permeability in nanopores and enhances ion transport in supercapacitors.[6] Emerging uses include atmospheric water harvesting via tailored porous structures and the controlled growth of nanomaterials, such as perovskite nanowires, underscoring its versatility in addressing challenges in sustainability and nanotechnology.[7][8]Fundamentals
Definition and Mechanism
Capillary condensation is the phenomenon in which a vapor phase transitions to a liquid phase 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 vapor pressure and P_sat is the saturation pressure of the bulk liquid.[2] This process is fundamentally driven by surface tension at the liquid-vapor interface and the associated Laplace pressure, which arises from the curvature of the interface in confined geometries.[2] The mechanism initiates with the adsorption of vapor molecules onto the pore walls, forming an initial thin liquid film due to attractive intermolecular forces, primarily van der Waals interactions.[2] As the relative vapor pressure rises, the film thickens until it reaches an instability point, where nucleation of discrete liquid clusters occurs, often described as a morphological transition from a uniform wetting layer to localized liquid phases.[9] These clusters develop into curved menisci, whose concave shape in wetting pores generates a capillary pressure that reduces the equilibrium vapor pressure needed for condensation, thereby promoting further liquid formation.[2] Key physical drivers include surface tension (γ), which dictates the meniscus curvature and pressure differential, and confinement effects that enhance intermolecular forces; in hydrophilic pores, strong wetting favors film stability and rapid filling, whereas hydrophobic pores exhibit delayed or altered condensation dynamics due to reduced adhesion.[2] The process unfolds step-by-step: vapor adsorption first establishes menisci at pore walls or entrances, creating localized regions of high curvature that draw additional vapor inward via capillary action.[2] This leads to the growth of liquid bridges connecting opposite pore surfaces, stabilized by the balance of surface tension and disjoining pressures from intermolecular forces.[9] As pressure increases, these bridges expand and coalesce, ultimately filling the pore volume with liquid condensate.[2] The quantitative foundation for this pressure reduction is captured by the Kelvin equation, relating vapor pressure to interface curvature.[10] 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.[10]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.[11] The Laplace pressure plays a central role, creating a pressure difference across the curved interface given by \Delta P = \frac{2\gamma \cos\theta}{r}, where \gamma is the liquid-vapor surface tension, \theta is the contact angle, and r is the pore radius. This pressure difference induces supersaturation within the pores, as the liquid phase experiences a lower pressure 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.[11] 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 condensation, leading to metastable states.[11] Unlike bulk condensation, which occurs at the saturation vapor pressure, confinement alters the phase diagram by the Laplace pressure and surface terms, enabling condensation at lower relative pressures and introducing hysteresis due to these energy barriers. This thermodynamic shift is particularly pronounced in nanopores, where surface effects dominate over bulk contributions.[12]Theoretical Framework
Kelvin Equation
The Kelvin equation 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 chemical potential balance between coexisting phases.[13] For the vapor phase, assumed to behave as an ideal gas, the chemical potential 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 chemical potential at temperature T, R is the universal gas constant, P_v is the equilibrium vapor pressure in the pore, and P_\mathrm{sat} is the saturation vapor pressure over a flat liquid-vapor interface.[13] For the liquid phase in the bulk at saturation, \mu_l^\mathrm{bulk} = \mu^0(T). In a capillary pore, the liquid experiences a pressure reduction due to the Laplace pressure across the curved meniscus at the liquid-vapor interface, given by \Delta P = \frac{2\gamma \cos\theta}{r} for a cylindrical geometry, where \gamma is the liquid-vapor surface tension, \theta is the contact angle, and r is the pore radius. The pressure in the liquid is thus P_l = P_v - \Delta P. Assuming the liquid is incompressible with constant molar volume V_m, the chemical potential of the confined liquid 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}). [13] Setting \mu_v = \mu_l at equilibrium 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 Kelvin equation:
\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 pressure 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 equilibrium to curved surfaces.[14][13] 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}. [13][15] 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}. [15] 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.[16][13]