One Dimensional Finite Depth Square Well
Michael Fowler, University of Virginia
We have considered in some detail a particle trapped between infinitely high walls a distance L apart, we have found the wave function solutions of the time independent Schrödinger equation, and the corresponding energies. The essential point was that the wave function had to go to zero at the walls, because there is zero probability of finding the particle penetrating an infinitely high wall. This meant that the lowest energy state couldn't have zero energy, that would give a constant nonzero wave function. Rather, the lowest energy state had to have the minimal amount of bending of the wave function necessary for it to be zero at the two walls but nonzero in between-this corresponds to half a period of a sine or cosine (depending on the choice of origin), these functions being the solutions of Schrödinger's equation in the zero potential region between the walls. The sequence of wave functions (eigenstates) as the energy increases have 0, 1, 2, … zeros (nodes) in the well.
Let us now consider how this picture is changed if the potential at the walls is not infinite. It will turn out to be convenient to have the origin at the center of the well, so we take
V(x) = V0 for x < -L/2
V(x) = 0 for -L/2 < x < L/2
V(x) = V0 for L/2 < x.
Having the potential symmetric about the origin makes it easier to catalog the wave functions. For a symmetric potential, the wave functions can always be taken to be symmetric or antisymmetric. (If a wave function ψ(x) is a solution of Schrödinger's equation with energy E, and the potential is symmetric, then ψ(-x) is a solution with the same energy. This means that ψ(x)+ ψ(-x) and ψ(x)- ψ(-x) are also solutions, since the equation is linear, and these are symmetric and antisymmetric respectively, and using them is completely equivalent to using the original ψ(x) and its reflection ψ(-x).)
How is the lowest energy state wave function affected by having finite instead of infinite walls? Inside the well, the solution to Schrödinger's equation is still of cosine form (it's a state symmetric about the origin). However, since the walls are now finite, ψ(x) cannot change slope discontinuously to a flat line at the walls. It must instead connect smoothly with a function which is a solution to Schrödinger's equation inside the wall.
The equation in the wall is
and has two exponential solutions (say, for x > L/2) one increasing to the right, the other decreasing,
(We are assuming here that E < V0, so the particle is bound to the well. We shall find this is always true for the lowest energy state.)
Let us try to construct the wave function for the energy E corresponding to this lowest bound state. From the equation with V0 = 0, the wavefunction inside the well (let's assume it's symmetric for now) is proportional to coskx, where
The wave function (and its derivative!) inside the well must match a sum of exponential terms—the wave function in the wall—at x = L/2, so
(By writing just a cosine term inside the well, we have left out the overall normalization constant. This can be put back in at the end.)
Solving these equations for the coefficients A, B in the usual way, we find that in general the cosine solution inside the well goes smoothly into a linear combination of exponentially increasing and decreasing terms in the wall. (By the symmetry of the problem, the same thing must happen for x < -L/2.) However, this cannot in general represent a bound state in the well. The increasing solution increases without limit as x goes to infinity, so since the square of the wave function is proportional to the probability of finding the particle at any point, the particle is infinitely more likely to be found at infinity than anywhere else. It got away! This clearly makes no sense—we're trying to find wave functions for particles that stay in, or at least close to, the well. We are forced to conclude that the only exponential wave function that makes sense is the one for which A is exactly zero, so that there is only a decreasing wave in the wall.
Requiring the decreasing wave function, A = 0, means that only a discrete set of values of k, or E, satisfy the boundary condition equations above. They are most simply found by taking A = 0 and dividing one equation by the other to give:
This cannot be solved analytically, but is easy to solve graphically by plotting the two sides as functions of k (recall , and) and finding where the curves intersect.
Solving Schrödinger's Equation Numerically
A more straightforward approach is to solve the Schrödinger equation numerically, using a spreadsheet, and explore how the wave function varies as E is changed. This method has the great advantage that the spreadsheet can be readily adapted to any reasonable potential. We present first the simplest, least sophisticated method, which gives an accuracy of around one percent for energy values using a few hundred rows of the spreadsheet.
To simplify adapting the spreadsheet to the equation, we choose units so that
From now on, in this section we shall denote the wave function by f(x) or just f instead of ψ(x) , to lessen headaches in setting up the spreadsheet.
Schrödinger's Equation is then:
To solve this using the spreadsheet, it is necessary to discretize the derivative. We replace the continuous variable x by a sequence of equally spaced points dx apart, such as 0, dx, 2dx, 3dx, … . Here dx is of course finite, we should really call it D x.
We approximate the second derivative by:
so the differential equation is replaced by the difference equation:
Notice that from this equation, if we know f(xj), f(xj-1) and V(xj) we can find f(xj+1). Of course, we do know V(xj), this is the potential we are trying to solve for. We can also specify E, the energy of the particle, although we shall find that for negative E, only a discrete set of bound state energies correspond to normalizable wave functions.
A second-order differential equation like this needs two boundary conditions for the solution to be defined. We can specify the initial value of f and its derivative df/dx. But this is just equivalent to giving the first two members of the discrete series, f(x0) and f(x1), say. We can then use the difference equation to find f(x2), then use it again to find f(x3) and so on. This is what the spreadsheet does for us: we find f(xj+1) row by row using the difference equation above written:
Using the Spreadsheet for the
Recall that by taking the origin in the center of the square well, we argued that we need only look at wave functions that were symmetric or antisymmetric about the origin. If you try to sketch such wave functions, you will find that symmetric wave functions must have zero slope at the origin, and antisymmetric wave functions must be zero at the origin.
Furthermore, if we know the wave function in the right-hand half, that is, for x > 0, we know it for all x, from the symmetry. Hence, we need only solve Schrödinger's equation going from the origin to the right - we can take x0 = 0, x1 = dx, etc.
As we discussed in the preceding section, on integrating out from the origin with energy E set less than V0, once we cross over into the wall the solution is some combination of an exponentially increasing function and an exponentially decreasing function,
where the exponential coefficient α is positive, and depends on E—it was defined in the preceding section.
Let us first look at the symmetric solutions for very low energies, so take f(0) = 1, f'(0) =0. (Note that we cannot find the correct overall normalization constant until we find the solution, then integrate its square over all space—this can always be done later, and is unnecessary for analyzing the properties of the state).
Let us begin with the trivial case E = 0. For zero energy, inside the well E - V will of course be identically zero, so from Schrödinger's equation the slope of f(x) can never change, consequently f(x) = 1 for all x < L. On reaching the wall, this wave function and its derivative connect smoothly from inside to outside if A = B = ½. It is clear that as we keep going to the right, the A term in the above equation dominates and f(x) diverges, signaling that there is no state localized in the well at E = 0.
At this point, you should create an Excel spreadsheet, with a graph showing f(x) and V(x)-E as functions of x. Details of how this might be done are given in the accompanying homework assignment. If you have a working spreadsheet, it will be much easier to understand the following discussion!
As we now increase E from zero, the symmetric wave function, having zero slope at the center of the well, will begin to curve downwards on moving away from the center, and as the energy increases so does the downward curvature. This naturally changes the mix of increasing and decreasing exponentials needed to connect smoothly at the wall, and in fact as a function of E, A(E) changes sign at a certain value we will call E0. For energies just below E0, f(x) diverges to plus infinity for large x. For energies above E0, it diverges to minus infinity for large x. Exactly at E0, f(x) goes to zero for large x. This is the wave function we are looking for: it corresponds to a particle localized close to the well, and in fact is the lowest possible energy—the ground state—for a particle in the well. E0 is called the ground state eigenvalue, the wave function is called an eigenstate.
Finding the Ground State Energy
The first problem is to find this value E0. If we just guess a value of E, the wave function will almost certainly diverge for large x. The way to find E0 is to notice that for E below E0, f(x) goes to large positive values on the far right, for E just above E0, f(x) goes to large negative values on the far right. So we take an increasing set of E values starting near zero, and watch for the tail to wag! When this happens, we back up half way, then back or forward as necessary, choosing a set of E values that bracket E0.
Actually, this rather tedious process can be automated in Excel, but it's worth going through it manually a time or two to get a feeling for how the wave functions behave. A general feature to check is that for E > V(x) they are oscillatory, sine or cosine like, although with changing wavelengths in general; whereas for E < V(x) they have exponential character. You can see why this is by looking at the equation, but it’s worth bearing in mind as you examine the curves.
Automating the Process
Now to a bit of automation. The strategy we are following is to adjust E until the function f(x) goes to zero for large x, that is, far down at the bottom of our column having several hundred entries. Since we would like to know how well we are doing, it’s convenient to have the spreadsheet copy that value to some cell near the top. So, if the end of our column of f(xi)'s is F600, say, then we enter =$F$600 in some convenient square near the top, say B22. Now we can keep an eye on the distant value as we adjust E. We can also get Excel to do this for us as follows (I assume the value of E is entered into square B15):
Click Tools, Goal Seek. Then fill out the dialog box: Set Cell B22 to value 0 by changing cell B15. (You can fill in the cell addresses by just clicking on the appropriate cells.) Click OK. It will then calculate and come up with a value. Click OK again.
The curve below was generated by this procedure. It is worth noting that there is a definite probability of finding the particle inside the wall, that is to say, the wave function is nonzero there (and, it turns out, for reasonable, physical values of the parameters):
All square well potentials in one dimension, however shallow, have a localized ground state with this general shape. Whether or not there are other eigenstates with other eigenvalues depends on the depth of the potential. For a sufficiently shallow potential, there is only one state. An infinitely deep well, as we discussed earlier, has an infinite number of bound states. As the well depth increases from zero, states are bound sequentially. These higher eigenstates are called excited states. A particle in one of them will usually decay to a lower state, emitting a photon, just as in the Bohr atom.
Finding Excited States
Let us now see how the wave function develops as E increases beyond E0 for a sufficiently deep well. As the energy is increased, the cosine term inside the well has tighter curvature, and the exponentially increasing term for large x is large and negative. However, continuing to increase the energy, eventually the cosine term bottoms out inside the well and begins to turn up again before reaching the wall. At a certain energy, it again becomes possible to match it to a decreasing function:
To find this new eigenvalue, we can get close to it then use: Tools, Goal Seek just as we did for the ground state—see above.
Creating a Macro
To really automate this procedure, it's worth creating a macro. So, go to Tools, Record Macro, Record New Macro. Call the Macro "Eigenvalue", and put in the shortcut Ctrl+e. In the Record Macro box, click OK to begin recording. Then Click Tools, Goal Seek. Then fill out the dialog box just as we did before, that is, Set Cell B22 to value 0 by changing cell B15. Click OK. It will then calculate and come up with a value. Click OK again. That is all we want to record, so click Tools, Macro, Stop Recording.
You can check to see if your macro is working by changing the value of E in B15, then hit Ctrl+e. (If you forgot to assign the shortcut to your macro, you can do it now by clicking Tools, Macro, Macros, click on your macro and click Options. You can fill in the shortcut key in the Options box.)
To get really fancy, it’s nice to run this macro by clicking on something. Insert a textbox somewhere in your spreadsheet, and write “Find Eigenvalue” on it. With the textbox selected, right click on one of the sizing handles and click “Assign Macro” . Choose your “Eigenvalue” macro. Now it will run if you click on the textbox.
A point to notice about the wave function pictured above is that has a node at about x = 2. (Actually, of course, this is just half the wave function for the complete well-there will be another node at x = -2.) As we continue to increase E, the high-x tail of the wave function wags up and down. Each time it crosses the axis there is an allowable wave function. Furthermore, each time it crosses the axis the wave function collects another node for x > 0. Thus, the complete wave functions generated by this method have 0, 2, 4, 6, … nodes. From the symmetry of the problem, the allowable wave functions with an odd number of nodes must have one node at the origin. They can be generated by taking as initial conditions that f is zero at the origin, and has finite slope.
Curvature of Wave Functions
Schrödinger's equation in the form
can be interpreted by saying that the left-hand side, the rate of change of slope, is the curvature—so the curvature of the function is equal to (V(x) - E)f(x). This means that if E > V(x), for f(x) positive f(x) is curving negatively, for f(x) negative f(x) is curving positively. In both cases, f(x) is always curving towards the axis. This means that for E > V(x), f(x) has a kind of stability: its curvature is always bringing it back towards the axis, so it has oscillatory character.
On the other hand, for V(x) > E, the curvature is always away from the axis. This means that f(x) tends to diverge to infinity. Only under exactly the right conditions will this curvature be just enough to bring the wave function to zero as x goes to infinity. (As f(x) tends to zero, the curvature tends to zero, too.)
It is worth examining the wave functions generated by the spreadsheet to see just how the curvature changes as V(x) - E, or for that matter f(x), changes sign.