Delta well + infinite barrier -> bound state

1. Jul 13, 2010

maverick280857

Hi,

I'm trying to understand the quantum mechanical solution to this potential:

$$V(x) = \left\{\begin{array}{cc}\infty & \mbox{ for } x < 0,\\-\lambda\delta(x-d) & \mbox { for } x > 0\end{array}\right.$$

A particle of mass m is constrained to move on the half straight line $\{x \in \mathbb{R}: x > 0\}$ under this potential. I want to examine the condition for existence of bound states, and its dependence on $\lambda$ and $d$, where both $\lambda$ and $d$ are positive.

For x < 0, $\psi(x) = 0$

For 0 < x < d, $\psi(x) = A\sin(kx)$ where $k = \sqrt{2mE}/\hbar$.

For x > d, $\psi(x) = Ce^{ikx} + De^{-ikx}$.

Enforcing continuity at x = d, $A sin(kd) = Ce^{ikd} + De^{-ikd}$.

Integrating Schrodinger's equation over a small interval around $x = d$

$$-\frac{\hbar^2}{2m}\int_{d-\epsilon}^{d+\epsilon}\frac{d^2 \psi}{dx^2} dx - \lambda\int_{d-\epsilon}^{d+\epsilon}\delta(x-d)\psi(x) dx = E\int_{d-\epsilon}^{d+\epsilon}\psi(x) dx$$

But I don't get a condition relating $\lambda$ and $d$ for the existence of a bound state. What have I missed?

Last edited: Jul 13, 2010
2. Jul 13, 2010

Jivesh

There might not be an explicit relation between \\lambda and d, if that what's you are trying to find. But you do have the relation that the difference of the derivatives to the left and right of the peak is equal to the negative of \\lambda times 2m, divided by hbar^2. Comparing the real and imaginary part, you might get a relation between the trigonometric ratios of d and \\lambda. I think the energy levels might have the discrete band structure, as is found for the case of a solid, with the potential represented by the dirac comb.

3. Jul 13, 2010

DrDu

For a bound state, the wavefunction has to decay exponentially for x>d.

4. Jul 13, 2010

maverick280857

Okay, I found the source of this problem and it asks me to find a relation between d and lambda for the existence of at least one bound state.

5. Jul 13, 2010

maverick280857

Right so that means k becomes imaginary.

6. Jul 13, 2010

DrDu

yes, and either C or D becomes 0.

7. Jul 13, 2010

maverick280857

Okay just to clarify matters...

For 0 < x < d, since the wavefunction is constrained to be 0 for x <= 0, it should have the form

$$\psi(x) = A\sin(kx)$$

where $k = \sqrt{\frac{2mE}{\hbar^2}}$ for E > 0. This makes k real.

For x > d, the wavefunction should have the form

$$\psi(x) = B e^{-kx}$$

for the same k. Now, what I want to understand is...what is the definition of a bound state here? Suppose the delta well didn't exist. Then, the "bound"
state would just be $\psi(x) = A\sin(kx)$ (or in other words, no bound state, as its just an infinite step at x = 0). With the delta well in place at x = d, what really is a bound state for the composite potential? In other words, how is an E > 0 solution a bound state?

8. Jul 13, 2010

DrDu

All solutions to the Schroedinger equation which are normalizable are bound states.
The other solutions, like sin(kx) in the free case are not normalizable and belong to the continuum. Whether one considers non-normalizable solutions as admissible solutions at all depends on the formalism. In the orthodox view a la von Neumann, they don't belong to the Hilbert space at all and are non-admissible. In the rigged Hilbert space approach of Bogoliubov which accomodates naturally the Dirac bra-ket formalism, they are ok.

9. Jul 13, 2010

maverick280857

Got it. So, are you saying that the solutions I listed -- the sin(kx) in (0, d) and the exponentially decaying wave for (d, \infty) are bound states for E > 0? The reason for my confusion in defining a "bound state" and distinguishing between E > 0 and E < 0 stems from the fact that we have a barrier which goes to +\infty and a well which goes to -\infty.

10. Jul 13, 2010

maverick280857

Here's my working.

$$\psi(x) = \left\{\begin{array}{cc}A\sin(kx) & 0 < x < d\\Ce^{-kx} & x > d\end{array}\right.$$

1. Continuity at x = d:

$$A\sin(kd) = Ce^{-kd}$$

2. Integration of Schrodinger's equation around x = d

$$\sin(kd) + \cos(kd) = \frac{2m\lambda}{\hbar^2 k}\sin(kd)$$

where $k = \sqrt{\frac{2mE}{\hbar^2}}$ and E > 0.

This gives (if my algebra is correct) the "condition relating $\lambda$ and $d$ for the existence of a bound state" as

$$\tan(kd) = \frac{1}{\frac{2m\lambda}{\hbar^2 k}-1}$$

Is this the desired condition, even though it involves another parameter (E)?

PS -- The question also says that an observer can detect the existence of the potential by bouncing back waves from the wall at x = 0, of the form sin(kx) which will acquire a phase shift when the delta well is present, and look like sin(kx + phi). If this is to be believed, the exponentially decaying solution should be replaced with the complex exponential wave, don't you think?

Last edited: Jul 13, 2010
11. Jul 13, 2010

DrDu

I didn't check your calculations but I'd rather expect that the bound state has E<0 in that problem as the continuum of free states still starts at E=0.

12. Jul 13, 2010

DrDu

Btw, instead of sin and cos you should also use sinh and cosh.

13. Jul 13, 2010

maverick280857

Why? Is that because E < 0? And why is that so?

PS -- The question also says that an observer can detect the existence of the potential by bouncing back waves from the wall at x = 0, of the form sin(kx) which will acquire a phase shift when the delta well is present, and look like sin(kx + phi). If this is to be believed, the exponentially decaying solution should be replaced with the complex exponential wave, don't you think?

14. Jul 13, 2010

DrDu

The ansatz in post #1 will yield the correct answer for scattering states.
For the bound state, as E<0, k will be imaginary and if you write exp(-kx) on the right, you should also write sinh (kx).

15. Jul 13, 2010

maverick280857

Ok, here is my revised working

$$k = \sqrt{-2mE/\hbar^2}$$

where E < 0

For 0 < x < d

$$\psi(x) = A\sinh(kx) + B\cosh(kx)$$

as $\psi(0) = 0$, B = 0.

for x > d,

$$\psi(x) = C\sinh(kx) + D\cosh(kx)$$

Integrating Schrodinger's equation,

$$k(C-A)\cosh(kd) + \left(kD + \frac{2m\lambda}{\hbar^2}A\right)\sinh(kd) = 0$$

which gives (please correct me if I'm wrong)

$$C = A$$
$$kD + \frac{2m\lambda}{\hbar^2}A = 0$$

So,

$$\psi(x) = \left\{\begin{array}{cc}A\sinh(kx) & \mbox{ for } 0 < x < d\\A\sinh(kx) -\frac{2m\lambda}{\hbar^2 k}A\cosh(kx) &\mbox{ for } x > d\end{array}\right.$$

But this doesn't give me any condition connecting $\lambda$ and $d$. Also, how can I write the wavefuntion as sin(kx + phi) when the delta barrier is present? The question suggests I should do so, to get to a transcendental equation for phi.

Last edited: Jul 13, 2010
16. Jul 13, 2010

DrDu

sigh, for x>d an independent combination of sinh and cosh is not allowed, as I thought I had explained before as these functions diverge when x goes to infinity. Only for D=-C, that is, for exp (-kx), a normalizable solution does result.

17. Jul 13, 2010

maverick280857

Oops, I'm sorry...the last step I wrote here is wrong...I figured out the problem btw, and was just going to post the correction. As usual, I had messed up because of 'k'. You're right of course, I wrote sinh() and cosh() but on my notebook I still had the old sin() and cos(). Everything works out now...including the condition on lambda and d, which is

$$\frac{2m\lambda d}{\hbar^2} > 1$$

(If this condition is satisfied, there is at least one bound state.)

As for the second part of the question, what they're actually proposing is that we bounce waves back from the wall at x = 0 (when lambda = 0). That is a scattering problem, and so has nothing to do with the bound states that we just computed. I was confused by the language because the question refers to the solution of the first part.

PS -- A bad day :-/

Last edited: Jul 13, 2010
18. Jul 15, 2010

maverick280857

DrDu, can you help me understand this second part properly. I am stuck somewhere.

An observer at $x = \infty$ can infer the existence of $V(x)$ by bouncing back waves from the wall at $x = 0$. For $\lambda = 0$, the wavefunction with energy $E = \hbar^2 k^2/2m$ would be $\psi(x) = \sin(kx)$. But with $\lambda \neq 0$, the wavefunction acquires a phase $\phi$ so that $\psi(x) = \sin(kx + \phi)$ at large x.

I have to find an equation determining $\phi$. In this case, the wavefunction is $\psi(x) = B\sin(kx + \phi)$ (for E > 0 and real k) for large x. Can this wavefunction be used as that near the delta well? And what is the wavefunction for $0 < x < d$ when the delta well is present? Can that just be $A\sin(kx)$?

19. Jul 15, 2010

DrDu

You can use that ansatz for x>d and for 0<x<d, A sin kx is ok.

20. Jul 15, 2010

maverick280857

Ok, so if I do that,

continuity at x = d gives

$$A\sin(kd + \phi) = B\sin(kd)$$

and integrating Schrodinger's equation around x = d gives

$$\frac{\hbar^2 k}{2m}\cos(kd + \phi) + \left(\lambda-\frac{\hbar^2 k}{2m}\cot(kd)\right)\sin(kd + \phi) = 0$$

This equation can be solved 'exactly' for $\tan(kd + \phi)$ and hence for $\phi$. But the questions asks for a 'transcendental equation' for $\phi$.