Conspectus
Microkinetic modeling based on density functional
theory (DFT)
energies plays an essential role in heterogeneous catalysis because
it reveals the fundamental chemistry for catalytic reactions and bridges
the microscopic understanding from theoretical calculations and experimental
observations. Microkinetic modeling requires building a set of ordinary
differential equations (ODEs) based on the calculation results of
thermodynamic properties of adsorbates and kinetic parameters for
the reaction elementary steps. Solving a microkinetic model can extract
information on catalytic chemistry, including critical reaction intermediates,
reaction pathways, the surface species distribution, activity, and
selectivity, thus providing vital guidelines for altering catalysts.
However, the quantitative reliability of traditional microkinetic
models is often insufficient to conclusively extrapolate the mechanistic
details of complex reaction systems. This can be attributed to several
factors, the most important of which is the limitation of obtaining
an accurate estimation of the energy inputs via traditional calculation
methods. These limitations include the difficulty of using static
DFT methods to calculate reaction energies of adsorption/desorption
processes, often rate-controlling or selectivity-determining steps,
and the inadequate consideration of surface coverage effects. In addition,
the robust microkinetic software is rare, which also complicates the
resolution of complex catalytic systems.
In this Account, we
review our recent works toward refining the
predictions of microkinetic modeling in heterogeneous catalysis and
achieving theory–experiment parity for activity and selectivity.
First, we introduce CATKINAS, a microkinetic software developed in
our group, and show how it disentangles the problem that traditional
microkinetic software has and how it can now be applied to obtain
kinetic results for more sophisticated reaction systems. Second, we
describe a molecular dynamics method developed recently to obtain
the free-energy changes for the adsorption/desorption process to fill
in the missing energy inputs. Third, we show that a rigorous consideration
of surface coverage effects is pivotal for building more realistic
models and obtaining accurate kinetic results. Following a series
of studies on acetylene hydrogenation reactions on Pd catalysts, we
demonstrate how this new approach can provide an improved quantitative
understanding of the mechanism, active site, and intrinsic structural
sensitivity. Finally, we conclude with a brief outlook and the remaining
challenges in this field.