The Harmonic Approximation¶

The exact Potential Energy Surface (PES) of a molecule is highly complex. As atoms move farther from their equilibrium positions, the energy changes in a nonlinear manner due to bond stretching, bond breaking, anharmonic interactions, and coupling between different vibrational modes.

Fortunately, molecular vibrations usually involve only very small displacements around the equilibrium geometry. Within this small region, the Potential Energy Surface can be approximated by a simple quadratic function.

This simplification is known as the Harmonic Approximation, and it forms the theoretical foundation of almost every vibrational frequency calculation performed by Gaussian.

Why is an Approximation Needed?¶

The true Potential Energy Surface is extremely complicated.

Energy
 ^
 |
 |                __
 |             __/  \__
 |          __/        \____
 |______ __/____________________

-------------------------------> Nuclear Coordinate

Its exact mathematical form is generally unknown and impossible to solve analytically for molecules containing many atoms.

To make vibrational calculations computationally feasible, Gaussian approximates only the small region surrounding the equilibrium geometry.

The Shape Near the Minimum¶

Near the equilibrium structure, the Potential Energy Surface becomes smooth and nearly parabolic.

Energy
 ^
 |
 |          *
 |        *   *
 |      *       *
 |_____*_________*________

        Equilibrium

-----------------------------> Nuclear Coordinate

Since atomic vibrations are very small, the molecule never moves far from this region during normal vibrational motion.

Therefore, the local parabolic approximation provides an excellent description of most molecular vibrations.

The Harmonic Oscillator Model¶

The Harmonic Approximation assumes that every chemical bond behaves like an ideal spring.

Atom A -----/////----- Atom B

When the atoms are displaced,

the bond stretches,
a restoring force appears,
the atoms oscillate around their equilibrium positions.

This model is identical to the classical mass–spring oscillator studied in introductory physics.

Characteristics of Harmonic Motion¶

Under the Harmonic Approximation,

the restoring force is directly proportional to the displacement.

Small Displacement
        │
        ▼
Small Restoring Force

Large Displacement
        │
        ▼
Large Restoring Force

This relationship produces smooth, periodic molecular vibrations.

Each vibrational mode behaves independently and possesses a well-defined vibrational frequency.

The Parabolic Potential¶

The Harmonic Approximation replaces the complicated Potential Energy Surface with a simple parabola.

Real Potential

Energy
 ^
 |
 |             __
 |          __/  \___
 |_______ _/___________

Approximate Harmonic Potential

Energy
 ^
 |
 |         *
 |       *   *
 |     *       *
 |___*___________*____

Near the equilibrium geometry,

the two curves are almost identical.

As the displacement becomes larger,

the approximation gradually becomes less accurate.

What Does Gaussian Assume?¶

During a standard Frequency Calculation,

Gaussian assumes that

all atomic displacements are small,
every vibrational mode behaves harmonically,
the Potential Energy Surface is locally quadratic.

Using these assumptions,

Gaussian can calculate

vibrational frequencies,
normal modes,
infrared intensities,
Zero-Point Energy,
thermodynamic properties.

without considering the full complexity of molecular motion.

Advantages of the Harmonic Approximation¶

The Harmonic Approximation offers several important advantages.

Greatly simplifies the mathematical treatment of vibrations.
Produces accurate frequencies for most molecules near equilibrium.
Makes vibrational calculations computationally efficient.
Allows analytical evaluation of many thermodynamic properties.
Forms the basis of nearly all routine quantum chemical frequency calculations.

Without this approximation, frequency calculations for large molecular systems would be prohibitively expensive.

Limitations of the Harmonic Approximation¶

Although highly successful, the Harmonic Approximation has several limitations.

It assumes that

chemical bonds never break,
restoring forces remain perfectly linear,
vibrational modes do not interact,
the Potential Energy Surface is exactly parabolic.

Real molecular vibrations deviate from these assumptions, particularly for

highly stretched bonds,
hydrogen-bonded systems,
floppy molecules,
low-frequency torsional motions,
highly excited vibrational states.

Consequently, calculated harmonic frequencies are often slightly higher than experimental values.

To compensate for this systematic error, empirical frequency scaling factors are commonly applied.

Harmonic vs Anharmonic Motion¶

Harmonic Approximation	Real (Anharmonic) Molecule
Parabolic potential	Asymmetric potential
Linear restoring force	Nonlinear restoring force
Constant force constant	Variable force constant
Independent vibrations	Coupled vibrations possible
Accurate near equilibrium	Accurate over the full PES

The harmonic model provides an excellent description of vibrations close to the equilibrium geometry but becomes less accurate for larger atomic displacements.

Role in Frequency Calculations¶

The Harmonic Approximation serves as the starting point for the mathematical treatment of molecular vibrations.

Optimized Geometry
        │
        ▼
Harmonic Approximation
        │
        ▼
Quadratic Potential
        │
        ▼
Hessian Matrix
        │
        ▼
Normal Modes
        │
        ▼
Vibrational Frequencies

Every subsequent step in Gaussian's frequency calculation is built upon this approximation.

Summary¶

The Harmonic Approximation assumes that atoms undergo only small oscillations about their equilibrium positions and that the local Potential Energy Surface can be represented by a quadratic (parabolic) function. Under this approximation, chemical bonds behave like ideal springs, allowing molecular vibrations to be described using the harmonic oscillator model. This simplification enables Gaussian to compute vibrational frequencies, infrared intensities, and thermodynamic properties efficiently while maintaining good accuracy for most equilibrium molecular structures.

Next Section¶

Having approximated the Potential Energy Surface as a quadratic function, the next chapter introduces the Hessian Matrix, which mathematically describes the curvature of this surface and provides the foundation for calculating molecular vibrational frequencies.