Modeling Addition and Substitution Reactions

by Milca Pierre and Ari Wagen · Mar 20, 2025

This blog was written with friend of Rowan Milca Pierre. Milca is an undergraduate studying chemistry at the University of Massachusetts Lowell.

In this blog post, we'll be looking at an exercise from The Molecular Modeling Workbook for Organic Chemistry by Warren J. Hehre, Alan J. Shusterman, and Janet E. Nelson (Internet Archive, Amazon). Exercise 1 of chapter 13 looks at addition and substitution reactions, asking:

Unsaturated hydrocarbons undergo a variety of reactions. Experimentally, alkenes & alkynes undergo addition reactions whereas aromatic molecules, such as benzene, undergo substitution reactions instead. Why?

To answer this question, the workbook directs us to calculate the thermodynamics of each reaction. We'll be looking at two different sets of products (cyclohexane and benzene) and calculating the energy of the addition of bromine and corresponding substitution reaction for each.

Our goal is to determine which reactions are exothermic, which are endothermic, and by how much.

Figure 1: Addition and substitution reactions of cyclohexene and benzene with Br2

Figure 1: Addition and substitution reactions of cyclohexene and benzene with Br2

To quickly explore the thermodynamics of these reactions, we'll calculate the energy of each set of reactants and products using the AIMNet2 (ωB97M-D3) level of theory.

Cyclohexene + Br2

Here are the involved Rowan calculations:

Running these optimizations and plotting the thermodynamic change using Rowan's graph builder utility shows us that, although both the addition and substitution reactions are exothermic, the addition reaction is more exothermic. The addition reaction is, therefore, thermodynamically favored, accounting for the observation that "alkenes & alkynes undergo addition reactions."

Figure 2: Thermodynamics of the cyclohexene + Br2 addition and substitution reactions (transition states not modeled)

Figure 2: Thermodynamics of the cyclohexene + Br2 addition and substitution reactions (transition states not modeled)

Benzene + Br2

Here are the involved Rowan calculations:

The results of these calculations show that the addition reaction is endothermic, while the substitution reaction is exothermic. This means that the substitution reaction is thermodynamically favored, explaining the observation that "aromatic molecules, such as benzene, undergo substitution reactions."

Figure 3: Thermodynamics of the benzene + Br2 addition and substitution reactions (transition states not modeled)

Figure 3: Thermodynamics of the benzene + Br2 addition and substitution reactions (transition states not modeled)

To model these reactions' kinetics and predict which pathways would be kinetically favored, we would need to find and optimize each involved transition state. This would be a research project too large for a single blog post—these reactions are very solvent-dependent, may directly involve multiple solvent molecules or catalysts, and the transition-state structures are hard to find. We encourage the interested reader to see other sources looking at these reactions, including:

Figure 4: Computed transition-state structures of the benzene + Br2 substitution and addition reactions. Reproduced from figure 2 in Kong et. al (2011).

Figure 4: Computed transition-state structures of the benzene + Br2 substitution and addition reactions. Reproduced from figure 2 in Kong et. al (2011).

If you're interested in using Rowan to model your own reactions, you can do so through Rowan's web platform (it's free to make an account and get started). Happy computing!

Banner background image

What to Read Next

Exploring Protein–Ligand Binding-Affinity Prediction

Exploring Protein–Ligand Binding-Affinity Prediction

Trying a few modern ML-based approaches for predicting protein–ligand binding affinity.
Aug 20, 2025 · Ishaan Ganti
What Ishaan and Vedant Learned This Summer

What Ishaan and Vedant Learned This Summer

Reflections from two of our interns on their time at Rowan and a few things they learned.
Aug 15, 2025 · Ishaan Ganti and Vedant Nilabh
Projects: Organization, Sharing, and Saving Structures

Projects: Organization, Sharing, and Saving Structures

better organization through projects; saving structures; usage tracking; new conf. search features; second-order SCF; ex. API repo; SMILES imports; a guide to the pKa-perplexed; our inaugural demo day
Aug 14, 2025 · Ari Wagen, Spencer Schneider, Corin Wagen, and Jonathon Vandezande
Macroscopic and Microscopic pKa

Macroscopic and Microscopic pKa

Two different ways to calculate acidity, what they mean, and when to use them.
Aug 11, 2025 · Corin Wagen
Computational Chemistry in the Classroom

Computational Chemistry in the Classroom

chemical modeling; Diels–Alder; call for more labs
Jul 31, 2025 · Jonathon Vandezande and Isaiah Sippel
Modeling Thia-Michael Reactions

Modeling Thia-Michael Reactions

In which the addition of a thiolate to an enone proves to be unexpectedly difficult to model.
Jul 25, 2025 · Corin Wagen
API v2, New BDE Methods, MCP, And More

API v2, New BDE Methods, MCP, And More

new API philosophy; streamlined interfaces for workflows; using NNPs and g-xTB to predict bond strength; an MCP server; .sdf files; benchmarking protein–ligand interactions; Diels–Alder visualizations
Jul 21, 2025 · Spencer Schneider, Corin Wagen, Ari Wagen, Jonathon Vandezande, Ishaan Ganti, and Isaiah Sippel
ExpBDE54: A Slim Experimental Benchmark for Exploring the Pareto Frontier of Bond-Dissociation-Enthalpy-Prediction Methods

ExpBDE54: A Slim Experimental Benchmark for Exploring the Pareto Frontier of Bond-Dissociation-Enthalpy-Prediction Methods

ExpBDE54 is a benchmark dataset of experimental homolytic bond-dissociation enthalpies (BDEs) for 54 small molecules, used for benchmarking DFT, semiempirical methods, and NNPs.
Jul 17, 2025 · Jonathon E. Vandezande, Corin C. Wagen
Benchmarking Protein–Ligand Interaction Energy

Benchmarking Protein–Ligand Interaction Energy

How new low-cost computational methods perform on the PLA15 benchmark.
Jul 11, 2025 · Ishaan Ganti
Efficient Black-Box Prediction of Hydrogen-Bond-Donor and Acceptor Strength

Efficient Black-Box Prediction of Hydrogen-Bond-Donor and Acceptor Strength

Here, we report a robust black-box workflow for predicting site-specific hydrogen-bond basicity and acidity in organic molecules with minimal computational cost.
Jul 1, 2025 · Corin C. Wagen