Double delta well separated by a delta barrier

[Flenzil]

Homework Statement

Consider a one-dimensional time-independent Schrodinger equation for an electron in a double quantum well separated by an additional barrier. The potential is modeled by:

V (x) = -γδ(x - a) -γδ(x + a) + βδ(x)​

Find algebraic equations which determine the energies (or k-values) of electron bound states for γ > 0 and arbitrary real β (positive or negative). Describe the symmetry of their wave functions in terms of even and odd
solutions. How many bound states do you expect for this system?

Homework Equations

You may nd it useful to work in units ħ = 1 and m = ½, and to introduce k defined as E = -k², where E < 0 is the energy of a bound state, so that k is real.

The Attempt at a Solution

I've tried to solve this for arbitrary E, since I don't understand what happens at the barrier for E<0. So I split the problem into 4 parts: x<-a, -a<x<0, 0<x<a and x>a. This results in the wavefunctions:

ψ₁ = Ae^iE½x + Be^-iE½x
ψ₂ = Ce^iE½x + De^-iE½x
ψ₃ = Fe^iE½x + Ge^-iE½x
ψ₄ = He^iE½x​

Using the conditions of continuity and differential continuity results in the following conditions for the coefficients:

Ce^-iE½a + De^iE½a = Ae^-iE½a + Be^iE½a
​

-Ce^-iE½a + De^iE½a = -(1+iγ/E^½)Ae^-iE½a + (1+iγ/E^½)Be^iE½a

F+G=C+D
​

-F+G=(iβ/E^½-1)C+(iβ/E^½+1)D

He^iE½a= Fe^-iE½a + Ge^iE½a
​

-He^-iE½a = (1+iγ/E^½)Fe^-iE½a - (1+iγ/E^½)Ge^iE½a
​

And now I'm not sure what to do. I can start to solve for transmission and reflection coefficients, but as far as I am aware bound states don't have these? It feels like at this point I need to specifiy the sign for E, but for E<0 I don't know what happens at the barrier since barriers don't have bound states.

 

[stevendaryl]

To simplify your life a little, you need to realize that if E &lt; 0, then the corresponding wave function looks like:

A e^{-kx} + B e^{+kx}

Where k = \sqrt{-E} (in the units with \hbar = 1, m = \frac{1}{2}. There is no "i" in the exponent. So when x \rightarrow -\infty, you have to have A = 0.

The continuity and derivative conditions across a delta function potential have the same form, whether the energy is positive or negative.

 

[Flenzil]

To simplify your life a little, you need to realize that if E &lt; 0, then the corresponding wave function looks like:

A e^{-kx} + B e^{+kx}

I know that but I was trying to solve it generically first because I don't understand what happens to an exponential at a barrier. I know that for a traveling wave, it would exponentially decrease throughout the barrier until it got to the other side with a lower amplitude. But if it's already exponential, and also below 0, and also width zero, then I have no idea.

 

[stevendaryl]

Schrodinger's equation is:

- \psi&#039;&#039;(x) + V(x) \psi(x) = E \psi(x)

Suppose that you have a delta-function at x=a, so V(x) = g \delta(x-a). Then we can integrate both sides of this equation from x=a - \epsilon to x=a + \epsilon.

\int_{a-\epsilon}^{a+\epsilon} - \psi&#039;&#039;(x) dx = - (\psi&#039;(a+\epsilon) - \psi&#039;(a-\epsilon))

\int_{a-\epsilon}^{a+\epsilon} g \delta(x-a) \psi(x) dx = g \psi(a)

\int_{a-\epsilon}^{a+\epsilon} E \psi(x) dx \approx E \epsilon \psi(a)

So when you take \epsilon very small, the right side is approximately zero, so we just have:

- (\psi&#039;(a+\epsilon) - \psi&#039;(a-\epsilon)) + g \psi(a) = 0

So that's the derivative condition. You have to apply that for all three delta-functions (which have different values for a and g). So that gives you three equations. You also have to apply continuity at all three delta-functions. That gives you three more equations. So you have 6 equations and 6 unknowns (the coefficients for B, C, D, F, G, H). If you solve them, that'll end up giving you an equation for k.

 

[Flenzil]

So doing that, I get
C=D\frac{2k}{\gamma} e^{2ka} (\frac{2k}{\gamma}-1) from the first well
G=F\frac{2k}{\gamma} e^{2ka} (\frac{2k}{\gamma}-1) from the second well
And then I use these with C+D=F+G from the barrier to get:
D=F and C=G

So that I have:
Be^{kx} for x<-a
Ce^{kx}+De^{-kx} for -a<x<0
De^{kx}+Ce^{-kx} for 0<x<a
He^{-kx} for x>a

Am I on the right track here?

Edit: I can also get B=H by using the D=F relation with the alternate solution from the wells:
D=B\frac{\gamma}{2k} e^{-2ka} from the first well
F=H\frac{\gamma}{2k} e^{-2ka} from the second well

 

[stevendaryl]

So doing that, I get
C=D\frac{2k}{\gamma} e^{2ka} (\frac{2k}{\gamma}-1) from the first well

I think that should be
C=De^{-2ka} (\frac{2k}{\gamma}-1)

 

[Flenzil]

I think that should be
C=De^{-2ka} (\frac{2k}{\gamma}-1)

Yes. Yes it should.

If I use the differential continuity at the barrier I come up with a different relation between C and D:
D=C\frac{2+\beta/k}{2-\beta/k}

These seems at odds with the relation from the first well:
C=D(\frac{2k}{\gamma}-1)e^{2ka}

Mathmatically they both seem correct but they don't agree with each other. Is this maybe because they're only valid at their respective potentials or are they just wrong?

 

[stevendaryl]

Let's recap: You have the four regions:

1. B e^{+kx}
2. C e^{kx} + D e^{-kx}
3. F e^{kx} + G e^{-kx}
4. H e^{-kx}

From the boundary between regions 1&2, you get:
B e^{-ka} = C e^{-ka} + D e^{+ka}
k (C e^{-ka} - D e^{+ka} - B e^{-ka}) = -\gamma B e^{-ka}

From the boundary between regions 2&3, you get:
C + D = F + G
k (F - G - C + D) = \beta (F+G)

From the boundary between regions 3&4, you get:
F e^{ka} + G e^{-ka} = H e^{-ka}
k (F e^{ka} - G e^{-ka} + H e^{-ka}) = -\gamma H e^{-ka}

I think you made a mistake in solving these. However, here is a trick to simplify things: Try two guesses about the parity:

First guess: \psi(-x) = \psi(x) This implies: B = H, C=G, D=F
Second guess: \psi(-x) = -\psi(x) This implies: B = -H, C=-G, D=-F

 

[Flenzil]

Ok, going through the calculations again, I've arrived at 3 equations for k, 1 for each boundary:

For even symmetry

Boundary 1: k=-k\frac{De^{ka}-Ce^{-ka}}{De^{ka}+Ce^{-ka}} +\gamma
Boundary 2: k=\frac{\beta}{2} \frac{D+C}{D-C}
Boundary 3: k=-k\frac{De^{ka}-Ce^{-ka}}{De^{ka}+Ce^{-ka}} +\gamma

And for odd symmetry:

Boundary 1: k=-k\frac{De^{ka}-Ce^{-ka}}{De^{ka}+Ce^{-ka}} +\gamma
Boundary 2: k is kind of undefined here because F-G-C+D =0, I assume because the wavefunction is zero here?
Boundary 3: k=-k\frac{De^{ka}-Ce^{-ka}}{De^{ka}+Ce^{-ka}} -\gamma

I'm a bit confused about this energy at the barrier, it doesn't depend on distance or anything.

 

[Flenzil]

I understand that if this were a non-delta well, I could compare k=\sqrt{-E} at the V=0 portion to \kappa =\sqrt{E+V_0} inside the well by forming an equation of k as a function of \kappa and finding the intersect of the two functions. But here that isn't possible here because there is no "inside the well".

 

[stevendaryl]

Ok, going through the calculations again, I've arrived at 3 equations for k, 1 for each boundary:

For even symmetry

Boundary 1: k=-k\frac{De^{ka}-Ce^{-ka}}{De^{ka}+Ce^{-ka}} +\gamma
Boundary 2: k=\frac{\beta}{2} \frac{D+C}{D-C}
Boundary 3: k=-k\frac{De^{ka}-Ce^{-ka}}{De^{ka}+Ce^{-ka}} +\gamma

I haven't checked whether those are right, but to get an eigenvalue for k, all three (or the first two, since the third is a repeat) equations must be true. So you have to solve the first equation for C in terms of D (which you already did, earlier). Plug that into the second equation. D should then cancel out, giving an equation involving k alone.

 

[Flenzil]

I haven't checked whether those are right, but to get an eigenvalue for k, all three (or the first two, since the third is a repeat) equations must be true. So you have to solve the first equation for C in terms of D (which you already did, earlier). Plug that into the second equation. D should then cancel out, giving an equation involving k alone.

Well, I get C=De^{2ka}(\frac{2k}{\gamma}-1) as before
Which I put into k=\frac{\beta}{2}\frac{D+C}{D-C}
Which comes out with k=\frac{\beta}{2}\frac{1+e^{2ka}(\frac{2k}{\gamma}-1)}{1-e^{2ka}(\frac{2k}{\gamma}-1)}

I'm having a hard time solving this for k, almost every term has a k in it. I've been checking the maths and it all seems to work, I'll keep checking though.

 

[stevendaryl]

Well, I get C=De^{2ka}(\frac{2k}{\gamma}-1) as before
Which I put into k=\frac{\beta}{2}\frac{D+C}{D-C}
Which comes out with k=\frac{\beta}{2}\frac{1+e^{2ka}(\frac{2k}{\gamma}-1)}{1-e^{2ka}(\frac{2k}{\gamma}-1)}

I'm having a hard time solving this for k, almost every term has a k in it. I've been checking the maths and it all seems to work, I'll keep checking though.

That's as good as you can get. It's a transcendental equation involving k. To solve for k, you can either do it numerically, or you can try expanding k in a power series (in powers of \gamma, maybe).

 

[Flenzil]

That's as good as you can get. It's a transcendental equation involving k. To solve for k, you can either do it numerically, or you can try expanding k in a power series (in powers of \gamma, maybe).

The question wants algebraic solutions for E though. Although I'm not sure if this forbids a power series.

 

[stevendaryl]

The question wants algebraic solutions for E though. Although I'm not sure if this forbids a power series.

Well, I don't think it has a closed-form solution. If \beta=0, then I think the equation for k turns out to be:

k = \frac{\gamma}{2}(1 + e^{-2ka})

There is no closed-form solution to that equation. What you can do is to graph both sides as a function of k, and see for what value of k the two curves intersect.

 

[Flenzil]

Well, I don't think it has a closed-form solution. If \beta=0, then I think the equation for k turns out to be:

k = \frac{\gamma}{2}(1 + e^{-2ka})

There is no closed-form solution to that equation. What you can do is to graph both sides as a function of k, and see for what value of k the two curves intersect.

I've been trying to solve it graphically but then I need values for \beta,\gamma and a and the functions change drastically depending on the values of these.

I'm also unsure about this solution since it's currently:

Be^{kx} for x<-a
De^{2ka}(\frac{2k}{\gamma}-1)e^{kx}+De^{-kx} for -a<x<0
De^{kx} + De^{2ka}(\frac{2k}{\gamma}-1)e^{-kx} for 0<x<a
Be^{-kx} for x>a

So going through the barrier has no effect on the wavefunction and going through the second well somehow rids the wavefunction of the (\frac{2k}{\gamma}-1)e^{2ka} term. I'm just worried that the question calls for some kind of equation for E and I'm unable to provide one.

 

[Flenzil]

So just in case people find this and want an answer. Yes, this is right and by "algebraic equations" the question didn't mean analytic so a relation that can be solved graphically is fine. Also, if you find that the odd solution ignores the central delta function completely, then you're headed in the right direction. Quantum mechanics is fun like that :).