Particle in a box

In physics, the particle in a box (also known as the infinite potential well or the infinite square well) is a problem consisting of a single particle inside of an infinitely deep potential well, 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.

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. The particle's confinement leads to a quantization of its allowed energy eigenvalues, where a vanishing energy eigenvalue is not allowed. In other words, the least energy state of the system has a non-zero energy. Furthermore, the probability of measuring particle in the box at a given energy level is non-uniform and probabilistic -- upon measurement, certain locations in the "box" where the particle may be found have a higher probability than others of being the location where the particle is actually found when a measurement is made. There are also places within the box where the particle cannot ever be found. Both of these results differ from the classical description of classical mechanics and everyday experience, but yet the results rest on principles that have been extensively experimentally verified.

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 particle 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 boundary 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:

# 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).
#Zero-point energy - The lowest possible energy level of the particle, called the zero-point energy, is nonzero.
#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).

1-dimensional box

olution

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\}\; ight)\}\; ,$::$E\_n\; =\; frac\{hbar^2\; pi^2\}\{2\; m\; L^2\}\; n^2\; ,$

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::$psileft(x\; ight),$ is the complex-valued stationary time-independent wavefunction that we want to find::$Vleft(x\; ight),$ 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:

::$-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 hand-waving 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")|², must be small in these regions and decreases with increasing potential. For the case of an infinite potential, |"ψ"("x")|² 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$ ($psi\; =\; 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 $psi\; =\; 0$ everywhere (i.e. the particle isn't in the box.) and can be thrown out. If $A\; eq\; 0$ then $sin(kL)\; =\; 0$ only when:

::$k\; =\; frac\{n\; pi\}\{L\}\; quad\; mathrm\{where\}\; quad\; nin\; mathbb\{Z\}^+\; quad\; (7)$

Note that $n\; =\; 0$ is ruled out because in this case, $psi=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)$, giving the same probability when squared. Now in order to find $A$ we must normalize the wavefunction. We recognize that the particle must exist somewhere in space. $|psi(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)\; ight|^2\; ,\; mathrm\{d\}x\; =\; left|\; A\; ight|^2\; int\_0^L\; sin^2\; left(kx\; ight)\; ,\; mathrm\{d\}x\; =\; left|\; A\; ight|^2\; frac\{L\}\{2\}$:or::$left|\; A\; ight|\; =\; 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.

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

:$psi\_n(x)\; =\; sqrt\{frac\{2\}\{L\; sinleft(frac\{n\; pi\; x\}\{L\}\; ight)\; quad\; (9)$

:$E\_n\; =\; frac\{n^2hbar^2\; pi\; ^2\}\{2mL^2\}\; =\; frac\{n^2\; h^2\}\{8mL^2\}\; quad\; (10)$with $n$ being a positive integer.

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")² both equal zero, indicating that for these energy levels, nodes exist in the box where the probability of finding the particle is zero. A distance between the particle and any point of the box or a distance between two particles in the box may be calculated using the probability metric.

Heuristic derivation of eigenenergies

The formula for the energy eigenvalues may be derived heuristically by considering that the wavefunction must have nodes at the boundary of the box. This is, in other words, a restatement of the continuity condition of the wavefunction at the boundaries of the box, which implies that a half-integer number of wavelengths must fit into the box::$nfrac\{lambda\}\{2\}\; =\; L$

Using the de Broglie relation, we may write the momentum of the particle as:$p\; =\; frac\{h\}\{lambda\}\; =\; frac\{h\}\{2L\}n.$

The classical formula relating energy to momentum may then be used to give the energy eigenvalues of the system.:$E\; =\; frac\{p^2\}\{2m\}\; =\; frac\{n^2\; h^2\}\{8mL^2\}$

2-dimensional box

olution

::$psi\_\{n\_x,n\_y\}\; =\; sqrt\{frac\{4\}\{L\_x\; L\_y\; sin\; left(\; frac\{n\_x\; pi\; x\}\{L\_x\}\; ight)\; sin\; left(\; frac\{n\_y\; pi\; y\}\{L\_y\}\; ight)$

::$E\_\{n\_x,n\_y\}\; =\; frac\{hbar^2pi^2\}\{2m\}\; left\; [\; left(\; frac\{n\_x\}\{L\_x\}\; ight)^2\; +\; left(\; frac\{n\_y\}\{L\_y\}\; ight)^2\; ight]$

Derivation

For the 2-dimensional case the particle is confined to a rectangular surface of length "L_x" in the "x"-direction and "L_y" 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^2psi\}\{partial\; x^2\}+frac\{partial^2\; psi\}\{partial\; y^2\}\; ight)\; =Epsi\; 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(\; Yfrac\{partial^2X\}\{partial\; x^2\}+Xfrac\{partial^2\; Y\}\{partial\; y^2\}\; ight)\; =E\; X\; Y\; quad(13)$

which upon dividing by XY and rewriting d²"X"/d"x"² as "X" and d²"Y"/d"y"² as "Y" becomes:

:$-frac\{hbar^2\}\{2m\}\; left(\; frac\{X"\}\{X\}+frac\{Y"\}\{Y\}\; ight)\; =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 "E"_x + "E"_y = "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\; Xquad\; (16)$:$-frac\{hbar^2\}\{2m\}frac\{partial^2Y\}\{partial\; y^2\}\; =\; E\_y\; Yquad\; (17)$

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\}\; ight)\; quad\; (18)$:$Y\_\{n\_y\}=sqrt\{frac\{2\}\{L\_y\; sin\; left(\; frac\{n\_y\; pi\; y\}\{L\_y\}\; ight)\; quad\; (19)$

Finally, since "ψ"="XY" and "E" = "E"_x + "E"_y, 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\}\; ight)\; sin\; left(\; frac\{n\_y\; pi\; y\}\{L\_y\}\; ight)\; quad\; (20)$

:$E\_\{n\_x,n\_y\}\; =\; frac\{hbar^2pi^2\}\{2m\}\; left\; [\; left(\; frac\{n\_x\}\{L\_x\}\; ight)^2\; +\; left(\; frac\{n\_y\}\{L\_y\}\; ight)^2\; ight]\; quad\; (21)$

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\}\; ight)\; sin\; left(\; frac\{n\_y\; pi\; y\}\{L\_y\}\; ight)\; sin\; left(\; frac\{n\_z\; pi\; z\}\{L\_z\}\; ight)\; quad\; (22)$

:$E\_\{n\_x,n\_y,n\_z\}\; =\; frac\{hbar^2pi^2\}\{2m\}\; left\; [\; left(\; frac\{n\_x\}\{L\_x\}\; ight)^2\; +\; left(\; frac\{n\_y\}\{L\_y\}\; ight)^2\; +\; left(\; frac\{n\_z\}\{L\_z\}\; ight)^2\; ight]\; 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. "L"_"x" = "L"_"y"), there are multiple wavefunctions corresponding to the same total energy. For example the wavefunction with "n"_"x" = 2, "n"_"y" = 1 has the same energy as the wavefunction with "n"_"x" = 1, "n"_"y" = 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.

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:

:$extrm\{in~1D\}:\; psi\_k(x)=\{1oversqrt\; L\}e^\{ikx\};\; k=\{2npiover\; L\};\; quad\; ninmathbb\{Z\}$:$extrm\{in~3D\}:\; psi\_\{mathbf\{k(mathbf\{r\})=\{1oversqrt\{L^3e^\{imathbf\{k\}cdotmathbf\{r;\; k\_x=\{2n\_xpiover\; L\};\; k\_y=\{2n\_ypiover\; L\};\; k\_z=\{2n\_zpiover\; L\};\; quad\; n\_x,\; n\_y,\; n\_zinmathbb\{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.

References

*cite book | author=Griffiths, David J.|title=Introduction to Quantum Mechanics (2nd ed.) | publisher=Prentice Hall |year=2004 |id=ISBN 0-13-111892-7

ee also

* finite potential well
* Delta function potential
* gas in a box
* particle in a ring
* Particle in a spherically symmetric potential
* Delta potential well (QM)
* Semicircle potential well
* Configuration integral (statistical mechanics)

External links

* [http://scienceworld.wolfram.com/physics/InfiniteSquarePotentialWell.html Scienceworld] (Infinite Potential Well)
* [http://scienceworld.wolfram.com/physics/FiniteSquarePotentialWell.html Scienceworld] (Finite Potential Well)
* [http://www.falstad.com/qm1d/ 1-D quantum mechanics java applet] simulates particle in a box, as well as other 1-dimensional cases.
* [http://www.falstad.com/qm2dbox/ 2-D particle in a box applet]