Collision theory

Reaction rate tends to increase with concentration phenomenon explained by collision theory

Collision theory is a principle of chemistry used to predict the rates of chemical reactions. It states that when suitable particles of the reactant hit each other with the correct orientation, only a certain amount of collisions result in a perceptible or notable change; these successful changes are called successful collisions. The successful collisions must have enough energy, also known as activation energy, at the moment of impact to break the pre-existing bonds and form all new bonds. This results in the products of the reaction. The activation energy is often predicted using the Transition state theory. Increasing the concentration of the reactant brings about more collisions and hence more successful collisions. Increasing the temperature increases the average kinetic energy of the molecules in a solution, increasing the number of collisions that have enough energy. Collision theory was proposed independently by Max Trautz in 1916^[1] and William Lewis in 1918.^{[2] [3]}

When a catalyst is involved in the collision between the reactant molecules, less energy is required for the chemical change to take place, and hence more collisions have sufficient energy for the reaction to occur. The reaction rate therefore increases.

Collision theory is closely related to chemical kinetics.

Collision theory was initially developed for the gas reaction system with no dilution. But most reactions involve solutions, for example, gas reactions in a carrying inert gas, and almost all reactions in solutions. The collision frequency of the solute molecules in these solutions is now controlled by diffusion or Brownian motion of individual molecules. The flux of the diffusive molecules follows Fick's laws of diffusion. For particles in a solution, an example model to calculate the collision frequency and associated coagulation rate is the Smoluchowski coagulation equation proposed by Marian Smoluchowski in a seminal 1916 publication.^[4] In this model, Fick's flux at the infinite time limit is used to mimic the particle speed of the collision theory. Jixin Chen proposed a finite-time solution to the diffusion flux in 2022 which significantly changes the estimated collision frequency of two particles in a solution.^[5]

Rate equations

The rate for a bimolecular gas-phase reaction, A + B → product, predicted by collision theory is^[6]

${\displaystyle r(T)=kn_{\text{A}}n_{\text{B}}=Z\rho \exp \left({\frac {-E_{\text{a}}}{RT}}\right)}$

where:

• k is the rate constant in units of (number of molecules)⁻¹⋅s⁻¹⋅m³.
• n_A is the number density of A in the gas in units of m⁻³.
• n_B is the number density of B in the gas in units of m⁻³. E.g. for a gas mixture with gas A concentration 0.1 mol⋅L⁻¹ and B concentration 0.2 mol⋅L⁻¹, the number of density of A is 0.1×6.02×10²³÷10⁻³ = 6.02×10²⁵ m⁻³, the number of density of B is 0.2×6.02×10²³÷10⁻³ = 1.2×10²⁶ m⁻³
• Z is the collision frequency in units of m⁻³⋅s⁻¹.
• ${\displaystyle \rho }$ is the steric factor.^[7]
• E_a is the activation energy of the reaction, in units of J⋅mol⁻¹.
• T is the temperature in units of K.
• R is the gas constant in units of J mol⁻¹K⁻¹.

The unit of r(T) can be converted to mol⋅L⁻¹⋅s⁻¹, after divided by (1000×N_A), where N_A is the Avogadro constant.

For a reaction between A and B, the collision frequency calculated with the hard-sphere model with the unit number of collisions per m³ per second is:

${\displaystyle Z=n_{\text{A}}n_{\text{B}}\sigma _{\text{AB}}{\sqrt {\frac {8k_{\text{B}}T}{\pi \mu _{\text{AB}}}}}=10^{6}N_{A}^{2}{\text{[A][B]}}\sigma _{\text{AB}}{\sqrt {\frac {8k_{\text{B}}T}{\pi \mu _{\text{AB}}}}}}$

where:

• n_A is the number density of A in the gas in units of m⁻³.
• n_B is the number density of B in the gas in units of m⁻³. E.g. for a gas mixture with gas A concentration 0.1 mol⋅L⁻¹ and B concentration 0.2 mol⋅L⁻¹, the number of density of A is 0.1×6.02×10²³÷10⁻³ = 6.02×10²⁵ m⁻³, the number of density of B is 0.2×6.02×10²³÷10⁻³ = 1.2×10²⁶ m⁻³.
• σ_AB is the reaction cross section (unit m²), the area when two molecules collide with each other, simplified to ${\displaystyle \sigma _{\text{AB}}=\pi (r_{\text{A}}+r_{\text{B}})^{2}}$, where r_A the radius of A and r_B the radius of B in unit m.
• k_B is the Boltzmann constant unit J⋅K⁻¹.
• T is the absolute temperature (unit K).
• μ_AB is the reduced mass of the reactants A and B, ${\displaystyle \mu _{\text{AB}}={\frac {{m_{\text{A}}}{m_{\text{B}}}}{{m_{\text{A}}}+{m_{\text{B}}}}}}$ (unit kg).
• N_A is the Avogadro constant.
• [A] is molar concentration of A in unit mol⋅L⁻¹.
• [B] is molar concentration of B in unit mol⋅L⁻¹.
• Z can be converted to mole collision per liter per second dividing by 1000N_A.

If all the units that are related to dimension are converted to dm, i.e. mol⋅dm⁻³ for [A] and [B], dm² for σ_AB, dm²⋅kg⋅s⁻²⋅K⁻¹ for the Boltzmann constant, then

${\displaystyle Z=N_{\text{A}}\sigma _{\text{AB}}{\sqrt {\frac {8k_{\text{B}}T}{\pi \mu _{\text{AB}}}}}[{\text{A}}][{\text{B}}]=k[A][B]}$

unit mol⋅dm⁻³⋅s⁻¹.

Quantitative insights

Derivation

Consider the bimolecular elementary reaction:

A + B → C

In collision theory it is considered that two particles A and B will collide if their nuclei get closer than a certain distance. The area around a molecule A in which it can collide with an approaching B molecule is called the cross section (σ_AB) of the reaction and is, in simplified terms, the area corresponding to a circle whose radius (${\displaystyle r_{AB}}$) is the sum of the radii of both reacting molecules, which are supposed to be spherical. A moving molecule will therefore sweep a volume ${\displaystyle \pi r_{AB}^{2}c_{A}}$ per second as it moves, where ${\displaystyle c_{A}}$ is the average velocity of the particle. (This solely represents the classical notion of a collision of solid balls. As molecules are quantum-mechanical many-particle systems of electrons and nuclei based upon the Coulomb and exchange interactions, generally they neither obey rotational symmetry nor do they have a box potential. Therefore, more generally the cross section is defined as the reaction probability of a ray of A particles per areal density of B targets, which makes the definition independent from the nature of the interaction between A and B. Consequently, the radius ${\displaystyle r_{AB}}$ is related to the length scale of their interaction potential.)

From kinetic theory it is known that a molecule of A has an average velocity (different from root mean square velocity) of ${\displaystyle c_{A}={\sqrt {\frac {8k_{\text{B}}T}{\pi m_{A}}}}}$, where ${\displaystyle k_{\text{B}}}$ is the Boltzmann constant, and ${\displaystyle m_{A}}$ is the mass of the molecule.

The solution of the two-body problem states that two different moving bodies can be treated as one body which has the reduced mass of both and moves with the velocity of the center of mass, so, in this system ${\displaystyle \mu _{AB}}$ must be used instead of ${\displaystyle m_{A}}$. Thus, for a given molecule A, it travels ${\displaystyle t=l/c_{A}=1/(n_{B}\sigma _{AB}c_{A})}$ before hitting a molecule B if all B is fixed with no movement, where ${\displaystyle l}$ is the average traveling distance. Since B also moves, the relative velocity can be calculated using the reduced mass of A and B.

Therefore, the total collision frequency,^[8] of all A molecules, with all B molecules, is

${\displaystyle Z=n_{\text{A}}n_{\text{B}}\sigma _{AB}{\sqrt {\frac {8k_{\text{B}}T}{\pi \mu _{AB}}}}=10^{6}N_{A}^{2}[A][B]\sigma _{AB}{\sqrt {\frac {8k_{\text{B}}T}{\pi \mu _{AB}}}}=z[A][B],}$

From Maxwell–Boltzmann distribution it can be deduced that the fraction of collisions with more energy than the activation energy is ${\displaystyle e^{\frac {-E_{\text{a}}}{RT}}}$. Therefore, the rate of a bimolecular reaction for ideal gases will be

${\displaystyle r=z\rho [A][B]\exp \left({\frac {-E_{\text{a}}}{RT}}\right),}$ in unit number of molecular reactions s⁻¹⋅m⁻³,

where:

• Z is the collision frequency with unit s⁻¹⋅m⁻³. The z is Z without [A][B].
• ${\displaystyle \rho }$ is the steric factor, which will be discussed in detail in the next section,
• E_a is the activation energy (per mole) of the reaction in unit J/mol,
• T is the absolute temperature in unit K,
• R is the gas constant in unit J/mol/K.
• [A] is molar concentration of A in unit mol/L,
• [B] is molar concentration of B in unit mol/L.

The product zρ is equivalent to the preexponential factor of the Arrhenius equation.

Validity of the theory and steric factor

Once a theory is formulated, its validity must be tested, that is, compare its predictions with the results of the experiments.

When the expression form of the rate constant is compared with the rate equation for an elementary bimolecular reaction, ${\displaystyle r=k(T)[A][B]}$, it is noticed that

${\displaystyle k(T)=N_{A}\sigma _{AB}\rho {\sqrt {\frac {8k_{\text{B}}T}{\pi \mu _{AB}}}}\exp \left({\frac {-E_{\text{a}}}{RT}}\right)}$

unit M⁻¹⋅s⁻¹ (= dm³⋅mol⁻¹⋅s⁻¹), with all dimension unit dm including k_B.

This expression is similar to the Arrhenius equation and gives the first theoretical explanation for the Arrhenius equation on a molecular basis. The weak temperature dependence of the preexponential factor is so small compared to the exponential factor that it cannot be measured experimentally, that is, "it is not feasible to establish, on the basis of temperature studies of the rate constant, whether the predicted T^⁠1/2⁠ dependence of the preexponential factor is observed experimentally".^[9]

Steric factor

If the values of the predicted rate constants are compared with the values of known rate constants, it is noticed that collision theory fails to estimate the constants correctly, and the more complex the molecules are, the more it fails. The reason for this is that particles have been supposed to be spherical and able to react in all directions, which is not true, as the orientation of the collisions is not always proper for the reaction. For example, in the hydrogenation reaction of ethylene the H₂ molecule must approach the bonding zone between the atoms, and only a few of all the possible collisions fulfill this requirement.

To alleviate this problem, a new concept must be introduced: the steric factor ρ. It is defined as the ratio between the experimental value and the predicted one (or the ratio between the frequency factor and the collision frequency):

${\displaystyle \rho ={\frac {A_{\text{observed}}}{Z_{\text{calculated}}}},}$

and it is most often less than unity.^[7]

Usually, the more complex the reactant molecules, the lower the steric factor. Nevertheless, some reactions exhibit steric factors greater than unity: the harpoon reactions, which involve atoms that exchange electrons, producing ions. The deviation from unity can have different causes: the molecules are not spherical, so different geometries are possible; not all the kinetic energy is delivered into the right spot; the presence of a solvent (when applied to solutions), etc.

Experimental rate constants compared to the ones predicted by collision theory for gas phase reactions

Reaction	A, s^−1M^−1	Z, s^−1M^−1	Steric factor
2ClNO → 2Cl + 2NO	9.4×10^9	5.9×10^10	0.16
2ClO → Cl2 + O2	6.3×10^7	2.5×10^10	2.3×10^−3
H2 + C2H4 → C2H6	1.24×10^6	7.3×10^11	1.7×10^−6
Br2 + K → KBr + Br	1.0×10^12	2.1×10^11	4.3

Collision theory can be applied to reactions in solution; in that case, the solvent cage has an effect on the reactant molecules, and several collisions can take place in a single encounter, which leads to predicted preexponential factors being too large. ρ values greater than unity can be attributed to favorable entropic contributions.

Experimental rate constants compared to the ones predicted by collision theory for reactions in solution^[10]

Reaction	Solvent	A, 10^11 s^−1⋅M^−1	Z, 10^11 s^−1⋅M^−1	Steric factor
C2H5Br + OH^−	ethanol	4.30	3.86	1.11
C2H5O^− + CH3I	ethanol	2.42	1.93	1.25
ClCH2CO2^− + OH^−	water	4.55	2.86	1.59
C3H6Br2 + I^−	methanol	1.07	1.39	0.77
HOCH2CH2Cl + OH^−	water	25.5	2.78	9.17
4-CH3C6H4O^− + CH3I	ethanol	8.49	1.99	4.27
CH3(CH2)2Cl + I^−	acetone	0.085	1.57	0.054
C5H5N + CH3I	C2H2Cl4	—	—	2.0 10×10^−6

Alternative collision models for diluted solutions

Collision in diluted gas or liquid solution is regulated by diffusion instead of direct collisions, which can be calculated from Fick's laws of diffusion. Theoretical models to calculate the collision frequency in solutions have been proposed by Marian Smoluchowski in a seminal 1916 publication at the infinite time limit,^[4] and Jixin Chen in 2022 at a finite-time approximation.^[5] A scheme of comparing the rate equations in pure gas and solution is shown in the right figure.

A scheme comparing direct collision and diffusive collision, with corresponding rate equations.

For a diluted solution in the gas or the liquid phase, the collision equation developed for neat gas is not suitable when diffusion takes control of the collision frequency, i.e., the direct collision between the two molecules no longer dominates. For any given molecule A, it has to collide with a lot of solvent molecules, let's say molecule C, before finding the B molecule to react with. Thus the probability of collision should be calculated using the Brownian motion model, which can be approximated to a diffusive flux using various boundary conditions that yield different equations in the Smoluchowski model and the JChen Model.

For the diffusive collision, at the infinite time limit when the molecular flux can be calculated from the Fick's laws of diffusion, in 1916 Smoluchowski derived a collision frequency between molecule A and B in a diluted solution:^[4]

${\displaystyle Z_{AB}=4\pi RD_{r}C_{A}C_{B}}$

where:

• ${\displaystyle Z_{AB}}$ is the collision frequency, unit #collisions/s in 1 m³ of solution.
• ${\displaystyle R}$ is the radius of the collision cross-section, unit m.
• ${\displaystyle D_{r}}$ is the relative diffusion constant between A and B, unit m²/s, and ${\displaystyle D_{r}=D_{A}+D_{B}}$.
• ${\displaystyle C_{A}}$ and ${\displaystyle C_{B}}$ are the number concentrations of molecules A and B in the solution respectively, unit #molecule/m³.

or

${\displaystyle Z_{AB}=1000N_{A}*4\pi RD_{r}[A][B]=k[A][B]}$

where:

• ${\displaystyle Z_{AB}}$ is in unit mole collisions/s in 1 L of solution.
• ${\displaystyle N_{\text{A}}}$ is the Avogadro constant.
• ${\displaystyle R}$ is the radius of the collision cross-section, unit m.
• ${\displaystyle D_{r}}$ is the relative diffusion constant between A and B, unit m²/s.
• ${\displaystyle [A]}$ and ${\displaystyle [B]}$ are the molar concentrations of A and B respectively, unit mol/L.
• ${\displaystyle k}$ is the diffusive collision rate constant, unit L mol⁻¹ s⁻¹.

There have been a lot of extensions and modifications to the Smoluchowski model since it was proposed in 1916.

In 2022, Chen rationales that because the diffusive flux is evolving over time and the distance between the molecules has a finite value at a given concentration, there should be a critical time to cut off the evolution of the flux that will give a value much larger than the infinite solution Smoluchowski has proposed.^[5] So he proposes to use the average time for two molecules to switch places in the solution as the critical cut-off time, i.e., first neighbor visiting time. Although an alternative time could be the mean free path time or the average first passenger time, it overestimates the concentration gradient between the original location of the first passenger to the target. This hypothesis yields a fractal reaction kinetic rate equation of diffusive collision in a diluted solution:^[5]

${\displaystyle Z_{AB}=(1000N_{A})^{4/3}*8\pi ^{-1}A\beta D_{r}([A]+[B])^{1/3}[A][B]=k([A]+[B])^{1/3}[A][B]}$

where:

• ${\displaystyle Z_{AB}}$ is in unit mole collisions/s in 1 L of solution.
• ${\displaystyle N_{\text{A}}}$ is the Avogadro constant.
• ${\displaystyle A}$ is the area of the collision cross-section in unit m².
• ${\displaystyle \beta }$ is the product of the unitless fractions of reactive surface area on A and B. ${\displaystyle \beta A}$ is the effective adsorption cross-section area.
• ${\displaystyle D_{r}}$ is the relative diffusion constant between A and B, unit m²/s, and ${\displaystyle D_{r}=D_{A}+D_{B}}$.
• ${\displaystyle [A]}$ and ${\displaystyle [B]}$ are the molar concentrations of A and B respectively, unit mol/L.
• ${\displaystyle k}$ is the diffusive collision rate constant, unit L^4/3 mol^-4/3 s⁻¹.

See also

• Two-dimensional gas
• Rate equation

References

1. Trautz, Max. Das Gesetz der Reaktionsgeschwindigkeit und der Gleichgewichte in Gasen. Bestätigung der Additivität von C_v − 3/2 R. Neue Bestimmung der Integrationskonstanten und der Moleküldurchmesser, Zeitschrift für anorganische und allgemeine Chemie, Volume 96, Issue 1, Pages 1–28, (1916).
2. IUPAC, Compendium of Chemical Terminology, 2nd ed. (the "Gold Book") (1997). Online corrected version: (2006–) "collision theory". doi:10.1351/goldbook.C01170
3. William Cudmore McCullagh Lewis, XLI.—Studies in catalysis. Part IX. The calculation in absolute measure of velocity constants and equilibrium constants in gaseous systems, J. Chem. Soc., Trans., 1918, 113, 471-492.
4. Smoluchowski, Marian (1916). "Drei Vorträge über Diffusion, Brownsche Molekularbewegung und Koagulation von Kolloidteilchen". Phys. Z. (in German). 17: 557–571, 585–599. Bibcode:1916ZPhy...17..557S.
5. Chen, Jixin (2022). "Why Should the Reaction Order of a Bimolecular Reaction be 2.33 Instead of 2?". J. Phys. Chem. A. 126: 9719–9725. doi:10.1021/acs.jpca.2c07500. PMC 9805503.
6. "6.1.6: The Collision Theory". 2 October 2013.
7. IUPAC, Compendium of Chemical Terminology, 2nd ed. (the "Gold Book") (1997). Online corrected version: (2006–) "steric factor". doi:10.1351/goldbook.S05998
8. IUPAC, Compendium of Chemical Terminology, 2nd ed. (the "Gold Book") (1997). Online corrected version: (2006–) "collision frequency". doi:10.1351/goldbook.C01166
9. Kenneth Connors, Chemical Kinetics, 1990, VCH Publishers.
10. E.A. Moelwyn-Hughes, The kinetics of reactions in solution, 2nd ed, page 71.

External links

• Introduction to Collision Theory