Catalytic resonance theory

Catalytic resonance theory is constructed on the Sabatier principle of catalysis developed by French chemistry Paul Sabatier. In the limit of maximum catalytic performance, the surface of a catalyst is neither too strong nor too weak. Strong binding results in an overall catalytic reaction rate limitation due to product desorption, while weak binding catalysts are limited in the rate of surface chemistry. Optimal catalyst performance is depicted as a 'volcano' peak using a descriptor of the chemical reaction defining different catalytic materials. Experimental evidence of the Sabatier principle was first demonstrated by Balandin in 1960.[1][2]

The concept of catalytic resonance was proposed on dynamic interpretation of the Sabatier volcano reaction plot.[3] As described, extension of either side of the volcano plot above the peak defines the timescales of the two rate-limiting phenomena such as surface reaction(s) or desorption.[4] For binding energy oscillation amplitudes that extend across the volcano peak, the amplitude endpoints intersect the transiently accessible faster timescales of independent reaction phenomena. At the conditions of sufficiently fast binding energy oscillation, the transient binding energy variation frequency matches the natural frequencies of the reaction and the rate of overall reaction achieves turnover frequencies greatly in excess of the volcano plot peak.[5]

The basis of catalytic resonance theory utilizes the transient behavior of adsorption, surface reactions, and desorption as surface binding energy and surface transition states oscillate with time. The binding energy of a single species, i, is described via a temporal functional including square or sinusoidal waves of frequency, f_i, and amplitude, dU_i:

${\displaystyle \Delta H_{\text{i,ads}}(t)=\Delta H_{\text{0}}+\Delta U\sin(f_{i}t),}$[6]

Other surface chemical species, j, are related to the oscillating species, i, by the constant linear parameter, gamma γ_i-j:

${\displaystyle \gamma _{\text{i-j}}={\frac {\Delta H_{\text{i,ads}}(t)}{\Delta H_{\text{j,ads}}(t)}}}$[7]

The two surface species also share the common enthalpy of adsorption, delta δ_i-j. Specification of the oscillation frequency and amplitude of species i and relating γ_i-j and δ_i-j for all other surface species j permits determination of all chemical surface species adsorption enthalpy with time. The transition state energy of a surface reaction between any two species i and j is predicted by the linear scaling relationship of the Bell–Evans–Polanyi principle which relates to the surface reaction enthalpy, ΔH_i-j, to the transition state energy, E_a, by parameters α and β with the following relationship:

${\displaystyle E_{a},_{i-j}(t)=\beta +\alpha \Delta H_{\text{i,j}}(t)}$

The oscillating surface and transition state energies of chemical species alter the kinetic rate constants associated with surface reaction, adsorption, and desorption. The surface reaction rate constant of species i converting to surface species j includes the dynamic activation energy:

${\displaystyle k_{i-j}(t,T)=Ae^{\frac {-E_{\rm {a,i-j}}(t)}{RT}},}$

The resulting surface chemistry kinetics are then described via a surface reaction rate expression containing dynamic kinetic parameters responding to the oscillation in surface binding energy:

${\displaystyle {\frac {d\theta _{i}(t)}{dt}}=k_{ads,i}P_{i}\theta _{*}(t)-k_{des,i}(t)\theta _{i}(t)+\sum _{k=1}^{N}{\nu _{i,k}r_{k}(t,T)}}$,

with k reactions with dynamic activation energy. The desorption rate constant also varies with oscillating surface binding energy by:

${\displaystyle k_{des,i}(t)=A_{des,i}e^{\left\lbrack {\frac {\Delta H_{ads,i}(t)}{RT}}\right\rbrack }}$.

Implementation of dynamic surface binding energy of a reversible A-to-B reaction on a heterogeneous catalyst in a continuous flow stirred tank reactor operating at 1% conversion of A produces a sinusoidal binding energy in species B as shown.[8] In the transition between surface binding energy amplitude endpoints, the instantaneous reaction rate (i.e., turnover frequency) oscillates over an order of magnitude as a limit cycle solution.

Oscillating binding energies of all surface chemical species introduces periodic instances of transient behavior to the catalytic surface. For slow oscillation frequencies, the transient period is only a small quantity of the oscillation time scale, and the surface reaction achieves a new steady state. However, as the oscillation frequency increases, the surface transient period approaches the timescale of the oscillation and the catalytic surface remains in a constant transient condition. A plot of the averaged turnover frequency of a reaction with respect to applied oscillation frequency identifies the 'resonant' frequency range for which the transient conditions of the catalyst surface match the applied frequencies.[9] The 'resonance band' exists above the Sabatier volcano plot maximum of a static system with average reaction rates as high as five orders of magnitude faster than that achievable by conventional catalysis.

Surface binding energy oscillation also occurs to different extent with the various chemical surface species as defined by the γ_i-j parameter. For any non-unity γ_i-j system, the asymmetry in the surface energy profile results in conducting work to bias the reaction to a steady state away from equilibrium.[10] Similar to the controlled directionality of molecular machines, the resulting ratchet (device) energy mechanism selectively moves molecules through a catalytic reaction against a free energy gradient.[11]

Application of dynamic binding energy to a surface with multiple catalytic reactions exhibits complex behavior derived from the differences in the natural frequencies of each chemistry; these frequencies are identified by the inverse of the adsorption, desorption, and surface kinetic rate parameters. Considering a system of two parallel elementary reactions of A-to-B and A-to-C that only occur on a surface, the performance of the catalyst under dynamic conditions will result in varying capability for selecting either reaction product (B or C).[12] For the depicted system, both reactions have the same overall thermodynamics and will produce B and C in equal amounts (50% selectivity) at chemical equilibrium. Under normal static catalyst operation, only product B can be produced at selectivities greater than 50% and product C is never favored. However, as shown, the application of surface binding dynamics in the form of a square wave at varying frequency and fixed oscillation amplitude but varying endpoints exhibits the full range of possible reactant selectivity. In the range of 1-10 Hertz, there exists a small island of parameters for which product C is highly selective; this condition is only accessible via dynamics. [13]

Catalytic rate enhancement via dynamic perturbation of surface active sites has been demonstrated experimentally with dynamic electrocatalysis and dynamic photocatalysis. Those results may be explained in the framework of catalytic resonance theory but conclusive evidence is still lacking:

• In 1978, the electro-oxidation of formic acid on a platinum electrode was studied under the application of constant potentials and square-wave pulsed potentials. The latter was found to enhance the current density (and thus catalytic activity) by up to 20 times compared to the potentiostatic conditions, with the optimal wave amplitude and frequency of 600 mV and 2000 Hz, respectively.[14] In 1988, the oxidation of methanol on a platinum electrode was conducted under pulsed potentials between 0.4 and 1.18 V, resulting in an average current almost 100 times higher than the steady-state current at 0.4 V.[15]

• Using the formic acid electro-oxidation reaction, oscillation of the applied electrodynamic potential between 0 and 0.8 volts accelerated the formation rate of carbon dioxide more than an order of magnitude higher (20X) than what was achievable on platinum, the best existing catalyst.[16] The maximum catalytic rate was experimentally observed at a frequency of 100 Hz; slower catalytic rates were observed at higher and lower electrodynamic frequencies. The resonant frequency was interpreted as the oscillation between conditions favorable to formic acid decomposition (0 V) and conditions favorable to form CO₂ (0.8 V).[17]
• The concept of implementing periodic illumination to improve the quantum yield of a typical photocatalytic reaction was first introduced in 1964 by Miller et al. In this work, they showed enhanced photosynthetic efficiency in the conversion of CO₂ to O₂ when the algal culture was exposed to periodic illumination in a Taylor vortex reactor.[18] Sczechowski et al. later implemented the same approach for heterogeneous photocatalysis in 1993, where they demonstrated 5-fold increment in photoefficiency of formate decomposition by cycling between light and dark conditions with periods of 72 ms and 1.45 s respectively.[19] They hypothesized that upon illumination of the catalyst, there is a critical illumination time during which absorbed photons generate oxidizing species (h_vb⁺) on the surface of the catalyst. The generated species or their intermediates go on to react with substrates on the surface or in the bulk. During dark period, adsorption, desorption, and diffusion generally occurs in the absence of photons. After a critical recovery period in the dark, the photocatalyst can efficiently use photons again when photons are reintroduced. A summary of work involving “dynamic” photocatalysis was provided by Tokode et al. in 2016.[20]

• Dynamic promotion of methanol decomposition was demonstrated on 2 nm Pt nanoparticles using pulsed light.[21] The rate acceleration to form H₂ relative to static illumination was attributed to the selective weakening of adsorbed carbon monoxide, thereby also increasing the quantum efficiency of applied light.

• In 2021, Sordello et al.[22] experimentally demonstrated a 50% increase of the quantum yield for the Hydrogen Evolution Reaction (HER) over Pt/TiO₂ nanoparticles via formic acid photoreforming under Controlled Period Illumination (CPI).

Implementation of catalyst dynamics has been proposed to occur by additional methods using oscillating light, electric potential, and physical perturbation.[23]