6.3 BET Theory and Surface Area Determination

BET Theory for Multilayer Adsorption

Principles and Assumptions

The Langmuir model assumes adsorption stops at a single monolayer, but real surfaces often accumulate multiple layers of adsorbed gas. BET theory (Brunauer, Emmett, and Teller) extends the Langmuir framework to account for this multilayer adsorption, making it the standard method for determining the specific surface area of porous materials.

BET theory rests on several key assumptions:

• Gas molecules can adsorb onto a solid surface in an infinite number of layers, with no interaction between layers (only vertical stacking, no lateral molecule-molecule interactions within a layer).
• The first layer adsorbs with a characteristic heat of adsorption that reflects the adsorbate-surface interaction. All subsequent layers adsorb with a heat equal to the heat of liquefaction of the adsorbate, since molecules in upper layers are essentially condensing onto other adsorbate molecules rather than onto the bare surface.
• Each layer reaches its own adsorption-desorption equilibrium independently.
• The concept of a statistical thickness of the adsorbed film emerges naturally: it's the ratio of the total amount adsorbed to the amount needed to complete a single monolayer.

Applicability and Adsorption Isotherms

BET theory works best for Type II and Type IV adsorption isotherms (IUPAC classification):

• Type II isotherms are characteristic of non-porous or macroporous materials (pore diameters > 50 nm). These show unrestricted monolayer-multilayer adsorption with a smooth, sigmoidal curve. The "knee" of the curve corresponds roughly to monolayer completion.
• Type IV isotherms are characteristic of mesoporous materials (pore diameters between 2 and 50 nm). They resemble Type II at low pressures but exhibit a hysteresis loop at higher pressures due to capillary condensation within the mesopores.

For microporous materials (pore diameters < 2 nm), pore filling dominates rather than layer-by-layer adsorption, so BET analysis becomes unreliable.

BET Equation Derivation and Limitations

Derivation of the BET Equation

The BET equation is derived by writing rate expressions for condensation and evaporation at each adsorption layer, then summing over all layers at equilibrium. The result is:

$\frac{1}{v\left(\frac{P_0}{P} - 1\right)} = \frac{c - 1}{v_m c} \cdot \frac{P}{P_0} + \frac{1}{v_m c}$

where:

• $v$ = volume of gas adsorbed at pressure $P$
• $P_0$ = saturation vapor pressure of the adsorbate at the experimental temperature
• $v_m$ = volume of gas required to form exactly one monolayer
• $c$ = BET constant, related to the difference between the heat of adsorption of the first layer and the heat of liquefaction

This equation has the form $y = mx + b$, which means you can linearize it:

1. Measure $v$ at several values of $P/P_0$.

2. Plot $\frac{1}{v\left(\frac{P_0}{P} - 1\right)}$ on the y-axis versus $\frac{P}{P_0}$ on the x-axis.

3. Fit a straight line through the linear region.

4. Extract the slope $s = \frac{c-1}{v_m c}$ and intercept $i = \frac{1}{v_m c}$.

5. Solve for the two unknowns:

   • $v_m = \frac{1}{s + i}$
   • $c = \frac{s}{i} + 1$

Assumptions and Limitations

• The surface is assumed to be energetically homogeneous, meaning every adsorption site on the bare surface has the same energy. Real surfaces have defects, edges, and varying crystal faces, so this is an approximation.
• No lateral interactions between adsorbed molecules within a given layer are considered.
• The BET equation is valid only over a limited range of relative pressures, typically $0.05 \leq P/P_0 \leq 0.35$. Below this range, surface heterogeneity effects distort the fit. Above it, capillary condensation or pore-filling effects take over.
• For microporous materials, the BET surface area can be physically meaningless because adsorption proceeds by pore filling rather than multilayer formation. Methods like Dubinin-Radushkevich are more appropriate in that case.

Surface Area Determination Using BET

Specific Surface Area Calculation

Once you've obtained $v_m$ from the BET plot, the specific surface area follows from a straightforward geometric argument: convert the monolayer volume to a number of molecules, then multiply by the area each molecule occupies.

$S_{\text{BET}} = \frac{v_m \cdot N_A \cdot \sigma}{V_m \cdot m}$

where:

• $N_A$ = Avogadro's number ($6.022 \times 10^{23}$ mol$^{-1}$)
• $\sigma$ = cross-sectional area of one adsorbate molecule
• $V_m$ = molar volume of the gas at STP (22,414 cm$^3$/mol for an ideal gas)
• $m$ = mass of the sample

Nitrogen ($N_2$) is the most commonly used probe gas, with $\sigma = 0.162$ nm$^2$ per molecule at 77 K. Other adsorbates (krypton, argon) are sometimes used for low-surface-area materials where higher sensitivity is needed.

Experimental Procedure

A typical BET measurement follows these steps:

1. Degas the sample by heating it under vacuum (or flowing inert gas) to remove adsorbed water, organics, and other contaminants. The temperature and duration depend on the material's thermal stability.
2. Cool the sample to the adsorption temperature. For $N_2$ adsorption, this is 77 K (liquid nitrogen bath).
3. Dose the sample incrementally with known quantities of $N_2$ gas using either a volumetric (manometric) or gravimetric apparatus.
4. Wait for equilibrium after each dose, then record the equilibrium pressure. This gives one data point on the adsorption isotherm ($v$ vs. $P/P_0$).
5. Repeat across the desired pressure range to build the full isotherm.
6. Construct the BET plot using data in the $0.05 \leq P/P_0 \leq 0.35$ window and extract $v_m$ and $c$.

Importance of Specific Surface Area

Specific surface area directly influences how a material performs in applications like catalysis, gas adsorption, and energy storage. More surface area means more accessible sites for reactions or adsorption events.

• Activated carbons can exceed 1000 m$^2$/g, making them effective adsorbents for water purification and gas separation.
• Zeolites and metal-organic frameworks (MOFs) can reach even higher values, sometimes above 5000 m$^2$/g for certain MOFs.
• Comparing BET surface areas across synthesis batches or different materials provides a quantitative basis for optimizing preparation conditions.

BET Isotherm Analysis and Surface Area Calculation

Adsorption Isotherm Interpretation

The raw adsorption isotherm plots $v$ (volume adsorbed, often reported at STP) against $P/P_0$ at constant temperature. The shape tells you about the material:

• A Type II curve with a clear knee near $P/P_0 \approx 0.05\text{–}0.15$ suggests monolayer completion on a non-porous or macroporous solid.
• A Type IV curve with a hysteresis loop at $P/P_0 \approx 0.4\text{–}0.8$ signals mesopores undergoing capillary condensation.
• A very steep uptake at low $P/P_0$ (Type I character) suggests significant microporosity, and BET analysis should be applied cautiously.

The linear region of the BET plot (within $0.05 \leq P/P_0 \leq 0.35$) is where you extract $v_m$ and $c$ from the slope and intercept.

Surface Area Calculation from Experimental Data

Here's a worked example for nitrogen adsorption at 77 K:

Given:

• $v_m = 10.0$ cm$^3$/g (from BET plot, at STP)
• $N_A = 6.022 \times 10^{23}$ mol$^{-1}$
• $\sigma = 0.162$ nm$^2$ = $0.162 \times 10^{-18}$ m$^2$
• $V_m = 22{,}414$ cm$^3$/mol

Calculation:

$S_{\text{BET}} = \frac{v_m \cdot N_A \cdot \sigma}{V_m} = \frac{10.0 \times 6.022 \times 10^{23} \times 0.162 \times 10^{-18}}{22{,}414}$

$S_{\text{BET}} = \frac{9.756 \times 10^{5} \times 10^{-18} \text{ m}^2}{22{,}414} \approx 43.5 \text{ m}^2/\text{g}$

Note that $v_m$ here is already per gram of sample, so the mass doesn't appear separately. If $v_m$ were reported as a total volume for a sample of mass $m$, you'd divide by $m$ as well.

Additional Information from BET Analysis

The BET constant $c$ carries physical meaning beyond being a fitting parameter:

• Large $c$ values (e.g., > 100) indicate strong adsorbate-surface interactions and a sharp knee in the isotherm. The monolayer is well-defined.
• Small $c$ values (e.g., < 20) suggest weak interactions, and the transition from monolayer to multilayer is gradual. The BET surface area becomes less reliable in this regime.

Deviations from linearity in the BET plot can guide you toward alternative methods:

• At low $P/P_0$: upward curvature may indicate micropore filling. Consider the Dubinin-Radushkevich or Horvath-Kawazoe methods.
• At high $P/P_0$: downward curvature often reflects capillary condensation in mesopores. The Barrett-Joyner-Halenda (BJH) method uses the Kelvin equation (assuming cylindrical pores) to extract a pore size distribution from the desorption branch of the isotherm.

These complementary methods, combined with BET surface area, give a more complete picture of a material's pore structure and adsorption characteristics.