Single-Route Linear Catalytic Mechanism: A New, Kinetico-Thermodynamic Form of the Complex Reaction Rate

For a complex catalytic reaction with a single-route linear mechanism, a new, kinetico-thermodynamic form of the steady-state reaction rate is obtained, and we show how its symmetries in terms of the kinetic and thermodynamic parameters allow better discerning their influence on the result. Its reciprocal is equal to the sum of n terms (n is the number of complex reaction steps), each of which is the product of a kinetic factor multiplied by a thermodynamic factor. The kinetic factor is the reciprocal apparent kinetic coefficient of the i-th step. The thermodynamic factor is a function of the apparent equilibrium constants of the i-th equilibrium subsystem, which includes the (n−1) other steps. This kinetico-thermodynamic form separates the kinetic and thermodynamic factors. The result is extended to the case of a buffer substance. It is promising for distinguishing the influence of kinetic and thermodynamic factors in the complex reaction rate. The developed theory is illustrated by examples taken from heterogeneous catalysis.


Introduction
How do we derive the reaction rate for a complex reaction, e.g., for a complex catalytic reaction? Based on graph theory, this general problem was posed by King and Atman [1] and, then, Volkenstein and Goldstein [2,3]. Yablonsky and Bykov, using graph theory as well, introduced the concept of relationships between the complex reaction rate equation and structure of the detailed mechanism [4,5]; see also the monographs [6,7]. The kinetic equation for a single-route catalytic reaction with a linear mechanism can always be presented in the form (see [1,2]): in which C c is the cycle characteristic. It is related to the overall reaction, C c = k + f + (C r ) − k − f − (C p ), in which k + = ∏ s k + s , k − = ∏ s k − s , k + /k − = K. The functions f + (C r ), f − (C p ) are the products of the reactant and product concentrations raised to certain powers, assuming that the rate of reaction follows the mass-action law formally, as if the overall reaction were an elementary reaction. The denominator, or using the phenomenological forms (see Equation (5.63) in [8]): Reaction Rate = Driving Force Kinetic Resistance (4) and: Kinetic Resistance = Driving Force Reaction rate (5) where the expression f + (C r ) − f − (C p )/K is a "driving force". Equations (4) and (5) are very similar to the well-known Ohm's Law considering the analogy between the reaction rate and current and between kinetic "driving force" and voltage as well.

Formulation of the Problem
Applying Equation (3) to this mechanism, the equation for the rate can also be written in the form: where is the quotient of the overall reaction and K is its equilibrium constant. The expression (1 − Q/K) is a driving force, which does not depend on the details of the mechanism, but only on the overall reaction. As fork i , it is the apparent forward rate constant of the i-th step, which can include the concentration of a gas substance as a parameter. Furthermore, T i is a function obtained under the assumption that all steps except for the i-th one are under equilibrium conditions. Therefore, every i-th step is characterized by the ensemble of equilibrium reactions. Previously, the concept of the ensemble of equilibrium subsystems was proposed as an efficient method for the mathematical analysis of nonlinear models (see the monographs [8][9][10] and the papers [11,12]). In fact, the concept of equilibrium subsystems is a generalization of the physico-chemical concept of the equilibrium step, which is well accepted in chemical kinetics. An equilibrium subsystem can be defined by two requirements: that (n − 1) steps are at equilibrium and that one step is rate limiting. In the papers [9,10], it was shown that the solutions of these n equilibrium subsystems (i.e., all roots, not just one root) determine the coefficients of the kinetic polynomial. In the theory of steady-state kinetic models, the kinetic polynomial is the general description of a single-route catalytic reaction.

Research Analysis
The advantage of Equation (6) is that it is fully symmetric in the mechanism's parameter values, every term presenting the product of a pure kinetic factor 1/k i by a purely thermodynamic term T i . If one or more steps are irreversible, Proof of Equation (6). For this proof, we assume that all gas concentrations are absorbed in the apparent rate constants, so that the mechanism is of the form Z 1 with global quotient Q = 1. The rate of the i-th step is denoted R i and the i-th intermediate as for the value in {1, . . . , n} that is equal to a modulo n. We then have that ; the equilibrium constant of the i-th step is K i = k + i /k − i , and the balance of intermediates is Z 1 + · · · + Z n = 1. In terms of (n + 1) × (n + 1) matrices, for instance, when n = 4, The determinant of the square matrix must therefore be zero. Expanding it along the rightmost column, a 1 R 1 + · · · + a n R n = 1 must hold for some values a 1 , . . . , a n . The equilibrium rate of the i-th step being denoted R i,e , by definition, one must have that a i = 1/R i,e ; consequently, Calculating the coefficients a i from the determinant, where . The concentrations then are determined by:

Single-Route Catalytic Reaction: All Reactions Are Irreversible
In this case, Equation (6) is extremely simplified;

Single-Route Catalytic Reaction: One
Step Is Irreversible and Rate Limiting.
Assuming that the rate limiting step's 1/k l is so large that the other terms with i = l can be neglected,

Single-Route Catalytic Reaction with a Buffer Step
Let Z 0 be a buffer attached to Z 1 . In the steady-state, the reaction step between Z 0 and Z 1 must be balanced; hence, Z 0 = Z 1 /K 0 always.
This is reflected in the formulas: when n = 2, and for general n,

Two Cycles that Have a Common Intermediate
For n = 2, see [8], and: (14)

Two
Step Single Route Mechanism (Water-Gas Shift Mechanism) where Z is the active catalytic center. The overall reaction is CO + H 2 O CO 2 + H 2 : wherek

Three
Step Mechanism of Sulfur Dioxide Oxidation The overall equation is 2 SO 2 + O 2 2 SO 3 . wherek

Discussion
The new kinetico-thermodynamic expressions for the complex catalytic reaction rate presented here, Equation (8) and similar ones, are convenient for analyzing different aspects of kinetic behavior.
First, it is possible to characterize the contributions of different steps because the proposed expression is a sum of terms that correspond to different steps. Based on this expression, the problem of distinguishing the domain of validity of a rate limiting assumption can be posed, at least qualitatively [13].
In a special case, if all steps are irreversible, the answer is obvious (see Equation (10)), and the i-th reaction will be limiting ask i tends to zero. Generally, two interwoven problems arise: Which step is rate limiting? What is the domain of the validity of the rate limiting assumption?
It is possible to simplify Equation (8) by reducing the number of terms that correspond to different steps. Our qualitative analysis shows that a step can be eliminated when both the apparent forward kinetic parameter and the apparent equilibrium constant are large enough. However, it is possible that some step may have a moderate (finite) equilibrium constant despite having both the kinetic parameters, forward and reverse, large. In this case, the corresponding rate term of this step can be eliminated, but the equilibrium parameter of this step will remain present in the other rate terms. In most cases, determining the rate-determining step requires an analysis based on concrete parameters. Every term of Equation (8) is an explicit function of the temperature and concentrations: the apparent kinetic coefficientk i is characterized by an Arrhenius dependence on the temperature and a linear dependence on the gas concentration, and the thermodynamic term T i is a function of the temperature as well: it is increasing or decreasing with temperature if all present steps are exothermic or endothermic, respectively. Regarding the concentration dependences, T i is proportional to the reactant concentrations and inversely proportional to the product concentrations.

Conclusions
A new, kinetico-thermodynamic form expressing the rate of a complex catalytic single route equation is obtained. According to this form, the reciprocal rate equals the sum of n terms (n is the number of steps), and every i-th term equals the product of the i-th kinetic factor, i.e., the reciprocal i-th apparent kinetic coefficient, multiplied by the i-th apparent thermodynamic term. The latter is a function of the apparent equilibrium constants of the (n − 1) other steps. Generally, the kinetico-thermodynamic form separates kinetic and thermodynamic factors. This form allows estimating the contribution of every step from both kinetic and thermodynamic aspects.