Heterogeneous Catalysis Simulator: Sabatier Principle & Volcano Plot

simulator advanced ~14 min
Loading simulation...
TOF = 0.68 s⁻¹ — near optimal catalytic activity

At E_b = 1.5 eV and 600 K with BEP coefficient α = 0.5, the turnover frequency is 0.68 s⁻¹ — near the peak of the volcano plot where adsorption and desorption rates are balanced.

Formula

TOF = A × exp(-E_a/RT) × θ × (1 - θ)
E_a = E₀ - α × E_b (Bronsted-Evans-Polanyi)
θ = K₀ × exp(E_b/RT) × P / (1 + K₀ × exp(E_b/RT) × P)

Catalysis at Surfaces

Over 90% of industrial chemical processes rely on heterogeneous catalysis — reactions accelerated by solid surfaces. From ammonia synthesis feeding billions to catalytic converters cleaning exhaust, these processes depend on the precise interaction between reactant molecules and catalyst surfaces. The key insight, formulated by Paul Sabatier in 1911, is that the best catalyst provides an intermediate level of interaction: strong enough to break bonds in reactants, weak enough to release the products.

The Volcano Plot

When catalytic activity is plotted against the binding energy of a key intermediate, a volcano-shaped curve emerges. On the ascending (left) side, stronger binding increases surface coverage and lowers activation barriers. On the descending (right) side, binding is so strong that products cannot desorb, poisoning the surface. The peak represents the Sabatier optimum — and remarkably, the same metals (Pt, Pd, Ru) sit near the peak for many different reactions.

The BEP Relation

The Bronsted-Evans-Polanyi (BEP) relation provides the theoretical backbone of volcano plots. It states that activation energy decreases linearly with increasing binding energy: stronger binding makes it easier to break reactant bonds. Combined with the Langmuir adsorption model, this creates the fundamental tradeoff: increasing binding energy simultaneously lowers the activation barrier (faster reaction) but increases surface coverage (fewer available sites), producing the volcano shape.

Designing Better Catalysts

Modern computational catalysis uses density functional theory (DFT) to calculate binding energies and screen thousands of candidate materials. Descriptor-based design identifies the optimal binding energy from volcano plots, then searches for alloys and nanostructures that achieve it. This approach has led to discoveries of non-precious metal catalysts for fuel cells, electrochemical CO₂ reduction, and nitrogen fixation — potentially replacing expensive platinum group metals with abundant alternatives.

FAQ

What is the Sabatier principle?

The Sabatier principle states that the optimal catalyst binds reactants with intermediate strength — strong enough to activate them but weak enough to release the products. This creates a volcano-shaped plot of catalytic activity vs. binding energy, with the best catalysts (Pt, Pd, Rh) sitting at the peak.

What is a volcano plot in catalysis?

A volcano plot graphs catalytic activity (turnover frequency) against a descriptor like adsorption energy. Metals on the left side bind too weakly (Au, Ag); those on the right bind too strongly (W, Mo). The peak metals (Pt, Pd, Ru) achieve optimal balance, explaining why they are the most effective catalysts.

What is the BEP relation?

The Bronsted-Evans-Polanyi (BEP) relation states that activation energy correlates linearly with reaction energy: E_a = E₀ - α × ΔE. In catalysis, stronger binding lowers the activation barrier for adsorption but raises it for desorption, creating the fundamental tradeoff captured by the volcano plot.

Why are platinum group metals the best catalysts?

Platinum group metals (Pt, Pd, Rh, Ru, Ir) sit near the peak of volcano plots for many reactions because their d-band electronic structure provides intermediate binding strength with common reactants (H₂, O₂, CO, NO). Their d-band center is neither too high (over-binding) nor too low (under-binding).

Sources

Embed

<iframe src="https://homo-deus.com/lab/surface-chemistry/heterogeneous-catalysis/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub