Skip to main content
Chemistry LibreTexts

5.3: The Harmonic Oscillator Approximates Vibrations

  • Page ID
    210814
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    The quantum harmonic oscillator is the quantum analog of the classical harmonic oscillator and is one of the most important model systems in quantum mechanics. This is due in partially to the fact that an arbitrary potential curve \(V(x)\) can usually be approximated as a harmonic potential at the vicinity of a stable equilibrium point. Furthermore, it is one of the few quantum-mechanical systems for which an exact, analytical solution exists. Solving other potentials typically require either approximations or numerical approaches to identify the corresponding eigenstates and eigenvalues (i.e., wavefunctions and energies).

    Anharmonic_oscillator.gif
    Figure \(\PageIndex{1}\): The \(\ce{HCl}\) molecule is really an anharmonic oscillator. Here \(HCl\) is vibrating at energy level \(E_3\). \(D_0\) is dissociation energy here, \(r_0\) bond length, \(V\) potential energy. Energy is expressed in wavenumbers. The hydrogen chloride molecule is attached to the coordinate system to show bond length changes on the curve. from Wikipedia.

    A general potential energy (\(V(x)\)) curve for a molecular vibration can be expanded as a Taylor series (Figure \(\PageIndex{2}\))

    \[V(x) = V(x_0) + \left. \dfrac {d V(x)}{d x} \right|_{x_0}^{x} (x - x_0) + \left. \dfrac {1}{2!} \dfrac {d^2 V(x)}{d x^2} \right|_{x_0}^{x} (x - x_0)^2 + \ldots + \left. \dfrac {1}{n!} \dfrac {d^n V(x)}{d x^n} \right|_{x_0}^{x} (x - x_0)^n \label{5.3.1}\]

    \(V(x)\) is often (but not always) shortened to the cubic term and can be rewritten as

    \[V(x) = \dfrac {1}{2} kx^2 + \dfrac {1}{6} \gamma x^3 \label{5.3.2}\]

    where \(V(x_0) = 0\), \(k\) is the harmonic force constant (harmonic term), and \(\gamma\) is the first (i.e., cubic) anharmonic term. It is important to note that this approximation is only good for \(x\) near \(x_0\), and that \(x_0\) stands for the equilibrium bond distance.

    altalt

    Figure \(\PageIndex{2}\): (left) Deviation from simple harmonic potential approximation (red curve) of true/exact potential (black curve) with cubic term (green). (right) Expansion with positive or negative quartic terms. Images used with permission from Heather Kulik (MIT) via OpenCourseWare.

    Almost all diatomics have experimentally determined \(\dfrac {d^2 V}{d x^2}\) for their lowest energy states. H2, Li2, O2, N2, and F2 have had terms up to \(n < 10\) determined of Equation \(\ref{5.3.1}\).

    Adding anharmonic perturbations to the harmonic oscillator (Equation \(\ref{5.3.2}\)) better describes molecular vibrations. Anharmonic oscillation is defined as the deviation of a system from harmonic oscillation, or an oscillator not oscillating in simple harmonic motion. Anharmonic oscillation is described as the restoring force is no longer proportional to the displacement. Figure \(\PageIndex{1}\) shows the the general potential with (numerically) calculated energy levels (\(E_0\), \(E_1\) etc.). \(D_o\) is the dissociation energy, which is different from the well depth \(D_e\). These vibrational energy levels of this plot can be calculated using the harmonic oscillator model (i.e., Equation \(\PageIndex{1}\) with the Schrödinger equation) and have the general form

    \[ E_v = \left(v + \dfrac{1}{2}\right) v_e - \left(v + \dfrac{1}{2}\right)^2 v_e x_e + \left(v + \dfrac{1}{2}\right)^3 v_e y_e + \text{higher terms} \label{5.3.7}\]

    where \( v \) is the vibrational quantum number and \( x_e\) and \( y_e\) are the first and second anharmonicity constants, respectively.

    The \(v = 0\) level is the vibrational ground state. Because this potential is less confining than a parabola used in the harmonic oscillator, the energy levels become less widely spaced at high excitation (Figure \(\PageIndex{1}\); top of potential).

    Limitations of the Harmonic Oscillator Model for Molecular Vibrations

    The harmonic oscillation is a great approximation of a molecular vibration, but has key limitations:

    • Due to equal spacing of energy, all transitions occur at the same frequency (i.e. single line spectrum). However experimentally many lines are often observed (called overtones).
    • The harmonic oscillator does not predict bond dissociation; you cannot break it no matter how much energy is introduced.
    Morse Potentials are better Approximations of Vibrational Motion

    A more powerful approach than just "patching up" the harmonic oscillator solution with anharmonic corrections is to adopt a different potential (\(V(x)\)). One such approach is the Morse potential, named after physicist Philip M. Morse, and a better approximation for the vibrational structure of the molecule than the harmonic oscillator because it explicitly includes the effects of bond breaking and accounts for the anharmonicity of real bonds (Figure \(\PageIndex{4}\)).

    alt
    Figure \(\PageIndex{4}\): The Morse potential (blue) and harmonic oscillator potential (green). Unlike the energy levels of the harmonic oscillator potential, which are evenly spaced by ħω, the Morse potential level spacing decreases as the energy approaches the dissociation energy. The dissociation energy \(D_e\) is larger than the true energy required for dissociation \(D_0\) due to the zero point energy of the lowest (\(v = 0\)) vibrational level. from Wikipedia.

    The Morse Potential is a good approximation to \(V(x)\) and is best when looking for a general formula for all \(x\) from 0 to \(\infty\), not just applicable for the local region around the \(x_o\):

    \[V(x) = D(1-e^{-\beta (x - x_0)})^2 \label{5.3.8}\]

    with \(V(x = x_0) = 0\) and \(V(x = \infty) = D\).

    The Morse Potential (Figure \(\PageIndex{4}\)) approaches zero at infinite \(r_e\) and equals \(-D_e\) at its minimum (i.e. \(r=r_e\)). It clearly shows that the Morse potential is the combination of a short-range repulsion term (small \(r\) values) and a long-range attractive term (large \(r\) values).

    Solving the Schrödinger Equation with the Morse Potential (Equation \ref{5.3.8}) is not trivial, but can be done analytically.

    \[ \hat{H}|\psi \rangle = E_n | \psi \rangle \]

    with

    \[ \begin{align} \hat{H} &= \hat{T} + \hat{V} \\[4pt] &= \dfrac{- \hbar ^2 d^2}{2m \;dx^2} + D(1-e^{-\beta (x - x_0)})^2 \label{5.3.9} \end{align} \]

    The solutions and energies for the Morse potential will not be used in this course and will not be discussed in more detail.

    Contributors

    • Peter Kelly (UCDavis)

    5.3: The Harmonic Oscillator Approximates Vibrations is shared under a not declared license and was authored, remixed, and/or curated by LibreTexts.