Time independent schrodinger equation: delta potential

[tim_lou]

I'm currently reading Griffith's book on time independent Schroedinger's equation about delta functions.

However, I complete dislike how the book deals with the delta distribution.

firstly, the book discusses how to solve:
-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2}-\alpha\delta\psi=E\psi

for E<0, by solving the equation piece-wise, and noticing that psi must be continuous:
\psi=\sqrt k e^{-k|x|}however, so far the alpha has not come into play, the book incorporate alpha by integrating the Schroedinger's equation from -epsilon to +epsilon, then take the limit of epsilon to zero.

However, in doing so, the first integration becomes (according to the book)
\left \lim_{\epsilon\rightarrow 0}\frac{d\psi}{dx} \right |_{-\epsilon}^{+\epsilon}that is complete bs to me. psi is not even twice differentiable in the usual sense. The whole idea is psi is solved to be a distribution, and the derivative is a distribution, how in the heck can the book just invoke the fundamental theorem of calculus?? at least if they do that, they should at least realize that the integral is improper (zero is a bigggg discontinuity) and make the integral:
\left \lim_{\epsilon\rightarrow 0}\frac{d\psi}{dx} \right |_{-\epsilon}^{0}+\left \frac{d\psi}{dx} \right |_{0}^{+\epsilon}<br />If anybody can provide additional insights and a more rigorous treatment of this thing... I'll be greatly appreciated. I want to understand psi in the distributional sense.

 

[StatMechGuy]

There's a discontinuity in the derivative of the function from the left and from the right, and the idea is to figure out that jump by looking at the delta function.

Also, don't worry about the wave function being a distribution. It's still a function that satisfies a certain differential equation.

 

[tim_lou]

ok, I get a better intuition out of it now... yes, the idea is to figure out the jumping of the discontinuity.

I guess I'm just a little bit uncomfortable treating psi by seeing it as a function after my math professor has taught us the "right" way in terms of distributions...

 

[StatMechGuy]

Never trust a math professor to teach you physics. They'll obscure a painfully obvious concept under layers of formalism just so that they can show that it does in fact work the way physicists know it works. The trick is to know when that level is needed, and when it isn't.

 

[lalbatros]

tim_lou,

Usually, when I deal with such things, I just replace the delta function by a very peaked Gaussian.
Then I try to get the results.
For elementary questions, this avoid me usually any trouble.
It the example here, everything becomes obvious by keeping that in mind.
Integrating the SE once around x=0 over a -sufficiently large compared to the Gaussian- is easy.

 

[jostpuur]

I'm not sure about rigor, but at least getting the same result in several different ways gives some confidence to believe that the calculations have some truth behind them. I know in total three different ways to deal with the delta potential in this problem. The epsilon trickery you have encountered is one.

Another way is to start with a square well potential, that has some fixed width R. The solution where psi is cosine in the well (without zeros in the well), and exponential outside, is the one that survives in the limit R->0. The boundary condition leads to some impossible equations. I don't remember precisly, but the are something like \cos(Ax)=Ax+B or similar, but this problem can be avoided but using parabola approximation to the cosine. Then the boundary conditions can be solved precisly. The parabola approximation leads to small error, but this error vanishes on the limit R->0, and you should be able to solve the wave function and it's energy level for the delta potential.

However, no doubt the easiest way is to just take the second derivative of the expression e^{-k|x|} by using identities \partial |x|=2\theta(x)-1 and \partial\theta(x)=\delta(x) blindly.

 

[George Jones]

Here's an attempt at a distributional treatment. The ideas were easier than the write-up seems to indicate. Informal reasoning usuall works well, but sometimes it leads ones astray. See https://www.physicsforums.com/showthread.php?p=998502". Here's Roger Penrose on the question of mathematical rigour in quantum theory: "Quantum mechanics is full of irritating issues of this kind. As the state of the art stands, one either be decidedly sloppy about such mathematical niceties and even pretend that position states and momentum sates are actually states, or else spend the whole time insisting on getting the mathematics right, in which case there is a contrasting danger of getting trapped in 'rigour mortis'. ... I am not at all sure what the correct answer is for making progress in the subject!''

A distribution is a continuous linear mapping (functional) from \mathcal{T} to \mathbb{C}, where \mathcal{T} is the space of test (i.e., sufficiently nice) functions. Any locally integrable function g naturally defines a distribution G:

G\left[ f\right] =\int_{-\infty }^{\infty }g\left( x\right) f\left( x\right) dx

for all test functions f. Not all distributions arise in such a fashion.

Even though the Dirac delta function, defined by \delta \left[ f\right] =f\left( 0\right), is an example of a distribution that doesn't arise in the above manner, it is convenient (and useful!) notationally to pretend that it does, i.e., that \delta is a function such that

f\left( 0\right) =\int_{-\infty }^{\infty }\delta \left( x\right) f\left( x\right) dx.

If g is any function such that

\int_{-\infty }^{\infty }g\left( x\right) dx=1,

then the family of functions g_{\varepsilon }\left( x\right) :=g\left( x/\varepsilon \right) /\varepsilon defines a family of distributions (as above) G_{\varepsilon }, with G_{\varepsilon }\rightarrow \delta (in the distributional or weak sense) as \varepsilon \rightarrow 0, i.e.,

f\left( 0\right) =\lim_{\varepsilon \rightarrow 0}\int_{-\infty }^{\infty }g_{\varepsilon }\left( x\right) f\left( x\right) dx.

The definitions of differentiation of distributions and multiplication of distributions by somewhat nice functions are both motivated by functions considered as distributions. Let g be a function that has locally integrable derivative g^{\prime }, and use g^{\prime } to define the distribution G^{\prime } in the usual way:

<br /> \begin{equation*}<br /> \begin{split}<br /> G^{\prime }\left[ f\right] &amp;= \int_{-\infty }^{\infty }g^{\prime }\left( x\right) f\left( x\right) dx \\<br /> &amp;= \left[ g\left( x\right) f\left( x\right) \right] _{-\infty }^{\infty }-\int_{-\infty }^{\infty }g\left( x\right) f^{\prime }\left( x\right) dx.<br /> \end{split}<br /> \end{equation*}<br />

Test functions die at \pm \infty, so the first term in the last line is zero, giving G^{\prime }\left[ f\right] =-G\left[ f^{\prime }\right]. This definition works for all distributions (including the Dirac delta function), not just those that correspond to functions.

Suppose v\left( x\right) =u\left( x\right) g\left( x\right). Then,

V\left[ f\right] =\int_{-\infty }^{\infty }v\left( x\right) f\left( x\right) dx=\int_{-\infty }^{\infty }g\left( x\right) \left[ u\left( x\right) f\left( x\right) \right] dx=G\left[ uf\right].

This motivates the definition of a distribution \left( uG\right) \left[ f\right] :=G\left[ uf\right] for G an arbitrary distribution and u a function.

Now to the problem at hand. Consider the standard textbook form of the Schrodinger equation with delta function potential,

-\frac{\hbar ^{2}}{2m}\frac{d^{2}\psi }{dx^{2}}-\alpha \delta \left( x\right) \psi =E\psi .

The second term on the left is a distribution, so all terms in the equation need to be distributions, i.e., try to find a function \psi that has corresponding distribution \Psi (not to be confused with the time-dependent wavefunction), and that satisfies the distributional Schrodinger equation

-\frac{\hbar ^{2}}{2m}\Psi ^{\prime \prime }\left[ f\right] -\left( \alpha \psi \delta \right) \left[ f\right] =E\Psi \left[ f\right]

for all test functions f.

The term on the right is a distribution that corresponds to a function, and therefore, if the first term also corresponds to a function, then (by rearrangement and linearity of distributions), then so does the delta distribution. Since the delta distribution doesn't correspond to a function, the first term on the left can't be a distribution that corresponds to a function. In other words, \psi is twice differentiable as a distribution, but not as a function. Clearly, the only place is a problem is at x=0, so provisionally, assume \psi is continuous everywhere and piecewise smooth on both sides of zero.

The two different notions of differentiation give two options: 1) turn \psi into the distribution \Psi and use distributional differentiation (defined above) to produce the distributions \Psi ^{\prime } and \Psi ^{\prime \prime }; 2) apply piecewise differentiation of functions to \psi (except at x=0) to produce functions \phi ^{\prime } and \phi ^{\prime \prime }, and use these function to define distributions \Phi ^{\prime } and \Phi ^{\prime \prime }. In order to distinguish between these two options, I've introduced somewhat awkward notation.

Interesting question: Does \Psi ^{\prime }=\Phi ^{\prime } and \Psi ^{\prime \prime }=\Phi ^{\prime \prime } ?

<br /> \begin{equation*}<br /> \begin{split}<br /> \Psi ^{\prime }\left[ f\right] &amp;= -\Psi \left[ f^{\prime }\right] \\<br /> &amp;= -\int_{-\infty }^{\infty }\psi \left( x\right) f^{\prime }\left( x\right) dx\\<br /> &amp;= -\left( \int_{-\infty }^{0^{-}}\psi \left( x\right) f^{\prime }\left( x\right) dx+\int_{0^{+}}^{\infty }\psi \left( x\right) f^{\prime }\left( x\right) dx\right)\\<br /> &amp;= -\left( \left[ \psi \left( x\right) f\left( x\right) \right] _{-\infty }^{0^{-}}-\int_{-\infty }^{0^{-}}\psi ^{\prime }\left( x\right) f\left( x\right) dx+\left[ \psi \left( x\right) f\left( x\right) \right] _{0^{+}}^{\infty }-\int_{0^{+}}^{\infty }\psi ^{\prime }\left( x\right) f\left( x\right) dx\right)\\<br /> &amp;= -\psi \left( 0^{-}\right) f\left( 0^{-}\right) +\psi \left( 0^{+}\right) f\left( 0^{+}\right) +\int_{-\infty }^{0^{-}}\phi ^{\prime }\left( x\right) f\left( x\right) dx+\int_{0^{+}}^{\infty }\phi ^{\prime }\left( x\right) f\left( x\right) dx<br /> \end{split}<br /> \end{equation*}<br />

Since both \psi and f are continuous,

<br /> -\psi \left( 0^{-}\right) f\left( 0^{-}\right) +\psi \left( 0^{+}\right) f\left( 0^{+}\right) =0,<br />

and hence \Psi ^{\prime }\left[ f\right] =\Phi ^{\prime }\left[ f\right].

<br /> \begin{equation*}<br /> \begin{split}<br /> \Psi ^{\prime \prime }\left[ f\right] &amp;= -\Psi ^{\prime }\left[ f^{\prime }\right]\\<br /> &amp;= -\Phi ^{\prime }\left[ f^{\prime }\right]\\<br /> &amp;= -\left( \int_{-\infty }^{0^{-}}\phi ^{\prime }\left( x\right) f^{\prime }\left( x\right) dx+\int_{0^{+}}^{\infty }\phi ^{\prime }\left( x\right) f^{\prime }\left( x\right) dx\right)\\<br /> &amp;= -\psi ^{\prime }\left( 0^{-}\right) f\left( 0^{-}\right) +\psi ^{\prime }\left( 0^{+}\right) f\left( 0^{+}\right) +\int_{-\infty }^{0^{-}}\phi ^{\prime \prime }\left( x\right) f\left( x\right) dx+\int_{0^{+}}^{\infty }\phi ^{\prime \prime }\left( x\right) f\left( x\right) dx\\<br /> &amp;= \left( \psi ^{\prime }\left( 0^{+}\right) -\psi ^{\prime }\left( 0^{-}\right) \right) f\left( 0\right) +\int_{-\infty }^{\infty }\phi ^{\prime \prime }\left( x\right) f\left( x\right) dx<br /> \end{split}<br /> \end{equation*}<br />

Consequently,

<br /> \Psi ^{\prime \prime }\left[ f\right] =\left( \psi ^{\prime }\left( 0^{+}\right) -\psi ^{\prime }\left( 0^{-}\right) \right) \delta \left[ f\right] +\Phi ^{\prime \prime }\left[ f\right]<br />

Use this in the distributional Schrodinger equation:

<br /> \begin{equation*}<br /> \begin{split}<br /> -\frac{\hbar ^{2}}{2m}\left( \left( \psi ^{\prime }\left( 0^{+}\right) -\psi ^{\prime }\left( 0^{-}\right) \right) \delta \left[ f\right] +\Phi ^{\prime \prime }\left[ f\right] \right) -\left( \alpha \psi \delta \right) \left[ f\right] &amp;= E\Psi \left[ f\right]\\<br /> -\left( \frac{\hbar^{2}}{2m}\left( \psi ^{\prime }\left( 0^{+}\right) -\psi ^{\prime }\left( 0^{-}\right) \right) +\alpha \psi \left( 0\right) \right) \delta \left[ f\right] &amp;= \left( E\Psi +\frac{\hbar ^{2}}{2m}\Phi ^{\prime \prime }\right) \left[ f\right]<br /> \end{split}<br /> \end{equation*}<br />

The distribution on the left side of this equation is a distribution that can't be produced by a function, while the distribution on the right is produced by a function. This is only possible if both sides equal zero. Hence, the left side gives the jump discontinuity in \psi ^{\prime }, while the right gives the standard differential equation satisfied by \psi (since f is an arbitrary test function) for x\neq 0.