An elastic sphere adheres to a rigid substrate in the presence of moisture. The adhesion–detachment trajectory is derived based on the Hertz contact theory that governs the contact mechanics and Laplace–Kelvin equation that governs the water meniscus at the interface. The intersurface attraction is solely provided by the Laplace pressure within the meniscus. Interrelation between the applied load, contact radius, and approach distance is derived based on a force balance. The resulting “pulloff” force to detach the sphere exceeds the critical load in the Derjaguin–Muller–Toporov (DMT) limit which only holds at saturated moisture. The new model accounts for the finite size of water molecules that is missing in virtually all classical models.

## Introduction

Adhesion of an elastic sphere onto another sphere or a rigid substrate is ubiquitous in many branches of science and technology such as colloidal particles, storage and transportation of glass powder, and micro-electromechanical systems especially when the typical dimension shrinks to submicron scale [1–3]. Exposure to moist air leads to meniscus formation at the contact interface and could post serious problems of operation hindrance and reliability.

Meticulous experimental work conducted by Christenson [4,5], Maugis and Gauthier-Manuel [6], and Xu et al. [7], and theoretical framework constructed by Fogden and White [8], and Grobelny et al. [9] demonstrated the effect of moisture on adhesion. In a relatively dry environment, the inevitable meniscus shrinks to a small dimension, and the resulting short range Laplace pressure leads to the Johnson–Kendall–Roberts (JKR) limit. Close to saturated moisture, the meniscus grows to such a large extent that adhesion mechanics approaches the Derjaguin–Muller–Toporov (DMT) limit. Wan and Lawn [10] and Wan et al. [11] investigated adhesion of two noninteracting elastic cantilevers in the presence of meniscus, and found that the critical energy release rate is twice the surface tension of water as if the crack front is right at the meniscus. The intermediate range of relative humidity falls into the JKR–DMT transition and is accounted for by the Tabor's parameter [12]. In such classical work, the “pulloff” force leading to spontaneous detachment is taken to be *F** = −*χ* × *πR* × 2*γ* cos*θ*, where *R* is the sphere radius, *γ* is the surface tension of water, *θ* is the contact angle at meniscus–solid interface, and *χ* is a numerical constant *χ*_{JKR} ≤ *χ* ≤ *χ*_{DMT} with *χ*_{JKR} = 3/2 and *χ*_{DMT} = 2. Therefore, *F** is a monotonic increasing function of relative humidity approaching the JKR and DMT limits in the extremes. Zeng and Streator [13] and Xue and Polycarpou [14,15] started from a long-range intersurface attraction and derived the sphere deformation. The numerical approach interestingly suggests that *F** is a monotonically decreasing function of humidity as it approaches the DMT limit.

In this paper, we will build an alternative model based on the Hertz contact theory [2] and cohesive zone model [16], rather similar to the DMT model. We will first present the theory and results, leaving the discussion and comparison with the classical models to Sec. 4. In the absence of moisture, the intrinsic intersurface forces such as van der Waals are present, but such interactions are taken as the baseline or ignored hereafter to avoid unnecessary mathematical complications. The total force acting on the sphere, therefore, becomes the sum of the applied load and the Laplace pressure within the meniscus. The compressive stress distribution within the contact and the deformed geometry strictly follow the classical Hertz contact theory [17]. The Dugdale–Barenblatt–Maugis cohesive zone model is adapted where the disjoining pressure is taken to be constant over a finite range determined by the Laplace–Kelvin equation [10,11,16,18]. The measurable relations of applied load, contact radius, and approach distance are derived.

## Theory

*F*, is applied to an elastic sphere with radius

*R*, elastic modulus,

*E*, and Poisson's ratio,

*v*, resulting in an approach distance,

*δ*, and contact radius,

*a*. In the absence of adhesion, the classical Hertz contact theory [17] gives the interrelation of (

*F*,

*a*,

*δ*)

*F*

_{H}=

*F*, governed by Eq. (2). In the presence of a general intersurface interaction potential ϕ(

*z*), a disjoining pressure

*p*=

*d*ϕ(

*z*)/

*dz*arises at the sphere–substrate interface, and the total adhesion force acting on the sphere becomes

*p**, with a range,

*z**, over a cohesive zone immediately outside the contact circle (Fig. 1), such that

*z*

_{0}is the equilibrium sphere–substrate separation to account for the finite dimension of molecular planes. The annular cohesive zone is bounded by

*a*<

*r*≤

*a*

_{1}with

*z*(

*r*=

*a*

_{1}) =

*z**. Substituting Eq. (6) into Eq. (4), the effective adhesion force becomes

*γ*=

*p** × (

*z** −

*z*

_{0}) ≈

*p** ×

*z**, since

*z** ≫

*z*

_{0}. It is apparent that

*a*

_{1}depends on both

*z** which depends on the interface chemistry and

*z*(

*r*) which depends on the elastic deformation. The measurable relations,

*a*(

*F*) and

*a*(

*δ*), can then be solved based on Eqs. (1)–(7).

*p** =

*γ*/

*r*

_{k}, with

*γ*the surface tension (or, equivalently, surface energy) of water. Rather than accounting for the sophisticated principal curvatures of the meniscus and the geometry of the deformed sphere, the force range is taken to be approximately the meniscus diameter or twice the Kelvin radius [10], $rk$, or

*θ*≈ 0 deg, molecular volume of water

*v*, Boltzmann's constant

_{m}*k*, temperature

*T*≈ 300 K, and relative humidity RH. At the saturation limit (RH ≈ 100%),

*p** → 0 and

*z** → ∞. Without the loss of generality,

*θ*is taken to be zero hereafter. In case of a nonzero

*θ*, all

*γ*term will be replaced by

*γ*cos

*θ*. In dry conditions with RH = 0%,

*p** → ∞ and

*z**→ 0. Using Eq. (7), the adhesive force is given by

## Results

Figure 2(a) shows *a*(*F*) for a range of RH for a sphere with *E* = 1.0 GPa. For RH = 5%, the detachment trajectory follows path ABCDO (dark curve). At A when *F* = 0, a nonzero contact radius is expected. Increasing tensile load (*F* < 0) reduces the contact circle along ABC until the sphere spontaneously detaches from the substrate or “pulloff” at C with a nonzero radius (*a*^{†} > 0). The “pulloff” at *F*^{†} occurs under fixed load when *dF/da* → ∞ or *da/dF* = 0. The branch CD with *dF/da* < 0 is physically inaccessible under fixed load, because *F*^{†} is already the minimum force along this curve, but possible under displacement control or fixed grips. The contact shrinks further along CD (*dark curve*). At *D*, the contact is reduced to a point with *a* = 0 and the solid resumes its spherical shape. Along DO, the sphere is readily out of intimate contact with the substrate, but the meniscus turns into a water bridge still exerting a traction on the sphere. Increasing tension reduces the bridge width until it breaks at O, the sphere then completely breaks free from the substrate at *F*^{††} = 0. Point O is coined “pulloff” under fixed grips. The detachment trajectory is similar for other moisture level. In the dry conditions with RH ≈ 0%, *r*_{k} → 0, *p** → ∞, the cohesive zone vanishes, and a large pulloff force |*F*^{†}| is predicted due to a constant *γ*. As RH increases, |*F*^{†}| approaches the limit of *F*^{†} = −4*πRγ* along CD (*gray curve*) and *a*^{†} approaches zero, which are summarized in Fig. 2. It is noted that the model fails when 2*r*_{k} is of the same order of magnitude as *R*.

*F*(

*δ*) for a range of RH. For RH = 5%,

*δ*is nonzero at

*F*= 0 due to adhesion (not shown). Applied tension (

*F*< 0) further reduces

*δ*until “pulloff” under fixed load occurs at C with maximum tensile load with

*dF*/

*dδ*= 0. Under fixed grips,

*δ*can decrease further under stable equilibrium. At D,

*δ*= 0 and the sphere resumes its spherical shape. Here the cohesive zone extends to the meniscus with

*a*

_{1}= 2 (

*R × r*

_{k})

^{1/2}, and the external force to maintain equilibrium becomes

*F*= −

*F*

_{ad}= −

*p** ×

*π a*

_{1}

^{2}= −4

*πRγ*. Further decrease in

*δ*(< 0) detaches the sphere from the substrate, but the sphere is linked by a bridge to the substrate. To maintain mechanical equilibrium along DH, the applied load is balanced by the adhesion force

where *a*_{1} = (2*R* (*δ* + 2*r*_{k}))^{1/2} is the bridge radius found by simple geometry. As the sphere moves further out from the substrate, the water bridge shrinks and the external load decreases linearly with *δ*. At *H*, *δ*^{††} = −2*r*_{k}, *a*_{1} = 0, the meniscus bridge collapses, the applied tension vanishes, and “pulloff” occurs at *F*^{††} = 0. The fixed grips “pulloff” locus is along the negative *δ*-axis, corresponding to O in Fig. 2(a). At saturation (RH → 100%), *F*^{††} → −4*πRγ* and *δ*^{††} → −∞ theoretically as *r _{k}* → ∞ in Eq. (8). However, the water pillar collapses before

*δ*reaches −

*R*due to the geometrical constraints.

Figure 3 shows *F*^{†}, *a*^{†}, and *δ*^{†} as functions of RH. In Fig. 3(a), the calculated *F*^{†}(*a* = 0, *δ* = 0) slightly deviates from the DMT limit of −4*πRγ* at very high RH as the cohesive zone now extends to the sphere dimension (*a*_{1} ≈ *R*), and the Hertz assumption breaks down. Figure 3(b) shows the expected monotonic decreasing function of *a*^{†} (RH) that approaches *a*^{†} = 0 in saturated moisture. Figure 3(c) shows *δ*^{†} (>0) under fixed load and *δ*^{††} (<0) under fixed grips as RH varies. A drastic decrease in *δ*^{††} is expected as RH exceeds 50% since the large meniscus (*r*_{k}) sustains a long water bridge (c.f. Eq. (8)).

## Discussion

The new model is rigorously compared with the classical adhesion models: (i) Bradley's model [1] where the sphere is rigid and nondeformable, (ii) DMT theory [19] for small but stiff spheres in the presence of a weak force with an ideal infinite range, and (iii) the JKR theory [20] for large but soft spheres in the presence of a strong surface attraction with an ideal zero range.

In Bradley's model, a rigid sphere makes only a single point contact with the substrate with *δ* = 0 and *a* = 0 independent of both RH and *F*. In the presence of a meniscus, the sum of the applied load and Laplace pressure is balanced by the reaction force at the point contact in the form of a delta function such that the adhesion force is a constant with *F*_{ad} = 4*πRγ* for *F* ≥ 0. When *δ* turns negative, the sphere is out of intimate contact with the substrate but is linked by a water pillar. Equation (10) holds, and the behavior is identical to the present model. Pulloff occurs when the water pillar breaks.

In the classical JKR model [12], the disjoining pressure is taken to be infinite with a vanishingly small force range though the energy to bring the two adhering surfaces together is finite. The tensile stress at the contact edge is therefore theoretically infinite and is capable of deforming the sufficiently soft sphere to the Griffith parabola with (∂*y*/∂*r*)* _{r = a}* → ∞. An energy balance predicts “pulloff” under fixed load at

*P*in Fig. 2, where

*F*

^{†}= −3

*πRγ*,

*a*

^{†}= (9

*πR*

^{2}

*γ*(1 −

*v*

^{2})/4

*E*)

^{1/3}, and

*δ*

^{†}= −(

*a*

^{†})

^{2}/3

*R*, whereas “pulloff” under fixed grips occurs at

*P′*where

*F*

^{††}= −(5/12)

*πRγ*,

*a*

^{††}= (

*πR*

^{2}

*γ*(1 −

*v*

^{2})/4

*E*)

^{1/3}, and

*δ*

^{††}= −3(

*a*

^{†})

^{2}/

*R*. A close examination of these underlying assumptions prompts one to check for limitations. Laplace pressure cannot increase indefinitely but has an upper bound of the hydrogen bonds of water in the molecular scale or vapor pressure in the macroscopic scale. It is unable to deform a sufficiently stiff sphere (e.g., silica glass) to the ideal parabolic shape at the contact edge, since hydrogen bonds will be broken and water molecules dissociate beforehand. The present model predicts

*a*

^{†}> 0 under fixed load, but since there is no local deformation at the contact edge as in JKR. Both

*F*

^{†}and

*a*

^{†}are distinctly different from the JKR values. It is noted that Hertz contact theory in Eqs. (1)–(3) is adopted in this work for stiff solid sphere. Should a large and soft sphere be considered, the deformed profile will significantly deviated from Hertz, and a better computational model (e.g., Refs. [6] and [12]) will be necessary for the thermodynamic energy balance.

The present model is similar to the DMT model to some extent; in that, the sphere deforms according to the Hertz theory, and the adhesion mechanics serves as the limit for RH approaching saturation as shown in Fig. 2(b). Maugis's version of DMT assumes an infinite force range [21] and *a*_{1} → ∞. “Pulloff” under fixed grips occurs at *δ* → ∞, and the external tension to maintain the infinite water pillar remains constant at *F*^{†} = −4*πRγ* for all *δ* ≤ 0. In the absence of sphere–substrate interaction, the “crack front” is located literally at the meniscus, which is consistent with the classical experiment in measuring surface tension of a liquid trapped between two cantilevers [10,11]. One interesting implication of DMT is that the water molecules must penetrate all the way from the meniscus *z*(*r* = *a*_{1}) = 2*r*_{k} to the contact edge at *z*(*r = a*) = *z*_{0} so that the Laplace pressure is present in the entire stretch of the cohesive zone. In reality, the physical dimension of two water molecules bonded by a typical hydrogen bond is roughly 8.4 Å. A *vacuum* thus exists between the contact edge and the molecular wedge, or, *p* = 0 in *r*(*z*_{0}) to *r*(*z* = 8.4 Å). In the limit of RH = 100%, *a*_{1} ≫ *r*(*z* = 8.4 Å), and our model approaches DMT.

In the intermediate range of RH, the new model requires the Laplace pressure to be bounded by the cohesive zone, *p* = *p** in *r*(*z* = 8.4 Å) to *r*(*z* = 2*r*_{k}). Based on Eq. (8), 2*r*_{k} = 8.4 Å at RH ≈ 30%, and no meniscus can exist below this critical RH. The corresponding Laplace pressure is too weak to deform the contact edge to the Griffith parabola. The contact edge can now be arbitrarily chosen to be in the small range between *r*(*z*_{0}) to *r*(*z* = 8.4 Å). Without loss of generality, we choose to assume the Hertzian compressive stress within the contact circle and the Laplace pressure without, regardless the exact location of the contact edge. To further justify the new model, Fig. 4 shows the deformed profile for *E* = 10 GPa and a range of RH according to Maugis's cohesive zone model along with the JKR and DMT limits. All curves show a characteristic kink, where the slope ∂*z*/∂*r* changes abruptly indicating the cohesive zone edge *r* = *a*_{1}. At RH < 30%, 2*r*_{k} < 8.4 Å and intersurface force does not exist. Therefore, meniscus alone is unable to lead to the JKR limit. In Maugis's calculation, the compressive contact stress gradually turns tensile radially outward and reaches the maximum tension matching the Laplace pressure. However, the absence of intrinsic solid–solid interaction and the finite water molecular dimension do *not* allow tension within the contact. Our model allows a discontinuity in stress at the contact edge where it abruptly drops from zero based on Hertz to the full strength of Laplace pressure.

It is worthwhile to mention a few existing models in the literature. Christenson [5] measured an *increasing* pulloff force as RH rises as shown in Fig. 5. The ostensible contradiction seems to be explainable by the van der Waals attraction at the solid–solid interface that is inevitable and overshadows the Laplace pressure. Xue and Polycarpou [14] used a numerical approach to show *F*^{†} as a monotonic *decreasing* function of RH > 70% that approaches the DMT limit at saturated moisture. Such results match with the present model quite well, though Polycarpou only focused on the wet environment with RH > 70% only to avoid the intrinsic solid–solid interaction. We calculated the data points for RH < 70% based Polycarpou's data and Maugis' model as suggested by Xue and Polycarpou [14], which do not match quite well with the present model. Based on the JKR model, Fogden and White [8] predicted the large and soft spheres at relatively low RH pull off at *F*^{†} → −3*πRγ*, shown in the Fig. 5. The size of water molecule was not considered in their model. In the presence of a strong sphere–substrate interaction, the meniscus becomes negligible, and the adhesion–detachment mechanics approaches the JKR limit, which is the beyond the scope of this work. Another set of data and theoretical model worth discussion is that due to Grobelny et al. [9]. The capillary force, equivalent to *F*_{ad} (*a*) in Eq. (7), is deduced based on the Kelvin radius, Laplace pressure, and Hertz contact theory, and serves as the critical threshold for pulloff. However, the *repulsive* reaction force at the contact is altogether ignored, or equivalently, *F*_{H} = 0 in Eq. (5) and *F*^{†} = *F*_{ad}. Since the force balance is offset, the mechanical response *F*(*a*,*δ*) and the subsequent *F*^{†} and *F*^{††} are different from the present work. Despite such inconsistencies, the experimental data in Ref. [9] are presented for gold microsphere with *E*_{gold} ≈ 78 GPa, which deviate significantly from trends in this work and others. However, the theoretical model for RH = 45% coincides with our model.

The classical Tabor parameter, *μ* ∝ (*R γ*^{2}/*E*^{2}*z*_{0}^{2})^{1/3} or (*R p**^{2}/*E*^{2})^{1/3}, and its various forms, are used to discuss the JKR to DMT transition in the literature such that small *μ* leads to DMT-limit and large *μ* to JKR-limit [6,7,22,23] including pulloff force as a function of RH. In this work, *μ* is not an appropriate parameter, because the deformed geometry of the sphere is established solely on the classical Hertz contact theory (c.f. Eq. (4)), and the local deformation at the contact edge similar to the JKR model is ignored. A comprehensive model will have to solve the mathematically involved governing solid mechanics equations, which is beyond the scope of this paper. For stiff sphere like silica and glass, a full-fledged computation unlikely leads to a deformed geometry significantly different from the Hertz model. This work aims to demonstrate how a proper description of the Laplace pressure leads to an adhesion mechanics distinctly different from the classical JKR and DMT model.

As a final remark, implications of the new model is significant. For instance, it is known that heating and drying helps preventing coagulation in powder (e.g., glass beads, dry spray) and stiction in microdevices. This work shows that once RH drops to roughly 30%, water condensation, and the consequent adhesion in glass become relatively unimportant.

## Conclusion

A simple adhesion model is derived for an elastic sphere adhering to a rigid planar substrate in the presence of moisture. The model is useful in discussing adhesion of microscopic objects, and the derived values of “pulloff” force and contact radius can be used to deduce the magnitude and range of intersurface forces. Contrary to the classical description based on Maugis's JKR–DMT transition, meniscus alone causes |*F*^{†}| to exceed 2*πRγ* for any relative humidity and the DMT limit only holds at saturation.

## Acknowledgment

This work was supported by National Science Foundation. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.

## Funding Data

Directorate for Engineering (CMMI No. #1333889).