Skip to content

Mathematics of Intrinsic Reaction Coordinate

The mathematics of an Intrinsic Reaction Coordinate (IRC) calculation builds upon concepts introduced in the previous Gaussian calculations. Unlike Geometry Optimization, which searches for a local minimum, or Transition State Search, which locates a first-order saddle point, an IRC calculation follows the minimum-energy pathway (MEP) connecting the transition state to the reactants and products.

Most of the mathematical framework—such as gradients, Hessians, and the Potential Energy Surface—has already been introduced in earlier chapters. This section therefore focuses only on the additional concepts required to understand how Gaussian traces a reaction pathway.

Prerequisite Mathematics

Before studying the IRC algorithm, readers should be familiar with the following topics.

Single Point Energy

Electronic energies are evaluated at every point along the reaction pathway.

📘 Single Point Energy Mathematics

Geometry Optimization

IRC points are obtained through constrained geometry optimizations using gradients and Hessian information.

📘 Geometry Optimization Mathematics

Frequency Calculation

The imaginary normal mode obtained from the frequency calculation defines the initial direction of the IRC path.

📘 Frequency Calculation Mathematics

The IRC calculation always begins from a verified transition state and follows the reaction coordinate identified during the Transition State Search.

📘 Transition State Search Mathematics

Mathematics Covered in this Section

Only the mathematical concepts unique to the IRC calculation are discussed.

  1. Minimum Energy Path (MEP)

    • Definition of the reaction pathway
    • Relationship to the Potential Energy Surface

  2. Reaction Path Following

    • Following the reaction coordinate from the transition state
    • Forward and reverse propagation

  3. Mass-Weighted Coordinates

    • Why Gaussian transforms Cartesian coordinates
    • Physical significance of mass weighting

  4. Predictor–Corrector Algorithm

    • Prediction of the next IRC point
    • Correction onto the true reaction path

  5. Step Size and Hessian Updates

    • Choosing appropriate IRC step lengths
    • Updating local curvature information

  6. Complete IRC Algorithm

    • Overall mathematical workflow used by Gaussian

Learning Objectives

After completing this section, you will understand

  • what defines the minimum-energy reaction pathway,
  • how Gaussian follows the reaction coordinate,
  • why mass-weighted coordinates are used,
  • how the Predictor–Corrector algorithm constructs the IRC,
  • how Gaussian maintains the reaction path during propagation,
  • and how the complete IRC algorithm combines these concepts to generate the full reaction pathway.

Summary

The mathematics of the IRC calculation extends the concepts developed in the previous Gaussian sections rather than introducing an entirely new optimization method. Building upon the Potential Energy Surface, gradients, Hessians, and transition-state theory, the IRC algorithm follows the minimum-energy pathway from the transition state to the reactants and products using specialized reaction-path following techniques. The chapters in this section focus exclusively on these additional mathematical ideas while referring to the earlier sections for the underlying theoretical framework.