Derivations for the Wenzel model and the Cassie-Baxter model.
Recently, Young’s equation, the Wenzel equation, and the Cassie-Baxter equation have been widely used with active research on superhydrophobic surfaces. However, experiments showed that the Wenzel equation and the Cassie-Baxter equation were not derived correctly. They should be reviewed on a firm physical ground. In this study, these equations are re-derived from a thermodynamic point of view by employing energy minimization and variational approach. The derivations provide a deeper understanding of these equations and the behavior of a contact angle. Also, in applying these equations, the limitations and considerations are discussed. It is expected that this study will provide a theoretical basis for the careful use of these equations on rough or chemically heterogeneous surfaces.
- Young’s equation
- Wenzel equation
- Cassie-Baxter equation
- contact angle
- energy minimization
- variational method
The easiest way to determine the wetting property is to drop a liquid drop on the surface. The drop on the surface forms a unique contact angle depending on the wetting property. By measuring the contact angle, it is easy to examine the surface wettability. Young’s equation on the ideal surface, the Wenzel equation on the surface with roughness, and the Cassie-Baxter equation on the surface with chemical heterogeneity have been widely used for the analysis of the contact angle. Although these equations were not derived correctly, they have been used without consideration of the limitations. Application of these equations to surfaces such as a surface with large contact angle hysteresis that do not meet the conditions for these equations can give errors inherently.
In this chapter, Young’s equation, the Wenzel equation, and the Cassie-Baxter equation will be re-derived by energy minimization and variational approach. From analyses of the derivations, properties of a contact angle will be reviewed. Also, the limitations and the considerations will be discussed in applying these equations to various surfaces. We expect that this study will help in the understanding of the nature of the contact angle and its application.
1.1. Young’s equation, Wenzel equation, and Cassie-Baxter equation
It is possible to quantify the wettability of a surface by simply measuring the contact angle of a drop resting on a surface. Young’s equation has been used as a basic model. Application of this equation is limited to an ideal surface that is rigid, perfectly flat, insoluble, non-reactive, and chemically homogenous. The surface is assumed to have no contact angle hysteresis. On the surface, a contact angle of liquid drop can be described by the following Young's equation:
where , , and are liquid/gas surface tension, solid/liquid interfacial energy, and solid/gas surface energy, respectively. The apparent contact angle is an equilibrium contact angle . However, since all the real surfaces are not ideal, models were developed to describe the contact angles on the real surfaces.
There are two models to describe the contact angle on a real surface, i.e. the Wenzel model and the Cassie-Baxter model. Contrary to the ideal surface, the real surface can have chemical heterogeneity and surface roughness. The Wenzel model considers the rough surface but with chemical homogeneity . The Cassie-Baxter model considers the flat surface but with chemical heterogeneity .
In the Wenzel model, the surface roughness
where is the apparent contact angle and is the equilibrium contact angle from Young’s equation on an ideal solid without roughness.
In the Cassie-Baxter model,
where is the apparent contact angle and is the equilibrium contact angles on the solid.
From the Wenzel model, it can be deduced that the surface roughness amplifies the wettability of the original surface. Hydrophilic surface becomes more hydrophilic and hydrophobic surface more hydrophobic. In the Cassie-Baxter model, the area fractions under the drop is important in that the larger the area fraction of air, the higher the contact angle. Although these two models were proposed half a century ago, these equations have been widely used recently with active research on superhydrophobic surface [3-6].
1.2. The fallacy of the Wenzel model and the Cassie-Baxter model
In the Wenzel model and the Cassie-Baxter model, the contact angles were obtained from the non-smooth or chemically heterogeneous state of the surface under the drop. However, Gao and McCarthy demonstrated the fallacy of these models experimentally . They prepared a surface with a hydrophilic spot on a hydrophobic surface, as shown in Fig. 1a. Fig. 1b shows a smooth hydrophobic surface with a superhydrophobic spot. D and d are mean diameters of the drop and the spot.
With various diameters of the drops and the spots, advancing and receding contact angles were measured. They proved that the state of internal surface inside the triple line does not affect the contact angles experimentally and the contact angles are determined only by the state of the surface at triple contact line. It means that the previous Wenzel model and Cassie-Baxter model should be revised for rigid physical meaning . Since then, an active discussion on them has been made [8-12]. Also, these models have been derived in a more rigorous way. We have summarized the derivations of these models studied to date in Table 1. All the derivations verify that a contact angle is determined at the triple line regardless of the external fields. Experiments also confirmed these findings [13-15]. Here, we will introduce the derivations by energy minimization using simple mathematics or calculus of variations.
|Homogenization approach||N/A||At triple line||Xu and Wang, 2010 |
|Fundamental calculus||N/A||At triple line||Seo et al
|Fundamental calculus||N/A||At triple line||Whyman et al
|Variational approach||Gravity||At triple line||Bormashenko, 2009 |
|Variational approach||Electric field||At triple line||Bormashenko, 2012 |
2. Derivation with simple mathematics
For the derivation of Young’s equation in a rigorous way, the following assumptions will be used. First, the surface is ideal and it has no contact angle hysteresis. Thus, the contact line can freely move around. Second, the drop is in zero gravity and the shape of the drop is always a section of sphere, i.e., spherical cap.
As shown in Fig. 2, when the shape of the drop is deformed by spreading or contracting, the solid/liquid interfacial area varies with a contact angle that is a one-to-one function of the interfacial area. By the free movement of the contact line on an ideal surface, the drop can change freely its shape in order to satisfy the minimum energy state of the system. When the drop is at the equilibrium state, there will be no residual force at the contact line. At this point, the contact line and the shape of the drop will be fixed.
2.1. Derivation of Young’s equation
With a thermodynamic approach, Young’s equation can be derived with simple mathematics. Fig. 3 shows a drop on an ideal surface. The volume of the spherical cap is or . The surface area of the cap is or . The total energy of the system can be written as
The variation of the energy is written as
The variation of the energy is equal to zero at the equilibrium state (
means that there is no energy variation by an infinitesimal change of the shape of the drop. At this point, the shape of the drop satisfies the minimum energy state. Here, in order to induce , the constant-volume condition of the drop is used (
The volume of the drop is given by
can be obtained from the condition of the constant volume.
Rearranging the above equation gives rise to the Young's equation.
2.2. Derivation of the Wenzel equation
Fig. 4 shows a drop in the Wenzel state. The radius of the drop is
From the figure, the total energy of the system can be written as
Rearranging above equation, it will be written as
It should be noted that the terms related to
where is the equilibrium contact angle on a smooth surface and is the apparent contact angle.
In the revised Wenzel equation, the definition on the surface roughness factor is different with the previous one. The roughness factor
2.3. Derivation of the Cassie-Baxter equation
Fig. 5 shows a drop on a composite substrate that consists of two kinds of ideal surfaces. The area fraction of the red region is f1 and the yellow region is f2.
From the figure, the total energy of the system can be written as
Assuming f2 to be a fraction for contacting with air and , the revised Cassie-Baxter equation can be derived from Eq. (17). There is an important difference in the previous Cassie-Baxter equation and the revised one. The definition on the surface fraction is different. While it came from the total area under the drop in the previous model, the surface fraction in Eq. (17) should be obtained from the local region near the contact line. From the derivation, it can be also deduced that the property or the state of internal surface inside the triple line does not affect the apparent contact angle.
3. Derivation with calculus of variations
In the previous chapter, an external field, such as gravity, was not considered for simple derivation and it was possible to deal with the shape of the drop as a part of sphere. In this chapter, variational approach is employed and the shape of the drop can be distorted by the external field. Bormashenko used this approach for the first time to derive and develop the contact angle models . To understand the variational approach, fundamental formulas in calculus of variations will be introduced briefly .
3.1. Calculus of variations
The basic concept of variational method is searching a function that has an extreme value (maximum or minimum) of a physical quantity, such as energy, length, area, time, and so on. In mathematical expression, objective function
Thus, the goal is to find function ‘
With the variational method, we can solve many problems involving the determination of maxima or minima of functionals, such as the shortest smooth curve joining two distinct points in the plane , the shape of solid of revolution moving in a flow of gas with least resistance , the plane curve down that a particle will slide without friction from one point to the other point in the shortest time , the curve passing through two given points to have minimum surface area by the rotation of the curve , and the shape of the flexible cable of given length suspended between two poles .
All of the above examples involve funtionals that can be written in the form,
There is a fundamental formula for solving the simple variational problems. This is the so-called Euler equation,
When a curve passes through two fixed end points,
However, as shown in Fig. 6, when both end points of the curve are always placed on
the transversality condition is as follows:
It is possible to obtain the curve that must simultaneously satisfy certain constraints as a subsidiary condition. When the form of a functional is the same as , and the curve ‘
where is the Lagrange multiplier to be deduced from the subsidiary condition.
Therefore, the curve ‘
These fundamental equations will be used to derive the contact angle models. How a drop takes its contact angle can be understood more clearly from the variational approach.
3.2. Derivation of Young’s equation
Fig. 7 shows a drop on an ideal surface in a three-dimensional (3D) system. The symmetrical 3D drop sitting on the surface subject to energy density , due to an external field such as gravity, is considered. The free energy for the drop, , is given by the integral.
The linear density of the additional energy with the dimension of (J m-1) is defined by the expression, , where is the volume energy density of the drop in the external field (J m-3),
The constant volume
and is the Lagrange multiplier to be deduced from Eq. (26). Transversality condition is the one that should be satisfied at the end points in order to minimize the total energy of the system. The transversality condition at the end point
where denotes . Considering and , substitution of Eq. (28) into the transversality condition gives
Rearranging Eq. (30) gives
Rearrangement of Eq. (32) gives Young’s equation:
where the apparent contact angle, , is the equilibrium contact angle .
Young’s equation was derived from the transversality condition. It means that the equilibrium contact angle (Young’s equation) must be satisfied at the contact line in order to minimize the total energy of the system. The variational approach assures that the contact angle is determined at the contact line. It should be noted that the external field, such as gravity, cannot affect the equilibrium contact angle, although it distorts the shape of the drop.
3.3. Derivation of the Wenzel equation
Fig. 8 shows a symmetrical 3D drop in the Wenzel state. The drop is placed on a rough surface with full contact with the solid surface (no air gap).
The above equation can be rearranged as
where the last two terms are constants, giving
Here, since the contact angle is independent of the absolute value of the total energy, the constant energy term has no effect on the contact angle . Thus, the free energy can be redefined as follows:
Now, considering the constant volume of a drop as a subsidiary condition and the transversality condition, the revised Wenzel equation is derived:
where is the equilibrium contact angle on a smooth solid surface and is an apparent contact angle.
3.4. Derivation of the Cassie-Baxter equation
Fig. 9 shows a symmetrical 3D drop in the Cassie-Baxter state. The drop is placed on a composite surface consisting of two surfaces. The radius of the drop is
Now, considering the constant volume of a drop as a subsidiary condition and the transversality condition, the revised Cassie-Baxter equation is derived.
where and are the equilibrium contact angles on the surface 1 and the surface 2. is an apparent contact angle.
Likewise, a general Cassie-Baxter equation for different composite surfaces can be derived . In the revised equation (Eq. (41)), the area fractions are related to only the local region at the contact line.
3.5. Contact angles under other conditions
By employing the variational approach, it is easy to understand the behavior of the contact angles under various conditions. As mentioned above, after obtaining the total energy of the system under each condition, the contact angles are obtained from the transversality condition. In this way, the contact angles of the drops on a gradient surface , a rotating surface , a curved surface , and a surface with an electric field  were studied. To conclude, it was demonstrated that if the factors or conditions to be considered do not affect the surface energy and the surface topography near the contact line, they cannot affect the contact angles.
The contact angle of a drop was considered from a thermodynamic point of view. The contact angle models (Young’s equation, Wenzel equation, and Cassie-Baxter equation) were re-derived by the energy minimization and the variational approach. It was clearly demonstrated from the derivations that the contact angle of a drop is a necessary condition that must be satisfied at the contact line in order to minimize the total energy of the system. In other words, a drop takes an optimal contact angle to have the lowest energy of the system. When the optimal contact angle is not satisfied, the total energy of the system is not at a local minimum. Thus, residual force exists at the contact line and changes the shape of the drop until it disappears, as shown in Fig. 2.
Two important points can be deduced from the derivations [19, 27]. Firstly, the contact angle is determined by the infinitesimal region at the contact line. The internal surface inside the contact line does not affect the contact angle. Thus, the roughness factor in the Wenzel equation or the area fraction in the Cassie-Baxter equation should be defined in the contact line region. Secondly, the contact angle is independent of the external factors that do not affect the surface energy. Fig. 10 shows the behavior of the contact angle on an ideal surface under various conditions. The contact angle is not affected by pressure, drop size, gravity, curvature of substrate, rotation of the substrate, and existence of a needle or defects.
During the derivations, it has been assumed that the contact line moves freely on the surfaces and there is no contact angle hysteresis. However, all of real surfaces are not ideal and have contact angle hysteresis . The contact line cannot move freely on them. The surfaces have a range of static contact angles between two extreme values of an advancing angle and a receding angle . So, the real surfaces are hard to describe with a single equation, while Young’s equation, the Wenzel equation, or the Cassie-Baxter equation yield a single contact angle. Especially, when a drop on a rough surface takes the Wenzel state, the contact angle hysteresis is very large violating the assumption of the free movement of the contact line [30, 31] and a specific static contact angle is hardly meaningful to describe the surface [32, 33]. Whereas a superhydrophobic surface has a very low contact angle hysteresis having a narrow range of static contact angles. A drop on the surface takes the Cassie-Baxter state. The contact line can easily move on it. Thus, the superhydrophobic surface can be considered as pseudo-ideal surface from the theoretical viewpoint in that the contact line can move freely. For this reason, the Cassie-Baxter equation has been widely used for the superhydrophobic surface.
In using the contact angle equations, careful attention is required. When the length scale of the pattern in the surface roughness is smaller than an order of a micrometer, additional correction factors should be considered including line tension and disjoining pressure [34-36]. However, the contact angle models do not contain these factors. In addition, they do not consider the shape and size of the pattern, which actually affect the contact angle . Therefore, these equations should be used with caution and more advanced equations should be developed in order to describe the contact angle on the real surface.
Young’s equation, the Wenzel equation, and the Cassie-Baxter equations were re-derived from a thermodynamic point of view. From the derivations, the behavior of the contact angle could be deduced. In an ideal situation, the contact angle is determined by the infinitesimal region in the vicinity of contact line, not by the internal surface inside the contact line. The contact angle is also independent of the external factors that do not affect the surface energy. Thus, it is not affected by pressure, drop size, gravity, curvature of the substrate surface, rotation of the substrate, and existence of a needle or defects. It was explained from the view point of the contact angle hysteresis why these equations are not proper to describe the real common surfaces although the Cassie-Baxter equation has been widely used for a superhydrophobic surface. Also, the limitations of the equations were discussed. It is expected that this study will provide a deeper understanding of the validity of the contact angle models and the nature of the contact angle.
Wenzel, R.N. (1936). Resistance of Solid Surfaces to Wetting by Water. Industrial & Engineering Chemistry,vol. 28, pp. 988-994, ISSN 0888-5885.
Cassie, A.B.D. & Baxter, S. (1944). Wettability of Porous Surfaces. Transactions of the Faraday Society,vol. 40, pp. 546-551, ISSN 0014-7672.
Xu, W., Shi X. & Lu S. (2011). Controlled Growth of Superhydrophobic Films without any Low-Surface-Energy Modification by Chemical Displacement on Zinc Substrates. Materials Chemistry and Physics, vol. 129, pp. 1042-1046, ISSN 0254 0584.
Tuteja, A., Choi, W., Ma, M., Mabry, J.M., Mazzella, S.A., Rutledge, G.C., McKinley, G.H. & Cohen, R.E. (2007). Designing Superoleophobic Surfaces. Science, vol. 318, pp. 1618-1622, ISSN 1095-9203.
Noh, J., Lee, J.-H., Na, S., Lim, H. & Jung, D.-H. (2010). Fabrication of Hierarchically Micro- and Nano-structured Mold Surfaces Using Laser Ablation for Mass Production of Superhydrophobic Surfaces. Japanese Journal of Applied Physics, vol. 49, 106502, ISSN 1347-4065.
Wang, C., Yao, T., Wu, J., Ma, C., Fan, Z., Wang, Z., Cheng, Y., Lin, Q. & Yang, B. (2009). Facile Approach in Fabricating Superhydrophobic and Superoleophilic Surface for Water and Oil Mixture Separation, ACS Applied Materials & Interfaces, vol. 1, pp. 2613-2617, ISSN 1944-8252.
Gao, L. & McCarthy, T.J. (2007). How Wenzel and Cassie were Wrong. Langmuir, vol. 23, pp. 3762-3765, ISSN 1520-5827.
McHale, G. (2007). Cassie and Wenzel: Were They Really So Wrong? Langmuir,vol. 23, pp. 8200-8205, ISSN 1520-5827.
Panchagnula, M.V. & Vedantam, S. (2007). Comment on How Wenzel and Cassie Were Wrong by Gao and McCarthy. Langmuir, vol. 23, pp. 13242, ISSN 1520-5827.
Nosonovsky, M. (2007). On the Range of Applicability of the Wenzel and Cassie Equations. Langmuir, vol. 23, pp. 9919-9920, ISSN 1520-5827.
Marmur, A. & Bittoun, E. (2009). When Wenzel and Cassie Are Right: Reconciling Local and Global Considerations. Langmuir, vol. 25, pp. 1277-1281, ISSN 1520-5827.
Erbil, H.Y. (2014). The Debate on the Dependence of Apparent Contact Angles on Drop Contact Area or Three-Phase Contact Line: A Review. Surface Science Reports, vol. 69, pp. 325-365, ISSN 0167-5729.
Choi, W., Tuteja, A., Mabry, J.M., Cohen, R.E. & McKinley, G.H. (2009). A Modified Cassie-Baxter Relationship to Explain Contact Angle Hysteresis and Anisotropy on non-Wetting Textured Surfaces. Journal of Colloid and Interface Science, vol. 339, pp. 208-216, ISSN 0021-9797.
Seo, K., Kim, M., & Kim, D.H. (2014). Candle-Based Process for Creating a Stable superhydrophobic surface. Carbon, vol. 68, pp. 583-596, ISSN 0008-6223.
ChangWei, Y., Feng, H. & PengFei, H. (2010). The Apparent Contact Angle of Water Droplet on the Micro-Structured Hydrophobic Surface. Science China Chemistry, vol. 53, pp. 912-916, ISSN 1869-1870.
Xu, X. & Wang, X. (2010). Derivation of the Wenzel and Cassie Equations from a Phase Field Model for Two Phase Flow on Rough Surface. SIAM Journal on Applied Mathematics, vol. 70, pp. 2929-2941, ISSN 0036-1399.
Seo, K., Kim, M., & Kim, D.H. (2013). Validity of the Equations for the Contact Angle on Real Surfaces. Korea-Australia Rheology Journal, vol. 25, pp.175-180, ISSN 2093-7660.
Whyman, G., Bormashenko, E. & Stein, T. (2008). The Rigorous Derivation of Young, Cassie-Baxter and Wenzel Equations and the Analysis of the Contact Angle Hysteresis Phenomenon. Chemical Physics Letters, vol. 450, pp. 355-359, ISSN 0009-2614.
Bormashenko, E. (2009). A Variational Approach to Wetting of Composite Surfaces: Is Wetting of Composite Surfaces a One-Dimensional or Two-Dimensional Phenomenon. Langmuir, vol. 25, pp. 10451-10454, ISSN 1520-5827.
Bormashenko, E. (2012). Contact Angles of Sessile Droplets Deposited on Rough and Flat Surfaces in the Presence of External Fields. Mathematical Modelling of Natural Phenomena, vol. 7, pp. 1-5, ISSN 0973-5348.
Gelfand, I.N. & Fomin, S.V. (2000). Calculus of Variations, Dover, ISBN-13: 978-0-486-41448-5; ISBN-10: 0-486-41448-5, New York.
Eggers, A.J., Resnikoff, M.M., & Dennis, D.H. (1957). Bodies of Revolutions Having Minimum Drag at High Supersonic Air Speeds, NACA Report No. 1306.
Filobello-Nino, U., Vazquaz-Leal, H., Pereyra-Diaz, D., Yildirim, A., Peres-Sesma, A., Castaneda-Sheissa, R., Sanchez-Orea, J. & Hoyos-Reyes, C. (2013) A Generalization of the Bernoulli’s Method Applied to Brachistochrone-like Problems. Applied Mathematics and Computation, vol. 219, pp. 6707-6718, ISSN 0096-3003.
Mareno, A. & English, L.Q. (2009). The Stability of the Catenary Shapes for a Hanging Cable of Unspecified Length. European Journal of Physics, vol. 30, pp. 97-108, ISSN 1361-6404.
Bormashenko, E. (2009). Wetting of Flat and Rough Curved Surfaces. The Journal of Physical Chemistry C, vol. 113, pp. 17275-17277, ISSN 1932-7455.
Bormashenko E. (2013). Contact Angles of Rotating Sessile Droplets. Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 432 pp. 38-41, ISSN 0927-7757.
Seo, K., Kim, M., & Kim, D.H. (2015). Effects of Drop Size and Measuring Condition on Static Contact Angle Measurement on a Superhydrophobic Surface with Goniometric Technique. Korean Journal of Chemical Engineering, in print, ISSN 1975-7220.
Gleiche, M., Chi, L., Geding, E. & Fuchs, H., (2001). Anisotropic Contact‐Angle Hysteresis of Chemically Nanostructured Surfaces. Chem Phys Chem, vol. 2, pp. 187-191, ISSN 1439-7641.
Marmur, A. (2006). Soft Contact: Measurement and Interpretation of Contact Angles. Soft Matter, vol. 2. pp. 12-17, ISSN 1744-6848.
Morra, M., Occhiello, E., & Garbassi, F. (1989). Contact Angle Hysteresis in Oxygen Plasma Treated Poly(tetrafluoroethylene). Langmuir, vol. 5, pp.872-876, ISSN 1520-5827.
Cao, L. & Gao, D. (2010). Transparent Superhydrophobic and Highly Oleophobic Coatings. Faraday Discussions, vol. 146, pp. 57-65, ISSN 1359-6640.
Gao, L. & McCarthy, T.J. (2006) Contact Angle Hysteresis Explained. Langmuir, vol. 22, pp. 6234-6237, ISSN 1520-5827.
Gao, L. & McCarthy, T.J. (2009). Wetting 101°. Langmuir, vol. 25, pp. 14105-14115, ISSN 1520-5827.
de Gennes, P.G., Brochard-Wyart, F., & Quéré, D. (2003). Capillarity and Wetting Phenomena, Springer, ISBN 0-387-00592-7, Berlin.
Wong, T.-S. & Ho, C.-M. (2009). Dependence of Macroscopic Wetting on Nanoscopic Surface Textures. Langmuir, vol. 25, pp. 12851-12854, ISSN 1520-5827.
Marmur, A. & Krasovitski, B. (2002). Line Tension on Curved Surfaces: Liquid Drops on Solid Micro- and Nanospheres . Langmuir, vol. 18, pp. 8919-8923, ISSN 1520-5827.
Cansoy, C.F., Erbil, H.Y., Akar, O., & Akin, T. (2011). Effect of Pattern Size and Geometry on the Use of Cassie-Baxter Equation for Superhydrophobic Surfaces. Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 386, pp. 116-124, ISSN 0927-7757.