Processing math: 100%

Divergent First-order Kinetics

Scheme:

Bk1Ak2C

Differential Equation:

A(t)=-k1A-k2A

B(t)=k1A

C(t)=k2A

Mass balance equation:

A0=A+B+C

Initial conditions:

A(0)=A0, B(0)=0 and C(0)=0

Solution:

A(t)=A0exp(-(k1+k2)t)

B=k1A0-(k1+k2)exp(-(k1+k2)t)+k1A0k1+k2

C=k2A0-(k1+k2)exp(-(k1+k2)t)+k2A0k1+k2

Graph:

Plotting A0=1Mol, B0=0, C0=0, k1=1e-2, k2=2e-2


Derivation:

The solution can be obtained directly by integration of the separable equation. First A(t),

A(t)=-k1A-k2A=-(k1+k2)A

AtA0dAA=-(k1+k2)t

lnAt=lnA0-(k1+k2)t

A(t)=A0exp(-(k1+k2)t)

With the solution for A(t) in terms of t only, we can solve B(t),

B(t)=k1A=k1A0exp(-(k1+k2)t)

We solve this integral with u-substitution. With u=-(k1+k2)t, then dudt=-(k1+k2); applying the chain rule we have,

dBdt=-(k1+k2)dBdu=k1A0exp(u)

which becomes,

-(k1+k2)dB=k1A0exp(u)du

We then integrate the expression,

-(k1+k2)dB=k1A0exp(u)du

...and obtain,

-(k1+k2)B=k1A0exp(u)+c1

Reversing the substitution for u,

-(k1+k2)B=k1A0exp(-(k1+k2)t)+c1

Making B the subject,

B=k1A0-(k1+k2)exp(-(k1+k2)t)+c1

We then solve the initial value problem by noting that at t=0, B=0,

0=k1A0-(k1+k2)+c1

Hence we can solve for c1,

c1=k1A0k1+k2

The solution for B, becomes,

B=k1A0-(k1+k2)exp(-(k1+k2)t)+k1A0k1+k2

We then perform the same procedure for C(t), however we note that this differential equation differs from that of B(t) only by replacement of the rate constant k1 with k2. Hence the solution is,

C=k2A0-(k1+k2)exp(-(k1+k2)t)+k2A0k1+k2