Modeling contact angle hysteresis on chemically patterned and superhydrophobic surfaces

We investigate contact angle hysteresis on chemically patterned and superhydrophobic surfaces, as the drop volume is quasistatically increased and decreased. We consider both two (cylindrical drops) and three (spherical drops) dimensions using analytical and numerical approaches to minimize the free energy of the drop. In two dimensions, we find, in agreement with other authors, a slip, jump, stick motion of the contact line. In three dimensions, this behavior persists, but the position and magnitude of the contact line jumps are sensitive to the details of the surface patterning. In two dimensions, we identify analytically the advancing and receding contact angles on the different surfaces, and we use numerical insights to argue that these provide bounds for the three-dimensional cases. We present explicit simulations to show that a simple average over the disorder is not sufficient to predict the details of the contact angle hysteresis and to support an explanation for the low contact angle hysteresis of suspended drops on superhydrophobic surfaces.     

Wetting on axially-patterned heterogeneous surfaces

Contact angle variability, leading to errors in interpretation, arises from various sources. Contact angle hysteresis (history-dependent wetting) and contact angle multiplicity (corrugation of three-phase contact line) are irrespectively the most frequent causes of this uncertainty. Secondary effects also derived from the distribution of chemical defects on solid surfaces, and so due to the existence of boundaries, are the known "stick/jump-slip" phenomena. Currently, the underlying mechanisms in contact angle hysteresis and their connection to "stick/jump-slip" effects and the prediction of thermodynamic contact angle are not fully understood. In this study, axial models of smooth heterogeneous surface were chosen in order to mitigate contact angle multiplicity. For each axial pattern, advancing, receding and equilibrium contact angles were predicted from the local minima location of the system free energy. A heuristic model, based on the local Young equation for spherical drops on patch-wise axial patterns, was fruitfully tested from the results of free-energy minimization. Despite the very simplistic surface model chosen in this study, it allowed clarifying concepts usually misleading in wetting phenomena.     

Simulation analysis of contact angles and retention forces of liquid drops on inclined surfaces

A simulation study of liquid drops on inclined surfaces is performed in order to understand the evolution of drop shapes, contact angles, and retention forces with the tilt angle. The simulations are made by means of a method recently developed for dealing with contact angle hysteresis in the public-domain Surface Evolver software. The results of our simulations are highly dependent on the initial contact angle of the drop. For a drop with an initial contact angle equal to the advancing angle, we obtain results similar to those of experiments in which a drop is placed on a horizontal surface that is slowly tilted. For drops with an initial contact angle equal to the mean between the advancing and the receding contact angles, we recover previous results of finite element studies of drops on inclined surfaces. Comparison with experimental results for molten Sn-Ag-Cu on a tilted Cu substrate shows excellent agreement.     

Modeling contact angle hysteresis of a liquid droplet sitting on a cosine wave-like pattern surface

A liquid droplet sitting on a hydrophobic surface with a cosine wave-like square-array pattern in the Wenzel state is simulated by using the Surface Evolver to determine the contact angle. For a fixed drop volume, multiple metastable states are obtained at two different surface roughnesses. Unusual and non-circular shape of the three-phase contact line of a liquid droplet sitting on the model surface is observed due to corrugation and distortion of the contact line by structure of the roughness. The contact angle varies along the contact line for each metastable state. The maximum and minimum contact angles among the multiple metastable states at a fixed viewing angle correspond to the advancing and the receding contact angles, respectively. It is interesting to observe that the advancing/receding contact angles (and contact angle hysteresis) are a function of viewing angle. In addition, the receding (or advancing) contact angles at different viewing angles are determined at different metastable states. The contact angle of minimum energy among the multiple metastable states is defined as the most stable (equilibrium) contact angle. The Wenzel model is not able to describe the contact angle along the three-phase contact line. The contact angle hysteresis at different drop volumes is determined. The number of the metastable states increases with increasing drop volume. Drop volume effect on the contact angles is also discussed.     

Assessing the accuracy of contact angle measurements for sessile drops on liquid-repellent surfaces

Gravity-induced sagging can amplify variations in goniometric measurements of the contact angles of sessile drops on super-liquid-repellent surfaces. The very large value of the effective contact angle leads to increased optical noise in the drop profile near the solid-liquid free surface and the progressive failure of simple geometric approximations. We demonstrate a systematic approach to determining the effective contact angle of drops on super-repellent surfaces. We use a perturbation solution of the Bashforth-Adams equation to estimate the contact angles of sessile drops of water, ethylene glycol, and diiodomethane on an omniphobic surface using direct measurements of the maximum drop width and height. The results and analysis can be represented in terms of a dimensionless Bond number that depends on the maximum drop width and the capillary length of the liquid to quantify the extent of gravity-induced sagging. Finally, we illustrate the inherent sensitivity of goniometric contact angle measurement techniques to drop dimensions as the apparent contact angle approaches 180°.     

Constitutive modeling of contact angle hysteresis

We introduce a phase field model of wetting of surfaces by sessile drops. The theory uses a two-dimensional non-conserved phase field variable to parametrize the Gibbs free energy of the three-dimensional system. Contact line tension and contact angle hysteresis arise from the gradient term in the free energy and the kinetic coefficient respectively. A significant advantage of this approach is in the constitutive specification of hysteresis. The advancing and receding angles of a surface, the liquid-vapor interfacial energy and three-phase line tension are the only required constitutive inputs to the model. We first simulate hysteresis on a smooth chemically homogeneous surface using this theory. Next we show that it is possible to study heterogeneous surfaces whose component surfaces are themselves hysteretic. We use this theory to examine the wetting of a surface containing a circular heterogeneous island. The contact angle for this case is found to be determined solely by the material properties at the contact line in accord with recent experimental data.     

Contact angles on spherical surfaces

In this paper, we explore the influence of curved surfaces on contact angles. Small liquid drops were deposited at the apex of spheres. Liquid was added to advance the contact line (or withdrawn to cause recession). As drop volume increased, the contact line advanced outward and downward. With the addition of each increment of liquid, the contact line encountered a steeper slope and showed progressively larger apparent advancing contact angles. Observed apparent contact angles could be explained in terms of intrinsic contact angles and surface orientation. We found that if curvature and geometry were correctly accounted for, the classic Gibbs relation held. The experimental approach and analysis used here for estimating intrinsic wettability from curved surfaces could easily be integrated into automated contact angle measurement systems.     

Anisotropic wetting characteristics on submicrometer-scale periodic grooved surface

Submicrometer-scale periodic structures consisting of parallel grooves were prepared on azobenzene-containing multiarm star polymer films by laser interference. The wetting characteristics on the patterned surfaces were studied by contact angle measurements. Macroscopic distortion of water drops was found on such small-scale surface structures, and the contact angles measured from the direction parallel to the grooves were larger than those measured from the perpendicular direction. A thermodynamic model was developed to calculate the change in the surface free energy as a function of the instantaneous contact angle when the three-phase contact line (TPCL) moves along the two orthogonal directions. It was found that the fluctuations, i.e., energy barriers, on the energy versus contact angle curves are crucial to the analysis of wetting anisotropy and contact angle hysteresis. The calculated advancing and receding contact angles from the energy versus contact angle curves were in good agreement with those measured experimentally. Furthermore, with the groove depth increasing, both the degree of wetting anisotropy and the contact angle hysteresis perpendicular to the grooves increased as a result of the increase in the energy barrier. The theoretical critical value of the groove depth, above which the anisotropic wetting appears, was determined to be 16 nm for the grooved surface with a wavelength of 396 nm. On the other hand, the effect of the groove wavelength on the contact angle hysteresis perpendicular to the grooves was also interpreted on the basis of the thermodynamic model. That is, with the wavelength decreasing, the contact angle hysteresis increased due to the increase in the number of energy barriers. These results may provide theoretical evidence for the design and application of anisotropic wetting surface.     

Water wetting transition parameters of perfluorinated substrates with periodically distributed flat-top microscale obstacles

Superhydrophobicity is obtained on photolithographically structured silicon surfaces consisting of flat-top pillars after a perfluorosilanization treatment. Systematic static contact angle measurements were carried out on these surfaces as a function of pillar parameters that geometrically determine the surface roughness, including pillar height, diameter, top perimeter, overall filling factor, and disposition. In line with thermodynamics models, two regimes of static contact angles are observed varying each parameter independently: the "Cassie" regime, in which the water drop sits suspended on top of the pillars (referred to as composite), corresponding to experimental contact angles greater than 140-150 degrees, and the "Wenzel" regime, in which water completely wets the asperities (referred to as wetted), corresponding to lower experimental contact angles. A transition between the Cassie and Wenzel regimes corresponds to a set of well-defined parameters. By smoothly depositing water drops on the surfaces, this transition is observed for surface parameter values far from the calculated ones for the thermodynamic transition, therefore offering evidence for the existence of metastable composite states. For all studied parameters, the position of the experimental transition correlates well with a rough estimation of the energy barrier to be overcome from a composite metastable state in order to reach the thermodynamically favored Wenzel state. This energy barrier is estimated as the surface energy variation between the Cassie state and the hypothetical composite state with complete filling of the surface asperities by water, keeping the contact angle constant.     

Three-dimensional equilibrium shapes of drops on hysteretic surfaces

In this paper, we study equilibrium three-dimensional shapes of drops on hysteretic surfaces. We develop a function coupled with the publicly available surface energy minimization code Surface Evolver to handle contact angle hysteresis. The function incorporates a model for the mobility of the triple line into Surface Evolver. The only inputs to the model are the advancing and receding contact angles of the surface. We demonstrate this model’s versatility by studying three problems in which parts of the triple line advance while other parts either recede or remain stationary. The first problem focuses on the three-dimensional shape of a static pendant drop on a vertical surface. We predict the finite drop volume when impending sliding motion is observed. In the second problem, we examine the equilibrium shapes of coalescing sessile drops on hysteretic surfaces. Finally, we study coalescing puddles in which gravity plays a leading role in determining the equilibrium puddle shape along with hysteresis.     

Directional self-cleaning superoleophobic surface

In this work, we report the creation of a grooved surface comprising 3 μm grooves (height ～4 μm) separated by 3 μm from each other on a silicon wafer by photolithography. The grooved surface was then modified chemically with a fluorosilane layer (FOTS). The surface property was studied by both static and dynamic contact angle measurements using water, hexadecane, and a polyethylene wax ink as the probing liquids. Results show that the grooved surface is both superhydrophobic and superoleophobic. Its observed contact angles agree well with the calculated Cassie-Baxter angles. More importantly, we are able to make a replica of the composite wax ink-air interface and study it by SEM. Microscopy results not only show that the droplet of the wax ink "sits" on air in the composite interface but also further reveal that the ink drop actually pins underneath the re-entrant structure in the side wall of the grooved structure. Contact angle measurement results indicate that wetting on the grooved surface is anisotropic. Although liquid drops are found to have lower static and advancing contact angles in the parallel direction, the drops are found to be more mobile, showing smaller hysteresis and lower sliding angles (as compared to the FOTS wafer surface and a comparable 3-μm-diameter pillar array FOTS surface). The enhanced mobility is attributable to the lowering of the resistance against an advancing liquid because 50% of the advancing area is made of a solid strip where the liquid likes to wet. This also implies that the contact line for advancing is no longer smooth but rather is ragged, having the solid strip area leading the wetting and the air strip area trailing behind. This interpretation is supported by imaging the geometry of the contact lines using molten ink drops recovered from the sliding angle experiments in both the parallel and orthogonal directions. Because the grooved surface is mechanically stronger against mechanical abrasion, the self-cleaning effect exhibited in the parallel direction suggests that groove texturing is a viable approach to create mechanically robust, self-cleaning, superoleophobic surfaces.     

A method for correcting the contact angle from the θ/2 method

The θ/2 method, a widely used technique on measuring the contact angle of a sessile drop, assumes that the drop profile is part of a sphere. However, the shape profile of a sessile drop is governed by the Young–Laplace equation and is different from a sphere, especially for drops with a large bound number (e.g. large volume or small surface tension). The spherical assumption, therefore, causes errors on evaluating the contact angles. The deviation of contact angle from the θ/2 method is evaluated from a theoretical calculation in this work. A simple means is given for correcting the measurement error. The corrected angle results from the drop volume, surface tension, liquid density and the contact angle from θ/2 method. An algorithm for finding the correct contact angle without knowing the density and surface tension is also given. At the end, two examples of pendant drops are given for the illustration.     

Wetting on the microscale: shape of a liquid drop on a microstructured surface at different length scales

Describing wetting of a liquid on a rough or structured surface is a challenge because of the wide range of involved length scales. Nano- and micrometer-sized textures cause pinning of the contact line, reflected in a hysteresis of the contact angle. To investigate contact angles at different length scales, we imaged water drops on arrays of 5 μm high poly(dimethylsiloxane) micropillars. The drops were imaged by laser scanning confocal microscopy (LSCM), which allowed us to quantitatively analyze the local and large-scale drop profile simultaneously. Deviations of the shape of drops from a sphere decay at two different length scales. Close to the pillars, the amplitude of deviations decays exponentially within 1-2 μm. The drop profile approached a sphere at a length scale 1 order of magnitude larger than the pillars' height. The height and position dependence of the contact angles can be understood from the interplay of pinning of the contact line, the principal curvatures set by the topography of the substrate, and the minimization of the air-water interfaces.     

Application of Young’s equation to two immiscible liquids resting on a solid substrate

The geometry of two immiscible liquid drops (water and mercury) resting on a solid substrate is reported. The problem under investigation is a very simple situation of a non-miscible two-phase liquid system which is properly documented. The simple experimental observation of water micro drop on top of a mercury drop on polypropylene substrate is discussed. The static situation of the drops is explored, and a new equation is derived for such situation by applying Young’s equation to the contact angles. The obtained equation simply represents the additive summation of Young equations for two independent drops, and may be applied for the case when the interfaces of two drops are close to each other and hence interact. The forces at the contact lines of each drop are treated separately as well as in correlation. This paper shows a schematic representation of two-liquid phase contact angle system. The value of this work may be appreciated from a pedagogical point of view.     

Sliding behavior of liquid droplets on tilted Langmuir-Blodgett surfaces

The sliding behavior of liquid droplets on inclined Langmuir-Blodgett surfaces was investigated. The critical sliding angle defined as the tilt angle of the surface at which the drop slides down as well as the advancing and receding contact angles was measured for five different liquids on five surfaces. In addition, the contact line geometry was analyzed at critical sliding angle. The experimental relationship between the surface tension forces resulting from contact angle hysteresis and the weight of the drop was compared to theoretical predictions. Even though the shape of the drop bases was found as skewed ellipses, a model assuming parallel-sided elongated drops is shown to describe reasonably the experimental values. This result probably indicates the main influence of the capillary forces at the rear and front edges of the drop with respect to that exerted on the lateral sides.     

Partial wetting of chemically patterned surfaces: the effect of drop size

Partial wetting of chemically heterogeneous substrates is simulated. Three-dimensional sessile drops in equilibrium with smooth surfaces supporting ordered chemical patterns are considered. Significant features are observed as a result of changing the drop volume. The number of equilibrated drops is found either to remain constant or to increase with growing drop volume. The shape of larger drops appears to approach that of a spherical cap and their three-phase contact line seems, on a larger scale, more circular in shape than that of smaller drops. In addition, as the volume is increased, the average contact angle of drops whose free energy is lowest among all equilibrium-shaped drops of the same volume appears to approach the angle predicted by Cassie. Finally, contrary to results obtained with two-dimensional drops, contact angle hysteresis observed in this system is shown to exhibit a degree of volume dependence in the advancing and receding angles. Qualitative differences in the wetting behavior associated with the two different chemical patterns considered here, as well as differences between results obtained with two-dimensional and three-dimensional drops, can possibly be attributed to variations in the level of constraint imposed on the drop by the different patterns and by the dimensionality of the system.     

The leading edge of evaporating droplets

New experiments on drops evaporating in normal atmosphere from smooth substrates in the situation of complete wetting are reported and compared with the available theoretical model. They are the continuation of previous work with alkane or water sessile drops, which is first briefly summarized. The model accounts very well for the dynamics of the drop radius, but the predictions are only qualitative for the contact angle, especially for small angles. Experiments with hanging drops allow us first to discard any influence of convection in the gas phase on the drops dynamics. Then the main part of the paper concerns new experiments with polydimethylsiloxane oligomers. These silicone oils are similar to alkanes as far as evaporation rate is concerned, but have lower surface tensions, and therefore smaller dynamic contact angles. The purity of the oils appears to be critical for the experiments, and requires a preliminary investigation. Then a systematic study of the drops dynamics is presented, as a basis for forthcoming theoretical work.     

Advancing contact lines on chemically patterned surfaces

We report an experimental investigation on advancing contact lines of large drops spreading on chemically patterned surfaces. The model substrates were prepared using microphotolithography allowing precise control of the position and the size of the wettability patterns. Experiments were performed exploring different surface geometries: from ordered to disordered fields of defects and from low to high surface densities. The shape of the contact line between two isolated defects was investigated as a function of the distance. Portions of the contact line on the defects and on the matrix were studied during spreading experiments and were related to the apparent contact angles measured from the final thickness of the drops. A modified Cassie equation based on the line fraction of defects is proposed.     

Force driven separation of drops by deterministic lateral displacement

We investigate the separation of drops in force-driven deterministic lateral displacement (f-DLD), a promising high-throughput continuous separation method in microfluidics. We perform scaled-up macroscopic experiments in which drops settle through a square array of cylindrical obstacles. These experiments demonstrate the separation capabilities-and provide insight for the design-of f-DLD for drops of multiple sizes, including drops that are larger than the gaps between cylinders and exhibit substantial deformation as they move through the array. We show that for any orientation of the driving force relative to the array of obstacles, the trajectories of the drops follow selected locking directions in the lattice. We also found that a simple collision model accurately describes the average migration angles of the drops for the entire range of sizes investigated here, and for all forcing directions. In addition, we found a difference of approximately 20° between the critical angles at which the smallest and largest drops first move across a line of obstacles (column) in the array, a promising result in terms of potential size resolution of this method. Finally, we demonstrate that a single line of cylindrical obstacles rotated with respect to the driving force is capable of performing binary separations. The critical angles obtained in such single line experiments, moreover, agree with those obtained using the full array, thus validating the assumption in which the trajectory (and average migration angle) of the drops is calculated from individual obstacle-drop collisions.     

Liquid drops on vertical and inclined surfaces; II. A method for approximating drop shapes

In this article, a new method is proposed to approximate the shapes of liquid drops on vertical and inclined surfaces. Based on observations from Part I, the profile of a drop at a given azimuthal angle is approximated by two circles sharing a common tangent at the maximum height. The drop volume is obtained by integrating all profiles over the circumference of the base. The volume is thus described as a function of the contact angles and the three-phase contact line. The new method accurately predicts the volumes of drops tested in Part I and independent measurements from the literature. Simplifying the drop shape to a spherical cap can lead to a 75% error in drop-volume prediction. The proposed method is used to study the effect of drop parameters on volume prediction. The two-circle geometry can also be used to measure contact angles from profile images.