4. Aromatic Chemistry

1. Introduction: Benzene and Aromaticity

Benzene ($\text{C}_6\text{H}_6$) was first isolated by Michael Faraday in 1825 and its molecular formula suggested a highly unsaturated structure. Yet benzene does not behave like a typical alkene: it resists addition reactions and instead undergoes substitution. This unusual stability is the hallmark of aromaticity.

The key experimental observations that define aromatic behavior include:

• ● Unusual thermodynamic stability: The heat of hydrogenation of benzene (208 kJ/mol) is 150 kJ/mol less than expected for three isolated double bonds (3 × 119 = 357 kJ/mol), giving a resonance energy of ~150 kJ/mol.
• ● Equal bond lengths: All C–C bonds in benzene are 1.40 Å, intermediate between a single bond (1.54 Å) and a double bond (1.34 Å).
• ● Planar geometry: All 12 atoms lie in one plane, allowing continuous overlap of p-orbitals.
• ● Ring current: In NMR spectroscopy, aromatic protons appear at unusually high chemical shifts ($\delta \approx 6.5$–$8.5$ ppm) due to the induced ring current from circulating $\pi$ electrons.

2. Derivation: Aromaticity and Hückel's Rule

2.1 Molecular Orbital Theory of Cyclic Systems

Consider a cyclic, planar, fully conjugated system of $N$ carbon atoms, each contributing one p-orbital. In the Hückel approximation, the secular determinant for a monocyclic ring yields the $\pi$ molecular orbital energies:

where $\alpha$ is the Coulomb integral (energy of an electron in an isolated p-orbital) and $\beta$ is the resonance integral (interaction energy between adjacent p-orbitals;$\beta < 0$ for bonding). The quantum number $k$ ranges over values that give$N$ distinct energy levels.

2.2 Deriving the Energy Spectrum

The Hückel Hamiltonian for a cyclic system of $N$ atoms is a circulant matrix:

Since this is a circulant matrix, its eigenvectors are the discrete Fourier modes$c_j = \frac{1}{\sqrt{N}} e^{2\pi i k j/N}$, and the eigenvalues are obtained by Fourier transforming the first row. Specifically, defining $x = (E - \alpha)/\beta$, the secular equation factorizes as:

giving $E_k = \alpha + 2\beta\cos(2\pi k/N)$. The lowest energy level ($k = 0$) is always non-degenerate with energy $E_0 = \alpha + 2\beta$. Higher levels come in degenerate pairs (for $k$ and $N - k$) except when $N$ is even, where$k = N/2$ is also non-degenerate.

2.3 The Frost Circle (Mnemonic)

A convenient mnemonic: inscribe a regular $N$-gon in a circle of radius $2|\beta|$ with one vertex pointing down. The heights of the vertices give the orbital energies (relative to$\alpha$). The horizontal midline corresponds to $E = \alpha$ (non-bonding level).

2.4 Application: Benzene (N = 6, Aromatic)

For benzene with $N = 6$, the energy levels are:

Benzene has 6 $\pi$ electrons. Filling from the bottom: 2 electrons in $E_0$ and 4 electrons in the degenerate pair $E_{\pm 1}$. All bonding orbitals are completely filled, giving a closed-shell configuration. The total $\pi$ energy is:

For three isolated ethylene units: $E_{\text{loc}} = 3 \times 2(\alpha + \beta) = 6\alpha + 6\beta$. The delocalization energy (resonance stabilization) is:

Since $\beta < 0$, this is a stabilization of $2|\beta| \approx 150$ kJ/mol.

2.5 Application: Cyclobutadiene (N = 4, Antiaromatic)

For cyclobutadiene with $N = 4$ and 4 $\pi$ electrons:

Filling 4 electrons: 2 in $E_0$ and 2 in the degenerate non-bonding pair $E_{\pm 1}$. By Hund's rule, these two electrons occupy separate orbitals with parallel spins, giving an open-shell (diradical) configuration. Total energy:$E_\pi = 2(\alpha + 2\beta) + 2\alpha = 4\alpha + 4\beta$. For two isolated ethylene units: $4\alpha + 4\beta$ — no delocalization energy. The system distorts from square to rectangular geometry (Jahn–Teller distortion) and is extremely reactive.

2.6 Application: Cyclooctatetraene (N = 8)

Cyclooctatetraene ($\text{C}_8\text{H}_8$, 8 $\pi$ electrons) has $4n$ electrons ($n = 2$). Like cyclobutadiene, the planar form would have a degenerate pair of non-bonding orbitals partially occupied. Instead, cyclooctatetraene adopts a tub-shaped, non-planar conformation with alternating single and double bonds — it behaves as a simple polyene, not an aromatic compound.

3. Derivation: Electrophilic Aromatic Substitution (EAS)

3.1 Why Substitution Instead of Addition?

Alkenes undergo electrophilic addition because the product retains $\sigma$ bonds that are comparably stable. For benzene, however, addition would destroy the aromatic$\pi$ system and its ~150 kJ/mol resonance stabilization. In contrast, substitution restores the aromatic ring after the reaction. Thermodynamically:

because the product of substitution is still aromatic, while the product of addition is not.

3.2 General Mechanism

The EAS mechanism proceeds in two steps:

3.3 Energy Profile

The reaction coordinate diagram for EAS shows:

• ● A high activation energy for Step 1 ($\Delta G_1^\ddagger$), reflecting the energetic cost of partially breaking aromaticity to form the $\sigma$ complex.
• ● The $\sigma$ complex sits in a shallow energy well — a reactive intermediate.
• ● A low activation energy for Step 2 ($\Delta G_2^\ddagger \ll \Delta G_1^\ddagger$), since the product recovers full aromaticity.

The rate law for EAS is:

This second-order kinetics confirms that the rate-determining step involves both the aromatic substrate and the electrophile.

3.4 Common EAS Reactions

4. Derivation: Directing Effects in EAS

4.1 The Origin of Regioselectivity

When a substituted benzene $\text{C}_6\text{H}_5\text{G}$ undergoes EAS, the incoming electrophile can attack at three distinct positions relative to the existing substituent $\text{G}$: ortho (1,2), meta (1,3), or para (1,4). The product distribution depends on the nature of $\text{G}$.

The directing effect is determined by the relative stability of the$\sigma$-complex intermediates for ortho, meta, and para attack. Since Step 1 is rate-determining, the pathway with the most stable (lowest energy) intermediate will be favored (Hammond's postulate).

4.2 Ortho/Para Directors (Electron-Donating Groups)

Electron-donating groups (EDGs) direct incoming electrophiles to the ortho and para positions. These groups stabilize the $\sigma$ complex through resonance donation of electrons into the ring.

Common ortho/para directors:

• ● Strong activators (lone pair donors): $\text{-NH}_2$, $\text{-NHR}$, $\text{-NR}_2$, $\text{-OH}$, $\text{-OR}$
• ● Moderate activators: $\text{-NHCOR}$, $\text{-OCOR}$
• ● Weak activators (hyperconjugation): $\text{-R}$ (alkyl groups)
• ● Deactivating but ortho/para directing: $\text{-F}$, $\text{-Cl}$, $\text{-Br}$, $\text{-I}$ (halogens: lone pair donation wins over inductive withdrawal for directing, but inductive effect dominates for rate)

4.3 Meta Directors (Electron-Withdrawing Groups)

Electron-withdrawing groups (EWGs) direct electrophiles to the meta position. These groups destabilize the $\sigma$ complex more strongly at ortho and para than at meta.

Common meta directors (all are also deactivating):

• ● $\text{-NO}_2$, $\text{-CN}$, $\text{-SO}_3\text{H}$
• ● $\text{-COOH}$, $\text{-COOR}$, $\text{-COR}$, $\text{-CHO}$
• ● $\text{-CF}_3$, $\overset{+}{\text{N}}\text{R}_3$

4.4 Activating vs. Deactivating Groups

Activating groups increase the rate of EAS relative to benzene by stabilizing the $\sigma$ complex (and the transition state leading to it). Deactivating groups decrease the rate by destabilizing the $\sigma$ complex. The key quantitative relationship uses the partial rate factor $f$:

where $\sigma$ is the Hammett substituent constant and $\rho$ is the reaction constant (negative for EAS, since electron density aids the reaction). A more negative $\sigma$ (EDG) gives a faster reaction; a more positive $\sigma$ (EWG) gives a slower one.

5. Derivation: Friedel-Crafts Reactions

5.1 Friedel-Crafts Alkylation

In Friedel-Crafts alkylation, a Lewis acid catalyst ($\text{AlCl}_3$) generates a carbocation (or polarized complex) from an alkyl halide, which then acts as the electrophile in EAS:

Limitations of Friedel-Crafts alkylation:

• ● Polyalkylation: The alkyl group is an activating, ortho/para-directing group. The product is more reactive than the starting material, leading to multiple substitutions.
• ● Carbocation rearrangement: Primary carbocations rearrange to more stable secondary or tertiary carbocations via 1,2-hydride or 1,2-methyl shifts.
• ● Does not work with deactivated rings: Rings bearing strong EWGs ($\text{-NO}_2$, $\text{-COR}$) or amines (which complex with $\text{AlCl}_3$) do not undergo Friedel-Crafts reactions.

5.2 Friedel-Crafts Acylation

Acylation uses an acyl chloride ($\text{RCOCl}$) with $\text{AlCl}_3$ to generate a stabilized acylium ion:

The acylium ion is resonance-stabilized (the triple bond contributes) and does not rearrange. The overall reaction:

5.3 Why Acylation Avoids Polyalkylation

The acyl group ($\text{-COR}$) introduced onto the ring is an electron-withdrawing, meta-directing, deactivating group. Therefore, the product is less reactive than the starting benzene — preventing a second substitution. Furthermore, the ketone product complexes with $\text{AlCl}_3$ (requiring >1 equivalent of the Lewis acid), further deactivating the ring.

5.4 Acylation-Reduction Strategy

To introduce a straight-chain alkyl group without rearrangement or polyalkylation, a two-step approach is used — Friedel-Crafts acylation followed by reduction of the carbonyl:

6. Derivation: Nucleophilic Aromatic Substitution

6.1 The SN Ar Mechanism (Addition-Elimination)

Unlike EAS, nucleophilic aromatic substitution (SNAr) involves attack by a nucleophile on an aromatic ring. This requires special activation: electron-withdrawing groups (typically ortho and/or para to the leaving group) that stabilize the anionic intermediate.

The requirements for SNAr:

• ● Strong EWGs ortho and/or para to the leaving group ($\text{-NO}_2$ is most effective; 2,4-dinitro and 2,4,6-trinitro substrates react readily)
• ● Good leaving group: $\text{F} > \text{Cl} > \text{Br} > \text{I}$ (note: opposite to SN2! Because bond breaking is not rate-determining; C-F bond is most polarized, facilitating nucleophilic attack)
• ● Strong nucleophile: $\text{OH}^-$, $\text{OR}^-$, $\text{NH}_2^-$, $\text{NH}_2\text{R}$

6.2 Why EWGs Are Required

The Meisenheimer complex bears a formal negative charge that must be stabilized. When EWGs are positioned ortho or para to the leaving group, the negative charge can be delocalized directly onto these groups through resonance:

Each additional nitro group ortho or para to the leaving group accelerates the reaction by roughly $10^4$–$10^7$ fold, because it provides additional resonance stabilization to the Meisenheimer complex.

6.3 The Benzyne Mechanism (Elimination-Addition)

When no EWGs are present, nucleophilic substitution on an aryl halide can still occur via a completely different pathway — the benzyne (aryne) mechanism:

7. Applications of Aromatic Chemistry

Summary: Key Aromatic Reactions at a Glance

All reactions above are examples of electrophilic aromatic substitution (EAS). Nucleophilic aromatic substitution (SNAr) and the benzyne mechanism are discussed in Section 6. Sulfonation is notable for being reversible — desulfonation occurs under dilute acid conditions, making the sulfonic acid group a useful blocking group in multi-step synthesis.

8. Historical Context

9. Python Simulation: Hückel MO Calculation

The following Python code performs Hückel molecular orbital calculations for benzene, naphthalene, and pyridine (a heterocyclic aromatic with nitrogen). It diagonalizes the Hückel Hamiltonian to obtain orbital energies and coefficients, computes delocalization energies, and generates energy level diagrams and coefficient maps. Cyclobutadiene and cyclooctatetraene (antiaromatic systems) are also analyzed for comparison.

For pyridine, the nitrogen atom is modeled with a modified Coulomb integral:$\alpha_{\text{N}} = \alpha + 0.5\beta$, reflecting the greater electronegativity of nitrogen. The MO energy formula for a general system is found by diagonalizing the Hückel matrix numerically using numpy: