
#1
Jan913, 03:19 AM

P: 991

Okay, I have solved the schrödinger equation numerically by making it dimensionless (though Im still confused about this proces). And then approximating it on a finite interval and solving the resulting eigenvalue equation. This allows me to solve for the wave function of different potentials.
I started with the harmonic oscillator but have no reached the Gaussian one: V = V_{0}exp(x^{2}) In one simulation I am asked to find the difference in energy between the ground state and the first excited state as a function of V_{0}. On the attached graph I have done this for V_{0}=1..2..3...10 Does it look right? I am then asked the following: Solve the problem analytically by taylorexpanding the potential. So I taylor expand around x=0 to second order and find: V(x) = V_{0} + V_{0}x^{2} Plugging this into my dimensionless Schrödinger equation I get: ½∂^{2}ψ/∂x^{2} + (V_{0} + V_{0}x^{2})ψ = Eψ I thought aha. The x^{2}term can just be put in the harmonic oscillator form if we pick k=2V_{0} and the V_{0} term will just shift the energy of the oscillator, not alter the difference between E1 and E0. But in thinking it over again there are some problems. With my dimensionless equation I just had V(x)=½x^{2} for the harmonic oscillator. Now I have 2V_{0} in front of that. How will this constant effect my energies? And is all this even the right procedure? 



#2
Jan913, 08:50 AM

Mentor
P: 10,767

The difference in energy levels scales with the square root of the prefactor in a harmonic potential  in the dimensionless shape, this is hidden in the parameter transformations. In a region where this harmonic approximation is good (probably large V_{0}), I would expect that the difference grows with the square root of V_{0}, and your graph roughly looks like that. For small V_{0}, you probably get additional effects from higher orders of the potential.




#3
Jan913, 01:19 PM

P: 991

okay but I am meant to show this analytically. How can I do that?




#4
Jan913, 05:05 PM

Mentor
P: 10,767

Gaussian potential
With the standard formulas for a harmonic oscillator. In your dimensionless version, you should have the conversion factors somewhere.
Or with the simple sqrtdependence if you like. 



#5
Jan1013, 02:14 AM

P: 991

well the conversion formula to dimensionless units is x' = x * √(mω/hbar)
So do I go back to formulate it all in terms of x? :S Im very confused sorry. 



#6
Jan1013, 05:14 AM

Mentor
P: 10,767

That is possible. Alternatively, you can extract the scale of V from that prefactor.




#7
Jan1013, 01:15 PM

P: 991

okay I did the problem and did indeed find that the difference went like √(2V0)  I am just curious  how is it you can see that the harmonic approximation is better for bigger v0?




#8
Jan1013, 03:49 PM

Mentor
P: 10,767

Intuition. Deeper wells of the same size tend to have more bound states, so the lowest states are "deeper" in the well and smaller in terms of their size in x.



Register to reply 
Related Discussions  
Time Evolution Of A 1D Gaussian Wave Packet Under The Gravitational Potential  Quantum Physics  4  
Potential due to anisotropic gaussian distrubution  General Physics  0  
scattering of a gaussian wave packet at a potential  Introductory Physics Homework  1  
Fourier transform of a gaussian, throw stuff away, close to a gaussian?.  Calculus  0  
How to prove the output of Linear Filtering a Gaussian Process is still Gaussian?  Set Theory, Logic, Probability, Statistics  1 