Intrinsic Reaction Coordinate (IRC)¶
The Intrinsic Reaction Coordinate (IRC) calculation follows the minimum-energy pathway connecting a transition state to its corresponding reactants and products on the Potential Energy Surface (PES). While a Transition State Search locates the highest-energy structure along a reaction pathway, an IRC calculation determines whether that transition state actually connects the expected molecular structures.
In Gaussian, an IRC calculation starts from a previously optimized transition state and follows the reaction coordinate in both the forward and reverse directions. At each point along the pathway, the molecular geometry is optimized while remaining on the minimum-energy path. The result is a continuous sequence of structures that illustrates how the reaction proceeds from reactants, through the transition state, to products.
Why Perform an IRC Calculation?¶
A successful Transition State Search confirms the existence of a first-order saddle point through the presence of one imaginary vibrational frequency. However, this alone does not prove that the transition state corresponds to the desired chemical reaction.
An IRC calculation answers the important question:
Does this transition state connect the intended reactants and products?
By tracing the reaction pathway downhill in both directions, the IRC calculation verifies the connectivity of the transition state and provides confidence that the identified reaction mechanism is correct.
What Does an IRC Calculation Compute?¶
During an IRC calculation, Gaussian repeatedly
- moves a small distance along the reaction coordinate,
- performs a constrained geometry optimization,
- calculates the electronic energy,
- records the optimized geometry,
- repeats the procedure until a stable minimum is reached.
The complete pathway is therefore represented as a sequence of optimized molecular structures connecting the transition state with the surrounding minima.
Typical Workflow¶
An IRC calculation is usually performed after a successful Transition State Search.
Reactants
│
▼
Geometry Optimization
│
▼
Transition State Search
│
▼
Frequency Calculation
│
▼
Exactly One Imaginary Frequency?
│
Yes
│
▼
IRC Calculation
│
▼
Reaction Pathway
│
├────────► Reactants
│
└────────► Products
Thus, the IRC calculation serves as the final verification of the proposed reaction mechanism.
Information Obtained from an IRC Calculation¶
An IRC calculation provides several important pieces of information:
- The minimum-energy reaction pathway (MEP).
- The sequence of molecular geometries along the reaction.
- Changes in bond lengths and bond angles during the reaction.
- The evolution of the electronic energy.
- Confirmation that the transition state connects the correct reactants and products.
- A detailed picture of the reaction mechanism.
These results are invaluable for understanding how chemical transformations occur at the molecular level.
When Should an IRC Calculation Be Used?¶
An IRC calculation is recommended whenever a transition state has been successfully located and the reaction mechanism needs to be verified.
Typical applications include
- proton transfer reactions,
- hydrogen migration,
- bond-breaking and bond-forming processes,
- molecular rearrangements,
- isomerization reactions,
- catalytic reaction mechanisms.
For most computational studies, an IRC calculation is considered the final step after Transition State optimization.
Prerequisites¶
Before performing an IRC calculation, the following steps should already be complete:
- Geometry Optimization of the reactant.
- Transition State Search.
- Frequency Calculation confirming exactly one imaginary frequency.
Without a valid transition state, an IRC calculation is generally not meaningful.
Related Sections¶
Transition State Search¶
Frequency Calculation¶
Learning Objectives¶
After completing this section, you will understand
- how to prepare an IRC calculation in Gaussian,
- the purpose of the important IRC keywords,
- how to interpret the output,
- how to verify the reaction pathway,
- how Gaussian follows the minimum-energy path,
- the mathematical principles behind the IRC algorithm.
Summary¶
An Intrinsic Reaction Coordinate (IRC) calculation follows the minimum-energy pathway connecting a transition state with its corresponding reactants and products. By tracing the reaction coordinate in both directions, Gaussian verifies that the optimized transition state describes the intended chemical transformation. Because it establishes the connectivity of the reaction pathway, the IRC calculation is widely regarded as the final validation step in computational studies of chemical reaction mechanisms.
-
Download Input
Complete Gaussian input used in this tutorial.
-
Download Output
Complete output file discussed below.