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

${\displaystyle V(x)={\begin{cases}0,&0

where ${\displaystyle L}$ is the length of the box and ${\displaystyle 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 ${\displaystyle \psi (x,t)}$ can be found by solving the Schrödinger equation for the system

${\displaystyle 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),}$

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]

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

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

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

which is known as the dispersion relation for a free particle.[1]

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 ${\displaystyle 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

${\displaystyle \psi _{n}(x,t)={\begin{cases}A\sin(k_{n}x)\mathrm {e} ^{-i\omega _{n}t},&0

where ${\displaystyle n}$ is a positive integer. The wavenumber is restricted to certain, specific values given by[5]

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

where ${\displaystyle L}$ is the size of the box.[7] Negative values of ${\displaystyle n}$ are neglected, since they give wavefunctions identical to the positive ${\displaystyle n}$ solutions except for a physically unimportant sign change.[6]

Finally, the unknown constant ${\displaystyle 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

${\displaystyle \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 ${\displaystyle 0}$ and ${\displaystyle L}$. It's expected that the eigenvalues, i.e., the energy ${\displaystyle E_{n}}$ of the box should be the same regardless of its position in space, but ${\displaystyle \psi _{n}(x,t)}$ changes. This is represented by a more general case of ${\displaystyle \psi _{n}(x,t)}$:

${\displaystyle \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}

Where ${\displaystyle x_{0}}$ is the initial position. Notice that ${\displaystyle {\frac {n\pi x_{0}}{L}}}$ represents a phase shift in the wave function, and ${\displaystyle \psi _{n}(x,t)}$ simplifies to the above case when ${\displaystyle 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 ${\displaystyle k=p/\hbar }$ and ${\displaystyle \omega _{n}={\frac {\pi hn^{2}}{4L^{2}m}}}$,

${\displaystyle \phi _{n}(p,t)={\frac {1}{\sqrt {2\pi \hbar }}}\int _{-\infty }^{\infty }\psi _{n}(x,t)e^{-ikx}\,dx={\sqrt {\frac {\pi L}{\hbar }}}\,\,{\frac {n\left(1-(-1)^{n}e^{-ikL}\right)e^{-i\omega _{n}t}}{\pi ^{2}n^{2}-k^{2}L^{2}}}}$

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 ${\displaystyle 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

${\displaystyle P_{n}(x)={\begin{cases}{\frac {2}{L}}\sin ^{2}\left({\frac {n\pi x}{L}}\right);&0

Thus, for any value of n greater than one, there are regions within the box for which ${\displaystyle 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

${\displaystyle \langle x\rangle =\int _{-\infty }^{\infty }xP_{n}(x)\,{\mathrm {d} }x.}$

For the steady state particle in a box, it can be shown that the average position is always ${\displaystyle \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 ${\displaystyle cos(\omega t)}$.

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

${\displaystyle {\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 ${\displaystyle 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

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

where, again, ${\displaystyle 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:

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

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

${\displaystyle \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 ${\displaystyle \hbar }$ which obeys the Heisenberg uncertainty principle, which states that the product will be greater than or equal to ${\displaystyle \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]

${\displaystyle E_{n}={\frac {n^{2}\hbar ^{2}\pi ^{2}}{2mL^{2}}}={\frac {n^{2}h^{2}}{8mL^{2}}}}$.

The energy levels increase with ${\displaystyle 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]

${\displaystyle 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

${\displaystyle \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 ${\displaystyle 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 nx=4 and ny=4

If a particle is trapped in a two-dimensional box, it may freely move in the ${\displaystyle x}$ and ${\displaystyle y}$-directions, between barriers separated by lengths ${\displaystyle L_{x}}$ and ${\displaystyle 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

${\displaystyle \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)}$,
${\displaystyle E_{n_{x},n_{y}}={\frac {\hbar ^{2}k_{n_{x},n_{y}}^{2}}{2m}}}$,

where the two-dimensional wavevector is given by

${\displaystyle {\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

${\displaystyle \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)}$,
${\displaystyle 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

${\displaystyle {\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

${\displaystyle \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. ${\displaystyle L_{x}=L_{y}}$), there are multiple wavefunctions corresponding to the same total energy. For example the wavefunction with ${\displaystyle n_{x}=2,n_{y}=1}$ has the same energy as the wavefunction with ${\displaystyle 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

Template:Expand section 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

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

References

1. Davies, p.4
2. Actually, any constant, finite potential ${\displaystyle V_{0}}$ can be specified within the box. This merely shifts the energies of the states by ${\displaystyle V_{0}}$.
3. Davies, p. 1
4. Bransden and Joachain, p. 157
5. Davies p. 5
6. Bransden and Joachain, p.158
7. The simplest solutions, ${\displaystyle k_{n}=0}$ or ${\displaystyle A=0}$ both yield the trivial wavefunction ${\displaystyle \psi (x)=0}$, which describes a particle that does not exist anywhere in the system.[6]
8. Bransden and Joachain, p. 159
9. Davies, p. 15
10. {{#invoke:Citation/CS1|citation |CitationClass=journal }}

Bibliography

• {{#invoke:citation/CS1|citation

|CitationClass=book }}

• {{#invoke:citation/CS1|citation

|CitationClass=book }}

• {{#invoke:citation/CS1|citation

|CitationClass=book }}