Chemical kinetics equations

Chemical kinetics equations for the closed heterogeneous gas-solid systems are of the type 

To discriminate between reactions on solid surfaces and in the gas phase, we have introduced different indices, i.e. s for the former reactions and a for the latter. 

Equations (17) and (18) describe the process of complex homogeneous-heterogeneous reactions. 

If the isomerization rate in the gas phase is w0 and that on the surface is w4, w2, w3, then the kinetic equations can be written as 

Kinetic equations can be reduced to a more compact form using a stoichiometric matrix and writing the rates for the various steps as a vector column. Then 

---

It turns out that there is another branch of mathematics, closely related to tire calculus of variations, although historically the two fields grew up somewhat separately, known as optimal control theory (OCT). Although the boundary between these two fields is somewhat blurred, in practice one may view optimal control theory as the application of the calculus of variations to problems with differential equation constraints. OCT is used in chemical, electrical, and aeronautical engineering where the differential equation constraints may be chemical kinetic equations, electrical circuit equations, the Navier-Stokes equations for air flow, or Newton s equations. In our case, the differential equation constraint is the TDSE in the presence of the control, which is the electric field interacting with the dipole (pemianent or transition dipole moment) of the molecule [53, 54, 55 and 56]. From the point of view of control theory, this application presents many new features relative to conventional applications perhaps most interesting mathematically is the admission of a complex state variable and a complex control conceptually, the application of control teclmiques to steer the microscopic equations of motion is both a novel and potentially very important new direction. 

Solving the chemical kinetic equations and comparison with the observed molecular abundance. 

In the case of classic chemical kinetics equations, one can get in a few cases analytical solution for the set of differential equations in the form of explicit expressions for the number or weight fractions of i-mcrs (cf. also treatment of distribution of an ideal hyperbranched polymer). Alternatively, the distribution is stored in the form of generating functions from which the moments of the distribution can be extracted. In the latter case, when the rate constant is not directly proportional to number of unreacted functional groups, or the mass action law are not obeyed, Monte-Carlo simulation techniques can be used (cf. e.g. [2,3,47-52]). This technique was also used for simulation of distribution of hyperbranched polymers [21, 51, 52],... 

Most electrode reactions of interest to the organic electrochemist involve chemical reaction steps. These are often assumed to occur in a homogeneous solution, that is, not at the electrode surface itself. They are described by the usual chemical kinetic equations, for example, first- or second-order reactions and may be reversible (chemical reversibility) or irreversible. 

To illustrate new effects which could be observed employing this technique, consider a simple model [67, 118]. Start from the chemical kinetics equation... 

Its obvious peculiarity as compared with the standard chemical kinetics, equation (2.1.10), is the emergence of the fluctuational second term in r.h.s. The stochastic reaction description by means of equation (2.2.37) permits us to obtain the equation for dispersions crjj which, however, contains higher-order momenta. It leads to the distinctive infinite set of deterministic equations describing various average quantities, characterizing the fluctuational spectrum. 

Fig. 5.18. The non-stationary part of the reaction rate K(t) for the transient kinetics of the A + A 0 reaction (curves 1 to 3) and A + B —> 0 reaction (curve 4), d = 3 (random initial distribution) [101]. The broken line shows prediction of the standard chemical kinetics, equation (5.3.10). The initial concentrations are given.

Only in the case when the reaction vessel is a long tube do the laws of chemical reaction in a jet approach the laws of chemical kinetics in a closed vessel uncomplicated by diffusion exchange, i.e., the laws which are obtained by integration of chemical kinetics equations of the form... 

Oregonator and "brusselator studied in detail by the Prigogine school were nevertheless extremely speculative schemes. A study of the behaviour of classical chemical kinetics equations assumed a high priority in order to select the structure responsible for the appearance of critical effects. The results of such a study, described in Chap. 3, can be applied to interpret critical effect experiments. 

A.I. Volpert, Qualitative Methods for Chemical Kinetics Equations, Institut Khimicheskoi Fiziki, Akad. Nauk SSSR, Chernogolovka, 1976 (in Russian). 

The phase space of eqn. (73) is the space of vectors c. Its points are specified by the coordinates cx,. . . , cn. The set of phase space points is the set of all possible states of the system. Phase space can be not only the whole vector space but also a certain part. Thus in chemical kinetic equations, variables are either concentrations or quantities of substances in the system. Their values cannot be negative. It is therefore natural to restrict ourselves to the set of those c all the components of which are not negative, i.e. Ci > 0. In what follows we shall refer to these d values as non-negative. Hence positive are those c values all the components of which are positive, i.e. Cj > 0. 

Chemical kinetic equations possess the following properties. For any non-negative initial conditions, c0 the only solution of eqn. (73) that exists is c(t, %, c0). At the initial moment it takes values of c0, i.e. c(0, d0) = c0. 

In what follows when discussing the general properties of the chemical kinetic equations, we will assume that the additional laws of conservation (if there are any) have b en discovered and the respective values of x are included in the matrix A as additional rows. 

So far (Sect. 1) we have discussed only approaches to derive chemical kinetic equations for closed systems, i.e. those having no exchange with the environment. Now let us study their dynamic properties. For this purpose let us formulate the basic property of closed chemical systems expressed by the principle of detailed equilibrium a rest point for the closed system is a point of detailed equilibrium (PDE), i.e. at this point the rate of every step equals zero... 

In the previous section we introduced the Lyapunov functions for chemical kinetic equations that are the dissipative functions G. The function RTG is treated as free energy. Since G < 0 and the equality is obtained only at PDE, and for the construction of G it suffices to know only the position of equilibrium N, there exist limitations on the non-steady-state behaviour of a closed system that are independent of the reaction mechanism. If in the initial composition N = N, the other composition N can be realized during the reaction only in the case when... 

This is a general fact. For monomolecular (or pseudo-monomolecular) reactions the graphs corresponding to compartments are acyclic. A similar property for the systems having either bi- or termolecular reactions is more complex. It can be formulated as follows. If every edge in the graph of predominant reaction directions for some compartment is ascribed to a positive "rate constant k and chemical kinetic equations are written with... 

The multitude constructed for point N°, and designated as J(29°), is co-invariant for all systems of chemical kinetic equations obtained in accordance with the mechanism prescribed and having N as a PDE. Moreover, it is minimal among the multitudes possessing this property, i.e. if a multitude that is co-invariant for all systems with a given reaction mechanism and an equilibrium point, contains 29°, it also has J(N°). In the general case, it is constructed as follows [33]. 

So far the quasi-steady-state hypothesis introduced in 1913 has remained the most favourable approach to operating with chemical kinetic equations. In short (and not quite strictly), its most applicable version can be formulated as follows. During the reaction, the concentrations of some (usually intermediate) compounds are the concentration functions of the other (usually observed) substances and "adapt to their values as if they were steady-state values. 

Specificity of a concrete system accounts for the source of the appearance of a small parameter and for its type. For homogeneous reactions, a small parameter is usually a ratio of rate constants for various reactions some reactions are much faster than the others. For just such a small parameter Vasiliev et al. [25] distinguished a class of chemical kinetic equations for which the application of the quasi-stationarity principle is correct (they considered a closed system). 

Studies of linear systems and systems without "intermediate interactions show that a positive steady state is unique and stable not only in the "thermodynamic case (closed systems). Horn and Jackson [50] suggested one more class of chemical kinetic equations possessing "quasi-ther-modynamic properties, implying that a positive steady state is unique and stable in a reaction polyhedron and there exist a global (throughout a given polyhedron) Lyapunov function. This class contains equations for closed systems, linear mechanisms, and intersects with a class of equations for "no intermediate interactions reactions, but does not exhaust it. Let us describe the Horn and Jackson approach. 

A.N. Gorban , Round the Equilibrium. Chemical Kinetics Equations and their Thermodynamic Analysis, Nauka, Novosibirsk, 1984 (in Russian). 

The results of the numerical experiment for system (20) necessitated a general mathematical investigation of slow relaxations in chemical kinetic equations. This study was performed by Gorban et al. [226-228] who obtained several theorems permitting them to associate the existence of slow relaxations in a system of chemical kinetic equations (and, in general, in dynamic systems) with the qualitative changes in the phase portrait with its parameters (see Chap. 7). 

Let the phase portrait of the system be characterized by some set of co-limit points. The concepts of an "co-limit point and an "co-limit set have been extensively used in the theory of dynamic systems. The thing is that the trajectory does not necessarily enter into a steady state. In the general case (as well as in the case of chemical kinetic equations), the existence of limit cycles is possible. The letter co is a symbol for the region of the phase space into which at t—>co the trajectory tends ("from a to co ). Let x0 be a vector... 

Finally, we can suggest a third explanation fast steps can compose a mechanism with slow relaxations. Indeed, nothing suggests that the relaxation time for a set of chemical kinetic equations is directly dependent on the characteristic times of the individual steps. But it cannot be treated as a reason for slow relaxations. It is only a simple indication for the possibility of finding such reasons here. Let us now indicate the reasons according to which fast steps can compose a mechanism with slow relaxations. 

If a system of chemical kinetic equations is non-linear and the reaction mechanism includes an interaction step between various substances, bifurcations are possible. They account for the effects of critical retardation. Let us illustrate this by the simplest (non-chemical) example. Consider the differential equation... 

Generally, the gas-phase concentrations Nk are described by a set of master equations of chemical kinetic (equations for active masses at the source/sink of the kth substance in the chemical reactions Wk)... 

The notion of thermodynamic tree (the graph, each point of which represents the set of thermodynamically equivalent states) was introduced by Gorban (1984) where he also revealed the possibilities of applying this notion for analysis of the chemical kinetics equations. In the work by Gorban et al. (2001, 2006) the authors consider the problems of employing thermodynamic tree to study the physicochemical systems using MEIS. 

Most property kinetic studies reported in the literature are conducted by analogy with the methodology of chemical kinetics. A physical property, P, observed to change monotonically with time is assumed to obey a differential expression similar to a rate law in chemical kinetics. Equation 1 is a general expression of this kind where k is a constant and... 

Equation (47) may be treated as the first term in a perturbation about equilibrium flow [27], [81]. Given a scheme for the chemical kinetics, equation (1-8) may be employed to generate an expansion of in the departures of T and of from their chemical-equilibrium dependences on Z. A general equation for the departures from equilibrium may then be derived from equation (45) [27], [81], and suitable modeling in this equation then provides the first correction to vv.. The procedure, never simple, is least complicated if the departure from chemical equilibrium can be characterized fully in terms of a single variable. More research remains to be done on the subject. 

---