The finite potential well (also known as the finite square well) is a concept from quantum mechanics. It is an extension of the infinite potential well, in which a particle is confined to a "box", but one which has finite potential "walls". Unlike the infinite potential well, there is a probability associated with the particle being found outside the box. The quantum mechanical interpretation is unlike the classical interpretation, where if the total energy of the particle is less than the potential energy barrier of the walls it cannot be found outside the box. In the quantum interpretation, there is a non-zero probability of the particle being outside the box even when the energy of the particle is less than the potential energy barrier of the walls (cf quantum tunnelling).
Particle in a 1-dimensional box
For the 1-dimensional case on the x-axis, the time-independent Schrödinger equation can be written as:
![{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi }{dx^{2}}}+V(x)\psi =E\psi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/8c2fa52f362e1ce6c9caeb6f319c7a5d63ba9376)
|
|
(1)
|
where
is the reduced Planck's constant,
is Planck's constant,
is the mass of the particle,
is the (complex valued) wavefunction that we want to find,
is a function describing the potential energy at each point x, and
is the energy, a real number, sometimes called eigenenergy.
For the case of the particle in a 1-dimensional box of length L, the potential is
outside the box, and zero for x between
and
. The wavefunction is considered to be made up of different wavefunctions at different ranges of x, depending on whether x is inside or outside of the box. Therefore, the wavefunction is defined such that:
![{\displaystyle \psi ={\begin{cases}\psi _{1},&{\text{if }}x<-L/2{\text{ (the region outside the box)}}\\\psi _{2},&{\text{if }}-L/2<x<L/2{\text{ (the region inside the box)}}\\\psi _{3},&{\text{if }}x>L/2{\text{ (the region outside the box)}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5275282c84323a27453ef1d3dd38093332cf50c3)
Inside the box
For the region inside the box, V(x) = 0 and Equation 1 reduces to
![{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi _{2}}{dx^{2}}}=E\psi _{2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aec75d9b4f6f0b361dc6de74b0211b1200bcce8b)
Letting
![{\displaystyle k={\frac {\sqrt {2mE}}{\hbar }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/90edd443def7d3e69c0eb8d4525a61b03afafd03)
the equation becomes
![{\displaystyle {\frac {d^{2}\psi _{2}}{dx^{2}}}=-k^{2}\psi _{2}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/df32011b45e2184a36f9de0ef6a16241b58e21b4)
This is a well-studied differential equation and eigenvalue problem with a general solution of
![{\displaystyle \psi _{2}=A\sin(kx)+B\cos(kx)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b88f62e3f17b480ff11d90886efe69e6e07b246d)
Hence,
![{\displaystyle E={\frac {k^{2}\hbar ^{2}}{2m}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab6e0ac7528fb684e630e935a53cd53737fcd4b9)
Here, A and B can be any complex numbers, and k can be any real number.
Outside the box
For the region outside of the box, since the potential is constant,
and equation 1 becomes:
![{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}\psi _{1}}{dx^{2}}}=(E-V_{0})\psi _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5a437dad26494d3aba7e2722eb79a006b4e9a6d1)
There are two possible families of solutions, depending on whether E is less than
(the particle is bound in the potential) or E is greater than
(the particle is free).
For a free particle,
, and letting
![{\displaystyle k'={\frac {\sqrt {2m(E-V_{0})}}{\hbar }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/edb900480575aa9229d4dff1bdc651f98113221b)
produces
![{\displaystyle {\frac {d^{2}\psi _{1}}{dx^{2}}}=-k'^{2}\psi _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/139ad8ded187dd8e6b6fdb535971d5a2b5f7f504)
with the same solution form as the inside-well case:
![{\displaystyle \psi _{1}=C\sin(k'x)+D\cos(k'x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/203ca4fbb1d44486f0d7dbbdfc0133c96f469849)
This analysis will focus on the bound state, where
. Letting
![{\displaystyle \alpha ={\frac {\sqrt {2m(V_{0}-E)}}{\hbar }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/db99d16c87906a0d52fe0f60185696b454c90529)
produces
![{\displaystyle {\frac {d^{2}\psi _{1}}{dx^{2}}}=\alpha ^{2}\psi _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4fb3569c7cf0bc092997cf6fbf9646ef5fc2fa1f)
where the general solution is exponential:
![{\displaystyle \psi _{1}=Fe^{-\alpha x}+Ge^{\alpha x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0ea399b82cab1c1ccdeaec40aaa5d85acedb4d14)
Similarly, for the other region outside the box:
![{\displaystyle \psi _{3}=He^{-\alpha x}+Ie^{\alpha x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8a0a6b1e65ad3331a189cd16b8ba94c1f4c7c48)
Now in order to find the specific solution for the problem at hand, we must specify the appropriate boundary conditions and find the values for A, B, F, G, H and I that satisfy those conditions.
Finding wavefunctions for the bound state
Solutions to the Schrödinger equation must be continuous, and continuously differentiable.[1] These requirements are boundary conditions on the differential equations previously derived, that is, the matching conditions between the solutions inside and outside the well.
In this case, the finite potential well is symmetrical, so symmetry can be exploited to reduce the necessary calculations.
Summarizing the previous sections:
![{\displaystyle \psi ={\begin{cases}\psi _{1},&{\text{if }}x<-L/2{\text{ (the region outside the box)}}\\\psi _{2},&{\text{if }}-L/2<x<L/2{\text{ (the region inside the box)}}\\\psi _{3}&{\text{if }}x>L/2{\text{ (the region outside the box)}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f57d414da936d753926e93c089cc3693777c28bb)
where we found
![{\displaystyle \psi _{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8cfdde1da54e02a016fe2a230c58b25dfcc014d6)
,
![{\displaystyle \psi _{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c5083a526766f85c4f39ab695791b0b739f06897)
, and
![{\displaystyle \psi _{3}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/47acb83b8aed59a815c7a6f3bd2eaed2c932336d)
to be:
![{\displaystyle {\begin{aligned}\psi _{1}&=Fe^{-\alpha x}+Ge^{\alpha x}\\\psi _{2}&=A\sin(kx)+B\cos(kx)\\\psi _{3}&=He^{-\alpha x}+Ie^{\alpha x}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ad6a4e878b1cd578061cac731a82390e384c4ad)
We see that as
goes to
, the
term goes to infinity. Likewise, as
goes to
, the
term goes to infinity. In order for the wave function to be square integrable, we must set
, and we have:
![{\displaystyle \psi _{1}=Ge^{\alpha x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/12f1cbf01a6c0d9bc275541153ecf38336c3309e)
and
![{\displaystyle \psi _{3}=He^{-\alpha x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/31ca6e16da98dc82c6b4186a617333f65252ffa9)
Next, we know that the overall
function must be continuous and differentiable. In other words, the values of the functions and their derivatives must match up at the dividing points:
![{\displaystyle \psi _{1}(-L/2)=\psi _{2}(-L/2)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49d29364e8f15c9aadac04e9377561b2022f78e6) |
|
|
![{\displaystyle \left.{\frac {d\psi _{1}}{dx}}\right|_{x=-L/2}=\left.{\frac {d\psi _{2}}{dx}}\right|_{x=-L/2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e1b808b873ddbbd682b61328f8d791afba509b77) |
|
|
These equations have two sorts of solutions, symmetric, for which
and
, and antisymmetric, for which
and
. For the symmetric case we get
![{\displaystyle He^{-\alpha L/2}=B\cos(kL/2)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/46617bb750d0a768b3955a8c3e34f1c50e79efd6)
![{\displaystyle -\alpha He^{-\alpha L/2}=-kB\sin(kL/2)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd34da57a97c0e2d28abff536b3d40d985b70990)
so taking the ratio gives
![{\displaystyle \alpha =k\tan(kL/2).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/56f8c4965f2afba4972d9a6b826f3eb2afacddb9)
Similarly for the antisymmetric case we get
![{\displaystyle \alpha =-k\cot(kL/2).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/96673c71a4025dc6ed2ce69e17f2b0358a02169c)
Recall that both
and
depend on the energy. What we have found is that the continuity conditions cannot be satisfied for an arbitrary value of the energy; because that's a result of the infinite potential well case. Thus, only certain energy values, which are solutions to one or either of these two equations, are allowed. Hence we find that the energy levels of the system below
are discrete; the corresponding eigenfunctions are bound states. (By contrast, for the energy levels above
are continuous.[2])
The energy equations cannot be solved analytically. Nevertheless, we will see that in the symmetric case, there always exists at least one bound state, even if the well is very shallow.[3]
Graphical or numerical solutions to the energy equations are aided by rewriting them a little. If we introduce the dimensionless variables
and
, and note from the definitions of
and
that
, where
, the master equations read
![{\displaystyle {\sqrt {u_{0}^{2}-v^{2}}}={\begin{cases}v\tan v,&{\text{(symmetric case) }}\\-v\cot v,&{\text{(antisymmetric case) }}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80fb76aa16930e7bb4ddd828abc40f7da976bcb6)
In the plot to the right, for
, solutions exist where the blue semicircle intersects the purple or grey curves (
and
). Each purple or grey curve represents a possible solution,
within the range
. The total number of solutions,
, (i.e., the number of purple/grey curves that are intersected by the blue circle) is therefore determined by dividing the radius of the blue circle,
, by the range of each solution
and using the floor or ceiling functions:[4]
![{\displaystyle N=\left\lfloor {\frac {2u_{0}}{\pi }}\right\rfloor +1=\left\lceil {\frac {2u_{0}}{\pi }}\right\rceil }](https://wikimedia.org/api/rest_v1/media/math/render/svg/27f30bf874aa197be5a6d18ec78eb7d09012347c)
In this case there are exactly three solutions, since
.
and
, with the corresponding energies
![{\displaystyle E_{n}={2\hbar ^{2}v_{n}^{2} \over mL^{2}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10768c1888bf3ea8158ce8f874f461af9a18d087)
If we want, we can go back and find the values of the constants
![{\displaystyle A,B,G,H}](https://wikimedia.org/api/rest_v1/media/math/render/svg/363c8e6830b77c7a7e2338af9070526ea409b97b)
in the equations now (we also need to impose the normalisation condition). On the right we show the energy levels and wave functions in this case (where
![{\textstyle x_{0}\equiv \hbar /{\sqrt {2mV_{0}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55cc04b7c6b99039ffd13071c0677eb7a27ef073)
):
We note that however small
is (however shallow or narrow the well), there is always at least one bound state.
Two special cases are worth noting. As the height of the potential becomes large,
, the radius of the semicircle gets larger and the roots get closer and closer to the values
, and we recover the case of the infinite square well.
The other case is that of a very narrow, deep well - specifically the case
and
with
fixed. As
it will tend to zero, and so there will only be one bound state. The approximate solution is then
, and the energy tends to
. But this is just the energy of the bound state of a Delta function potential of strength
, as it should be.
A simpler graphical solution for the energy levels can be obtained by normalizing the potential and the energy through multiplication by
. The normalized quantities are
![{\displaystyle {\tilde {V}}_{0}=V_{0}{\frac {8m}{h^{2}}}L^{2}\qquad {\tilde {E}}=E{\frac {8m}{h^{2}}}L^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67de5dad02c23cf741682c4164e9d29353a78aa3)
giving directly the relation between the allowed couples
![{\displaystyle (V_{0},E)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/128f3bab14f8729aa384a32763c8881cdfc7f11e)
as
[5]
![{\displaystyle {\sqrt {{\tilde {V}}_{0}}}={\sqrt {\tilde {E}}}\,\left|{\sec({\sqrt {\tilde {E}}}\,{\pi }/{2})}\right|,\qquad {\sqrt {{\tilde {V}}_{0}}}={\sqrt {\tilde {E}}}\,\left|{\csc({\sqrt {\tilde {E}}}\,{\pi }/{2})}\right|}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dbca7892f4710900ac20c6a6cdcba382f536d898)
for the even and odd parity wave functions, respectively. In the previous equations only the positive derivative parts of the functions have to be considered. The chart giving directly the allowed couples
![{\displaystyle (V_{0},E)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/128f3bab14f8729aa384a32763c8881cdfc7f11e)
is reported in the figure.
Unbound states
If we solve the time-independent Schrödinger equation for an energy
, the solutions will be oscillatory both inside and outside the well. Thus, the solution is never square integrable; that is, it is always a non-normalizable state. This does not mean, however, that it is impossible for a quantum particle to have energy greater than
, it merely means that the system has continuous spectrum above
. The non-normalizable eigenstates are close enough to being square integrable that they still contribute to the spectrum of the Hamiltonian as an unbounded operator.[6]
Asymmetrical well
Consider a one-dimensional asymmetrical potential well given by the potential[7]
![{\displaystyle V(x)={\begin{cases}V_{1},&{\text{if }}-\infty <x<0{\text{ (the region outside the box)}}\\0,&{\text{if }}0<x<a{\text{ (the region inside the box)}}\\V_{2}&{\text{if }}a<x<\infty {\text{ (the region outside the box)}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b6e883e61e16c29d3cda8b0447c7c1b7d836341e)
with
. The corresponding solution for the wave function with
is found to be
![{\displaystyle \psi (x)={\begin{cases}c_{1}e^{k_{1}x},&{\text{for }}x<0,{\text{where }}k_{1}={\sqrt {(2m/\hbar ^{2})(V_{1}-E)}}\\c\sin(kx+\delta ),&{\text{for }}0<x<a,{\text{where }}k={\sqrt {2mE/\hbar ^{2}}}\\c_{2}e^{-k_{2}x},&{\text{for }}x>a,{\text{where }}k_{2}={\sqrt {(2m/\hbar ^{2})(V_{2}-E)}}\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/10e3a7a93d989777c24b2f5ac9ebb9c1c6408c9d)
and
![{\displaystyle \sin \delta ={\frac {k\hbar }{\sqrt {2mV_{1}}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/330452f8f0a4af997a6f61a6d11a2d9d818bc233)
The energy levels
are determined once
is solved as a root of the following transcendental equation
![{\displaystyle ka=n\pi -\sin ^{-1}\left({\frac {k\hbar }{\sqrt {2mV_{1}}}}\right)-\sin ^{-1}\left({\frac {k\hbar }{\sqrt {2mV_{2}}}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4cbdd827612c9e01641ba9ae50cffdfdc7936e01)
where
Existence of root to above equation is not always guaranteed, for example, one can always find a value of
so small, that for given values of
and
, there exists no discrete energy level. The results of symmetrical well is obtained from above equation by setting
.
Spherical cavity
The results above can be used to show that, contrary to the one-dimensional case, there is not always a bound state in a spherical cavity.
The ground state (n = 1) of a spherically symmetric potential will always have zero orbital angular momentum (ℓ = n−1), and the reduced wave function
satisfies the equation
![{\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {d^{2}U}{dr^{2}}}+V(r)U(r)=EU(r)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/902d68cbe579d61406782d7bd1c32855f1f48e27)
This is identical to the one-dimensional equation, except for the boundary conditions. As before,
![{\displaystyle U(r)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/354eeb1c2804f4b2bbc06122ce45472a61682b09)
and its first derivative must be continuous at the edge of the well
![{\displaystyle r=R}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a088bb7b6b4be9ff6b0fdb33fc1dd53af91e356)
. However, there is another condition, that
![{\displaystyle \psi (0)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fcfdf3fefa6feee8b2192815ac15a2650d8a8db2)
must be finite, and that requires
![{\displaystyle U(0)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/89b1fc1082774b08c5e4ab211a7bfbaafb5eb82c)
.
By comparison with the solutions above, we can see that only the antisymmetric ones have nodes at the origin. Thus only the solutions to
are allowed. These correspond to the intersection of the semicircle with the grey curves, and so if the cavity is too shallow or small, there will be no bound state.
See also
References
Further reading