# Surface energy

In surface science, surface energy (also interfacial free energy or surface free energy) quantifies the disruption of intermolecular bonds that occurs when a surface is created. In solid-state physics, surfaces must be intrinsically less energetically favorable than the bulk of the material (that is, the atoms on the surface must have more energy than the atoms in the bulk), otherwise there would be a driving force for surfaces to be created, removing the bulk of the material by sublimation. The surface energy may therefore be defined as the excess energy at the surface of a material compared to the bulk, or it is the work required to build an area of a particular surface. Another way to view the surface energy is to relate it to the work required to cut a bulk sample, creating two surfaces. There is "excess energy" as a result of the now-incomplete, unrealized bonding between the two created surfaces.
Cutting a solid body into pieces disrupts its bonds and increases the surface area, and therefore increases surface energy. If the cutting is done reversibly, then conservation of energy means that the energy consumed by the cutting process will be equal to the energy inherent in the two new surfaces created. The unit surface energy of a material would therefore be half of its energy of cohesion, all other things being equal; in practice, this is true only for a surface freshly prepared in vacuum. Surfaces often change their form away from the simple "cleaved bond" model just implied above. They are found to be highly dynamic regions, which readily rearrange or react, so that energy is often reduced by such processes as passivation or adsorption.

## Assessment
### Measurement
#### Contact angle

The most common way to measure surface energy is through contact angle experiments. In this method, the contact angle of the surface is measured with several liquids, usually water and diiodomethane. Based on the contact angle results and knowing the surface tension of the liquids, the surface energy can be calculated. In practice, this analysis is done automatically by a contact angle meter.
There are several different models for calculating the surface energy based on the contact angle readings. The most commonly used method is OWRK, which requires the use of two probe liquids and gives out as a result the total surface energy as well as divides it into polar and dispersive components.
Contact angle method is the standard surface energy measurement method due to its simplicity, applicability to a wide range of surfaces and quickness. The measurement can be fully automated and is standardized.
In general, as surface energy increases, the contact angle decreases because more of the liquid is being "grabbed" by the surface. Conversely, as surface energy decreases, the contact angle increases, because the surface doesn't want to interact with the liquid.

#### Other methods

The surface energy of a liquid may be measured by stretching a liquid membrane (which increases the surface area and hence the surface energy). In that case, in order to increase the surface area of a mass of liquid by an amount, δA, a quantity of work, γ δA, is needed (where γ is the surface energy density of the liquid). However, such a method cannot be used to measure the surface energy of a solid because stretching of a solid membrane induces elastic energy in the bulk in addition to increasing the surface energy.
The surface energy of a solid is usually measured at high temperatures. At such temperatures the solid creeps and even though the surface area changes, the volume remains approximately constant. If γ is the surface energy density of a cylindrical rod of radius r and length l at high temperature and a constant uniaxial tension P, then at equilibrium, the variation of the total Helmholtz free energy vanishes and we have
    $\delta F=-P~\delta l+\gamma ~\delta A=0\quad \implies \quad \gamma =P{\frac {\delta l}{\delta A}}$
where F is the Helmholtz free energy and A is the surface area of the rod:
    $A=2\pi r^{2}+2\pi rl\quad \implies \quad \delta A=4\pi r\delta r+2\pi l\delta r+2\pi r\delta l$
Also, since the volume (V) of the rod remains constant, the variation (δV) of the volume is zero, that is,
    $V=\pi r^{2}l{\text{ is constant}}\quad \implies \quad \delta V=2\pi rl\delta r+\pi r^{2}\delta l=0\quad \implies \quad \delta r=-{\frac {r}{2l}}\delta l~.$
Therefore, the surface energy density can be expressed as
    $\gamma ={\frac {Pl}{\pi r(l-2r)}}~.$
The surface energy density of the solid can be computed by measuring P, r, and l at equilibrium.
This method is valid only if the solid is isotropic, meaning the surface energy is the same for all crystallographic orientations. While this is only strictly true for amorphous solids (glass) and liquids, isotropy is a good approximation for many other materials. In particular, if the sample is polygranular (most metals) or made by powder sintering (most ceramics) this is a good approximation.
In the case of single-crystal materials, such as natural gemstones, anisotropy in the surface energy leads to faceting. The shape of the crystal (assuming equilibrium growth conditions) is related to the surface energy by the Wulff construction. The surface energy of the facets can thus be found to within a scaling constant by measuring the relative sizes of the facets.

### Calculation
#### Deformed solid

In the deformation of solids, surface energy can be treated as the "energy required to create one unit of surface area", and is a function of the difference between the total energies of the system before and after the deformation:
    $\gamma ={\frac {1}{A}}\left(E_{1}-E_{0}\right)$.
Calculation of surface energy from first principles (for example, density functional theory) is an alternative approach to measurement. Surface energy is estimated from the following variables: width of the d-band, the number of valence d-electrons, and the coordination number of atoms at the surface and in the bulk of the solid.

#### Surface formation energy of a crystalline solid

In density functional theory, surface energy can be calculated from the following expression:
    $\gamma ={\frac {E_{\text{slab}}-NE_{\text{bulk}}}{2A}}$
where
    E$_{slab}$ is the total energy of surface slab obtained using density functional theory.
    N is the number of atoms in the surface slab.
    E$_{bulk}$ is the bulk energy per atom.
    A is the surface area.
For a slab, we have two surfaces and they are of the same type, which is reflected by the number 2 in the denominator. To guarantee this, we need to create the slab carefully to make sure that the upper and lower surfaces are of the same type.
Strength of adhesive contacts is determined by the work of adhesion which is also called relative surface energy of two contacting bodies. The relative surface energy can be determined by detaching of bodies of well defined shape made of one material from the substrate made from the second material. For example, the relative surface energy of the interface "acrylic glass – gelatin" is equal to 0.03 N/m. Experimental setup for measuring relative surface energy and its function can be seen in the video.

### Estimation from the heat of sublimation

To estimate the surface energy of a pure, uniform material, an individual region of the material can be modeled as a cube. In order to move a cube from the bulk of a material to the surface, energy is required. This energy cost is incorporated into the surface energy of the material, which is quantified by:
    $\gamma ={\frac {\left(z_{\sigma }-z_{\beta }\right){\frac {1}{2}}W_{\text{AA}}}{a_{0}}}$
where z$_σ$ and z$_β$ are coordination numbers corresponding to the surface and the bulk regions of the material, and are equal to 5 and 6, respectively; a$_0$ is the surface area of an individual molecule, and W$_{AA}$ is the pairwise intermolecular energy.
Surface area can be determined by squaring the cube root of the volume of the molecule:
    $a_{0}=V_{\text{molecule}}^{\frac {2}{3}}=\left({\frac {\bar {M}}{\rho N_{\text{A}}}}\right)^{\frac {2}{3}}$
Here, M̄ corresponds to the molar mass of the molecule, ρ corresponds to the density, and N$_A$ is the Avogadro constant.
In order to determine the pairwise intermolecular energy, all intermolecular forces in the material must be broken. This allows thorough investigation of the interactions that occur for single molecules. During sublimation of a substance, intermolecular forces between molecules are broken, resulting in a change in the material from solid to gas. For this reason, considering the enthalpy of sublimation can be useful in determining the pairwise intermolecular energy. Enthalpy of sublimation can be calculated by the following equation:
    $\Delta _{\text{sub}}H=-{\frac {1}{2}}W_{\text{AA}}N_{\text{A}}z_{b}$
Using empirically tabulated values for enthalpy of sublimation, it is possible to determine the pairwise intermolecular energy. Incorporating this value into the surface energy equation allows for the surface energy to be estimated.
The following equation can be used as a reasonable estimate for surface energy:
    $\gamma \approx {\frac {-\Delta _{\text{sub}}H\left(z_{\sigma }-z_{\beta }\right)}{a_{0}N_{\text{A}}z_{\beta }}}$

## Interfacial energy

The presence of an interface influences generally all thermodynamic parameters of a system. There are two models that are commonly used to demonstrate interfacial phenomena: the Gibbs ideal interface model and the Guggenheim model. In order to demonstrate the thermodynamics of an interfacial system using the Gibbs model, the system can be divided into three parts: two immiscible liquids with volumes V$_α$ and V$_β$ and an infinitesimally thin boundary layer known as the Gibbs dividing plane (σ) separating these two volumes.
The total volume of the system is:
    $V=V_{\alpha }+V_{\beta }$
All extensive quantities of the system can be written as a sum of three components: bulk phase α, bulk phase β, and the interface σ. Some examples include internal energy U, the number of molecules of the ith substance n$_i$, and the entropy S.
    ${\begin{aligned}U&=U_{\alpha }+U_{\beta }+U_{\sigma }\\N_{i}&=N_{i\alpha }+N_{i\beta }+N_{i\sigma }\\S&=S_{\alpha }+S_{\beta }+S_{\sigma }\end{aligned}}$
While these quantities can vary between each component, the sum within the system remains constant. At the interface, these values may deviate from those present within the bulk phases. The concentration of molecules present at the interface can be defined as:
    $N_{i\sigma }=N_{i}-c_{i\alpha }V_{\alpha }-c_{i\beta }V_{\beta }$
where c$_{iα}$ and c$_{iβ}$ represent the concentration of substance i in bulk phase α and β, respectively.
It is beneficial to define a new term interfacial excess Γ$_i$ which allows us to describe the number of molecules per unit area:
    $\Gamma _{i}={\frac {N_{i\alpha }}{A}}$

## Wetting
### Spreading parameter

Surface energy comes into play in wetting phenomena. To examine this, consider a drop of liquid on a solid substrate. If the surface energy of the substrate changes upon the addition of the drop, the substrate is said to be wetting. The spreading parameter can be used to mathematically determine this:
    $S=\gamma _{\text{s}}-\gamma _{\text{l}}-\gamma _{\text{s-l}}$
where S is the spreading parameter, γ$_s$ the surface energy of the substrate, γ$_l$ the surface energy of the liquid, and γ$_{s-l}$ the interfacial energy between the substrate and the liquid.
If S < 0, the liquid partially wets the substrate. If S > 0, the liquid completely wets the substrate.

### Contact angle

A way to experimentally determine wetting is to look at the contact angle (θ), which is the angle connecting the solid–liquid interface and the liquid–gas interface (as in the figure).
    If θ = 0°, the liquid completely wets the substrate.
    If 0° < θ < 90°, high wetting occurs.
    If 90° < θ < 180°, low wetting occurs.
    If θ = 180°, the liquid does not wet the substrate at all.
The Young equation relates the contact angle to interfacial energy:
    $\gamma _{\text{s-g}}=\gamma _{\text{s-l}}+\gamma _{\text{l-g}}\cos \theta$
where γ$_{s-g}$ is the interfacial energy between the solid and gas phases, γ$_{s-l}$ the interfacial energy between the substrate and the liquid, γ$_{l-g}$ is the interfacial energy between the liquid and gas phases, and θ is the contact angle between the solid–liquid and the liquid–gas interface.

### Wetting of high- and low-energy substrates

The energy of the bulk component of a solid substrate is determined by the types of interactions that hold the substrate together. High-energy substrates are held together by bonds, while low-energy substrates are held together by forces. Covalent, ionic, and metallic bonds are much stronger than forces such as van der Waals and hydrogen bonding. High-energy substrates are more easily wetted than low-energy substrates. In addition, more complete wetting will occur if the substrate has a much higher surface energy than the liquid.

## Modification techniques

The most commonly used surface modification protocols are plasma activation, wet chemical treatment, including grafting, and thin-film coating. Surface energy mimicking is a technique that enables merging the device manufacturing and surface modifications, including patterning, into a single processing step using a single device material.
Many techniques can be used to enhance wetting. Surface treatments, such as corona treatment, plasma treatment and acid etching, can be used to increase the surface energy of the substrate. Additives can also be added to the liquid to decrease its surface tension. This technique is employed often in paint formulations to ensure that they will be evenly spread on a surface.

## The Kelvin equation

As a result of the surface tension inherent to liquids, curved surfaces are formed in order to minimize the area. This phenomenon arises from the energetic cost of forming a surface. As such the Gibbs free energy of the system is minimized when the surface is curved.
The Kelvin equation is based on thermodynamic principles and is used to describe changes in vapor pressure caused by liquids with curved surfaces. The cause for this change in vapor pressure is the Laplace pressure. The vapor pressure of a drop is higher than that of a planar surface because the increased Laplace pressure causes the molecules to evaporate more easily. Conversely, in liquids surrounding a bubble, the pressure with respect to the inner part of the bubble is reduced, thus making it more difficult for molecules to evaporate. The Kelvin equation can be stated as:
    $RT\ln {\frac {P_{0}^{K}}{P_{0}}}=\gamma V_{m}\left({\frac {1}{R_{1}}}+{\frac {1}{R_{2}}}\right)$
where P$^K$ $_0$ is the vapor pressure of the curved surface, P$_0$ is the vapor pressure of the flat surface, γ is the surface tension, V$_m$ is the molar volume of the liquid, R is the universal gas constant, T is temperature (in kelvin), and R$_1$ and R$_2$ are the principal radii of curvature of the surface.

## Surface modified pigments for coatings

Pigments offer great potential in modifying the application properties of a coating. Due to their fine particle size and inherently high surface energy, they often require a surface treatment in order to enhance their ease of dispersion in a liquid medium. A wide variety of surface treatments have been previously used, including the adsorption on the surface of a molecule in the presence of polar groups, monolayers of polymers, and layers of inorganic oxides on the surface of organic pigments.
New surfaces are constantly being created as larger pigment particles get broken down into smaller subparticles. These newly-formed surfaces consequently contribute to larger surface energies, whereby the resulting particles often become cemented together into aggregates. Because particles dispersed in liquid media are in constant thermal or Brownian motion, they exhibit a strong affinity for other pigment particles nearby as they move through the medium and collide. This natural attraction is largely attributed to the powerful short-range van der Waals forces, as an effect of their surface energies.
The chief purpose of pigment dispersion is to break down aggregates and form stable dispersions of optimally sized pigment particles. This process generally involves three distinct stages: wetting, deaggregation, and stabilization. A surface that is easy to wet is desirable when formulating a coating that requires good adhesion and appearance. This also minimizes the risks of surface tension related defects, such as crawling, cratering, and orange peel. This is an essential requirement for pigment dispersions; for wetting to be effective, the surface tension of the pigment's vehicle must be lower than the surface free energy of the pigment. This allows the vehicle to penetrate into the interstices of the pigment aggregates, thus ensuring complete wetting. Finally, the particles are subjected to a repulsive force in order to keep them separated from one another and lowers the likelihood of flocculation.
Dispersions may become stable through two different phenomena: charge repulsion and steric or entropic repulsion. In charge repulsion, particles that possess the same like electrostatic charges repel each other. Alternatively, steric or entropic repulsion is a phenomenon used to describe the repelling effect when adsorbed layers of material (such as polymer molecules swollen with solvent) are present on the surface of the pigment particles in dispersion. Only certain portions (anchors) of the polymer molecules are adsorbed, with their corresponding loops and tails extending out into the solution. As the particles approach each other their adsorbed layers become crowded; this provides an effective steric barrier that prevents flocculation. This crowding effect is accompanied by a decrease in entropy, whereby the number of conformations possible for the polymer molecules is reduced in the adsorbed layer. As a result, energy is increased and often gives rise to repulsive forces that aid in keeping the particles separated from each other.

## Surface energies of common materials

| Material | Orientation | Surface energy (mJ/m$^2$) |
| Polytetrafluoroethylene (PTFE) | | 19 |
| Glass | | 83.4 |
| Gypsum | | 370 |
| Copper | | 1650 |
| Magnesium oxide | (100) plane | 1200 |
| Calcium fluoride | (111) plane | 450 |
| Lithium fluoride | (100) plane | 340 |
| Calcium carbonate | (1010) plane | 230 |
| Sodium chloride | (100) plane | 300 |
| Sodium chloride | (110) plane | 400 |
| Potassium chloride | (100) plane | 110 |
| Barium fluoride | (111) plane | 280 |
| Silicon | (111) plane | 1240 |