Variational Approximation for the Harmonic Oscillator II (Worksheet)
- Page ID
- 40085
\( \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}\)Name: ______________________________
Section: _____________________________
Student ID#:__________________________
Template:HideTOCWork 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.
Constructing the Variational Energy
The variational method is one way of finding approximations to the lowest energy eigenstate or ground state. The method consists of constructing a "trial wavefunction" depending on one or more parameters (e.g., \(|\phi(\alpha, \beta, \gamma, \ldots) \rangle\) and then evaluating the "trial energy" (variational energy)
\[E_\phi (\alpha, \beta, \gamma, \ldots)=\dfrac{\langle\phi (\alpha, \beta, \gamma, \ldots) |\hat H|\phi (\alpha, \beta, \gamma, \ldots) \rangle}{\langle\phi (\alpha, \beta, \gamma, \ldots) |\phi (\alpha, \beta, \gamma, \ldots) \rangle} \label{W1}\]
where the Hamiltonian for the system. The wavefunction obtained by fixing the parameters to such values is then an approximation to the true wavefunction.
Q1
What are the limits on the number of parameters that a trial wavefunction (\(|\phi \rangle\) can have?
Q2
Identify which is hardest and why:
- Constructing the Hamiltonian for the system,
- Constructing the trial wavefunction for the system, or
- Evaluating the variational energy.
Q3
What is the origin of the denominator in Equation \ref{W1}? Does it always have to be in the equation? If not, when can you ignore it?
Variational Theorem
The variational theorem argues that this trial energy, \(E_\phi\) associated with the trial wavefunction for the known Hamiltonian is always greater than or equal to the true energy (\(E_\psi\)). Proof is not given.
\[ E_\phi (\alpha, \beta, \gamma, \ldots) \ge E_\psi \label{VM}\]
The variation method approximates the lowest energy eigenvalue, \(E_\psi\), and eigenfunction, \(\psi\), for a quantum mechanical system by guessing a function that is well-behaved over the limits of the system and minimizing the energy.
Q4
Under what condition(s) will this equation be true?
\[ E_\phi (\alpha, \beta, \gamma, \ldots) = E_\psi \label{VM1}\]
Q5
Explain the power (utility) of Equation \ref{VM}. If this equation were not true, would we be able to approximate the true solutions to the Hamiltonian of the unsolveable system using the variational method?
Minimizing the Variational Energy
The variational theory argues that when the energy is minimized, then
\[E_\phi (\alpha, \beta, \gamma, \ldots) \approx E_{actual}\]
and
\[| \phi (\alpha, \beta, \gamma, \ldots) \rangle \approx | \psi \rangle\]
The better the trial wavefunction resembles the true wavefunction, the more accurate these approximations are.
Q6
Since one does not know the true eigenstate \(| \psi \rangle\), how would one conclude that the optimized trial wavefunction is a good approximation to the true wavefunction?
Continued Application to the Harmonic Oscillator Hamilitonian
From the previous groupwork, we explore the trial wavefunction
\[| \phi(x) (\beta) \rangle =\dfrac{1}{1+\beta x^2} \label{trial}\]
where \(|\phi \rangle \) is the wavefunction that we guess and \(\hat H\) is the Hamiltonian for the system. The variational energy (from solving Equation \ref{W1}) was shown to be
\[E_\phi (\beta) =\dfrac{\langle\phi|\hat H|\phi\rangle}{\langle\phi |\phi\rangle}=\dfrac{\hbar^2\beta}{4\mu}+\dfrac{k}{2\beta}\]
Some equations you may find useful for the following discussion:
\[\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^{-ax^2/2} \]
Q7
Is the trial wavefunction (Equation \ref{trial}) normalized? And does it matter? Why or why not?
Q8
What is \(\dfrac{dE_\phi (\beta)}{d\beta}=0\) and why do you need to evaluate this derivative?
Q9
What value for \(\beta\) fulfills the minimized \(E_\phi\)?
Q10
What is \(E_\phi(\beta)\) for the value for \(\beta\) that fulfills the minimization?
Q11
The variation method approximates the ground state energy for the system. What is the expression for the exact energy of the harmonic oscillator?
Q12
What is the value for the quantum number for the ground state of the harmonic oscillator?
Q13
What is the exact energy for the ground state of the harmonic oscillator?
Q14
Considering only energies, how well does the optimizing trial wavefunction (Equation \ref{trial}) approximate the lowest energy harmonic oscillator eigenstate?
Q15
How could you improve this approximation?
Q16
Given the knowledge you have of the true harmonic oscillator wavefunction, how well would a different trial wavefunction, \(|\phi(x) (\beta) \rangle =e^{-\beta x^2}\) approximate the solution for the lowest energy state of the harmonic oscillator?