A Guide to Surface Energy


Written by Emma Spooner, a PhD student in Fullerene-Free Photovoltaic Devices at the University of Sheffield in collaboration with Ossila Ltd.

Surface free energy is a measure of the excess energy present at the surface of a material, in comparison to at its bulk. It can be used to describe wetting and adhesion between materials, but is not often used quantitatively. This article will explain the fundamental definition of surface energy, how to calculate it directly using contact angles, and how to tune it to fit your needs. For further information on contact angles, see Ossila’s guides to the background theory of contact angle measurements, and how to take them using the Ossila Contact Angle Goniometer.




What is Surface Energy?

In the bulk form of a material, atoms are generally stable and have a balanced set of bonds/interactions. In contrast, the surface atoms will have an incomplete, unbalanced set of interactions, and therefore have unrealised bonding energy. ‘Surface energy’ is a relative measurement of the energy at the surface (which is a result of this incomplete bonding). It is positively correlated to the strength of bulk interactions and the level of surface exposure. Therefore, surface energy will be higher if the bulk interactions are stronger, or if the surface exposure is greater.

A surface will always try and minimise its energy. This can be done by adsorbing a material with a lower energy onto its surface. Through the adsorption process, the number of exposed surface atoms with high surface energy are minimised and replaced with lower energy atoms or molecules. Liquids will often have lower surface energies than solids (due to the weak forces interacting between molecules), which is why liquids usually spread out. Surface energy can be defined as the energy required per unit area to increase the size of the surface, and as such is often quoted in units of mN/m.

 

Bonding in a solid, which impacts its surface energy
Figure 1: A basic depiction of bonding in a solid, showing balanced interactions in the bulk and unbalanced interactions at the surface.

Although numerically equivalent to liquid surface tension, surface energy is mostly used to describe a solid. The tendency of a solution to spread out (or ‘wet’) on a solid surface depends on several factors. Here it is assumed that the liquid and solid are capable of similar polarity interactions.

Generally, a surface with a ‘low’ surface energy will cause poor wetting, and therefore result in a high contact angle. This is because the surface is not capable of forming strong bonds, so there is little energetic reward for the liquid to break bulk bonding in favour of interacting with the surface. Typical surfaces with low surface energy include hydrocarbons, as these are held together with weak van der Waals forces.

 

A depiction of how contact angle changes with surface energy and surface tension
Figure 2: The basic relationship between contact angle and surface energy. Surface tension of the liquid will also play a role, as will the type of interactions between the liquid and solid.

The opposite is true for a ‘high’ surface energy, which will generally cause good wetting and a low contact angle. Examples of surfaces with high surface energy include glasses, ceramics, metals, and materials held together by stronger bonds (e.g. ionic, covalent and metallic).

 Ossila Contact Angle Goniometer



How can Surface Energy be calculated?

Surface energy is usually measured indirectly with contact angle measurements, which can be a useful qualitative measure. Quantitative data can also be obtained by using a variety of different models. Most of these are based on Young’s equation,1 a form of which is seen in Equation 1:

Young's equation describing surface energy
Equation 1: Young’s equation relating surface energy of a solid, surface tension of a liquid and the interfacial tension between the liquid and solid

Where σs is the surface energy of the solid, σsl is the interfacial tension between liquid and solid, σl is the surface tension of the liquid, and θ is the contact angle of the liquid on the solid. 

Dupré’s equation can be used to relate interfacial tension and the interactions between the solid and liquid. Interactions are generally described by the work of adhesion, Wsl, which represents the work which must be done to separate the two phases, or the energy released in wetting.  This is seen Equation 2:

Dupre's equation relating adhesion and interfacial tension
Equation 2: Dupré’s equation relating interfacial tension between a solid and liquid and the interactions between the phases.

This relationship is similar to (but not exactly the same as) that between the interfacial tension and spreading coefficient, S:

Equation describing spreading coefficient
Equation 3: Equation relating the interfacial tension and spreading coefficient between a liquid and solid

The spreading coefficient measures the tendency for a liquid to wet on a solid, where wetting occurs for positive values. 

Equation 2 and Equation 4 can be combined into the Young-Dupré equation, which is the basic form used by many of the surface energy models:

Young-Dupre equation relating adhesion and contact angle
Equation 4: The Young-Dupré equation, relating the interactions between the liquid and solid to surface tension and contact angle.

Different models will include different interactions in their calculation of surface energy, meaning values, even for the same sample, are not always directly comparable between them.



Zisman Model

One of the most basic and widely used methods of calculating surface energy is the Zisman model, published in 1964.2 This model assumes that the surface energy of a solid is equal to the maximum surface tension liquid that will give a contact angle of 0°. This is known as the critical surface tension.

Contact angle on the surface (in the form cosθ) for a series of liquids is plotted against surface tension. This is then extrapolated to find the surface tension for  cosθ = 1, where θ = 0°, and surface tension is equal to surface energy of the solid.

A Zisman plot used to calculate surface energy from contact angle
Figure 3: A basic Zisman plot, where cosθ is plotted against surface tension. The maximum surface tension that achieves θ=0° is equal to surface energy. There may be liquids with lower surface tension that also achieve complete wetting (as seen above).

As the Zisman model ignores the impact of polar interactions, it is only accurate for non-polar surfaces (such as polyethylene). For polar surfaces (including those that that contain heteroatoms), a model that includes polar interactions must be used.



Fowkes & Extended Fowkes Model

The most common method of calculating surface energy for polar surfaces is the Fowkes model, published in 1964.3 Fowkes' theory assesses the interactions between a liquid and solid in terms of ‘dispersive’ (van der Waals) and ‘polar’ interactions. These interactions will sum to form overall energy:

Fowkes describes surface energy as a sum of dispersive and polar components
Equation 5: Surface energy for Fowkes theory is a sum of dispersive and polar components. The same is true for surface tension of the liquid (not shown).

Where σsD and σsP are the dispersive and polar components of the surface energy of the solid respectively.

Initial Fowkes theory is shown in Equation 6:

Fowkes theory in terms of work of adhesion
Equation 6: Fowkes theory in terms of work of adhesion.

Where σlD and σlP are the dispersive and polar components of the surface tension of the liquid respectively. This theory uses the geometric mean for each type of interaction, in contrast to the similar Wu model,4 which uses the harmonic mean.

By combining Equation 4 and Equation 6, we can obtain Fowkes' primary equation, seen in Equation 7:

Fowkes equation relating interactions and young's equation
Equation 7: Fowkes equation relating the solid and liquid phase interactions with surface tension and contact angle.

The first step for this surface energy calculation is to measure the contact angle for a purely dispersive liquid, meaning σlP = 0 and σl = σlD, reducing Equation 7 to Equation 8:

Fowkes equation for a purely dispersive liquid
Equation 8: Reduced form of Fowkes equation using a purely dispersive liquid

Using this, σsD can be calculated directly, when the surface tension for the liquid is known. A common material used here is diiodomethane, which has effectively no polar component to its surface tension (due to molecular symmetry), meaning σl = σlD = 50.8 mN/m.5

The structure of diiodomethane
Figure 4: The chemical structure of diiodomethane, which is often used to calculate surface energy due to having only dispersive components to its interactions.

The next step is to measure the contact angle for a liquid with known dispersive and polar components. A common material used here is water, where σlP = 51.0 mN/m and σlD = 21.80 mN/m.5 By inserting this into Equation 7, along with surface tension and previously calculated σsD, the value of σsP can be calculated. The final surface energy value is then a sum of the surface energy components, as seen in Equation 5.

This model is best for low-charge surfaces that are moderately polar, such as polymers featuring heteroatoms.

Fowkes' model can also be extended6 to include a third interaction component, σH, which describes the hydrogen bonding within the phase. This calculation therefore requires three liquids with known components, and thus is not used as often as the standard Fowkes model.



Owens-Wendt-Rabel & Kaelble Model

The Owens-Wendt-Rabel & Kaelble model (OWRK) was published in 19697 and 1970.8 It is mathematically equivalent to the Fowkes model, but derived from different principles.

The OWRK model is shown in Equation 9:

The equation for the Owens-Wendt-Rabel & Kaelble model
Equation 9: Owens-Wendt-Rabel & Kaelble model, which is mathematically equivalent to the Fowke’s equation seen in Eq.7, but has been expressed in the form y = mx + c.

OWRK requires at least two liquids with known dispersive and polar interactions. If these are not known they can be determined with a reference solid (such as untreated PTFE).  This can be assumed9 to have a surface energy of ~18.0 mN/m and be incapable of polar interactions, meaning σs = σsD = 18 and σsP = 0. By inserting this into Equation 9, we can achieve Equation 10:

OWRK equation for calculating dispersive interactions from contact angle on PTFE
Equation 10: A modification of the OWRK model that can be used to calculate the dispersive interaction of liquids from their contact angle on PTFE.

As a result, a contact angle on PTFE can be used to calculate σlD for any liquid where the overall surface tension is known. σlP can then be calculated using the difference between the overall surface tension and σlD.

Once dispersive and polar interactions of the liquids are known, they can be used to calculate surface energy of a new surface, by plotting each liquid on a graph in the form of Equation 9. The slope between the points of the liquids will equate to √σsP, and the intercept will equate to √σsD. The overall surface energy can then be calculated using Equation 5

An OWRK plot, used to calcaulate surface energy from contact angle
Figure 5: An OWRK plot, where contact angle is plotted against surface tension in the form of Eq.9, and the components of surface energy are found from the intercept and gradient of a best fit line.

OWRK can be used for similar materials to Fowkes, but fits better to slightly lower energy surfaces, and requires much more experimental work.



Van Oss-Good Model

Whilst dispersive and polar interactions represent many surfaces well, other significant interactions (such as hydrogen bonds) are not considered by either Fowkes or OWRK. The van Oss-Good model was published in 199210 and considers acid-base interactions. This model combines polar and dispersive interactions into one term, alongside acid (σ+) and base (σ-) components. The acid component describes the ability of a surface to interact with a basic liquid (i.e. one that can donate electron density) through polar interactions, such as dipole-dipole bonding and hydrogen bonding. The base component describes the opposite of this- e.g. interactions with an electron-accepting, acidic liquid.

The primary equation for the Oss-Good model can be seen in Equation 11:

The van-oss-good model, which uses dispersive, acidic and basic interactions
Equation 11: The Oss-Good model, where the interactions between phases are expressed in the form of dispersive, acidic and basic interactions.

By choosing a liquid with only a dispersive component, such as cyclohexane, the equation reduces to Equation 12:

The van-oss-good equation for a purely dispersive liquid
Equation 12: The reduced form of the Oss-Good model where a purely dispersive liquid is used.

This means the contact angle can be used to calculate σsD. Repeating this with a liquid with no acid component, such as tetrahydrofuran, gives Equation 13:

The reduced form of the van-oss-good model where the liquid has no acidic component
Equation 13: The reduced form of the Oss-Good model where a liquid with no acidic component is used.

This can be used to calculate σs+. Repeating this with a liquid with no base component (e.g. chloroform), will give σs-, which can be used to calculate σs by summing all three components. Alternatively, once σsD and either σs+ or σs- is known, the remaining component can be calculated from Equation 11 by using a liquid with known acid and base components (e.g. water). This is useful if a liquid with no acid/base component is difficult to find.

The Oss-Good model works best for polar surfaces, such as organometallics. However, it can be difficult to establish unknown acid and base components for liquids because there is no defined set of reference solids.



Which Model Should I Use?

  • Non-polar surfaces with low surface energies: Zisman
  • Moderately polar surfaces: Fowkes (better for slightly higher energies) or OWRK
  • Polar surfaces (inorganic, organometallic and ionic): Oss-Good
  • To find polar and dispersive interactions of liquids: Reference surface & Fowkes/OWRK

There are also additional surface energy models that are not discussed here, such as:

  • The Wu model:4 Good for materials with low surface energy up to 40 mN/m
  • The Schultz model:11 Used for high energy surfaces like bare metals 
  • The Neumann equation of state:12 Approximates surface energy with only a single liquid by ignoring the impact of different interaction components


How can Surface Energy be changed?

A high surface energy is important for solution wetting, especially in processes such as spin coating. Most solids with high surface energy will not maintain a high-energy surface when exposed to atmospheric conditions. Hydrocarbon contaminants present within the air will adsorb onto the solid's surface, reducing the surface energy. The most common method of tuning surface energy is through surface treatment, which is typically designed to raise the energy by removing these contaminants or forming high surface energy functional groups. Many of these techniques only produce temporary changes to the surface energy. This is because adsorption of low surface energy molecules will occur over time slowly reducing the average surface energy.



UV Ozone Treatment

An in-depth description of UV Ozone can be seen on Ossila's UV Ozone Cleaner page. UV radiation will ionise atmospheric oxygen to form oxygen radicals, which will go on to react and form ozone. Radiation of a different wavelength will hit organic molecules, exciting them or forming free radicals, which will allow reaction with the ozone. This forms volatile molecules which can easily desorb from the surface, leaving it clean from contaminants, and with an increased surface energy.

Ossila UV Ozone Cleaner



Plasma & Corona Treatment

Plasma treatment is also commonly used to clean surfaces and increase surface energy. Here a gas or gas mixture is ionised by a high frequency voltage to form a reactive gas plasma. This is a mixture of free radicals, ions, electrons and gas molecules. The plasma will interact with the surface in several ways. Contaminants can be removed via ablation, by bombardment with electrons and ions, or through reactions with species in the plasma. This can form volatile organic molecules that will desorb (and sometimes also replace) contaminants with high surface energy functional groups, such as C=O and HO-. In some cases, plasma treatment can etch the surface and affect surface roughness. The exact impact of plasma treatment will depend on the gas used, for example oxygen plasma will react with and oxidise organic contaminants, whereas argon and other inert gases may only remove contaminants through ablation.

A depiction of oxygen plasma treatment
Figure 6: A depiction of plasma treatment using oxygen as the gas.

Corona treatment is similar to atmospheric plasma, but uses a lower plasma density. Here, the ionised gas is created by discharging a high frequency voltage across an electrode over a grounded surface. The film is passed underneath the electrode, and the surface is oxidised to remove contaminants and raise surface energy.

Beyond UV ozone and plasma treatment, other common methods of raising surface energy include flame treatment, etching and coating a high surface energy interfacial layer of a different material.



Lowering Surface Energy

In some cases (e.g. to generate a hydrophobic coating), surface energy needs to be lowered instead of increased. This is typically done with a coating of a lower surface energy material (e.g. a wax), or through specific surface microstructures that are designed to minimise solution wetting.13



References

  1. A Thermodynamic Derivation of Wenzel’s Modification of Young’s Equation for Contact Angles; Together with a Theory of Hysteresis, R. J. Good, J. Am. Chem. Soc. (74), 5041–5042 (1952); DOI: 10.1021/ja01140a014.
  2. Relation of the Equilibrium Contact Angle to Liquid and Solid Constitution in Contact Angle, Wettability, and Adhesion, W. A. Zisman (ed. R. F. Gould), American Chemical Society, Advances in Chemistry series (43), Ch1., 1-51 (1964); DOI: 10.1021/ba-1964-0043.fw001.
  3. Attractive Forces At Interfaces, F. M. Fowkes, Ind. Eng. Chem. (56), 40–52 (1964); DOI: 10.1021/ie50660a008.
  4. Calculation of Interfacial Tension in Polymer Systems, S. Wu, J. Polym. Sci. (30), 19–30 (1971); DOI: 10.1002/polc.5070340105.
  5. Surface free-energy components of liquids and low energy solids and contact angles, B. Jańczuk et al., J. Colloid Interface Sci. (127), 189–204 (1989); DOI: 10.1016/0021-9797(89)90019-2.
  6. Studies on the surface wettability and surface structure of polyester treated with low temperature plasma, C. Jie-Rong et al., J. Appl. Polym. Sci. (63), 1733–1739 (1997); DOI: 10.1002/(SICI)1097-4628(19970328)63:13<1733::AID-APP4>3.0.CO;2-H.
  7. Estimation of the surface free energy of polymers, D. K. Owens et al., J. Appl. Polym. Sci. (13), 1741–1747 (1969); DOI: 10.1002/app.1969.070130815.
  8. Dispersion-Polar Surface Tension Properties of Organic Solids, D. H. Kaelble, J. Adhes. (2), 66–81 (1970); DOI: 10.1080/0021846708544582.
  9. Wettability of Halogenated Organic Solid Surfaces, A. H. Ellison et al., J. Phys. Chem. (58), 260–265 (1954), DOI: 10.1021/j150513a020.
  10. The Modern Theory of Contact Angles and the Hydrogen Bond Components of Surface Energies in Modern Approaches to Wettability: Theory and Applications, R. J. Good et al. (ed. M. E. Schrader), Springer US, Ch.1, 1–27 (1992); DOI: 10.1007/978-1-4899-1176-6_1.
  11. Fibre Surface Energy Characterization, J. Schultz et al, J. Adhes. (12), 221–231 (1981); DOI: 10.1080/00218468108071202.
  12. An equation-of-state approach to determine surface tensions of low-energy solids from contact angles, A. W. Neumannn et al., J. Colloid Interface Sci. (49), 291–304 (1974); DOI: 10.1016/0021-9797(74)90365-8.
  13. Super-Hydrophobic Surfaces: From Natural to Artificial, L. Feng et al., Adv. Mater. (14), 1857–1860 (2002); DOI: 10.1002/adma.200290020.