How a Frequency Calculation Works¶
A Frequency Calculation determines how a molecule vibrates about its equilibrium geometry. Unlike a Geometry Optimization, which searches for the lowest-energy structure, a Frequency Calculation assumes that the molecular geometry has already been optimized and investigates how the atoms respond to small displacements from this equilibrium position.
The calculation provides information about molecular vibrations, infrared (IR) spectra, and thermodynamic properties. It also confirms whether the optimized structure corresponds to a true minimum or a transition state on the Potential Energy Surface (PES).
Why Perform a Frequency Calculation After Optimization?¶
Once a geometry optimization has converged, the forces acting on the atoms are essentially zero. However, this alone does not guarantee that the structure is a true minimum on the Potential Energy Surface.
A Frequency Calculation examines the curvature of the energy surface around the optimized geometry.
Depending on this curvature, the optimized structure may be
- a stable minimum,
- a transition state,
- or a higher-order saddle point.
This is determined by analyzing the molecular vibrational frequencies.
The Basic Idea¶
Imagine that each atom in a molecule is connected to its neighboring atoms by tiny springs.
If one atom is displaced slightly from its equilibrium position, the surrounding atoms respond, causing the entire molecule to vibrate.
A Frequency Calculation determines
- how the atoms move,
- how fast they vibrate,
- and how much energy is associated with each vibration.
Starting from the Optimized Geometry¶
The calculation begins with an optimized molecular structure.
Since the geometry is already at (or very near) a stationary point on the Potential Energy Surface, only very small atomic displacements are required for the vibrational analysis.
Small Atomic Displacements¶
Gaussian moves each atom by a very small distance in the positive and negative directions.
The displacement is extremely small—typically only a fraction of an angstrom—so that the molecule remains close to its equilibrium geometry.
Measuring the Change in Energy¶
After each displacement, Gaussian performs an electronic structure calculation to determine how the molecular energy changes.
By repeating this process for all atomic coordinates, Gaussian determines how the energy varies with nuclear motion.
Constructing the Hessian Matrix¶
The energy changes obtained from these small displacements are used to construct the Hessian matrix.
The Hessian contains the second derivatives of the molecular energy with respect to the nuclear coordinates.
Conceptually, it describes the curvature of the Potential Energy Surface around the optimized geometry.
A steep curvature corresponds to a stiff vibration, while a shallow curvature corresponds to a softer vibration.
Determining the Normal Modes¶
Once the Hessian matrix has been constructed, Gaussian analyzes the collective motions of the atoms.
Instead of considering the movement of each atom independently, the program identifies combinations of atomic motions that occur together.
These collective vibrations are called normal modes.
Each normal mode has
- a characteristic vibrational frequency,
- a specific pattern of atomic motion,
- and an associated vibrational energy.
Calculating Vibrational Frequencies¶
The vibrational frequencies are obtained from the curvature of the Potential Energy Surface.
Generally,
- steeper potential wells produce higher frequencies,
- shallower potential wells produce lower frequencies.
For example,
- O–H bond stretching vibrations occur at relatively high frequencies,
- bending and torsional motions occur at much lower frequencies.
Infrared Intensities¶
For each normal mode, Gaussian also evaluates how the molecular dipole moment changes during the vibration.
If the vibration changes the dipole moment of the molecule, the mode is infrared active and can absorb infrared radiation.
If the dipole moment remains unchanged, the vibration is infrared inactive.
Thermodynamic Properties¶
The vibrational frequencies are then used to calculate several thermodynamic quantities.
These include
- Zero-Point Energy (ZPE),
- Thermal correction to energy,
- Enthalpy,
- Entropy,
- Gibbs free energy,
- Heat capacity.
These properties are evaluated using statistical thermodynamics at the specified temperature and pressure (298.15 K and 1 atm by default).
Verifying the Optimized Structure¶
One of the most important purposes of a Frequency Calculation is to determine the nature of the optimized structure.
The interpretation is straightforward:
| Imaginary Frequencies | Interpretation |
|---|---|
| 0 | Stable minimum |
| 1 | Transition state |
| 2 or more | Higher-order saddle point or unconverged structure |
Therefore, a successful geometry optimization is usually followed by a Frequency Calculation to verify that the optimized structure corresponds to the desired stationary point.
Overall Workflow¶
The complete process can be summarized as
Optimized Geometry
│
▼
Small Atomic Displacements
│
▼
Electronic Structure Calculations
│
▼
Construct Hessian Matrix
│
▼
Determine Normal Modes
│
▼
Calculate Vibrational Frequencies
│
▼
Compute IR Intensities
│
▼
Evaluate Thermodynamic Properties
│
▼
Identify Stationary Point
Each step builds upon the previous one to provide a comprehensive description of the molecule's vibrational behavior.
Summary¶
A Gaussian Frequency Calculation analyzes the behavior of a molecule near its optimized geometry by introducing small atomic displacements and evaluating how the molecular energy changes. From these energy variations, Gaussian constructs the Hessian matrix, determines the normal modes of vibration, and calculates the corresponding vibrational frequencies. These frequencies are then used to predict infrared spectra, compute thermodynamic properties, and verify whether the optimized structure is a stable minimum or a transition state. As a result, Frequency Analysis is an essential component of most quantum chemical investigations.
Next Section¶
The next part of this documentation presents the Mathematics of Frequency Calculations, beginning with the Harmonic Approximation, which forms the theoretical foundation of molecular vibrational analysis.