Skip to content

Mass-Weighted Coordinates

When following a reaction pathway, Gaussian does not work directly with ordinary Cartesian coordinates. Instead, the molecular coordinates are transformed into mass-weighted coordinates, where each atomic displacement is scaled according to the mass of the corresponding atom.

This transformation ensures that the reaction pathway reflects the true motion of the nuclei and is independent of the arbitrary choice of Cartesian coordinates.

Mass-weighted coordinates are therefore fundamental to the mathematical definition of the Intrinsic Reaction Coordinate (IRC).

Why are Mass-Weighted Coordinates Needed?

Atoms in a molecule have very different masses.

For example,

  • Hydrogen is very light.
  • Carbon is much heavier.
  • Oxygen is heavier still.

If Cartesian coordinates were used directly, a small displacement of a hydrogen atom would contribute just as much as the same displacement of an oxygen atom.

Physically, however, moving a heavy atom requires much more energy than moving a light atom.

Mass-weighted coordinates account for this difference automatically.

Mathematical Definition

Suppose an atom has Cartesian coordinate

\[ x_i \]

and atomic mass

\[ m_i. \]

The corresponding mass-weighted coordinate is

\[ q_i=\sqrt{m_i}\,x_i. \]

Similarly,

\[ q_y=\sqrt{m_i}\,y_i, \]
\[ q_z=\sqrt{m_i}\,z_i. \]

Every Cartesian coordinate is therefore multiplied by the square root of the corresponding atomic mass.

Physical Interpretation

The transformation changes the coordinate system so that nuclear motion is measured according to its physical importance.

Cartesian Coordinates

x, y, z

        │
        ▼

Multiply by √mass

        │
        ▼

Mass-Weighted Coordinates

As a result,

  • heavy atoms contribute more to the coordinate system,
  • light atoms contribute less,
  • the reaction coordinate becomes physically meaningful.

Effect on the Reaction Path

The IRC is defined in mass-weighted coordinate space.

Instead of simply connecting molecular geometries,

Gaussian follows the path corresponding to the smallest change in the mass-weighted nuclear motion.

Transition State
        │
        ▼
Mass-Weighted Coordinates
        │
        ▼
Reaction Path
        │
        ▼
Reactants / Products

This produces a reaction pathway that better represents the actual dynamics of the nuclei.

Connection with Vibrational Analysis

Mass-weighted coordinates are also used in vibrational frequency calculations.

During frequency analysis,

the Hessian matrix is transformed into mass-weighted coordinates before being diagonalized to obtain

  • normal modes,
  • vibrational frequencies.

The IRC calculation uses the same coordinate system because the reaction coordinate originates from the imaginary normal mode of the transition state.

Advantages of Mass-Weighted Coordinates

Using mass-weighted coordinates provides several important benefits.

  • Nuclear motion has a clear physical meaning.
  • Heavy and light atoms contribute appropriately.
  • The reaction coordinate becomes independent of the Cartesian coordinate system.
  • Numerical stability of the IRC algorithm is improved.
  • The reaction path is consistent with vibrational analysis.

Relation to Previous Mathematics

The use of mass-weighted coordinates was introduced during the derivation of molecular vibrations in the Frequency Calculation section.

Frequency Calculation

📘 Frequency Calculation Mathematics

The present chapter extends that concept by showing how the same coordinate transformation is used to define and follow the Intrinsic Reaction Coordinate.

Summary

Mass-weighted coordinates transform ordinary Cartesian coordinates by scaling each atomic displacement with the square root of the atomic mass. This transformation ensures that the Intrinsic Reaction Coordinate represents the physically correct motion of the nuclei rather than simple geometric displacements. Because the IRC is defined in mass-weighted coordinate space, the resulting reaction pathway accurately reflects the minimum-energy motion connecting the transition state with the reactants and products.

Next Section

The next chapter introduces the Predictor–Corrector Algorithm, which explains how Gaussian predicts successive IRC points and then corrects them so that they remain on the true Minimum Energy Path.