Chemistry:Rate equation
In chemistry, the rate law or rate equation for a reaction is an equation that links the initial or forward reaction rate with the concentrations or pressures of the reactants and constant parameters (normally rate coefficients and partial reaction orders).^{[1]} For many reactions, the initial rate is given by a power law such as
- [math]\displaystyle{ v_0\; =\; k[\mathrm{A}]^x[\mathrm{B}]^y }[/math]
where [math]\displaystyle{ [\mathrm{A}] }[/math] and [math]\displaystyle{ [\mathrm{B}] }[/math] express the concentration of the species [math]\displaystyle{ \mathrm{A} }[/math] and [math]\displaystyle{ \mathrm{B}, }[/math] usually in moles per liter (molarity, [math]\displaystyle{ M }[/math]). The exponents [math]\displaystyle{ x }[/math] and [math]\displaystyle{ y }[/math] are the partial orders of reaction for [math]\displaystyle{ \mathrm{A} }[/math] and [math]\displaystyle{ \mathrm{B} }[/math] and the overall reaction order is the sum of the exponents. These are often positive integers, but they may also be zero, fractional, or negative. The constant [math]\displaystyle{ k }[/math] is the reaction rate constant or rate coefficient of the reaction. Its value may depend on conditions such as temperature, ionic strength, surface area of an adsorbent, or light irradiation. If the reaction goes to completion, the rate equation for the reaction rate [math]\displaystyle{ v\; =\; k[\ce{A}]^x[\ce{B}]^y }[/math] applies throughout the course of the reaction.
Elementary (single-step) reactions and reaction steps have reaction orders equal to the stoichiometric coefficients for each reactant. The overall reaction order, i.e. the sum of stoichiometric coefficients of reactants, is always equal to the molecularity of the elementary reaction. However, complex (multi-step) reactions may or may not have reaction orders equal to their stoichiometric coefficients. This implies that the order and the rate equation of a given reaction cannot be reliably deduced from the stoichiometry and must be determined experimentally, since an unknown reaction mechanism could be either elementary or complex. When the experimental rate equation has been determined, it is often of use for deduction of the reaction mechanism.
The rate equation of a reaction with an assumed multi-step mechanism can often be derived theoretically using quasi-steady state assumptions from the underlying elementary reactions, and compared with the experimental rate equation as a test of the assumed mechanism. The equation may involve a fractional order, and may depend on the concentration of an intermediate species.
A reaction can also have an undefined reaction order with respect to a reactant if the rate is not simply proportional to some power of the concentration of that reactant; for example, one cannot talk about reaction order in the rate equation for a bimolecular reaction between adsorbed molecules:
- [math]\displaystyle{ v_0=k \frac{K_1K_2C_AC_B}{(1+K_1C_A+K_2C_B)^2}. }[/math]
Definition
Consider a typical chemical reaction in which two reactants A and B combine to form a product C:
- [math]\displaystyle{ \ce{{A} + {2B} -\gt {3C}}. }[/math]
This can also be written
- [math]\displaystyle{ 0 = -\mathrm{A} - 2\mathrm{B} + 3\mathrm{C}. }[/math]
The prefactors −1, −2 and 3 (with negative signs for reactants because they are consumed) are known as stoichiometric coefficients. One molecule of A combines with two of B to form 3 of C, so if we use the symbol [X] for the number of moles of chemical X,^{[2]}
- [math]\displaystyle{ -\frac{d[\mathrm{A}]}{dt} = - \frac{1}{2} \frac{d[\mathrm{B}]}{dt} = \frac{1}{3} \frac{d[\mathrm{C}]}{dt}. }[/math]
If the reaction takes place in a closed system at constant temperature and volume, without a build-up of reaction intermediates, the reaction rate [math]\displaystyle{ v }[/math] is defined as
- [math]\displaystyle{ v = \frac{1}{\nu_i} \frac{d[\mathrm{X}_i]}{dt}, }[/math]
where ν_{i} is the stoichiometric coefficient for chemical X_{i}, with a negative sign for a reactant.^{[3]}
The initial reaction rate [math]\displaystyle{ v_0 = v(t = 0) }[/math] has some functional dependence on the concentrations of the reactants,
- [math]\displaystyle{ v_0 = f\left([\mathrm{A}], [\mathrm{B}], \ldots\right), }[/math]
and this dependence is known as the rate equation or rate law.^{[4]} This law generally cannot be deduced from the chemical equation and must be determined by experiment.^{[5]}
Power laws
A common form for the rate equation is a power law:^{[5]}
- [math]\displaystyle{ v_0 = k[\ce A]^x[\ce B]^y \cdots }[/math]
The constant [math]\displaystyle{ k }[/math] is called the rate constant. The exponents, which can be fractional,^{[5]} are called partial orders of reaction and their sum is the overall order of reaction.^{[6]}
In a dilute solution, an elementary reaction (one having a single step with a single transition state) is empirically found to obey the law of mass action. This predicts that the rate depends only on the concentrations of the reactants, raised to the powers of their stoichiometric coefficients.^{[7]}
Determination of reaction order
Method of initial rates
The natural logarithm of the power-law rate equation is
- [math]\displaystyle{ \ln v_0 = \ln k + x\ln[\ce A] + y\ln[\ce B] + \cdots }[/math]
This can be used to estimate the order of reaction of each reactant. For example, the initial rate can be measured in a series of experiments at different initial concentrations of reactant [math]\displaystyle{ \rm A }[/math] with all other concentrations [math]\displaystyle{ [\rm B], [\rm C], \dots }[/math] kept constant, so that
- [math]\displaystyle{ \ln v_0 = x\ln[\ce A] + \textrm{constant}. }[/math]
The slope of a graph of [math]\displaystyle{ \ln v }[/math] as a function of [math]\displaystyle{ \ln [\ce A] }[/math] then corresponds to the order [math]\displaystyle{ x }[/math] with respect to reactant [math]\displaystyle{ \rm A }[/math].^{[8]}^{[9]}
However, this method is not always reliable because
- measurement of the initial rate requires accurate determination of small changes in concentration in short times (compared to the reaction half-life) and is sensitive to errors, and
- the rate equation will not be completely determined if the rate also depends on substances not present at the beginning of the reaction, such as intermediates or products.
Integral method
The tentative rate equation determined by the method of initial rates is therefore normally verified by comparing the concentrations measured over a longer time (several half-lives) with the integrated form of the rate equation; this assumes that the reaction goes to completion.
For example, the integrated rate law for a first-order reaction is
- [math]\displaystyle{ \ln{[\ce A]} = -kt + \ln{[\ce A]_0}, }[/math]
where [math]\displaystyle{ [\rm A] }[/math] is the concentration at time [math]\displaystyle{ t }[/math] and [math]\displaystyle{ [\rm A]_0 }[/math] is the initial concentration at zero time. The first-order rate law is confirmed if [math]\displaystyle{ \ln{[\ce A]} }[/math] is in fact a linear function of time. In this case the rate constant [math]\displaystyle{ k }[/math] is equal to the slope with sign reversed.^{[10]}^{[11]}
Method of flooding
The partial order with respect to a given reactant can be evaluated by the method of flooding (or of isolation) of Ostwald. In this method, the concentration of one reactant is measured with all other reactants in large excess so that their concentration remains essentially constant. For a reaction a · A + b · B → c · C with rate law [math]\displaystyle{ v_0 = k \cdot [{\rm A}]^x \cdot [{\rm B}]^y, }[/math] the partial order [math]\displaystyle{ x }[/math] with respect to [math]\displaystyle{ \rm A }[/math] is determined using a large excess of [math]\displaystyle{ \rm B }[/math]. In this case
[math]\displaystyle{ v_0 = k' \cdot [{\rm A}]^x }[/math] with [math]\displaystyle{ k' = k \cdot [{\rm B}]^y, }[/math]
and [math]\displaystyle{ x }[/math] may be determined by the integral method. The order [math]\displaystyle{ y }[/math] with respect to [math]\displaystyle{ \rm B }[/math] under the same conditions (with [math]\displaystyle{ \rm B }[/math] in excess) is determined by a series of similar experiments with a range of initial concentration [math]\displaystyle{ [\rm B]_0 }[/math] so that the variation of [math]\displaystyle{ k' }[/math] can be measured.^{[12]}
Zero order
For zero-order reactions, the reaction rate is independent of the concentration of a reactant, so that changing its concentration has no effect on the rate of the reaction. Thus, the concentration changes linearly with time. This may occur when there is a bottleneck which limits the number of reactant molecules that can react at the same time, for example if the reaction requires contact with an enzyme or a catalytic surface.^{[13]}
Many enzyme-catalyzed reactions are zero order, provided that the reactant concentration is much greater than the enzyme concentration which controls the rate, so that the enzyme is saturated. For example, the biological oxidation of ethanol to acetaldehyde by the enzyme liver alcohol dehydrogenase (LADH) is zero order in ethanol.^{[14]}
Similarly reactions with heterogeneous catalysis can be zero order if the catalytic surface is saturated. For example, the decomposition of phosphine (PH
3) on a hot tungsten surface at high pressure is zero order in phosphine, which decomposes at a constant rate.^{[13]}
In homogeneous catalysis zero order behavior can come about from reversible inhibition. For example, ring-opening metathesis polymerization using third-generation Grubbs catalyst exhibits zero order behavior in catalyst due to the reversible inhibition that occurs between pyridine and the ruthenium center.^{[15]}
First order
A first order reaction depends on the concentration of only one reactant (a unimolecular reaction). Other reactants can be present, but their concentration has no effect on the rate. The rate law for a first order reaction is
- [math]\displaystyle{ -\frac{d[\ce A]}{dt} = k[\ce A], }[/math]
Although not affecting the above math, the majority of first order reactions proceed via intermolecular collisions. Such collisions, which contribute the energy to the reactant, are necessarily second order. The rate of these collisions is, however, masked by the fact that the rate determining step remains the unimolecular breakdown of the energized reactant.
The half-life is independent of the starting concentration and is given by [math]\displaystyle{ t_{1/2} = \frac{\ln{(2)}}{k} }[/math].
Examples of such reactions are:
- [math]\ce{ 2N2O5 -> 4NO2 + O2 }[/math] ^{[16]}^{[17]}
- [math]\ce{ [CoCl(NH3)5]^2+ + H2O -> [Co(H2O)(NH3)5]^3+ + Cl- }[/math]
- [math]\ce{ H2O2 -> H2O + 1/2O2 }[/math] ^{[18]}
In organic chemistry, the class of S_{N}1 (nucleophilic substitution unimolecular) reactions consists of first-order reactions. For example, in the reaction of aryldiazonium ions with nucleophiles in aqueous solution, ArN+
2 + X−
→ ArX + N
2, the rate equation is [math]\displaystyle{ v_0 = k[\ce{ArN2+}], }[/math] where Ar indicates an aryl group.^{[19]}
Second order
A reaction is said to be second order when the overall order is two. The rate of a second-order reaction may be proportional to one concentration squared, [math]\displaystyle{ v_0 = k [\ce A]^2, }[/math] or (more commonly) to the product of two concentrations, [math]\displaystyle{ v_0 = k[\ce A][\ce B]. }[/math] As an example of the first type, the reaction NO
2 + CO → NO + CO
2 is second-order in the reactant NO
2 and zero order in the reactant CO. The observed rate is given by [math]\displaystyle{ v_0 = k [\ce{NO2}]^2, }[/math] and is independent of the concentration of CO.^{[20]}
For the rate proportional to a single concentration squared, the time dependence of the concentration is given by
- [math]\displaystyle{ \frac{1}{\ce{[A]}} = \frac{1}{\ce{[A]0}} + kt. }[/math]
The time dependence for a rate proportional to two unequal concentrations is
- [math]\displaystyle{ \frac{\ce{[A]}}{\ce{[B]}} = \frac{\ce{[A]0}}{\ce{[B]0}} e^{\left(\ce{[A]0} - \ce{[B]0}\right) kt}; }[/math]
if the concentrations are equal, they satisfy the previous equation.
The second type includes nucleophilic addition-elimination reactions, such as the alkaline hydrolysis of ethyl acetate:^{[19]}
- [math]\ce{ CH3COOC2H5 + OH- -> CH3COO- + C2H5OH }[/math]
This reaction is first-order in each reactant and second-order overall:
- [math]\displaystyle{ v_0 = k [\ce{CH3COOC2H5}][\ce{OH-}] }[/math]
If the same hydrolysis reaction is catalyzed by imidazole, the rate equation becomes^{[19]}
- [math]\displaystyle{ v_0 = k[\text{imidazole}][\ce{CH3COOC2H5}]. }[/math]
The rate is first-order in one reactant (ethyl acetate), and also first-order in imidazole, which as a catalyst does not appear in the overall chemical equation.
Another well-known class of second-order reactions are the S_{N}2 (bimolecular nucleophilic substitution) reactions, such as the reaction of n-butyl bromide with sodium iodide in acetone:
- [math]\ce{ CH3CH2CH2CH2Br + NaI -> CH3CH2CH2CH2I + NaBr(v) }[/math]
This same compound can be made to undergo a bimolecular (E2) elimination reaction, another common type of second-order reaction, if the sodium iodide and acetone are replaced with sodium tert-butoxide as the salt and tert-butanol as the solvent:
- [math]\ce{ CH3CH2CH2CH2Br + NaO\mathit{t}-Bu -> CH3CH2CH=CH2 + NaBr + HO\mathit{t}-Bu }[/math]
Pseudo-first order
If the concentration of a reactant remains constant (because it is a catalyst, or because it is in great excess with respect to the other reactants), its concentration can be included in the rate constant, leading to a pseudo–first-order (or occasionally pseudo–second-order) rate equation. For a typical second-order reaction with rate equation [math]\displaystyle{ v_0 = k[\ce A][\ce B], }[/math] if the concentration of reactant B is constant then [math]\displaystyle{ v_0 = k[\ce{A}][\ce{B}] = k'[\ce{A}], }[/math] where the pseudo–first-order rate constant [math]\displaystyle{ k' = k[\ce B]. }[/math] The second-order rate equation has been reduced to a pseudo–first-order rate equation, which makes the treatment to obtain an integrated rate equation much easier.
One way to obtain a pseudo-first order reaction is to use a large excess of one reactant (say, [B]≫[A]) so that, as the reaction progresses, only a small fraction of the reactant in excess (B) is consumed, and its concentration can be considered to stay constant. For example, the hydrolysis of esters by dilute mineral acids follows pseudo-first order kinetics, where the concentration of water is constant because it is present in large excess:
- [math]\ce{ CH3COOCH3 + H2O -> CH3COOH + CH3OH }[/math]
The hydrolysis of sucrose (C
12H
22O
11) in acid solution is often cited as a first-order reaction with rate [math]\displaystyle{ v_0 = k[\ce{C12H22O11}]. }[/math] The true rate equation is third-order, [math]\displaystyle{ v_0 = k[\ce{C12H22O11}][\ce{H+}][\ce{H2O}]; }[/math] however, the concentrations of both the catalyst H+
and the solvent H
2O are normally constant, so that the reaction is pseudo–first-order.^{[21]}
Summary for reaction orders 0, 1, 2, and n
Elementary reaction steps with order 3 (called ternary reactions) are rare and unlikely to occur. However, overall reactions composed of several elementary steps can, of course, be of any (including non-integer) order.
Zero order | First order | Second order | nth order (g = 1-n) | |
---|---|---|---|---|
Rate Law | [math]\displaystyle{ -{d[\ce A]}/{dt} = k }[/math] | [math]\displaystyle{ -{d[\ce A]}/{dt} = k[\ce A] }[/math] | [math]\displaystyle{ -{d[\ce A]}/{dt} = k[\ce A]^2 }[/math]^{[22]} | [math]\displaystyle{ -{d[\ce A]}/{dt} = k[\ce A]^n }[/math] |
Integrated Rate Law | [math]\displaystyle{ \ce{[A] = [A]0} - kt }[/math] | [math]\displaystyle{ \ce{[A] = [A]0} e^{-kt} }[/math] | [math]\displaystyle{ \frac{1}\ce{[A]} = \frac{1}\ce{[A]0} + kt }[/math]^{[22]} | [math]\displaystyle{ [\ce A]^{g} = {\ce{[A]0}^{g}} - gkt }[/math]
[Except first order] |
Units of Rate Constant (k) | [math]\displaystyle{ \rm\frac{M}{s} }[/math] | [math]\displaystyle{ \rm\frac{1}{s} }[/math] | [math]\displaystyle{ \rm\frac{1}{M \cdot s} }[/math] | [math]\displaystyle{ \frac{{\rm M}^{g}}{\rm s} }[/math] |
Linear Plot to determine k | [A] vs. t | [math]\ce{ \ln ([A]) }[/math] vs. t | [math]\ce{ \frac{1}{[A]} }[/math] vs. t | [math]\ce{ {\rm [A]}^{g} }[/math] vs. t
[Except first order] |
Half-life | [math]\displaystyle{ t_\frac{1}{2} = \frac\ce{[A]0}{2k} }[/math] | [math]\displaystyle{ t_\frac{1}{2} = \frac{\ln (2)}{k} }[/math] | [math]\displaystyle{ t_\frac{1}{2} = \frac{1}{k\ce{[A]0}} }[/math]^{[22]} | [math]\displaystyle{ t_\frac{1}{2} = \frac{\ce{[A]0}^{g}(1-2^{-g})}{gk} }[/math]
[Limit is necessary for first order] |
Here [math]\displaystyle{ \rm M }[/math] stands for concentration in molarity (mol · L^{−1}), [math]\displaystyle{ t }[/math] for time, and [math]\displaystyle{ k }[/math] for the reaction rate constant. The half-life of a first order reaction is often expressed as t_{1/2} = 0.693/k (as ln(2)≈0.693).
Fractional order
In fractional order reactions, the order is a non-integer, which often indicates a chemical chain reaction or other complex reaction mechanism. For example, the pyrolysis of acetaldehyde (CH
3CHO) into methane and carbon monoxide proceeds with an order of 1.5 with respect to acetaldehyde: [math]\displaystyle{ v_0 = k[\ce{CH3CHO}]^{3/2}. }[/math]^{[23]} The decomposition of phosgene (COCl
2) to carbon monoxide and chlorine has order 1 with respect to phosgene itself and order 0.5 with respect to chlorine: [math]\displaystyle{ v_0 = k \ce{[COCl2] [Cl2]}^{1/2}. }[/math]^{[24]}
The order of a chain reaction can be rationalized using the steady state approximation for the concentration of reactive intermediates such as free radicals. For the pyrolysis of acetaldehyde, the Rice-Herzfeld mechanism is
- Initiation
- [math]\ce{ CH3CHO -> .CH3 + .CHO }[/math]
- Propagation
- [math]\ce{ .CH3 + CH3CHO -> CH3CO. + CH4 }[/math]
- [math]\ce{ CH3CO. -> .CH3 + CO }[/math]
- Termination
- [math]\ce{ 2 .CH3 -> C2H6 }[/math]
where • denotes a free radical.^{[23]}^{[25]} To simplify the theory, the reactions of the *CHO to form a second *CH
3 are ignored.
In the steady state, the rates of formation and destruction of methyl radicals are equal, so that
- [math]\displaystyle{ \frac{d[\ce{.CH3}]}{dt} = k_i[\ce{CH3CHO}]-k_t[\ce{.CH3}]^2 = 0 , }[/math]
so that the concentration of methyl radical satisfies
- [math]\ce{ [.CH3] \quad\propto \quad[CH3CHO]^{1/2}. }[/math]
The reaction rate equals the rate of the propagation steps which form the main reaction products CH
4 and CO:
- [math]\displaystyle{ v_0 = \frac{d[\ce{CH4}]}{dt}|_0 = k_p\ce{[.CH3][CH3CHO]} \quad\propto \quad\ce{[CH3CHO]^{3/2}} }[/math]
in agreement with the experimental order of 3/2.^{[23]}^{[25]}
Complex laws
Mixed order
More complex rate laws have been described as being mixed order if they approximate to the laws for more than one order at different concentrations of the chemical species involved. For example, a rate law of the form [math]\displaystyle{ v_0 = k_1[A]+k_2[A]^2 }[/math] represents concurrent first order and second order reactions (or more often concurrent pseudo-first order and second order) reactions, and can be described as mixed first and second order.^{[26]} For sufficiently large values of [A] such a reaction will approximate second order kinetics, but for smaller [A] the kinetics will approximate first order (or pseudo-first order). As the reaction progresses, the reaction can change from second order to first order as reactant is consumed.
Another type of mixed-order rate law has a denominator of two or more terms, often because the identity of the rate-determining step depends on the values of the concentrations. An example is the oxidation of an alcohol to a ketone by hexacyanoferrate (III) ion [Fe(CN)_{6}^{3−}] with ruthenate (VI) ion (RuO_{4}^{2−}) as catalyst.^{[27]} For this reaction, the rate of disappearance of hexacyanoferrate (III) is [math]\displaystyle{ v_0 = \frac{\ce{[Fe(CN)6]^2-}}{k_\alpha + k_\beta\ce{[Fe(CN)6]^2-}} }[/math]
This is zero-order with respect to hexacyanoferrate (III) at the onset of the reaction (when its concentration is high and the ruthenium catalyst is quickly regenerated), but changes to first-order when its concentration decreases and the regeneration of catalyst becomes rate-determining.
Notable mechanisms with mixed-order rate laws with two-term denominators include:
- Michaelis-Menten kinetics for enzyme-catalysis: first-order in substrate (second-order overall) at low substrate concentrations, zero order in substrate (first-order overall) at higher substrate concentrations; and
- the Lindemann mechanism for unimolecular reactions: second-order at low pressures, first-order at high pressures.
Negative order
A reaction rate can have a negative partial order with respect to a substance. For example, the conversion of ozone (O_{3}) to oxygen follows the rate equation [math]\displaystyle{ v_0 =k\ce{[O_3]^2}\ce{[O_2]^{-1}} }[/math] in an excess of oxygen. This corresponds to second order in ozone and order (−1) with respect to oxygen.^{[28]}
When a partial order is negative, the overall order is usually considered as undefined. In the above example, for instance, the reaction is not described as first order even though the sum of the partial orders is [math]\displaystyle{ 2 + (-1) = 1 }[/math], because the rate equation is more complex than that of a simple first-order reaction.
Opposed reactions
A pair of forward and reverse reactions may occur simultaneously with comparable speeds. For example, A and B react into products P and Q and vice versa (a, b, p, and q are the stoichiometric coefficients):
- [math]\ce{ {\mathit{a}A} + {\mathit{b}B} <=> {\mathit{p}P} + {\mathit{q}Q} }[/math]
The reaction rate expression for the above reactions (assuming each one is elementary) can be written as:
- [math]\displaystyle{ v = k_1 [\ce A]^a[\ce B]^b - k_{-1} [\ce P]^p[\ce Q]^q }[/math]
where: k_{1} is the rate coefficient for the reaction that consumes A and B; k_{−1} is the rate coefficient for the backwards reaction, which consumes P and Q and produces A and B.
The constants k_{1} and k_{−1} are related to the equilibrium coefficient for the reaction (K) by the following relationship (set v=0 in balance):
- [math]\displaystyle{ k_1 [\ce A]^a[\ce B]^b = k_{-1} [\ce P]^p[\ce Q]^q }[/math]
- [math]\displaystyle{ K = \frac{[\ce P]^p[\ce Q]^q}{[\ce A]^a[\ce B]^b} = \frac{k_1}{k_{-1}} }[/math]
Simple example
In a simple equilibrium between two species:
- [math]\ce{ A <=> P }[/math]
where the reaction starts with an initial concentration of reactant A, [math]\ce{ [A]0 }[/math], and an initial concentration of 0 for product P at time t=0.
Then the equilibrium constant K is expressed as:
- [math]\displaystyle{ K \ \stackrel{\mathrm{def}}{=}\ \frac{k_{1}}{k_{-1}} = \frac{\left[\ce P\right]_e} {\left[\ce A\right]_e} }[/math]
where [math]\displaystyle{ [\ce A]_e }[/math] and [math]\displaystyle{ [\ce P]_e }[/math] are the concentrations of A and P at equilibrium, respectively.
The concentration of A at time t, [math]\displaystyle{ [\ce A]_t }[/math], is related to the concentration of P at time t, [math]\displaystyle{ [\ce P]_t }[/math], by the equilibrium reaction equation:
- [math]\ce{ [A]_\mathit{t} = [A]0 - [P]_\mathit{t} }[/math]
The term [math]\ce{ [P]0 }[/math] is not present because, in this simple example, the initial concentration of P is 0.
This applies even when time t is at infinity; i.e., equilibrium has been reached:
- [math]\ce{ [A]_\mathit{e} = [A]0 - [P]_\mathit{e} }[/math]
then it follows, by the definition of K, that
- [math]\displaystyle{ [\ce P]_e = \frac{k_{1}}{k_1+k_{-1}}\ce{[A]0} }[/math]
and, therefore,
- [math]\displaystyle{ \ [\ce A]_e = \ce{[A]0} - [\ce P]_e = \frac{k_{-1}}{k_1+k_{-1}}\ce{[A]0} }[/math]
These equations allow us to uncouple the system of differential equations, and allow us to solve for the concentration of A alone.
The reaction equation was given previously as:
- [math]\displaystyle{ v = k_1 [\ce A]^a[\ce B]^b - k_{-1} [\ce P]^p[\ce Q]^q }[/math]
For [math]\ce{ A <=> P }[/math] this is simply
- [math]\displaystyle{ -\frac{d[\ce A]}{dt} = k_1 [\ce A]_t - k_{-1} [\ce P]_t }[/math]
The derivative is negative because this is the rate of the reaction going from A to P, and therefore the concentration of A is decreasing. To simplify notation, let x be [math]\displaystyle{ [\ce A]_t }[/math], the concentration of A at time t. Let [math]\displaystyle{ x_e }[/math] be the concentration of A at equilibrium. Then:
- [math]\displaystyle{ \begin{align} -\frac{d[\ce A]}{dt} &= {k_1 [\ce A]_t} - {k_{-1} [\ce P]_t}\\ -\frac{dx}{dt} &= {k_1 x} - {k_{-1} [\ce P]_t}\\ &= {k_1 x} - {k_{-1} (\ce{[A]0} - x)}\\ &= {(k_1 + k_{-1})x} - {k_{-1} \ce{[A]0}} \end{align} }[/math]
Since:
- [math]\displaystyle{ k_1 + k_{-1} = k_{-1} \frac{\ce{[A]0}}{x_e} }[/math]
the reaction rate becomes:
- [math]\displaystyle{ \frac{dx}{dt} = \frac{k_{-1}\ce{[A]0}}{x_e} (x_e - x) }[/math]
which results in:
- [math]\displaystyle{ \ln \left(\frac{\ce{[A]0} - [\ce A]_e}{[\ce A]_t-[\ce A]_e}\right) = (k_1 + k_{-1})t }[/math].
A plot of the negative natural logarithm of the concentration of A in time minus the concentration at equilibrium versus time t gives a straight line with slope k_{1} + k_{−1}. By measurement of [A]_{e} and [P]_{e} the values of K and the two reaction rate constants will be known.^{[29]}
Generalization of simple example
If the concentration at the time t = 0 is different from above, the simplifications above are invalid, and a system of differential equations must be solved. However, this system can also be solved exactly to yield the following generalized expressions:
- [math]\displaystyle{ \left[ \ce A \right]=\ce{[A]0}\frac{1}{k_{1}+k_{-1}}\left( k_{-1}+k_{1}e^{-\left( k_{1}+k_{-1} \right)t} \right)+\ce{[P]0}\frac{k_{-1}}{k_{1}+k_{-1}}\left( 1-e^{-\left( k_{1}+k_{-1} \right)t} \right) }[/math]
- [math]\displaystyle{ \left[ \ce P \right]=\ce{[A]0}\frac{k_{1}}{k_{1}+k_{-1}}\left( 1-e^{-\left( k_{1}+k_{-1} \right)t} \right)+\ce{[P]0}\frac{1}{k_{1}+k_{-1}}\left( k_{1}+k_{-1}e^{-\left( k_{1}+k_{-1} \right)t} \right) }[/math]
When the equilibrium constant is close to unity and the reaction rates very fast for instance in conformational analysis of molecules, other methods are required for the determination of rate constants for instance by complete lineshape analysis in NMR spectroscopy.
Consecutive reactions
If the rate constants for the following reaction are [math]\displaystyle{ k_1 }[/math] and [math]\displaystyle{ k_2 }[/math]; [math]\ce{ A -> B -> C }[/math], then the rate equation is:
- For reactant A: [math]\displaystyle{ \frac{d[\ce A]}{dt} = -k_1 [\ce A] }[/math]
- For reactant B: [math]\displaystyle{ \frac{d[\ce B]}{dt} = k_1 [\ce A] - k_2 [\ce B] }[/math]
- For product C: [math]\displaystyle{ \frac{d[\ce C]}{dt} = k_2 [\ce B] }[/math]
With the individual concentrations scaled by the total population of reactants to become probabilities, linear systems of differential equations such as these can be formulated as a master equation. The differential equations can be solved analytically and the integrated rate equations are
- [math]\displaystyle{ [\ce A]=\ce{[A]0} e^{-k_1 t} }[/math]
- [math]\displaystyle{ \left[ \ce B \right]= \begin{cases} \ce{[A]0}\frac{k_{1}}{k_{2}-k_{1}}\left( e^{-k_{1}t}-e^{-k_{2}t} \right)+\ce{[B]0}e^{-k_{2}t}& k_{1}\ne k_{2} \\ \ce{[A]0}k_{1}te^{-k_{1}t}+\ce{[B]0}e^{-k_{1}t} & \text{otherwise} \\ \end{cases} }[/math]
- [math]\displaystyle{ \left[ \ce C \right]= \begin{cases} \ce{[A]0}\left( 1+\frac{k_{1}e^{-k_{2}t}-k_{2}e^{-k_{1}t}}{k_{2}-k_{1}} \right)+\ce{[B]0}\left( 1-e^{-k_{2}t} \right)+\ce{[C]0} & k_{1}\ne k_{2} \\ \ce{[A]0}\left( 1-e^{-k_{1}t}-k_{1}te^{-k_{1}t} \right)+\ce{[B]0}\left( 1-e^{-k_{1}t} \right)+\ce{[C]0} & \text{otherwise} \\ \end{cases} }[/math]
The steady state approximation leads to very similar results in an easier way.
Parallel or competitive reactions
When a substance reacts simultaneously to give two different products, a parallel or competitive reaction is said to take place.
Two first order reactions
[math]\ce{ A -> B }[/math] and [math]\ce{ A -> C }[/math], with constants [math]\displaystyle{ k_1 }[/math] and [math]\displaystyle{ k_2 }[/math] and rate equations [math]\displaystyle{ -\frac{d[\ce A]}{dt}=(k_1+k_2)[\ce A] }[/math]; [math]\displaystyle{ \frac{d[\ce B]}{dt}=k_1[\ce A] }[/math] and [math]\displaystyle{ \frac{d[\ce C]}{dt}=k_2[\ce A] }[/math]
The integrated rate equations are then [math]\displaystyle{ [\ce A] = \ce{[A]0} e^{-(k_1+k_2)t} }[/math]; [math]\displaystyle{ [\ce B] = \frac{k_1}{k_1+k_2}\ce{[A]0} \left(1 - e^{-(k_1+k_2)t}\right) }[/math] and [math]\displaystyle{ [\ce C] = \frac{k_2}{k_1+k_2}\ce{[A]0} \left(1 - e^{-(k_1+k_2)t}\right) }[/math].
One important relationship in this case is [math]\displaystyle{ \frac{\ce{[B]}}{\ce{[C]}}=\frac{k_1}{k_2} }[/math]
One first order and one second order reaction
This can be the case when studying a bimolecular reaction and a simultaneous hydrolysis (which can be treated as pseudo order one) takes place: the hydrolysis complicates the study of the reaction kinetics, because some reactant is being "spent" in a parallel reaction. For example, A reacts with R to give our product C, but meanwhile the hydrolysis reaction takes away an amount of A to give B, a byproduct: [math]\ce{ A + H2O -> B }[/math] and [math]\ce{ A + R -> C }[/math]. The rate equations are: [math]\displaystyle{ \frac{d[\ce B]}{dt}=k_1\ce{[A][H2O]}=k_1'[\ce A] }[/math] and [math]\displaystyle{ \frac{d[\ce C]}{dt}=k_2\ce{[A][R]} }[/math], where [math]\displaystyle{ k_1' }[/math] is the pseudo first order constant.^{[30]}
The integrated rate equation for the main product [C] is [math]\displaystyle{ \ce{[C]=[R]0} \left [ 1-e^{-\frac{k_2}{k_1'}\ce{[A]0} \left(1-e^{-k_1't}\right)} \right ] }[/math], which is equivalent to [math]\displaystyle{ \ln \frac{\ce{[R]0}}{\ce{[R]0-[C]}}=\frac{k_2\ce{[A]0}}{k_1'}\left(1 - e^{-k_1't}\right) }[/math]. Concentration of B is related to that of C through [math]\displaystyle{ [\ce B]=-\frac{k_1'}{k_2} \ln \left ( 1 - \frac{\ce{[C]}}{\ce{[R]0}} \right ) }[/math]
The integrated equations were analytically obtained but during the process it was assumed that [math]\displaystyle{ \ce{[A]0} - \ce{[C]}\approx \ce{[A]0} }[/math]. Therefore, previous equation for [C] can only be used for low concentrations of [C] compared to [A]_{0}
Stoichiometric reaction networks
The most general description of a chemical reaction network considers a number [math]\displaystyle{ N }[/math] of distinct chemical species reacting via [math]\displaystyle{ R }[/math] reactions.^{[31]} ^{[32]} The chemical equation of the [math]\displaystyle{ j }[/math]-th reaction can then be written in the generic form
- [math]\displaystyle{ s_{1j} \ce X_1 + s_{2j} \ce X_2 + \cdots + s_{Nj} \ce X_{N} \ce{-\gt [k_j]} \ r_{1j} \ce X_{1} + \ r_{2j} \ce X_{2} + \cdots + r_{Nj} \ce X_{N}, }[/math]
which is often written in the equivalent form
- [math]\displaystyle{ \sum_{i=1}^{N} s_{ij} \ce X_i \ce{-\gt [k_j]} \sum_{i=1}^{N}\ r_{ij} \ce X_{i}. }[/math]
Here
- [math]\displaystyle{ j }[/math] is the reaction index running from 1 to [math]\displaystyle{ R }[/math],
- [math]\displaystyle{ \ce X_i }[/math] denotes the [math]\displaystyle{ i }[/math]-th chemical species,
- [math]\displaystyle{ k_j }[/math] is the rate constant of the [math]\displaystyle{ j }[/math]-th reaction and
- [math]\displaystyle{ s_{ij} }[/math] and [math]\displaystyle{ r_{ij} }[/math] are the stoichiometric coefficients of reactants and products, respectively.
The rate of such a reaction can be inferred by the law of mass action
- [math]\displaystyle{ f_j([\vec{\ce X}])= k_j \prod_{z=1}^N [\ce X_z]^{s_{zj}} }[/math]
which denotes the flux of molecules per unit time and unit volume. Here [math]\ce{ [\vec{X}]=([X1], [X2], \ldots ,[X_\mathit{N}]) }[/math] is the vector of concentrations. This definition includes the elementary reactions:
- zero order reactions
- for which [math]\displaystyle{ s_{zj}=0 }[/math] for all [math]\displaystyle{ z }[/math],
- first order reactions
- for which [math]\displaystyle{ s_{zj}=1 }[/math] for a single [math]\displaystyle{ z }[/math],
- second order reactions
- for which [math]\displaystyle{ s_{zj}=1 }[/math] for exactly two [math]\displaystyle{ z }[/math]; that is, a bimolecular reaction, or [math]\displaystyle{ s_{zj} = 2 }[/math] for a single [math]\displaystyle{ z }[/math]; that is, a dimerization reaction.
Each of these is discussed in detail below. One can define the stoichiometric matrix
- [math]\displaystyle{ S_{ij}=r_{ij}-s_{ij}, }[/math]
denoting the net extent of molecules of [math]\displaystyle{ i }[/math] in reaction [math]\displaystyle{ j }[/math]. The reaction rate equations can then be written in the general form
- [math]\displaystyle{ \frac{d [\ce X_i]}{dt} =\sum_{j=1}^{R} S_{ij} f_j([\vec{\ce X}]). }[/math]
This is the product of the stoichiometric matrix and the vector of reaction rate functions. Particular simple solutions exist in equilibrium, [math]\displaystyle{ \frac{d [\ce X_i]}{dt}=0 }[/math], for systems composed of merely reversible reactions. In this case the rate of the forward and backward reactions are equal, a principle called detailed balance. Detailed balance is a property of the stoichiometric matrix [math]\displaystyle{ S_{ij} }[/math] alone and does not depend on the particular form of the rate functions [math]\displaystyle{ f_j }[/math]. All other cases where detailed balance is violated are commonly studied by flux balance analysis, which has been developed to understand metabolic pathways.^{[33]}^{[34]}
General dynamics of unimolecular conversion
For a general unimolecular reaction involving interconversion of [math]\displaystyle{ N }[/math] different species, whose concentrations at time [math]\displaystyle{ t }[/math] are denoted by [math]\displaystyle{ X_1(t) }[/math] through [math]\displaystyle{ X_N(t) }[/math], an analytic form for the time-evolution of the species can be found. Let the rate constant of conversion from species [math]\displaystyle{ X_i }[/math] to species [math]\displaystyle{ X_j }[/math] be denoted as [math]\displaystyle{ k_{ij} }[/math], and construct a rate-constant matrix [math]\displaystyle{ K }[/math] whose entries are the [math]\displaystyle{ k_{ij} }[/math].
Also, let [math]\displaystyle{ X(t) = (X_1(t), X_2(t), \ldots,X_N(t))^T }[/math] be the vector of concentrations as a function of time.
Let [math]\displaystyle{ J=(1,1,1,\ldots,1)^T }[/math] be the vector of ones.
Let [math]\displaystyle{ I }[/math] be the [math]\displaystyle{ N \times N }[/math] identity matrix.
Let [math]\displaystyle{ \operatorname{Diag} }[/math] be the function that takes a vector and constructs a diagonal matrix whose on-diagonal entries are those of the vector.
Let [math]\displaystyle{ \mathcal{L}^{-1} }[/math] be the inverse Laplace transform from [math]\displaystyle{ s }[/math] to [math]\displaystyle{ t }[/math].
Then the time-evolved state [math]\displaystyle{ X(t) }[/math] is given by
- [math]\displaystyle{ X(t) = \mathcal{L}^{-1}[(sI + \operatorname{Diag}(KJ)-K^T)^{-1}X(0)], }[/math]
thus providing the relation between the initial conditions of the system and its state at time [math]\displaystyle{ t }[/math].
See also
- Michaelis–Menten kinetics
- Molecularity
- Petersen matrix
- Reaction–diffusion system
- Reactions on surfaces: rate equations for reactions where at least one of the reactants adsorbs onto a surface
- Reaction progress kinetic analysis
- Reaction rate
- Reaction rate constant
- Steady state approximation
- Gillespie algorithm
- Balance equation
- Belousov-Zhabotinsky reaction
- Lotka–Volterra equations
References
- ↑ IUPAC Gold Book definition of rate law. See also: According to IUPAC Compendium of Chemical Terminology.
- ↑ Atkins & de Paula 2006, p. 794
- ↑ IUPAC, Compendium of Chemical Terminology, 2nd ed. (the "Gold Book") (1997). Online corrected version: (2006–) "Rate of reaction". doi:10.1351/goldbook.R05156
- ↑ Atkins & de Paula 2006, p. 795
- ↑ ^{5.0} ^{5.1} ^{5.2} Atkins & de Paula 2006, p. 796
- ↑ Connors 1990, p. 13
- ↑ Connors 1990, p. 12
- ↑ Atkins & de Paula 2006, pp. 797–8
- ↑ Espenson 1987, pp. 5–8
- ↑ Atkins & de Paula 2006, pp. 798–800
- ↑ Espenson 1987, pp. 15–18
- ↑ Espenson 1987, pp. 30–31
- ↑ ^{13.0} ^{13.1} Atkins & de Paula 2006, p. 796
- ↑ Tinoco Jr. & Wang 1995, p. 331
- ↑ Walsh, Dylan J.; Lau, Sii Hong; Hyatt, Michael G.; Guironnet, Damien (2017-09-25). "Kinetic Study of Living Ring-Opening Metathesis Polymerization with Third-Generation Grubbs Catalysts" (in EN). Journal of the American Chemical Society 139 (39): 13644–13647. doi:10.1021/jacs.7b08010. ISSN 0002-7863. PMID 28944665.
- ↑ Atkins & de Paula 2006, p. 813-4
- ↑ Keith J. Laidler, Chemical Kinetics (3rd ed., Harper & Row 1987), p.303-5 ISBN:0-06-043862-2
- ↑ R.H. Petrucci, W.S. Harwood and F.G. Herring, General Chemistry (8th ed., Prentice-Hall 2002) p.588 ISBN:0-13-014329-4
- ↑ ^{19.0} ^{19.1} ^{19.2} Connors 1990
- ↑ Whitten K. W., Galley K. D. and Davis R. E. General Chemistry (4th edition, Saunders 1992), pp. 638–9 ISBN:0-03-072373-6
- ↑ Tinoco Jr. & Wang 1995, pp. 328–9
- ↑ ^{22.0} ^{22.1} ^{22.2} NDRL Radiation Chemistry Data Center. See also: Capellos, Christos; Bielski, Benon H. (1972). Kinetic systems: mathematical description of chemical kinetics in solution. New York: Wiley-Interscience. ISBN 978-0471134503. OCLC 247275.
- ↑ ^{23.0} ^{23.1} ^{23.2} Atkins & de Paula 2006, p. 830
- ↑ Laidler 1987, p. 301
- ↑ ^{25.0} ^{25.1} Laidler 1987, pp. 310–311
- ↑ Espenson 1987, pp. 34,60
- ↑ Mucientes, Antonio E.; de la Peña, María A. (November 2006). "Ruthenium(VI)-Catalyzed Oxidation of Alcohols by Hexacyanoferrate(III): An Example of Mixed Order" (in EN). Journal of Chemical Education 83 (11): 1643. doi:10.1021/ed083p1643. ISSN 0021-9584. http://jchemed.chem.wisc.edu/Journal/Issues/2006/Nov/abs1643.html.
- ↑ Laidler 1987, p. 305
- ↑ Rushton, Gregory T.; Burns, William G.; Lavin, Judi M.; Chong, Yong S.; Pellechia, Perry; Shimizu, Ken D. (September 2007). "Determination of the Rotational Barrier for Kinetically Stable Conformational Isomers via NMR and 2D TLC" (in EN). Journal of Chemical Education 84 (9): 1499. doi:10.1021/ed084p1499. ISSN 0021-9584. http://jchemed.chem.wisc.edu/Journal/Issues/2007/Sep/abs1499.html.
- ↑ Manso, José A.; Pérez-Prior, M. Teresa; García-Santos, M. del Pilar; Calle, Emilio; Casado, Julio (2005). "A Kinetic Approach to the Alkylating Potential of Carcinogenic Lactones". Chemical Research in Toxicology 18 (7): 1161–1166. doi:10.1021/tx050031d. PMID 16022509.
- ↑ Heinrich, Reinhart; Schuster, Stefan (2012). The Regulation of Cellular Systems. Springer Science & Business Media. ISBN 9781461311614.
- ↑ Chen, Luonan; Wang, Ruiqi; Li, Chunguang; Aihara, Kazuyuki (2010). Modeling Biomolecular Networks in Cells. doi:10.1007/978-1-84996-214-8. ISBN 978-1-84996-213-1.
- ↑ Szallasi, Z. and Stelling, J. and Periwal, V. (2006) System modeling in cell biology: from concepts to nuts and bolts. MIT Press Cambridge.
- ↑ Iglesias, Pablo A.; Ingalls, Brian P. (2010). Control theory and systems biology. MIT Press. ISBN 9780262013345.
Books cited
- Atkins, Peter; de Paula, Julio (2006). "The rates of chemical reactions". Atkins' Physical chemistry (8th ed.). W.H. Freeman. pp. 791–823. ISBN 0-7167-8759-8. https://archive.org/details/atkinsphysicalch00pwat/page/791.
- Connors, Kenneth Antonio (1990). Chemical kinetics : the study of reaction rates in solution. John Wiley & Sons. ISBN 9781560810063.
- Espenson, James H. (1987). Chemical kinetics and reaction mechanisms (2nd ed.). McGraw Hill. ISBN 9780071139496.
- Laidler, Keith James (1987). Chemical kinetics (3rd ed.). Harper & Row. ISBN 9780060438623.
- Tinoco Jr., Ignacio; Wang, James C. (1995). Physical chemistry : principles and applications in biological sciences (3rd ed.). Prentice Hall. ISBN 9780131865457.
External links
- Chemical kinetics, reaction rate, and order (needs flash player)
- Reaction kinetics, examples of important rate laws (lecture with audio).
- Rates of Reaction
cy:Cyfradd adwaith#Hafaliadau cyfradd
Original source: https://en.wikipedia.org/wiki/Rate equation.
Read more |