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

Click For Summary

Homework Help Overview

The discussion revolves around the formation of a bound state in a one-dimensional delta-function potential well, particularly when the well is near an impenetrable potential wall. Participants are exploring the implications of the potential energy defined by V(x) and the conditions for bound states based on the Schrödinger equation.

Discussion Character

  • Exploratory, Conceptual clarification, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • Participants are attempting to derive a transcendental equation for the energy of the bound state and are discussing the implications of the neighboring wall on the potential strength C. There are questions about the boundary conditions that must be satisfied at the delta function location and the approach to integrating the Schrödinger equation around that point.

Discussion Status

Some participants have suggested alternative methods for treating the delta potential directly rather than approximating it with a finite square well. There is ongoing exploration of the boundary conditions and how they affect the wave function's behavior at the delta function location. The discussion is active, with multiple interpretations being considered.

Contextual Notes

Participants are considering the limit as the width of the delta function approaches zero and the implications of the potential's strength on the existence of bound states. There is a focus on the mathematical formulation and the physical reasoning behind the conditions for bound states.

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 Schrödinger 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   Reactions: Wavefunction
vela said:
The idea is to integrate the Schrödinger 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   Reactions: Wavefunction
Awesome thanks for help everyone! :) Cheers
 

Similar threads

Replies
5
Views
3K
Replies
2
Views
2K
Replies
7
Views
2K
Replies
1
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 20 ·
Replies
20
Views
3K
  • · Replies 0 ·
Replies
0
Views
2K
  • · Replies 2 ·
Replies
2
Views
4K
Replies
7
Views
3K
  • · Replies 16 ·
Replies
16
Views
3K