Explore Formation of Bound State in a Delta-Function Potential Well

Wavefunction
Messages
99
Reaction score
4

Homework Statement



We learned in class that a particle exposed to a 1D delta-function potential well would
always have a single bound state. Let us now explore this question for the case where the
delta-function potential well is situated in the vicinity of the impenetrable potential wall:
Given the potential energy:

V(x)=\begin{cases}<br /> +\infty,\quad\text{for}\quad x&lt;0\\<br /> -C\delta(x-a),\quad\text{for}\quad 0&lt;x&lt;\infty<br /> \end{cases}

Explore the formation of a bound state for a particle of mass m by solving the time independent
Schrödinger equation. Here a &gt; 0 specifies the position of the potential well
and C &gt; 0 its strength.
(a) Find the (transcendental) equation that would allow you to determine the energy of the
bound state. Indicate how you would be solving this equation graphically.
(b) In the absence of the neighboring wall (think if this may correspond to the limit a\rightarrow\infty),
the bound state is formed for any value of strength C. But now, do we have a restriction
on the value of C for a bound state to exist? If so, what is the minimal possible value
of C in terms of other given parameters of the problem?
(c) Can you bring simple physical arguments to rationalize the results in item (b)?

Homework Equations



Eqn 1):\frac{d^2\psi(x)}{dx^2}-\frac{2m}{\hbar^2}[V(x)-E]\psi(x)=0

The Attempt at a Solution



Part a):

First since the well is centered at x=a I'll split the x-axis into two regions:

Region 1: 0&lt;x&lt;a

Region 2: x&gt;a

Next since I'll setup my finite well of strength C:

Let the well have a width of 2\varepsilon. Then it's depth is \frac{C}{2\varepsilon}. So C must have dimensions of (J)(m). Let C=V_0(2\varepsilon) Where the well has a maximum depth of -V_0. The well has its walls at x=a\pm\varepsilon. In order for a bound state to exist E&lt;0. So there are two cases to consider within the well:

Case I: A particle with energy E

\frac{d^2\psi(x)}{dx^2}+\frac{2m}{\hbar^2}[E]\psi(x)=0

Let \kappa^2=-\frac{2mE}{\hbar^2}

\frac{d^2\psi(x)}{dx^2}-\kappa^2\psi(x)=0

Which has solutions of C_1e^{\kappa x}andC_2e^{-\kappa x}

Case II: A particle with energy -V_0-E

\frac{d^2\psi(x)}{dx^2}-\frac{2m}{\hbar^2}[-V_0-E]\psi(x)=0

Let k^2=\frac{2m}{\hbar^2}[V_0+E]

\frac{d^2\psi(x)}{dx^2}+k^2\psi(x)=0

which has solutions of C_3\cos(kx)andC_4\sin(kx)

Now I can impose my boundary conditions:

In Region 1; \psi_{I}(0)=0
\psi_{I}(a-\varepsilon)=\psi_{II}(a-\varepsilon)
\psi&#039;_{I}(a-\varepsilon)=\psi&#039;_{II}(a-\varepsilon)

In Region 2; \lim_{x\rightarrow\infty}\psi_{I}(x)=0
\psi_{I}(a+\varepsilon)=\psi_{II}(a+\varepsilon)
\psi&#039;_{I}(a+\varepsilon)=\psi&#039;_{II}(a+\varepsilon)

Let a&#039;=a-\varepsilon and a&#039;&#039;=a+\varepsilon

Now to solve the TISE in Region 1:

\psi_{I}(x)=C_1e^{\kappa x}+C_2e^{-\kappa x}

\psi_{II}(x)=C_3\cos(kx)+C_4\sin(kx)

\psi_{I}(0)=0=C_1+C_2\Rightarrow C_2=-C_1

\psi_{I}(x) = C_1[e^{\kappa x}-e^{-\kappa x}]

Now I'll select an even function from \psi_{II}(x)(Letting C_4=0). Such that \psi_{II}(x)=C_3\cos(kx)

Using my B.Cs:

\psi_{I}(a&#039;)=\psi_{II}(a&#039;)

Eqn 2:C_1[e^{\kappa a&#039;}-e^{-\kappa a&#039;}]=C_3\cos(ka&#039;)

\psi&#039;_{I}(a&#039;)=\psi&#039;_{II}(a&#039;)

Eqn 3:C_1\kappa[e^{\kappa a&#039;}+e^{-\kappa a&#039;}] = -C_3k\sin(ka&#039;)

Dividing Eqn 3 by Eqn 2:

Eqn 4: \kappa\frac{e^{\kappa a&#039;}+e^{-\kappa a&#039;}}{e^{\kappa a&#039;}-e^{-\kappa a&#039;}}=-k\tan(ka&#039;)

or equivalently

Eqn 4:\kappa a&#039;\tanh(\kappa a&#039;)=-ka&#039;\tan(ka&#039;)

Now I'll select an odd function from \psi_{II}(x)(Letting C_3=0). Such that \psi_{II}(x)=C_4\sin(kx)

\psi_{I}(a&#039;)=\psi_{II}(a&#039;)

Eqn 5:C_1[e^{\kappa a&#039;}-e^{-\kappa a&#039;}]=C_4\sin(ka&#039;)

\psi&#039;_{I}(a&#039;)=\psi&#039;_{II}(a&#039;)

Eqn 6:C_1\kappa[e^{\kappa a&#039;}+e^{-\kappa a&#039;}] = C_4k\cos(ka&#039;)

Dividing Eqn 3 by Eqn 2:

Eqn 7: \kappa\frac{e^{\kappa a&#039;}+e^{-\kappa a&#039;}}{e^{\kappa a&#039;}-e^{-\kappa a&#039;}}=k\cot(ka&#039;)

or equivalently

Eqn 7:\kappa a&#039;\tanh(\kappa a&#039;)=ka&#039;\cot(ka&#039;)

Now to solve the TISE in Region 2:

\psi_{I}(x)=C_5e^{\kappa x}+C_6e^{-\kappa x}

\psi_{II}(x)=C_7\cos(kx)+C_8\sin(kx)

\lim_{x\rightarrow\infty}\psi_{I}(x)=0\Rightarrow C_5=0

\psi_{I}(x)=C_6e^{-\kappa x}

Again I'll start with the even solution \psi_{II}(x)=C_7\cos(kx)

\psi_{I}(a&#039;&#039;)=\psi_{II}(a&#039;&#039;)

Eqn 8:C_6e^{-\kappa a&#039;&#039;}=C_7\cos(ka&#039;&#039;)

\psi&#039;_{I}(a&#039;&#039;)=\psi&#039;_{II}(a&#039;&#039;)

Eqn 9: -C_6\kappa e^{-\kappa a&#039;&#039;}=-C_7k\sin(ka&#039;&#039;)

Dividing Eqn 9 by Eqn 8:

Eqn 10: \kappa a&#039;&#039;=ka&#039;&#039;\tan(ka&#039;&#039;)

Now for the odd portion of \psi_{II}(x).

\psi_{I}(a&#039;&#039;)=\psi_{II}(a&#039;&#039;)

Eqn 11:C_6e^{-\kappa a&#039;&#039;}=C_8\sin(ka&#039;&#039;)

\psi&#039;_{I}(a&#039;&#039;)=\psi&#039;_{II}(a&#039;&#039;)

Eqn 12: -C_6\kappa e^{-\kappa a&#039;&#039;}=C_8k\cos(ka&#039;&#039;)

Dividing Eqn 12 by Eqn 11:

Eqn 13: -\kappa a&#039;&#039;=ka&#039;&#039;\cot(ka&#039;&#039;)

I also have that:

Eqn 14: \kappa^2+k^2=\frac{2mV_0}{\hbar^2}

\kappa=\sqrt{\frac{2mV_0}{\hbar^2}-k^2}

Eqn 15: \kappa a&#039;=\sqrt{\frac{2mV_0(a&#039;)^2}{\hbar^2}-(ka&#039;)^2}

Eqn 16: \kappa a&#039;&#039;=\sqrt{\frac{2mV_0(a&#039;&#039;)^2}{\hbar^2}-(ka&#039;&#039;)^2}

At this point I'm a little lost as to what I need to do to get a single transcendental equation; however, here is how I think I should proceed:

1) for all equations sub in the expressions for a&#039; and a&#039;&#039; in terms of a
2) use eqns 15 and 16 to get an expression for \kappa a
3) for regions 1 and 2 plot odd and even solutions in terms of \kappa a
(Edit: actually I just realized since \varepsilon is small I can take \varepsilon^2 to be 0 so what I'll end up with is just \kappa a = \sqrt{\frac{V_0}{\frac{\hbar^2}{2ma^2}}-(ka)^2}) I'm still not sure how to resolve the transcendental equations in regions 1 and 2 though. Thanks in advance!
 
Last edited:
Physics news on Phys.org
Instead of treating the delta distribution as a limit of a square potential, I suggest using it as it is and infer from this the boundary condition that must be fulfilled (in addition to continuity) by the wave function at x = a. You will then only have one matching to perform.
 
Orodruin said:
Instead of treating the delta distribution as a limit of a square potential, I suggest using it as it is and infer from this the boundary condition that must be fulfilled (in addition to continuity) by the wave function at x = a. You will then only have one matching to perform.

I'm not sure how to deal directly with the delta function potential. The way we dealt with such a potential in class was to construct such a finite square well. Though I'm curious to hear more about your approach. Thank you for the reply :)
 
The idea is to integrate the Schrodinger equation from ##a-\varepsilon## to ##a+\varepsilon## and take the limit as ##\varepsilon \to 0^+##. That'll give you the condition ##\psi'(x)## has to satisfy at the boundary.
 
  • Like
Likes Wavefunction
vela said:
The idea is to integrate the Schrodinger equation from ##a-\varepsilon## to ##a+\varepsilon## and take the limit as ##\varepsilon \to 0^+##. That'll give you the condition ##\psi'(x)## has to satisfy at the boundary.
Thank you I'll reformulate my problem and post my solution later :)
 
\lim_{\varepsilon\rightarrow0}\int_{x=a-\varepsilon}^{a+\varepsilon}\frac{-\hbar^2}{2m}\frac{d^2\psi}{dx^2}-C\delta(x-a)\psi(x)dx=\lim_{\varepsilon\rightarrow0}\int_{x=a-\varepsilon}^{a+\varepsilon}E\psi(x)dx

\lim_{\varepsilon\rightarrow0}[\psi&#039;(a+\varepsilon)-\psi&#039;(a-\varepsilon)]=\frac{2mC}{\hbar^2}\psi(a)

Now in region 1: \psi_{I}(x)=C_1e^{\kappa x}+C_2e^{-\kappa x}

applying the B.C in region 1 of \psi_{I}(0)=0\Rightarrow C_2=-C_1

Then \psi_{I}(x)=C_1(e^{\kappa x}-e^{-\kappa x})

Now for region 2: \psi_{I}(x)=C_3e^{-\kappa x}

which satisfies in region 2: \lim_{x\rightarrow\infty}\psi_{I}(x)=0

Now \psi_{I}(x) must be continuous at x=a so

C_1(e^{\kappa a}-e^{-\kappa a})=C_3e^{-\kappa a}

C_1=C_3\frac{e^{-\kappa a} }{e^{\kappa a}-e^{-\kappa a}}
 
Last edited:
In region 1: \psi&#039;_{I}(a-\varepsilon)=C_3\frac{e^{-\kappa a}}{e^{\kappa a}-e^{-\kappa a}}\kappa(e^{\kappa(a-\varepsilon)}+e^{-\kappa(a-\varepsilon)})

In region 2: \psi&#039;_{I}(a+\varepsilon)=-C_3\kappa e^{-\kappa a}

So \lim_{\varepsilon\rightarrow0}[\psi&#039;_{I}(a+\varepsilon)-\psi&#039;_{I}(a-\varepsilon)]=-\frac{2mC}{\hbar^2}C_3e^{-\kappa a}

Which becomes \lim_{\varepsilon\rightarrow0}[(-C_3\kappa e^{-\kappa (a+\epsilon)})-(C_3\frac{e^{-\kappa a}}{e^{\kappa a}-e^{-\kappa a}}\kappa(e^{\kappa(a-\varepsilon)}+e^{-\kappa(a-\varepsilon)}))]=-\frac{2mC}{\hbar^2}C_3e^{-\kappa a}

[(-C_3\kappa e^{-\kappa a})-(C_3\frac{e^{-\kappa a}}{e^{\kappa a}-e^{-\kappa a}}\kappa(e^{\kappa a}+e^{-\kappa a}))]=-\frac{2mC}{\hbar^2}C_3e^{-\kappa a}

Finally,

[(\kappa)+(\frac{1}{e^{\kappa a}-e^{-\kappa a}}\kappa(e^{\kappa a}+e^{-\kappa a}))]=\frac{2mC}{\hbar^2}

Which as a\rightarrow\infty gives me \kappa = \frac{mC}{\hbar^2}
 
Looks good. You could simplify the algebra a bit if you used
\begin{eqnarray*}
\psi_I(x) &= A\sinh ka \\
\psi_{II}(x) &= Be^{-ka}
\end{eqnarray*} In terms of hyperbolic trig functions, your transcendental equation is
$$k(1+\coth ka) = \frac{2mC}{\hbar^2}.$$
 
  • Like
Likes Wavefunction
Awesome thanks for help everyone! :) Cheers
 
Back
Top