Langmuir adsorption model

Fig 1. A schematic showing equivalent sites, occupied (blue) and unoccupied (red) clarifying the basic assumptions used in the model. The adsorption sites (heavy dots) are equivalent and can have unit occupancy. Also, the adsorbates are immobile on the surface.

The Langmuir adsorption model explains adsorption by assuming an adsorbate behaves as an ideal gas at isothermal conditions. At these conditions the adsorbate's partial pressure, $p_A$, is related to the volume of it, $V$, adsorbed onto a solid adsorbent. The adsorbent, as indicated in Figure 1, is assumed to be an ideal solid surface composed of series of distinct sites capable of binding the adsorbate. The adsorbate binding is treated as a chemical reaction between the adsorbate molecule $A_{g}$ and an empty site, $S$. This reaction yields an adsorbed complex $A_{ad}$ with an associated equilibrium constant $K_{eq}$

$A_{g}\,\,\, + \,\,\, S \,\,\, \rightleftharpoons \,\,\, A_{ad}$

From these assumptions the Langmuir isotherm can be derived (see below), which states that:

$\theta_A = \frac{V}{V_m} = \frac{K_{eq}^A\,p_A }{1+K_{eq}^A\,p_A}$

where $\theta_A$ is the fractional occupancy of the adsorption sites and $V_m$ is the volume of the monolayer. A continuous monolayer of adsorbate molecules surrounding a homogeneous solid surface is the conceptual basis for this adsorption model.^[1]

The Langmuir isotherm is formally equivalent to the Hill equation in biochemistry.

Background and experiments

In 1916, Irving Langmuir presented his model for the adsorption of species onto simple surfaces. Langmuir was awarded the Nobel Prize in 1932 for his work concerning surface chemistry. He hypothesized that a given surface has a certain number of equivalent sites that a species can “stick”, either by physisorption or chemisorption. His theory began when he postulated that gaseous molecules do not rebound elastically from a surface, but are held by it in a similar way to groups of molecules in solid bodies.^[2]

Langmuir published two papers that proved the assumption that adsorbed films do not exceed one molecule in thickness. The first experiment involved observing electron emission from heated filaments in gases.^[3] The second, a more direct proof, examined and measured the films of liquid on an adsorbent surface layer. He also noted that generally the attractive strength between the surface and the first layer of adsorbed substance is much greater than the strength between the first and second layer. However, there are instances where the subsequent layers may condense given the right combination of temperature and pressure.^[4]

The most important empirical data came from a set of experiments that Langmuir ran to test the adsorption of several gases on mica, glass and platinum. The experiments began at very low pressures (~100 bar) in order to more easily measure the change in quantities of free gas and also to avoid condensation. He then ran the experiments at different temperatures and pressures, which proved the pressure dependence demonstrated below.

Basic assumptions of the model

Inherent within this model, the following assumptions^[5] are valid specifically for the simplest case: the adsorption of a single adsorbate onto a series of equivalent sites on the surface of the solid.

1. The surface containing the adsorbing sites is a perfectly flat plane with no corrugations (assume the surface is homogeneous) .
2. The adsorbing gas adsorbs into an immobile state.
3. All sites are equivalent.
4. Each site can hold at most one molecule of A (mono-layer coverage only).
5. There are no interactions between adsorbate molecules on adjacent sites.

Derivations of the Langmuir Adsorption Isotherm

Kinetic derivation

This section^[5] provides a kinetic derivation for a single adsorbate case. The multiple adsorbate case is covered in the Competitive adsorption sub-section. The model assumes adsorption and desorption as being elementary processes, where the rate of adsorption r_ad and the rate of desorption r_d are given by:

$r_{ad} = k_{ad} \, p_A \, [S]$

$r_{d} = k_d \, [A_{ad}]$

where P_A is the partial pressure of A over the surface, [S] is the concentration of bare sites in number/m², [A_ad] is the surface concentration of A in molecules/m², and k_ad and k_d are constants of forward adsorption reaction and backward desorption reaction in the above reactions.

At equilibrium, the rate of adsorption equals the rate of desorption. Setting r_ad=r_d and rearranging, we obtain:

$\frac {[A_{ad}]}{p_A[S]} = \frac{k_{ad}}{k_d} = K_{eq}^A$

The concentration of all sites [S₀] is the sum of the concentration of free sites [S] and of occupied sites:

$[S_0] = [S] + [A_{ad}]\,$

Combining this with the equilibrium equation, we get:

$[S_0] = \frac {[A_{ad}]}{K_{eq}^A\,p_A} + [A_{ad}] = \frac{1+K_{eq}^A\,p_A}{K_{eq}^A\,p_A}\,[A_{ad}]$

We define now the fraction of the surface sites covered with A, θ_A, as:

$\theta_A = \frac{[A_{ad}]}{[S_0]}$

This, applied to the previous equation that combined site balance and equilibrium, yields the Langmuir adsorption isotherm:

$\theta_A = \frac{K_{eq}^A\,p_A }{1+K_{eq}^A\,p_A}$

Statistical mechanical derivation

This derivation^[6][7] was originally provided by Volmer and Mahnert^[8] in 1925.

The partition function of the finite number of adsorbents adsorbed on a surface, in a canonical ensemble is given by

$Z(N) \, = \, \frac {\zeta^{N}_{L}}{N!} \, \frac { N_{S}!} { (N_{S}-N)!}$

where $\zeta_{L}$ is the partition function of a single adsorbed molecule, $N_{S}$ are the number of sites available for adsorption. Hence, N, which is the number of molecules that can be adsorbed, can be less or equal to N_s. The first term of Z(n) accounts the total partition function of the different molecules by taking a product of the individual partition functions (Refer to Partition function of subsystems). The latter term accounts for the overcounting arising due to the indistinguishable nature of the adsorption sites. The grand canonical partition function is given by

$\mathcal{Z}(\mu) \, = \, \sum_{N=0}^{N_{S}} \, \exp \left(\, \frac {N\mu}{k_{B}T} \right)\frac {\zeta^{N}_{L}}{N!} \, \frac { N_{S}!} { (N_{S}-N)!} \,$

As it has the form of binomial series, the summation is reduced to

$\mathcal{Z} \, = \, (1+x)^{N_{S}}$

where $x \, = \, \zeta_{L}\exp \left( \frac {\mu}{k_{B}T} \right)$

The Landau free energy, which is generalized Helmholtz free energy is given by

$L \, = \, -k_{B}T\ln(\mathcal{Z})\, = \, -k_{B}TN_{S}\ln(1+x)$

According to the Maxwell relations regarding the change of the Helmholtz free energy with respect to the chemical potential,

$\left (\frac {\partial L}{\partial \mu}\right)_{T, V, Area} = \, -N$

which gives

$\theta_{A} \,= \, \frac{N}{N_{S}} \, = \, \frac{x}{1+x}$

Now, invoking the condition that the system is in equilibrium, the chemical potential of the adsorbates is equal to that of the gas surroundings the absorbent.

$\mu_{A} \, = \, \mu_{g}$

An example plot of the surface coverage θ_A = P/(P+P₀) with respect to the partial pressure of the adsorbate. P₀ = 100mtorr. The graph shows levelling off of the surface coverage at pressures higher than P₀.

$\mu_{g} = \left( \frac{\partial A_{g}}{\partial N_i} \right)_{T,V, N_{j \ne i}} \,\, = k_{B}T\ln \frac{N^{3D}}{Z^{3D}}$

where N^3D is the number of gas molecules, Z^3D is the partition function of the gas molecules and A_g=-k_BT ln Z_g. Further, we get

$x \, = \, \frac {\theta_A}{1- \theta_A} \, = \, \zeta_{L} \frac{N^{3D}}{\zeta^{3D}} \, = \, \zeta_L \left ( \frac {h^2}{2 \pi mk_BT} \right)^{3/2} \frac{P}{k_BT} \, = \, \frac{P}{P_0}$

where
$P_0 = \frac {k_BT}{\zeta_L} \left ( \frac {2 \pi mk_BT}{h^2} \right)^{3/2}$

Finally, we have

$\theta_A = \frac {P}{P+P_0}$

It is plotted in the figure alongside demonstrating the surface coverage increases quite rapidly with the partial pressure of the adsorbants but levels off after P reaches P₀.

Competitive Adsorption

The previous derivations assumes that there is only one species, A, adsorbing onto the surface. This section^[9] considers the case when there are two distinct adsorbates present in the system.Consider two species A and B that compete for the same adsorption sites. The following assumptions are applied here:

1. All the sites are equivalent.
2. Each site can hold at most one molecule of A or one molecule of B, but not both.
3. There are no interactions between adsorbate molecules on adjacent sites.

As derived using kinetical considerations, the equilibrium constants for both A and B are given by

$\frac {[A_{ad}]}{p_A\,[S]} = K^A_{eq}$

and

$\frac {[B_{ad}]}{p_B\,[S]} = K^B_{eq}$

The site balance states that the concentration of total sites [S₀] is equal to the sum of free sites, sites occupied by A and sites occupied by B:

$[S_0] = [S] + [A_{ad}] + [B_{ad}]\,$

Inserting the equilibrium equations and rearranging in the same way we did for the single-species adsorption, we get similar expressions for both θ_A and θ_B:

$\theta_A = \frac {K^A_{eq}\,p_A}{1+K^A_{eq}\,p_A+K^B_{eq}\,p_B}$

$\theta_B = \frac {K^B_{eq}\,p_B}{1+K^A_{eq}\,p_A+K^B_{eq}\,p_B}$

Dissociative Adsorption

The other case of special importance is when a molecule D₂ dissociates into two atoms upon adsorption.^[9] Here, the following assumptions would be held to be valid:

1. D₂ completely dissociates to two molecules of D upon adsorption.
2. The D atoms adsorb onto distinct sites on the surface of the solid and then move around and equilibrate.
3. All sites are equivalent.
4. Each site can hold at most one atom of D.
5. There are no interactions between adsorbate molecules on adjacent sites.

Using similar kinetic considerations, we get:

$\frac {[D_{ad}]}{p^{1/2}_{D_2}[S]} = K^D_{eq}$

The 1/2 exponent on p_D2 arises because one gas phase molecule produces two adsorbed species. Applying the site balance as done above:

$\theta_D = \frac {(K^D_{eq}\,p_{D_2})^{1/2}}{1 + (K^D_{eq}\,p_{D_2})^{1/2}}$

Entropic considerations

The formation of Langmuir monolayers by adsorption onto a surface dramatically reduces the entropy of the molecular system. This conflicts with the second law of thermodynamics, which states that entropy will increase in an isolated system. This implies that either another locally active force is stronger than the thermodynamic potential, or that our expression of the entropy of the system is incomplete.

To find the entropy decrease, we find the entropy of the molecule when in the adsorbed condition.^[10]

$S \, = \, S_\text{configurational} \, + \, S_\text{vibrational}$

$S_\text{conf} = k_B \, \ln \Omega_\text{conf}$

$\Omega_\text{conf} \, = \, \frac {N_S!}{N! (N_S -N)!}$

Using Stirling's approximation, we have,

$\ln N! \, = \, N \ln N - N$

$S_\text{conf}/k_B = -\theta_A \, \ln (\theta_A) - (1-\theta_A) \, \ln (1- \theta_A)$

On the other hand, the entropy of a molecule of an ideal gas is
$\frac {S_\text{gas}}{Nk_B} \, = \, \ln \left (\frac {k_BT}{P \lambda^3} \right) + 5/2$

where $\lambda$ is the Thermal de Broglie wavelength of the gas molecule.

Disadvantages of the model

The Langmuir adsorption model deviates significantly in many cases, primarily because it fails to account for the surface roughness of the adsorbent. Rough inhomogeneous surfaces have multiple site-types available for adsorption, and some parameters vary from site to site, such as the heat of adsorption.

The model also ignores adsorbate/adsorbate interactions. Experimentally, there is clear evidence for adsorbate/adsorbate interactions in heat of adsorption data. There are two kinds of adsorbate/adsorbate interactions: direct interaction and indirect interaction. Direct interactions are between adjacent adsorbed molecules, which could make adsorbing near another adsorbate molecule more or less favorable and greatly affects high-coverage behavior. In indirect interactions, the adsorbate changes the surface around the adsorbed site, which in turn affects the adsorption of other adsorbate molecules nearby.

Modifications of the Langmuir Adsorption Model

The modifications try to account for the points mentioned in above section like surface roughness, inhomogeneity, and adsorbate-adsorbate interactions.

The Freundlich Adsorption Isotherm

The Freundlich isotherm is the most important multisite adsorption isotherm for rough surfaces.

$\theta_A = \alpha_F\,p^{C_F}$

where α_F and C_F are fitting parameters.^[11] This equation implies that if one makes a log-log plot of adsorption data, the data will fit a straight line. The Freundlich isotherm has two parameters while Langmuir's equations has only one: as a result, it often fits the data on rough surfaces better than the Langmuir's equations.

A related equation is the Toth equation. Rearranging the Langmuir equation, one can obtain:

$\theta_A = \frac{p_A}{\frac{1}{K_{eq}^A} + p_A}$

Toth^[12] modified this equation by adding two parameters, α_T0 and C_T0 to formulate the Toth equation:

$\theta^{C_{T_0}} = \frac {\alpha_{T_0}\,p_A^{C_{T_0}}}{\frac{1}{K_{eq}^A} + p_A^{C_{T_0}}}$

The Temkin Adsorption Isotherm

This isotherm takes into accounts of indirect adsorbate-adsorbate interactions on adsorption isotherms. Temkin^[13] noted experimentally that heats of adsorption would more often decrease than increase with increasing coverage.

The heat of adsorption ΔH_ad is defined as:

$\frac{[A_{ad}]}{p_A\,[S]} = K^A_{eq} \propto \mathrm{e}^{-\Delta G_{ad}/RT} = \mathrm{e}^{\Delta S_{ad}/R}\,\mathrm{e}^{-\Delta H_{ad}/RT}$

He derived a model assuming that as the surface is loaded up with adsorbate, the heat of adsorption of all the molecules in the layer would decrease linearly with coverage due to adsorbate/adsorbate interactions:

$\Delta H_{ad} = \Delta H^0_{ad}\,(1-\alpha_T\,\theta)$

where α_T is a fitting parameter. Assuming the Langmuir Adsorption isotherm still applied to the adsorbed layer, $K_{eq}^A$ is expected to vary with coverage, as follows:

$K^A_{eq} = K^{A,0}_{eq} \mathrm{e}^{(\Delta H^0_{ad}\,(1-\alpha_T\,\theta)/ k\,T)}$

Langmuir's isotherm can be rearranged to this form:

$K^A_{eq}\,p_A = \frac{\theta }{1-\theta}$

Substituting the expression of the equilibrium constant and taking the natural logarithm:

$\ln (K^{A,0}_{eq}\,p_A) = \frac{-\Delta H^0_{ad} \, \alpha_T \, \theta}{k\,T} + \ln \left( \frac{\theta}{1-\theta}\right)$

BET equation

Brunauer's model of multilayer adsorption, that is, a random distribution of sites covered by one, two, three, etc., adsorbate molecules.

Brunauer, Emmett and Teller^[14] derived the first isotherm for multilayer adsorption. It assumes a random distribution of sites that are empty or that are covered with by one monolayer, two layers and so on, as illustrated alongside. The main equation of this model is:

$\frac{[A]}{S_0} = \frac{c_B \, x_B}{(1-x_B)\,[1 + (c_B - 1)\,x_B]}$

where

$x_B = p_A\,K_m, \qquad c_B = \frac{K_1}{K_m}$

and [A] is the total concentration of molecules on the surface, given by:

$[A] = \sum^{\infty}_{i=1} i\,[A]_i = \sum^{\infty}_{i=1}i \, K_1 \, K^{i-1}_m \, p^i_A \, [A]_0$

where

$K_i = \frac{[A]_i}{p_A\,[A]_{i-1}}$

in which [A]₀ is the number of bare sites, and [A]_i is the number of surface sites covered by i molecules.

Adsorption of binary liquid adsorption on solids

This section describes the surface coverage when the adsorbate is in liquid phase and is a binary mixture^[15]

For ideal both phases - no lateral interactions, homogeneous surface - the composition of a surface phase for a binary liquid system in contact with solid surface is given by a classic Everett isotherm equation (being a simple analogue of Langmuir equation), where the components are interchangeable (i.e. "1" may be exchanged to "2") without change of eq. form:

$x_1^s \, = \, \frac{Kx_1^l}{1+(K-1)x_1^l}$

where the normal definition of multicomponent system is valid as follows :

$\sum_{i=1}^{k} x^s_i =1 \,\,\sum_{i=1}^{k} x^l_i =1$

By simple rearrangement, we get

$x_1^s \, = \, \frac{K[x_1^l/(1-x_1^l)]}{1+K[x_1^l/(1-x_1^l)]}$

This equation describes competition of components "1" and "2".

References

1. Hanaor, D.A.H.; Ghadiri, M.; Chrzanowski, W.; Gan, Y. (2014). "Scalable Surface Area Characterization by Electrokinetic Analysis of Complex Anion Adsorption" (PDF). Langmuir 30 (50): 15143–15152. doi:10.1021/la503581e.
2. Langmuir, Irving (June 1918). "The Adsorption of Gases on Plane Surface of Glass, Mica and Platinum". The Research Laboratory of The General Electric Company: 1361–1402. doi:10.1021/ja02242a004. Retrieved 11 June 2013.
3. Langmuir, Irving (1916). "Part I". The Research Laboratory of The General Electric Company: 2221.
4. Langmuir, Irving (1918). "Part II". The Research Laboratory of The General Electric Company: 1848.
5. Masel, Richard (1996). Principles of Adsorption and Reaction on Solid Surfaces. Wiley Interscience. p. 240. ISBN 0-471-30392-5.
6. Masel, Richard (1996). Principles of Adsorption and Reaction on Solid Surfaces. Wiley Interscience. p. 242. ISBN 0-471-30392-5.
7. Cahill, David (2008). "Lecture Notes 5 Page 2" (pdf). University of Illinois, Urbana Champaign. Retrieved 2008-11-09.
8. Volmer, M.A., and P. Mahnert, Z. Physik. Chem 115, 253
9. Masel, Richard (1996). Principles of Adsorption and Reaction on Solid Surfaces. Wiley Interscience. p. 244. ISBN 0-471-30392-5.
10. Cahill, David (2008). "Lecture Notes 5 Page 13" (pdf). University of Illinois, Urbana Champaign. Retrieved 2008-11-09.
11. Freundlich, H. (1909). "eine darstellung der chemie der kolloide und verwanter gebiete.". Kapillarchemie (Leipzig: Academishe Bibliotek).
12. Toth, J., Acta. Chim. Acad. Sci. Hung 69, 311(1971)
13. Temkin, M. I.; Pyzhev, V. (1940). Acta Physicochima USSR 12: 327. Missing or empty |title= (help)
14. Brunauer, Stephen; Emmett, P. H.; Teller, Edward (1938). "Adsorption of Gases in Multimolecular Layers". Journal of the American Chemical Society 60 (2): 309–319. doi:10.1021/ja01269a023. ISSN 0002-7863.
15. Marczewski, A. W (2002). "Basics of Liquid Adsorption". Retrieved 2008-11-24.

• The constitution and fundamental properties of solids and liquids. part i. solids. Irving Langmuir; J. Am. Chem. Soc. 38, 2221-95 1916

External links

• Langmuir isotherm from Queen Mary, University of London
• LMMpro, Langmuir equation-fitting software