Quantum mechanics, potential well

In summary, the homework statement states that a particle of mass is found to be inside a uni-dimensional potential well of the form: V(x)=0 for x \leq -a and a\leq x and V(x)=-V_0 for -a <x<a. The equation that is solved to determine the contour conditions and conditions of continuity is Schrödinger's equation. The eigenfunctions have a definite parity and are real for the "linked states". The energy of the linked states only take discrete values. When E>0, the wave function of the particle is found to be a function of x that takes discrete values.
  • #1
fluidistic
Gold Member
3,953
265

Homework Statement


A particle of mass m is find to be inside a uni-dimensional potential well of the form: [itex]V(x)=0[/itex] for [itex]x \leq -a[/itex] and [itex]a\leq x[/itex] and [itex]V(x)=-V_0[/itex] for [itex]-a <x<a[/itex].
1)Write down the corresponding Schrödinger's equation.
2)Consider the case [itex]-V_0<E<0[/itex]. Determine the contour conditions and conditions of continuity that an eigenfunction must satisfy inside this potential well.
3)Show that the eigenfunctions have a definite parity and that they are real for the "linked states". I'm not really sure about the word "linked state", this is my own translation of the original problem.
4)Show that the energy of the linked states only take discrete values. (Maybe linked states means eigenstates.)
5)Consider now that E>0. Find and graph the wave function of the particle. Calculate the transmission coefficient when the particle moves along the potential well.

Homework Equations



[itex]-\frac{\hbar ^2 }{2m} \frac{d ^2 \Psi (x)}{d x^2}+[V(x)-E]\Psi (x)=0[/itex].

The Attempt at a Solution


Part 1) is solved if you take the equation above and replace V(0) by its values, for each region of the potential well.
I won't put all my work for 2), since it's extremely long.
I reached that in the first region [itex](x \leq -a)[/itex], [itex]\Psi _I (x)=De^{\alpha x}[e^{x(\alpha+\beta )}-e^{x(\alpha - \beta)}][/itex].
For the region of the well, region II, [itex]\Psi _{II} (x)=D(e^{-\beta x}-e^{\beta x})[/itex].
For the third region, when [itex]x \geq a[/itex], [itex]\Psi _{III}(x)=De^{-\alpha x} [e^{x(\alpha - \beta)}-e^{x(\alpha + \beta)}][/itex].
Where [itex]\alpha = \frac{\sqrt {-2mE}}{\hbar}[/itex] and [itex]\beta = \frac{\sqrt{2m(-V_0-E)}}{\hbar}[/itex].
I never met any complex exponential since despite its negative look, all the arguments I met under the square roots were positive.
For the contour condition, I assumed that Psi should not diverge when x tends to negative and positive infinity. From this, I could "drop" 2 constants, or more precisely, equal them to 0.
For the continuity, I assumed that the function Psi and its first derivative must be continuous over the whole 3 regions. This is how I could reduce the number of unknown constant to 1 (I started with A, B, C and D).
My question is now... how do I obtain D? Should I normalize something?
I don't think I'm done for part 2).

For the first part of part 3), it's easy. In previous exercise I showed that if V(x)=V(-x) then Psi is either odd or even. And this is the case in this exercise so by the previous exercise, the eigenfunctions have a definite parity. I'll have to show that they are real and take discrete values.
But first, I want to deal with D... but have no idea how to.
Any help will be appreciated! Thanks!
 
Physics news on Phys.org
  • #2
Yes, you determine D by normalizing the wave function.
 
  • #3
vela said:
Yes, you determine D by normalizing the wave function.

Ah ok thanks! So I must do [itex]\int _{-\infty} ^{\infty} |\Psi (x)|^2dx=1[/itex], right? Of course I'll have to split the integral according to the 3 regions...
 
  • #4
Yes, exactly.
 
  • #5
I'm depressed. I made at least an error somewhere. If you look at my Psi II and Psi III, you can see that they are equal and this shouldn't be so... I realized this when trying to normalize the Psi function.
 
  • #6
Yeah, I was wondering why V0 makes an appearance in your solutions for regions I and III.
 
  • #7
Ok let's retake this.
Do you buy that the solution to S. equation for the region I is of the form [itex]\Psi _I (x)=Ae^{\alpha x}+Be^{-\alpha x}[/itex] (but B is worth 0 from the contour conditions), thus [itex]\Psi ' _I (x)=\alpha A e^{\alpha x}[/itex]?
For the other 2 regions I get [itex]\Psi _{II}(x)= C e^{\beta x}+De^{- \beta x}[/itex], thus [itex]\Psi ' _{II}(x)= \beta C e^{\beta x}-\beta D e^{- \beta x}[/itex].
And [itex]\Psi _{III}(x)=Fe^{\alpha x }+Ge^{-\alpha x}[/itex] (G=0 for the later mentioned reason), thus [itex]\Psi ' _{III}(x)=-\alpha G e^{-\alpha x}[/itex].
Now I say that F and B are worth 0 as conditions of no divergence of Psi when x tends to infinities.
I get the continuity equations:
(1) [itex]Ae^{- \alpha a}=Ce^{-\beta a}+D e^{\beta a}[/itex]
(2) [itex]Ce^{\beta a}+D e^{-\beta a}=Ge^{-\alpha a}[/itex]
(3) [itex]\alpha A e^{-\alpha a}=\beta C e^{-\beta a}-\beta D e^{\beta a}[/itex]
(4) [itex]\beta C e^{\beta a}-\beta D e^{-\beta a}=G e^{-\alpha a}[/itex].
Is what I've done so far okay?
After this, what I've done of my lengthy draft is playing with algebra to reach that, assuming C=-D, G must equal -A. I also isolated A to be worth [itex]D[e^{a (\alpha + \beta)}-e^{a(\alpha - \beta)}][/itex].
So that I have A, G, and C in terms of D... the remaining unknown constant to be determined with the normalization of the Psi function.
 
  • #8
You should get sinusoidal solutions inside the well when E > -V0

Probably just a typo, but you're missing an alpha on the righthand side of (4). Your approach is fine.
 
Last edited:
  • #9
vela said:
You should get sinusoidal solutions inside the well when E > -V0
Hmm ok, good to know.
Probably just a typo, but you're missing an alpha on the righthand side of (4). Your approach is fine.
A typo that costs me a lot! Originally a typo on my draft. I wrote [itex]\Psi ' _{II}(a)= \Psi _{III}(a)[/itex] instead of [itex]\Psi ' _{II}(a)= \Psi' _{III}(a)[/itex]. I must redo the algebra.
By the way, is it ok to assume C=-D?
Thank you very much for your help!
 
  • #10
From the symmetry of the potential, you should expect to get even and odd solutions, which will correspond to C=D and C=-D or, respectively, cosines and sines. Those two conditions should come out of the math somehow.
 
  • #11
vela said:
From the symmetry of the potential, you should expect to get even and odd solutions, which will correspond to C=D and C=-D or, respectively, cosines and sines. Those two conditions should come out of the math somehow.

Ah I see. Well I'll try to reach this mathematically. I've tried to deal with matrices but after about 1 hour of algebra, I reach non sense so I made at least 1 error as usual. I even tried wolfram alpha (failed to calculate the coefficient A, C, D and G.)
Even mathematica 7:
Code:
Solve[{A*e^(-k*a) - C*e^(-b*a) - D*e^(b*a) == 0, 
  C*e^(b*a) + D*e^(-b*a) - G*e^(-a*k) == 0, 
  k*A*e^(-a*k) - C*b*e^(-b*a) + b*D*e^(b*a) == 0, 
  b*C*e^(b*a) - b*D*e^(-b*a) - k*G*e^(-k*a) == 0}, {A, C, D, G}]
returns
Code:
{{G -> 0, A -> 0, C -> 0, D -> 0}}
...
It's almost 2 am, time to sleep for me. I hope I'll wake up and attack this when I just wake up!
 
  • #12
Being lazy, I tried the same thing in Mathematica and got the same result. The problem is if you have four independent equations, the only solution is the trivial solution. You have to find the appropriate values for [itex]\alpha[/itex] and [itex]\beta[/itex] so the equations are no longer independent. These values happen to correspond to the allowed energies of the bound state.

I'm going to suggest you switch to sines and cosines for region II, just to simplify the algebra. Your four equations become:
\begin{align*}
C\cos k_2 a - D \sin k_2 a &= A e^{-k_1 a} \\
C \sin k_2 a + D \cos k_2 a &= \frac{k_1}{k_2} A e^{-k_1 a} \\
C\cos k_2 a + D \sin k_2 a &= G e^{-k_1 a} \\
C \sin k_2 a - D \cos k_2 a &= \frac{k_1}{k_2} G e^{-k_1 a}
\end{align*}where I replaced your [itex]\alpha[/itex] and [itex]\beta[/itex] by k1 and k2 because they're shorter to type.

If you add the first and third equations together, you get
[tex]2C\cos k_2 a = (A+G)e^{-k_1 a}[/tex]Similarly, if you add the second and fourth equations together, you get
[tex]2C\sin k_2 a = \frac{k_1}{k_2}(A+G)e^{-k_1 a}[/tex]If you divide these two equations, a bunch of factors cancel, and you get
[tex]\tan k_2 a = \frac{k_1}{k_2}[/tex]This is one of the conditions you are looking for. You should be able to show that when this condition is satisfied, D has to vanish. Then it follows that A=G, and you have your wave function up to the normalization constant. If you do the same thing except this time isolating the terms with D in it, you'll get another condition, which leads to the rest of the solutions.
 
  • #13
Ok thank you very much vela. I do have a few questions though, I don't really get what you get.
From your first equation it seems like [itex]e^{- \beta a}[/itex] transforms into [itex]\cos (\beta a)[/itex] and [itex]e^{\beta a}[/itex] transforms into [itex]-\sin (\beta a)[/itex]. However if I transform all equations according to this rule, I get different equations than yours, but for equation 1. So I didn't really understand how you simplified the exponential function for the region II.
 
  • #14
Your C and D aren't the same as my C and D. I should have made that clear. Remember that [itex]\beta=ik_2[/itex] will be complex, so [itex]e^{\beta a} = \cos k_2a + i\sin k_2a[/itex] and [itex]e^{-\beta a} = \cos k_2a - i\sin k_2a[/itex]. Then you can collect terms and rename the constants. Alternately, you can write the solution to region II down as [itex]\psi_\mathrm{II}(x) = C \cos k_2 x + D \sin k_2 x[/itex] and then apply the boundary conditions.
 
  • #15
vela said:
Your C and D aren't the same as my C and D. I should have made that clear. Remember that [itex]\beta=ik_2[/itex] will be complex, so [itex]e^{\beta a} = \cos k_2a + i\sin k_2a[/itex] and [itex]e^{-\beta a} = \cos k_2a - i\sin k_2a[/itex]. Then you can collect terms and rename the constants. Alternately, you can write the solution to region II down as [itex]\psi_\mathrm{II}(x) = C \cos k_2 x + D \sin k_2 x[/itex] and then apply the boundary conditions.

Ok, I'm going to retake all. I've tried to keep exponentials and reached [itex]C=D \frac{\frac{\beta }{\alpha }+1}{1-\frac{\beta}{\alpha}}[/itex] but I think it's wrong because [itex]C \neq \pm D[/itex] in this case I think.
So basically if I restart the whole exercise, I should define [itex]\beta = \frac{\sqrt{2m(-V_0-E)}}{\hbar}= \frac{i \sqrt {2m (V_0+E)}}{\hbar}[/itex] for convenience?
 
  • #16
I would actually pull the i out and define β as a real quantity, so your solution in region II would be
[tex]\psi(x) = C e^{i\beta x} + D e^{-i\beta x}[/tex]or
[tex]\psi(x) = C \cos \beta x + D \sin \beta x[/tex]I find it's simpler to work with real-valued variables if you can.
 
  • #17
vela said:
I would actually pull the i out and define β as a real quantity, so your solution in region II would be
[tex]\psi(x) = C e^{i\beta x} + D e^{-i\beta x}[/tex]or
[tex]\psi(x) = C \cos \beta x + D \sin \beta x[/tex]I find it's simpler to work with real-valued variables if you can.

Hmm I don't really understand. My beta was supposed to be real. It is/was worth [itex]\frac{\sqrt {2m(-V_0-E)}}{\hbar}[/itex] where both [itex]V_0[/itex] and [itex]E[/itex] are negative, so that [itex]-V_0-E>0[/itex]. Thus the square root is real and so is beta.
If I factorizes by "i", I get a negative argument in the square root, which is worth a complex number (I can't define the new "beta" to be real in this case, which is what you've done I think. I don't really understand this part). In all cases beta is real valued.
 
  • #18
The way you originally stated the problem, the potential is -V0 at the bottom of the well, which implied V0>0 so that -V0<0. Since E<0 and it needs to be greater than -V0, you have E > -V0 so 0 > -V0-E. (These signs always mess me up too.)
 
  • #19
vela said:
The way you originally stated the problem, the potential is -V0 at the bottom of the well, which implied V0>0 so that -V0<0. Since E<0 and it needs to be greater than -V0, you have E > -V0 so 0 > -V0-E. (These signs always mess me up too.)

Oh... I totally missed this! Well thank you very much for pointing this out. I'm going to redo everything with this in mind!
 
  • #20
vela said:
I would actually pull the i out and define β as a real quantity, so your solution in region II would be
[tex]\psi(x) = C e^{i\beta x} + D e^{-i\beta x}[/tex]or
[tex]\psi(x) = C \cos \beta x + D \sin \beta x[/tex]I find it's simpler to work with real-valued variables if you can.

Okay your beta is worth [itex]\frac{ \sqrt{2m(V_0+E)}}{\hbar}[/itex] which is indeed real.
But I don't understand how you get such values for Psi.
[itex]\Psi _{II}(x)=Ce^{i kx}+De^{-ikx}[/itex] where my k is worth your beta.
Now, this is worth [itex]C[\cos (kx)-i \sin (kx)]+D[\cos (kx)+i \sin (kx)][/itex]. But from here, I don't undestand how you simplified to get [itex]C \cos (kx)+D \sin (kx)[/itex]. Could you explain a bit more please?
Thanks for all your help and time.
 
  • #21
Again, the first C and D isn't the same as the second C and D. :) They're arbitrary constants, so I just keep reusing the names. I should stop doing that.

After you expand the exponentials, you collect the terms and you'll get
[tex]\psi_\mathrm{II} = (C+D)\cos kx + i(-C+D)\sin kx[/tex]Since C and D are arbitrary, we can replace the combinations with two new arbitrary constants, so I just renamed them as
\begin{align*}
C+D &\rightarrow C \\
i(-C+D) &\rightarrow D
\end{align*}
 
  • #22
Ah I see, thank you once again.
I think I still reach non sense. :frown:
The continuity equations are:
(1) [itex]Ae^{\alpha a}=C \cos (ka)-D \sin (ka)[/itex]
(2)[itex]-\alpha Ae^{\alpha a}=Ca \sin (ka)-Da \cos (ka)[/itex]
(3) [itex]Ge^{-\alpha a}=C \cos (ka) + D \sin (ka)[/itex]
(4) [itex]- \alpha G e^{-\alpha a}=-Ca \sin (ka) + Da \cos (ka)[/itex]
From (1) and (2) (divide (2) by - alpha and then equate both equations) I reach that [itex]C=D \frac{[\frac{a}{\alpha}\cos (ka)+sin (ka)]}{[\cos (ka)+ \frac{a}{\alpha}\sin (ka)]}[/itex].
Now [itex]C=D[/itex] + the condition that [itex]a\neq \alpha[/itex] implies that [itex]\tan (ka)=1[/itex] while [itex]C=-D[/itex] + the condition that [itex]a\neq \alpha[/itex] implies that [itex]\tan (ka)=-1[/itex]. I'm more than sure that I made at least 1 error since [itex]C= \pm D[/itex] didn't pop up as it should.
 
  • #23
There are a couple of mistakes in your equations. Go back and reread what I wrote in post #12.

If you use sines and cosines, you get either C=0 or D=0, not C=D or C=-D.
 
  • #24
vela said:
There are a couple of mistakes in your equations. Go back and reread what I wrote in post #12.

If you use sines and cosines, you get either C=0 or D=0, not C=D or C=-D.

Ah I see. I'm digging into this now. I reached your result. I'll try to show that D vanishes.
 
  • #25
I must say I'm getting lost in the algebra. I reached [itex]D=\frac{\alpha}{k} \left [ \frac{Ae^{-\alpha a}}{\cos (ka)}-C \right ]=C \frac{k}{\alpha}-\frac{A e^{-\alpha a}}{\sin (ka)}[/itex].
I isolated C in function of A from this expression. But I couldn't show that D was worth 0 yet.
 
  • #26
Use the fact that [itex]\alpha/k = \tan ka[/itex] (or cot ka, I don't remember which).
 
  • #27
vela said:
Use the fact that [itex]\alpha/k = \tan ka[/itex] (or cot ka, I don't remember which).

In fact I used it to reach the expression of D I gave in my last post. If I use them once again I go backward one step.
[itex]D= \tan (ka) \left [ \frac{Ae^{-\alpha a}}{\cos (ka)}-C \right ]=C \cot (ka) - \frac{Ae^{-\alpha a}}{\sin (ka)}[/itex].
 
  • #28
What I did was solve for Ae-αa in two of the equations and then set the results equal to each other:
\begin{align*}
C\cos ka - D \sin ka &= \frac{k}{\alpha}(C\sin ka+D\cos ka) \\
&= (\cot ka)(C\sin ka+D\cos ka)
\end{align*}Can you take it from there?
 
  • #29
vela said:
What I did was solve for Ae-αa in two of the equations and then set the results equal to each other:
\begin{align*}
C\cos ka - D \sin ka &= \frac{k}{\alpha}(C\sin ka+D\cos ka) \\
&= (\cot ka)(C\sin ka+D\cos ka)
\end{align*}Can you take it from there?

Oh this helped a lot, thanks for sharing your approach.
I reach [itex]D[\cot (ka) \cos (ka)+\sin (ka)]=0[/itex]. It's not very obvious to me that this imply that D=0. But I guess it does?
 
  • #30
Try multiplying through by sin ka.
 
  • #31
vela said:
Try multiplying through by sin ka.

Wow. I totally missed it. Thank you very much for all your help. Time for me to sleep. Hopefully I'll solve the problem before tomorrow bed-time.
Have a nice day/night.
 
  • #32
So this gives me:
[itex]\Psi _I (x)=C \cos (ka) e^{\alpha (a+x)}[/itex]
[itex]\Psi _{II}(x)=C \cos (kx)[/itex]
[itex]\Psi _{III}(x)=C \cos (ka)e^{\alpha (a-x)}[/itex].
Now I know I must normalize but... I think I made some error(s).
The question 3 says I must show that Psi is either odd or even. But the exponentials are a problem to me...
Psi is real valued as I must have showed though, so this work.
 
  • #33
That looks fine. Remember [itex]\alpha a[/itex] term in the exponent merely results in a multiplicative constant because [itex]e^{\alpha(a+x)} = e^{\alpha a} e^{\alpha x}[/itex].

If you wrote the eigenfunction in terms of A instead of C, it becomes even more plain to see:
\begin{align*}
\psi_\mathrm{I}(x) &= Ae^{\alpha x} \\
\psi_\mathrm{II}(x) &= A(e^{-\alpha a}\sec ka)\cos kx \\
\psi_\mathrm{III}(x) &= Ae^{-\alpha x}
\end{align*}
 
  • #34
Ah I see now, even with Psi written in terms of C.
I tried to normalize it. I reached that [itex]\int _{-\infty} ^{-a}|\Psi _I (x)|^2dx=\frac{C^2 \cos ^2 (ka)}{2\alpha}=\int _{a}^{\infty} |\Psi _{III}(x)|^2dx[/itex] (and that's a good sign that they are equal!) while [itex]\int _{-a}^a |\Psi _{II}|dx=C^2 \left [ \frac{\sin (2ak)}{2k}+a \right ][/itex].
So that [itex]C^2=\left [ \frac{\cos ^2 (ka)}{\alpha}+\frac{\sin(2ak)}{2k}+a \right ] ^{-1}[/itex]. Now I suppose I must take the positive square root since the probability must be positive.
Is that okay? This doesn't look very nice on my eyes though.
 
  • #35
Do you need to normalize the wave functions for this problem?
 
Back
Top