Easy variational principle question that I can't integrate

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
8 replies · 3K views
mrausum
Messages
44
Reaction score
0

Homework Statement



Use trial wavefunction exp(-bx^2) to get an upper limit for the groundstate energy of the 1-d harmonic oscillator

The Attempt at a Solution



This is always going to give an integral of x^2*exp(-x^2). How do you do it? :/
 
Physics news on Phys.org
It's a Gaussian integral that comes out to be

[tex]\int_{-\infty}^{\infty} x^2 e^{-b x^2} = \frac{\sqrt{\pi}}{2b^{3/2}}[/tex]
 
yeah, I know that from wolfram alpha too. But that can't be right, because then my upper limit for the ground energy would depend on 1/b^3/2, which would mean if b was very large it would violate the uncertainty principle.
 
Did you normalize your wave test function?
 
Feldoh said:
Did you normalize your wave test function?

Yeah that's the one, thanks.
 
Wait..that still doesn't work. The norm.const = b^.5/(pi/2)^.5, which doesn't quite cancel with the b^3/2 in the denominator from the integral?
 
Can you show me some of your work? When I do it I get the appropriate contributions from the kinetic and potential expectation values.
 
sure:

[tex]<H> = <(b/\pi).exp(-bx^2)|H|(b/\pi).exp(-bx^2)>[/tex]


[tex]<H> = b/\pi\int^{inf}_{-inf}exp(-2bx^2).\left(\frac{2h^2x}{m}+\frac{1}{2}mw^2x^2\right)dx[/tex]
 
Last edited:
The normalization constant is [tex]A = (\frac{2b}{\pi})^{1/4}[/tex]

From there the Hamiltonian of a 1D harmonic oscillator:

[tex]-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+\frac{1}{2}m \omega x^2[/tex]

I'm not sure if you all ready dealt with the second partial with respect to x do your expectation value of the Hamiltonian looks wrong so far.

Remember it's [tex]|A|^2\int_{-\infty}^{\infty} (e^{-b x^2})^* (-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}+\frac{1}{2}m \omega x^2) e^{-b x^2} dx[/tex]