Low Pressure Hysteresis in Materials with Narrow Slit Pores

Humidity-dependent closing and reopening slit pores can produce hysteresis loops in sorption diagrams even at low relative vapor pressure. Pore closing is supported by adhesion of the slit wall surfaces. In a macroscopic model for sorption hysteresis in narrow slits, the adhesion energy jumps by a finite value when touching slit walls are separated from each other. We consider a more realistic adhesion model by introducing a smoothly-varying adhesion force, which depends on the distance between the slit walls. The range of the attraction between the slit walls is found to have a pronounced influence on the shape of hysteresis loops at low vapor pressure. A large interaction range avoids an extraordinarily small relative vapor pressure necessary for pore closing, which is a precondition for low pressure hysteresis. Our extended model allows us to describe a discontinuity, which can appear in the desorption branch of swelling/shrinkage diagrams for hardened cement paste.


Introduction
Cracking of porous solids like concrete or ceramic green bodies due to capillary or solvation forces is a serious problem in many technological processes.Variations of moisture content in porous materials are associated with compression and swelling of the solid skeleton [1].If there is a thermodynamic equilibrium between adsorbed water and vapor in the environment, the capillary force depends on the relative air humidity.An additional effect accompanied by the shrinkage and swelling of porous materials is sorption hysteresis at low vapor pressure.
In most porous materials, a precondition of sorption hysteresis is supposed to be capillary condensation in mesopores.Mesopores have diameters or thicknesses larger than a few nanometers.For example, in rigid materials such as Vycor glass, water sorption hysteresis appears at a relative vapor pressure larger than about 0.4.Generally, in the case of rigid porous solids, the IUPAC classification of adsorption isotherms [2] only considers hysteresis at large relative vapor pressure m.Actually, experimental evidence and theoretical reasoning supports the notion that sorption hysteresis should be restricted to the region of high relative vapor pressure.For non-rigid solids, however, the IUPAC classification for sorption isotherms appreciates the existence of low pressure hysteresis in materials containing micropores.As also mentioned in this IUPAC report [2], low pressure hysteresis is associated with swelling of non-rigid porous structures.The observed hysteresis loops may extend to the lowest attainable vapor pressure.
Low pressure hysteresis in some porous solids such as hardened cement paste and montmorillonite clay is thoroughly investigated experimentally and well illustrated in a large number of papers [3,4].These materials have a layered structure and a large fraction of slit-like micropores with a slit width comparable to or lower than a nanometer.Experimental results for hardened cement paste suggest that condensed water in very narrow slit-like pores produces the experimentally-observed pronounced sorption hysteresis, which is accompanied by swelling and shrinkage of this material [5][6][7][8].Surprisingly, although there is a widely-accepted plausible model for explaining low pressure hysteresis in hardened cement paste (Calcium-Silicate-Hydrate, C-S-H phase) by Feldman and Sereda [5], a corresponding mathematical description for water vapor sorption and swelling/shrinkage does not seem to exist.The approach of Pinson et al. [8] can well describe mesopore hysteresis at higher humidity, but the low pressure hysteresis has been considered only in the framework of a semi-empirical model.Some difficulties in modeling hardened cement pastes result from the complex pore network structure, which is comprised of micropores, mesopores and macropores.There are many proposals for structure models of the C-S-H phase [3,5,9], but the existence of layer stacks that have a local structure similar to the mineral tobermorite is generally accepted [7,10].Absorbed water layers (interlayer water) between the C-S-H sheets of the stacks can be desorbed and reabsorbed by changing the humidity.If the relative water vapor pressure or the humidity m is low, say 0.01 < m < 0.4, almost all condensed water (which can be desorbed without heating the cement paste) is stored as interlayer water between the C-S-H sheets [8].However, at higher humidity (m > 0.4), water condensed in mesopores and adsorbed on surfaces modifies the experimentally-obtained sorption and swelling/shrinkage isotherms.
Recently, we proposed a simple mathematical model to describe sorption, as well as swelling and shrinkage in porous solids with narrow slit pores [11][12][13].According to this model, narrow slits can be completely compressed by condensate tension at low relative vapor pressure of the ambient vapor (Figure 1).This compressed and closed configuration of a slit pore is expected to be stabilized by short-range adhesion forces.At high relative vapor pressure, when the compressive tension is weak, elastic restoring forces of the solid skeleton could be sufficiently strong to separate touching slit pore walls and recover the open pore configuration.Thus, compression of slits along the desorption branch and hampered recovery along the adsorption branch may be responsible for sorption hysteresis at low relative vapor pressure.In the framework of a simple macroscopic theory, adhesion is described by a surface energy that is non-zero for touching solid surfaces and equal to zero for separated surfaces [14].This simple macroscopic adhesion model requires extremely low relative vapor pressure for slit closing if the porous solid is not very soft.In this paper, we consider a more realistic hysteresis model by introducing an adhesion potential that is a smooth function of the slit width.The corresponding distance-dependent adhesion force supports pore closing, which becomes possible at a moderately low relative vapor pressure.In accordance with experimental results, extremely low relative vapor pressure is not needed for observing low pressure hysteresis.In comparison to a simple macroscopic theory, which only uses a single parameter for describing adhesion, the sorption isotherms based on the extended adhesion potential turn out to be modified, but the essential physical results remain unchanged.
Experimental data for swelling and shrinkage of hardened cement pastes can be compared with predictions of the theoretical model.An interesting feature is visible in some swelling/shrinkage isotherms, which display the humidity-dependent change of the mean distance between C-S-H sheets (basal spacing) of hardened cement pastes: upon drying, the average value of the basal spacing has been found to decrease discontinuously at a certain threshold of humidity [15].It will be shown that such a discontinuity is also predicted by the extended adhesion model with smooth adhesion potential.

Simple Macroscopic Model
Let us consider a narrow circular slit with radius R and thickness H embedded in a soft solid.The slit is assumed to be completely filled with condensate and connected to the ambient vapor outside the soft solid.This connection (which is not shown in Figure 1) could be a channel between the micropore and an adjacent mesopore, which itself is part of a mesopore system connected to the ambient vapor outside the porous solid.It is useful to introduce the relative vapor pressure m = p v /p s where p v is the vapor pressure of the ambient vapor and p s the saturation pressure of the bulk fluid.m may vary between a value above the threshold for capillary condensation and one.In this case, a slit pore is either filled with condensate or completely compressed.According to the classical theory of capillarity, the condensate exerts a compressive stress Π cap (m) = ρk B T ln (m) on pore walls, where ρ is the particle number density of the condensate, T the temperature and k B the Boltzmann constant.As a surprising result, Israelachvili [14] pointed out that this formula for Π cap (m) leads to formally correct results for adhesion forces between solid surfaces even at rather low relative vapor pressure (m = 0.1, ..., 0.2), when the corresponding meniscus radius turns out to be equal to a molecular diameter.In a previous publication [16], we have used the capillary pressure Π cap (m) (Kelvin's equation) to derive the compression of slits in soft materials.In very soft solids, complete slit compression can also occur if its thickness H is relatively large.Hardened cement pastes are not very soft materials.For these porous pastes, length change measurements [17] support the assumption that the thickness of very narrow water-filled slit-like spaces (called interlayer space in the cement literature) and the number of water molecules absorbed in these slits N are directly proportional (N ∝ H).Based on this proportionality, a more general approach leads to a solvation pressure, which is the sum of a humidity-dependent and a humidity-independent term (Appendix A): The humidity-dependent term: is proportional to a dimensionless constant factor q, which should not deviate much from one.
Here, we always assume that q = 1, which leads to the formal coincidence If the solid is soft and Young's modulus E is low, tension Π T is capable of compressing narrow slits.We define H 0 as the value of the slit width H at saturation pressure (m = 1), which is accompanied by vanishing condensate tension (Π (1) = 0).When the relative vapor pressure is reduced from m = 1 to a lower value, the tension can change the pore thickness from H 0 to H < H 0 .Reduction of the pore thickness produces elastic deformations in the solid skeleton of the soft solid.The elastic free energy of the deformations resulting from a compressed slit can be approximately evaluated by using the dislocation theory of solids [18,19].
A closed dislocation loop is formally produced by a Volterra process [18,20,21], which is used in the theory of edge dislocations.In the monograph [18], a circular dislocation loop is considered.The production of a circular dislocation loop by a Volterra process is comprised of two steps, where the first step produces the open pore.After the second step, the circular dislocation loop appears at the edge of the former circular slit pore: Step 1: An open circular slit is produced by removing all atoms of the solid in a disk-like region.The resulting structure appears as a circular slit with a thickness equal to the length b of Burgers vector, which corresponds to an open slit pore with thickness b.
Step 2: The opposite slit walls are pressed together, and the touching slit interfaces are glued.In the physics of solids, chemical bonds bridge both sides of the former slit-like cavity, but we assume that a compressive adhesion force hampers reopening of the pore.
This procedure can also be applied to planar circular slits in amorphous materials, where the slit thickness cannot be identified with the length of a Burgers vector (vector of a crystal lattice) because the elasticity theory is a macroscopic description of matter.
Analogously to the first step of the Volterra process, we can produce the open slit that is shown in the upper part of Figure 1.The slit with thickness H 0 = b and radius R is created by removing all molecules of the disk-shaped region.Then, the slit is filled with water at humidity m = 1.The contribution f (H) to the solvation force (Equation ( 1)) is zero for H = 0 because at saturation pressure (m = 1), adhesion appears only in the borderline case H = 0 (touching surfaces), as we only consider the simplest macroscopic description in this section (a force f (H) with a finite range ξ will be considered in Section 2.2).Hence, Π (m) (Equation ( 2)) defines the force (per unit area) between separated walls.As Π (m = 1) = 0, the water-filled cavity retains its thickness H 0 , and there are no deformations of the solid skeleton caused by slit compression.For this state (H = H 0 ), the elastic free energy F ela should be equal to zero, and the water-filled slit is in mechanical equilibrium ((∂F ela /∂H) H 0 ,R = 0).If humidity m is reduced, the tension Π (m) compresses the slit.Assuming a Hookean behavior, the elastic free energy (divided by the area A = πR 2 ) for the deformed solid is expected to be a quadratic function of the slit wall displacement H 0 − H.The ansatz: becomes quite plausible if F ela is approximately evaluated by using a Taylor series (Appendix B, Equation (A11)).An approximate expression for k e f f (H 0 , R) can be obtained by using the free energy per slit area F ela /A = k e f f (H 0 , R) H 2 0 /2 resulting from inserting H = 0 into Equation (3).This completely compressed slit configuration (H = 0) should appear when the compressive tension Π (m) is strong enough.On the other hand, the same closed slit configuration appears after the second step of the Volterra process to create a circular edge dislocation loop with b = H 0 for the magnitude of the Burgers vectors (Figure 1).The second step of the Volterra process is equivalent to slit closing.Hence, the elastic free energy for a circular edge dislocation loop F loop is equal to the elastic free energy F ela of a completely compressed slit.Equation F loop /A = F ela /A for H = 0 allows us to evaluate an approximate expression for the coefficient k e f f (H 0 , R) in Equation ( 3).As is demonstrated in the Appendix B, coefficient k e f f (H 0 , R) turns out to be independent of H 0 .Thus, we use the notation k (R) instead of k e f f (H 0 , R) in the subsequent equations of this paper.Then, the elastic free energy divided by the slit area A = πR 2 can be written as: where (Appendix B, Equation (A15)): E is Young's modulus, ν Poisson's ratio and radius R 0 is comparable to the core radius of the dislocation.For simplicity, we assume that R 0 H 0 .If the condensate tension is strong enough to completely compress the slit (H = 0), we obtain the free energy F ela /A = k (R) H 2 0 /2, which is equal to the elastic free energy of a circular dislocation loop (Figure 1).The closed slit configuration H = 0 is stabilized by adhesion of the touching slit walls.In previous papers [11][12][13], we have considered the simple macroscopic model for the adhesion potential where w is equal to the work per unit area to separate two touching planar surfaces at saturation pressure (m = 1).The difference of the grand potentials (divided by the slit area A) for the open and the closed slit configuration ω (m, can be written as [11]: The last term The potential difference ω (m, H) is equal to zero if there is an equilibrium between the open (H > 0) and the closed pore configuration (H = 0).If ω (m, H) < 0, the open slit is absolutely stable, and if ω (m, H) > 0, the closed configuration is the absolutely stable state.
Let us first briefly discuss the mechanism that leads to sorption hysteresis.Figure 2 illustrates the humidity-dependent width H of a narrow circular slit pore with radius R = 5 nm.At saturation pressure (m = 1), the value of the pore width H 0 is chosen to be 0.7 nm.The slit is embedded in a soft solid with Young's modulus E = 5 GPa, Poisson's ratio ν = 0.33 and adhesion energy per unit area w = 0.1 N/m.The number density of water is ρ = 33.5 nm −3 at the chosen temperature T = 300 K.The narrow slit is supposed to be always filled with condensate if H > 0. However, instead of evaporation, the amount of condensate in the pore can decrease by compression of the pore space.Condensate can be squeezed out of the slit pore, when the compression force (Equation ( 1)) increases with decreasing humidity m.This process occurs along the desorption branch in Figure 2 (Curve A → B → C).According to Equation (8), the slit thickness H is equal to H 0 at saturation pressure (m = 1), which refers to Point A in Figure 2.
After decreasing m to m C = 0.1, the slit is completely closed (Point C), and the onset of the short-range adhesion stabilizes the closed slit configuration (H = 0).This configuration remains closed along the line C → D as long as the difference of grand potentials ω(m, H) (Equation ( 7)) is positive.When the humidity m is increased to values exceeding m D = m B , the potential difference ω (m, H) is negative, and the closed slit pore can reopen.We evaluate the threshold for pore reopening from the condition ω (m D , H) = 0, which refers to the thermodynamic equilibrium between the open and the closed configuration.Thus, the adsorption branch starts at Point C, moves to Point D, jumps to Point B and finally arrives at Point A when m is further increased.The desorption-adsorption cycle (A → B → C → D → B → A) can be repeated many times without serious material damage, since mainly elastic deformations occur.In a region close to the dislocation line, however, the large strain may cause some plastic deformations.
For hardened cement pastes, experimental data for the average distance between adjacent C-S-H sheets (basal spacing) are available, and thus, the relation between m and H can be tested.Let us denote changes of the basal spacing by ∆H.Humidity-dependent changes of the slit width H are equal to ∆H.According to Equation ( 8), the change of the basal spacing ∆H by varying humidity m should be proportional to ln (m).The lower diagram in Figure 2 shows a semi-log plot for illustrating this logarithmic humidity dependence of H.In the plots for experimental data, the absolute vapor pressure p v is often chosen instead of the relative humidity m = p v /p s .Because ln (m) = ln (p v ) − ln (p s ), this different choice is not important.Using p v instead of m only causes a uniform shift of data points parallel to the humidity axis in semi-log plots.In Figure 2, the desorption branch A-C appears as a sloped straight line.Furthermore, the absorption branch consists of a straight line C-D parallel to the humidity axis and a perpendicular line D-B describing a discontinuity appearing when the closed slit reopens.This plot can be compared with analogous semi-log plots for experimental data obtained from X-ray measurements.For the branches A-C and C-D, Smith [7] found straight lines as depicted in Figure 2. Some deviations from the linear behavior at high relative humidity along the desorption branch can readily be explained by water-filled mesopores above the threshold of capillary condensation, because the condensed water in mesopores exerts a compression force on the solid skeleton [1,16].Finally, the experimentally-observed data reveal a rapid change of the basal spacing after pore reopening.The experimental curves are relatively steep, but not perpendicular to the humidity axis, while in Figure 2, the depicted straight line D-B is perpendicular to the humidity axis.This deviation could be explained by a slit size distribution, as illustrated in a previous paper [11].It should be noted, however, that some publications report on the appearance of jumps along the desorption branch in semi-log plots of the basal spacing versus humidity [15,22].The macroscopic model that uses the discontinuous adhesion potential (Equation ( 6)) is not capable of describing a discontinuity of the basal spacing for the desorption branch.In this paper, a more realistic adhesion potential that has a finite interaction range is used.We want to check if this extension is capable of describing a discontinuous dependence of H on m.

Extended Theoretical Model
The considered model for the adhesion potential (Equation ( 6)) is only a proper approximation in cases where the interaction force f (H) = −dW (H) /dH per unit area decays rapidly with increasing slit wall distance H.For elucidating the impact of a finite interaction range on sorption diagrams, we can use the force: where the length ξ is a measure of the interaction distance and A is the magnitude of f (H).
The choice A > 0 leads to a compression force accompanied by a cohesion of slit walls.The constants A and ξ in Equation ( 9) could be fitted to results of Monte Carlo simulations.In hardened cement paste, the cohesion force between C-S-H sheets is supposed to originate from calcium ions in the interlayer water [23].Equation ( 7) is replaced by (Appendix A, Equation (A10)): where h 0 is a reference width of the slit at vanishing solvation pressure Π T introduced in the Appendix A. It is expedient to define the reference width as the value of H at vapor saturation (m = 1).This width would be equal to h 0 if condition f (H = h 0 ) = 0 were satisfied.Usually, condition f (H = h 0 ) = 0 is not satisfied because the attraction force between walls is not equal to zero, even at vapor saturation (m = 1).Actually, the condition for mechanical stability of an open slit ∂ω (m, h) /∂H = 0 leads to: At saturation pressure of the vapor, we have Π (m = 1) = 0, and according to Equation ( 11), the equilibrium slit thickness deviates from h 0 if f (H = h 0 ) = 0.It is expedient to introduce a slightly shifted slit thickness H 0 by the equation: It should be noted, however, that the difference between H 0 and h 0 is rather small, as the force magnitude | f (H)| is a strongly decreasing function of the wall distance H. Thus, the condition h 0 holds.Replacing h 0 in Equation ( 10) by using Equation (12), ω (m, H) is transformed into: where an adhesion potential: is introduced, and f (H 0 ) is a constant, independent of m.Now, the condition for mechanical stability ∂ω (m, H) /∂H = 0 leads to: This equation allows us to evaluate how the width H of an open slit depends on the relative vapor pressure m.Obviously, Equation ( 15) satisfies the condition H (m = 1) = H 0 , and thus, the value of the wall width H at saturation pressure (m = 1) is equal to H 0 .In a similar way, as has been done previously, we define a work per unit area w necessary to separate two touching planar surfaces and to restore the original slit width H 0 at saturation pressure (m = 1).According to Equation (14), this work can be evaluated by w = W (H = H 0 ) or: Combining Equations ( 9) and ( 16) yields: where: 17)) into Equation ( 14), we obtain the adhesion potential: Figure 3 illustrates the adhesion potential W (H) for three different interaction distances ξ.If ξ → 0, the force f (H) deviates from zero only in the case of touching wall surfaces (H = 0), as in the simplest macroscopic adhesion model (Equation ( 6)).

Mechanical Stability of the Open and the Closed Slit Configuration
In the previous section, we have introduced the difference ω (m, H) of grand potentials for deciding whether an open (H > 0) or closed slit configuration (H = 0) has a lower grand potential.Checking whether or not ω (m, H) has a negative value provides us the preferred state of the slit, but mechanical stability is not guaranteed.Let us first check the stability of the open slit configuration (H > 0).Mechanical stability requires that the grand potential Ω as a function of H has a local minimum.Apart from the equilibrium condition A −1 ∂Ω/∂H = 0, the condition for a local minimum of the grand potential A −1 ∂ 2 Ω/∂H 2 > 0 must be satisfied.Hence, taking into account that the grand potential (divided by the constant slit area A = πR 2 ) can be written as: we obtain the stability condition for an open slit where the width H > 0 is evaluated by Equation ( 15).The point on the desorption branch where the open slit configuration loses its stability is the critical point of a fold in the framework of catastrophe theory [24].At this point, the open pore collapses on the desorption branch, and the closed slit configuration (H = 0) is the only stable state at low humidity (Figure 4).Unstable branches of the fold catastrophe are not depicted in the diagrams shown in Figure 4.In the borderline case ξ → 0, the extended model for sorption hysteresis based on a continuous function W (H) (Equation ( 18)) reduces to the simpler model with a discontinuous adhesion potential (Equation ( 6)).Taking into account that k(R) > 0 and: for H > 0 and arbitrary values of n, the stability condition for open pores (Equation ( 20) turns out to be always satisfied in the case of the discontinuous adhesion model (Equation ( 6)).For confirming the consistency of the considered model, we also check the mechanical stability of the closed slit configuration (H = 0).Stability considerations for the closed slit state are more subtle because the width H must be always positive (H ≥ 0).This constraint for H has the consequence that the condition A −1 ∂ 2 Ω/∂H 2 H=0 > 0 is not necessary for the closed slit (H = 0) [24].Instead, the closed slit configuration is stable if condition A −1 (∂Ω/∂H) H=0 > 0 holds.From a physical point of view, this condition allows for a compressive stress P = −A −1 (∂Ω/∂H) H=0 < 0 in the plane bounded by the circular dislocation loop (Figure 1).Using Equation ( 19), condition A −1 (∂Ω/∂H) H=0 > 0 yields: Obviously, for sufficiently small values of the interaction length ξ, the condition for local stability (Equation ( 22)) of closed slits is satisfied in the whole range of relative vapor pressure (or relative humidity) 0 < m ≤ 1.Using the stability conditions (Equations ( 20) and ( 22)), the mechanical stability of all states defined by the isotherms illustrated in this paper can also be checked.

Results and Discussion
Sorption hysteresis at moderate and low humidity can be explained by strong cohesion forces, which attract the walls of narrow slit pores [11][12][13].The presented model suggests a possible mechanism, which leads to a sorption hysteresis loop accompanied by a related loop of the pore volume change.Low pressure hysteresis can be based on pore closing by compressive stress (solvation force) and a stabilization of the closed pore configuration by short-range adhesion.In previous papers [11][12][13], we have used a primitive discontinuous adhesion potential, which is zero if slit walls are touching each other and jumps to a constant non-zero value if the slit is open and the walls are separated from each other.For moderately rigid materials with slit pores, such as hardened cement paste, this simple model leads to a very low humidity threshold necessary for pore closing, which is a precondition for the appearance of large hysteresis loops.However, experimentally-observed sorption isotherms for hardened cement paste show large hysteresis loops even in those cases where the humidity is not very small [3].Thus, the considered simple model for sorption hysteresis seems to be not applicable to porous materials that are not very soft.This disadvantage of the sorption model can be remedied by introducing a smooth adhesion potential, which produces a compression force between separated slit walls and supports pore closing.As an example, we have used an exponentially-decaying adhesion potential, which is suitable to define an interaction length ξ.If ξ tends to zero, the extended model with continuous adhesion potential reduces to the simple model with discontinuous adhesion potential.Figure 4 shows sorption isotherms with the hysteresis loop for four different interaction lengths ξ.The first diagram (a) depicts the case of a very low value for the interaction length (ξ = 0.001 nm).In this case, isotherms evaluated by the continuous and the discontinuous adhesion model are indistinguishable.The slit wall distance H (and thus, the adsorbate volume) is a monotone continuous function of the relative humidity m (Figure 4a).However, if ξ is distinctly larger than zero, the slit width H jumps from a finite value to zero on the desorption branch.As shown in Diagrams b, c and d, the humidity threshold for pore closing, where H jumps from a finite value to zero on the desorption branch, strongly increases with increasing ξ.In contrast to the threshold for slit closing on the desorption branch, the humidity threshold for pore reopening on the adsorption branch is almost constant if ξ increases.Obviously, the more general continuous adhesion model can lead to pore closing at a much higher relative vapor pressure than the discontinuous model.A smoothly-varying adhesion potential avoids extraordinarily small values of relative vapor pressure required for slit pore closing.Therefore, we expect that extended models with a smooth adhesion potential are capable of providing us a better description of the sorption hysteresis for porous solids if they are not very soft.
Tests of the model can be done by comparing theoretical results with available experimental results for hardened cement pastes.Some theoretical and experimental sorption isotherms have been compared in a previous publication [11].The comparison of sorption isotherms is difficult because hardened cement pastes have a complex pore structure.Apart from micropores, the pastes also have mesopores, which contain a large amount of water at high humidity m.Thus, analyzing the humidity dependence of the basal spacing may provide us a better interpretation, since X-ray experiments only probe micropores.The basal spacing changes due to humidity-dependent absorption and desorption of water layers between the C-S-H sheets.The change of the basal spacing ∆H is equal to the change of the slit width H.According to Equation (8), the change ∆H of the basal spacing should obey a logarithmic law (∆H ∝ ln (m)) on the desorption branch, as shown in the lower diagram of Figure 2.This logarithmic dependence is confirmed by semi-log plots of experimental data by Smith [7].In other presumably differently-prepared hardened cement pastes, the change of the basal spacing ∆H is discontinuous at a certain water vapor pressure [15].Our extended model allows us to describe such a jump in the desorption branch, as illustrated in Figure 4.However, more experimental data are needed for a comprehensive test of the simple mesoscopic model.the slit.At vanishing solvation force, if Π T were equal to zero, the slit thickness H would be equal to the reference value h 0 .A slit can be compressed by the solvation tension Π T exerted on the slit walls.The elastic free energy F ela of the solid depends on the deviation H− h 0 of the slit width H from its reference value h 0 at Π T = 0. Hook's law in elasticity theory suggests that a Taylor series of F ela up to second order terms should be an accurate approximation.Thus, we obtain: where k e f f (h 0 , R) = ∂ 2 F ela /∂H 2 h 0 ,R and the Landau symbol O (H − h 0 ) 3 indicates that higher order terms are neglected.In the reference state H = h 0 for Π T = 0, the elastic free energy F ela must have its minimum value (mechanical equilibrium).Hence, condition (∂F ela /∂H) h 0 ,R = 0 is satisfied, and thus, Equation (A11) can be simplified to: where we have omitted the constant term F ela (h 0 , R). Coefficient k e f f (h 0 , R) could be considered as an empirical quantity that must be fitted to experiments.An approximate expression for k e f f (h 0 , R) can be obtained from the dislocation theory.A strong solvation force may squeeze out the water of a narrow slit, so that its thickness H becomes zero.As discussed in Section 2.1, the completely compressed slit pore (H = 0) in the solid material can be considered as a circular edge dislocation loop with Burgers vector − → b satisfying the relation h 0 = − → b = b.According to Hull [18] (Chapter 5, Equation (5.10)), an approximate expression for the elastic free energy of a circular edge dislocation loop is given by: where G denotes the shear modulus and R 0 is a fixed radius comparable to the core radius of the dislocation.This formula was also published in a monograph on rheophysics [19].Taking into account the relation E = 2G/ (1 + ν) between Young's modulus and shear modulus, the elastic free energy per loop area A = πR 2 is written as: Now, the unknown coefficient k e f f (h 0 , R) in the elastic free energy F ela = k e f f (h 0 , R) (h 0 − H) 2 /2 can be adjusted by comparing F ela /A = k e f f (h 0 , R) h 2 0 /2 (closed slit, H = 0) and F loop /A (circular edge dislocation with magnitude b = h 0 of the Burgers vector).If H = 0 and b = h 0 , equation F ela /A = F loop /A leads to a function k e f f (h 0 , R) = k(R), which does not depend on h 0 : Finally, we note that Equation (A15) for coefficient k (R) is not changed if Equation ( 12) is used to replace h 0 by H 0 (reference slit width H 0 for humidity m = 1 instead of h 0 for vanishing solvation tension Π T = 0).In the case of the macroscopic description of adhesion (Section 2.1), the humidity-independent part of the solvation force f (H) is non-zero only if H = 0 (completely closed slit).Then, Equation ( 12) leads to h 0 = H 0 and rather simple equations (Equations ( 4), (7) and ( 8)) in Section 2.1.

Figure 1 .
Figure 1.Cross-section of a narrow circular slit pore with radius R and width H 0 is embedded in a soft solid filled with condensate at saturation pressure (m = 1).At low values of m, the slit is completely compressed (H = 0) and degenerates to a closed dislocation loop.The dark triangles symbolize the cross-section of the dislocation core, where the solid material is strongly compressed.Connections to larger pores or to the ambient vapor outside the porous solid are not depicted.
2 is equal to the elastic free energy F ela /A for H = 0.The condition for mechanical stability ∂Ω open /∂H = ∂ω (m, H) /∂H = 0 for open pores (H > 0) yields: