Next Article in Journal
Synthesis of 2-Benzylidene-3-Pyrrolines and Their Synthetic Transformation
Next Article in Special Issue
Heterogeneously Catalyzed γ-Valerolactone Hydrogenation into 1,4-Pentanediol in Milder Reaction Conditions
Previous Article in Journal / Special Issue
An Iron(III) Complex with Pincer Ligand—Catalytic Water Oxidation through Controllable Ligand Exchange
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Reactions
Volume 1
Issue 1
10.3390/reactions1010004
reactions-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Requiem for the Rate-Determining Step in Complex Heterogeneous Catalytic Reactions?

by
Dmitry Yu. Murzin
Process Chemistry Centre, Åbo Akademi University, 20500 Turku/Åbo, Finland
Reactions 2020, 1(1), 37-46; https://0-doi-org.brum.beds.ac.uk/10.3390/reactions1010004
Submission received: 1 September 2020 / Revised: 8 September 2020 / Accepted: 9 September 2020 / Published: 10 September 2020
(This article belongs to the Special Issue Mechanism and Kinetics of Complex Reactions— 2020 Selected Papers from Reactions’ Editorial Board Members)
Browse Figures
Versions Notes

Abstract

:
The concept of the rate determining step, i.e., the step having the strongest influence on the reaction rate or even being the only one present in the rate equation, is often used in heterogeneous catalytic reactions. The utilization of this concept mainly stems from a need to reduce complexity in deriving explicit rate equations or searching for a better catalyst based on the theoretical insight. When the aim is to derive a rate equation with eventual kinetic modelling for single-route mechanisms with linear sequences, the analytical rate expressions can be obtained based on the theory of complex reactions. For such mechanisms, a single rate limiting step might not be present at all and the common practice of introducing such steps is due mainly to the convenience of using simpler expressions. For mechanisms with a combination of linear and nonlinear steps or those just comprising non-linear steps, the reaction rates are influenced by several steps depending on reaction conditions, thus a reduction in complexity to a single rate limiting step can lead to misinterpretations. More widespread utilization of a microkinetic approach when the reaction rate constants can be computed with reasonable accuracy based on the theoretical insight, and availability of software for kinetic modelling, when a system of differential equations for reactants and products will be solved together with differential equations for catalytic species and the algebraic conservation equation for the latter, will eventually make the concept of the rate limiting step obsolete.

1. Introduction

The kinetics of heterogeneous catalytic reactions has been the focus of numerous studies because of its theoretical and practical importance [1,2,3,4,5,6,7,8,9,10]. Kinetic analysis was applied for decades to assist reactor design and scaling up of various reactions relevant for oil refining and synthesis of basis and specialty chemicals. More recently, because of a widespread application of so-called flow chemistry to streamline pharmaceutical research and products, kinetic analysis also started to be more widely utilized in the field of fine chemicals [11,12,13]. It was argued [11] that the framework of a mechanistically based rate equation combined with in situ experimental data is useful for the interpretation of complex reaction networks. Increasing complexity of the catalytic reactions, either relevant for oil refining or pharmaceutical synthesis, poses some restrictions regarding incorporation of such complex mechanisms into tractable rate equations.
Thus, another trend is to simplify complex reaction networks and assume that there is a limiting step, the rate of which is a good approximation of the reaction rate of the whole network [14,15]. This approach of a rate-determining step (RDS) became very popular in chemical reaction engineering, allowing the derivation of kinetic equations for reaction mechanisms of differing complexity. Such an assumption of a rate-determining step makes the rate expressions more tractable, facilitating the kinetic modelling of mainly steady-state kinetics even at the expense of losing important mechanistic information. Moreover, in several instances, without a substantial reduction in complexity, kinetic modelling can be computationally very demanding, especially when coupled with transport phenomena [16].
On the other hand, the concept of RDS was suggested to give fruitful ideas for modifying the reaction conditions or design the catalyst. Subsequently, there have been many theoretical papers addressing how to determine the step which is the rate determining [17,18]. The most popular is probably the degree of rate control proposed by Campbell, based on the general principles of differential sensitivity analysis to determine sensitivity of the output to one input parameter [19,20,21,22,23,24,25]. It was argued [16] that an approach analogous to electrical networks is even better suited for the identification of the key steps determining the reaction rate allowing the derivation of accurate rate expressions even for mechanisms with nonlinear steps. Such nonlinear steps are present not only in heterogeneous catalytic reactions, such as, for example, a water gas shift reaction [16], decomposition of ammonia [26,27] or methanol [28], but also in organometallic catalysis. In the latter case, the reaction order in catalysts deviates from unity [5]. The formation of dimeric catalytic species, on the contrary, results in the reaction order being below one, which is rather rare [29,30,31].
Contrary to the case of linear (Christiansen) sequences of a single-route reaction [5,32,33] for single-route mechanisms with nonlinear steps, a general expression for the reaction rate in steady-state and quasi-steady-state conditions cannot be derived [34]. When such derivation is possible, the rate expressions become rather complicated to handle [5,35,36,37].
To overcome difficulties in deriving explicit rate expressions, a semi-empirical kinetic polynomial was proposed [38] as a substitute of analytical rate equations.
However, in catalysis by solids or organometallic complexes, quite often either there are nonlinear steps involved or the rate equations for even single-route mechanisms and obviously multi route reactions do not have rate limiting steps [34].
It was recognized a long time ago that for multi-route (or composite) reactions, the concept of the rate-controlling steps is not required [34,39], while in other cases it should be considered with care, as misinterpretation is easily possible.
In the current work, several cases of complex single- and multi-route heterogeneous catalytic reactions will be presented, highlighting an apparent danger of falling into a trap of using the rate determining step, which leads to confusion rather than clarification of reaction mechanisms.

2. Three Step Sequence

Let us consider first a three-step linear sequence (Figure 1)
The general equation for this three-step sequence in the case of all reversible steps can be written in a form of frequencies of steps [6]
r = ω + 1 ω + 2 ω + 3 − ω − 1 ω − 2 ω − 3 ω + 2 ω + 3 + ω − 3 ω + 2 + ω − 3 ω − 2 + ω + 3 ω + 1   +   ω − 1 ω + 3 + ω − 1 ω − 3 + ω + 1 ω + 2 +   ω − 2 ω + 1 + ω − 2 ω − 1 C c a t
where r is the reaction rate, ω + 1 is the frequency of step 1, etc., and Ccat is the catalyst concentration. The frequencies of steps are obtained by dividing the rate expressions of a particular step by the concentration of catalytic species in that expression. For example, in the case of the reaction mechanism presented in Figure 2, the general form of the Equation (1) gives [6]
r = ( k + 1 C H 2 k + 2 C C 2 H 4 k + 3 − k − 1 k − 2 k − 3 C C 2 H 6 ) C c a t k + 2 k + 3 + k − 3 C C 2 H 6 k + 2 + k − 3 C C 2 H 6 k − 2 + k + 3 k + 1   C H 2 +   k − 1 k + 3 + k − 1 k − 3 C C 2 H 6 + k + 1 C H 2 k + 2 C C 2 H 4 +   k − 2 k + 1 C H 2 + k − 2 k − 1
Thus, the frequency of the first step in the forward direction is k + 1 C H 2 , while for the reverse step it would be just k − 1 . When all steps are irreversible, Equation (2) can be easily simplified, giving
r = ω + 1 ω + 2 ω + 3 ω + 2 ω + 3 + ω + 3 ω + 1   + ω + 1 ω + 2 C c a t = k + 1 C H 2 k + 2 C C 2 H 4 k + 3 C c a t k + 2 k + 3 + k + 3 k + 1   C H 2 + k + 1 C H 2 k + 2 C C 2 H 4 = k + 1 C H 2 C C 2 H 4 C c a t 1 + k + 1   k + 2 C H 2 + k + 1 k + 3 C H 2 C C 2 H 4
It follows from Equation (3) that depending on the values of partial pressures of, for example, hydrogen, the reaction order can change from unity at low pressures, to zero at high hydrogen pressures. The relative contribution of the last term in the denominator of Equation (3) k + 1 C H 2 C C 2 H 4 / k + 3 in comparison with two other terms defines the reaction order towards ethylene.
This example illustrates an apparent danger of defining the rate and even more turnover frequency through just one rate limiting step without taking into account kinetic regularities. To highlight the need to consider all steps in a linear catalytic sequence even further, let us consider the apparent activation energy for the three-step mechanism (Equation (3)). In what follows, the term “apparent” corresponds to the observed activation energy of a complex multistep reaction in the kinetic regime, which should not be confused with the case when the observed kinetics is influenced by the transport phenomena.
A more simple case of the two-step sequence was addressed in the literature [41] using the definition of the apparent activation energy
E a , app = − R ∂ ln r ∂ ( 1 / T )
and giving
E a , app = ω + 2 ω + 1   + ω + 2 E a , 1 + ω + 1 ω + 1   + ω + 2 E a , 2
As follows from Equation (5), the apparent activation energy depends not only on the step with the highest activation energy, but also on the contribution of the other step with a lower activation energy. For the three-step sequence using Equations (3) and (4), one gets
E a , app = R T 2 ∂ ∂ T ( ln k + 1 + ln C H 2 + ln k + 2 + ln C C 2 H 4 + ln k + 3 + ln C c a t − ln ( k + 2 k + 3 + k + 3 k + 1   C H 2 + k + 1 C H 2 k + 2 C C 2 H 4 ) )
Some terms in Equation (6) can be easily differentiated
∂ ∂ T ( ln k 0 + 1 exp ( − E + 1 / R T ) = ∂ ∂ T ( ln k 0 + 1 − E + 1 / R T ) = E + 1 R T 2
Taking Equation (7) into account and also using a standard approach for the differentiation of a complex function
∂ ln U ∂ T = ∂ ln U ∂ U ∂ U ∂ T = 1 T ∂ U ∂ T
one gets the final expression for the activation energy of a three step sequence
E a , app = ω + 2 ω + 3 E a , 1 + ω + 1 ω + 3 E a , 2 + ω + 1 ω + 2 E a , 3 ω + 2 ω + 3 + ω + 3 ω + 1   + ω + 1 ω + 2
It can be shown that in only exceptional cases the apparent activation energy of the overall reaction is determined by the activation energy of only one rate-limiting step in a steady-state sequence.

3. Four Step Sequence

As an example of such a sequence, the liquid-phase hydrogenation of aromatic compounds will be considered. The reaction kinetics was extensively investigated previously in [42,43,44,45,46]. The following generic scheme for hydrogenation of an aromatic compound A to the product B was proposed
1.*A + H2 ⇔ *AH2
2.*AH2 + H2 ⇔ *AH4
3.*AH4 ⇒ *Y
3′.*Y + H2 ⇔ *B (fast)
4.*B + A Ξ *A + B
A + 3H2 = B
where Y is cyclohexene or its derivative, AH2 and AH4 are intermediate complexes, step 3′ is fast because hydrogenation of cycloalkenes is much faster than aromatic compounds and step 4 is at quasi-equilibrium.
A general equation for the four step sequence is [6]
r = ω + 1 ω + 2 ω + 3 ω + 4 − ω − 1 ω − 2 ω − 3 ω − 4 D C c a t
where
D = ω + 2 ω + 3 ω + 4 +   ω − 1 ω + 3 ω + 4 + ω − 1 ω − 2 ω + 4 + ω − 1 ω − 2 ω − 3 + + ω + 1 ω + 2 ω + 3 +   ω − 4 ω + 2 ω + 3 + ω − 4 ω − 1 ω + 3 + ω − 4 ω − 1 ω − 2 + + ω + 4 ω + 1 ω + 2 +   ω − 3 ω + 1 ω + 2 + ω − 3 ω − 4 ω + 2 + ω − 3 ω − 4 ω − 1 + + ω + 3 ω + 4 ω + 1 +   ω − 2 ω + 4 ω + 1 + ω − 2 ω − 3 ω + 1 + ω − 2 ω − 3 ω − 4
Because of irreversibility of step (3), Equations (11) and (12) can be simplified
r = ω + 1 ω + 2 ω + 3 ω + 4 D C c a t
D = ω + 2 ω + 3 ω + 4 +   ω − 1 ω + 3 ω + 4 + ω − 1 ω − 2 ω + 4 + ω + 1 ω + 2 ω + 3 +   ω − 4 ω + 2 ω + 3 + ω − 4 ω − 1 ω + 3 + ω − 4 ω − 1 ω − 2 + + ω + 4 ω + 1 ω + 2 +   ω + 3 ω + 4 ω + 1 +   ω − 2 ω + 4 ω + 1
Step 4 is at quasi-equilibrium, and therefore the rates are fast in both directions. Further simplifications of Equation (13) are thus possible by dividing the numerator and the denominator by ω + 4 , neglecting subsequently the terms ω + 1 ω + 2 ω + 3 / ω + 4 and ω − 3 ω + 1 ω + 2 / ω + 4 . Then the following holds
r = ω + 1 ω + 2 ω + 3 ( ω + 2 ω + 3 +   ω − 1 ω + 3 + ω − 1 ω − 2 ) ( 1 + ω − 4 ω + 4 ) + ω + 1 ω + 2 + ω + 3 ω + 1 +   ω − 2 ω + 1 C c a t
Replacing frequencies in Equation (15) with the explicit expressions, one gets
r = k + 1 C H 2 k + 2 C H 2 k + 3 ( k + 2 C H 2 k + 3 +   k − 1 k + 3 + k − 1 k − 2 ) ( 1 + C B K 4 C A ) + k + 1 C H 2 k + 2 C H 2 + k + 3 k + 1 C H 2 +   k − 2 k + 1 C H 2 C c a t
where K 4 = k + 4 / k − 4 . For a binary mixture of the reactant and substrate, it is possible to use mole fractions instead of concentrations. Moreover, quite often partial pressures of hydrogen are applied instead of concentrations, resulting in another form of the rate equation
r = k ′ + 1 P H 2 k ′ + 2 P H 2 k + 3 ( k ′ + 2 P H 2 k + 3 +   k − 1 k + 3 + k − 1 k − 2 ) ( 1 + 1 − N A K 4 N A ) + k ′ + 1 P H 2 k ′ + 2 P H 2 + k + 3 k ′ + 1 P H 2 +   k − 2 k ′ + 1 P H 2 C c a t
where N A is the mole fraction of the reactant. Equation (16) contains the modified rate constants k ′ + 1 ; k ′ + 2 as instead of hydrogen concentration in the liquid phase, the pressure of hydrogen is used.
Equation (17) presented in a slightly different form
r = k ′ + 1 k ′ + 2 k + 3 k + 3 k ′ + 1 + k − 2 k ′ + 1 + k ′ + 2 k + 3 P H 2 N A C c a t k − 1 k + 3 + k − 1 k − 2 k + 3 k ′ + 1 + k − 2 k ′ + 1 + k ′ + 2 k + 3 N A P − 1 H 2 + (   k − 1 k + 3 + k − 1 k − 2 ) ( 1 − N A ) ( k + 3 k ′ + 1 + k − 2 k ′ + 1 + k ′ + 2 k + 3 ) K 4 ) P − 1 H 2 + k ′ + 2 k + 3 ( 1 − N A ) ( k + 3 k ′ + 1 + k − 2 k ′ + 1 + k ′ + 2 k + 3 ) K 4 + k ′ + 1 k ′ + 2 k + 3 k ′ + 1 + k − 2 k ′ + 1 + k ′ + 2 k + 3 P H 2 N A + N A
was derived in [43] using a more tedious procedure, namely applying directly the steady-state conditions
r + 1 − r − 1 = r + 2 − r − 2 = r + 3
as well as the quasi-equilibrium for step 4 and finally the balance equation (i.e., the sum of all coverages is equal to unity).
All steps in Equation (10) are needed to explain the observed kinetic regularities, because it was demonstrated that at low hydrogen partial pressures, the reaction order in hydrogen can exceed unity, which is possible when k−1 ≠ 0, while at high pressures the reaction kinetics obeys frequently observed zero orders in the substrate and hydrogen. In the latter case, the overall rate is determined by the isomerization of the adsorbed AH4 species. The surface of the catalyst is completely covered with this complex and the reaction rate ceases to depend on either hydrogen pressure or the concentration of the substrate.
Apparently, the concept of a single rate-determining step is not able to account for a rich kinetic behavior, observed experimentally in hydrogenation of aromatics.

4. Oxidative Dehydrogenation of Ethanol: A Multi-Route Mechanism with Nonlinear Steps

The mechanism and kinetics of oxidative dehydrogenation of ethanol over gold catalysts was recently reported [47]. DFT calculations indicated that the activation of molecular oxygen is facilitated in the presence of ethanol acting as a hydrogen donor. As a result of hydrogen abstraction from the donor leading to the formation of an OOH intermediate, the dioxygen undergoes dissociation on the catalyst surface. The scheme of the reaction mechanism is given below
N ( 1 a ) N ( 1 b ) N ( 1 c ) 1 .   O 2 + * Ξ O 2 * 1 1 1 2 .   EtOH + * Ξ EtOH * 2 2 2 3 .   EtOH *   +   O 2 * ↔ EtO * + OOH * 2 1 2 4 .   OOH * + * → O * + OH * 2 1 1 5 .   O * + EtO * → AcH + * + OH * 2 2 2 6 .   2 OH * ↔ O * + H 2 O + * 2 2 2 7 a .   2 O * → O 2 * + * 1 0 0 7 b .   EtOH * + O * → EtO * + OH * 0 1 0 7 c .   O * + OOH * → O 2 * + OH * 0 0 1 N ( 1 ) 2 EtOH + O 2 → 2 AcH + 2 H 2 O
In (20), on the right hand side of the equations for the steps (i.e., 1 to 7c), the stoichiometric numbers (i.e., 0, 1, 2) along alternative routes N(1a) - N(1c) are given [32]. These numbers can be equal to zero, if a particular step is not involved in the reaction mechanism. In Equation (20), EtOH stands for ethanol, while AcH represents acetaldehyde and * denotes a surface site.
According to the DFT calculations [47] illustrated in Figure 3, the largest activation barrier corresponds to step 3 in Equation (20), and thus superficially it can be considered in the conventional approaches as the rate limiting step. However, the activation energy of the same step in the reverse direction is very low and step 3 can even be viewed to be close to the quasi-equilibrium.
From the considerations above on the apparent activation energy of the three-step sequence with linear steps, it is clear that the apparent activation energy and the reaction rate depend on the contribution of all steps in the mechanism, which could be treated using the steady-state approximation. This approach is possible because of the specific features of the reaction mechanism comprising only nonlinear steps and thus allowing derivation of an analytic rate expression. Such derivation was done in [47] based on the DFT computed reaction energies assuming quasi-equilibria for steps 1 and 2, irreversibility for steps 4 and 5 and reversibility for steps 3 and 6.
From the viewpoint of the global reaction kinetics, not only steps 3–6 are important but also the final step in the mechanism, step 7, reflecting the fate of atomically adsorbed oxygen. Therefore, three options are possibly influencing the stoichiometric numbers of other steps. Step 7a in the route N(1a) corresponds to recombination of atomic oxygen, while such recombination in step 7c (route N(1c)) involves a reaction with a peroxy species OOH. In the route N(1b), atomically adsorbed oxygen assists in the abstraction of hydrogen from the substrate.
Derivation of the rate expression (per mole of oxygen) for the route N(1b) was described in [47], giving
r N ( 1 b ) = k + 3 k + 4 K 1 K 2 C O 2 C E t O H ( k + 4 + k − 3 k + 7 b k + 5 C E t O H ) ( 1 + K 2 C E t O H + K 1 C O 2 + 2 k + 7 b K 2 C E t O H k + 5 + k + 3 k + 4 K 1 C O 2 k + 7 b ( k + 4 + k − 3 k + 7 b k + 5 C E t O H ) + 2 k + 6 ( k + 3 k + 4 K 1 K 2 C O 2 C E t O H k + 4 + k − 3 k + 7 b k + 5 C E t O H + k − 6 k + 3 k + 4 K 1 C O 2 C H 2 O k + 7 b ( k + 4 + k − 3 k + 7 b k + 5 C E t O H ) ) + k + 3 K 1 K 2 C O 2 C E t O H k + 4 + k − 3 k + 7 b k + 5 C E t O H ) 2
Just for illustration purposes, let us consider a case when step 3 is at quasi-equilibrium and step 6 is irreversible. Equation (21) is then transformed to
r N ( 1 b ) = K 3 k + 4 k + 5 K 1 K 2 C O 2 C E t O H / k 7 b ( 1 + K 2 C E t O H + K 1 C O 2 + 2 k + 7 b K 2 C E t O H k + 5 + K 3 k + 4 k + 5 K 1 C O 2 k 2 + 7 b C E t O H + 2 k + 6 ( k + 4 k + 5 K 1 K 2 K 3 C O 2 k + 7 b ) + K 3 K 1 K 2 k + 5 C O 2 C E t O H k + 7 b ) 2
For two other options also considering the quasi-equilibria of the step 3, the rates are for the routes N(1b)
r N ( 1 c ) = k + 4 k + 5 / 2 k + 7 c 2 k + 7 c K 3 k + 5 K 1 K 2 C O 2 C E t O H ( 1 + K 2 C E t O H + K 1 C O 2 + k + 4 k + 7 c + ( 1 + k + 5 / 2 k + 7 c ) 2 k + 7 c K 3 k + 5 K 1 K 2 C O 2 C E t O H + k + 4 k + 5 k + 6 k + 7 c ( 2 k + 7 c K 3 k + 5 K 1 K 2 C O 2 C E t O H ) 0.25 ) 2
and N(1a), respectively
r N ( 1 c ) = k 2 + 5 / k + 7 a ( 2 k + 7 a k + 4 k 2 + 5 K 1 K 2 K 3 C O 2 C E t O H ) 2 / 3 ( 1 + K 2 C E t O H + K 1 C O 2 + ( 2 k + 7 a k + 4 k 2 + 5 K 1 K 2 K 3 C O 2 C E t O H ) 1 / 3 ( 1 + k + 5 2 k + 7 a + k + 4 k + 6 ( k 2 + 5 2 k + 7 a k + 4 ) 1 / 6 ) + ( k 2 + 5 2 k + 7 a k + 4 ) 1 / 3 ( K 1 K 2 K 3 C O 2 C E t O H ) 2 / 3 ) 2
illustrating that the rate expression and even the kinetic regularities which follow from the corresponding rate expressions are very sensitive to the chemistry of step 7 and the rate constant of this step. The mechanism N(1a) is less probable compared to other alternatives considering a high activation barrier anticipated for this step [47], while other options in principle are possible. For kinetic modelling in [47], a feasible channel of atomically adsorbed oxygen consumption with a low activation energy barrier was considered to be step 7b (Equation (20)).

5. Conclusions and Outlook

Apparently, the concept of the rate limiting will continue to be utilized in the future. This concept is useful as a part of the education of future specialists in chemical reaction engineering and catalyst development. For simple single-route reaction mechanisms, the application of the rate-determining step along with the quasi-equilibria of other steps results in easily derived rate equations. For more complex cases with several linear steps, the general approach based on the theory of complex reactions of Horiuti–Temkin allows the derivation of the corresponding rate equations. Apparent difficulties in doing this, as illustrated above for the four-step sequence, often result in a reductionistic approach assuming a rate-determining step either based on chemical intuition or a more rigorous theoretical analysis based on quantum chemical calculations.
There are cases when the rate equation for linear sequences with a small number of steps was not derived and a system of differential equations was solved even if the rate equation for the steady -state system can be readily derived. This was done, for example in [48], when a solver of differential equations was used for a four step sequence of the Fujiwara–Moritani reaction (Figure 4), unfortunately not reporting how the balance equation for catalytic species was taken into account,
For a particular case of the mechanism in Figure 4, the rate equation can be easily derived as this mechanism is a special case of Equations (11) and (12) with irreversible steps 2 and 3 simply resulting in
r = ω + 1 ω + 2 ω + 3 ω + 4 C c a t ω 2 ω 3 ω 4 +   ω − 1 ω 3 ω 4 + ω 1 ω 2 ω 3 +   ω − 4 ω 2 ω 3 + ω − 4 ω − 1 ω 3 + + ω 4 ω 1 ω 2 +   ω 3 ω 4 ω 1
Or, after introducing the explicit expressions for the frequencies of steps
r = k + 1 k + 2 k + 3 k + 4 C a c e tan i l i d e C a c r y l a t e C B Q C c a t k + 2 k + 3 k + 4 C a c r y l a t e C B Q + k − 1 k + 3 k + 4 C a c r y l a t e C B Q + k + 1 k + 2 k + 3 C a c e tan i l i d e C a c r y l a t e +   k − 4 k + 2 k + 3 C a c r y l a t e + k − 4 k − 1 k + 3 C a c r y l a t e + k + 1 k + 2 k + 4 C a c e tan i l i d e C B Q +   k + 1 k + 3 k + 4 C a c e tan i l i d e C a c r y l a t e C B Q
In a more general case, when the reaction mechanisms comprise linear and nonlinear steps without any rate-limiting ones, derivation of an explicit rate equation can be too tedious. Instead, a comparison between the experimental and calculations should be done in an implicit form using a system of differential equations considering also an algebraic balance equation for the catalytic species. Availability of computer programs with a customer written subroutine comprising a set of elementary reactions will pave the way for eventual marginalization of the concept of the rate-determining step in the kinetic modelling of complex heterogeneous catalytic reactions.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Tsukamoto, M.; Gopalaiah, K.; Kagan, H.B. Equilibrium of homochiral oligomerization of a mixture of enantiomers. Its relevance to nonlinear effects in asymmetric catalysis. Phys. Chem. 2008, 112, 15361–15368. [Google Scholar] [CrossRef]
  2. Blackmond, D.G. Kinetic profiling of catalytic organic reactions as a mechanistic tool. Am. J. Chem. Soc. 2015, 137, 10852–10866. [Google Scholar] [CrossRef]
  3. Helfferich, F.G. Kinetics of Homogeneous Multistep Reactions. In Comprehensive Chemical Kinetics; Compton, R.G., Hancock, G., Eds.; Elsevier: Amsterdam, The Netherlands, 2001; Volume 38. [Google Scholar]
  4. Leskovac, V. Comprehensive Enzyme Kinetics; Kluwer Academic/Plenum Publishers: New York, NY, USA, 2003. [Google Scholar]
  5. Temkin, O.N. Homogeneous Catalysis with Metal Complexes: Kinetic Aspects and Mechanisms; Wiley: Hoboken, NJ, USA, 2012. [Google Scholar]
  6. Murzin, D.Y.; Salmi, T. Catalytic Kinetics, Chemistry and Engineering, 2nd ed.; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
  7. Salmi, T.O.; Mikkola, J.-P.; Wärnå, J. Chemical Reaction Engineering and Reactor Technology; CRC Press: Boca Raton, FL, USA, 2010. [Google Scholar]
  8. Fogler, H.S. Elements of Chemical Reaction Engineering, 3rd ed.; Prentice Hall: Upper Saddle River, NJ, USA, 1998. [Google Scholar]
  9. Marin, G.; Yhablonsky, G.S. Kinetics of Chemical Reactions; Wiley: Hoboken, NJ, USA, 2011. [Google Scholar]
  10. Kapteijn, F.; Berger, R.J.; Moulijn, J.A. Rate Procurement and Kinetic Modelling, Handbook of Heterogeneous Catalysis; Wiley: Hoboken, NJ, USA, 2008. [Google Scholar]
  11. Blackmond, D.G. Requiem for the reaction rate equation? Catal. Lett. 2002, 83, 133–136. [Google Scholar] [CrossRef]
  12. Ball, L.T.; Corrie, T.J.A.; Cresswell, A.J.; Lloyd-Jones, G.C. Kinetic analysis of domino catalysis: A case study on gold-catalyzed arylation. ACS Catal. 2020. [Google Scholar] [CrossRef]
  13. Schmidt, O.P.; Dechert, A.-M.; Garnsey, M.; Wisniewska, H.; Blackmond, D.G. Kinetic analysis of catalytic organic reactions using a temperature scanning protocol. ChemCatChem 2019, 11, 3808–3813. [Google Scholar] [CrossRef] [Green Version]
  14. Gorban, A. Model reduction in chemical dynamics: Slow invariant manifolds, singular perturbations, thermodynamic estimates, and analysis of reaction graph. arXiv 2018, arXiv:1802.05745. [Google Scholar] [CrossRef] [Green Version]
  15. Constales, D.; Yablonsky, G.S.; D’Hooge, D.R.; Thybaut, J.W.; Marin, G.B. Physicochemical Principles of Simplifications of Complex Models, Advanced Data Analysis and Modelling in Chemical Engineering; Springer Science+Business Media: Berlin, Germany, 2017; Chapter 4; pp. 83–103. [Google Scholar]
  16. O’Malley, P.D.; Datta, R.; Vilekar, S.A. Ockham’s razor for paring microkinetic mechanisms: Electric analogy vs. Campbell’s degree of rate control. AIChE J. 2015, 61, 4332–4346. [Google Scholar] [CrossRef]
  17. Motagamwala, A.H.; Dumesic, J.A. Analysis of reaction schemes using maximum rates of constituent steps. Proc. Natl. Acad. Sci. USA 2016, 113, E2879–E2888. [Google Scholar] [CrossRef] [Green Version]
  18. Motagamwala, A.H.; Ball, M.R.; Dumesic, J.A. Microkinetic analysis and scaling relations for catalyst design. Ann. Rev. Chem. Biomol. Eng. 2018, 9, 413–450. [Google Scholar] [CrossRef]
  19. Steigelmann, C.; Andreasen, A.; Campbell, C.T. Degree of rate control: How much the energies of intermediates and transition states control rates. JACS 2009, 131, 8077–8082. [Google Scholar] [CrossRef]
  20. Campbell, C.T. The degree of rate control: A powerful tool for catalysis research. ACS Catal. 2017, 7, 2770–2779. [Google Scholar] [CrossRef] [Green Version]
  21. Wolcott, C.A.; Medford, A.J.; Studt, F.; Campbell, C.T. Degree of rate control approach to computation catalyst screening. J. Catal. 2015, 330, 197–207. [Google Scholar] [CrossRef] [Green Version]
  22. Mao, Z.; Campbell, C.T. Apparent activation energies in complex reaction mechanisms: A simple relationship via degrees of rate control. ACS Catal. 2019, 9, 9465–9473. [Google Scholar] [CrossRef]
  23. Mao, Z.; Campbell, C.T. Kinetic isotope effects: Interpretation and prediction using degrees of rate control. ACS Catal. 2020, 10, 4181–4192. [Google Scholar] [CrossRef]
  24. Mao, Z.; Campbell, C.T. The degree of rate control of catalyst-bound intermediates in catalytic reaction mechanisms: Relationships to site coverage. J. Catal. 2020, 381, 53–62. [Google Scholar] [CrossRef]
  25. Foley, B.L.; Bhan, A. Degree of rate control and De Donder relations—An interpretation based on transition state theory. J. Catal. 2020, 384, 231–251. [Google Scholar] [CrossRef]
  26. Vilekar, S.A.; Fishtik, I.; Datta, R. The steady-state kinetics of a catalytic reaction sequence. Chem. Eng. Sci. 2009, 64, 1968–1979. [Google Scholar] [CrossRef]
  27. Vilekar, S.A.; Fishtik, I.; Datta, R. The peculiar catalytic sequence of the ammonia decomposition reaction and its steady-state kinetics. Chem. Eng. Sci. 2012, 71, 333–344. [Google Scholar] [CrossRef]
  28. Vilekar, S.A.; Fishtik, I.; Datta, R. Topological analysis of catalytic reaction networks: Methanol decomposition on Pt (111). J. Catal. 2007, 252, 258–270. [Google Scholar] [CrossRef]
  29. Konsler, R.G.; Karl, J.; Jacobsen, E.N. Coperative asymmetric catalysis with dimeric salen complexes. J. Am. Chem. Soc. 1998, 120, 10780–10781. [Google Scholar] [CrossRef]
  30. Rosner, T.; le Bars, J.; Pfaltz, A.; Blackmond, D.G. Kinetic studies of Heck coupling reactions using palladacycle catalysts: Experimental and kinetic modelling of the role of dimer species. J. Am. Chem. Soc. 1998, 123, 1848–1855. [Google Scholar] [CrossRef]
  31. Noble-Teran, M.E.; Buhse, T.; Cruz, J.-M.; Coudret, C.; Micheau, J.-C. Nonlinear effects in asymmetric synthesis: A practical tool for the discrimination between monomer and dimer catalysis. ChemCatChem 2016, 8, 1836–1845. [Google Scholar] [CrossRef]
  32. Temkin, M. The kinetics of some industrial heterogeneous catalytic reactions. Adv. Catal. 1979, 28, 173–291. [Google Scholar]
  33. Horiuti, J.; Nakamura, T. On the theory of heterogeneous catalysis. Adv. Catal. 1967, 17, 1–74. [Google Scholar]
  34. Temkin, O.N. On Various Interconnections between Kinetics and Thermodynamics; LAP Lambert Academi Publishing: Saarbrueken, Germany, 2016; p. 119. [Google Scholar]
  35. Murzin, D.Y.; Touroude, R. On the rate of heterogeneous catalytic reactions with ionic intermediates. Catal. Lett. 1993, 20, 185–190. [Google Scholar] [CrossRef]
  36. Carucci, J.R.H.; Kurman, A.; Karhu, H.; Arve, K.; Eränen, K.; Wärnå, J.; Salmi, T.; Murzin, D.Y. Kinetics of the biofuels-assisted SCR of NOx over Ag/alumina-coated microchannels. Chem. Eng. J. 2009, 154, 34–44. [Google Scholar] [CrossRef]
  37. Murzin, D.Y. Interpretation of rate optima vs reaction parameters in steady state catalytic kinetics: Molecular aspects beyond concentration dependences. Mol. Catal. 2017, 433, 321–333. [Google Scholar] [CrossRef]
  38. Lazman, M.Z.; Yablonsky, G.S. Overall reaction rate equation of single-route complex catalytic reaction in terms of hypergeometric series. Adv. Chem. Eng. 2008, 34, 47–102. [Google Scholar]
  39. Laidler, K.J. Rate-controlling step: A necessary or useful concept? J. Chem. Educ. 1988, 65, 250–254. [Google Scholar] [CrossRef]
  40. Kozuch, S.; Shaik, S. How to conceptualize catalytic cycles? The energetic span model. Acc. Chem. Res. 2011, 44, 101–110. [Google Scholar] [CrossRef]
  41. Marin, G.B.; Yablonsky, G.S.; Constales, D. Kinetics of Chemical Reactions: Decoding Complexity, 2nd ed.; Wiley: Hoboken, NJ, USA, 2019. [Google Scholar]
  42. Temkin, M.I.; Murzin, D.Y.; Kul’kova, N.V. Kinetics and mechanism of liquid-phase hydrogenation. Doklady Akademii Nauk SSSR 1988, 303, 659–662. [Google Scholar]
  43. Murzin, D.Y.; Sokolova, N.A.; Kul’kova, N.V.; Temkin, M.I. Kinetics of liquid-phase hydrogenation of benzene and toluene on a nickel catalyst. Kinet. Catal. 1989, 30, 1352. [Google Scholar]
  44. Temkin, M.I.; Murzin, D.Y.; Kul’kova, N.V. Mechanism of liquid-phase hydrogenation of a benzene ring. Kinet. Catal. 1989, 30, 637. [Google Scholar]
  45. Murzin, D.Y.; Kul’kova, N.V. Kinetics and mechanism of the liquid-phase hydrogenation. Kinet. Catal. 1995, 36, 70–76. [Google Scholar]
  46. Murzin, D.Y.; Kul’kova, N.V. Non-equilibrium effects in the liquid-phase catalytic hydrogenation. Catal. Today 1995, 24, 35. [Google Scholar] [CrossRef]
  47. Behravesh, E.; Melander, M.M.; Wärnå, J.; Salmi, T.; Honkala, K.; Murzin, D.Y. Oxidative dehydrogenation of ethanol on gold: Combination of kinetic experiments and computation approach to unravel the reaction mechanism. J. Catal. 2020. [Google Scholar] [CrossRef]
  48. Mulligan, C.J.; Parker, J.S.; Hii, K.K. Revisiting the mechanism of the Fujiwa-Moritani reaction. React. Chem. Eng. 2020, 5, 1104–1111. [Google Scholar] [CrossRef]
Reactions 01 00004 g001 550
Figure 1. Potential energy diagram using Gibbs energy and Gibbs activation energy for the reaction S->P. T1–T3 represent transition states, and I1–I2 intermediates, I0 is the free form of catalyst. Adapted from [40].
Figure 1. Potential energy diagram using Gibbs energy and Gibbs activation energy for the reaction S->P. T1–T3 represent transition states, and I1–I2 intermediates, I0 is the free form of catalyst. Adapted from [40].
Reactions 01 00004 g001
Reactions 01 00004 g002 550
Figure 2. Three-step direct ethylene hydrogenation mechanisms.
Figure 2. Three-step direct ethylene hydrogenation mechanisms.
Reactions 01 00004 g002
Reactions 01 00004 g003 550
Figure 3. The minimum energy path for ethanol (EtOH) dehydrogenation. The black and orange bars correspond to minima and transition states, respectively. * denotes adsorbed species, while (g) refers to gas phase species. Square brackets enclose a transition state structure where ‘ - -’ marks the bond being broken. Reproduced with permission from [47].
Figure 3. The minimum energy path for ethanol (EtOH) dehydrogenation. The black and orange bars correspond to minima and transition states, respectively. * denotes adsorbed species, while (g) refers to gas phase species. Square brackets enclose a transition state structure where ‘ - -’ marks the bond being broken. Reproduced with permission from [47].
Reactions 01 00004 g003
Reactions 01 00004 g004 550
Figure 4. Reaction mechanism for the Fujiwara–Moritani reaction [48]. Published by The Royal Society of Chemistry.
Figure 4. Reaction mechanism for the Fujiwara–Moritani reaction [48]. Published by The Royal Society of Chemistry.
Reactions 01 00004 g004

Share and Cite

MDPI and ACS Style

Murzin, D.Y. Requiem for the Rate-Determining Step in Complex Heterogeneous Catalytic Reactions? Reactions 2020, 1, 37-46. https://0-doi-org.brum.beds.ac.uk/10.3390/reactions1010004

AMA Style

Murzin DY. Requiem for the Rate-Determining Step in Complex Heterogeneous Catalytic Reactions? Reactions. 2020; 1(1):37-46. https://0-doi-org.brum.beds.ac.uk/10.3390/reactions1010004

Chicago/Turabian Style

Murzin, Dmitry Yu. 2020. "Requiem for the Rate-Determining Step in Complex Heterogeneous Catalytic Reactions?" Reactions 1, no. 1: 37-46. https://0-doi-org.brum.beds.ac.uk/10.3390/reactions1010004

Article Metrics

Back to TopTop