Catalyst Development
Surface area directly correlates with catalytic activity. BET monitors catalyst deactivation through surface area loss. DFT reveals micropore blockage in zeolite catalysts.
Last reviewed: May 8, 2026
Pillar II: Surface Area & Micropore/Mesopore Analysis
Gas Adsorption BET & BJH
Gas adsorption analysis quantifies surface area via BET theory and pore size distribution through BJH, t-plot, and DFT methods. By measuring physisorption isotherms with nitrogen (77 K), argon (87 K), or krypton, this technique characterizes micropores (<2 nm) and mesopores (2-50 nm) with surface areas from 0.01 to 3,000+ m²/g.
0.35 nm
Min Pore Size
300 nm
Max Pore Size
3000+ m²/g
Max Surface Area
ISO 9277
Key Standard
On This Page
- 1. BET Theory & Surface Area Determination
- 2. IUPAC Isotherm Types (I-VI) & Interpretation
- 3. Hysteresis Loops (H1-H5) & Pore Geometry
- 4. BJH Method for Mesopore Analysis
- 5. Micropore Analysis: t-Plot & Dubinin Methods
- 6. DFT/NLDFT: Modern Pore Size Analysis
- 7. Adsorbate Selection: N₂ vs. Ar vs. Kr
- 8. BET Validation & C-Constant Interpretation
- 9. Sample Preparation & Degassing
- 10. Applications & Case Studies
Fundamental Theory
BET Theory & Surface Area Determination
The Brunauer-Emmett-Teller (BET) theory, developed in 1938, extends Langmuir's monolayer adsorption to multilayer physisorption. It remains the most widely used method for determining specific surface area, despite its simplifying assumptions.
BET Equation (Linear Form)
1/[W(P₀/P - 1)] = 1/(Wₘ·C) + (C-1)/(Wₘ·C) · (P/P₀)
Where:
W = weight of gas adsorbed
Wₘ = monolayer capacity
P/P₀ = relative pressure
C = BET constant (~exp(E₁-E_L)/RT)
Key Assumptions:
• Uniform surface
• No lateral interactions
• Infinite layer formation
• E₁ = heat of first layer
• E_L = heat of liquefaction
Surface Area Calculation
From Monolayer to Surface Area
S_BET = (Wₘ × N_A × A_m) / M
• N_A = Avogadro's number (6.022 × 10²³)
• A_m = molecular cross-sectional area
- N₂ at 77 K: 16.2 Ų
- Ar at 87 K: 13.8 Ų (or 14.2 Ų on heterogeneous surfaces)
- Kr at 77 K: 20.2 Ų
• M = molecular weight of adsorbate
BET Plot Construction
- Collect isotherm data: Measure adsorption at multiple P/P₀ points (typically 5-10 points)
-
Select linear range: Usually P/P₀ = 0.05-0.30 for mesoporous materials
- Microporous: 0.01-0.10
- Mesoporous: 0.05-0.30
- Macroporous: 0.05-0.35
- Plot BET transform: 1/[W(P₀/P-1)] vs. P/P₀
- Linear regression: ISO 9277 asks for a correlation coefficient of at least 0.999 over the chosen points
-
Extract parameters:
• Slope = (C-1)/(Wₘ·C)
• Intercept = 1/(Wₘ·C)
• Wₘ = 1/(slope + intercept)
⚠️ Critical Validation:
The C-constant must be positive (typically 50-300). Negative or very low C values (<20) indicate invalid BET application. Very high C values (>1000) suggest strong microporosity where BET may overestimate surface area.
Classification System
IUPAC Isotherm Types & Interpretation
The 2015 IUPAC classification defines six physisorption isotherm types, each revealing specific pore structures and surface interactions. Understanding these patterns is essential for proper material characterization.
TYPE I - Microporous
Langmuir-type isotherm with rapid initial uptake followed by plateau. Subdivisions:
- • Type I(a): Narrow micropores < 1 nm
- • Type I(b): Wider micropores & narrow mesopores < 2.5 nm
Examples: Activated carbons, zeolites, MOFs
TYPE II - Non-porous/Macroporous
Unrestricted monolayer-multilayer adsorption. Point B marks monolayer completion. Typical S-shaped curve.
BET analysis most reliable for this type.
Examples: Non-porous powders, macroporous materials
TYPE III - Weak Interactions
No identifiable monolayer formation. Convex curve throughout. Adsorbate-adsorbate > adsorbent-adsorbate interactions.
BET not applicable.
Examples: Water on hydrophobic surfaces
TYPE IV - Mesoporous
Capillary condensation with hysteresis. Subdivisions:
- • Type IV(a): With hysteresis (pores > 4 nm)
- • Type IV(b): No hysteresis (pores < 4 nm)
Examples: MCM-41, SBA-15, mesoporous oxides
TYPE V - Weak + Pores
Similar to Type III initially but with pore filling at higher P/P₀. Weak adsorbent-adsorbate interactions.
Rare; water on hydrophobic microporous carbons.
TYPE VI - Layer-by-Layer
Stepwise multilayer adsorption on uniform surface. Each step represents a completed layer.
Examples: Ar or Kr on graphitized carbon black
Isotherm Analysis Guidelines
- Initial slope (P/P₀ < 0.05): Steep = micropores; gradual = weak interaction
- Knee sharpness: Sharp = strong adsorbate-surface interaction (high C-constant)
- Middle section (0.1-0.8): Linear = monolayer-multilayer; curved = capillary condensation
- High P/P₀ behavior: Plateau = limited pore volume; continuous rise = macropores/interparticle voids
- Hysteresis presence: For N₂ at 77 K, hysteresis appears when mesopores are wider than about 4 nm (cylindrical pores); its shape reflects pore geometry and connectivity
Pore Structure Analysis
Hysteresis Loops & Pore Geometry
Hysteresis between adsorption and desorption branches reveals pore shape, size distribution, and network connectivity. The 2015 IUPAC classification defines five hysteresis types (H1-H5), each associated with specific pore geometries. The dedicated hysteresis loops H1–H5 reading guide walks through each loop class, the pore geometry it indicates, and which branch (adsorption or desorption) is the safer input to a Kelvin / BJH analysis.
H1 - Uniform Pores
Narrow distribution of uniform mesopores. Steep, narrow hysteresis loop with parallel branches. Condensation is delayed on the adsorption branch; desorption is close to equilibrium evaporation.
Geometry: Cylindrical, uniform channels
Examples: MCM-41, SBA-15, controlled pore glass
H2 - Complex Networks
Network effects with pore blocking. Subdivisions based on desorption branch steepness:
- • H2(a): Steep desorption - narrow neck distribution
- • H2(b): Gradual desorption - wide neck distribution
Geometry: Ink-bottle pores
Examples: Silica gels, porous glasses
H3 - Slit-Shaped Pores
No limiting adsorption at high P/P₀. Non-rigid aggregates forming slit-like pores. The lower closure point is usually at the cavitation-controlled p/p₀ ≈ 0.4–0.45 (N₂, 77 K).
Geometry: Plate-like particles
Examples: Clays, layered materials
H4 - Narrow Slit Pores
Similar to H3 but with more pronounced uptake at low P/P₀. Indicates micropores in mesoporous materials.
Geometry: Narrow slits
Examples: Activated carbons, zeolites
H5 - Partial Pore Blocking
Distinctive step in desorption branch. Open pores partially blocked but not completely. Unusual loop shape.
Geometry: Plugged hexagonal templates
Examples: Certain mesoporous silicas
Cavitation and Tensile Strength Effects
For nitrogen at 77 K, cavitation occurs at P/P₀ ≈ 0.42-0.45, causing forced closure of hysteresis loops. This represents the tensile strength limit of liquid nitrogen, not a pore size effect.
Implications:
• Pores with necks < 5-6 nm show cavitation-controlled desorption
• True pore size cannot be determined from desorption branch
• Argon at 87 K shows cavitation at P/P₀ ≈ 0.38
• Consider using argon or alternative methods for narrow-necked pores
Classical Method
BJH Method for Mesopore Analysis
The Barrett-Joyner-Halenda (BJH) method, developed in 1951, calculates pore size distribution from the desorption branch using the Kelvin equation. While widely used, it has known limitations for pores below 10 nm.
Modified Kelvin Equation
r_k = -2γV_m cosθ / RT ln(P/P₀)
Where r_k is the Kelvin radius (core radius), to which the statistical thickness t must be added:
r_p = r_k + t
BJH Calculation Procedure
- Start from highest pressure: Work backwards along desorption branch
- Calculate Kelvin radius: r_k from current P/P₀
-
Add statistical thickness: t = thickness of adsorbed layer
- Harkins-Jura: t = [13.99/(0.034 - log(P/P₀))]^0.5
- Halsey: t = 3.54 × [5/ln(P₀/P)]^(1/3)
- Calculate pore volume: From desorbed volume minus film thinning
- Correct for area change: Account for newly exposed surface area
BJH Limitations & Corrections
⚠️ Known BJH Limitations:
- • Underestimates pore size by 20-30% for pores < 10 nm
- • Assumes cylindrical pore geometry
- • Neglects fluid-wall interactions in small pores
- • Statistical thickness models may be inaccurate
- • Tensile strength effects below P/P₀ = 0.42
Micropore Methods
Micropore Analysis: t-Plot & Dubinin Methods
Micropores (<2 nm) require specialized analysis methods since BET and BJH fail at these scales. The t-plot and Dubinin-Radushkevich methods provide micropore volume and surface area determination.
t-Plot Method
The t-plot compares sample adsorption to a reference non-porous material with similar surface chemistry. Deviations from linearity indicate microporosity.
Plot Construction
- Calculate statistical thickness t for each P/P₀
- Plot V_ads vs. t
- Identify linear regions
- Extrapolate to t = 0
Interpretation
- • Positive intercept = micropore volume
- • Slope = external surface area
- • Downward deviation = micropore filling
- • Upward deviation = capillary condensation
Dubinin-Radushkevich (DR) Equation
log W = log W₀ − D·[log(P₀/P)]², with D = 2.303·B·(RT)²/β²
• W = amount adsorbed at P/P₀
• W₀ = micropore capacity
• B = structural constant related to pore size; T = analysis temperature
• β = affinity coefficient relative to benzene (N₂ ≈ 0.33, CO₂ ≈ 0.35)
• Applied to the low-pressure (micropore-filling) part of the isotherm; the linear range must be checked for each sample
αs-Plot Method
Superior to t-plot as it uses reduced adsorption (αs = n/n₀.₄), so no statistical-thickness equation is needed; the reference should match the surface chemistry of the sample.
- • Plot V_ads vs. αs
- • Reference: non-porous material at P/P₀ = 0.4
- • More robust for heterogeneous surfaces
- • Standard references available for carbons, silicas, aluminas
Modern Methods
DFT/NLDFT: State-of-the-Art Pore Analysis
Density Functional Theory (DFT) methods represent the current gold standard for pore size analysis, providing molecular-level accuracy across the entire micro-mesopore range. Unlike classical methods, DFT accounts for fluid-wall interactions and confined space effects.
DFT Method Variants
NLDFT - Non-Local DFT
Accounts for fluid density oscillations near pore walls. Assumes smooth, chemically homogeneous surfaces.
- • Best for: Ordered mesoporous materials (MCM-41, SBA-15)
- • Limitations: Overestimates narrow microporosity in heterogeneous materials
- • Kernels available: N₂@77K, Ar@87K on silica, carbon
QSDFT - Quenched Solid DFT
Accounts for surface roughness and chemical heterogeneity through quenched solid density profile.
- • Best for: Activated carbons, heterogeneous materials
- • Advantages: More realistic for real materials
- • Introduced by Ravikovitch & Neimark (2006); now widely used for carbons
2D-NLDFT - Two-Dimensional
Two-dimensional NLDFT (Jagiello & Olivier, 2013) treats energetic heterogeneity and geometric corrugation of carbon pore walls.
- • Best for: microporous and micro-mesoporous carbons
- • Kernel-specific: results are only as good as the model assumptions
- • Often combined with CO₂ (273 K) data for ultramicropores
DFT Kernel Selection Guide
Choosing and reporting a DFT kernel
- Match the model to the system: adsorbate, temperature, wall material (silica, carbon, zeolite) and pore shape (slit, cylinder, sphere) must all match the kernel.
- Check the fit: the fitted isotherm should overlay the measured one; a poor fit means the kernel is wrong for the sample.
- Pick the right branch: equilibrium kernels need the desorption branch; metastable-adsorption kernels use the adsorption branch.
- Report it: name the kernel, the branch, and the regularisation used, so others can reproduce the distribution.
Method Optimization
Adsorbate Selection: N₂ vs. Ar vs. Kr
Choosing the correct probe molecule is critical for accurate surface area and pore size determination. Each adsorbate offers unique advantages depending on material properties and analysis goals.
Selection Guidelines
Nitrogen (77 K) - Standard
✓ Most common, extensive literature
✓ Good for mesopores and general use
✗ Quadrupole moment causes orientation
✗ Slow diffusion in narrow micropores
✗ May miss ultramicropores < 0.7 nm
Argon (87 K) - Recommended
✓ No quadrupole, spherical molecule
✓ Better for micropore analysis
✓ Fills micropores at higher P/P₀
✓ Cleaner isotherms for zeolites
✗ Less common, fewer references
Krypton (77 K) - Low SA
✓ Low P₀ gives high sensitivity
✓ Ideal for SA < 1 m²/g
✓ Reduced thermal transpiration
✗ Expensive gas
✗ Limited pore size analysis
CO₂ (273 K) - Ultramicropores
✓ Room temperature analysis
✓ Fast diffusion in < 1 nm pores
✓ No cryogenics needed
✗ Limited to < 1.5 nm pores
✗ Cannot determine mesopores
Quality Assurance
BET Validation & C-Constant Interpretation
Proper BET analysis requires careful validation to ensure physically meaningful results. ISO 9277:2022 provides strict criteria for determining the validity of BET surface area measurements.
Rouquerol Criteria for BET Range Selection
- C must be positive (a negative intercept means the range is wrong)
- V(1 - P/P₀) must increase continuously with P/P₀
- The p/p₀ corresponding to nm must lie inside the selected range
- 1/(√C + 1) should agree with that p/p₀ to within about 10 %
- Also check: linear fit with a correlation coefficient ≥ 0.999 (ISO 9277)
C-Constant Interpretation
Common BET Errors
⚠️ Common Pitfalls to Avoid:
- • Using default 0.05-0.30 range without validation
- • Ignoring negative intercept (indicates micropores)
- • Applying BET to Type I isotherms above P/P₀ = 0.1
- • Not checking if n_m falls within selected range
- • Using wrong cross-sectional area for adsorbate
- • Insufficient equilibration time during measurement
Methodology
Sample Preparation & Degassing
Proper sample preparation is critical for reproducible gas adsorption measurements. Complete removal of physisorbed species without altering the pore structure requires careful control of temperature, vacuum, and time. The cross-method sample preparation and degassing guide covers the same constraints across MIP and CFP and includes a generic protocol with a common-mistakes checklist.
Degassing Protocols by Material Type
Sample Mass Guidelines
Optimal sample mass depends on expected surface area:
- • > 1000 m²/g: 30-50 mg
- • 100-1000 m²/g: 50-100 mg
- • 10-100 m²/g: 100-300 mg
- • 1-10 m²/g: 300-500 mg
- • < 1 m²/g: 1-2 g (use Kr)
Target: Total surface area in tube > 10 m² for accurate measurement
Industry Applications
Gas Adsorption Applications Across Industries
BET surface area and pore size distribution are critical quality parameters across diverse industries, from catalyst development to pharmaceutical formulation.
Battery Materials
Electrode surface area affects rate capability. Typical values: LFP 10-20 m²/g, graphite 2-5 m²/g. Micropore analysis critical for SEI formation studies.
MOFs & Zeolites
Record surface areas >7000 m²/g achieved. Ar@87K preferred for accurate micropore analysis. QSDFT essential for heterogeneous pore networks.
Carbon Capture
CO₂ capacity correlates with micropore volume <0.7 nm. BET determines optimal activation conditions. Isosteric heat from multi-temperature isotherms.
Pharmaceuticals
API surface area affects dissolution rate (Noyes-Whitney). Excipient porosity controls tablet disintegration. BET validates milling and granulation.
Activated Carbons
Surface areas 500-3000 m²/g typical. Combined N₂ and CO₂ analysis for full pore spectrum. QSDFT-slit model for accurate PSD.
Nanomaterials
Particle size calculation from BET: d = 6000/(ρ × S_BET). Aggregation assessment through t-plot external surface area.
Cement & Construction
N₂ surface areas of hydrated cement paste depend strongly on how the sample was dried, so drying protocols must be reported. N₂ adsorption complements MIP for complete pore structure.
Next Steps
Continue with gas adsorption
The multi-point BET calculator turns isotherm points into a specific surface area with Rouquerol checks; the BJH guide, t-plot guide and 4V/A guide walk through the other standard results. The nitrogen-versus-argon comparison covers when to choose each adsorbate. For materials whose performance depends on the macropore range, mercury intrusion porosimetry is the appropriate companion technique.