Short news


See System

About us



Which model do you use for calculation of surface free energy?
Surface Free Energy PDF Tisk


The determination of surface free energy of a solid is of a great importance in a wide range of problems in applied science. Because of the difficulties to measure directly the surface free energy of a solid phase, indirect approaches have been used. Among the indirect methods used for estimation of the surface free energy of the solids, the contact angle measurements are believed to be the simplest. The contact angle measurement could be easily performed by establishing the tangent angle of a liquid drop with a solid surface. The contact angle of a liquid drop on a solid surface is defined by the mechanical equilibrium of the drop under the action of three interfacial tensions solid/vapor, solid/liquid and liquid/vapor. The determination of the surface free energy of solids from contact angles relies on a relation, which has been proposed by Young in 1805 and the equilibrium relation is known as Young’s equation. In the literature there is wide range of different methods for surface free energy calculation which are based on the Young’s relation [1-13]. There are various graphical or numerical methods (Zisman Theory, Fowkes Theory, The Extended Fowkes Theory, Wu Equation of State, Owens–Wendt-Raeble-Kaeble Theory, Lifshitz/van der Waals-Lewis Acid/Base Theory etc.). One of the most often used methods is the Owens-Wendt-Rable-Kaeble (OWRK) method.

Surface Tension and Surface Free Energy

In case of liquid surfaces the term “surface tension” and the term “surface free energy” are often interchanged. The term surface tension arose from historical concept that the liquid has a “skin” which can exert a force per unit length.

Generally, the surface tension γ can be defined as the reversible work, W done in creating a unit surface area.

W = γΔA                                                                                      (1)


Here ΔA is the newly created surface area. γ may be associated with units of either Jm-2 (surface free energy) or Nm-1 (surface tension). These two units are dimensionally equivalent.

The total Helmholtz free energy Ftot of the system of may be written in the following form

Ftot = NFa +AFsurf                                                                                      (2)

Here Fa is the Helmholtz free energy per atom, N is the number of atoms of the bulk material and Fsurf is the Helmholtz free energy per unit area of the surface:

dFtot = -SdT – PdV + mdN +FsurfdA  ›  dFtot(constant T,V and N) = FsurfdA               (3)

Here S is the entropy, P is the pressure and m is the chemical potential. The reversible work done in creating of unit area of surface is the Helmholtz free energy for that area at constant temperature T and volume V

.                                                         (4)


From (3) and (4) follows:

γ = Fsurf                                                          (5)

The surface tension equals the Helmoltz free energy for the surface, which is generally termed the surface free energy.

In case of solid surfaces the interatomic/intermolecular forces are stronger than in case of liquids, the solid material is not able to rearrange itself in order to minimise its energy. There is another contribution to the surface free energy, contribution of the surface stress, resulting from the work required to deform the surface. The relationship between the surface tension, surface stress and surface free energy is too complex to deal with here in detail.

Fig 1.1 Illustration of cohesion between identical substances i and adhesion between different substance i and j.

The surface free energy of solids can be described as a measure of disruption of chemical bonds that occurs when a new surface is created. Surfaces must be energetically less favourable than the bulk of the material otherwise there would be a driving force for creation of new surfaces. Cutting a piece of a solid material in half (see Fig 1.1) breaking its bonds consumes energy.  If the cutting is done reversibly, then the energy W consumed by the cutting process is equal to the energy of the two newly created surfaces 2γi. The unit surface energy would be therefore half of the cohesion energy Wc /2= γi (see Fig 1.1) of the solid material. In practice it is true only for surfaces newly prepared in ultra-high vacuum conditions. All surfaces are energetically unfavourable in that they have a positive energy of formation. Due to positive energy of newly created surface there are driving forces leading to creation of surface layers of natural oxides or contaminations. Therefore clean surfaces in ultra-high vacuum conditions have higher free energy, than oxidised or contamined surfaces.

Energy of adhesive bonds Wa is energy released at creation of adhesive bonding between two different materials i and j. The energy of adhesion can be also defined as the energy which is needed for the separation of two different materials i and j.

Wa may be determined according to Dupre formula [14]:

Wi + Wj = Wij + Wa,                                                             (6)

where Wi  is the cohesion energy of the substance i, Wj  is the cohesion energy of the substance j and Wij is the energy of the solid-solid or solid-liquid interface.

Direct measurements of surface free energy of solid surface Wj and of the interfacial energy Wij are rather difficult, so the equation in the above described form could not be used easily for determination of energy of adhesion Wa.

Young’s Equation

If a drop of a liquid l is brought into contact with a flat horizontal solid surface s in a gas atmosphere v, it can either spread over the surface or assume the shape of a spherical segment from a flat lens up to an almost complete ball. The most well known description of wetting state was derived by Young in 1805 [13]. The Young’s equation describes the equilibrium of forces between the surface tension at the three-phase boundary:

 ,                                                     (7)


where γsv is the solid-vapor interfacial free energy, γsl is the solid-liquid interfacial free energy, and γlv is the liquid-vapor interfacial tension.



Fig 1.2 Description of the wetting state by Young [13]:γsv is the solid-vapor interfacial free energy, γsl is the solid-liquid interfacial free energy, and γlv is the liquid-vapor interfacial tension.


If the contact angle is equal to 0° then the liquid will completely spread over the surface of the solid phase. If the contact angle is equal to 180°, liquid will not wet the solid surface. Generally, the solid surface is called wettable if the contact angle of the liquid drop is less than 90° and is not wettable if the contact angle is higher than 90°.  The solid surfaces may be divided into two groups, they may be hydrophilic (wettable with water) or hydrophobic (not wettable with water). The first case concerns the high-energy surfaces and the second case concerns the low-energy surfaces.

High-energy surfaces such as metal, oxide or ceramic surfaces should be completely wetted by fluids. However, the complete wetting (q =0°) may be obtained only in the case of completely clean, highly polished surfaces and pure liquids in inert gas atmospheres that are completely free of contaminants. Even a short-time contact with a gas atmosphere of less than 1% relative atmospheric humidity substantially changes the surface properties of the high-energetic materials.


From known γsv , γsl  and γlv the following thermodynamic quantities may be defined:


Work of adhesion                                                              (8)

Spreading coefficient                                                              (9)

Wetting energy                                                              (10)

Harkins and Feldman [15] defined the initial spreading coefficient So as

,                                                             (11)

where both, γso and Wao are associated with vapor-free “dry” surface, Wc is the work of cohesion and is the change in the Helmholtz free energy of spreading per unit area.

The liquid spreads spontaneously over the solid surface, when the initial spreading coefficient So > 0 and forms a drop resulting in a partial wetting, when So < 0.

The equilibrium spreading coefficient of a system may be expressed in term of the difference between the work of adhesion and cohesion according to


The difference in the initial and final spreading coefficient is the so called spreading pressure


The Young’s equation can be expressed also in terms of γso and P:


The spreading pressure may be determined from the Gibbs adsorption isotherm equation determining the surface coverage of the adsorbate G at several values of vapour pressures p

,                                                (15)

where R is the gas constant and T is the temperature in Kelvins, G is the amount of the vapour adsorbed in moles per unit area, po is the initial value of the pressure  ps is the saturated vapour pressure. The spreading pressure is negligible for contact angles higher than 10o.

Lifshitz -Van der Waals Interactions

The non-covalent and non-electrostatic molecular forces are referred to as van der Waals forces [23, 24]. These intermolecular forces originate from:

  • Randomly orienting dipole-dipole interactions (Keesom)
  • Randomly orienting dipole-induced dipole (Debye) interactions
  • Dispersion  interactions (fluctuating dipole-induced dipole interactions) described by London

Keesom’s dipole-dipole interaction

According to Keesom [25], two molecules with permanent dipole orient themselves in such a way as to attract each other. Interaction of freely rotating dipoles depends on the thermal energy kT and is referred to as Keesom‘s energy:

.                                       (16)

The Keesom’s energy between two different molecules i and j is expressed as


where V is the potential energy, mi and mj dipole moments, eo dielectric permitivity, k Boltzman constant, T absolute temperature and r is the distance between two interacting molecules.

Debye’s dipole-induced dipole interactions

According to Debye [26] a molecule with permanent dipole induces a dipole in the neighboring neutral molecule.  The corresponding potential energy is referred to as Debye energy:

,                                                  (18)



where a is the polarizability. Vii is the formula for the interaction between molecules of the same types of molecules, Vij is the formula for interaction between different molecules. The interaction energy is independent of the temperature because the induced dipole follows the motion of the permanent dipole independently on the thermal motion.

London’s dispersive interactions

Random fluctuation in a polarizable molecule leads to a temporary dipole, which induces a corresponding dipole in a nearby molecule. According to London [27, 28], the dispersive interaction between the like molecules may be described using the following formula

,                                         (20)

where h is the Planck’s constant, ν is the frequency of fluctuation and C is the London’s constant. In case of unlike molecules the dispersion interaction can be expressed as


where I  = hν is the ionisation energy.

The London’s interaction is universal and is present in atom-atom interactions.  Of the three components of intermolecular forces, only van der Waals – London interactions have significant importance between the macroscopic bodies in the condensed systems. The dispersion forces play significant role in such important processes as adhesion, polymer conformation, physical adsorption etc.

Keesom’s and Debye’s interactions are found only among the molecules which have permanent dipole moments. The Keesom and Debye theory explain the interactions between the polar molecules.

The total van der Waals interaction energy can be expressed as the sum of Keesom, Debye and London interactions in the following form:


Van der Waals forces between two parallel flat surfaces

The total attractive dispersion energy, V between two flat parallel surfaces separated by a distance r is given by [29, 30]


Here Aii is the so called Hamaker constant and qi is the number of atoms per unit volume.

Berthelot’s geometric mean combining rule

According to Berthelot’s principle [31], the interaction constant between two materials equals to the geometric mean of the interaction constant of individual materials.

  and                                           (24)

For identical materials i in environment k the combining rule for Hamaker constants can be written as


For materials i and j in an environment k, the combining rule for Hamaker constants is given by


Fowkes approach

According to Fowkes [29, 32] the Hamaker constant can be expressed using the following formulae


where  is the dispersion component of the surface free energy of the material. Fowkes found, that the value of  is 1.44x10-14 cm2 for most materials.

Lifshitz approach

Lifshitz [39] calculated the interaction energy between two macroscopic materials using quantum mechanical and electromagnetic approach and obtained for Aijk the following expression:



where ei, ej and ek are the dielectric constants of the three media and ni, nj and nk are the refractive indexes of the same.

In each of these expressions the first term on the rigth side gives the zero frequency part of the van der Waals energy, while the latter term represents the dispersion energy which includes London’s energy contribution. ne is the absorption energy in the UV region (~ 3x1015 s-1).

Interaction parameter

According to Berthelot’s rule, the free energy of adhesion is equal to the geometric mean of free energies of cohesion of the separate phases. For cases where both phases are apolar


and apolar systems is valid the following equation

.                                                       (29b)

In case of polar phases, the experimental results are often in conflict with Eq. (1.29b).  In order to correct these discrepancies the interaction parameter Φ was defined in the following way:


The free energy of adhesion is the negative work of the adhesion and the free energy of cohesion is the negative work of cohesion:


Using the Dupré‘s equation  (6) and the geometric combining rule , the following equation will be obtained


From Young’s equation it is possible to obtain the relationship between the equilibrium contact angle and the interaction parameter:



.                                             (35)

The interaction parametr Φ depends on the mutual properties of the liquid and the solid


Φ = 1 for so called “regular” interfaces, i.e. systems for which the cohesive forces of the two phases and the adhesive forces across the interface are of the same type.

Φ < 1 the dominant forces within the separate phases are unlike, e.g. London–van der Waals vs. metallic or ionic or dipolar.

Φ >1 there are specific interactions between the molecules forming the two phases.

Interaction parameter obtained in formula (31) is closely related to the molecular interaction  parameters D.  For substance i the following expressions are valid:

,                                                  (36)

,                                                  (37)

,                                                    (38)

where Dip is the polar component of molecular interactions, DiL is the dispersion (London) component of molecular interactions, Dii is the component of dipole induced dipole interactions .

.                                                               (39)

The interaction parameter may be expressed as


Substituting values of interaction constants from equations (1.16-1.21) into (39) the following formula is obtained:

.                        (41)

The equation for polar components is valid only if the dipoles are maximum of about 2-2,5 Debyes. Then is possible to write

.                                                     (42)

This equation is not valid, when one of the dipoles is too big or its rotation is limited.

If Is and Il are approximately of the similar value and asms2 and alml2 are also almost the same,

then the above described relationships can be simplified in the following way:



We can express F  in terms of the particular components of the molecular interactions

.                                            (45)

Acid-base interactions

The new term, Lewis acid/base interaction was introduced in 1986 by van Oss et. al. in order to explain the interaction between biopolymers and low–energetic solids. Acid/base interactions occure when a base (electron donor or proton acceptor) and an acid (electron acceptor or proton donor ) are given together. This can be described by the general equation:


The acid-base interactions are the following:

  1. Hydrogen bonding compounds treated under Bronsted proton donor/proton acceptor theory:
    • Substances that are both proton donors (acids) and proton acceptors (bases), for example water.
    • Substances that are much more effective as proton donors than as proton acceptors, e.g. CHCl3
    • Substances that are much more effective as proton acceptors  than as proton donors, e.g. ketones
  2. Compounds that interact as Lewis acids (electron acceptors) and bases (electron donors).
    • Substances that have both kinds of functionalities, electron acceptors and donors, termed as bipolar.
    • Substances that are much more effective as electron donors than as electron acceptors.
    • Substances that are much more effective as electron acceptors than as electron donors.


[1] Good, R. J., and van Oss, C. J., in Modern Approaches to Wettability: Theory and Applications, Eds: M. E. Schrader, and G. Loeb, Plenum Press, New York, pp. 1-27, 1992.

[2] Padday, J. F., (Ed), Wetting, Spreading, and Adhesion, Academic Press, New York, 1978.

[3] Zisman, W. A., in: Handbook of Adhesives, Skeist, I., Ed: van Nostrand, NewYork, Chapter 3, 1977.

[4] Fowkes, F.M., Industrial and Engineering Chemistry, 56/12, 40, 1964.

[5] Fowkes, F.M., J. of Adhesion, 4, 153, 1972.

[6] Van Oss, C. J., Chaudhury, M. K., and Good, R. J., J. Separation Sci. and Technology, 22, 1, 1987.

[7] Van Oss, C. J., Good, R. J., and Chaudhury, M. K., Langmuir, 4, pp. 884-891, 1988.

[8] Van Oss, C. J., Chaudhury, M. K., and Good, R. J., J. Colloid Interface Science,

[9] Neumann, A. W. and Good, R.J., Techniques of Measuring Contact Angles, Surface and Colloid Science, Volume II, Experimental Methods, Eds: R. J. Good and R. R. Stromberg, Plenum Pres, New York, 1979.

[10] Van Oss, C. J., Giese, R. F., and Good, R. J., Langmuir, 6, pp. 1711-1713, 1990.

[11] Spelt, J. K., Li, D., and Neumann, A. W., in Modern Approaches to Wettability: Theory and Applications, Eds: M. E. Schrader, and G. Loeb, Plenum Press, New York, pp. 101-142, 1992.

[12] Van Oss, C. J., Interfacial Forces in Aqueous Media, Marcel Decker Inc., New York, 1994.

[13] Young, T., 1805. in Ref. 12.

[14] Dupré, A., 1869. in Ref. 12.

[15] Bruil, H. G., and van Aartsen, J. J., Colloid & Polymer Science, 252, pp. 32-38, 1974.

[16] Crawford, R., Koopal, L. K. and Ralston, J., Colloids and Surfaces, 27, 57-64, 1987.

[17]. Van Oss C.J., Griese, R.F. Jr., Li, Z., Murphy, K., Norris, J., Chaudhury, M. K., and Good, R. J., J. Adhesion Sci. and Technology, Vol. 6, No.4, 413-428, 1992.

[18] Wu, W., Griese, R.F. Jr., and van Oss C.J., Powder Technology, 89, 129-132, 1996.

[19] Young, G. J., Chessick, J. J., Healey, F. H., and Zettlemoyer, A. C., J. Phys. Chem., 58, 313, 1954.

[20] Fowkes, F.M., Ind. and Eng. Chem., 56/12, 40, 1967.

[21] Fowkes, F. M., and Joslin, S. T., see Ref. 75, 1983.

[22] Zettlemoyer, A. C., in Surfaces and Coatings Related to Paper and Wood, Eds: R. H. Marchessault and C. Skaar, Syracuse University Press, 1967.

[23] Van der Waals, J. D., 1873, in Ref. 12.

[24] Van der Waals, 1899, in Ref. 12.

[25] Keesom, W. H., Phys. Z., 22, 643, 1921.

[26] Debye, P., Phys. Z., 21, 178, 1920.

[27] London, F., Z. Phys. Chem., 11, 222, 1930.

[28] London, F., Trans. Faraday Soc., 33, 8, 1937.

[29] Fowkes, F.M., in: Physicochemical Aspects of Polymer Surfaces, Ed: K.L. Mittal, vol. 2, p. 583, Plenum Press, New York, 1983.

[30]. Israelachvili, J. N., Intermolecular and Surface Forces, Second Edition, Academic Press, London, 1991.

[31] Berthelot, D., Comp. Rend. Acad. Sci., 126, 1703, 1898.

[32] Hamaker, H. C., Physica, 4, 1058, 1937.

[33] Israelachvili, J. N., Intermolecular and Surface Forces: With Applications to Colloidal and Biological Systems, Academic Press, London, 1985

[34] Fowkes, F. M., and Pugh, R. J., in Polymer Adsorption and Dispersion Stability,Eds: E. D. Goddard and B. Vincent, American Chemical Society, 1984.

[35] Hough, D. B., and White, L. R., Adv. Colloid Interface Sci., 14, 3, 1980.

[36] Van Oss, C. J., and Good, R. J., Colloids and Surfaces, 8, 373, 1984.

[37] Van Oss, C. J., Chaudhury, M. K., and Good, R. J., Chem. Revs., 88, 927, 1988.

[38] Padday, J. F., and Uffindell, N. D., J. Phys. Chem., 72, 1407, 1968.

[39] Lifshitz, E. M., J. Exp. Theor. Phys., 29, 94, 1955.

[40] Rytov, S. M., Theory of Electric Fluctuations and Thermal Radiations, Moscow Academy of Science Press, Moscow, 1953.

[41] Israelachvili, J. N., Q. Rev. Biophys., 6, 341-387, 1974.

[42] Garifalco, L. A., and Good, R. J., J. Phys. Chem., 61, 904, 1957.

[43] Fowkes, F. M., J. Phys. Chem., 67, pp. 2538-2541, 1963.

[44] Fowkes, F. M., Ronay, G. S., and Schick, M. J., J. Phys. Chem., 63, 1684, 1959.

[45] Fowkes, F. M., J. Phys. Chem., 64, 726, 1960.

[46] Bolger, J. C., and Michaels, A. S., Interface Conversion for Polymer Coatings, Ed: P. Weiss, pp. 3-60, Elsevier, New York, 1968

[47] Lipatov, Y. S., and Sergeeva, L. M., Adsorption of Polymers, Halsted, New York, 1974.

[48] Sorensen, P., J. Paint Technol.,47, 31, 1975.

[49] Fowkes, F. M., and Mostafa, M. A., Ind. Eng. Chem. Prod. Res. Dev., 71, pp. 3-7, 1978.

[50] Fowkes, F. M., J. Adhes. Sci., 1, pp. 7-27, 1987.

[51] Van Oss, C. J., Good, R. J., and Chaudhury, M. K., J. of Colloid and Interface Science, 111, pp. 378, 1986.

[52] Van Oss, C. J., Chaudhury, M. K., and Good, R. J., Adv. Coll. Interface Sci., 28, pp. 35-64, 1987.

[53] Van Oss, C. J., Chaudhury, M. K., and Good, R. J., J. of Colloid and Interface Science, 128, 313, 1988.

[54] Fowkes, F. M., J. Phys. Chem., 66, 382, 1962.

[55] Morra, M. J., J. Colloids Interface Sci.,182, 312, 1996.

[56] Kwok, D. Y., Lee, Y. and Neumann, A. W., Langmuir, 14, pp. 2548-2553, 1998.

[57] Van Oss, C. J., Chaudhury, M. K., and Good, R. J., Adv. Colloid Interface Sci., 28, 35, 1987.

[58] Van Oss, C. J. and Good, R. J., J. Dispersion Sci. Tech., 1988.

[59] Fowkes, F. M., in Surfaces and Interfaces, Eds: J. J. Burke, Syracuse University Press, 1967.

[60] Laskowski, J. S., and Kitchener, J. A., J. of Colloid and Interface Science, 29, 670, 1969.

[61] Bangham, D. H., and Razouk, R. I., Trans. Fraday Soc., 33, 1459, 1937.

[62] Harkins, W. D., and Livingston, H. K., J. Chem. Phys., 10, 342, 1942.

[63] Fowkes, F. M., McCarthy, D. C., and Mostafa, M. A., J. of Colloid and Interface Science, 78, 200, 1980.

[64] Busscher, H. J., Kim, K. G. A., Van Silfout, G. A. M., and Arends, J., J. of Colloid and Interface Science, 114, 307, 1986.

[65] Bellon-Fontaine, M.-N. and Cerf, O., J. Adhesion Sci. and Technology, 4, 475, 1990.

[66] Good, R. J., Surface Colloid Sci., 11, 1, 1979.

[67] Neumann, A. W., Adv. Colloids Interface Sci., 4, 105, 1974.

[68] Israelachvili, J. N., and McGee, M. L., Langmuir, 5, 288, 1989.

[69] Drelich, J., and Miller, J. D., Langmuir, 9, 619, 1993.

[70] Cassie, A. B. D., and Baxter, S., Trans. Faraday Soc., 40, 456, 1944.

[71] Wenzel, R. N., J. Phys. Chem., 53, 1466, 1949.


From Bursikova V et al. (2004): Surface Energy Evaluation of Plasma Treated Materials by Contact Angle Measurement. Masaryk University, Brno. ISBN 80-210-3563-3.