Part II: Reaction Mechanisms | Chapter 1

Nucleophilic Substitution Reactions (S_N1 & S_N2)

A rigorous treatment of bimolecular and unimolecular nucleophilic substitution at sp³ carbon, including mechanistic derivations, stereochemical outcomes, solvent effects, nucleophilicity parameters, and leaving group correlations.

1. Introduction: Substitution at sp³ Carbon

Nucleophilic substitution is one of the most fundamental reaction types in organic chemistry. In these reactions, a nucleophile (an electron-rich species) replaces a leaving group bonded to an sp³-hybridized carbon. The general transformation is:

$\text{Nu}^{-} + \text{R-LG} \longrightarrow \text{Nu-R} + \text{LG}^{-}$

Two limiting mechanisms describe this transformation. The S_N2 mechanism is a concerted, one-step process where bond formation and bond breaking occur simultaneously. The S_N1 mechanism is a two-step process involving prior ionization to form a carbocation intermediate, followed by nucleophilic capture.

The "S" stands for substitution, "N" for nucleophilic, and the numeral (1 or 2) denotes the molecularity of the rate-determining step — that is, how many species appear in the transition state of the slow step.

Understanding which mechanism operates for a given substrate requires analyzing steric effects, electronic effects, solvent polarity, nucleophile strength, and leaving group ability. This chapter derives the kinetic rate laws, stereochemical predictions, and thermodynamic factors that govern each pathway.

2. Derivation 1: The S_N2 Mechanism

2.1 Bimolecular Rate Law

The S_N2 reaction is bimolecular: both the nucleophile and the substrate appear in the rate-determining (and only) transition state. Experimentally, the rate depends on the concentration of both species:

$\text{rate} = k_2 [\text{Nu}^{-}][\text{R-X}]$

This second-order rate law (first order in each reactant) is the defining kinetic signature of S_N2. Doubling either the nucleophile or substrate concentration doubles the rate. Doubling both quadruples it.

2.2 Transition State Theory Derivation

From Eyring's transition state theory (TST), the rate constant for any elementary reaction is:

$k = \frac{k_B T}{h} \exp\!\left(-\frac{\Delta G^{\ddagger}}{RT}\right)$

where $k_B$ is Boltzmann's constant, $T$ is temperature, $h$ is Planck's constant, and $\Delta G^{\ddagger}$ is the Gibbs free energy of activation. For the S_N2 reaction, the activated complex is a pentacoordinate species:

$\text{Nu}^{-} + \text{R-X} \rightleftharpoons{K^{\ddagger}} [\text{Nu} \cdots \text{R} \cdots \text{X}]^{\ddagger} \longrightarrow \text{Nu-R} + \text{X}^{-}$

The pseudo-equilibrium constant for forming the transition state is:

$K^{\ddagger} = \frac{[\text{TS}]}{[\text{Nu}^{-}][\text{R-X}]} = \exp\!\left(-\frac{\Delta G^{\ddagger}}{RT}\right)$

The rate of product formation equals the frequency of barrier crossing times the transition state concentration:

$\text{rate} = \frac{k_B T}{h} [\text{TS}] = \frac{k_B T}{h} K^{\ddagger} [\text{Nu}^{-}][\text{R-X}]$

Therefore the second-order rate constant is:

$k_2 = \frac{k_B T}{h} \exp\!\left(-\frac{\Delta G^{\ddagger}}{RT}\right)$

Expanding $\Delta G^{\ddagger} = \Delta H^{\ddagger} - T \Delta S^{\ddagger}$ gives the Eyring equation. For S_N2, the activation entropy $\Delta S^{\ddagger}$ is typically large and negative (around $-80$ to $-120 \;\text{J mol}^{-1}\text{K}^{-1}$) because two particles combine into one transition state, losing translational freedom.

2.3 Backside Attack and Walden Inversion

The S_N2 mechanism proceeds through backside attack: the nucleophile approaches the electrophilic carbon from the side opposite the leaving group. This geometry is dictated by orbital symmetry — the nucleophile's HOMO (highest occupied molecular orbital) must overlap with the $\sigma^{*}$ antibonding orbital of the C–X bond.

Orbital Picture

The HOMO of the nucleophile (e.g., a lone pair on $\text{HO}^{-}$) donates electron density into the $\sigma^{*}_{\text{C-X}}$ orbital, which has its larger lobe on the side opposite the leaving group. This HOMO → $\sigma^{*}$ interaction simultaneously forms the new Nu–C bond and breaks the C–X bond:

$\text{HOMO}(\text{Nu}^{-}) \;\longrightarrow\; \sigma^{*}(\text{C-X})$

The stereochemical consequence is Walden inversion (also called "inversion of configuration"). If the substrate carbon is a stereocenter, the product has the opposite configuration:

$(R)\text{-substrate} \xrightarrow{\text{S}_N 2} (S)\text{-product}$

This is analogous to an umbrella inverting in the wind — the three substituents flip to the opposite face as the nucleophile pushes through.

2.4 Steric Effects on S_N2 Rates

Because the S_N2 transition state is pentacoordinate, steric crowding at the reaction center dramatically raises $\Delta G^{\ddagger}$ and slows the reaction. The relative rates for alkyl bromides reacting with a common nucleophile are:

Substrate	Structure	Relative Rate
Methyl	$\text{CH}_3\text{Br}$	30
Primary (1°)	$\text{CH}_3\text{CH}_2\text{Br}$	1.0
Secondary (2°)	$(\text{CH}_3)_2\text{CHBr}$	0.025
Tertiary (3°)	$(\text{CH}_3)_3\text{CBr}$	≈ 0

The trend is clear: methyl > 1° > 2° >> 3°. Tertiary substrates effectively do not undergo S_N2 because the three bulky substituents block backside approach of the nucleophile. Each additional alkyl group on the reaction center introduces 1,3-diaxial-type steric interactions with the incoming nucleophile.

2.5 Energy Diagram: S_N2

The S_N2 energy profile has a single transition state and no intermediate. The reaction coordinate diagram shows reactants climbing directly to the pentacoordinate transition state, then descending to products. The height of the barrier ($\Delta G^{\ddagger}$) determines the rate. For an exergonic S_N2 reaction, the products lie lower in energy than the reactants ($\Delta G^{\circ} < 0$).

3. Derivation 2: The S_N1 Mechanism

3.1 Unimolecular Rate Law

The S_N1 reaction is unimolecular: only the substrate appears in the rate-determining step. The experimental rate law is first order:

$\text{rate} = k_1 [\text{R-X}]$

The rate does not depend on nucleophile concentration. This is the defining kinetic signature of S_N1 and tells us the nucleophile is not involved in the rate-determining step.

3.2 Two-Step Mechanism

The S_N1 mechanism proceeds in two distinct steps:

Step 1 (slow, rate-determining): Ionization

$\text{R-X} \xrightarrow{k_1} \text{R}^{+} + \text{X}^{-}$

Step 2 (fast): Nucleophilic capture

$\text{R}^{+} + \text{Nu}^{-} \xrightarrow{k_2} \text{Nu-R}$

3.3 Steady-State Derivation of the Rate Law

We can rigorously derive the S_N1 rate law using the steady-state approximation on the carbocation intermediate $\text{R}^{+}$. The full mechanism includes the reverse of step 1:

$\text{R-X} \underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}} \text{R}^{+} + \text{X}^{-}$

$\text{R}^{+} + \text{Nu}^{-} \xrightarrow{k_2} \text{Nu-R}$

Applying the steady-state approximation to $[\text{R}^{+}]$:

$\frac{d[\text{R}^{+}]}{dt} = k_1[\text{R-X}] - k_{-1}[\text{R}^{+}][\text{X}^{-}] - k_2[\text{R}^{+}][\text{Nu}^{-}] = 0$

Solving for $[\text{R}^{+}]$:

$[\text{R}^{+}] = \frac{k_1[\text{R-X}]}{k_{-1}[\text{X}^{-}] + k_2[\text{Nu}^{-}]}$

The rate of product formation is:

$\text{rate} = k_2[\text{R}^{+}][\text{Nu}^{-}] = \frac{k_1 k_2 [\text{R-X}][\text{Nu}^{-}]}{k_{-1}[\text{X}^{-}] + k_2[\text{Nu}^{-}]}$

Now we apply the S_N1 limiting condition: the nucleophilic capture step is much faster than recombination with $\text{X}^{-}$, i.e., $k_2[\text{Nu}^{-}] \gg k_{-1}[\text{X}^{-}]$. The denominator simplifies:

$\text{rate} \approx \frac{k_1 k_2 [\text{R-X}][\text{Nu}^{-}]}{k_2[\text{Nu}^{-}]} = k_1[\text{R-X}]$

The nucleophile concentration cancels, yielding the characteristic first-order rate law. This is the mathematical proof that when ionization is rate-determining, the observed rate depends only on substrate concentration.

3.4 Carbocation Stability and Substrate Effects

Since the rate-determining step of S_N1 is carbocation formation, anything that stabilizes the carbocation lowers $\Delta G^{\ddagger}$ and accelerates the reaction. By Hammond's postulate, the transition state for an endergonic ionization step resembles the carbocation intermediate. The stability order is:

$\text{3}^{\circ} > \text{2}^{\circ} > \text{1}^{\circ} > \text{methyl}$

This ordering arises from hyperconjugation: adjacent C–H $\sigma$ bonding orbitals donate electron density into the empty p orbital on the carbocation center. More alkyl groups provide more hyperconjugative stabilization:

$\sigma_{\text{C-H}} \;\longrightarrow\; p^{+}_{\text{empty}}$

Additionally, alkyl groups are weakly electron-donating by induction, further stabilizing the positive charge. The combined effect means tertiary substrates are approximately $10^3$–$10^6$ times more reactive in S_N1 than primary substrates.

3.5 Stereochemistry: Racemization

The carbocation intermediate is sp²-hybridized and planar. The nucleophile can attack from either face of the empty p orbital with roughly equal probability:

$\text{R}^{+}\text{ (planar)} + \text{Nu}^{-} \longrightarrow \begin{cases} (R)\text{-product} & (50\%) \\ (S)\text{-product} & (50\%) \end{cases}$

In practice, S_N1 reactions at a stereocenter give predominantly racemization, though slight excess inversion is sometimes observed due to ion-pair effects — the departing leaving group partially shields one face of the carbocation (front-side shielding).

3.6 Energy Diagram: S_N1

The S_N1 energy profile has two transition states and a carbocation intermediate between them. The first transition state (ionization, the rate-determining step) has the highest energy. The intermediate sits in a local energy minimum. The second transition state (nucleophilic capture) is lower in energy. The overall barrier $\Delta G^{\ddagger}_1$ determines the observed rate.

4. Derivation 3: Solvent Effects (Hughes–Ingold Rules)

4.1 The Hughes–Ingold Framework

Edward D. Hughes and Christopher K. Ingold developed a systematic framework in the 1930s for predicting how solvent polarity affects substitution and elimination rates. The key principle is:

"An increase in solvent polarity accelerates reactions that develop charge in the transition state and decelerates reactions that disperse charge in the transition state."

4.2 Derivation: Solvation Free Energy Analysis

The observed rate constant depends on the difference in solvation free energy between the transition state and the reactants:

$\Delta G^{\ddagger}_{\text{solution}} = \Delta G^{\ddagger}_{\text{gas}} + \Delta G_{\text{solv}}(\text{TS}) - \Delta G_{\text{solv}}(\text{reactants})$

Case 1: S_N1 (charge generation)

In S_N1, the rate-determining step converts a neutral substrate into a transition state with developing positive and negative charges:

$\underbrace{\text{R-X}}_{\text{neutral}} \longrightarrow \underbrace{[\delta^{+}\text{R} \cdots \text{X}^{\delta^{-}}]^{\ddagger}}_{\text{charge developing}} \longrightarrow \text{R}^{+} + \text{X}^{-}$

The transition state is more polar than the reactants. Therefore:

$|\Delta G_{\text{solv}}(\text{TS})| > |\Delta G_{\text{solv}}(\text{reactants})|$

In a more polar solvent, the transition state is preferentially stabilized relative to the reactants, lowering $\Delta G^{\ddagger}_{\text{solution}}$ and increasing the rate. Polar protic solvents (water, methanol, formic acid) are particularly effective because they stabilize both the cation and the anion through hydrogen bonding and electrostatic interactions.

Case 2: S_N2 (charge dispersal)

In the typical S_N2 case (anionic nucleophile + neutral substrate), the charge is localized on the nucleophile in the reactants but dispersed across two atoms in the transition state:

$\underbrace{\text{Nu}^{-}}_{\text{localized charge}} + \text{R-X} \longrightarrow \underbrace{[\text{Nu}^{\delta^{-}} \cdots \text{R} \cdots \text{X}^{\delta^{-}}]^{\ddagger}}_{\text{dispersed charge}} \longrightarrow \text{Nu-R} + \text{X}^{-}$

Here the reactant anion has a larger solvation energy than the diffuse transition state. In polar protic solvents, the nucleophile is tightly solvated by hydrogen bonds, which must be partially stripped away to reach the transition state. This raises the effective barrier:

$\Delta G^{\ddagger}_{\text{protic}} = \Delta G^{\ddagger}_{\text{gas}} + \underbrace{\Delta G_{\text{solv}}(\text{TS})}_{\text{smaller}} - \underbrace{\Delta G_{\text{solv}}(\text{Nu}^{-})}_{\text{larger (more negative)}}$

Switching to a polar aprotic solvent (DMSO, DMF, acetonitrile, acetone) changes the picture. These solvents dissolve salts but cannot hydrogen-bond to anions. The nucleophile is poorly solvated — essentially "naked" — making it far more reactive:

$k_{\text{SN2}}(\text{DMSO}) \gg k_{\text{SN2}}(\text{MeOH})$

Rate enhancements of $10^3$ to $10^6$ are typical when switching from protic to aprotic solvents for S_N2 reactions.

4.3 Summary of Hughes–Ingold Rules

Mechanism	Charge Change	Polar Protic	Polar Aprotic
S_N1	Generation	Faster	Moderate
S_N2 (Nu¹⁻)	Dispersal	Slower	Faster

5. Derivation 4: Nucleophilicity vs Basicity

5.1 Kinetic vs Thermodynamic Measures

Nucleophilicity is a kinetic property: it measures the rate at which a Lewis base attacks an electrophilic carbon. It is quantified by how fast the species reacts in an S_N2 displacement.

Basicity is a thermodynamic property: it measures the equilibrium affinity for a proton, quantified by $\text{p}K_a$ of the conjugate acid.

$\text{Nucleophilicity:} \quad \text{Nu}^{-} + \text{CH}_3\text{-X} \xrightarrow{k_{\text{Nu}}} \text{products} \quad (\text{rate-based})$

$\text{Basicity:} \quad \text{B}^{-} + \text{H}^{+} \rightleftharpoons \text{BH} \quad (K_a\text{-based})$

These two properties often correlate but can diverge. For example, $\text{I}^{-}$ is a weak base ($\text{p}K_a(\text{HI}) \approx -10$) but an excellent nucleophile because of its large, polarizable electron cloud.

5.2 The Swain–Scott Equation

In 1953, C. Gardner Swain and Carleton B. Scott proposed a linear free-energy relationship to quantify nucleophilicity. They defined a nucleophilicity parameter $n$ for each nucleophile by measuring its rate relative to water attacking methyl bromide:

$\log\!\left(\frac{k_{\text{Nu}}}{k_{\text{H}_2\text{O}}}\right) = s \cdot n$

where:

●$k_{\text{Nu}}$ = rate constant for the nucleophile of interest
●$k_{\text{H}_2\text{O}}$ = rate constant for water (the reference nucleophile, $n = 0$)
●$n$ = nucleophilicity parameter (characteristic of the nucleophile)
●$s$ = substrate sensitivity parameter ($s = 1.00$ for CH₃Br by definition)

5.3 Derivation from Free Energy

The Swain–Scott equation is a linear free-energy relationship (LFER). Starting from transition state theory:

$\ln k_{\text{Nu}} = \ln\!\left(\frac{k_B T}{h}\right) - \frac{\Delta G^{\ddagger}_{\text{Nu}}}{RT}$

$\ln k_{\text{H}_2\text{O}} = \ln\!\left(\frac{k_B T}{h}\right) - \frac{\Delta G^{\ddagger}_{\text{H}_2\text{O}}}{RT}$

Subtracting:

$\ln\!\left(\frac{k_{\text{Nu}}}{k_{\text{H}_2\text{O}}}\right) = -\frac{\Delta G^{\ddagger}_{\text{Nu}} - \Delta G^{\ddagger}_{\text{H}_2\text{O}}}{RT} = -\frac{\delta\Delta G^{\ddagger}}{RT}$

Converting to log base 10 and identifying $s \cdot n = -\delta\Delta G^{\ddagger} / (2.303\,RT)$:

$\log\!\left(\frac{k_{\text{Nu}}}{k_{\text{H}_2\text{O}}}\right) = s \cdot n$

The parameter $n$ thus encodes the differential stabilization of the transition state by the nucleophile relative to water. High $n$ values (e.g., $n \approx 5$ for I⁻, CN⁻, RS⁻) indicate powerful nucleophiles. The substrate sensitivity $s$ measures how responsive a given electrophilic center is to nucleophile quality.

5.4 Factors Affecting Nucleophilicity

Charge

Anions are stronger nucleophiles than their neutral counterparts:$\text{HO}^{-} \gg \text{H}_2\text{O}$ and $\text{NH}_2^{-} \gg \text{NH}_3$.

Polarizability

Larger atoms with more diffuse electron clouds form better overlap with $\sigma^{*}_{\text{C-X}}$. Down a column: $\text{I}^{-} > \text{Br}^{-} > \text{Cl}^{-} > \text{F}^{-}$ in protic solvents.

Electronegativity

Less electronegative atoms hold their electrons more loosely, making them better nucleophiles. Across a row: C > N > O > F.

Steric Effects

Bulky nucleophiles react more slowly even if strongly basic.$t$-BuO⁻ is a strong base but poor nucleophile due to steric bulk.

6. Derivation 5: Leaving Group Ability

6.1 Correlation with Conjugate Acid pK_a

A good leaving group must be able to depart with the bonding electrons and stabilize itself as a free anion or neutral molecule. The key thermodynamic predictor is the pK_a of the conjugate acid of the leaving group:

Good leaving groups are the conjugate bases of strong acids (low pK_a).

The reasoning is straightforward: if HX is a strong acid (low pK_a), then X⁻ is a stable, weak base that is thermodynamically comfortable bearing a negative charge. This same stability makes X⁻ a good leaving group because the C–X bond heterolysis produces a stable anion.

6.2 Derivation: Why I⁻ > Br⁻ > Cl⁻ > F⁻

For the halide leaving groups, we can derive the ordering from two physical factors: bond dissociation energies and solvation energies.

Factor 1: C–X Bond Strength

The enthalpy required to break the C–X bond heterolytically (in the gas phase) decreases as we go down the group:

Bond	BDE (kJ/mol)	pK_a(HX)
C–F	485	3.2
C–Cl	339	−7
C–Br	285	−9
C–I	222	−10

Weaker bonds (lower BDE) are easier to break. The C–I bond is the weakest because iodine's large atomic radius gives poor overlap with carbon's 2p orbital, producing a long, weak bond.

Factor 2: Solvation of X⁻

In solution, the leaving group anion must be stabilized by solvation. The total energy for leaving group departure is:

$\Delta G_{\text{departure}} = \text{BDE}(\text{C-X}) - \Delta G_{\text{solv}}(\text{X}^{-}) - \text{EA}(\text{X})$

where EA is the electron affinity of X. Although F⁻ has the largest solvation energy in protic solvents (due to its high charge density), this does not compensate for the very strong C–F bond. The net effect across the series is:

$\Delta G_{\text{departure}}(\text{I}^{-}) < \Delta G_{\text{departure}}(\text{Br}^{-}) < \Delta G_{\text{departure}}(\text{Cl}^{-}) < \Delta G_{\text{departure}}(\text{F}^{-})$

Therefore: I⁻ > Br⁻ > Cl⁻ > F⁻ as leaving groups. The bond strength factor dominates the trend. Fluoride is an exceptionally poor leaving group because the C–F bond is the strongest single bond to carbon.

6.3 Common Leaving Groups Ranked

Leaving Group	Conjugate Acid	pK_a	Quality
N₂ (diazonium)	—	—	Superb
OTs, OMs	TsOH, MsOH	−1 to −2	Excellent
I⁻	HI	−10	Excellent
Br⁻	HBr	−9	Good
Cl⁻	HCl	−7	Moderate
H₂O	H₃O⁺	−1.7	Moderate
F⁻	HF	3.2	Poor
HO⁻	H₂O	15.7	Very poor

7. Applications

7.1 Williamson Ether Synthesis

The Williamson ether synthesis is a classic S_N2 reaction used to prepare ethers from an alkoxide nucleophile and a primary alkyl halide:

$\text{R-O}^{-} + \text{R}'\text{-X} \xrightarrow{\text{S}_N 2} \text{R-O-R}' + \text{X}^{-}$

The alkoxide (strong nucleophile, strong base) attacks the primary or methyl halide via backside attack. The reaction must use primary or methyl halides to avoid elimination (E2) as a competing pathway with bulky alkoxides. For example, diethyl ether is synthesized from sodium ethoxide and ethyl bromide:

$\text{CH}_3\text{CH}_2\text{O}^{-}\text{Na}^{+} + \text{CH}_3\text{CH}_2\text{Br} \longrightarrow \text{CH}_3\text{CH}_2\text{-O-CH}_2\text{CH}_3 + \text{NaBr}$

7.2 Alkyl Halide Transformations

Nucleophilic substitution is the principal method for converting alkyl halides into other functional groups. By choosing different nucleophiles, a single alkyl halide can be transformed into alcohols, ethers, amines, nitriles, thiols, and more:

$\text{R-Br} + \text{NaOH} \longrightarrow \text{R-OH} + \text{NaBr} \quad \text{(alcohol)}$

$\text{R-Br} + \text{NaCN} \longrightarrow \text{R-CN} + \text{NaBr} \quad \text{(nitrile, extends C chain)}$

$\text{R-Br} + \text{NaN}_3 \longrightarrow \text{R-N}_3 + \text{NaBr} \quad \text{(azide)}$

$\text{R-Br} + \text{NaSH} \longrightarrow \text{R-SH} + \text{NaBr} \quad \text{(thiol)}$

The cyanide reaction is particularly valuable because it creates a new C–C bond, extending the carbon chain by one atom. The azide reaction is used extensively in "click chemistry" for bioconjugation.

7.3 Drug Metabolism: Glutathione Conjugation

In biochemistry, nucleophilic substitution plays a critical role in Phase II drug metabolism. The tripeptide glutathione (GSH, $\gamma$-Glu-Cys-Gly) acts as a biological nucleophile through its cysteine thiol group:

$\text{GS}^{-} + \text{R-X} \xrightarrow{\text{GST}} \text{GS-R} + \text{X}^{-}$

The enzyme glutathione S-transferase (GST) catalyzes this S_N2 reaction, detoxifying electrophilic metabolites (including drug intermediates and epoxides) by conjugating them with glutathione. The thiolate $\text{GS}^{-}$ is an excellent nucleophile due to sulfur's high polarizability ($n \approx 5.1$ on the Swain–Scott scale).

This detoxification pathway is essential for clearing reactive electrophiles that would otherwise damage DNA and proteins. Genetic polymorphisms in GST enzymes contribute to individual variation in drug sensitivity and cancer susceptibility.

8. Historical Context

Ingold & Hughes (1930s)

Christopher K. Ingold and Edward D. Hughes at University College London established the mechanistic framework for nucleophilic substitution and elimination reactions. They introduced the S_N1, S_N2, E1, and E2 nomenclature that remains standard today, and developed the Hughes–Ingold solvent rules based on extensive kinetic studies of hundreds of reactions.

Saul Winstein (1950s–60s)

Winstein at UCLA discovered neighboring group participation and the concept of ion pairs as intermediates in S_N1 reactions. He showed that the simple free carbocation model was incomplete: contact ion pairs, solvent-separated ion pairs, and free ions represent a continuum that explains the partial racemization and partial inversion often observed in S_N1.

Swain & Scott (1953)

C. Gardner Swain and Carleton B. Scott at MIT introduced the quantitative nucleophilicity scale that bears their name. Their linear free-energy relationship provided the first systematic way to predict S_N2 rates from nucleophile and substrate parameters, paralleling the Hammett equation for aromatic reactivity.

Paul von Ragué Schleyer (1960s–2000s)

Schleyer used computational chemistry to challenge and refine the S_N1/S_N2 dichotomy. His work showed that many reactions have "borderline" mechanisms with characteristics of both pathways, and helped establish the modern view that S_N1 and S_N2 represent extremes on a mechanistic continuum.

9. Python Simulation: Energy Diagrams & Rate Comparisons

The simulation below generates three sets of plots: (1) S_N1 vs S_N2 reaction coordinate energy diagrams showing the contrasting energy profiles, (2) bar charts comparing relative rates for different substrate types under each mechanism, and (3) nucleophilicity parameters (Swain–Scott) and leaving group ability correlated with conjugate acid pK_a.

Uses numpy only. No scipy.

Simulation

Python

script.py231 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# ============================================================
# SN1 vs SN2 Reaction Coordinate Energy Diagrams
# ============================================================

fig, axes = plt.subplots(1, 2, figsize=(14, 6))
fig.suptitle('Reaction Coordinate Energy Diagrams: SN1 vs SN2', fontsize=14, fontweight='bold', color='white')
fig.patch.set_facecolor('#0f0a1a')

# --- SN2 Energy Diagram (single transition state) ---
ax1 = axes[0]
ax1.set_facecolor('#1a0f2e')

# SN2: reactants -> single TS -> products
x_sn2 = np.linspace(0, 10, 500)
# Smooth curve: reactants at E=2, TS at E=7, products at E=1
E_sn2 = 2.0 - 1.0 * (x_sn2 / 10.0) + 5.0 * np.exp(-((x_sn2 - 5.0)**2) / 1.5)

ax1.plot(x_sn2, E_sn2, 'c-', linewidth=2.5, label='SN2 pathway')
ax1.axhline(y=2.0, color='gray', linestyle='--', alpha=0.4, linewidth=0.8)
ax1.axhline(y=1.0, color='gray', linestyle='--', alpha=0.4, linewidth=0.8)

# Annotations
ax1.annotate('Reactants\nNu$^-$ + R-X', xy=(0.5, 2.2), fontsize=9, color='#f0f0f0', ha='center')
ax1.annotate('Transition State\n[Nu---R---X]$^{\\ddagger}$', xy=(5.0, 7.4), fontsize=9, color='#ff6b6b', ha='center')
ax1.annotate('Products\nNu-R + X$^-$', xy=(9.2, 1.3), fontsize=9, color='#f0f0f0', ha='center')

# Activation energy arrow
ax1.annotate('', xy=(2.5, 7.0), xytext=(2.5, 2.0),
            arrowprops=dict(arrowstyle='<->', color='#ff9999', lw=1.5))
ax1.text(3.1, 4.5, '$\\Delta G^{\\ddagger}$', fontsize=11, color='#ff9999')

ax1.set_xlabel('Reaction Coordinate', color='#cccccc', fontsize=11)
ax1.set_ylabel('Free Energy (G)', color='#cccccc', fontsize=11)
ax1.set_title('SN2 Mechanism', color='cyan', fontsize=13, fontweight='bold')
ax1.tick_params(colors='#999999')
ax1.set_xlim(0, 10)
ax1.set_ylim(-0.5, 9)
for spine in ax1.spines.values():
    spine.set_color('#444444')

# --- SN1 Energy Diagram (two transition states, carbocation intermediate) ---
ax2 = axes[1]
ax2.set_facecolor('#1a0f2e')

x_sn1 = np.linspace(0, 10, 500)
# Two humps: TS1 at x=3, intermediate at x=5, TS2 at x=7
E_sn1 = (2.0
         + 6.0 * np.exp(-((x_sn1 - 3.0)**2) / 1.0)
         + 3.0 * np.exp(-((x_sn1 - 7.0)**2) / 1.0)
         - 1.5 * (x_sn1 / 10.0))

ax2.plot(x_sn1, E_sn1, color='#ff6bb5', linewidth=2.5, label='SN1 pathway')
ax2.axhline(y=2.0, color='gray', linestyle='--', alpha=0.4, linewidth=0.8)

# Annotations
ax2.annotate('Reactants\nR-X', xy=(0.5, 2.3), fontsize=9, color='#f0f0f0', ha='center')
ax2.annotate('TS1\n[R---X]$^{\\ddagger}$', xy=(3.0, 8.3), fontsize=9, color='#ff6b6b', ha='center')
ax2.annotate('Carbocation\nR$^+$', xy=(5.0, 3.0), fontsize=9, color='#ffdd57', ha='center')
ax2.annotate('TS2\n[R---Nu]$^{\\ddagger}$', xy=(7.0, 5.3), fontsize=9, color='#ff6b6b', ha='center')
ax2.annotate('Products\nNu-R', xy=(9.2, 1.0), fontsize=9, color='#f0f0f0', ha='center')

# Rate-determining step arrow
ax2.annotate('', xy=(1.2, 7.8), xytext=(1.2, 2.0),
            arrowprops=dict(arrowstyle='<->', color='#ff9999', lw=1.5))
ax2.text(1.6, 4.8, '$\\Delta G^{\\ddagger}_1$\n(RDS)', fontsize=10, color='#ff9999')

ax2.set_xlabel('Reaction Coordinate', color='#cccccc', fontsize=11)
ax2.set_ylabel('Free Energy (G)', color='#cccccc', fontsize=11)
ax2.set_title('SN1 Mechanism', color='#ff6bb5', fontsize=13, fontweight='bold')
ax2.tick_params(colors='#999999')
ax2.set_xlim(0, 10)
ax2.set_ylim(-0.5, 9.5)
for spine in ax2.spines.values():
    spine.set_color('#444444')

plt.tight_layout(rect=[0, 0, 1, 0.93])
plt.savefig('output.png', dpi=130, bbox_inches='tight', facecolor='#0f0a1a')
plt.close()

# ============================================================
# Relative Rates vs Substrate Structure
# ============================================================

fig2, axes2 = plt.subplots(1, 2, figsize=(14, 6))
fig2.suptitle('Relative Reaction Rates vs Substrate Structure', fontsize=14, fontweight='bold', color='white')
fig2.patch.set_facecolor('#0f0a1a')

substrates = ['Methyl', 'Primary\n(1\u00b0)', 'Secondary\n(2\u00b0)', 'Tertiary\n(3\u00b0)']
x_pos = np.arange(len(substrates))

# SN2 relative rates (decreasing with steric bulk)
sn2_rates = [30.0, 1.0, 0.025, 0.0]  # tertiary essentially zero
# SN1 relative rates (increasing with carbocation stability)
sn1_rates = [0.0, 0.001, 1.0, 1200.0]  # methyl essentially zero

ax_sn2 = axes2[0]
ax_sn2.set_facecolor('#1a0f2e')
bars1 = ax_sn2.bar(x_pos, sn2_rates, color=['#00e5ff', '#00bcd4', '#0097a7', '#455a64'],
                     edgecolor='#00e5ff', linewidth=0.8, width=0.6)
ax_sn2.set_xticks(x_pos)
ax_sn2.set_xticklabels(substrates, fontsize=10, color='#cccccc')
ax_sn2.set_ylabel('Relative Rate', color='#cccccc', fontsize=11)
ax_sn2.set_title('SN2 Relative Rates', color='cyan', fontsize=13, fontweight='bold')
ax_sn2.tick_params(colors='#999999')
for spine in ax_sn2.spines.values():
    spine.set_color('#444444')

# Add rate labels on bars
for bar, rate in zip(bars1, sn2_rates):
    label = f'{rate:.1f}' if rate >= 0.01 else ('~0' if rate == 0 else f'{rate:.3f}')
    ax_sn2.text(bar.get_x() + bar.get_width()/2, bar.get_height() + 0.5,
               label, ha='center', va='bottom', color='white', fontsize=10, fontweight='bold')

ax_sn1 = axes2[1]
ax_sn1.set_facecolor('#1a0f2e')
bars2 = ax_sn1.bar(x_pos, sn1_rates, color=['#455a64', '#e91e63', '#f06292', '#ff6bb5'],
                     edgecolor='#ff6bb5', linewidth=0.8, width=0.6)
ax_sn1.set_xticks(x_pos)
ax_sn1.set_xticklabels(substrates, fontsize=10, color='#cccccc')
ax_sn1.set_ylabel('Relative Rate', color='#cccccc', fontsize=11)
ax_sn1.set_title('SN1 Relative Rates', color='#ff6bb5', fontsize=13, fontweight='bold')
ax_sn1.tick_params(colors='#999999')
for spine in ax_sn1.spines.values():
    spine.set_color('#444444')

# Add rate labels on bars
for bar, rate in zip(bars2, sn1_rates):
    label = f'{rate:.0f}' if rate >= 1 else ('~0' if rate == 0 else f'{rate:.3f}')
    ax_sn1.text(bar.get_x() + bar.get_width()/2, bar.get_height() + 20,
               label, ha='center', va='bottom', color='white', fontsize=10, fontweight='bold')

plt.tight_layout(rect=[0, 0, 1, 0.93])
plt.savefig('output2.png', dpi=130, bbox_inches='tight', facecolor='#0f0a1a')
plt.close()

# ============================================================
# Nucleophilicity & Leaving Group Trends
# ============================================================

fig3, axes3 = plt.subplots(1, 2, figsize=(14, 5))
fig3.suptitle('Nucleophilicity & Leaving Group Ability Trends', fontsize=14, fontweight='bold', color='white')
fig3.patch.set_facecolor('#0f0a1a')

# Swain-Scott nucleophilicity parameters
nucleophiles = ['H$_2$O', 'F$^-$', 'Cl$^-$', 'Br$^-$', 'I$^-$', 'HO$^-$', 'CN$^-$', 'HS$^-$']
n_params = [0.0, 2.0, 3.04, 3.89, 5.04, 4.20, 5.1, 5.1]

ax3a = axes3[0]
ax3a.set_facecolor('#1a0f2e')
colors_nuc = plt.cm.plasma(np.linspace(0.2, 0.9, len(nucleophiles)))
ax3a.barh(range(len(nucleophiles)), n_params, color=colors_nuc, edgecolor='white', linewidth=0.5, height=0.6)
ax3a.set_yticks(range(len(nucleophiles)))
ax3a.set_yticklabels(nucleophiles, fontsize=10, color='#cccccc')
ax3a.set_xlabel('Swain-Scott n parameter', color='#cccccc', fontsize=11)
ax3a.set_title('Nucleophilicity (Swain-Scott)', color='#ff9ff3', fontsize=13, fontweight='bold')
ax3a.tick_params(colors='#999999')
for spine in ax3a.spines.values():
    spine.set_color('#444444')

# Leaving group ability vs conjugate acid pKa
leaving_groups = ['I$^-$', 'Br$^-$', 'Cl$^-$', 'F$^-$', 'HO$^-$']
pka_values = [-10.0, -9.0, -7.0, 3.2, 15.7]
lg_ability = [5.0, 4.0, 3.0, 1.0, 0.01]  # relative leaving group ability

ax3b = axes3[1]
ax3b.set_facecolor('#1a0f2e')
scatter = ax3b.scatter(pka_values, lg_ability, c=['#00e5ff', '#00bcd4', '#0097a7', '#ff6b6b', '#ff6b6b'],
                        s=120, edgecolors='white', linewidths=1, zorder=5)
for i, lg in enumerate(leaving_groups):
    ax3b.annotate(lg, (pka_values[i], lg_ability[i]), textcoords="offset points",
                 xytext=(10, 8), fontsize=10, color='#f0f0f0')

# Trend line
pka_line = np.linspace(-12, 18, 200)
lg_line = 5.0 * np.exp(-0.12 * (pka_line + 10))
ax3b.plot(pka_line, lg_line, '--', color='#ff6bb5', alpha=0.6, linewidth=1.5)

ax3b.set_xlabel('Conjugate Acid pK$_a$', color='#cccccc', fontsize=11)
ax3b.set_ylabel('Relative Leaving Group Ability', color='#cccccc', fontsize=11)
ax3b.set_title('Leaving Group vs pK$_a$', color='#ff6bb5', fontsize=13, fontweight='bold')
ax3b.tick_params(colors='#999999')
for spine in ax3b.spines.values():
    spine.set_color('#444444')

plt.tight_layout(rect=[0, 0, 1, 0.93])
plt.savefig('output3.png', dpi=130, bbox_inches='tight', facecolor='#0f0a1a')
plt.close()

print("=" * 70)
print("NUCLEOPHILIC SUBSTITUTION REACTIONS: COMPUTATIONAL ANALYSIS")
print("=" * 70)

print("\n--- SN2 Relative Rates by Substrate ---")
print(f"  {'Substrate':<15} {'Relative Rate':<15} {'Explanation'}")
print(f"  {'Methyl':<15} {'30.0':<15} {'No steric hindrance, fastest SN2'}")
print(f"  {'Primary':<15} {'1.0':<15} {'Reference substrate'}")
print(f"  {'Secondary':<15} {'0.025':<15} {'Significant steric hindrance'}")
print(f"  {'Tertiary':<15} {'~0':<15} {'Too hindered, SN2 does not occur'}")

print("\n--- SN1 Relative Rates by Substrate ---")
print(f"  {'Substrate':<15} {'Relative Rate':<15} {'Explanation'}")
print(f"  {'Methyl':<15} {'~0':<15} {'Cannot form stable carbocation'}")
print(f"  {'Primary':<15} {'0.001':<15} {'Very unstable primary carbocation'}")
print(f"  {'Secondary':<15} {'1.0':<15} {'Moderately stable carbocation'}")
print(f"  {'Tertiary':<15} {'1200':<15} {'Highly stable 3-degree carbocation'}")

print("\n--- Swain-Scott Nucleophilicity Parameters ---")
for nuc, n_val in zip(['H2O', 'F-', 'Cl-', 'Br-', 'I-', 'HO-', 'CN-', 'HS-'], n_params):
    print(f"  {nuc:<8}  n = {n_val:.2f}")

print("\n--- Leaving Group Ability (correlates with conjugate acid pKa) ---")
for lg, pka in zip(['I-', 'Br-', 'Cl-', 'F-', 'HO-'], pka_values):
    quality = 'Excellent' if pka < -8 else ('Good' if pka < 0 else ('Poor' if pka < 10 else 'Very poor'))
    print(f"  {lg:<8}  pKa(HX) = {pka:>6.1f}   Quality: {quality}")

print("\n--- Hughes-Ingold Solvent Rules ---")
print("  SN1 (charge generation): favored by polar protic solvents")
print("    Reactant (neutral) -> TS (charge developing) -> ions")
print("    Polar protic solvents stabilize developing charges via H-bonding")
print("  SN2 (charge dispersal): favored by polar aprotic solvents")
print("    Concentrated charge (Nu-) -> TS (charge dispersed) -> products")
print("    Polar aprotic solvents do not solvate anions strongly")
print("    Result: naked, more reactive nucleophile")

print("\n[3 plots saved: energy diagrams, rate comparisons, nucleophilicity/LG trends]")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← Part II Overview Elimination Reactions (E1/E2) →