Contact angle dependence of solid probe–liquid drop forces in AFM measurements
Introduction
Following the early work of Ducker et al. [4], atomic force microscopy (AFM) has been applied to direct measurement of forces between gas bubbles [5], [6], [7] or liquid drops [8], [9], [10], [11], [12] and solid colloidal probe particles.
The theoretical challenge posed by these measurements is to allow, in the force calculation, for deformation of the interface, itself dependent on the induced pressure distribution. The problem has received attention along numerical [12], [13] and analytical [1], [2], [3], [14], [15] lines.
This paper addresses the interaction of a solid colloidal probe with an incompressible liquid drop, however as has been noted elsewhere [2], [13], the finite compressibility correction for gas bubbles is negligible in typical experimental situations. Thus the treatment here applies equally to bubbles.
Numerical solutions of the augmented Young–Laplace equation were used by Bhatt et al. [13] to obtain interface deformation and force–distance curves. An iteration on the central drop height and the pressure change in the drop was used to simultaneously implement fixed contact line and constant volume constraints. Similar numerical solutions were examined by Aston and Berg [12], in conjunction with experimental results.
Difficulties arise, however, in interpreting these numerical results, since the central drop height is sensitive to the small change in drop liquid pressure, requiring evaluation of the drop volume to very high numerical resolution [15].
Earlier numerical work focussed on bubble–bubble [16] and liquid drop–planar solid [17], [18] interactions. Bachmann and Mikavcic [18] assume the drop contact line is pinned to a capillary. Miklavcic et al. [17] make a similar assumption, but also discuss an example where the drop extends across a substrate and the contact angle θ is governed by Young’s equation,i.e. the constant θ condition, where γSL and γSA are respectively the surface free energy densities of the substrate–liquid drop interface and the substrate interface with the bathing aqueous solution, and γ is the droplet surface tension.
Young’s equation (1) will apply to real surfaces (i.e. in the presence of physical and chemical defects), provided the defects are sufficiently weak [19]. The presence of strong defects will lead to pinning [19] of the contact line, corresponding to contact angle hysteresis, and requires the use of fixed contact line boundary conditions.
Miklavcic [14] considered the interaction via electrical double-layer forces of a deformable droplet (pinned to a substrate) with a solid sphere, using perturbation theory to obtain the first order deformations away from sphericity. This approach was found to be appropriate for cases where the Debye length scale is at most a few orders of magnitude smaller than the unperturbed droplet radius of curvature, but it was noted that for cases where the droplet radius is many orders of magnitude larger, a matched asymptotic expansion would be more suitable.
This suggestion was implemented in Chan et al. [1] for purely repulsive potentials and fixed contact line, for the drop geometry and also for a spherical probe interacting with an (initially) planar deformable substrate. The matched expansion has the advantage that the Young–Laplace equation is solved analytically over the vast majority of the drop profile, alleviating the difficulties associated with purely numerical calculations. The method has been extended to potentials with attractive components by Dagastine and White [15].
Attard and Miklavcic [2] consider the problem under fixed contact angle assumptions in the weak force regime, also using a matched expansion approach.
In analyzing colloidal interactions between solid AFM probe particles and liquid drops, we wish to obtain the total force F(X) as a function of distance between the probe and the rigid substrate (i.e. the AFM stage). A crucial element in the calculation is the small change δR in the mean radius of curvature of the drop, which is proportional to the small change in liquid pressure inside the dropThese perturbation quantities depend on the total force F, the contact angle θ made by the unperturbed drop at the substrate, and the nature of the boundary conditions imposed at the contact line.
In experiments such as those of Hartley et al. [11], where the drop is deposited onto a uniform substrate, the appropriate boundary conditions will depend on the surface defect structure at the contact line, which can of course be controlled [20]. We will see that the motion of a free contact line (implied by fixing θ) would typically occur on a submicron scale in AFM experiments.
The probe-AFM stage separation can be expressed [1] as a function of the central liquid drop–probe separation D0,where a is the probe particle radius, H and G are certain integrals of the disjoining pressure over the (numerically calculated) drop profile in the inner (near to probe) region, z0∞ is the central height of the unperturbed drop and θ is the contact angle at the base of the unperturbed drop, which has curvature R0−1. The total force exerted by the probe is 2πγG, so that (3) provides the force–distance relationship. Typical length scales for AFM experiments [11], [12] are a∼5–10 μm and R0∼0.25–0.5 mm. From the maximum force values and surface tensions [11], [12], we see that Gmax∼0.01–0.4 μm, so that at the highest experimental loads the ratio G/R0 is typically of order 10−4 to 2×10−3, and is therefore a useful expansion parameter.
Under the fixed contact line assumption, where the radial extent r1 of the drop on the stage remains constant, the perturbation to R0, derived in [1], iswhich then leads [1] to a contact angle dependence in (3) ofWe show here, however that the correct perturbation to R0 isresulting in the contact angle dependence,Note that in addition to the probe–drop calculation, a calculation was also carried out in [1] for a spherical probe interacting with a planar substrate. The planar results are unaffected by the erroneous δR calculation. In a follow-up paper, Dagastine and White [15] extend the treatment in [1] to include a general colloidal force, with both attractive and repulsive components. The expression (7) developed here is employed in [15] for P(θ).
The constant θ (free contact line) condition was treated in [2], with the perturbation resultwhich is positive for G>0 (i.e. when the total force is repulsive), implying via (2) that the effect of a repulsive probe on the drop is to reduce the liquid pressure in the drop. Since this is physically implausible, (8) clearly cannot be correct. We find that the correct expression isand that the corresponding function P(θ) for constant θ is then
We show that the perturbation treatment breaks down for small θ (i.e. small drop volume), as would be expected physically. This breakdown is related to the singularity of δR and thus of the liquid pressure and surface curvature perturbations, evident in , , but absent in the relevant expressions (4) and (8) from [1] and [2] respectively.
A detailed analysis is given of the angular domain of validity of the perturbation theory results for δR.
The θ-dependent perturbation treatment presented here, in common with those of [1] and [2], is invalid near θ=π/2, due to the structure of the small G/R0 expansion. Since this region is obviously experimentally important, we conduct separate analyses for θ=π/2 for both constant θ and constant r1, leading in the latter case to a second order expansion. We show that the δR and P(θ) formulae, derived under the assumption that θ is not too close to π/2, do in fact give the correct θ=π/2 limit.
A nonperturbative matched asymptotic expansion is developed, avoiding inessential assumptions on θ. This method, while still invoking a small G/R expansion, does not assume that the change in curvature of the perturbed drop is small (i.e. of order G), in contrast with the perturbative treatment. We apply the resulting numerical method to determine the accuracy of the perturbative formulae for P(θ), confirming that the analytical formulae for P(θ) are useful for a large angular region, and, in particular, remain accurate in the region θ∼π/2.
Numerical results for F(X) are presented in some cases of interest, for both fixed r1 and fixed θ boundary conditions. We show that boundary condition effects on F(X) are greatest for obtuse drops with θ∼160°.
Section snippets
Formalism
The development in this section is based on that of [1], with some additional details provided for clarity. Throughout the paper, r1 denotes the contact line radius of the perturbed drop and θp is the perturbed contact angle.
We consider a drop or bubble of density ρ1, situated on an AFM stage, and immersed in a bathing medium of density ρ2. The unperturbed drop is taken to be a portion of a sphere of radius R0, with R0≪λ, where λ=(γg−1|ρ1−ρ2|−1)1/2 is the capillary length (g is the
Small G/R0 expansion
The solutions to (30) can be expressed in terms of the incomplete elliptic integrals [21] E and F,withand the upper sign now refers to acute θp, and the lower to obtuse θp. The parameter α is used to specify the solutions corresponding to the upper and lower sections of the drop. For the acute case (where only the upper part is physically relevant) and for the upper part of the obtuse case, α
Fixed contact line
If we assume a fixed contact line, r1 remains constant, and is related to the unperturbed contact angle θ,where the upper sign refers to acute, and the lower to obtuse unperturbed contact angles. Recasting (43) as an expansion in G/R0 and using the above relationship (45) between r1 and θ, we findThis result differs drastically from the corresponding result (Eq. (A.51)) in [1]. Setting Vp−Vu to zero yields, from (46), the
Fixed contact angle
When the contact angle θ remains fixed while the contact line is free to move, the perturbed contact radius r1 is a more complicated function of θ than in the fixed contact line case. A constant contact angle is equivalent to constant profile derivative at the contact line, i.e.where we have used (30). The upper sign is for acute contact angles and the lower for obtuse. Solving for r12, we findExpanding (49) in G/R0, we
Nonperturbative matched expansion in G/R
The nonsingular nature of both , at π/2 suggests that these formulae for P(θ) may be accurate beyond the angular range in which their derivation remains valid. To assess the true range of validity of these expressions, we pursue an approach which requires only that G/R≪1, for both fixed contact angle and fixed contact line cases. That is, we avoid expansions which introduce inessential conditions on the contact angle. This leads to a nonperturbative approach: we are still considering
The region θ∼π/2
Prior to applying the above numerical approach for θ∼π/2, we briefly digress to utilize , as the basis for extending the analytical (perturbation) treatment to special case θ=π/2, for both fixed r1 and fixed θ boundary conditions. We first expand the matching region terms in (56) in G/R,retaining in full the terms in r1 and using the estimate for the matching coordinate.
The region θ∼0
We now examine the numerical solutions of , as θ is reduced to values where the perturbation analysis is expected to break down, but large enough that the nonperturbative treatment is valid.
The nonperturbative approach requires in the small θ region (see Appendix A), which will be satisfied for , since R<R0 (i.e. the perturbed curvature is greater). The G/R expansion in the matching region will lose accuracy for small θ, as R decreases with θ. This ratio must therefore
Numerical F(X) results
Having assembled analytical expressions for P(θ)|r1 and P(θ)|θ and determined their domains of validity, we are now in a position to calculate some indicative force curves F(X), using numerical solutions of the inner equation (20) to calculate the integrals G(D0)=(2πγ)−1F(D0) and H(D0). Then from (3), X(D0) is evaluated as a function of central separation D0, yielding F(X), parametrically in D0.
As in [1], a purely repulsive colloidal interaction is used for the disjoining pressure, calculated
Conclusion
We have investigated the contact angle dependence, P(θ), of the AFM observable, X, for the solid probe–liquid drop configuration, pursuing both an analytical perturbation theory approach and a nonperturbative numerical approach, for both fixed contact line (constant r1) and fixed contact angle (constant θ) conditions at the liquid boundary with the substrate.
The analytical results given here, which we confirm by extensive numerical investigation, differ radically from those given in Refs. [1],
Acknowledgements
The author was introduced to Ref. [1] by Derek Chan, who suggested the fixed θ calculation in order to examine the effect of boundary conditions on experimental results, and provided the Mathematica code used in [1] for numerical evaluation of F(X). This code was then modified for the present work.
References (21)
- et al.
Forces between a rigid probe particle and a liquid interface. I. The repulsive case
J. Colloid Interface Sci.
(2001) - et al.
Effective spring description of a bubble or a droplet interacting with a particle
J. Colloid Interface Sci.
(2002) A technique for measuring the force between a colloidal particle in water and a bubble
J. Colloid Interface Sci.
(1994)- et al.
Measurement of critical disjoining pressure for dewetting of solid surfaces
J. Colloid Interface Sci.
(1996) - et al.
The direct measurement of the forces of interaction between a colloid particle and an oil droplet
J. Colloid Interface Sci.
(1996) - et al.
Quantitative analysis of fluid interface––atomic force microscopy
J. Colloid Interface Sci.
(2001) - et al.
Forces between a rigid probe particle and a liquid interface. II. The general case
J. Colloid Interface Sci.
(2002) - et al.
Effective spring constant of bubbles and droplets
Langmuir
(2001) - et al.
Measurements of hydrophobic and DLVO forces in bubble–surface interactions in aqueous solutions
Langmuir
(1994) - et al.
Surface and capillary forces affecting air bubble–particle interactions in aqueous electrolyte
Langmuir
(1996)
Cited by (28)
Extreme reversal in mechanical anisotropy in liquid-liquid interfaces reinforced with self-assembled protein nanosheets
2021, Journal of Colloid and Interface ScienceQuantifying the force between mercury and mica across an ionic liquid using white light interferometry
2019, Journal of Colloid and Interface ScienceCitation Excerpt :However, these devices are not yet in widespread use. Concurrent theoretical studies, beginning with the work of Miklavcic, Horn and colleagues [31,32] to model the equilibrium LSFA measurements and culminating in the accurate and robust numerical models developed by Chan, Manica, Klaseboer and coworkers [33–35], have been pursued to model both particle-drop, flat surface-drop, as well as more general configurations and situations, both dynamic and equilibrium [16,34–41]. Overall, through the combined theoretical and experimental effort an understanding has been established of the equilibrium colloidal interaction between a solid particle and a fluid drop or between two fluid drops across an aqueous electrolyte medium.
Interactions between micro-scale oil droplets in aqueous surfactant solution determined using optical tweezers
2018, Journal of Colloid and Interface ScienceExperimental and theoretical study of wettability alteration during low salinity water flooding-an state of the art review
2017, Colloids and Surfaces A: Physicochemical and Engineering AspectsReview and perspectives of AFM application on the study of deformable drop/bubble interactions
2015, Advances in Colloid and Interface ScienceInteraction forces between a deformable air bubble and a spherical particle of tuneable hydrophobicity and surface charge in aqueous solutions
2012, Journal of Colloid and Interface ScienceCitation Excerpt :For all the cases of R0 shown in Fig. 12, the bubble interface profile (z(r) − z(0)) changes from flattening of the original bubble surface at relatively large D0 values (∼50 nm) to the formation of a “dimple” at smaller D0 values close to the point of jump-in (D0j = 8.9 nm). The “dimple” formation as schematically shown in Fig. 1 was predicted theoretically as reported in literature for interactions of a rigid sphere with a deformable interface [44,45,24]. The effect of particle radius (a) on the normalised force (F/a) profile between a probe particle and an air bubble is shown in Fig. 13.