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 Meta's Open Molecules 2025 (OMol25) & Universal Models for Atoms (UMA)

Exploring Meta's Open Molecules 2025 (OMol25) & Universal Models for Atoms (UMA)

A close look at the OMol25 dataset, the pre-trained eSEN and UMA models, and some thoughts about the future of NNP-accelerated atomistic simulation.
May 23, 2025 · Corin Wagen and Ari Wagen
Protein–Ligand Co-Folding

Protein–Ligand Co-Folding

folding vs co-folding; free open-source models; running Boltz-1 and Chai-1 through Rowan; decentralized data generation with Macrocosmos
May 9, 2025 · Spencer Schneider, Ari Wagen, and Corin Wagen
Rowan Research Spotlight: Turki Alturaifi

Rowan Research Spotlight: Turki Alturaifi

How Rowan helps researchers understand and optimize complex catalytic reactions.
May 7, 2025 · Corin Wagen
Partnering with Macrocosmos to Accelerate Next-Generation NNP Development

Partnering with Macrocosmos to Accelerate Next-Generation NNP Development

Starting today, Rowan is teaming up with Macrocosmos to accelerate the development of the next generation of NNPs through Bittensor Subnet 25 - Mainframe.
May 1, 2025 · Ari Wagen
Introducing Egret-1

Introducing Egret-1

trusting computation; speed vs accuracy; Egret-1, Egret-1e, and Egret-1t; benchmarks; speed on CPU and GPU; download Egret-1 or use it through Rowan
Apr 30, 2025 · Eli Mann, Corin Wagen, Jonathon Vandezande, Ari Wagen, and Spencer Schneider
Egret-1: Pretrained Neural Network Potentials For Efficient and Accurate Bioorganic Simulation

Egret-1: Pretrained Neural Network Potentials For Efficient and Accurate Bioorganic Simulation

Here, we present Egret-1, a family of large pre-trained NNPs based on the MACE architecture with general applicability to main-group, organic, and biomolecular chemistry.
Apr 30, 2025 · Elias L. Mann, Corin C. Wagen, Jonathon E. Vandezande, Arien M. Wagen, Spencer C. Schneider
Introducing Egret-1

Introducing Egret-1

Today, we're releasing Egret-1, a family of open-source NNPs for bioorganic simulation.
Apr 30, 2025 · Eli Mann, Corin Wagen, Jonathon Vandezande, Ari Wagen, and Spencer Schneider
Starling: Macroscopic pKa, logD, and Blood–Brain-Barrier Permeability

Starling: Macroscopic pKa, logD, and Blood–Brain-Barrier Permeability

microscopic vs. macroscopic pKa; Uni-pKa and Starling; microstate ensembles; logD and Kp,uu predictions
Apr 25, 2025 · Corin Wagen
Physics-Informed Machine Learning Enables Rapid Macroscopic pKa Prediction

Physics-Informed Machine Learning Enables Rapid Macroscopic pKa Prediction

Here we introduce Starling, a physics-informed neural network based on the Uni-pKa architecture trained to predict per-microstate free energies and compute macroscopic pKa values via thermodynamic ensemble modeling.
Apr 25, 2025 · Corin C. Wagen
Predicting Infrared Spectra and Orb-v3

Predicting Infrared Spectra and Orb-v3

light and its manifold interactions with matter; why IR spectroscopy is useful; predicting IR spectra through Rowan; Orb-v3
Apr 17, 2025 · Ari Wagen, Corin Wagen, and Jonathon Vandezande