Integrating the time-indepdent Schrodinger equation

wrldt
Messages
10
Reaction score
0

Homework Statement



A particle of mass m is confined to move in a one-dimensional and Diract delta-function attractive potential V(x)=-\frac{\hbar^2}{m}\alpha\delta(x) where \alpha is positive.

Integrate eh time-independent Schrodinger equation between -\epsilon and \epsilon. Let \epsilon\to 0, show that the derivative of the eigenfunction \phi(x) is discontinuous at x=0, and determine its jump at the point.

Homework Equations


\frac{\hbar^2}{2m}\frac{\partial^2 \psi}{\partial x^2}-V(x)\psi=E\psi

The Attempt at a Solution



First I want to determine if I am approaching the problem correctly.
First we have:
\int\limits_{-\epsilon}^{\epsilon}-\frac{\hbar^2}{2m}\frac{\partial^2 \psi}{\partial x^2}+\int\limits_{-\epsilon}^{\epsilon}\frac{\hbar^2}{m}\alpha\delta(x)\psi=\int\limits_{-\epsilon}^{\epsilon}E\psi.
Then this simplifies to
\int\limits_{-\epsilon}^{\epsilon}\frac{\hbar^2}{m}\alpha\delta\psi=\int\limits_{-\epsilon}^{\epsilon}\frac{\hbar^2}{2m}\frac{\partial^2 \psi}{\partial x^2}
Further simplifying
2\alpha\int\limits_{-\epsilon}^{\epsilon}\delta(x)\psi(x)=\frac{\partial \psi}{\partial x}\bigg|_{-\epsilon}^{\epsilon}

The last equation shows that as $\epsilon \to 0$, the right hand side will tend to zero but the left hand side will explode. I'm wondering if my work sufficiently shows this.

The part I am stuck on is showing what its jump is at x=0.
 
Physics news on Phys.org
wrldt said:
2\alpha\int\limits_{-\epsilon}^{\epsilon}\delta(x)\psi(x)=\frac{\partial \psi}{\partial x}\bigg|_{-\epsilon}^{\epsilon}

Your work up to here looks very good.

The last equation shows that as ##\epsilon \to 0##, the right hand side will tend to zero but the left hand side will explode. I'm wondering if my work sufficiently shows this.

Here's where you need to think some more about what you're getting. The left hand side does not explode. What does integrating over the delta function give you for any finite ##\epsilon##? Then let ##\epsilon## → 0. Similarly for the right hand side.
 
Last edited:
2\alpha\int\limits_{-\epsilon}^{\epsilon}\delta(x)\psi(x)=\frac{\partial \psi}{\partial x}\bigg|_{-\epsilon}^{\epsilon}

This simplies to something like 2\alpha\psi(0)=\Delta\left(\frac{\partial \psi}{\partial x}\right).

After some work I believe this simplies down to k=\alpha

So now our E=-\frac{\hbar^2\alpha^2}{2m}.

I'm basically following the process that is going on on p. 71-72 of Introduction to Quantum Mechanics.
 
Last edited:
To solve this, I first used the units to work out that a= m* a/m, i.e. t=z/λ. This would allow you to determine the time duration within an interval section by section and then add this to the previous ones to obtain the age of the respective layer. However, this would require a constant thickness per year for each interval. However, since this is most likely not the case, my next consideration was that the age must be the integral of a 1/λ(z) function, which I cannot model.
Back
Top