Science Journaling Club Founded 2024

INTERACTIVE COMPANION · VOLUME 1, ISSUE 3 · SPRING 2025

The Kinetics Bench

A live model accompanying “The Reaction That Gives You the Wrong Product Because You Were Impatient”

← Read the full article

Two reactions share one starting material, \(A \rightleftharpoons B\) and \(A \rightleftharpoons C\), integrated exactly (this page solves the same linear equations as the article's code, in closed form, so nothing here is stepped through time one tick at a time). Drag the barriers and stabilities below and watch the kinetic product \(B\) rise, peak, and fade while the thermodynamic product \(C\) slowly overtakes it. The defaults reproduce the article's headline numbers: a crossover at 111.7 days, a kinetic-limit ratio of 0.01770, and a \(B\) peak of 98.0% at 5,809 s.

Model 1. Watch the crossover happen

k₁ = s⁻¹ k₂ = s⁻¹ K₂/K₁ = crossover:

Model 2. The regime map: temperature against time

Same four barrier and stability sliders (Model 1's settings carry over); this map fixes nothing but temperature and time. The curve is the crossover time itself, recomputed at every temperature from 250 to 400 K. Below the curve, B still leads. Above it, C has taken over. The open circle marks the current temperature slider and the crossover time it implies.

at this T, the thermodynamic product wins in:
Three numbers to check by hand, on the default settings:
1. \(k_1 = 10^{13}\exp(-90000/(8.314\times298.15)) \approx 1.709\times10^{-3}\,\text{s}^{-1}\) — matches the k₁ readout above.
2. Kinetic limit \(k_2/k_1 = 3.025\times10^{-5} / 1.709\times10^{-3} \approx 0.01770\) — drag both barrier sliders to their minimum step size apart and read the earliest part of the B/C traces.
3. Thermodynamic limit \(K_2/K_1 = \exp((\Delta G_B-\Delta G_C)/RT) = \exp(25000/2478.8) \approx 2.40\times10^4\) — matches the ratio the two curves settle into on the right edge of Model 1's plot.

What to try

Set Ea2 equal to Ea1 and watch the crossover collapse to almost nothing: with no kinetic head start, the more stable product wins almost immediately. Then reset Ea2 and instead push ΔGC up toward ΔGB: the crossover time grows without bound and the regime map's boundary curve runs off the top of the chart, then the whole map turns to the "never" colour once C stops being more stable at all. Finally leave the barriers and stabilities alone and just drag temperature: the boundary curve in Model 2 is this article's Figure 4, live.