Reactions on surfaces

Reactions on surfaces are reactions in which at least one of the steps of the reaction mechanism is the adsorption of one or more reactants. The mechanisms for these reactions, and the rate equations are of extreme importance for heterogeneous catalysis. Via scanning tunneling microscopy, it is possible to observe reactions at the solid gas interface in real space, if the time scale of the reaction is in the correct range. Reactions at the solid–gas interface are in some cases related to catalysis.

Simple decomposition
If a reaction occurs through these steps:


 * A + S ⇌ AS → Products

where A is the reactant and S is an adsorption site on the surface and the respective rate constants for the adsorption, desorption and reaction are k1, k−1 and k2, then the global reaction rate is:


 * $$r=k_2 C_\mathrm{AS}=k_2 \theta C_\mathrm{S} $$

where:
 * r is the rate, mol·m−2·s−1
 * $$C_A$$is the concentration of adsorbate, mol·m−3
 * $$C_\mathrm{AS}$$ is the surface concentration of occupied sites, mol·m−2
 * $$C_\mathrm{S}$$ is the concentration of all sites (occupied or not), mol·m−2
 * $$\theta$$ is the surface coverage, (i.e. $$C_{AS}/C_{S}$$) defined as the fraction of sites which are occupied, which is dimensionless
 * $$t$$ is time, s
 * $$k_2$$ is the rate constant for the surface reaction, s−1.
 * $$k_1$$is the rate constant for surface adsorption, m3·mol−1·s−1
 * $$k_{-1}$$is the rate constant for surface desorption, s−1

$$C_\mathrm{S}$$ is highly related to the total surface area of the adsorbent: the greater the surface area, the more sites and the faster the reaction. This is the reason why heterogeneous catalysts are usually chosen to have great surface areas (in the order of a hundred m2/gram)

If we apply the steady state approximation to AS, then:


 * $$\frac {dC_\mathrm{AS}}{dt}= 0 = k_1 C_\mathrm{A} C_\mathrm{S} (1-\theta)- k_2 \theta C_\mathrm{S} -k_{-1}\theta C_\mathrm{S} $$ so $$\theta =\frac {k_1 C_\mathrm{A}}{k_1 C_\mathrm{A} + k_{-1}+k_2}$$

and


 * $$r=\frac {k_1 k_2 C_\mathrm{A} C_\mathrm{S}}{k_1 C_\mathrm{A} + k_{-1}+k_2}.$$

The result is equivalent to the Michaelis–Menten kinetics of reactions catalyzed at a site on an enzyme. The rate equation is complex, and the reaction order is not clear. In experimental work, usually two extreme cases are looked for in order to prove the mechanism. In them, the rate-determining step can be:
 * Limiting step: adsorption/desorption
 * $$k_2 \gg \ k_1C_\mathrm{A}, k_{-1},\text{ so }r \approx k_1 C_\mathrm{A} C_\mathrm{S}.$$

The order respect to A is 1. Examples of this mechanism are N2O on gold and HI on platinum
 * Limiting step: reaction of adsorbed species


 * $$k_2 \ll \ k_1C_\mathrm{A}, k_{-1}\text{ so }\theta =\frac {k_1 C_\mathrm{A}}{k_1 C_\mathrm{A} + k_{-1}} =\frac {K_1 C_\mathrm{A}}{K_1 C_\mathrm{A} + 1}$$

The last expression is the Langmuir isotherm for the surface coverage. The adsorption equilibrium constant $$K_1=\frac{k_1}{k_{-1}}$$, and the numerator and denominator have each been divided by $$k_{-1}$$. The overall reaction rate becomes $$r= \frac {K_1 k_2 C_\mathrm{A} C_\mathrm{S}}{K_1 C_\mathrm{A}+1}$$.

Depending on the concentration of the reactant the rate changes:
 * Low concentrations, then $$r= K_1 k_2 C_\mathrm{A} C_\mathrm{S}$$, that is to say a first order reaction in component A.
 * High concentration, then $$r= k_2 C_\mathrm{S}$$. It is a zeroth order reaction in component A.

Langmuir–Hinshelwood mechanism
In this mechanism, suggested by Irving Langmuir in 1921 and further developed by Cyril Hinshelwood in 1926, two molecules adsorb on neighboring sites and the adsorbed molecules undergo a bimolecular reaction:


 * A + S ⇌ AS


 * B + S ⇌ BS


 * AS + BS → Products

The rate constants are now $$k_1$$,$$k_{-1}$$,$$k_2$$,$$k_{-2}$$ and $$k$$ for adsorption/desorption of A, adsorption/desorption of B, and reaction. The rate law is: $$r=k \theta_\mathrm{A} \theta_\mathrm{B} C_\mathrm{S}^2 $$

Proceeding as before we get $$\theta_\mathrm{A}=\frac{k_1C_\mathrm{A}\theta_E}{k_{-1}+kC_\mathrm{S}\theta_\mathrm{B}}$$, where $$\theta_E$$ is the fraction of empty sites, so $$\theta_\mathrm{A}+\theta_\mathrm{B}+\theta_E=1$$. Let us assume now that the rate limiting step is the reaction of the adsorbed molecules, which is easily understood: the probability of two adsorbed molecules colliding is low. Then $$\theta_\mathrm{A}=K_1C_\mathrm{A}\theta_E$$, with $$K_i=k_i/k_{-i}$$, which is nothing but Langmuir isotherm for two adsorbed gases, with adsorption constants $$K_1$$ and $$K_2$$. Calculating $$\theta_E$$ from $$\theta_\mathrm{A}$$ and $$\theta_\mathrm{B}$$ we finally get
 * $$r=k C_\mathrm{S}^2 \frac{K_1K_2C_\mathrm{A}C_\mathrm{B}}{(1+K_1C_\mathrm{A}+K_2C_\mathrm{B})^2}$$.

The rate law is complex and there is no clear order with respect to either reactant, but we can consider different values of the constants, for which it is easy to measure integer orders: That means that $$1 \gg K_1C_\mathrm{A}, K_2C_\mathrm{B}$$, so $$r=k C_\mathrm{S}^2 K_1K_2C_\mathrm{A}C_\mathrm{B}$$. The order is one with respect to each reactant, and the overall order is two. In this case $$K_1C_\mathrm{A}, 1 \gg K_2C_\mathrm{B}$$, so $$r=k C_\mathrm{S}^2 \frac{K_1K_2C_\mathrm{A}C_\mathrm{B}}{(1+K_1C_\mathrm{A})^2}$$. The reaction order is 1 with respect to B. There are two extreme possibilities for the order with respect to A:
 * Both molecules have low adsorption
 * One molecule has very low adsorption
 * At low concentrations of A, $$r=k C_\mathrm{S}^2 K_1K_2C_\mathrm{A}C_\mathrm{B}$$, and the order is one with respect to A.
 * At high concentrations, $$r=k C_\mathrm{S}^2 \frac{K_2C_\mathrm{B}}{K_1C_\mathrm{A}}$$. The order is minus one with respect to A. The higher the concentration of A, the slower the reaction goes, in this case we say that A inhibits the reaction.

One of the reactants has very high adsorption and the other one doesn't adsorb strongly.
 * One molecule has very high adsorption

$$K_1C_\mathrm{A} \gg 1, K_2C_\mathrm{B}$$, so $$r=k C_\mathrm{S}^2 \frac{K_2C_\mathrm{B}}{K_1C_\mathrm{A}}$$. The reaction order is 1 with respect to B and −1 with respect to A. Reactant A inhibits the reaction at all concentrations.

The following reactions follow a Langmuir–Hinshelwood mechanism:
 * 2 CO + O2 → 2 CO2 on a platinum catalyst.
 * CO + 2H2 → CH3OH on a ZnO catalyst.
 * C2H4 + H2 → C2H6 on a copper catalyst.
 * N2O + H2 → N2 + H2O on a platinum catalyst.
 * C2H4 + $1⁄2$ O2 → CH3CHO on a palladium catalyst.
 * CO + OH → CO2 + H+ + e− on a platinum catalyst.

Eley–Rideal mechanism
In this mechanism, proposed in 1938 by D. D. Eley and E. K. Rideal, only one of the molecules adsorbs and the other one reacts with it directly from the gas phase, without adsorbing ("nonthermal surface reaction"):


 * A(g) + S(s) ⇌ AS(s)


 * AS(s) + B(g) → Products

Constants are $$k_1, k_{-1}$$ and $$k$$ and rate equation is $$r = k C_\mathrm{S} \theta_\mathrm{A} C_\mathrm{B}$$. Applying steady state approximation to AS and proceeding as before (considering the reaction the limiting step once more) we get $$r=kC_\mathrm{S} C_\mathrm{B}\frac{K_1C_\mathrm{A}}{K_1C_\mathrm{A}+1}$$. The order is one with respect to B. There are two possibilities, depending on the concentration of reactant A:


 * At low concentrations of A, $$r=kC_\mathrm{S} K_1C_\mathrm{A}C_\mathrm{B}$$, and the order is one with respect to A.


 * At high concentrations of A, $$r=kC_\mathrm{S} C_\mathrm{B}$$, and the order is zero with respect to A.

The following reactions follow an Eley–Rideal mechanism:
 * C2H4 + $1⁄2$ O2 (adsorbed) → (CH2CH2)O The dissociative adsorption of oxygen is also possible, which leads to secondary products carbon dioxide and water.
 * CO2 + H2 (ads.) → H2O + CO
 * 2 NH3 + $1 1⁄2$ O2 (ads.) → N2 + 3H2O on a platinum catalyst
 * C2H2 + H2 (ads.) → C2H4 on nickel or iron catalysts