How IRC Works¶
An Intrinsic Reaction Coordinate (IRC) calculation determines the minimum-energy pathway (MEP) connecting a transition state with its corresponding reactants and products. Rather than searching for a stationary point, as in Geometry Optimization or Transition State Search, the IRC calculation begins from an already verified transition state and follows the reaction coordinate downhill in both directions.
The result is a continuous sequence of molecular structures that describes how the chemical reaction proceeds on the Potential Energy Surface (PES).
The Basic Idea¶
A transition state represents the highest-energy point along the reaction pathway.
Once this point has been located, there are two possible directions in which the system can evolve:
- toward the reactants,
- toward the products.
The IRC calculation follows both directions by taking many small steps away from the transition state while remaining on the minimum-energy path.
Instead of finding a new stationary point, Gaussian traces the pathway connecting existing stationary points.
Step 1 — Start from the Transition State¶
The calculation begins with a fully optimized transition-state geometry.
Before an IRC calculation is performed,
the transition state should satisfy two conditions:
- the optimization has converged,
- exactly one imaginary vibrational frequency is present.
The imaginary frequency defines the reaction coordinate that Gaussian will follow.
Step 2 — Determine the Initial Direction¶
Gaussian identifies the imaginary normal mode obtained from the frequency calculation.
This normal mode indicates the direction of negative curvature on the Potential Energy Surface.
The IRC calculation starts by moving a small distance along this direction.
This displacement creates the first IRC point.
Step 3 — Optimize the New Geometry¶
The displaced geometry is not yet on the exact reaction pathway.
Gaussian therefore performs a constrained optimization to locate the nearest point on the minimum-energy path.
At this stage,
Gaussian calculates
- the electronic energy,
- the gradient,
- the local curvature,
and adjusts the geometry accordingly.
Step 4 — Continue Along the Path¶
After one IRC point has been located,
Gaussian repeats the same procedure.
Each cycle consists of
- moving a small distance along the reaction coordinate,
- optimizing the displaced geometry,
- storing the new structure,
- calculating the electronic energy,
- repeating the process.
In this way,
the reaction pathway is constructed point by point.
Step 5 — Reach the Reactant Minimum¶
As the calculation proceeds,
the molecular energy gradually decreases.
Eventually,
Gaussian reaches a geometry where
- the energy no longer decreases,
- the forces become negligible,
- a stable minimum is obtained.
This structure corresponds to one side of the reaction.
Step 6 — Repeat in the Opposite Direction¶
After completing one direction,
Gaussian automatically returns to the transition state.
The procedure is then repeated in the opposite direction.
The result is a complete reaction pathway connecting both minima.
Energy Changes Along the IRC¶
Unlike a Geometry Optimization,
an IRC calculation does not simply minimize the energy.
Instead,
the energy decreases smoothly away from the transition state in both directions.
The transition state is therefore the highest point,
while the reactants and products correspond to lower-energy minima.
Structural Changes Along the IRC¶
Each IRC point corresponds to an optimized molecular geometry.
As the calculation progresses,
the molecular structure changes continuously.
Typical structural changes include
- bond breaking,
- bond formation,
- proton transfer,
- hydrogen migration,
- molecular rearrangement.
Because every geometry is stored,
the entire reaction can be animated using molecular visualization software such as GaussView.
Reaction Path Verification¶
The primary purpose of an IRC calculation is verification.
After the calculation finishes,
the user confirms that
correspond to the expected molecular structures.
If either endpoint is incorrect,
the transition state does not describe the intended reaction.
Overall Workflow¶
The complete IRC procedure can be summarized as
Optimized Transition State
│
▼
Identify Imaginary Mode
│
▼
Small Displacement
│
▼
Optimize Geometry
│
▼
IRC Point
│
▼
Repeat
│
▼
Reactant Minimum
│
│
Return to Transition State
│
▼
Repeat in Opposite Direction
│
▼
Product Minimum
Why is IRC Important?¶
A successful Transition State Search identifies a first-order saddle point.
However,
only an IRC calculation can demonstrate that the transition state actually connects the desired reactants and products.
Without an IRC calculation,
the reaction mechanism remains uncertain.
For this reason,
IRC calculations are widely regarded as the final verification step in computational reaction mechanism studies.
Relation to Previous Sections¶
An IRC calculation combines concepts developed throughout the previous Gaussian calculations.
Single Point Energy¶
The electronic energy is evaluated at every point along the reaction pathway.
Geometry Optimization¶
Each IRC point is obtained through a constrained optimization on the Potential Energy Surface.
Frequency Calculation¶
The imaginary vibrational mode obtained from the frequency calculation defines the initial reaction coordinate.
Transition State Search¶
The IRC calculation always starts from a previously verified transition state.
Summary¶
An Intrinsic Reaction Coordinate calculation begins with an optimized transition state and follows the reaction coordinate toward both the reactants and products. Gaussian repeatedly takes small steps along the reaction coordinate, performs constrained geometry optimizations, and records the corresponding molecular structures and energies. The resulting sequence of IRC points defines the minimum-energy pathway of the reaction and confirms that the transition state connects the correct molecular minima. This makes the IRC calculation the final and most reliable verification of a proposed reaction mechanism.