Skip to main content
Chemistry LibreTexts

Worksheet 8 Solutions

  • Page ID
    95490
  • \( \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}\)

    Solution

    Name: ______________________________

    Section: _____________________________

    Student ID#:__________________________

    Work in groups on these problems. You should try to answer the questions without referring to your textbook. If you get stuck, try asking another group for help.

    The variation method approximates the lowest energy eigenvalue, \(E_\phi\), and eigenfunction, \(\phi\), for a quantum mechanical system by guessing a function that is well-behaved over the limits of the system and minimizing the energy. For the trial wavefunction \(| \phi(x)\rangle =\dfrac{1}{1+\beta x^2}\), the energy is

    \[E_\phi=\dfrac{\langle\phi|\hat H|\phi\rangle}{\langle\phi |\phi\rangle}\]

    where \(|\phi \rangle \) is the wavefunction that we guess and \(\hat H\) is the Hamiltonian for the system. The variational theory argues that when the energy is minimized, then \(E_\phi\geq E_{actual}\). Some equations you may find useful for a harmonic oscillator:

    \[\hat H_{HO}=-\dfrac{\hbar^2}{2\mu}\dfrac{d^2}{dx^2}+\dfrac{kx^2}{2}\]

    \[E_{n,HO}=h\nu \left(n+\dfrac{1}{2} \right)=\hbar \omega \left(n+\dfrac{1}{2}\right)\]

    \[ \underset{\text{lowest energy eigenstate}}{| \psi(x) \rangle}=\left(\dfrac{a}{\pi}\right)^{1/4}e^{-\dfrac{ax^2}{2}} \]

    Q1

    From the previous groupwork, we got

    \[E_\phi=\dfrac{\langle\phi|\hat H|\phi\rangle}{\langle\phi |\phi\rangle}=\dfrac{\hbar^2\beta}{4\mu}+\dfrac{k}{2\beta}\]

    To get the best value for the energy, we must minimize \(E_\phi(\beta)\) with respect to \(\beta\).

    What is, \(\dfrac{dE_\phi}{d\beta}=0\).

    \[\dfrac{dE_\phi}{d\beta}= \dfrac{\hbar^2}{4\mu}-\dfrac{k}{2\beta^2}\]

    What value for \(\beta\) fulfills the minimized \(E_\phi\)?

    \[\beta_{min} = \sqrt{\dfrac{4\mu k}{2\hbar^2}} =\sqrt{\dfrac{2\mu k}{\hbar^2}}\]

    Using the value for \(\beta\) that fulfills the minimization, what is \(E_\phi(\beta)\)?

    \[E_\phi ^{min}=\dfrac{\hbar^2\beta_{min}}{4\mu}+\dfrac{k}{2\beta_{min}}=\dfrac{\hbar^2\sqrt{\dfrac{2\mu k}{\hbar^2}}}{4\mu}+\dfrac{k}{2\sqrt{\dfrac{2\mu k}{\hbar^2}}} = \dfrac{\hbar\sqrt{2\mu k}}{4\mu}+\dfrac{k \hbar }{2\sqrt{2\mu k}} = \dfrac{\hbar\sqrt{2}}{4} \sqrt {\dfrac{k}{\mu}}+ \dfrac{\hbar}{2\sqrt{2}} \sqrt {\dfrac{k}{\mu}} = \dfrac{\hbar}{2} \sqrt {\dfrac{k}{\mu}} ( \dfrac{\sqrt{2}}{2} + \dfrac{1}{\sqrt{2}}) = \dfrac{\hbar}{2} \sqrt {\dfrac{k}{\mu}} ( \dfrac{\sqrt{2}}{2} + \dfrac{\sqrt{2}}{2})= \dfrac{\hbar}{2} \sqrt {\dfrac{k}{\mu}} \sqrt{2}\]

    The variation method approximates the ground state energy for the system. What is the expression for the exact energy of the harmonic oscillator?

    \[\hbar \sqrt {\dfrac{k}{\mu}} (\frac{1}{2} +v)\]

    where v = 0,1,2,... is a quantum number for a vibrational state.

    Q2

    What is the quantum number for the ground state of the harmonic oscillator?

    \[v=0\]

    What is the exact energy for the ground state of the harmonic oscillator?

    \[\frac{1}{2} \hbar \sqrt {\dfrac{k}{\mu}}\]

    How well does the guess of \(| \phi(x) \rangle =\dfrac{1}{1+\beta x^2}\) approximate the lowest energy harmonic oscillator eigenstate?

    \[\frac {E_{trial} - E_{true}}{E_{true}} = \frac {\dfrac{\hbar}{2} \sqrt {\dfrac{k}{\mu}} \sqrt{2} - \dfrac{1}{2} \hbar \sqrt {\dfrac{k}{\mu}}}{\frac{1}{2} \hbar \sqrt {\dfrac{k}{\mu}}} = \frac {\dfrac{\hbar}{2} \sqrt {\dfrac{k}{\mu}} (\sqrt{2} - 1) }{\dfrac{\hbar}{2} \sqrt {\dfrac{k}{\mu}}} = 0.41 \]

    The estimate with 41% error is not very good. The error of 5-10% would be considered acceptable.

    Given the knowledge you have of the true harmonic oscillator wavefunction, how well would the new trial wavefunction, \(|\phi(x) \rangle =e^{-\beta x^2}\) approximate the solution for the lowest energy state of the harmonic oscillator?

    Very well, since the functional form is the same for trial wavefunction and true wavefunction of harmonic oscillator.


    This page titled Worksheet 8 Solutions is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Nancy Levinger.

    • Was this article helpful?