Skip to content

Complete Intrinsic Reaction Coordinate Algorithm

An Intrinsic Reaction Coordinate (IRC) calculation determines the Minimum Energy Path (MEP) connecting a transition state with its corresponding reactants and products. Unlike a Transition State Search, which locates a first-order saddle point, the IRC calculation follows the reaction pathway away from that saddle point until stable molecular minima are reached.

The IRC algorithm combines concepts from electronic structure theory, geometry optimization, vibrational analysis, and reaction-path following to construct a continuous description of the chemical reaction.

Overall Workflow

The complete IRC procedure can be summarized as

Verified Transition State
           │
           ▼
Frequency Calculation
           │
           ▼
Imaginary Normal Mode
           │
           ▼
Choose Initial Direction
           │
           ▼
Transform to Mass-Weighted Coordinates
           │
           ▼
Predict Next IRC Point
           │
           ▼
Correct onto Minimum Energy Path
           │
           ▼
Calculate Energy & Gradient
           │
           ▼
Update Hessian
           │
           ▼
Converged to Minimum?
      │              │
     No             Yes
      │              │
      └──────────────┘
           │
           ▼
Reactant Minimum

Return to Transition State

           │
           ▼

Repeat Entire Procedure

           │
           ▼
Product Minimum

Step 1 — Verified Transition State

An IRC calculation always begins from a previously verified transition state.

Before starting the calculation,

Gaussian requires that

  • the Transition State optimization has converged,
  • exactly one imaginary frequency is present.

The transition state serves as the starting point for the reaction pathway.

📘 Transition State Search Mathematics

Step 2 — Determine the Reaction Coordinate

The imaginary vibrational mode obtained from the frequency calculation defines the initial reaction coordinate.

This normal mode indicates the direction of negative curvature on the Potential Energy Surface.

Gaussian follows this direction during the first IRC step.

Frequency Calculation

📘 Frequency Calculation Mathematics

Step 3 — Transform to Mass-Weighted Coordinates

Before following the reaction path,

Gaussian transforms the Cartesian coordinates into mass-weighted coordinates.

This transformation

  • accounts for differences in atomic masses,
  • produces a physically meaningful reaction coordinate,
  • improves numerical stability.

The IRC is therefore defined in mass-weighted coordinate space rather than ordinary Cartesian space.

Step 4 — Predictor Step

Using the local gradient,

the Hessian,

and the current reaction direction,

Gaussian predicts the location of the next IRC point.

Current IRC Point
        │
        ▼
Predict Next Point

This predicted geometry is only an approximation.

Step 5 — Corrector Step

The predicted geometry is then corrected so that it lies on the true Minimum Energy Path.

During this stage,

Gaussian minimizes the energy in all directions perpendicular to the reaction coordinate while preserving movement along the reaction path.

Predicted Geometry
        │
        ▼
Geometry Correction
        │
        ▼
Correct IRC Point

The corrected geometry becomes the next point on the IRC.

Step 6 — Electronic Structure Calculation

At every IRC point,

Gaussian performs a complete electronic structure calculation.

This includes

  • SCF calculation,
  • electronic energy,
  • molecular orbitals,
  • electron density,
  • energy gradient.

The electronic energy is therefore evaluated continuously along the reaction pathway.

Single Point Energy

📘 Single Point Energy Mathematics

Step 7 — Hessian Update

As the molecular geometry changes,

the local curvature of the Potential Energy Surface also changes.

Rather than computing a completely new Hessian at every IRC point,

Gaussian updates the existing Hessian.

This updated curvature information improves the prediction of subsequent IRC points while reducing computational cost.

Step 8 — Continue Along the Reaction Path

The Predictor–Corrector cycle is repeated many times.

IRC Point
      │
      ▼
Predict
      │
      ▼
Correct
      │
      ▼
SCF Calculation
      │
      ▼
Update Hessian
      │
      ▼
Next IRC Point

Each cycle generates another optimized structure along the Minimum Energy Path.

Step 9 — Reach a Stable Minimum

Eventually,

the energy no longer decreases,

and the molecular geometry converges to a stable minimum.

This minimum corresponds to either

  • the reactant,
  • or the product.

One branch of the IRC calculation is now complete.

Step 10 — Repeat in the Opposite Direction

Gaussian returns to the transition state.

The entire procedure is repeated along the opposite branch of the reaction coordinate.

Reactants
      ▲
      │
Transition State
      │
      ▼
Products

Once both directions have converged,

the complete reaction pathway has been determined.

Complete Mathematical Flow

The complete mathematical framework can be summarized as

Verified Transition State
           │
           ▼
Imaginary Frequency
           │
           ▼
Reaction Coordinate
           │
           ▼
Mass-Weighted Coordinates
           │
           ▼
Predictor Step
           │
           ▼
Corrector Step
           │
           ▼
SCF Calculation
           │
           ▼
Energy & Gradient
           │
           ▼
Hessian Update
           │
           ▼
Next IRC Point
           │
           ▼
Repeat Until Minimum
           │
           ▼
Reactants

Return to Transition State

           │
           ▼

Repeat Entire Procedure

           │
           ▼
Products

Mathematical Foundations Used

The IRC algorithm combines concepts developed throughout the previous Gaussian documentation.

Single Point Energy

Electronic energies are evaluated at every IRC point.

📘 Single Point Energy Mathematics

Geometry Optimization

Each IRC point is obtained through a constrained optimization using gradients and Hessian information.

📘 Geometry Optimization Mathematics

Frequency Calculation

The imaginary vibrational mode defines the initial reaction coordinate.

📘 Frequency Calculation Mathematics

Transition State Search

The starting point of every IRC calculation is a verified transition state.

📘 Transition State Search Mathematics

Together, these topics provide the complete theoretical foundation for understanding how Gaussian constructs an Intrinsic Reaction Coordinate.

Summary

The Intrinsic Reaction Coordinate algorithm begins from a verified transition state and follows the Minimum Energy Path toward both the reactants and products. At each step, Gaussian transforms the molecular coordinates into mass-weighted coordinates, predicts the next point along the reaction path, corrects the geometry onto the Minimum Energy Path, performs an electronic structure calculation, and updates the local curvature information. Repeating this Predictor–Corrector cycle produces a continuous sequence of optimized molecular structures that describes the complete reaction mechanism.

Conclusion

The Intrinsic Reaction Coordinate calculation is the final stage in the characterization of a chemical reaction. While a Transition State Search identifies the highest-energy structure along the reaction pathway, the IRC algorithm demonstrates how the system moves from reactants to products through that transition state. By combining electronic structure theory, geometry optimization, vibrational analysis, and reaction-path following, Gaussian provides a complete and physically meaningful description of the reaction mechanism.