Curious problem using the variational method to find the ground state energy of the Harmonic oscillator.
Posted by peeterjoot on October 3, 2011
Motivation
Recap. Variational method to find the ground state energy.
Problem 3 of section 24.4 in the text [1] is an interesting one. It asks to use the variational method to find the ground state energy of a one dimensional harmonic oscillator Hamiltonian.
Somewhat unexpectedly, once I take derivatives equate to zero, I find that the variational parameter beta becomes imaginary?
I tried this twice on paper and pencil, both times getting the same thing. This seems like a noteworthy problem, and one worth reflecting on a bit.
Recap. The variational method.
Given any, not necessarily normalized wavefunction, with a series representation specified using the energy eigenvectors for the space
where
and
We can perform an energy expectation calculation with respect to this more general state
This allows us to form an estimate of the ground state energy for the system, by using any state vector formed from a superposition of energy eigenstates, by simply calculating
One of the examples in the text is to use this to find an approximation of the ground state energy for the Helium atom Hamiltonian
This calculation is performed using a trial function that was a solution of the interaction free Hamiltonian
This is despite the fact that this is not a solution to the interaction Hamiltonian. The end result ends up being pretty close to the measured value (although there is a pesky error in the book that appears to require a compensating error somewhere else).
Part of the variational technique used in that problem, is to allow Z to vary, and then once the normalized expectation is computed, set the derivative of that equal to zero to calculate the trial wavefunction as a parameter of Z that has the lowest energy eigenstate for a function of that form. We find considering the Harmonic oscillator that this final variation does not necessarily produce meaningful results.
The Harmonic oscillator variational problem.
The problem asks for the use of the trial wavefunction
to perform the variational calculation above for the Harmonic oscillator Hamiltonian, which has the one dimensional position space representation
We can find the normalization easily
After a bit of calculation, using integration by parts, we find for the energy expectation
Note that evaluating this integral requires the origin to be avoided where the derivative of becomes undefined. So one has to evaluate this in the limit as (which is an easy way to do it anyways since the absolute value can be eliminated by doubling the integral). Because of the curious end result, I also verified my calculation using Mathematica.
So, our ground state energy estimation, parametrized by becomes
(NOTE: TYPO FIXED ABOVE … THERE WERE TWO MINUS SIGNS, BUT THIS ERROR DIDN’T MAKE IT INTO THE DERIVATIVE BELOW).
Observe that if we set the derivative of this equal to zero to find the “best” beta associated with this trial function
we find that the parameter beta that best minimizes this ground state energy function is complex with value
So, it appears that we can’t minimize 4.10 to find a best ground state energy estimate associated with the trial function 4.7. We do however, know the exact ground state energy for the Harmonic oscillator. Is is possible to show that for all we have
? This inequality would be expected if we can assume that the trial wavefunction has a Fourier series representation utilizing the actual energy eigenfunctions for the system.
Some reflection.
I think that it is notable that I don’t believe the trial wave function for this problem lies in the span of the Hilbert space that describes the solutions to the Harmonic oscillator. Another thing of possible interest is the trouble near the origin for this wave function, when operated on by . Observe that Mathematica also required some special hand holding to deal with the origin.
I had initially thought that part of the value of this variational method was that we can use it despite not even knowing what the exact solution is (and in the case of the Helium atom, I believe it was stated in class that an exact closed form solution is not even known). This makes me wonder what restrictions must be imposed on the trial solutions to get a meaningful answer from the variational calculation?
Having started with a wavefunction that is probably not representable in the solution space is likely the bigger problem here. We probably need to adjust the treatment to account for that. Suppose we have
where is unknown, and presumed not orthogonal to any of the energy eigenkets. We can still calculate the norm of the trial function
Similarly we can calculate the energy expectation for this unnormalized state and find
Our normalized energy expectation is therefore the considerably messier
With a requirement to include the perpendicular cross terms the norm doesn’t just cancel out, leaving us with a clean estimation of the ground state energy. In order to utilize this variational method, we implicitly have an assumption that the and terms in the denominator are sufficiently small that they can be neglected.
For this harmonic oscillator problem, it would be interesting to calculate that Fourier remainder explicitly.
References
[1] BR Desai. Quantum mechanics with basic field theory. Cambridge University Press, 2009.
Mike Bantegui said
Peeter
I ran through your calculations and I just wanted to point out a useful for trick integrating quantities such as absolute values within Mathematica.
As you know, an expression of the form Integrate[Abs[f[x]], {x, -Infinity, +Infinity}] usually needs to be split up into Integrate[f[-x], {x, -Infinity, 0} + Integrate[f[x], {x, 0, Infinity}].
To get around this, you can simply use Sqrt[f[x]^2] for the absolute value. This has all the properties of the absolute value and Mathematica can handle this much more cleanly than Abs[f[x]]. All the analytic derivatives of this slightly modified square root are well defined.
peeterjoot said
Excellent. I’ll give that a try.
Curious problem using the variational method to find the ground state energy of the Harmonic oscillator (cont.) « Peeter Joot's Blog. said
[…] Curious problem using the variational method to find the ground state energy of the Harmonic os… […]
Ken said
Hey Peeter
I think your solution is in-correct. The double derivative of the trial wave function has a delta function. So, when you work out the expectation value of the kinetic part of the hamiltonian, you get a positive answer, which is turn makes your beta positive.
So, the derivative of the wave function gives you a step function and the derivative of the step function is a delta function…
peeterjoot said
Yes, I eventually figured out something along those lines (after talking to Professor Sipe about the problem). I put the corrected solution in my compilation of notes from when I took that class:
Click to access phy456.pdf
See: section: 2.3.3.6 Correcting, treating the origin this way.
Ken said
I see! Thanks!
Tom said
I know this is way after the fact but following your link to section 2.3.3.6 was very helpful after I was trying to solve Sakurai 5.20. Heavisides…Who would have thought! Thank-you!
peeterjoot said
Glad that this was useful. fyi, the current version of those notes, which may have some updates, is linked from here: http://peeterjoot.com/writing/graduate-quantum-mechanics-phy1520h/ (this particular section is “2.3.2.6 Correcting, treating the origin this way” in phy456.V0.1.9.pdf ). Also note that I have notes for a grad quantum course available (which followed Sakurai).