Matrix formulation of chemical reaction rates: A ... - ACS Publications

Pharmaceutical companies employ teams of lawyers and intellectual property experts to make sure synthetic... BUSINESS ... American Chemical Society ...
1 downloads 0 Views 1MB Size
Matrix Formulation of Chemical Reaction Rates A Mathematical Chemical Exercise Lionello poglianil and Mario Terenzi Dipartimento di Chimica, Universita della Calabria, 87030 Rende (CS), Italy Some years ago, L. Papula2presented a method to obtain the kinetic equations of a consecutive reaction scheme involving only irreversible fwst-order steps. The method was based on matrix algebra. In fact, matrix algebra offers an elegant way of solving the differential equations for simple kinetic systems. The use of matrix algebra has an advantage in that it allows the simultaneous solution for the rates of the reaction for the different chemical species. Accordingly, we present here an application of matrix algebra to first-order rate equations. One of the cases treated (side opposing reactions) is also treated in a different way by N a q ~ i who ,~ presents a more general treatment on the similarity between chemical kineties and vibration modes. First-Order Reactions Consider the reaction

yields the eigenvalues, w{s wl = 0 and wz = --k

To the eigenvalue wl = 0, belongs the following system of homogeneous linear equations

solutions of which are Al = 0 and Az that can be regarded as a free parameter and be chosen equal to A, withA some nonzero constant. Eigenvector C1, which belongs to eigenvalue w,, then has the following form

.. To the eigenvalue wz = -k, belongs the following system. X-+Y

Let the molecular concentrations of X and Y at time t be Cx(t)= Cx and Cdt) = Cy

whose solution is A,=-Az=B

respectively The rate of formation of X and Y are thus and C;

where B is a nonzero constant. Thus, eigenvedor Cz,which belongs to eigenvalue w2,is

where

Then this reaction obeys the following system of homogeneuos linear differential equations of the first order with constant coefficient.

k)= e qq

(I)

This first-order rate equation can be rewritten in the following compact matrix form. C'=AC

(2)

The general solution of the fxst-order kinetic problem is then the followine linear combination of the two indeoendent C1 and C2 eigenvedors.

-

With the aid of the initial conditions, Cx(0)=a and Cy(0) = b, we obtain B = a and A = a + b The final chemical kinetics equations will be

which has a general solution

with C' = wC. By the aid ofthis solution, eq 2 leadus to the following matrix eigenvalue problem. wC=AC (4) Since wC = wIC, where I is the unit matrix (A-wI)C=O

Opposing First-Order Reactions Consider the reaction

(5)

The matrix form of the differential rate equations has the following form.

For nontrivial solutions to this equation to exist, the determinant I A - wI I must be zero. Solution of

(13)

'Author to whom correspondence should be addressed. 'Papula, L. Mathernatik fur Chemiker; Ferdinand Enke Verlag: Stuttgattl975. 3Naqvi,K. R. J. Chem. Ed. 1989, 66(9),703.

Here the coefficient square matrix is the result of the sum of two coefficient square matrices of the fust-order reaction type. Note that the two diagonals have opposite signs. This equation lead us t o the same kind of matrix

278

Journal of Chemical Education

eigenvalue problem of the first-order reaction, with the same kind of nontrivial solutions. The wi eigenvalues are found by solving the following determinantal equation. (14)

Solving, we obtain

Diverging Side Reactions

Consider the reaction

e x + z

Y

The coefficient square matrix that describes this system has the following form

w1=0

and w2 = i k , + k,)

k, 0 0

The C1 and C z eigenvetors calculated from

(25)

with k=k,+kz

With the aid of the same kind of mathematical approach and with the initial conditions

have the following form

Cx(0)= a

CdO) = b

CZ(0)= e

we obtain the following rate equations Cx = a exp-kt)

.

.

where

where A and B are some nonzero constant and where k K = A and h = k f + k , k,

and

For the general solution to this problem c = R ] = (KA A +-BBexp(-fit) exp-ht)

CZ= e +K" ( 1- exp-kt))

where (18)

with the aid of the initial conditions

Converging Side Reactions

Cx(0)= a and CdO) = b

Consider the reaction

and with

Y(akf- bk,) m= h

X

e-z

The coefficient square matrix that describes this system has the following form.

we get the final chemical rates (29)

The final rate equations are Cy=-K(a b, - m exp-ht) +

(l+K)

where

Consecutive First-Order Reactions

Consider a simple, consecutive reaction involving only irreversible fust-order steps

and Cy = b exp(-k',t)

X + Y + Z

Papula's original treatment starts with the following matrix.

where Cx(0)= a

(21)

(31)

CdO) = b

Cz(0)= c

Opposing Side Reactions

Consider the following hypothetical reaction.

Then, with The resulting coefficient square matrix is the sum of the two previous square matrices.

it yields the following final kinetic equations

c, = a kl exp(-kit) k

+ (bk - akl)k exp-k2t)

(23)

-k

k',

k,

0 -k;

Volume 69 Number 4 April 1992

(33)

279

For kl=k,=kd

(d stands far divergent)

and k', = k', = k c

(C

A rapid comparison between Cx and Naqvi's3 Bz shows the similarity of the two different approaches.

stands for convergent)

with CAO) = a

Cy(0)= b

Conclusion

CZ(0)= c

and with

the fmal rate equations will be (34)

Ks

b-c

+ 2 -P(-@) cv = ZK+~

280

a-b-c - X2(2K+ 1)

Journal of Chemical Education

We notice that the number of rows in our matrices eauals the number of chemical species in the chemical ecpacons. Also. the coefficientmatrix is a sauare matrix that can be used to fully characterize the chekcal reaction. While the simple formalism presented applies fairly well for chemical equations with first-order steps only, a more elaborated matrix method was developped very recently for kinetic systems composed by steps of any order? Berberan-Santos. M. N.; Maltinho, J. M. G. J. Chem. E h c .

-(Fht'

(35)

1990, 67(5),375.