Particle in a box

Some trajectories of a particle in a box according to Newton's laws of classical mechanics (A), and according to the Schrödinger equation of quantum mechanics (B-F). In (B-F), the horizontal axis is position, and the vertical axis is the real part (blue) and imaginary part (red) of the wavefunction. The states (B,C,D) are energy eigenstates, but (E,F) are not.

In quantum mechanics, the particle in a box model (also known as the infinite potential well or the infinite square well) describes a particle free to move in a small space surrounded by impenetrable barriers. The model is mainly used as a hypothetical example to illustrate the differences between classical and quantum systems. In classical systems, for example a ball trapped inside a large box, the particle can move at any speed within the box and it is no more likely to be found at one position than another. However, when the well becomes very narrow (on the scale of a few nanometers), quantum effects become important. The particle may only occupy certain positive energy levels. Likewise, it can never have zero energy, meaning that the particle can never "sit still". Additionally, it is more likely to be found at certain positions than at others, depending on its energy level. The particle may never be detected at certain positions, known as spatial nodes.

The particle in a box model provides one of the very few problems in quantum mechanics which can be solved analytically, without approximations. This means that the observable properties of the particle (such as its energy and position) are related to the mass of the particle and the width of the well by simple mathematical expressions. Due to its simplicity, the model allows insight into quantum effects without the need for complicated mathematics. It is one of the first quantum mechanics problems taught in undergraduate physics courses, and it is commonly used as an approximation for more complicated quantum systems. See also: the history of quantum mechanics.

One-dimensional solution

The barriers outside a one-dimensional box have infinitely large potential, while the interior of the box has a constant, zero potential.

The simplest form of the particle in a box model considers a one-dimensional system. Here, the particle may only move backwards and forwards along a straight line with impenetrable barriers at either end.^[1] The walls of a one-dimensional box may be visualised as regions of space with an infinitely large potential energy. Conversely, the interior of the box has a constant, zero potential energy.^[2] This means that no forces act upon the particle inside the box and it can move freely in that region. However, infinitely large forces repel the particle if it touches the walls of the box, preventing it from escaping. The potential energy in this model is given as

V(x) = \begin{cases} 0, & 0 < x < L,\\ \infty, & \text{otherwise,} \end{cases},

where $L$ is the length of the box and $x$ is the position of the particle within the box.

Wavefunctions

In quantum mechanics, the wavefunction gives the most fundamental description of the behavior of a particle; the measurable properties of the particle (such as its position, momentum and energy) may all be derived from the wavefunction.^[3] The wavefunction $\psi(x,t)$ can be found by solving the Schrödinger equation for the system

i\hbar\frac{\partial}{\partial t}\psi(x,t) = -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2}\psi(x,t) +V(x)\psi(x,t),

where $\hbar$ is the reduced Planck constant, $m$ is the mass of the particle, $i$ is the imaginary unit and $t$ is time.

Inside the box, no forces act upon the particle, which means that the part of the wavefunction inside the box oscillates through space and time with the same form as a free particle:^[1]^[4]

$\psi(x,t) = [A \sin(kx) + B \cos(kx)]\mathrm{e}^{-i\omega t},$

(1)

where $A$ and $B$ are arbitrary complex numbers. The frequency of the oscillations through space and time are given by the wavenumber $k$ and the angular frequency $\omega$ respectively. These are both related to the total energy of the particle by the expression

E = \hbar\omega = \frac{\hbar^2 k^2}{2m},

which is known as the dispersion relation for a free particle.^[1] Here one must notice that now, since the particle is not entirely free but under the influence of a potential (the potential V described above), the energy of the particle given above is not the same thing as $\frac{p^2}{2m}$ where p is the momentum of the particle, and thus the wavenumber k above actually describes the energy states of the particle, not the momentum states (id est, it turns out that the momentum of the particle is not given by $p = \hbar k$ ). In this sense, it is quite dangerous to call the number k a wavenumber, since it is not related to momentum like "wavenumber" usually is. The rationale for calling k the wavenumber is that it enumerates the number of crests that the wavefuntion has inside the box, and in this sense it is a wavenumber. This indiscrepancy can be seen more clearly below, when we find out that the energy spectrum of the particle is discrete (only discrete values of energy is allowed) but the momentum spectrum is continuous (momentum can vary continuously) and in particular, the relation $E = \frac{p^2}{2m}$ for the energy and momentum of the particle does not hold. As said above, the reason this relation between energy and momentum does not hold is that the particle is not free, but there is a potential V in the system, and the energy of the particle is $E = T + V$ , where T is the kinetic and V the potential energy.

Initial wavefunctions for the first four states in a one-dimensional particle in a box

The size (or amplitude) of the wavefunction at a given position is related to the probability of finding a particle there by $P(x,t) = |\psi(x,t)|^2$ . The wavefunction must therefore vanish everywhere beyond the edges of the box.^[1]^[4] Also, the amplitude of the wavefunction may not "jump" abruptly from one point to the next.^[1] These two conditions are only satisfied by wavefunctions with the form

\psi_n(x,t) = \begin{cases} A \sin(k_n x)\mathrm{e}^{-i\omega_n t}, & 0 < x < L,\\ 0, & \text{otherwise,} \end{cases}

where $n$ is a positive integer. Usually in quantum mechanics it is also demanded that the derivative of the wavefunction in addition to the wavefunction itself be continuous; here this demand would lead to the only solution being the constant zero function, which is not what we desire, so we give up this demand (as this system with infinite potential can be regarded as a nonphysical abstract limiting case, we can treat it as such and "bend the rules"). Note that giving up this demand means that the wavefuntion is not differentiable function at the boundary of the box, and thus it can be said that the wavefunction does not solve the Schrödinger equation at the boundary points $x = 0$ and $x = L$ (but does solve everywhere else). The wavenumber above is restricted to certain, specific values given by^[5]

k_n = \frac{n \pi}{L}, \quad \mathrm{where} \quad n = \{1,2,3,4,\ldots\},

where $L$ is the size of the box.^[6] Negative values of $n$ are neglected, since they give wavefunctions identical to the positive $n$ solutions except for a physically unimportant sign change.^[7] Here one sees that only a discrete set of energy values and wavenumbers k are allowed for the particle.

Finally, the unknown constant $A$ may be found by normalizing the wavefunction so that the total probability density of finding the particle in the system is 1. It follows that

\left| A \right| = \sqrt{\frac{2 }{L}}.

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

The above solution is for the specific case of a box situated between $0$ and $L$ . It's expected that the eigenvalues, i.e., the energy $E_n$ of the box should be the same regardless of its position in space, but $\psi_n(x,t)$ changes. This is represented by a more general case of $\psi_n(x,t)$ :

\psi_n(x,t) = \begin{cases} \sqrt{\frac{2 }{L}} \sin(k_n x - \frac{n \pi {x_0}}{L})\mathrm{e}^{-i\omega_n t}, & x_0 < x < x_0 + L,\\ 0, & \text{otherwise,} \end{cases}

Where $x_0$ is the initial position. Notice that $\frac{n \pi x_0}{L}$ represents a phase shift in the wave function, and $\psi_n(x,t)$ simplifies to the above case when $x_0 = 0$ . Also, the phase shift has no effect when solving the Schrödinger equation, thus not affecting the eigenvalue.

The momentum wavefunction is proportional to the Fourier transform of the position wavefunction. With $k'=p/\hbar$ (note that this $k'$ , describing the momentum states, is not the same thing as the $k$ above, which described the energy states) and $\omega_n=\frac{\pi h n^2}{4 L^2 m}$ ,

\phi_n(p,t)=\frac{1}{\sqrt{2\pi\hbar}}\int_{-\infty}^\infty \psi_n(x,t)e^{-ik'x}\,dx = \sqrt{\frac{\pi L}{\hbar}}\,\,\frac{n\left(1-(-1)^ne^{-ik'L}\right) e^{-i \omega_n t}}{\pi ^2 n^2-(k')^2 L^2}

Here one sees that the momentum spectrum is continuous, and one can conclude that for energy state described by the $k_n$ wavenumber, the momentum can, when measured, attain also other values than $p = \pm \hbar k_n$ , although these are the most likely momentum values (the momentum wavefunction attains its greatest values at these momentum values). Hence one also sees that since the energy is $E_n = \frac{\hbar^2 k_n^2}{2m}$ (since we are talking about the n:th energy eigenstate) the relation $E = \frac{p^2}{2m}$ does not necessarily hold; here p is the measured momentum (since the energy eigenstate $\psi_n$ is not a momentum eigenstate, and actually not even a superposition of two momentum eigenstates, as one could imagine from equation (1) above, it has no well-defined momentum before measurement).

Position and momentum

In classical physics, the particle can be detected anywhere in the box with equal probability. In quantum mechanics, however, the probability density for finding a particle at a given position is derived from the wavefunction as $P(x) = |\psi(x)|^2.$ For the particle in a box, the probability density for finding the particle at a given position depends upon its state, and is given by

P_n(x) = \begin{cases} \frac{2 }{L}\sin^2\left(\frac{n\pi x}{L}\right); & 0 < x < L \\ 0; & \text{otherwise}. \end{cases}

Thus, for any value of n greater than one, there are regions within the box for which $P(x)=0$ , indicating that spatial nodes exist at which the particle cannot be found.

In quantum mechanics, the average, or expectation value of the position of a particle is given by

\langle x \rangle = \int_{-\infty}^{\infty} x P_n(x)\,\mathrm{d}x.

For the steady state particle in a box, it can be shown that the average position is always $\langle x \rangle = L/2$ , regardless of the state of the particle. For a superposition of states, the expectation value of the position will change based on the cross term which is proportional to $cos(\omega t)$ .

The variance in the position is a measure of the uncertainty in position of the particle:

\mathrm{Var}(x)=\int_{-\infty}^\infty (x-\langle x\rangle)^2 P_n(x)\,dx = \frac{L^2}{12}\left(1-\frac{6}{n^2\pi^2}\right)

The probability density for finding a particle with a given momentum is derived from the wavefunction as $P(x) = |\phi(x)|^2.$ . As with position, the probability density for finding the particle at a given momentum depends upon its state, and is given by

P_n(p)=\frac{2 \pi L}{\hbar}\, \frac{n^2 \left(1-(-1)^n \cos (k L)\right)}{\left(k^2 L^2-\pi ^2 n^2\right)^2}

where, again, $k=p/\hbar$ . The expectation value for the momentum is then calculated to be zero, and the variance in the momentum is calculated to be:

\mathrm{Var}(p)=\left(\frac{\hbar n\pi}{L}\right)^2

The uncertainties in position and momentum ( $\Delta x$ and $\Delta p$ ) are defined as being equal to the square root of their respective variances, so that:

\Delta x \Delta p = \frac{\hbar}{2} \sqrt{\frac{n^2\pi^2}{3}-2}

This product increases with increasing n, having a minimum value for n=1. The value of this product for n=1 is about equal to 0.568 $\hbar$ which obeys the Heisenberg uncertainty principle, which states that the product will be greater than or equal to $\hbar/2$

Energy levels

The energy of a particle in a box (black circles) and a free particle (grey line) both depend upon wavenumber in the same way. However, the particle in a box may only have certain, discrete energy levels.

The energies which correspond with each of the permitted wavenumbers may be written as^[5]

E_n = \frac{n^2\hbar^2 \pi ^2}{2mL^2} = \frac{n^2 h^2}{8mL^2}

The energy levels increase with $n^2$ , meaning that high energy levels are separated from each other by a greater amount than low energy levels are. The lowest possible energy for the particle (its zero-point energy) is found in state 1, which is given by^[8]

E_1 = \frac{\hbar^2\pi^2}{2mL^2}.

The particle, therefore, always has a positive energy. This contrasts with classical systems, where the particle can have zero energy by resting motionlessly. This can be explained in terms of the uncertainty principle, which states that the product of the uncertainties in the position and momentum of a particle is limited by

\Delta x\Delta p \geq \frac{\hbar}{2}

It can be shown that the uncertainty in the position of the particle is proportional to the width of the box.^[9] Thus, the uncertainty in momentum is roughly inversely proportional to the width of the box.^[8] The kinetic energy of a particle is given by $E=p^2/(2m)$ , and hence the minimum kinetic energy of the particle in a box is inversely proportional to the mass and the square of the well width, in qualitative agreement with the calculation above.^[8]

Higher-dimensional boxes

The wavefunction of a 2D well with n_x=4 and n_y=4

If a particle is trapped in a two-dimensional box, it may freely move in the $x$ and $y$ -directions, between barriers separated by lengths $L_x$ and $L_y$ respectively. Using a similar approach to that of the one-dimensional box, it can be shown that the wavefunctions and energies are given respectively by

\psi_{n_x,n_y} = \sqrt{\frac{4}{L_x L_y}} \sin \left( k_{n_x} x \right) \sin \left( k_{n_y} y\right)

E_{n_x,n_y} = \frac{\hbar^2 k_{n_x,n_y}^2}{2m}

where the two-dimensional wavevector is given by

\mathbf{k_{n_x,n_y}} = k_{n_x}\mathbf{\hat{x}} + k_{n_y}\mathbf{\hat{y}} = \frac{n_x \pi }{L_x} \mathbf{\hat{x}} + \frac{n_y \pi }{L_y} \mathbf{\hat{y}}

For a three dimensional box, the solutions are

\psi_{n_x,n_y,n_z} = \sqrt{\frac{8}{L_x L_y L_z}} \sin \left( k_{n_x} x \right) \sin \left( k_{n_y} y \right) \sin \left( k_{n_z} z \right)

E_{n_x,n_y,n_z} = \frac{\hbar^2 k_{n_x,n_y,n_z}^2}{2m}

where the three-dimensional wavevector is given by

\mathbf{k_{n_x,n_y,n_z}} = k_{n_x}\mathbf{\hat{x}} + k_{n_y}\mathbf{\hat{y}} + k_{n_z}\mathbf{\hat{z}} = \frac{n_x \pi }{L_x} \mathbf{\hat{x}} + \frac{n_y \pi }{L_y} \mathbf{\hat{y}} + \frac{n_z \pi }{L_z} \mathbf{\hat{z}}

In general for an n-dimensional box, the solutions are

\psi = \sqrt{\frac{2^n}{\prod_{i} L_i}} \prod_{i}\sin(k_i i)

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.

Applications

Because of its mathematical simplicity, the particle in a box model is used to find approximate solutions for more complex physical systems in which a particle is trapped in a narrow region of low electric potential between two high potential barriers. These quantum well systems are particularly important in optoelectronics, and are used in devices such as the quantum well laser, the quantum well infrared photodetector and the quantum-confined Stark effect modulator. It is also used to model a lattice in the Kronig-Penny model and for a finite metal with the free electron approximation.

Relativistic Effects

The probability density does not go to zero at the nodes if relativistic effects are taken into account.^[10]

References

↑ 1.0 1.1 1.2 1.3 1.4 Davies, p.4
↑ Actually, any constant, finite potential $V_0$ can be specified within the box. This merely shifts the energies of the states by $V_0$ .
↑ Davies, p. 1
↑ 4.0 4.1 Bransden and Joachain, p. 157
↑ 5.0 5.1 Davies p. 5
↑ The simplest solutions, $k_n=0$ or $A=0$ both yield the trivial wavefunction $\psi(x)=0$ , which describes a particle that does not exist anywhere in the system.^[7]
↑ 7.0 7.1 Bransden and Joachain, p.158
↑ 8.0 8.1 8.2 Bransden and Joachain, p. 159
↑ Davies, p. 15
↑ Alberto, P; Fiolhais, C; Gil, V M S (1996). "Relativistic particle in a box" (PDF). European Journal of Physics 17: 19–24. Bibcode:1996EJPh...17...19A. doi:10.1088/0143-0807/17/1/004.

Bibliography

Bransden, B. H.; Joachain, C. J. (2000). Quantum mechanics (2nd ed.). Essex: Pearson Education. ISBN 0-582-35691-1.
Davies, John H. (2006). The Physics of Low-Dimensional Semiconductors: An Introduction (6th reprint ed.). Cambridge University Press. ISBN 0-521-48491-X.
Griffiths, David J. (2004). Introduction to Quantum Mechanics (2nd ed.). Prentice Hall. ISBN 0-13-111892-7.

External links

Scienceworld (Infinite Potential Well)
1-D quantum mechanics java applet simulates particle in a box, as well as other 1-dimensional cases.
2-D particle in a box applet