Particle in a box

From Wikipedia, the free encyclopedia

In physics, the particle in a box (also known as the infinite potential well or the infinite square well) is a very simple problem consisting of a single particle bouncing around inside of an immovable box, from which it cannot escape, and which loses no energy when it collides with the walls of the box. In classical mechanics, the solution to the problem is trivial: The particle moves in a straight line, always at the same speed, until it reflects from a wall. When it reflects from a wall, it always reflects at an equal but opposite angle to its angle of approach, and its speed does not change.

The problem becomes very interesting when one attempts a quantum-mechanical solution, since many fundamental quantum mechanical concepts need to be introduced in order to find the solution. Nevertheless, it remains a very simple and solvable problem. This article will only be concerned with the quantum mechanical solution.

The problem may be expressed in any number of dimensions, but the simplest problem is one dimensional, while the most useful solution is the particle in the three dimensional box. In one dimension this amounts to the particle existing on a line segment, with the "walls" being the endpoints of the segment.

In physical terms, the particle in a box is defined as a single point particle, enclosed in a box inside of which it experiences no force whatsoever, i.e. it is at zero potential energy. At the walls of the box, the potential rises to infinity, forming an impenetrable wall. Using this description in terms of potentials allows the Schrödinger equation to be used to determine the solution.

As mentioned above, if we were studying this system under the rules of classical mechanics we would apply Newton's laws of motion to the initial conditions and the result would seem reasonable and intuitive. In quantum mechanics, when the Schrödinger equation is applied to the proposed system, the results are not intuitive. In the first place, the particle can only have certain specific energy levels, and the zero energy level is not one of them. Secondly, the chances of detecting the particle in the box at any specific energy level is not uniform - there are certain locations in the box where the particle might be found, but there are also places where it can never be found. Both of these results differ from the usual way we perceive the world, yet rest on principles that have been extensively experimentally verified.

Contents

[edit] Formal introduction

The particle in a box (or the infinite potential well or infinite square well) is a simple idealized system that can be completely solved within quantum mechanics. The infinite potential well is a finite sized region in space (the box) with an infinite potential at its boundaries (the walls). A particle experiences no forces while inside the box, but as the walls are 'infinitely high', it is constrained to remain in the box. This is similar to the situation of a gas confined in a non-porous container. For simplicity we start with the 1-dimensional case, where the motion of the electron is considered only in a single dimension - the direction of the quantum confinement. Later we extend the discussion to the 2 and 3 dimensional cases.

As we shall see, the solution of the Schrödinger equation for the particle in a box problem reveals some decidedly quantum behavior of the particle that agrees with observation but contrasts sharply with the predictions of classical mechanics. This is a particularly useful illustration because this behaviour is not "forced" on the system, it arises naturally from the initial conditions. It neatly demonstrates that quantum behaviour is a natural outcome of any wave-like system, contrary to the common concept of a "quantum leap" where the behavior is almost magical.

The quantum behavior in the box includes:

  1. Energy quantization - It is not possible for the particle to have any arbitrary definite energy. Instead only discrete definite energy levels are allowed (if the state is not a steady state, however, any energy past zero-point energy is allowed on average).
  2. Zero-point energy - The lowest possible energy level of the particle, called the zero-point energy, is nonzero.
  3. Spatial nodes - In contrast to classical mechanics the Schrödinger equation predicts that for some energy levels there are nodes, implying positions at which the particle can never be found.

One can solve the Schrödinger equation analytically for such a simple potential. However trivial, this case is both of great technical value for the insights it allows, and of paramount physical importance. Depending on the boundary conditions, one can use the solutions to describe two important systems. If one considers real valued solutions (of which detailed derivation is given below), one describes actual potentials of heterostructures called quantum wells which spatially confine particles - typically electrons and holes. If one considers complex valued solutions, one describes conveniently a particle propagating freely in a constrained volume (like a solid).

[edit] 1-dimensional box

[edit] Solution

The eigenfunction and eigenenergies of a particle of mass m in a one dimensional box of length L are

\psi_n = \sqrt{\frac{2}{L}} \sin{\left(\frac{n \pi x}{L} \right)} \,
E_n = \frac{\hbar^2 \pi^2}{2 m L^2} n^2 \,

[edit] Derivation

The problem is defined as the case of a particle situated in a 1D infinite square well with momentum only in the direction of quantum confinement (the x direction).

The 1D time-independent Schrödinger equation can be written as:

-\frac{\hbar^2}{2m} \frac{\mathrm{d}^2 \psi(x)}{\mathrm{d}x^2} + V(x) \psi(x) = E \psi(x) \quad (1)
where
\hbar is the Reduced Planck Constant
m \, is the mass of the particle
\psi\left(x\right)\, is the complex-valued stationary time-independent wavefunction that we want to find
V\left(x\right)\, is the spatially varying potential and
E\, is the energy, a real number.

For the case of the particle in a 1-dimensional box of length L, the potential is zero inside the box, but rises abruptly to infinity at x = 0 and x = L. Thus for the region inside the box V(x) = 0 and Equation 1 reduces to:

The Potential is 0 inside the box, and infinite elsewhere
The Potential is 0 inside the box, and infinite elsewhere
-\frac{\hbar^2}{2m} \frac{\mathrm{d}^2 \psi(x)}{\mathrm{d} x^2} = E \psi(x) \quad (2)

This is a well studied differential equation and eigenvalue problem with a general solution of:

\psi(x) = A \sin(kx) + B \cos(kx)\quad
E = \frac{k^2 \hbar^2}{2m} \quad (3)

Here, A and B can be any complex numbers and the wavenumber, k can be any real number (k must be real because E is real).

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 and B that satisfy those conditions. One usually resorts to one of the following two choices, describing two kinds of systems. The first case, with which we shall pursue our derivation, demands that ψ(x) equal zero at x = 0 and x = L. A handwaving argument to motivate these boundary conditions is that the particle is unlikely to be found at a location with a high potential (the potential repulses the particle), thus the probability of finding the particle, |ψ(x)|2, must be small in these regions and decreases with increasing potential. For the case of an infinite potential, |ψ(x)|2 must be infinitesimally small or 0, thus ψ(x) must also be zero in this region. In summary,

\psi(0)=\psi(L)=0 \quad (4)

Substituting the general solution from Equation 3 into Equation 2 and evaluating at x = 0 (ψ = 0), we find that B = 0 (since sin(0) = 0 and cos(0) = 1). It follows that the wavefunction must be of the form:

\psi(x) = A \sin(kx) \quad (5)

and at x = L we find:

\psi(L) = A \sin(kL) = 0 \quad (6)

One solution for Equation 6 is A = 0, however, this "trivial solution" would imply that ψ = 0 everywhere (i.e. the particle isn't in the box.) and can be thrown out. If A \neq 0 then sin(kL) = 0 only when:

k = \frac{n \pi}{L} \quad \mathrm{where} \quad n\in \mathbb{N}^+ \quad (7)

Note that n = 0 is ruled out because in this case, ψ = 0 everywhere, corresponding to the case where the particle is not in the box. Negative values of n are also neglected, since they merely change the sign of sin(nx). Now in order to find A we must normalise the wavefunction. We recognize that the particle must exist somewhere in space. | ψ(x) | 2 is the probability of finding the particle at a particular point in space, so the integral of this value over all x must be equal to 1:

1 = \int_{-\infty}^{\infty} \left| \psi(x) \right|^2 \, \mathrm{d}x = \left| A \right|^2 \int_0^L \sin^2 \left(kx\right) \, \mathrm{d}x = \left| A \right|^2 \frac{L}{2}
or
\left| A \right| = \sqrt{\frac{2}{L}} \quad (8)

Thus, A may be any complex number with absolute value √(2/L); these different values of A yield the same physical state, so we choose A = √(2/L) to simplify.

Energy Levels (dashed lines) and Wavefunctions (solid lines) of the one-dimensional particle in a box.
Energy Levels (dashed lines) and Wavefunctions (solid lines) of the one-dimensional particle in a box.

Finally, substituting the results from Equations 7 and 8 into Equation 3, the complete set of energy eigenfunctions for the 1-dimensional particle in a box problem is:

\psi_n(x) = \sqrt{\frac{2}{L}} \sin\left(\frac{n \pi x}{L}\right) \quad (9)
E_n = \frac{n^2\hbar^2 \pi ^2}{2mL^2} = \frac{n^2 h^2}{8mL^2} \quad (10)

with

n\in\mathbb{N}^+

Note, that as mentioned previously, only "quantized" energy levels are possible. Also, since n cannot be zero, the lowest energy from Equation 10 is also non-zero. This zero-point energy, as it is called, can be explained in terms of the uncertainty principle. Because the particle is constrained within a finite region, the variance in its position is upper-bounded. Thus due to the uncertainty principle the variance in the particle's momentum cannot be zero, so the particle must contain some amount of energy that increases as the length of the box, L, decreases.

Also, since ψ consists of sine waves, for any value of n greater than one, there are regions within the box for which ψ(x) and thus ψ(x)2 both equal zero, indicating that for these energy levels, nodes exist in the box where the probability of finding the particle is zero.

[edit] 2-dimensional box

[edit] Solution

\psi_{n_x,n_y} = \sqrt{\frac{4}{L_x L_y}} \sin \left( \frac{n_x \pi x}{L_x} \right) \sin \left( \frac{n_y \pi y}{L_y} \right)
E_{n_x,n_y} = \frac{\hbar^2\pi^2}{2m} \left[ \left( \frac{n_x}{L_x} \right)^2 + \left( \frac{n_y}{L_y} \right)^2 \right]

[edit] Derivation

For the 2-dimensional case the particle is confined to a rectangular surface of length Lx in the x-direction and Ly in the y-direction. Again the potential is zero inside the "box" and infinite at the walls. For the region inside the box, where the potential is zero, the two dimensional analogue of Equation 2 applies:

-\frac{\hbar^2}{2m} \left( \frac{\partial^2\psi}{\partial x^2}+\frac{\partial^2 \psi}{\partial y^2} \right) =E\psi \quad(11)

In this case ψ is a function of both x and y, so ψ=ψ(x,y). In order to solve Equation 11, we use the method of separation of variables. First, we assume that ψ can be expressed as the product of two independent functions, the first depending only on x and the second depending only on y; i.e.:

\psi(x,y) = X(x) Y(y) \quad (12)

Substituting Equation 12 into Equation 11 and evaluating the partial derivatives gives:

-\frac{\hbar^2}{2m} \left( Y\frac{\partial^2X}{\partial x^2}+X\frac{\partial^2 Y}{\partial y^2} \right) =E X Y \quad(13)

which upon dividing by XY and rewriting d2X/dx2 as X" and d2Y/dy2 as Y" becomes:

-\frac{\hbar^2}{2m} \left( \frac{X''}{X}+\frac{Y''}{Y} \right) =E \quad(14)

Now we note that since X"/X is independent of y, varying y can only change the Y"/Y term. However, from Equation 14 we see that changing Y"/Y without varying X"/X, would also change E, but E is a constant, so Y"/Y must also be a constant, independent of y. The same argument can be applied to show that X"/X is independent of x. Since X"/X and Y"/Y are constants, we can write:

-\frac{\hbar^2}{2m}\frac{X''}{X} = E_x \quad and \quad -\frac{\hbar^2}{2m}\frac{Y''}{Y} = E_y \quad (15)

where Ex + Ey = E. Expanding X" and Y" in terms of the derivatives and rearranging gives:

-\frac{\hbar^2}{2m}\frac{\partial^2X}{\partial x^2} = E_x X\quad (16)
-\frac{\hbar^2}{2m}\frac{\partial^2Y}{\partial y^2} = E_y Y\quad (17)
The wavefunction of a 2D well with nx=4 and ny=4
The wavefunction of a 2D well with nx=4 and ny=4

each of which are of the same form as the 1-dimensional Schrödinger equation (Equation 2) we solved in the previous section. Thus, adapting the results from the previous section gives:

X_{n_x}=\sqrt{\frac{2}{L_x}} \sin \left( \frac{n_x \pi x}{L_x} \right) \quad (18)
Y_{n_y}=\sqrt{\frac{2}{L_y}} \sin \left( \frac{n_y \pi y}{L_y} \right) \quad (19)

Finally, since ψ=XY and E = Ex + Ey, we obtain the solutions:

\psi_{n_x,n_y} = \sqrt{\frac{4}{L_x L_y}} \sin \left( \frac{n_x \pi x}{L_x} \right) \sin \left( \frac{n_y \pi y}{L_y} \right) \quad (20)
E_{n_x,n_y} = \frac{\hbar^2\pi^2}{2m} \left[ \left( \frac{n_x}{L_x} \right)^2 + \left( \frac{n_y}{L_y} \right)^2 \right] \quad (21)

[edit] 3 dimensional box and higher

The same separation of variables technique can be applied to the three dimensional case to give the energy eigenfunctions:

\psi_{n_x,n_y,n_z} = \sqrt{\frac{8}{L_x L_y L_z}} \sin \left( \frac{n_x \pi x}{L_x} \right) \sin \left( \frac{n_y \pi y}{L_y} \right) \sin \left( \frac{n_z \pi z}{L_z} \right) \quad (22)
E_{n_x,n_y,n_z} = \frac{\hbar^2\pi^2}{2m} \left[ \left( \frac{n_x}{L_x} \right)^2 + \left( \frac{n_y}{L_y} \right)^2 + \left( \frac{n_z}{L_z} \right)^2 \right] \quad (23)

with

n_i=1,2,3,\ldots

An interesting feature of the above solutions is that when two or more of the lengths are the same (e.g. Lx = Ly), there are multiple wavefunctions corresponding to the same total energy. For example the wavefunction with nx = 2, ny = 1 has the same energy as the wavefunction with nx = 1, ny = 2. This situation is called degeneracy and for the case where exactly two degenerate wavefunctions have the same energy that energy level is said to be doubly degenerate. Degeneracy results from symmetry in the system. For the above case two of the lengths are equal so the system is symmetric with respect to a 90° rotation.

[edit] Free propagation

If the potential is zero (or constant) everywhere, one describes a free particle. This leads to some difficulties of normalization of the momentum or energy eigenfunctions. One way around is to constrain the particle in a finite volume V of arbitrary (large) extension, in which it is free to propagate. It is expected that in the limit of V→ ∞ we recover the free particle while allowing in the intermediate calculations the use of properly normalized states. Also, when describing for instance a particle propagating in a solid, one does not expect spatially localized states but instead completely delocalized states (within the solid), meaning that the particles propagates inside it (since it can be everywhere with the same probability, conversely to the sine solutions we encountered where the particle has favored locations). This understanding follows from the solutions of the Schrödinger equation for zero potential following from the so-called Born-von Karman boundary conditions; i.e., the wavefunction assumes same values on opposite sides of the box but it is not required to be zero here. One can then check that the following solutions obey eq. 1:

\textrm{in~1D}:\ \psi_k(x)={1\over\sqrt L}e^{ikx}; k={2n\pi\over L}; \quad n\in\mathbb{Z}
\textrm{in~3D}:\ \psi_{\mathbf{k}}(\mathbf{r})={1\over\sqrt{L^3}}e^{i\mathbf{k}\cdot\mathbf{r}}; k_x={2n_x\pi\over L}; k_y={2n_y\pi\over L}; k_z={2n_z\pi\over L}; \quad n_x, n_y, n_z\in\mathbb{Z}

The energy remains \hbar^2 k^2/2m (cf. eq. 3) but interestingly, now the k are twice as before (cf. eq. 7). This is because in the previous case, n was strictly positive whereas now it can be negative or zero (the ground state). The solutions where the sine does not superpose to itself after a translation of L can not be recovered with exponentials, since in this propagating particle interpretation, the derivative is discontinuous at the border, meaning that the particle acquires infinite velocity here. This shows how the two interpretations bear intrinsically differing behaviours.

[edit] References

  • Griffiths, David J. (2004). Introduction to Quantum Mechanics (2nd ed.). Prentice Hall. ISBN 0-13-111892-7. 

[edit] See also

[edit] External links

In other languages