Unit 8: Chain Reactions, Branched Kinetics & Thermal/Branching Explosions

Chain reactions are complex chemical networks in which a reactive intermediate (chain carrier, typically a free radical or atom) is consumed in a reaction step that regenerates one or more new chain carriers, enabling a single initiation event to trigger hundreds or thousands of product-forming cycles.

The Four Essential Stages

1. Initiation: Generation of active chain carriers from stable closed-shell molecules, driven by thermal activation, photolysis, or chemical initiators:

2. Propagation: Elementary reactions between chain carriers and stable molecules that yield final products while regenerating active chain carriers:

3. Inhibition / Retardation: Reversible or irreversible scavenging of chain carriers by products or added inhibitors:

4. Termination: Destruction of active chain carriers, halting the propagation cycle:

• Homogeneous (Gas-Phase): Binary or termolecular radical recombination:

• Heterogeneous (Wall): Diffusion to the vessel wall followed by adsorption and recombination:

Chain Length ($\Lambda_{\text{chain}}$)

The kinetic chain length $\Lambda_{\text{chain}}$ is defined as the average number of propagation cycles executed per initiation event:

For highly exothermic chain systems like $H_2 + Cl_2$, $\Lambda_{\text{chain}} \sim 10^5\text{--}10^6$, meaning a tiny flash of light can cause explosive transformation.

Hydrocarbon Combustion & Explosion Limits Reference Matrix

| Fuel / Oxidizer System | Stoichiometric Composition | Lower Explosion Limit (LEL, $\%$) | Upper Explosion Limit (UEL, $\%$) | Autoignition Temp ($T_{\text{auto}}$, $^\circ\text{C}$) | Laminar Flame Speed $S_L$ ($\text{cm/s}$) | Adiabatic Flame Temp ($T_b$, $\text{K}$) | |---|---|---|---|---|---|---| | Hydrogen / Air | $29.6\%\text{ }H_2$ | $4.0\%$ | $75.0\%$ | $560^\circ\text{C}$ | $210$ | $2380$ | | Methane / Air | $9.5\%\text{ }CH_4$ | $5.0\%$ | $15.0\%$ | $580^\circ\text{C}$ | $38$ | $2220$ | | Propane / Air | $4.0\%\text{ }C_3H_8$ | $2.1\%$ | $9.5\%$ | $470^\circ\text{C}$ | $43$ | $2260$ | | Ethylene / Air | $6.5\%\text{ }C_2H_4$ | $2.7\%$ | $36.0\%$ | $450^\circ\text{C}$ | $68$ | $2375$ | | Acetylene / Air | $7.7\%\text{ }C_2H_2$ | $2.5\%$ | $100.0\%$ (pure decomposition) | $305^\circ\text{C}$ | $155$ | $2540$ | | Carbon Monoxide / Air| $29.6\%\text{ }CO$ | $12.5\%$ | $74.0\%$ | $609^\circ\text{C}$ | $45$ (moist) | $2385$ |

Max Bodenstein (1906) experimentally discovered that the gas-phase reaction $H_2 + Br_2 \longrightarrow 2 HBr$ follows an extraordinarily intricate empirical rate law:

The Christiansen-Kramers-Polanyi Mechanism (1919)

1. Initiation:

2. Propagation 1:

3. Propagation 2:

4. Inhibition (Product Retardation):

5. Termination:

Derivation via Bodenstein SSA

Apply SSA to the two chain carriers, $[Br^\bullet]$ and $[H^\bullet]$:

Adding both steady-state equations:

From the first equation:

The rate of product formation is:

Substituting $[H^\bullet]$:

This derivation mathematically reproduces Bodenstein's empirical law, proving that $k = 2 k_2 (k_1/k_5)^{1/2}$ and $m = k_4 / k_3$.

Chain-growth free-radical polymerization represents an industrial manifestation of stationary chain kinetics.

Kinetic Steps

1. Initiation: Thermal decomposition of initiator $I$ (e.g., AIBN, benzoyl peroxide) followed by addition to monomer $M$:

where $f$ is initiator efficiency ($0.5 < f < 0.9$).

2. Propagation: Sequential addition of monomer units:

assuming rate constant $k_p$ is independent of radical chain length $n$.

3. Termination: Bimolecular combination or disproportionation:

Steady-State Radical Concentration

Applying the Bodenstein SSA to total active radicals $[M^\bullet]$:

Overall Polymerization Rate ($r_p$)

The rate is strictly first-order in monomer and half-order in initiator.

Number-Average Degree of Polymerization ($\bar{X}_n$)

Increasing initiator concentration increases polymerization rate but decreases average polymer molecular weight.

In a branched chain reaction, a propagation step generates more chain carriers than it consumes. Nikolai Semenov (1934; Nobel Prize 1956) formulated the mathematical theory of branching explosion kinetics.

Semenov Differential Equation

Let $n(t)$ be the concentration of active chain carriers.

1. Rate of initiation: $w_0$
2. Rate of chain branching: $f \, n$ (where $f$ is the branching probability coefficient per second)
3. Rate of chain termination: $g \, n$ (where $g$ is the termination probability coefficient per second)

The differential equation for radical population growth is:

where $\phi = f - g$ is the net branching factor.

Solution and Stability Criteria

Integrating with initial condition $n(0) = 0$:

Three distinct regimes emerge:

1. Stationary Sub-Critical Regime ($\phi < 0$, $g > f$):

Termination exceeds branching. As $t \to \infty$, $e^{\phi t} \to 0$:

The radical population reaches a finite, stable steady state. The reaction proceeds smoothly and slowly.

2. Critical Limit ($\phi = 0$, $f = g$):

Branching exactly balances termination:

Linear growth marking the boundary of explosion.

3. Explosive Super-Critical Regime ($\phi > 0$, $f > g$):

Branching exceeds termination. The exponential term $e^{\phi t}$ diverges:

The active radical concentration multiplies exponentially within milliseconds, producing an isothermal chain-branching explosion.

The reaction $2 H_2 + O_2 \longrightarrow 2 H_2O$ exhibits an iconic "explosion peninsula" in a pressure-temperature ($P-T$) phase diagram, governed by competition between branching and termination.

Elementary Branching Mechanism

• Initiation:

• Branching Step 1:

• Branching Step 2:

• Propagation:

Net result of the branching cycle:

From one $H$ atom, three active radicals are produced (net generation of $+2$ radicals).

The Three Explosion Limits

1. First (Lower) Limit ($P_1$):

At very low pressures (few Torr), mean free path is long. Radicals diffuse rapidly to the reactor walls where they are destroyed ($g_{\text{wall}} \propto D / d^2 \propto 1 / (P d^2)$). Explosion occurs when branching overcomes wall loss:

2. Second (Upper) Limit ($P_2$):

As pressure rises, termolecular gas-phase termination becomes dominant:

The hydroperoxyl radical $HO_2^\bullet$ is unreactive at moderate temperatures and diffuses to the wall without branching. Condition for explosion:

Because $k_2$ has high activation energy ($E_a \approx 70\text{ kJ/mol}$) and $k_4$ has near-zero activation energy, $P_2$ increases exponentially with temperature.

3. Third Limit ($P_3$):

At high pressures ($P > 1\text{ bar}$), $HO_2^\bullet$ begins to react via $HO_2^\bullet + H_2 \longrightarrow H_2O_2 + H^\bullet$, regenerating radicals while massive exothermic heat release triggers thermal runaway.

University Honors Research Monograph: Cool Flame Oscillations & Low-Temperature Hydrocarbon Autoignition

In advanced Homogeneous Charge Compression Ignition (HCCI) engines, hydrocarbon combustion exhibits two-stage autoignition accompanied by faint blue chemiluminescent cool flames ($T = 550 - 750\text{ K}$):

• Radical Addition to Dioxygen:

Alkyl radicals formed by initial $H$-atom abstraction react reversibly with molecular oxygen:

• Intramolecular Radical Isomerization:

The alkylperoxy radical ($RO_2^\bullet$) undergoes intramolecular hydrogen abstraction via a cyclic six-membered transition state to form a hydroperoxyalkyl radical ($^\bullet QOOH$):

• Chain Branching via Second $O_2$ Addition:

At intermediate temperatures, $^\bullet QOOH$ adds a second oxygen molecule:

The fragile $O-O$ peroxide bond of the ketohydroperoxide decomposes into two additional radicals ($OH^\bullet + \text{alkoxy radical}$), yielding net degenerate chain branching that drives first-stage autoignition.

• Negative Temperature Coefficient (NTC) Regime:

As temperature rises above $750\text{ K}$, the equilibrium $R^\bullet + O_2 \rightleftharpoons RO_2^\bullet$ shifts back toward reactants, while $RO_2^\bullet$ decomposes to non-branching conjugate alkene $+ HO_2^\bullet$. Overall oxidation rate paradoxically decreases with increasing temperature, creating the NTC phenomenon.

The oxidation of hydrocarbons ($RH + O_2$) exhibits complex kinetic behavior including two-stage ignition, cool flames, and negative temperature coefficient (NTC) behavior.

Degenerate Chain Branching (Semenov)

Unlike $H_2 + O_2$ where branching is instantaneous via unstable atoms, hydrocarbon oxidation forms relatively stable molecular intermediates (hydroperoxides $ROOH$, aldehydes $RCHO$) that slowly decompose to produce radicals:

The hydroperoxide undergoes occasional unimolecular homolysis:

Because $ROOH$ has a lifetime of seconds, radical multiplication is delayed: degenerate branching.

Cool Flame Phenomena and NTC Regime

Between $300^\circ\text{C}$ and $400^\circ\text{C}$, hydrocarbons exhibit a pale bluish luminescence called a cool flame, emitted by electronically excited formaldehyde ($HCHO^* \to HCHO + h\nu$). In this temperature window, the reaction rate slows down as temperature increases (Negative Temperature Coefficient, NTC).

• Kinetic Origin: The peroxy radical equilibrium $R^\bullet + O_2 \rightleftharpoons RO_2^\bullet$ shifts backward at higher temperatures, favoring non-branching alkene formation ($R^\bullet + O_2 \to \text{alkene} + HO_2^\bullet$). This kinetic competition underpins engine knock in internal combustion engines and defines fuel octane ratings.

In contrast to isothermal branched-chain explosions, a thermal explosion occurs when the rate of exothermic chemical heat generation exceeds the rate of heat dissipation to the surroundings.

Semenov Thermal Explosion Model (Uniform Temperature)

Consider a reaction vessel of volume $V$, surface area $S$, and heat transfer coefficient $\chi$, containing an exothermic reaction ($\Delta H_r < 0$) with rate $r = k_0 e^{-E_a / R T} c^n$:

1. Rate of Heat Generation ($q_{\text{gen}}$):

Increases exponentially with temperature $T$.

2. Rate of Heat Removal ($q_{\text{loss}}$):

Increases linearly with temperature $T$ above ambient wall temperature $T_0$ (Newton's law of cooling).

Critical Semenov Condition

Steady-state heat balance requires $q_{\text{gen}} = q_{\text{loss}}$. The boundary of stability occurs when the heat generation curve is tangent to the heat loss line:

Evaluating derivatives:

Substituting $q_{\text{gen}} = S \chi (T_{\text{crit}} - T_0)$:

For typical activation energies ($E_a \approx 100\text{ kJ/mol}$, $T_0 \approx 300\text{ K}$):

If the self-heating exceeds this critical temperature rise $\Delta T_{\text{crit}}$, heat generation outstrips heat removal, triggering catastrophic thermal runaway.

Elucidating the intricate multi-step branching radical networks of hydrocarbon combustion and atmospheric ozone depletion requires isomeric identification of transient radicals and reactive peroxide intermediates.

1. The Challenge of Isomer Differentiation in Complex Radical Chains

In high-temperature oxidation of hydrocarbons ($RH + O_2$), dozens of isomeric radical intermediates with identical molecular weights coexist (e.g., resonance-stabilized propenyl vs. allyl radicals, $C_3H_5$, $m/z = 41$; cyclic QOOH vs. peroxy radicals). Conventional electron impact mass spectrometry ($70\text{ eV}$) causes extensive fragmentation, obscuring parent ion identities and destroying fragile organic peroxides ($ROOH$, $OOQOOH$).

2. Synchrotron Vacuum Ultraviolet Photoionization Mass Spectrometry (SVUV-PIMS)

Synchrotron Vacuum Ultraviolet Photoionization Mass Spectrometry couples laminar flow reactors, low-pressure flat-flame burners, or shock tubes to a tunable synchrotron radiation beamline ($7 - 15\text{ eV}$) with an orthogonal time-of-flight mass spectrometer.

Principle of Threshold Photoionization

Synchrotron VUV radiation provides continuously tunable photon energies with narrow bandwidth ($\Delta E \approx 0.01\text{ eV}$). When photon energy $h\nu$ exceeds the adiabatic ionization potential of a molecule ($h\nu \ge AIE$), single-photon "soft" ionization occurs:

Because the excess energy ($h\nu - AIE$) is negligible, zero dissociative fragmentation occurs, preserving fragile radical cations intact as parent molecular ions.

Photoionization Efficiency (PIE) Curves

By scanning photon energy across an energy interval and recording ion intensity at a fixed mass-to-charge ratio $m/z$, a Photoionization Efficiency (PIE) curve is measured:

Each chemical isomer has a distinct ionization threshold (AIE) and Franck-Condon structural overlap spectrum:

• Allyl radical ($\text{H}_2\text{C=CH-CH}_2^\bullet$): $AIE = 8.13\text{ eV}$.
• 2-Propenyl radical ($\text{H}_2\text{C=C}^\bullet\text{-CH}_3$): $AIE = 8.68\text{ eV}$.
• 1-Propenyl radical ($\text{CH}_3\text{-CH=CH}^\bullet$): $AIE = 7.70\text{ eV}$.

Fitting experimental PIE curves to absolute calibrated photoionization cross-sections deconvolutes isomeric mole fraction profiles $[X_i](T, z)$ through flame fronts with sub-millimeter spatial resolution.

3. Direct Discovery of the Elusive Hydroperoxyalkyl ($QOOH$) Intermediates

The low-temperature oxidation regime ($500 - 800\text{ K}$) that governs automotive knock and autoignition relies on the central chain-branching pathway:

For decades, $^\bullet QOOH$ radicals eluded detection because the isomerization barrier is high and subsequent reaction is rapid. In 2012, Taatjes, Welz, Osborn, and co-workers used SVUV-PIMS at the Advanced Light Source to directly detect the simplest $QOOH$ radical ($\cdot CH_2CH_2OOH$) in photolysis flow reactors, confirming the mechanistic core of combustion chemical modeling.