Required math: calculus
Reference: Griffiths, David J. (2005), Introduction to Quantum Mechanics, 2nd Edition; Pearson Education – Sec 2.2, Problem 2.7.
In the solution of the one-dimensional particle in a box problem we find that there are an infinite number of wave functions that satisfy the spatial part of the Schrödinger equation, and that each such solution corresponds to a discrete energy. To summarize, the problem was to solve the spatial equation
for the potential
We found that the solutions, properly normalized, are
with outside the box. The energy corresponding to solution is
But this solves only the spatial part of the Schrödinger equation. The full solution requires bringing back the function containing the time dependence that arises from solving the equation using separation of variables. Thus the full, time-dependent solution for a particular energy is
where and are given above.
The probability of finding a particle that is in state at location at time is therefore . That is, if a particle is in one of the states with a definite energy , its probability density is independent of time, since the only time dependence comes from the complex exponential function whose modulus is always 1. The time dependence always disappears when the square modulus is calculated. For this reason, the states are called stationary states.
This isn’t the end of the story, however. Since the Schrödinger equation is linear, once we have found two or more solutions, any linear combination of these solutions gives another solution. Thus the most general solution of the time-independent Schrödinger equation for the particle in a box is
where the are arbitrary complex constants. Solutions that are linear combinations of two or more of the stationary states are, however, not stationary states themselves, since each term in the sum contains its own dependence on time in the form of a complex exponential, and these exponentials will not cancel out when the square modulus is calculated.
To explore the consequences of this general solution, we need first to demonstrate some of the properties of the stationary states. We know that is normalized, in the sense that
However the set of functions is also orthogonal, in the sense that if
(This can be checked by direct integration, using the trigonometric identity .) A set of functions that are both orthogonal and normalized is called an orthonormal set.
This result may surprise you, but you might also be thinking ‘so what?’. The normalization is clearly important if is to be used as a probability density, but the orthogonality seems pretty much irrelevant. To see why this property is important, suppose you want a full dynamical description of a particle starting from time when the particle is in some initial state . In general, we should be able to specify this initial state to be anything at all. In the classical case, an analogous situation might be something like this: standing on the Earth’s surface we have a rock of mass which we can throw from any height above the ground. We can aim the rock in any direction and at any speed. All of these things are independent of the law of gravity which takes over after the rock has been thrown. But the theory of motion that we use to describe the motion of the rock must allow us to specify the initial conditions to be anything we like.
It is a similar situation in quantum mechanics. The Schrödinger equation tells us how the particle behaves, but only after we have specified the initial conditions. In this case, we can ‘throw’ the particle into the potential well by specifying its initial state, and the theory must be able to accept any (within the constraints of the problem) initial state. For the particle in a box, since we have an infinite potential outside the box, there is no way the particle could ever be found outside, so any initial condition must confine the particle to the region , but apart from that, can be anything at all.
From the general solution above, we can get an expression for :
The problem is therefore, given and , find the constants . It is here that the orthonormal nature of the set of stationary solutions comes into play. Suppose we want to find for some particular value of . If we multiply this expression by and integrate over the range of the box (from 0 to ), then all integrals where are zero because is orthogonal to , so we are left with just the one term we want:
where the symbol is called the Kronecker delta and is a shorthand way of writing a quantity that is zero if , and 1 if .
So for any initial condition, we now have a way of writing it as a sum over the spatial part of the stationary states at time . The general solution is obtained merely by restoring the time dependence:
Note that I am not saying that any function is a solution of the Schrödinger equation – that is clearly absurd. What is true, however, is that any function of can be expressed as a linear combination of the solutions of the spatial part of the Schrödinger equation. It is important to note that any of the functions on its own is not a solution of the Schrödinger equation unless it is multiplied by , so the linear combination is not in general a solution of the Schrödinger equation. It is only the final form that is a solution.
Example: A particle in the infinite square well starts off with a wave function as follows:
The initial wave function is thus a triangle with its peak at and a height of . We can normalize the wave function first
Next, we can find in terms of the . We need to calculate the , so:
(The integral is straightforward, if a bit tedious. I used mathematical software to do it on a computer.) Thus we get
The term is zero if is even, and for odd values of , giving the sums shown in the last line. The series consists of odd only, with half the terms being positive and the other half negative. The full solution is found by replacing the exponentials containing the time dependence:
The probability that the energy is is
[Incidentally, since we know that we get from 20 that
The average energy can be found by using the energy levels for the infinite square well:
We then have
is not commonly found in tables, but a link (live at the time of writing) showing a derivation is here. The derivation requires knowledge of either the Riemann zeta function or the residue theorem from complex variable theory.
To summarize, the situation with the particle in a box is:
- Using separation of variables, we can split the solution into a spatial function and a temporal factor .
- When we solve the spatial equation with appropriate boundary conditions and normalization, we get the set of functions . Each of these functions is associated with a specific discrete energy .
- The set of functions is orthonormal.
- The full solution of the Schrödinger equation for a particular energy is . The square modulus of this solution is independent of time, so it is a stationary state.
- The general solution of the Schrödinger equation is a linear combination of the set of stationary states. The square modulus of a general solution will not, in general, be stationary (time-independent).
- The constants in the general solution can be obtained using the initial condition and the orthonormality of the .