A potential energy surface (PES) maps the potential energy of a molecular system as a function of its geometric coordinates (bond lengths, bond angles, dihedral angles). Think of it as a landscape where every point corresponds to a specific arrangement of atoms, and the "elevation" at that point is the system's potential energy.

The topography of this landscape tells you a great deal about the chemistry:

Valleys (minima) correspond to stable species: reactants, products, or intermediates.
Mountain passes (saddle points) correspond to transition states.
Ridges and slopes reveal energy barriers and alternative reaction pathways.

For a system with $N$ atoms, the full PES has $3N - 6$ internal degrees of freedom ( $3N - 5$ for linear molecules). That means even a triatomic system like $H_2 + H$ already requires a three-dimensional surface. For larger molecules, you can't visualize the full PES directly, so we work with slices or projections onto the most chemically relevant coordinates.

Reaction pathways and transition states

A reaction pathway is the minimum energy path (MEP) connecting reactants to products on the PES. It traces the route a system is most likely to follow as it transforms, because deviations from this path cost extra energy.

The transition state sits at the saddle point along this path. It's the highest-energy point the system must pass through, and it has a very specific mathematical character: a maximum in one direction (the reaction coordinate) and a minimum in all perpendicular directions.

The energy gap between the reactants and the transition state defines the activation energy $E_a$ . This quantity controls the reaction rate through the Arrhenius equation:

$k = A \exp\!\left(-\frac{E_a}{RT}\right)$

where $A$ is the pre-exponential factor, $R$ is the gas constant, and $T$ is the absolute temperature. A higher $E_a$ means a slower reaction at a given temperature.

Potential energy diagrams for reactions

Graphical representations and insights, physical chemistry - How should I think about potential energy surfaces for reactions of complex ...

Simplified representations of potential energy surfaces

A potential energy diagram is a 2D cross-section of the full PES, plotting energy along the vertical axis against a single reaction coordinate on the horizontal axis. The reaction coordinate tracks the progress of the reaction and is typically defined by the geometric parameter that changes most during the transformation (a bond length being broken or formed, a bond angle opening, etc.).

On this diagram:

Reactants and products appear as minima (valleys).
Transition states appear as maxima (peaks).
Intermediates, if present, appear as local minima between two maxima.

These diagrams sacrifice the multi-dimensional detail of the full PES but make it much easier to read off the key energetic quantities at a glance.

Constructing potential energy diagrams

To build a potential energy diagram, you calculate (or measure) the system's potential energy at a series of points along the chosen reaction coordinate. Here's the general process:

Choose the reaction coordinate. Pick the geometric parameter most relevant to the transformation (e.g., the breaking bond length in a dissociation).
Fix the reaction coordinate at a series of values spanning reactants to products.
Optimize all other coordinates at each fixed value, so you're tracing the minimum energy path.
Plot energy vs. reaction coordinate.

The energy difference between reactants and products gives you the reaction enthalpy: negative $\Delta H$ for exothermic reactions, positive $\Delta H$ for endothermic ones.

Some classic examples where this construction is straightforward:

Dissociation of $H_2$ : the reaction coordinate is the H–H bond length. The curve is essentially a Morse potential, with a minimum at the equilibrium bond length (~0.74 Å) and an asymptote at the dissociation energy (~436 kJ/mol).
Rotation around the C–C bond in ethane: the reaction coordinate is the dihedral angle. The diagram shows periodic minima at staggered conformations and maxima at eclipsed conformations, with a rotational barrier of about 12 kJ/mol.
Formation of $H_2O$ from $H_2 + O$ : the diagram shows the energy drop as new O–H bonds form, reflecting the strongly exothermic nature of the reaction.

Potential energy surfaces and reaction coordinates

Graphical representations and insights, Activation energy - Wikipedia

Relationship between potential energy surfaces and reaction coordinates

The reaction coordinate is a one-dimensional simplification of a multi-dimensional PES. Choosing it well is critical: a poor choice can obscure important features of the reaction or make the energy profile misleading.

In practice, the reaction coordinate is selected as the geometric parameter that changes most significantly between reactants and products. For a bond-breaking reaction, that's usually the length of the bond being broken. For a conformational change, it might be a dihedral angle.

The full PES contains all the information, but it lives in $3N - 6$ dimensions. The reaction coordinate projects this high-dimensional surface down to a single dimension, capturing the essential energetic story of the reaction while discarding perpendicular degrees of freedom.

Analyzing potential energy surfaces along reaction coordinates

When you trace the minimum energy path along the reaction coordinate, the resulting energy profile reveals:

The number of transition states the reaction passes through (each appears as a maximum).
The presence of intermediates (local minima between transition states).
Which step is rate-determining (the step with the highest-energy transition state relative to the preceding minimum).

Some instructive examples:

$S_N2$ reaction (e.g., $Cl^- + CH_3Br \rightarrow ClCH_3 + Br^-$ ): The reaction coordinate is roughly the distance between the nucleophile and the electrophilic carbon. The PES shows a single saddle point (the pentacoordinate transition state), with no intermediate. This is a concerted, one-step mechanism.
Diels-Alder cycloaddition (1,3-butadiene + ethene $\rightarrow$ cyclohexene): The reaction coordinate involves the simultaneous formation of two new C–C bonds. The PES shows a single, concerted transition state, consistent with a pericyclic mechanism.
Cyclohexane conformational analysis: The reaction coordinate is a combination of dihedral angles. The PES reveals the chair conformation as the global minimum, the twist-boat as a local minimum (~23 kJ/mol higher), and the half-chair as the transition state (~45 kJ/mol higher) connecting them.

Stationary points on potential energy surfaces

Types of stationary points

A stationary point is any point on the PES where the gradient of the energy with respect to all internal coordinates vanishes:

$\nabla E = \mathbf{0}$

This means the forces on all atoms are zero at that geometry. But not all stationary points are the same. To classify them, you examine the Hessian matrix (the matrix of second derivatives of the energy, $\frac{\partial^2 E}{\partial q_i \partial q_j}$ ). The eigenvalues of the Hessian tell you the curvature in each direction:

Minimum (0 negative eigenvalues): All eigenvalues are positive. The surface curves upward in every direction. This corresponds to a stable species (reactant, product, or intermediate).
First-order saddle point (1 negative eigenvalue): The surface curves downward in exactly one direction and upward in all others. This is a transition state.
Higher-order saddle points (2+ negative eigenvalues): These have multiple downward-curving directions. They don't typically correspond to simple transition states but can appear in complex PES landscapes with multiple competing pathways.

Significance of stationary points in reaction dynamics

Locating and characterizing stationary points is central to computational chemistry and reaction dynamics. Here's why each type matters:

Minima define the species you can actually observe. Their vibrational frequencies (obtained from the Hessian eigenvalues) determine IR and Raman spectra and contribute to thermodynamic quantities like zero-point energy and entropy.

First-order saddle points (transition states) control reaction rates. The curvature of the PES at the saddle point determines:

The imaginary frequency along the reaction coordinate (exactly one, corresponding to the single negative Hessian eigenvalue). This is the hallmark of a true transition state in a computational search.
The vibrational frequencies perpendicular to the reaction coordinate, which enter into transition state theory rate expressions.

A few examples that highlight the practical importance of stationary point analysis:

In the $S_N2$ reaction, the transition state saddle point is the rate-determining point. Its energy relative to reactants gives $E_a$ , and its geometry (trigonal bipyramidal at carbon) confirms the backside-attack mechanism.
In hydroboration-oxidation of alkenes, the PES has an intermediate minimum (the alkylborane) between two transition states, revealing a two-step mechanism.
In racemization of chiral molecules, the transition state geometry is typically a planar or symmetric configuration. Comparing $E_a$ for the normal substrate versus a deuterium-substituted analog reveals kinetic isotope effects, since the heavier isotope changes the zero-point energy at both the minimum and the saddle point.

2,589 studying →