Steady state (chemistry)

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

<templatestyles src="Module:Hatnote/styles.css"></templatestyles>

In chemistry, a steady state is a situation in which all state variables are constant in spite of ongoing processes that strive to change them. For an entire system to be at steady state, i.e. for all state variables of a system to be constant, there must be a flow through the system (compare mass balance). A simple example of such a system is the case of a bathtub with the tap running but with the drain unplugged: after a certain time, the water flows in and out at the same rate, so the water level (the state variable Volume) stabilizes and the system is in a steady state.

The steady state concept is different from chemical equilibrium. Although both may create a situation where a concentration does not change, in a system at chemical equilibrium, the net reaction rate is zero (products transform into reactants at the same rate as reactants transform into products), while no such limitation exists in the steady state concept. Indeed, there does not have to be a reaction at all for a steady state to develop.

The term steady state is also used to describe a situation where some, but not all, of the state variables of a system are constant. For such a steady state to develop, the system does not have to be a flow system. Therefore such a steady state can develop in a closed system where a series of chemical reactions take place. Literature in chemical kinetics usually refers to this case, calling it steady state approximation.

In simple systems the steady state is approached by state variables gradually decreasing or increasing until they reach their steady state value. In more complex systems state variable might fluctuate around the theoretical steady state either forever (a limit cycle) or gradually coming closer and closer. It theoretically takes an infinite time to reach steady state, just as it takes an infinite time to reach chemical equilibrium.

Both concepts are, however, frequently used approximations because of the substantial mathematical simplifications these concepts offer. Whether or not these concepts can be used depends on the error the underlying assumptions introduce. So, even though a steady state, from a theoretical point of view, requires constant drivers (e.g. constant inflow rate and constant concentrations in the inflow), the error introduced by assuming steady state for a system with non-constant drivers may be negligible if the steady state is approached fast enough (relatively speaking).

Steady state approximation in chemical kinetics

The steady state approximation,[1] occasionally called the stationary-state approximation, involves setting the rate of change of a reaction intermediate in a reaction mechanism equal to zero.

It is important to note that steady state approximation does not assume the reaction intermediate concentration to be constant (and therefore its time derivative being zero), it assumes that the variation in the concentration of the intermediate is almost zero: the concentration of the intermediate is very low, so even a big relative variation in its concentration is small, if considered quantitatively.

Its use facilitates the resolution of the differential equations that arise from rate equations, which lack an analytical solution for most mechanisms beyond the most simple ones. The steady state approximation is applied, for example in Michaelis-Menten kinetics.

As an example, the steady state approximation will be applied to two consecutive, irreversible, homogeneous first order reactions in a closed system. (For heterogeneous reactions, see reactions on surfaces.) This model corresponds, for example, to a series of nuclear decompositions like <ce> ^{239}U -> ^{239}Np -> ^{239}Pu\! </ce>.

If the rate constants for the following reaction are k_1 and k_2; <ce> A -> B -> C </ce>, combining the rate equations with a mass balance for the system yields:

Reaction rates

For species A: Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): \frac{d[\ce A]}{dt} = -k_1 [\ce A]


For species B: Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): \frac{d[\ce B]}{dt} = k_1 [\ce A] - k_2 [\ce B] , Here the first (positive) term represents the formation of B by the first step <ce> A -> B</ce>, whose rate depends on the initial reactant A. The second (negative) term represents the consumption of B by the second step <ce> B -> C</ce>, whose rate depends on B as the reactant in that step.

For species C: Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): \frac{d[\ce C]}{dt} = k_2 [\ce B] , the rate of formation of C by the second step.

Analytical Solutions

The analytical solutions for these equations (supposing that initial concentrations of every substance except for A are zero) are:[2]

Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): [\ce A]=[\ce A]_0 e^{-k_1 t}
Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): \left[ \ce B \right]= \begin{cases} \left[ \ce A \right]_{0}\frac{k_{1}}{k_{2}-k_{1}}\left( e^{-k_{1}t}-e^{-k_{2}t} \right);&k_{1}\ne k_{2} \\\\ \left[ \ce A \right]_{0}k_{1}te^{-k_{1}t};&\text{otherwise} \\ \end{cases}


Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): \left[ \ce C \right]= \begin{cases} \left[ \ce A \right]_{0}\left( 1+\frac{k_{1}e^{-k_{2}t}-k_{2}e^{-k_{1}t}}{k_{2}-k_{1}} \right);&k_{1}\ne k_{2} \\\\ \left[ \ce A \right]_{0}\left( 1-e^{-k_{1}t}-k_{1}te^{-k_{1}t} \right);&\text{otherwise} \\ \end{cases}


Steady State

If the steady state approximation is applied, then the derivative of the concentration of the intermediate is set to zero.

Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): \frac{d[\ce B]}{dt} = 0 = k_1 [\ce A] - k_2 [\ce B] \Rightarrow \; [\ce B] = \frac{k_1}{k_2} [\ce A]
therefore Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.):  \frac{d[\ce C]}{dt} =  k_1 [\ce A]
so Failed to parse (Missing <code>texvc</code> executable. Please see math/README to configure.): [\ce C]=[\ce A]_0 \left (1- e^{-k_1 t} \right )

.

Validity

Concentration vs. time. Concentration of intermediate in green, product in blue and substrate in red
(k_2/k_1=0.5)
Concentration vs. time. Concentration of intermediate in green, product in blue and substrate in red
(k_2/k_1=10)

The analytical and approximated solutions should now be compared in order to decide when it is valid to use the steady state approximation. The analytical solution transforms into the approximate one when  k_2 \gg k_1 , because then e^{-k_2t} \ll e^{-k_1t} and k_2-k_1 \approx \; k_2. Therefore it is valid to apply the steady state approximation only if the second reaction is much faster than the first one (k2/k1 > 10 is a common criterion), because that means that the intermediate forms slowly and reacts readily so its concentration stays low.

The graphs show concentrations of A (red), B (green) and C (blue) in two cases, calculated from the analytical solution:

  • When the first reaction is faster it is not valid to assume that the variation of [B] is very small, because [B] is neither low or close to constant: first A transforms into B rapidly and B accumulates because it disappears slowly. As the concentration of A decreases its rate of transformation decreases, at the same time the rate of reaction of B into C increases as more B is formed, so a maximum is reached when t=\begin{cases}
   \frac{\ln \left( \frac{k_{1}}{k_{2}} \right)}{k_{1}-k_{2}} & \, k_{1}\ne k_{2}  \\\\
   \frac{1}{k_{1}} & \, \text{otherwise}  \\
\end{cases}.
    From then on the concentration of B decreases.
  • When the second reaction is faster, after a short induction period, concentration of B remains low (and more or less constant) because its rate of formation and disappearance are almost equal and the steady state approximation can be used.

The equilibrium approximation can be used sometimes in chemical kinetics to yield similar results as the steady state approximation (Michaelis-Menten kinetics can be derived assuming equilibrium instead of steady state): it consists in assuming that the intermediate is at chemical equilibrium. Normally the requirements for applying the steady state approximation are laxer: the concentration of the intermediate is only needed to be low and more or less constant (as seen, this has to do only with the rates at which it appears and disappears) but it is not needed to be at equilibrium, which is usually difficult to prove and involves heavier assumptions.

External links

Notes and references

  1. IUPAC Gold Book definition of steady state
  2. P. W. Atkins "Physical Chemistry" For a complete derivation of the exact and approximate solution and a comparison between them

See also