
==== Front
Proc Natl Acad Sci U S A
Proc Natl Acad Sci U S A
PNAS
Proceedings of the National Academy of Sciences of the United States of America
0027-8424
1091-6490
National Academy of Sciences

39264743
202318386
10.1073/pnas.2318386121
research-articleResearch ArticlevideoVideoapp-physApplied Physical Sciences405
Physical Sciences
Applied Physical Sciences
Deformation dynamics of nanopores upon water imbibition
Sanchez Juan a 1
Dammann Lars a b c 1 https://orcid.org/0000-0001-5759-3896

Gallardo Laura a d https://orcid.org/0000-0002-5851-0105

Li Zhuoqing a b https://orcid.org/0009-0007-6367-5063

Fröba Michael e
Meißner Robert H. c f https://orcid.org/0000-0003-1926-114X

Stone Howard A. g https://orcid.org/0000-0002-9670-0639

Huber Patrick patrick.huber@tuhh.de
a b 2 https://orcid.org/0000-0002-2126-9100

aInstitute for Materials and X-ray Physics, Hamburg University of Technology, Hamburg 21073, Germany
bCenter for X-ray and Nano Science, Deutsches Elektronen-Synchrotron, Hamburg 22607, Germany
cInstitute of soft Matter Modeling, Hamburg University of Technology, Hamburg 21073, Germany
dCentre for the Study of Manuscript Cultures, Hamburg University, Hamburg 20354, Germany
eInstitute of Inorganic and Applied Chemistry, University of Hamburg, Hamburg 20146, Germany
fInstitute of Surface Science, Helmholtz-Zentrum Hereon, Geesthacht 21502, Germany
gDepartment of Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ 08540
2To whom correspondence may be addressed. Email: patrick.huber@tuhh.de.
Edited by Abraham D. Stroock, Cornell University, Ithaca, NY; received November 17, 2023; accepted July 26, 2024, by Editorial Board Member John A. Rogers

1J.S. and L.D. contributed equally to this work.

12 9 2024
17 9 2024
12 9 2024
121 38 e231838612117 11 2023
26 7 2024
Copyright © 2024 the Author(s). Published by PNAS.
2024
https://creativecommons.org/licenses/by/4.0/ This open access article is distributed under Creative Commons Attribution License 4.0 (CC BY).

Significance

Capillarity-driven flows in nanometer-sized pores play a dominant role in many natural and technological processes, ranging from water transport and transpiration in trees, clay swelling, and catalysis to transport through microfluidic structures and fabrication of battery materials. Here, we show by a combination of experiments and computer simulations of water imbibition in nanopores that the competition between expansive, surface stress release at pore walls and negative, contractile Laplace pressures of nanoscale menisci lead to an unusual macroscopic behavior of the porous medium, which is generic for any liquid/nanoporous solid combination. The results allow one to quantify surface and Laplace stresses and to monitor nanoscale flow and infiltration states by relatively simple length measurements of the porous medium.

Capillarity-driven transport in nanoporous solids is widespread in nature and crucial for modern liquid-infused engineering materials. During imbibition, curved menisci driven by high negative Laplace pressures exert an enormous contractile load on the porous matrix. Due to the challenge of simultaneously monitoring imbibition and deformation with high spatial resolution, the resulting coupling of solid elasticity to liquid capillarity has remained largely unexplored. Here, we study water imbibition in mesoporous silica using optical imaging, gravimetry, and high-resolution dilatometry. In contrast to an expected Laplace pressure-induced contraction, we find a square-root-of-time expansion and an additional abrupt length increase when the menisci reach the top surface. The final expansion is absent when we stop the imbibition front inside the porous medium in a dynamic imbibition-evaporation equilibrium, as is typical for transpiration-driven hydraulic transport in plants, especially in trees. These peculiar deformation behaviors are validated by single-nanopore molecular dynamics simulations and described by a continuum model that highlights the importance of expansive surface stresses at the pore walls (Bangham effect) and the buildup or release of contractile Laplace pressures as menisci collectively advance, arrest, or disappear. Our model suggests that these observations apply to any imbibition process in nanopores, regardless of the liquid/solid combination, and that the Laplace contribution upon imbibition is precisely half that of vapor sorption, due to the linear pressure drop associated with viscous flow. Thus, simple deformation measurements can be used to quantify surface stresses and Laplace pressures or transport in a wide variety of natural and artificial porous media.

nanoporous materials
water
nanofluidics
capillarity
strain
Deutsche Forschungsgemeinschaft (DFG) 501100001659 CRC 986 Juan SanchezLars DammannLaura GallardoPatrick Huber Deutsche Forschungsgemeinschaft (DFG) 501100001659 TRG 2462 Juan SanchezLars DammannLaura GallardoPatrick Huber Deutsche Forschungsgemeinschaft (DFG) 501100001659 SPP 2171 Juan SanchezLars DammannLaura GallardoPatrick Huber Deutsche Forschungsgemeinschaft (DFG) 501100001659 EXC 2176 Juan SanchezLars DammannLaura GallardoPatrick Huber
==== Body
pmcLiquid-infused nanoporous solids play an increasingly important role as functional materials. One can find them as soft-hard hybrids with superior properties in terms of adaptive wettability (1), mechanical actuation (2, 3), adjustable photonics (4, 5), in nanofluidics (6), and in energy storage and energy conversion (7) applications. Also, the development of structural materials that are capable of adapting to changing environmental conditions can be achieved by combining soft, dynamic liquid phases with static, solid phases that act as mechanically robust scaffold structures. The resulting hybrid materials have demonstrated unprecedented properties of stability, adaptability, and stimuli-responsiveness, as desired in many applications such as 3D printing, soft robotics, omniphobic surfaces, microfluidics, and multiphase separation (1, 8, 9).

It has also long been known, starting with the pioneering work of Meehan (10) and Bangham and Fakhoury (11), that stresses and strains are generated in porous solids during the adsorption and desorption of gases and liquids (12–22). Two effects play a major role in the resulting sorption-induced deformation (18, 23). First, the adsorption of matter onto the solid pore wall surface typically leads to a reduction in surface stress. The reduced stress results in an expansion of the porous solid (Bangham effect) (24). Second, the formation of menisci in the pore space during capillary condensation can create a negative pressure in the liquid, the magnitude of which is inversely proportional to the pore size. The effect of the negative pressure on the solid surface causes the solid to contract (Laplace effect) (25).

Similar to fluid adsorption processes, deformation of porous solids is also expected to occur in capillarity-driven spontaneous imbibition experiments with wetting liquids (12, 26, 27) as well as in forced intrusion of nonwetting liquids (28, 29). In fact, studies highlight the importance of Laplace pressure-induced pore space deformations and their impact on imbibition dynamics for soft (macro)porous media such as cellulose sponges (30–32). However, the influence of surface stress on the deformation behavior, which is particularly important for nanoporous media with their large internal surface-to-volume ratio and which is predicted to lead to pore space expansion, has been little studied. Thus, in contrast to gas adsorption-induced deformations, a comprehensive understanding of imbibition-induced deformations as a result of the interplay of Bangham and Laplace effects is still lacking for the important class of liquid-infused nanoporous materials.

Here, we aim to achieve such a mechanistic understanding of imbibition-induced deformation in nanoporous solids by performing combined experiments on water imbibition and the corresponding deformation dynamics of a model nanoporous solid (monolithic Vycor glass). Furthermore, by comparison with atomistic simulations of water imbibition on a single silica nanopore, we trace the distinct deformation regimes observed at the macroscopic scale to the single-pore behavior.

Results and Discussion

Imbibition Experiments with Water: The Interplay of Front Movement, Front Broadening, and Deformation Dynamics.

We study the deformation of nanoporous Vycor glass upon imbibition of liquid water. Vycor is an almost pure fused silica glass permeated by a 3D network of interconnected pores (33, 34). It is formed by a leaching process after spinodal decomposition of a borosilicate glass. Therefore, its geometric structure can be well represented by clipped Gaussian random fields (35, 36).

To monitor the capillarity-driven spontaneous water imbibition, we perform optical imaging experiments. They allow us to determine the mean front position and its broadening, respectively, due to a transparency change induced by the capillary filling of the pore space and light scattering typical of the inhomogeneously filled pores near the imbibition front; see Movie S1 (37–40). A set of selected snapshots is shown in Fig. 1A aligned with the optically determined mean imbibition front position (Fig. 1B), mass uptake of the porous monolith (Fig. 1C), and the measured sample strain as a function of time (Fig. 1D).

Fig. 1. Elastocapillarity dynamics upon water imbibition in nanoporous Vycor glass on the porous-medium scale. (A) Optical snapshots of water imbibition in a rod-like Vycor monolith. The time between snapshots is 15 min. (B) Mean front position, (C) mass of water uptake, and (D) sample strain measured as a function of time by optical imaging, gravimetry, and dilatometry, respectively. (E) Schematic illustration of the macroscopic deformation of the glass monolith as a function of sample height (Left), caused by the interplay of sample contraction, negative hydrostatic pressure (purple), and sample expansion due to surface stress release at the water–silica wall interface per unit sample length (red) as a function of sample height for one selected time t during imbibition. The strain sensor or pushrod is represented as a cylindrical glass piece in contact with the Top face of the sample, while the Bottom face is in contact with a liquid reservoir. (F) Schematic deformation states for selected times during imbibition (regime I) and for complete filling (regime III); also indicated are the resulting changes in the sample length (green line) as well as the imbibition front position (blue line).

Since the relative humidity (RH) and thus potential adsorption of water via the vapor phase before and during the infiltration experiments is expected to affect both the deformation and imbibition dynamics (37, 41), we perform the experiments in a closed cell with controlled humidity. In the following, we first discuss the experiment with a RH of 50%. According to a water vapor sorption isotherm, this corresponds to approximately two preadsorbed monolayers of water molecules at the silica pore walls (37).

Evaporation from the porous block was prevented by shielding the side surface of the Vycor monolith, leaving the top and bottom surfaces unshielded to allow unaltered contact with the dilatometer pushrod and liquid water inlet, respectively. Notice that this configuration allows for air escape through the top facet. Therefore we do not expect any air-entrapment-induced impact on the imbibition and deformation dynamics as has been observed for dead-end single nanopores (42). After initiating capillary rise in the porous glass by raising the level of the bulk liquid reservoir until it reaches the bottom of the glass cylinder, we observe deformation kinetics that include three distinct regimes. Regime I lasts from the moment of contact with the bulk liquid reservoir until the imbibition front reaches the top of the matrix; see Fig. 1D. Here, the sample expands in accordance with the mean front motion according to a square-root-of-time dynamics. This regime is terminated by a relatively abrupt, additional expansion of the matrix that occurs when the advancing part of the white imbibition front (Fig. 1A) reaches the top of the matrix. This step-like deformation continues until the receding part of the white front reaches the top surface, and thus until complete filling of the pore space is achieved. This can be seen in both the optical imaging (Fig. 1B) and mass uptake (Fig. 1C) experiments. In agreement with a constant final mass uptake and the complete disappearance of the white front, the subsequent regime III after this “strain jump” indicates the complete filling of the pore space and is characterized by a deformation plateau. The total strain measured after complete filling of the matrix is 0.09%, which corresponds to a length increase of 10.79 μm of our sample (original length 11.99 mm).

Next, we propose a model for the time-dependent strain evolution. The model requires a macroscopic description of the imbibition dynamics, as well as a characterization of the strain-inducing elastocapillary phenomena occurring at the single-nanopore scale and how they translate to the macroscopic porous-medium scale.

Elastocapillarity Dynamics: From the Single-Nanopore to the Porous-Medium Scale.

The imbibition dynamics of the Lucas-Washburn front motion can be derived from a Darcy description of the spontaneous imbibition process (37, 43). Since at these scales inertial and gravitational effects are negligible compared to capillarity, the imbibition kinetics are governed by a dynamic equilibrium between the approximately constant Laplace pressure at the highly curved menisci of the imbibition front, which drives the flow, and the increasing viscous drag in the liquid-saturated part of the sample, which opposes the flow. As outlined in ref. 37, for pore-space geometries such as those present in Vycor glass, the following equation is obtained for the imbibition height h(t) as a function of time t:[1] h(t)=γ2ητξ(rp−d)3rp2t,

where rp is the pore radius, γ is the water-vapor surface tension, η is the viscosity of water, and τ is the hydraulic tortuosity. The preadsorbed water film on the pore walls has a thickness d. This pore wall coverage results in a reduced, initial porosity Φi accessible to the imbibing water, compared to the original host material porosity Φ0, measured in vacuum. Therefore, a scaling factor ξ=ΦiΦ0 must be included in Eq. 1. We will refer to this equation as the modified Lucas Washburn equation (MLW).

Eq. 1 can be extended to describe the mass uptake of water by multiplying the host material area A, the reduced porosity Φi, and the water density ρ, giving Δm(t)=ρΦiAh(t). Assuming two immobile water layers (d= 0.5 nm), one obtains an excellent fit to the observed mean imbibition position and mass uptake; see Fig. 1 B and C. This reduction in the hydraulic radius can be attributed to the strongly bound, highly viscous interfacial water layers at the silica surfaces, in agreement with previous experiments and computer simulations (37, 44). Note, that the preadsorbed water layers lead to an initial Bangham effect, i.e., adsorption-induced expansion of the Vycor glass relative to the completely dry sample, even before the imbibition experiment is started.

Upon contact with the bulk reservoir, two processes occur in parallel. In the part of the pore space filled with liquid, i.e., the part of the sample up to the capillary rise height h(t), the Bangham expansion takes place due to the corresponding release of surface stress. In addition, the pressure in the liquid increases linearly from the tensile Laplace pressure of the menisci in the advancing imbibition front to the atmospheric pressure at the bulk water reservoir. This results in a contractile stress contribution that decreases linearly from the imbibition front to the bulk reservoir. Thus, the total deformation of the matrix must be considered as a superposition of two opposing deformation effects, the expansive surface stress release (Bangham effect) and the contractile Laplace pressure effect; see Fig. 1E.

Since the induced strain and thus pore size changes are small, on the order of fractions of a percent, the resulting time-dependent changes in the hydraulic permeabilities and thus capillarity-driven flow dynamics due to the changes in the pore-sizes are negligible. Therefore, both deformation effects scale, to good approximation, linearly with the capillary-filling dynamics. More precisely, the strain dynamics follow the classical square-root-of-time LW dynamics, as proven by our experimental strain observation in regime I; see Fig. 1D.

Since we observe the matrix to be expanding, the expanding Bangham contribution obviously dominates the Laplace pressure effect. The abrupt expansion observed in regime II is explained by the sudden vanishing of the contractile Laplace pressure contribution. Upon reaching the top surface of the matrix the highly curved menisci and the corresponding Laplace pressure contribution vanish. The disappearance of the contractile contribution leads to an additional sample expansion, which we refer to as the “Laplace expansion jump.” The duration of the observed expansion is dictated by the broadening of the imbibition front upon complete filling of the sample and thus by the time characteristic of its vanishing, when the front reaches the top of the sample. This can be seen by a comparison of the vanishing of the white scattering front in the optical images with the times characteristic of the beginning and end of deformation regime II, respectively (Fig. 1).

Water vapor transport through the pore network is very slow compared to capillarity-driven liquid flow, so the effects of evaporation through the top surface during imbibition are negligible. Additionally, the upper dilatometer glass pushrod effectively sealed the upper sample surface. In the absence of a mechanism to counteract water uptake, the liquid will continue to rise according to MLW dynamics until the pore space is completely filled (37, 45). In Fig. 2, Top we show such dynamics and qualitatively illustrate the flattening of the menisci once the liquid reaches the upper porous block surface and thus the upper monolithic glass pushrod.

Fig. 2. Impact of evaporation on imbibition-induced deformation. Imbibition-induced deformation of a nanoporous Vycor glass monolith (rp= 3.4 nm and RH≈ 5%) measured by in situ dilatometry in two different scenarios: the first without evaporation (NEv) in (A and B) and the second with evaporation (Ev) in (C and D). In both cases, we present our dilatometry data (NEv: green circles; EV: red circles), combined with optical imaging. We select a vertical pixel line centered in the sample and integrate over time. The visible white region reveals the broadening of the imbibition front separating the imbibed and dry portions of the porous host. The dashed lines in (A and C) correspond to fitting curves based on a t scaling. The fitting curves in (C) do not start at t=0 because they do not account for the early times when unavoidable lateral imbibition resulted from the insertion of the Bottom of the unsealed into the liquid reservoir. The deformation dynamics are compatible with the expected exponential relaxation of the imbibition process, supported by the blue solid fitting line (see SI Appendix for derivation). In the nonevaporation scenario, such a condition is achieved by sealing the monolith lateral surface as shown in panel (B). Here, we show a representative 2-D slice of the completely filled porous matrix in contact with a liquid reservoir. The Top of the sample is effectively sealed by the dilatometer pushrod. The orange arrows represent the water flow from the reservoir to the pore space (imbibition). Zooming into the Upper region, we see the absence of curved menisci since there is no evaporating outflow to counterbalance the capillarity-driven imbibition. This contrasts with the presence of menisci at the saturated–unsaturated interface, as shown in panel (D). Here, we illustrate the dynamic equilibrium reached for an unsealed sample when both the evaporation and imbibition flow rates are equal. The orange arrows at the curved water–gas interface symbolize the evaporation of water molecules.

Impact of Evaporation on Elastocapillarity Dynamics: The Extreme Case of Reaching an Evaporation-Imbibition Equilibrium.

We performed a second imbibition experiment in which, unlike the first, the side facets are not sealed against evaporation. Then, the evaporation fluxes increase in proportion to the wetted, vapor-phase exposed part of the side facets and thus in proportion to the capillary rise height. The imbibition front then stops at a height well below the upper monolith surface as the capillary rise reaches dynamic equilibrium; see Movie S2. The inflow from the bulk liquid reservoir equals the evaporative outflow to the bulk vapor phase, so that an evaporation-driven water transport reminiscent of an (artificial) tree is achieved (46–48); see Fig. 2A.

In this case, the dynamics of capillary rise deviates from the classical t scaling. As shown in Fig. 2, Bottom and Left, the water reaches a final capillary rise height, where the arrested imbibition front consists of highly curved menisci. From the imbibition front to the water reservoir there exists a linear pressure gradient that continues to exert a tensile contraction on the wetted porous material. Therefore, no “Laplace expansion jump” can be observed in the sample even though the imbibition front stopped rising.

As outlined in SI Appendix, the capillary-rise dynamics for this scenario can be described by[2] h(t)=heq·1−e−4qRΦit,

where q is the evaporation rate of water per surface unit and heq=KΔpR2qμ represents the equilibrium capillary rise or saturation height, reached at sufficiently long times. The hydraulic permeability K and every variable involved in heq except for q can be obtained for the case of water imbibition in Vycor glass from refs. 37 and 45.

Since the deformation dynamics evolve proportionally to the capillary rise, the inclusion of a proportionality factor C allows the use of Eq. 2 as a fitting function, ϵ(t)=C·1−e−βt, for the dilatometry data. Our experiment leads to β=2.1×10−4 s−1, which translates to an evaporation rate q= 0.05 μm s−1. This value is in reasonable agreement with values reported in literature q≈ 0.045 μm s−1 for bulk deionized water at similar conditions of humidity and temperature (49). The determined evaporation rate q leads to an equilibrium capillary rise of heq= 14 mm or a filling height of 74% of the total length (19 mm), which is in reasonable agreement with the observed averaged position of the arrested front visible in the snapshots in Fig. 2C.

Note that a first-order Taylor series expansion of Eq. 2 around t=0 leads back to our MLW law or Eq. 1. This means that for early times (t≪tE=Rφi4q), where tE=4700s±600s corresponds to the time when the evaporative flow rate has reached 1 − 1/e = 63% of the imbibition flow rate, imbibition dominates over evaporation and the rising dynamics can be described exclusively in terms of the classical LW square-root-of-time dependent law (SI Appendix). This is evidenced by the overlapping of the fitting curves in Fig. 2C, for t<tE, where the continuous line indicates the fitting model based on Eq. 2 and the dashed line denotes the t-fitting based on Eq. 1. The latter is in reasonable agreement with the deformation dynamics shown in Fig. 2A, corresponding to a nonevaporation scenario (sealed sample), and is very well described by the MLW model. The fitting coefficients are Cunsealed= (8.53 ± 0.50) ×10−4 s−0.5 and Csealed= (10.4 ± 0.2) ×10−4 s−0.5.

In the derivation of Eq. 2, we consider only the liquid flow in the vertical imbibition direction, so that a one-dimensional (1-D) problem results for h=h(t) and no transversal imbibition front height variation is considered. However, a closer look at the shape of the imbibition front reveals that it becomes increasingly curved with increasing imbibition time, with h being maximal in the center of the sample and decreasing symmetrically in the radial direction toward the sample sides, as evident in the snapshot sequence in SI Appendix, Fig. S1. This suggests an increasing importance of the transversal flow components in the imbibition phenomenology with increasing imbibition time. In fact, the slight but systematic overestimation of the strain by our simple 1-D model compared to the measured strain in the intermediate time range from about 5,000 s to the onset of saturation at about 15,000 s remains unexplained. Perhaps it results from neglecting the curvature of the imbibition front and thus the lateral flow components, since the average height of the curved front and thus the averaged Bangham effect is smaller than expected in the simplified model with a flat imbibition front. Furthermore, we neglect the Kelvin effect, i.e. a reduction of the vapor pressure at the highly curved menisci, which could also lead to an underestimation of the equilibrium height heq, since it results in reduced evaporative fluxes compared to the planar bulk interface.

Quantitative Modeling of the Deformation Dynamics on the Porous-Medium Scale.

To achieve a quantitative description of the deformation kinetics we consider the monolithic porous glass as a 3-D network of interconnected cylindrical channels with an isotropic orientation distribution. We neglect the broadening of the imbibition front and therefore assume that the imbibed part of the network is completely filled with liquid up to the height h(t).

As outlined in detail in ref. 24, given the narrow pore-size distribution of Vycor glass and the disordered channel-like structure of the pores, its deformation can be considered uniform, isotropic, and a single pore-size approximation can be used. The deformations are determined by a surface stress contribution and by the hydrostatic fluid pressure in the pore pF, both acting as an effective loading pressure on the solid and normal to the pore walls. These contributions result in the equations:[3] ϵBangham=−1MPLΔfrandϵLaplace=1MPLpF,

where MPL is the hydrostatic pressure-associated pore-load modulus (13), Δf is the released surface stress, i.e., the difference in surface stress before and after imbibition, and r=rp−d is the corrected pore radius, which appears in Eq. 1.

For isotropically distributed cylindrical pores, the total strain of a liquid-filled representative volume ϵ∗ is calculated then as a combination of the two strains:[4] ϵ∗=1MPL−Δfr+pF.

On the porous-medium scale, the fluid pressure increases linearly from the negative, contractile Laplace pressure pL at the position of the advancing imbibition front h(t) up to the bulk reservoir pressure p0, while the change in surface stress per unit length for the imbibed part of the pore network is constant along the sample height at any imbibition time t. Thus, the Laplace pressure-induced deformation is maximal at the imbibition front and vanishes at the bulk reservoir. By contrast, the Bangham effect is constant along the imbibed part of the matrix, so that the water-filled part of the originally dry, cylindrical monolith deforms to a frustum (see Fig. 1E for the corresponding trapezoidal deformation in side view). The resulting overall deformation and thus relative sample length change ΔL/L0 in the imbibition direction is thus proportional to the fraction of imbibed sample (h(t)/hfinal). And because of the linearly increasing magnitude of the fluid pressure from the bulk reservoir to the imbibition front, it is only half as large as it would be for a porous solid homogeneously experiencing the tensile pressure pL. Additionally considering that for nanoporous materials the atmospheric pressure is much smaller than the Laplace pressure, p0≪pL, we arrive at:[5] ϵ(t)=−1MPLγ+Δfrh(t)hfinal.

At imbibition time tc the average position of the imbibition front is equal to the height of the sample L or, in other words, using Eq. 1, h(tc)=L. For a single nanocapillary, the tensile Laplace pressure vanishes suddenly at tc, but for our porous matrix such a vanishing is not instantaneous but gradual, since the vanishing dynamics is governed by the broadened imbibition front. The Laplace pressure contribution starts to decrease when the upper limit of the broadening front reaches h(t1)=L and continues to do so until its disappearance when the whole matrix is filled at t2. This time difference t2−t1 over which the Laplace pressure vanishes is indicated in all of our figures as regime II (light blue background). The extent of the Laplace jump is[6] Δϵ=ϵ(t2)−ϵ(t1)=pFMPL=−pL2MPL=−γ(rp−d)MPL.

The individual contributions to the strain and consequently the total strain are time dependent and evolve proportional to the MLW dynamics of Eq. 1. Fig. 3A shows a model of both strain contributions and, how the combination of both allows reconstructing the observed deformation dynamics represented by the green curve. This figure panel illustrates the competitive deformation process and does not represent an actual experiment. In this way, we understand the imbibition-induced strain curve of water in nanoporous Vycor glass and why, upon filling, we observe a dynamic expansion proportional to t; the green curve acceptably describes the first stage of the process, prior to a strain jump that lasts until time tc of complete filling.

Fig. 3. Contributions to imbibition-induced deformation and influence of preimbibition humidity. (A) Contributions to deformation of the porous medium along the imbibition direction that originate from contractile Laplace pressure (purple) and expansive surface stress contributions (red) along with the resulting overall deformation (green line), as calculated by a linear superposition according to our elastocapillary model. The maximum values of both contributions are represented by the dotted lines. (B) Strain isotherm measurement for a Vycor glass rod rp= 3.4 nm, T= 25.8 °C and p0= 3.3 kPa, with the logarithmic vanishing of the Laplace pressure highlighted with a black line (adapted from ref. 50). Overlapping with such dataset, we include our own experimental data (green points) obtained by subtracting the maximum strain (Bangham contribution) in our imbibition experiments to the maximum strain value in the isotherm (0.2%) represented by the green dashed line. The red dashed lines represent the measured strain in our dilatometry experiment and the position along the Y-axis of the influence of the adsorption-induced strain due to the humidity. (C) Water imbibition-induced deformation in a nanoporous Vycor sample rp= 3.4 nm with preimbibition humidity RH2≈ 50%) (dark green points) compared to one (light green dots) performed with lower preimbibition humidity (RH1≈ 5%). The strain jump and maximum strain are given by the red dashed and dotted lines for RH1 and RH2, respectively. (D) Schematic illustrations of the interplay of negative hydrostatic pressure (purple) and a sample-expanding surface stress release at the water–silica wall interface per unit of sample length (red) for specific rise heights (h): (h1), (h2), (h3), and (h4). Analogously, for a selected time ts, in (E), this interplay is shown at the single-pore scale for two distinct preimbibition humidity values. The white arrows represent the parabolic velocity profile for pressure-driven flow in a cylindrical channel.

The strain jump in our dilatometry measurement can be quantified according to Eq. 6 provided the pore-load modulus MPL of the material is known. MPL is extractable through a sorption strain isotherm measurement (24). The water adsorption–desorption-induced strain on nanoporous Vycor rods was experimentally studied in the 1950s by McIntosh and Amberg (50). We work with nanoporous Vycor hosts with comparable pore size and porosity, and at identical temperature. The pore-load modulus of Vycor is accessible via the Kelvin-Laplace equation upon reaching bulk liquid–vapor coexistence. The McIntosh and Amberg data show the logarithmic vanishing of the Laplace pressure upon approaching bulk liquid coexistence and thus vanishing menisci curvature (Fig. 3B). It results in MPL= 39.45 gPa. Introducing this value to Eq. 6 while assuming the bulk surface tension of water (72 mN m−1) we get a value for the maximum Laplace pressure-induced strain of 0.036%, which is in reasonable agreement with the experimental value of 0.032% (Fig. 3C). Here the final strain value after the jump (dotted red line) corresponds to the Bangham effect contribution to the strain.

An identical experiment is performed for a relative humidity of ∼5% and also shown in Fig. 3C. Again, the three different deformation regimes are observed. Note, however, that the total relative expansion after complete filling is 0.16%, almost twice as large as in the experiment at RH= 50%, while the strain jump at complete filling is the same in both experiments.

Whereas the effect of the Laplace pressure on the strain upon imbibition is predictable, the surface stress contribution is highly dependent on the preimbibition atmosphere. Higher humidity values in our experimental environment, prior to the start of the imbibition experiment, will induce larger preimbibition adsorption of water molecules on the pore walls, leading to expansion of the porous host (as illustrated in Fig. 3D). Such dependency is evidenced by comparing the dilatometry experiments at different preimbibition humidities. Again in Fig. 3B we show how, for a certain relative humidity precondition RH1, the measurable imbibition-induced deformation (c) is comparable to the strain difference between equilibrium states at RH1 and RHf=p/p0=1 (fully filled sample).

From the final strain values, it is possible to calculate Δf using Eq. 4. We find Δf(RH=50%)=−116 mJ m−2 and Δf(RH=5%)=−206 mJ m−2. Our simple semiempirical approach to calculate Δf from ϵ is in reasonable agreement with an effective medium-based description of surface stress-induced mechanical deformation as proposed by Weissmüller et al. (51). For a network of interconnected elongated fibers with circular cross-section they propose[7] ϵ=−αΔf3K1−ν1−2ν,

with α being the volume-specific surface area and K being the bulk modulus. The Δf values obtained in this way are −116 mJ m−2 and −190 mJ m−2 for RH=50% and RH=5%, respectively.

Considering the common case where the surface stress release equals the difference between the surface energy of the dry solid and the wetted solid, Δf=γsl−γsv (sl: solid–liquid, sv: solid–vapor) (52), Eq. 4 can be written as[8] ϵ∗=−1MPLγ+γsl−γsvr=1MPLSr,

where S=γsv−(γsl+γ), is the spreading coefficient at the liquid meniscus contact line of the wetting theory (53). Hence, for the complete wetting case, Eq. 8 always predicts a positive strain, and we arrive at the important conclusion that for any capillarity-driven, spontaneous imbibition process, and for any liquid/porous solid combination, the Bangham effect prevails over the Laplace effect, i.e., in regime I, an expansive behavior will always be observed. Note that the complete wetting case assumes a zero liquid/solid contact angle θ. The inclusion of a finite contact angle (θ>0) implies the replacement of γ by γcosθ both for the calculation of the Laplace pressure and the spreading coefficient.

Unfortunately, the surface energy values for dry amorphous silica and hydrous amorphous silica interfaces reported in the literature are quite inconsistent, presumably because of the complex water–silica interaction so that they depend sensitively on any contamination, the thickness of the adsorbed water layers, and the surface chemistry (silanol density and hydroxylation of the surface) (54, 55). However, for the surface energies (γ) reported in ref. 54 (γsl= 130 mJ m−2 and γsv= 260 mJ m−2, respectively) we arrive at a predicted surface stress release of Δf=−130 mJ m−2 for our RH= 50% measurements, in reasonably good agreement with our finding of −116 mJ m−2. This also corroborates the conclusion in ref. 24 that for water–silica interfaces the surface stress corresponds to the surface energy changes.

Molecular Dynamics Simulations of Water Imbibition-Induced Deformation on the Single-Nanopore Scale.

As demonstrated above, dilatometry experiments allow the study of the macroscopic effects of imbibition-induced strain on porous host materials and a quantification of the total strain. However, the individual strain contributions originating from the superposition of the Bangham and Laplace pressure effects on the single-pore scale can not be investigated with this technique.

To better understand how these microscopic strain contributions accumulate to yield the total macroscopically measured deformation and how they are related to individual menisci movements in the pore space, Molecular Dynamics (MD) computer simulations of the imbibition of water into silica nanopores on the single-pore scale are conducted. Since in MD simulations atom–atom interactions for each individual atom in the system are modeled they are well suited to investigate the microscopic imbibition-induced strain contributions on a fundamental level. The porous silica system used is based on the original work of refs. 56 and 57 and has a diameter of around 3 nm and a wall thickness of around 1 nm. The total length the water can imbibe is about 12 nm.

Fig. 4 shows the results of an ensemble of 20 MD simulations with independent starting conditions; see Movie S3 of one of the simulations in SI Appendix. Fig. 4A displays snapshots of slices through a nanopore during the imbibition process. The snapshots are arranged from left to right in chronological order at times 0 ps, 183 ps, 585 ps, and 780 ps as indicated by the dotted vertical lines in Fig. 4B–D. In the third snapshot, a zoom-in is provided that focuses on the water meniscus formed at the imbibing water column. Note that with a 3 nm pore diameter, the silica pore size is only around one order of magnitude larger than the size of a water molecule as visible from the Inset.

Fig. 4. Elastocapillary dynamics upon water imbibition on the single-nanopore scale as inferred from molecular dynamics computer simulations. (A) Snapshots from an exemplary simulation at selected imbibition times as indicated in the figure. Silicon atoms are depicted in yellow, oxygen in red, and hydrogen in white. The imbibing water is represented through the blue-colored surface. (B) Imbibition length h of the water column in the silica pores from one pore side in comparison with a t fit (black). (C and D) Observed lateral ϵ⊥ and longitudinal ϵ‖ strain in the pore matrix.

The imbibition length h of the water column into the silica pore from one pore side is shown as a function of imbibition time t in Fig. 4B. The error bars are omitted since the plot markers are larger than the SD of the simulation ensemble. In the simulations, the imbibition length is measured indirectly through the extent of the simulation box in the longitudinal pore direction. The initial imbibition dynamics of around 3 ps follows a linear trend, which is consistent with an initial water inertia-governed dynamics during which the influence of the water viscosity can be neglected (see SI Appendix for further information). After the initial 3 ps the imbibition dynamics exhibits a t dependence as shown by the black fit, which confirms the validity of the LW-model even for pore sizes that only fit a few adjacent water molecule layers. These results are in agreement with previous MD simulations of nanopore imbibition for simple Lennard-Jones fluids and polymer melts (58). The pore is completely filled with water after 780 ps. The measurement of the imbibition length through the extent of the simulation box in the longitudinal pore direction assumes a constant pore space volume. For the process of imbibition, this is a good approximation since the pore space volume only changes gradually. However, at the end of the imbibition process, a relatively strong lateral pore expansion occurs (Fig. 4C). This lateral expansion leads to a significantly increased pore space volume that triggers an additional water inflow from the water reservoir into the pore space, even though the pore is already entirely filled by the imbibing water column. The inflow of the water in turn reduces the reservoir volume and leads to a shrinkage of the simulation box in the longitudinal pore direction. This effect explains the additional increase of measured imbibition length even after the water columns fill the entire pore space.

The green markers in Fig. 4 C and D display the dynamics of the lateral ϵ⊥(t) (perpendicular to the long pore axis) and longitudinal ϵ‖(t) (parallel to the long pore axis) strains, respectively. Error bars indicate the SDs over the averaged 20 individual simulations. The use of different analysis methods to determine the lateral and the longitudinal strains explains the different magnitudes of SD in the measured strains (Materials and Methods). The strong changes in strain during the first time steps of the simulation are in response to the initial inertia-dominated imbibition and initial system equilibration dynamics.

The temporal strain evolution, laterally ϵ⊥(t) and longitudinally ϵ‖(t), can be divided in three regimes. Regime I from 0 ps to about 780 ps displays an expansion of the host material with t dynamics as indicated by the black curve fit. In agreement with our macroscopic experiments on the porous-medium scale, the expansive Bangham effect dominates in comparison to the contractile Laplace pressure effect on the single-pore scale. However, the influence of the Laplace pressure effect is different for the two strain directions. In the lateral direction, the negative Laplace pressure leads to a pronounced contractile effect. In contrast, the longitudinal strain is only influenced by the Laplace pressure through the Poisson effect. More precisely, the negative Laplace pressure leads to an expansion in pore wall thickness, which, in turn through the Poisson effect, leads to a contractile effect of small magnitude in the longitudinal pore direction. At the end of regime I the menisci at the water columns vanished upon complete filling of the pore space, which leads to a disappearance of the contractile Laplace pressure and the resulting strain effects in the lateral and longitudinal directions. Consequently the strain increases during the relaxation of the porous material from 780 ps to 1,200 ps in regime II. Regime III occurs after the expansion around 1,200 ps when the porous material exhibits the final observed strain created by the Bangham effect acting on the entire inner surface of the nanopore.

Thus, overall the strain evolution at the single-pore scale is analogous to that observed at the macroscopic porous-medium scale. In particular, by comparison with the single-pore simulation, the Laplace pressure strain jump observed when the imbibition front reaches the top of the Vycor matrix can be attributed to a collective disappearance of menisci in Vycor.

Analogous to the macroscopic case, for the simulated material, MPL can be calculated from Eq. 6 with the Laplace jump Δϵ= 1%. This results in a pore-load modulus value of MPL= 4.8 GPa, an order of magnitude smaller than the 39.45 GPa calculated for Vycor glass. This is due to the significantly smaller silica pore wall thicknesses in the simulation compared to the experiment and also explains why the magnitude of the imbibition-induced strains is an order of magnitude larger than the experimental one. From the total strain we calculated Δfsimulation=Δγsimulation=−177 mJ m−2 provided that pL=(γcosθ)/rp and the contact angle is about θ= 45°, as can be inferred from the meniscus snapshots in Fig. 4. We find that this value (starting from vacuum) is the same order of magnitude as the calculated values of Δγ at low humidity in the macroscopic case (between −190 mJ m−2 and −206 mJ m−2). Thus, the single-pore simulations corroborate not only qualitatively, but also (semi)quantitatively the experimental findings on the porous-medium scale.

Summary

We have studied experimentally and by computer simulation, the deformation behavior of a mesoporous silica monolith (Vycor) and a single cylindrical nanopore, respectively, as a function of water imbibition. In contrast to the expected Laplace pressure-induced contraction, we observed two distinct imbibition-induced expansion kinetics, which can be directly attributed to the advance, arrest, and disappearance of the imbibition front or single nanoscale menisci, respectively. In particular, the surface stress release-induced expansion known as the Bangham effect was found to dominate the Laplace pressure-induced contraction in both cases. As we derive from general wetting theory, this is the generic case for any spontaneous imbibition process and any liquid/solid combination, provided that the surface stress release is equal to the surface energy difference between the wetted and dry solid.

Notably, there is also an important difference in the Laplace contribution to the deformation and thus in the total deformation between filling by classical, quasi-static water vapor adsorption and dynamic liquid water imbibition. In the case of gas sorption, the entire filled part is under a homogeneous negative Laplace pressure (18, 25). In the case of spontaneous imbibition, the tensile pressure in the porous medium decreases linearly from the advancing menisci to the bulk reservoir due to viscous liquid flow. This results in exactly half the Laplace pressure-induced contraction of the porous medium compared to the situation with gas adsorption (capillary condensation/evaporation). Since the contractile Laplace effect counteracts the dominant expansive Bangham effect, the total deformation is always smaller for gas sorption than for liquid imbibition. Thus, for fragile porous solids, filling by gas sorption should generally be advantageous, since smaller mechanical strains are expected.

The observed competition between surface stress and Laplace pressure should be relevant for any imbibition phenomenon in nanoporous media. For significantly stiff systems, where deformation leads to negligible changes in pore size and thus hydraulic permeability, simple LW deformation kinetics are expected. For soft solids, on the other hand, the deformation could lead to linear hydraulic resistance gradients and thus to deviations of both the liquid advancement kinetics and the deformation dynamics from the classical LW behavior. Also, for extremely narrow, subnanometer pores, repulsive hydration layer forces could have a nontrivial influence on the elastic imbibition response.

Our results show that relative length changes can be used to study the liquid filling kinetics and thus the fabrication of the emerging class of liquid-infused functional materials in detail, but also under in operando conditions a study of the deformation behavior allows for monitoring the appropriate filling state of a liquid-infused materials system.

Materials and Methods

Humidity-Dependent Dilatometry Experiments.

As porous host, we work with cylindrical monoliths of porous Vycor glass (Corning Glass, 7930). All of the samples are identical, only differing in length, and were previously used by Simon Grüner. An extensive description of such samples can be found in ref. 37. Important properties such as the mean pore radius rp= 4.9 nm and porosity Φ0=0.3 are obtained from a water sorption isotherm measurement. Gruener et al. (37) assumed that there are two preadsorbed monolayers of water molecules of thickness d= 0.5 nm attached to the pore walls prior to the sorption isotherm measurements. Since we measure in comparable atmospheric conditions, we will assume the same value for our measurements. This magnitude of water layer thickness is consistent with the reduced porosity Φ0=0.27 that we estimate by measuring the sample mass difference between the beginning and end of the experiment. The tortuosity value of τ=3.6 is also taken from the above mentioned study.

The samples are precleaned by immersing them in a mixture of concentrated sulfuric acid (H2SO4) and a 30% hydrogen peroxide solution (H2O2) with a 3:1 ratio (Piranha etch solution). The strongly oxidizing strength of this solution removes most organic matter from the pore walls while hydroxylating most surfaces (i.e., adding OH groups) in the process. The porous Vycor monoliths are immersed in the solution for 3 wk and subsequently for 3 d in ultrapure water. After drying the process is repeated one more time. Before the experiment, the samples are dried by introducing them in a vacuum chamber for 48 h. Using pure water for every experiment avoids the necessity to repeat the whole cleaning process between different measurements, hence in this case only the drying is to be performed multiple times.

For the measurements, we use a vertically oriented dilatometer (Linseis L75 VS500) with a precision up to 0.1 μm and a temporal resolution of 1 s.

Humidity control during our experiments is achieved by using a custom-built humidity chamber that can be integrated into the dilatometry setup. The desired relative humidity can be maintained by controlling the incoming air flow, which is a mixture of pure dry nitrogen and water vapor. The samples are exposed to the desired preimbibition humidity, where they remain until equilibrium is reached (monitored by dilatometry). After equilibration, the initially empty liquid reservoir is filled by injecting liquid through an external syringe until the water level reaches the sample, initiating the capillary filling process.

To seal the side facets of our samples, we used pressure-sensitive adhesive (PSA) tape provided by Tesa®. The adhesive tape consists of a polypropylene–polyethylene backing film for transparency and flexibility, paired with a polymeric adhesive layer. Both the adhesive and the liner are hydrophobic, which, combined with the viscoelastic properties of the tape (very low elastic modulus), makes it an excellent solution for sealing our specimens without compromising the mechanical stiffness of the monoliths. The adhesive mechanism of the PSA is based on the spreading of the polymer over the applied surface under pressure and the bonding of the polymers to the surface by van der Waals forces. Since the strain-driven pressures involved in the deformation process are on the order of MPa, the effect of the sealing on the monolith deformation is neglected. The Tesa® tape used is optimized for easy removal with minimal residue. Using a magnifying glass, we have observed that small traces of polymer remain on the surface after the PSA is removed in between measurements. We therefore treat the surface with isopropanol to remove the remaining residues and to facilitate evaporation after washing the samples with ultrapure water.

Molecular Dynamics Simulations.

Average quantities from the imbibition process are obtained from 20 individual Molecular Dynamics simulations with alternating starting velocities of the water molecules conducted with the Large-scale Atomic/Molecular Massively Parallel Simulator (59). As host material, a periodically repeated amorphous silica pore initially introduced in refs. 56 and 57 is used. The silica–silica interactions are modeled by a modified Demiralp potential (60), water–silica interactions are modeled as described in refs. 61 and 62, and water–water interactions are modeled with the TIP4P/2005 water model (63). The simulated pore structure consists of amorphous silica SiO2. The surface of the porous structure is saturated with hydroxyl groups of 7.23 nm−2 density. The total longitudinal pore length is around 13 nm along the z axis of the simulation box. However, in the simulation, the water column starts imbibing from just inside the pore leaving around 12 nm of possible imbibition length. In the lateral direction, that is the x and y direction of the simulation box, the pores are hexagonally shaped and have a diameter of around 3 nm while the thickness of the pore walls is around 1 nm. The imbibition process is simulated with periodic boundary conditions in all directions in an isothermal–isobaric ensemble (NpT) at 1 bar and 300 K employing a Nosé-Hover barostat (64–67). To trigger the imbibition process a water reservoir outside of the nanopores at the longitudinal pore ends is created. Due to the periodic boundary conditions the water can enter the pore from both sides. The x and y axes are coupled such that the rate of expansion or contraction of the simulation box is equal for both axes. The total simulated time is 2,000 ps employing a 1 fs timestep (see SI Appendix for further information).

The measured strains are defined as Cauchy strains ϵ(t)=[L(t)−L(t=t0)]/L(t=t0), where L(t=t0) constitutes the pore extensions at the start of the simulation. The radial strain ϵ⊥(t) is measured indirectly through the change of the simulation box extension Lbox in the lateral direction. In contrast, the longitudinal strain ϵ‖(t) is measured by the distances between the groups of atoms forming the last 0.5 nm of bulk material of the pores on each end (see SI Appendix for further information). The thermal motions of the atoms lead to strongly varying distances during the simulation. These fluctuations account for the high SD in the measured longitudinal strain represented in the error bars in Fig. 4.

Supplementary Material

Appendix 01 (PDF)

Movie S1. Video of our optical imaging experiment on water imbibition in nanoporous Vycor glass, in absence of evaporation. The sample length is ~ 12 mm, and the side and top of the sample are sealed using transparent adhesive film (TESA film®). The real time duration of the experiment is 2.05 hours.

Movie S2. Video of our optical imaging experiment on water imbibition in nanoporous Vycor glass, with evaporation. The sample length is ~ 19 mm, and the side and top of the sample are not sealed allowing evaporation. The real time duration of the experiment is 27.87 hours.

Movie S3. Video of snapshots originating from the first simulated 780 ps (regime I) of one exemplary imbibition simulation trajectory. The simulation ensemble used to investigate the imbibition dynamics and the imbibition induced strain consists of 20 trajectories with independent starting conditions.

We are indebted to Klaus Knorr (Mainz, Germany) and Jörg Weißmüller (Hamburg, Germany) for a critical reading of the manuscript and for their helpful comments. J.S., P.H., and R.H.M. acknowledge funding by the Deutsche Forschungsmeinschaft (DFG, German Research Foundation) within the Research Training Group 2462 “Processes in natural and technological Particle-Fluid-Systems” (Project No. 390794421) and within the Collaborative Research Centre 986 “Tailor-Made Multi-Scale Materials Systems” (Project No.192346071). L.D. was supported by the Data Science in Hamburg HELMHOLTZ Graduate School DASHH, Helmholtz Association Grant-No. HIDSS-0002. L.G. profited from the DFG research grant “Dynamic Electrowetting at Nanoporous Surfaces: Switchable Spreading, Imbibition, and Elastocapillarity,” Project number 422879465 (Schwerpunktprogramm 2171), and from the Cluster of Excellence “Understanding Written Artefacts” (Exzellenzcluster 2176). This research was supported in part through the Maxwell computational resources operated at Deutsches Elektronen-Synchrotron, Hamburg, Germany. M.F and P.H. thank the Free and Hanseatic City of Hamburg for funding within the research group “Control of the special properties of water in nanopores.” We also acknowledge the scientific exchange and support of the Centre for Molecular Water Science, Hamburg (Germany).

Author contributions

J.S., L.D., Z.L., and P.H. designed research; J.S., L.D., L.G., Z.L., and P.H. performed research; J.S., L.D., L.G., Z.L., R.H.M., H.A.S., and P.H. analyzed data; and J.S., L.D., L.G., Z.L., M.F., R.H.M., H.A.S., and P.H. wrote the paper.

Competing interests

The authors declare no competing interest.

Data, Materials, and Software Availability

Experimental dilatometry (68) and Molecular Dynamics (69) simulation data have been deposited in TORE.

Supporting Information

This article is a PNAS Direct Submission. A.D.S. is a guest editor invited by the Editorial Board.
==== Refs
1 X. Yao , Adaptive fluid-infused porous films with tunable transparency and wettability. Nat. Mater. 12 , 529–534 (2013).23563739
2 W. Wang , Multifunctional ferrofluid-infused surfaces with reconfigurable multiscale topography. Nature 559 , 77–82 (2018).29942075
3 M. Brinker, P. Huber, Wafer-scale electroactive nanoporous silicon: Large and fully reversible electrochemo-mechanical actuation in aqueous electrolytes. Adv. Mater. 34 , 2105923 (2022).
4 K. Sentker , Self-assembly of liquid crystals in nanoporous solids for adaptive photonic metamaterials. Nanoscale 11 , 23304–23317 (2019).31788679
5 L. G. Cencha, G. Dittrich, P. Huber, C. L. A. Berli, R. Urteaga, Precursor film spreading during liquid imbibition in nanoporous photonic crystals. Phys. Rev. Lett. 125 , 234502 (2020).33337190
6 T. Emmerich , Enhanced nanofluidic transport in activated carbon nanoconduits. Nat. Mater. 21 , 696–702 (2022).35422506
7 P. Simon, Y. Gogotsi, Materials for electrochemical capacitors. Nat. Mater. 7 , 845–854 (2008).18956000
8 J. Zhang, B. Chen, X. Chen, X. Hou, Liquid-Based Adaptive Structural Materials. Adv. Mater. 2005664 , 33 (2021).
9 T. S. Wong , Bioinspired self-repairing slippery surfaces with pressure-stable omniphobicity. Nature 477 , 443–447 (2011).21938066
10 F. T. Meehan, The expansion of charcoal on sorption of carbon dioxide. Proc. R. Soc. Lond. Ser. A Containing Pap. Math. Phys. Character 115 , 199–207 (1927).
11 D. H. Bangham, N. Fakhoury, The expansion of charcoal acompanying sorption of gases and vapours. Nature 122 , 681–682 (1928).
12 G. W. Scherer, Dilatation of porous glass. J. Am. Ceram. Soc. 69 , 473–480 (1986).
13 J. Prass, D. Mueter, P. Fratzl, O. Paris, Capillarity-driven deformation of ordered nanoporous silica. Appl. Phys. Lett. 95 , 83121 (2009).
14 G. Y. Gor, A. V. Neimark, Adsorption-induced deformation of mesoporous solids. Langmuir 26 , 13021–13027 (2010).20695535
15 K. Schappert, R. Pelster, Unexpected sorption-induced deformation of nanoporous glass: Evidence for spatial rearrangement of adsorbed argon. Langmuir 30 , 14004–14013 (2014).25358117
16 C. Balzer , Relationship between pore structure and sorption-induced deformation in hierarchical silica-based monoliths. Z. Phys. Chem. 229 , 1189–1209 (2014).
17 A. Grosman, J. Puibasset, E. Rolley, Adsorption-induced strain of a nanoscale silicon honeycomb. Europhys. Lett. 109 , 56002 (2015).
18 G. Y. Gor, P. Huber, N. Bernstein, Adsorption-induced deformation of nanoporous materials—A review. Appl. Phys. Rev. 4 , 011303 (2017).
19 M. Chen, B. Coasne, R. Guyer, D. Derome, J. Carmeliet, Molecular simulation of sorption-induced deformation in atomistic nanoporous materials. Langmuir 35 , 7751–7758 (2019).31117732
20 Q. Yang , Capillary condensation under atomic-scale confinement. Nature 588 , 250–253 (2020).33299189
21 S. G. Harrellson , Hydration solids. Nature 619 , 500–505 (2023).37286609
22 G. Y. Gor, G. W. Scherer, H. A. Stone, Bacterial spores respond to humidity similarly to hydrogels. Proc. Natl. Acad. Sci. U.S.A. 121 , e2320763121 (2024).38416680
23 G. Y. Gor , Elastic response of mesoporous silicon to capillary pressures in the pores. Appl. Phys. Lett. 106 , 1–13 (2015).
24 G. Y. Gor, N. Bernstein, Revisiting Bangham’s law of adsorption-induced deformation: Changes of surface energy and surface stress. Phys. Chem. Chem. Phys. 18 , 9788–9798 (2016).27001041
25 G. Y. Gor , Elastic response of mesoporous silicon to capillary pressures in the pores. Appl. Phys. Lett. 106 , 261901 (2015).
26 C. Duprat, J. Aristoff, H. A. Stone, Dynamics of elastocapillary rise. J. Fluid Mech. 679 , 641–654 (2011).
27 J. H. Weijs, B. Andreotti, J. H. Snoeijer, Elasto-capillarity at the nanoscale: On the coupling between elasticity and surface energy in soft solids. Soft Matter 9 , 8494–8503 (2013).
28 M. Tortora , Giant negative compressibility by liquid intrusion into superhydrophobic flexible nanoporous frameworks. Nano Lett. 21 , 2848–2853 (2021).33759533
29 L. Michel , Bowtie-shaped deformation isotherm of superhydrophobic cylindrical mesopores. Langmuir 38 , 211–220 (2022).34964631
30 T. B. Hoberg, E. Verneuil, A. E. Hosoi, Elastocapillary flows in flexible tubes. Phys. Fluids 26 , 122103 (2014).
31 J. Ha , Poro-elasto-capillary wicking of cellulose sponges. Sci. Adv. 4 , 1–7 (2018).
32 J. I. Siddique, D. M. Anderson, A. Bondarev, Capillary rise of a liquid into a deformable porous material. Phys. Fluids 21 , 013106 (2009).
33 P. Levitz, G. Ehret, S. K. Sinha, J. M. Drake, Porous Vycor glass: The microstructure as probed by electron microscopy, direct energy transfer, small-angle scattering, and molecular adsorption. J. Chem. Phys. 95 , 6151 (1991).
34 P. Huber, K. Knorr, Adsorption-desorption isotherms and X-ray diffraction of AR condensed into a porous glass matrix. Phys. Rev. B 60 , 12657 (1999).
35 P. Levitz, Off-lattice reconstruction of porous media: Critical evaluation, geometrical confinement and molecular transport. Adv. Colloid Interface Sci. 76–77 , 71–106 (1998).
36 C. J. Gommes, Stochastic models of disordered mesoporous materials for small-angle scattering analysis and more. Microporous Mesoporous Mater. 257 , 62–78 (2018).
37 S. Gruener, T. Hofmann, D. Wallacher, A. V. Kityk, P. Huber, Capillary rise of water in hydrophilic nanopores. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 79 , 067301 (2009).
38 S. Gruener , Anomalous front broadening during spontaneous imbibition in a matrix with elongated pores. Proc. Natl. Acad. Sci. U.S.A. 109 , 10245–10250 (2012).22689951
39 S. Gruener, H. E. Hermes, B. Schillinger, S. U. Egelhaaf, P. Huber, Capillary rise dynamics of liquid hydrocarbons in mesoporous silica as explored by gravimetry, optical and neutron imaging: Nano-rheology and determination of pore size distributions from the shape of imbibition fronts. Colloids Surf. A Physicochem. Eng. Asp. 496 , 13–27 (2015).
40 J. H. Page, J. Liu, B. Abeles, H. W. Deckman, D. A. Weitz, Pore-space correlations in capillary condensation in Vycor. Phys. Rev. Lett. 71 , 1216–1219 (1993).10055479
41 C. H. Amberg, R. McIntosh, A study of adsorption hysteresis by means of length changes of a rod of porous glass. Can. J. Chem. 30 , 1012–1032 (1952).
42 L. G. L. Cencha , Nondestructive high-throughput screening of nanopore geometry in porous membranes by imbibition. Appl. Phys. Lett. 115 , 113701 (2019).
43 S. Gruener, P. Huber, Imbibition in mesoporous silica: Rheological concepts and experiments on water and a liquid crystal. J. Phys. Condens. Matter 23 , 184109 (2011).21508488
44 A. Schlaich, J. Kappler, R. R. Netz, Hydration friction in nanoconfinement: From bulk via interfacial to dry friction. Nano Lett. 17 , 5969–5975 (2017).28910108
45 B. Maillet, G. Dittrich, P. Huber, P. Coussot, Diffusionlike drying of a nanoporous solid as revealed by magnetic resonance imaging. Phys. Rev. Appl. 10 , 1 (2022).
46 E. Seker, M. R. Begley, M. L. Reed, M. Utz, Kinetics of capillary wetting in nanoporous films in the presence of surface evaporation. Appl. Phys. Lett. 92 , 013128 (2008).
47 T. D. Wheeler, A. D. Stroock, The transpiration of water at negative pressures in a synthetic tree. Nature 455 , 208 (2008).18784721
48 C. K. Camplisson, K. M. Schilling, W. L. Pedrotti, H. A. Stone, A. W. Martinez, Two-ply channels for faster wicking in paper-based microfluidic devices. Lab Chip 15 , 4461–4466 (2015).26477676
49 Z. Huang , Fast water evaporation from nanopores. Adv. Mater. Interfaces 8 , 2100660 (2021).
50 C. H. Amberg, R. Mcinthosh, A study of adsorption hysteresis by means of length changes of a rod of porous glass. Can. J. Chem. 30 , 1012–1032 (1952).
51 J. Weissmüller, H. L. Duan, D. Farkas, Deformation of solids with nanoscale pores by the action of capillary forces. Acta Mater. 58 , 1–13 (2010).
52 G. Y. Gor, A. V. Neimark, Adsorption-induced deformation of mesoporous solids. Langmuir ACS J. Surf. Colloids 26 , 13021–13027 (2010).
53 P. G. de Gennes, F. Brochard-Wyart, D. Quere, Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves (Springer, New York, 2004).
54 S. Brunauer, D. L. Kantro, C. H. Weise, The surface energies of amorphous silica and hydrous amorphous silica. Can. J. Chem. 34 , 1483–1496 (1956).
55 D. Díaz , How water wets and self-hydrophilizes nanopatterns of physisorbed hydrocarbons. J. Colloid Interface Sci. 606 , 57–66 (2022).34388573
56 P. Ugliengo , Realistic models of hydroxylated amorphous silica surfaces and MCM-41 mesoporous material simulated by large-scale periodic B3LYP calculations. Adv. Mater. 20 , 4579–4583 (2008).
57 M. Delle Piane, M. Corno, A. Pedone, R. Dovesi, P. Ugliengo, Large-scale B3LYP simulations of ibuprofen adsorbed in MCM-41 mesoporous silica as drug delivery system. J. Phys. Chem. C 118 , 26737–26749 (2014).
58 D. I. Dimitrov, A. Milchev, K. Binder, Capillary rise in nanopores: Molecular dynamics evidence for the Lucas-Washburn equation. Phys. Rev. Lett. 99 , 54501 (2007).
59 A. P. Thompson , LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comput. Phys. Commun. 271 , 108171 (2022).
60 R. H. Meißner, J. Schneider, P. Schiffels, L. Colombi Ciacchi, Computational prediction of circular dichroism spectra and quantification of helicity loss upon peptide adsorption on silica. Langmuir 30 , 3487–3494 (2014).24627945
61 D. J. Cole, M. C. Payne, G. Csányi, S. M. Spearing, L. C. Ciacchi, Development of a classical force field for the oxidized Si surface: Application to hydrophilic wafer bonding. J. Chem. Phys. 127 , 204704 (2007).18052443
62 A. Butenuth , Ab initio derived force-field parameters for molecular dynamics simulations of deprotonated amorphous-SiO2/water interfaces. Phys. Stat. Solidi 249 , 292–305 (2012).
63 J. L. Abascal, C. Vega, A general purpose model for the condensed phases of water: TIP4P/2005. J. Chem. Phys. 123 , 234505 (2005).16392929
64 G. J. Martyna, D. J. Tobias, M. L. Klein, Constant pressure molecular dynamics algorithms. J. Chem. Phys. 101 , 4177–4189 (1994).
65 M. Parrinello, A. Rahman, Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 52 , 7182–7190 (1981).
66 L. N. Coregliano , A Liouville-operator derived measure-preserving integrator for molecular dynamics simulations in the isothermal-isobaric ensemble. J. Phys. A Math. Gen. 39 , 5629 (2006).
67 W. Shinoda, M. Shiga, M. Mikami, Rapid estimation of elastic constants by molecular dynamics simulation under constant stress. Phys. Rev. B 69 , 134103 (2004).
68 J. Sanchez Calzado, Deformation dynamics of nanopores upon water imbibition: Optical, gravimetric and dilatometry experiments. TORE. 10.15480/882.13233. Deposited 21 August 2024.
69 L. Dammann, Deformation dynamics of nanopores upon water imbibition: Molecular dynamics simulations. TORE. 10.15480/882.13223. Deposited 21 August 2024.
