Unit 2: Quantum Chemistry-II: Simple Harmonic Oscillator and Rigid Rotator Models

1. The Simple Harmonic Oscillator (SHO) Model

The Simple Harmonic Oscillator (SHO) model is a fundamental quantum mechanical system used to describe the vibrational motion of molecules. Classically, a harmonic oscillator is a system governed by Hooke's Law, where the restoring force is directly proportional to the displacement from the equilibrium position.

Hooke's Law:
F = -kx

Where:

• F is the restoring force.
• k is the force constant (a measure of the stiffness of the chemical bond).
• x is the displacement from the equilibrium position (x = r - r_e).

The potential energy (V) associated with this restoring force is obtained by integrating the force over the displacement:

V(x) = (1/2)kx²

Setting up the Schrödinger Equation for SHO

To transition from classical mechanics to quantum mechanics, we substitute the potential energy term V(x) into the one-dimensional time-independent Schrödinger equation:

- (ħ² / 2μ) (d²ψ / dx²) + V(x)ψ = Eψ

For a diatomic molecule, we replace the mass with the reduced mass (μ) to describe the relative motion of the two nuclei about their center of mass:

μ = (m₁ × m₂) / (m₁ + m₂)

Substituting the potential energy term, we obtain the definitive Schrödinger equation for the quantum harmonic oscillator:

- (ħ² / 2μ) (d²ψ / dx²) + (1/2)kx²ψ = Eψ

Discussion of Wave Functions and Solutions

To solve this second-order differential equation, we introduce dimensionless variables. Let:

α = (kμ / ħ²)^1/2
y = α^1/2x

Using these substitutions, the wave equation is transformed and solved. The acceptable, physically well-behaved wave functions (which must be single-valued, continuous, and square-integrable so that ψ → 0 as x → ±∞) are expressed as:

ψ_v(x) = N_v × H_v(α^1/2x) × e^-αx2/2

Where:

• N_v is the normalization constant, given by:
  N_v = [ (α / π)^1/2 / (2^v v!) ]^1/2
• H_v(y) represents the Hermite Polynomials of degree v. These are mathematical polynomials that satisfy the boundary conditions.
• e^-αx2/2 is the Gaussian decay function, which ensures the wave function asymptotically approaches zero at infinite displacement.

The first three Hermite polynomials are:

• H₀(y) = 1
• H₁(y) = 2y
• H₂(y) = 4y² - 2

Key Features of the SHO Wave Functions

• Symmetry (Parity): The wave functions are either symmetric (even parity) or antisymmetric (odd parity) about the origin. For even quantum numbers (v = 0, 2, ...), the wave function is even; for odd quantum numbers (v = 1, 3, ...), the wave function is odd.
• Nodes: The wave function ψ_v(x) has exactly v nodes (points where the probability density is zero, excluding infinity).
• Classical Penetration (Tunneling): In quantum mechanics, there is a non-zero probability of finding the oscillator in the classically forbidden region (where the total energy E is less than the potential energy V), a phenomenon known as quantum tunneling.

2. Vibrational Energy of Diatomic Molecules and Zero-Point Energy

The solution to the SHO Schrödinger equation yields quantized energy levels. The allowed vibrational energy states for a diatomic molecule are given by:

E_v = (v + 1/2)hν

Where:

• v is the vibrational quantum number, which can take integer values: v = 0, 1, 2, 3, ...
• h is Planck's constant.
• ν is the classical vibrational frequency of the bond, defined as:
  ν = (1 / 2π) × √(k / μ)

Alternatively, in spectroscopic units (wavenumbers, cm^-1), the energy expression is written as:

G(v) = E_v / hc = (v + 1/2)ω_e

Where ω_e is the equilibrium vibrational frequency in cm^-1.

Significance of Zero-Point Energy

One of the most important consequences of quantum mechanics is that the minimum energy of the quantum harmonic oscillator is not zero. When the system is in its lowest possible energy level (v = 0), the energy is:

E₀ = (1/2)hν

This minimum, non-zero energy is called the Zero-Point Energy (ZPE).

Physical Significance of ZPE:

• Heisenberg Uncertainty Principle: If the ground-state energy were zero, the molecule would be completely at rest at the bottom of the potential well. This would mean its position (x = 0) and momentum (p = 0) would both be known with absolute precision, which directly violates the Uncertainty Principle (Δx × Δp ≥ ħ/2).
• Absolute Zero Behavior: Even at absolute zero (T = 0 K), the atoms in a diatomic molecule continue to vibrate with this characteristic residual energy.
• Chemical Bond Stability: ZPE must be factored into thermodynamic calculations when estimating bond dissociation energies. The actual energy required to dissociate a bond is less than the classical depth of the potential well (D_e) by an amount equal to the zero-point energy.

3. The Rigid Rotator Model

To analyze the rotational motion of a diatomic molecule, we use the Rigid Rotator model. In this idealized model, the molecule is treated as two point masses (m₁ and m₂) separated by a fixed, unchanging distance (r₀), rotating about their common center of mass.

Moment of Inertia

The rotation of the system is characterized by its moment of inertia (I), which acts as the rotational equivalent of mass:

I = μr₀²

Where μ is the reduced mass of the diatomic system.

Setting up the Schrödinger Equation

Since the distance between the nuclei is fixed, there is no potential energy associated with rotation (V = 0). The Hamiltonian operator consists solely of the kinetic energy operator. In spherical polar coordinates (r, θ, φ), with r fixed at r₀, the Schrödinger equation is expressed as:

- (ħ² / 2I) [ (1 / sinθ) (∂ / ∂θ)(sinθ ∂ψ/∂θ) + (1 / sin²θ) (∂²ψ / ∂φ²) ] = Eψ

Here, the wave function ψ depends only on the angular coordinates: ψ = Y(θ, φ).

Discussion of Application and Wave Solutions

The solution to this angular differential equation yields a set of functions known as Spherical Harmonics, denoted as Y_J^M(θ, φ). These functions are acceptable only when two quantum numbers are introduced to satisfy the cyclic boundary conditions:

1. J: The rotational quantum number (J = 0, 1, 2, ...).
2. M_J: The magnetic/rotational projection quantum number, which can take any of the 2J + 1 integer values: -J, -(J-1), ..., 0, ..., (J-1), J.

Allowed Rotational Energy Levels

Solving the equation yields quantized rotational energy levels:

E_J = J(J + 1) × (ħ² / 2I) = J(J + 1) × [ h² / (8π²I) ]

In rotational spectroscopy, these levels are typically expressed in terms of the rotational constant (B):

E_J = BJ(J + 1)hc

Where:

B = h / (8π²Ic)

B is typically expressed in wavenumbers (cm^-1).

Key Features of the Rigid Rotator Model

• Degeneracy: Because the energy expression E_J depends only on J and is independent of M_J, each rotational energy level J is degenerate. The degeneracy (g_J) of a level J is given by:
  g_J = 2J + 1
• Energy Spacing: Unlike the harmonic oscillator, which has equally spaced energy levels, the rotational energy levels get progressively further apart as J increases. The spacing between consecutive levels J and J+1 is:
  ΔE = E_J+1 - E_J = 2B(J + 1)
• No Zero-Point Rotational Energy: For J = 0, the rotational energy E₀ = 0. This is physically permissible because a non-rotating molecule does not violate the uncertainty principle in a way that prevents J = 0 from being a valid state (as angular orientation is unconstrained without rotation).

4. Summary and Comparison of Models

The following table summarizes the fundamental differences between the quantum mechanical models of vibration (SHO) and rotation (Rigid Rotator) for diatomic molecules.