In this article, we provide a quick introduction to adsorption science isotherms. We begin with a derivation of the Gibbs adsorption equation. Then, we introduce two of the most common adsorption isotherm models: the Langmuir isotherm and the BET isotherm. Each section is accompanied by an example problem.
1. The Gibbs adsorption equation (Shaw, 1992; Pashley & Karaman, 2004; Cosgrove, 2005)
The Gibbs adsorption equation enables the extent of adsorption at a liquid surface to be estimated from surface tension data. Before deriving this expression, it is important to realize the limitations associated with this model – and, for that matter, all of the adsorption isotherms studied in this article. It is convenient to regard the interface between two phases as a mathematical plane, such as SS in Figure 1. This approach, however, is unrealistic, especially if an adsorbed film is present. Not only will such a film have a certain thickness, but also its presence may influence the nearby structure (for example, by dipole-dipole interaction, especially in an aqueous phase) and result in an interfacial region of varying composition with an appreciable thickness in terms of molecular dimensions. Nevertheless, theories based on discrete layer separation offer a good degree of accuracy relatively to experimental measurements. One of the earliest such models, proposed by the legendary American thermodynamicist, Josiah Gibbs, is briefly described below.
Figure 1. Interfaces in ideal and real systems.
If a mathematical plane is indeed taken to represent the interface between the two phases, adsorption can be described conveniently in terms of surface excess concentrations. Let us consider the interface between two phases, say between a liquid and a vapor, where a solute (i) is dissolved in the liquid phase. When the solute increases in concentration near the interface between the two phases, there must be a surface excess of solute , compared with the bulk value continued right up to the interface. We can define a surface excess concentration, in moles per unit area, as
where A is interfacial area. Bear in mind that may be positive or negative, and its magnitude depends on the location of SS, which is chosen somewhat arbitrarily as discussed above.
Let us now examine the effect of adsorption on the interfacial energy (). If a solute i is positively adsorbed with a surface density of , we would expect the surface energy to decrease on increasing the bulk concentration of this component (and vice versa). This situation is illustrated in Figure 2, where the total free energy of the system and are both increased by addition of component i, but, because this component is favorably adsorbed at the surface (only relative to the solvent) since both have a higher energy state at the surface, the work required to create a new surface (i.e., ) is reduced. Thus, although the total free energy of the system increases with the creation of new surface, this process is made easier as the chemical potential of the selectively adsorbed component increases (i.e., with concentration). This reduction in surface energy must be directly related to the change in chemical potential of the solute and to the amount adsorbed, and is therefore given by the simple relationship
or, in the case of several components,
The change in is caused by the change in bulk solute concentration. The result above is the Gibbs surface tension equation. Basically, this equation describes the fact that increasing the chemical potential of the adsorbing species reduces the energy required to produce new surface (i.e., ). This is the main physicochemical feature of surface-active materials, more commonly known as surfactants.
Figure 2. Change in surface energy caused by addition of a solute.
Gibbs’ equation can be readily extended to a solution of two components, in which case we may write
and hence the adsorption excess for one of the components is given by
Thus, in principle, we could determine the adsorption excess of one of the components from surface tension measurements, if we could vary independently of . But this appears to not be possible, because the chemical potentials are dependent on the concentration of each component. However, for dilute solutions the change in for the solvent is negligible compared to that of the solute. Hence, the change for the solvent can be ignored, leading to the straightforward result
differentiation with respect to concentration gives
so that, combining this result with the expression for , we ultimately obtain
This relation is the Gibbs adsorption isotherm in its usual formulation. Experimental measurements of over a range of concentrations allows us to plot against and thereby obtain , the adsorption density at the surface. The validity of this fundamental representation of adsorption phenomena has been corroborated by comparison with direct adsorption measurements.
It is important to bear in mind that the Gibbs equation in the form above applies to non-dissociating materials (e.g., non-ionic surfactants). For ionic surfactants of the form , we can write, assuming ideal behavior,
If no electrolyte is added, electroneutrality of the surface requires that = = . Substituting above brings to
Cosgrove (2005) recommends using the equation above with the mean ionic activity = .
For a M aqueous solution of n-butanoic acid, = –0.08 N-m²/mol at 25ºC. Using the Gibbs adsorption equation, determine the surface excess of butanoic acid and calculate the surface area available to each molecule.
Note that mol/dm³ = 0.1 mol/m³. Substituting in the Gibbs equation gives
The surface area available for each molecule is then
2. Adsorption isotherms – the Langmuir isotherm (Modified from Tien, 2019)
The Langmuir isotherm is perhaps the most widely used isotherm expression for representing physical adsorption equilibrium data. Although it may be derived in several different ways, the kinetic approach is used in most textbooks. The basic assumptions are:
1. Adsorption of adsorbates takes place at well-defined sites.
2. All adsorption sites are the same (energetically), and each site accommodates only one adsorbate molecule.
3. There is no lateral interaction between adsorbed molecules.
The adsorption rate is assumed to be
where is the adsorption rate constant, is adsorption, is the maximum adsorption (when all sites are filled), and is the adsorbate pressure. The desorption rate, in turn, is given by
where is the desorption rate constant. At equilibrium, we must have
Let / = and / = . The latter is the Langmuir adsorption constant. Making these substitutions brings to
This is the basic form of the Langmuir adsorption isotherm.
At 273 K and 10 bar, the Langmuir adsorption of a gas on a solid surface gave the fraction of surface coverage as 0.01. In this case, what is the value of the Langmuir adsorption constant?
From the Langmuir isotherm equation, we have
We were given the fraction of surface coverage = 0.01 and the pressure = 10 bar. Substituting and solving for K, it follows that
3. Adsorption isotherms – the Brunauer-Emmet-Teller (BET) isotherm (Ghosh, 2009)
The theory of chemical adsorption proposed by Langmuir is based on the formation of a monolayer. However, multilayers can form in physical adsorption. The isotherm in the case of a monolayer assumes the shape shown as Type I in Figure 3. The extent of adsorption increases with pressure and ultimately reaches a limiting value, as predicted by the Langmuir theory. However, for adsorption involving multilayer formation at low temperatures, at least four other types of adsorption isotherm can be observed; see below. (In the graphs shown, y denotes the amount of gas adsorbed by m kg of adsorbent at pressure P.)
Figure 3. Adsorption isotherms.
To explain such varied adsorption isotherms, Brunauer, Emmet and Teller proposed, in 1938, a theory of multilayer adsorption. They posited that, even after the formation of the monolayer, other layers of gas may condense on it. Langmuir’s idea of fixed adsorption sites was, however, retained. They suggested that there is equilibrium between the first layer and the adsorbent, and similar dynamic criteria were assumed to exist for the successive molecular layers. The Van der Waals forces provide the binding energy in these successive layers. The equation derived by them is
where is the volume of gas adsorbed under a pressure , is the saturated vapor pressure, is the volume of gas adsorbed when the surface is covered with a monolayer, and is a constant that can be related to the heat of adsorption of the first layer () and the enthalpy of liquefaction of the gas (),
From the BET equation, it is clear that a plot of versus should yield a straight line with the following characteristics:
Accordingly, can be determined with the simple relation
Thus, the BET equation permits us to extract from multilayer adsorption data the volume of adsorbed gas that would saturate the surface if the adsorption were limited to a monolayer. The most immediate application of this quantity is the calculation of specific surface area. The volume of gas adsorbed when the surface is covered by a monolayer, , is recorded at standard temperature and pressure, STP, at which the volume occupied by one mole of gas is m³. The number of moles of gas adsorbed is /(), where is Avogadro’s number. If the cross-section of each adsorbed gas molecule is and if m kg of adsorbent is present, then the surface area per kg of adsorbent is given by
Here, is expressed in m³ and in m². The following example illustrates use of the BET isotherm in calculations of specific surface area.
The following data were obtained for adsorption of nitrogen on 0.88 kg of nitrogen-activated alumina. The data are adapted from Tien (2019).
A) Estimate the volume of nitrogen adsorbed over 1 kg of alumina when a monolayer is formed.
B) Calculate the specific surface area of this alumina sample. The area occupied by one nitrogen molecule is m².
According to the BET model, a plot of versus should yield a straight line. The data we need are processed below.
Then, the data in the red column are plotted against the data in the blue column, as shown.
Clearly, there is a linear trend in this dataset. Fitting these data to a first-degree equation yields a line with slope
Solving both equations simultaneously gives = 0.0739 m³ and = 108.2. The volume of nitrogen adsorbed over 1 kg of alumina is then
Lastly, the specific surface area is
• COSGROVE, T. (Ed.) (2005). Colloid Science: Principles, Methods and Applications. Oxford: Blackwell Publishing.
• GHOSH, P. (2009). Colloid and Surface Science. New Delhi: PHI Learning.
• PASHLEY, R. and KARAMAN, M. (2004). Applied Colloid and Surface Chemistry. Hoboken: John Wiley and Sons.
• SHAW, D. (1992). Introduction to Colloid and Surface Chemistry. 4th edition. Oxford: Butterworth-Heinemann.
• TIEN, C. (2019). Introduction to Adsorption. Amsterdam: Elsevier.