Revision 4025068 of "Particle in a box" on simplewiki

[[File:Particle in a box.svg|thumb|300px| A classical particle, such as a ball, bouncing between two high walls. The particle moves forwards and backwards, and can not go out.]]

[[File:particle in a box wavefunctions 2.svg|thumb|right|upright|A '''particle in a box''' is free to move in a space surrounded by impenetrable barriers). When the barriers lie very close together, quantum effects are observed.  For example, the particle is more likely to be found at certain positions than others and it may only occupy specific energy levels.]]

In [[quantum mechanics]], the '''particle in a box''' model (also known as the '''infinite potential well''' or the '''infinite square well''') describes a particle that can move freely between two walls.  This is used to show the differences between [[classical physics|classical]] and quantum systems.  In classical systems, for example a ball trapped inside a heavy 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 box becomes very narrow (on the scale of a few [[nanometer]]s), quantum effects become important.  The particle may only have certain [[energy levels]].  Likewise, it can never have zero energy, meaning that the particle can never "sit still". It is more likely to be found at certain places than at others, depending on its energy level.  The particle may never be found at certain places, known as spatial nodes.

The particle in a box model provides one of the very few problems in quantum mechanics which can be solved without approximations.  This means that the energy and position of the particle 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 ==
[[File:Infinite potential well-en.svg|thumb|right|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.<ref name="Davies4">Davies, p.4</ref>
The walls of a one-dimensional box may be visualised as regions of space with an infinitely large [[potential]].  Conversely, the interior of the box has a constant, zero potential.<ref>Actually, any constant, finite potential <math>V_0</math> can be specified within the box.  This merely shifts the energies of the states by <math>V_0</math>.</ref>  This means that no forces act upon the particle inside the box and it can move freely in that region.  However, infinitely large [[Force (physics)|forces]] repel the particle if it touches the walls of the box, preventing it from escaping.  The potential in this model is given as
:<math>V(x) =
\begin{cases}
0, & 0 < x < L,\\
\infty, & \text{otherwise,}
\end{cases},
</math>
where <math>L</math> is the length of the box and <math>x</math> 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.<ref name="Davies1">Davies, p. 1</ref>
The wavefunction <math>\psi(x,t)</math> can be found by solving the [[Schrödinger equation]] for the system
:<math>\mathrm{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)</math>

where <math>\hbar</math> is the [[reduced Planck constant]], <math>m</math> is the [[mass]] of the particle, <math>\mathrm{i}</math> is the [[imaginary unit]] and <math>t</math> is time.


To solve for <math>\Psi(x,t)\;</math> the first thing we do is assume that <math>\Psi(x,t)\;</math> is the product of two functions, <math>\Psi(x,t)=\psi(x)\phi(t)\;</math>

Substituting into <math>\mathrm{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)</math> gives us:
:<math>\mathrm{i}\hbar\frac{\partial}{\partial t}\psi(x)\phi(t)=-\frac{\hbar^{2}}{2m}\frac{\partial^{2}}{\partial x^{2}}\psi(x)\phi(t)+V(x)\psi(x)\phi(t)</math>

Since the <math>\phi(t)</math> is a constant relative to x, and <math>\psi(x)</math> is constant relative to t, we can pull out the 'constants' out of the partial derivatives:
:<math>\mathrm{i}\hbar\psi(x)\frac{\partial}{\partial t}\phi(t)=-\frac{\hbar^{2}}{2m}\phi(t)\frac{\partial^{2}}{\partial x^{2}}\psi(x)+V(x)\psi(x)\phi(t)</math>

We can then divide the entire equation by <math>\psi(x)\phi(t)</math>
:<math>\mathrm{i}\hbar\frac{\frac{\partial}{\partial t}\phi(t)}{\phi(t)}=-\frac{\hbar^{2}}{2m}\frac{\frac{\partial^{2}}{\partial x^{2}}\psi(x)}{\psi(x)}+V(x)</math>

Now, the entire equation on the left is dependent on t, and the entire equation on the right is dependent on x. However, the only way this could be true would be if both sides are a constant. We can call this constant C. So now we have two equations:
:<math>\mathrm{i}\hbar\frac{\frac{\partial}{\partial t}\phi(t)}{\phi(t)}=C</math>
and
:<math>-\frac{\hbar^{2}}{2m}\frac{\frac{\partial^{2}}{\partial x^{2}}\psi(x)}{\psi(x)}+V(x)=C</math>

==== Solving the Time-Dependent Portion ====
Now we'll focus on solving for <math>\phi(t)</math> the time dependent part of the equation:
:<math>\mathrm{i}\hbar\frac{\frac{\partial}{\partial t}\phi(t)}{\phi(t)}=C</math>

Some algebraic manipulation gives us:
:<math>\frac{\partial\phi(t)}{\phi(t)}=\frac{-\mathrm{i}C}{\hbar}\partial t</math> (because <math>(\frac{1}{i}=-i)</math>)

Solving this differential equation yields:
:<math>ln[\phi(t)]=\frac{-\mathrm{i}C}{\hbar}t</math>

We then raise e to the powers of both sides to isolate <math>\phi(t)</math>:
:<math>\phi(t)=e^{\frac{-\mathrm{i}C}{\hbar}t}</math>

We can then expand <math>\phi(t)</math> with Eulers equation, and we get:
:<math>\phi(t)=\cos(\frac{C}{\hbar}t)-i\sin(\frac{C}{\hbar}t)</math>

So we see that it's a wave with a frequency <math>\omega=\frac{C}{\hbar}</math> . Therefore, <math>C=\hbar\omega</math>

However, we can take it one step further, and say that in this case, <math>E=\hbar\omega</math>

Substituting back into <math>\phi(t)=e^{\frac{-\mathrm{i}C}{\hbar}t}</math> gives us our final solution for the time dependent part of the function:
:<math>\phi(t)=e^{-\mathrm{i(}\frac{E}{\hbar})t}</math>

==== Solving the Spatial Portion ====
Continuing from before, we have:
:<math>-\frac{\hbar^{2}}{2m}\frac{\frac{\partial^{2}}{\partial x^{2}}\psi(x)}{\psi(x)}+V(x)=E</math>

As mentioned earlier, the definition of the 'particle in a box' scenario is that there is an infinite potential outside of the box, and no potential within the box. Therefore:
:<math>V(x) =
\begin{cases}
0, & 0 < x < L,\\
\infty, & \text{otherwise,}
\end{cases},
</math>

Substituting <math>V(x)=0\;</math> within the box gives us:
:<math>-\frac{\hbar^{2}}{2m}\frac{\frac{\partial^{2}}{\partial x^{2}}\psi(x)}{\psi(x)}=E</math>

Doing some algebra gives us:
:<math>\frac{\partial^{2}}{\partial x^{2}}\psi(x)+\frac{2mE}{\hbar^{2}}\psi(x)=0</math>
The solutions for a differential equation of the form <math>ay''+by'+cy=0\;</math> where <math>b^{2}-4ac<0\;</math> are given to be
:<math>y(t)=c_{1}e^{\lambda t}cos(\mu t)+c_{2}e^{\lambda t}sin(\mu t)\;</math>, where <math>\lambda=\frac{-b}{2a}\;</math>, and <math>\mu=\frac{\sqrt{b^{2}-4ac}}{2a*i}\;</math>
In our case, <math>b=0\;</math>, and <math>c = \frac{2m}{\hbar^{2}}E\;</math>, so <math>\lambda=0\;</math> and <math>\mu=\sqrt{\frac{2m}{\hbar^{2}}E}\;</math>. For ease in writing, we'll set <math>k=\mu=\sqrt{\frac{2m}{\hbar^{2}}E}\;</math>, so our solution is <math>y(t)=c_{1}cos(kt)+c_{2}sin(kt)\;</math>

Therefore, we can say the form of <math>\psi(x)\;</math> is:
<math>\psi(x) = Asin(kx) + Bcos(kx)\;</math>

The wave function needs to be continuous, and because it is 0 outside of the walls of the well, it needs to be 0 at the point <math>x = 0\;</math>. Therefore <math>B=0\;</math>, and <math>\psi(x) = Asin(kx)\;</math>

On the same note, the wave function needs to be continuous, and as it is 0 outside of the walls of the well, it needs to be 0 at the point <math>x = L\;</math>. In order for that to be true, the sin must be 0. <math>\sin\;</math> is 0 at integer values of <math>\pi\;</math>. Therefore <math>kL=n\pi\;</math>, and <math>k=\frac{n\pi}{L}\;</math>where <math>n\;</math> is an integer.

Therefore, <math>\psi(x) = Asin(\frac{n\pi x}{L})\;</math>

One of the conditions of the wave function is that the probability of finding it somewhere the enclosed well must be a total of 1 (100%).
The size (or [[amplitude]]) of the wavefunction at a given position is related to the probability of finding a particle there by <math>P(x) = |\psi(x)|^2\;</math>.
Therefore, <math>\int_{0}^{L}|\psi(x)|^{2}=1\;</math>, or
:<math>\int_{0}^{L}A^{2}\sin^{2}(kx)dx=1\;</math>. 
:<math>A^{2}\int_{0}^{L}\frac{1-\cos(2kx)}{2}dx=1\;</math>
:<math>\frac{A^{2}}{2}\int_{0}^{L}dx-\frac{A^{2}}{2}\int_{0}^{L}\cos(2kx)dx=1\;</math>

However, integrating the cos portion gives us:<math>\int_{0}^{L}\cos(2kx)=\frac{L\sin(2kL)}{2k}\;</math>. Substituting in <math>k=\frac{n\pi}{L}\;</math> gives us <math>\frac{L\sin(2\frac{n\pi}{L}L)}{2\frac{n\pi}{L}}=\frac{L\sin(2n\pi)}{2\frac{n\pi x}{L}}=0\;</math>, because <math>\sin(2n\pi)\;</math>, where n is an integer is always 0.

Therefore, the only part remaining is the <math>\frac{A^{2}}{2}\int_{0}^{L}dx\;</math>, which gives us <math>\frac{A^{2}L}{2}=1\;</math>, and after some algebraic manipulation, we reach our result of:
:<math>A=\sqrt{\frac{2}{L}}\;</math>

Now we have our solution for the spatial component:
<math>\psi(x) = \sqrt{\frac{2}{L}}sin(\frac{n\pi x}{L})\;</math>, where n is an integer.

[[File:particle in a box wavefunctions 2.svg|thumb|right|upright|Initial wavefunctions for the first four states in a one-dimensional particle in a box]]

==== Putting them together ====
The first thing we did was assume that <math>\Psi(x,t)\;</math> is the product of two functions, <math>\Psi(x,t)=\psi(x)\phi(t)\;</math>

:<math>\phi(t)=e^{-\mathrm{i(}\frac{E}{\hbar})t}</math>
:<math>\psi(x) = \sqrt{\frac{2}{L}}sin(\frac{n\pi x}{L})\;</math>, where n is an integer.

Therefore, our result for the time dependent Schrödinger equation is:
:<math>\Psi(x,t)=\psi(x)\phi(t)=\sqrt{\frac{2}{L}}sin(\frac{n\pi x}{L})e^{-\mathrm{i(}\frac{E}{\hbar})t}\;</math>

=== Energy levels ===
[[File:Confined particle dispersion - positive.svg|thumb|upright|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.]]

Because <math>k=\sqrt{\frac{2m}{\hbar^{2}}E}\;</math> and <math>kL=n\pi\;</math>, we can combine the two: <math>\sqrt{\frac{2m}{\hbar^{2}}E}L=n\pi\;</math>

After doing some algebra, we can solve for E (Energy) to be:
:<math>E=\frac{n^{2}\pi^{2}\hbar^{2}}{2mL^{2}}\;</math>

We can get an alternate form by substituting <math>\hbar=\frac{h}{2\pi}\;</math>
:<math>E_n = \frac{n^2 h^2}{8mL^2}</math>.

The energy levels increase with <math>n^2</math>, 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 (when <math>n=1</math>), which is given by<ref name="Bransden159">Bransden and Joachain, p. 159</ref>
:<math>E_1 = \frac{\hbar^2\pi^2}{2mL^2}.</math>
The particle, therefore, always has a positive energy.  This contrasts with classical systems, where the particle can have zero energy by resting motionless at the bottom of the box.  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
:<math>\Delta x\Delta p \geq \frac{\hbar}{2}</math>
It can be shown that the uncertainty in the position of the particle is proportional to the width of the box.<ref name="Davies15">Davies, p. 15</ref>  Thus, the uncertainty in momentum is roughly inversely proportional to the width of the box.<ref name="Bransden159" />  The kinetic energy of a particle is given by <math>E=p^2/(2m)</math>, 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.<ref name="Bransden159" />

=== Spatial location ===
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 <math>P(x) = |\psi(x)|^2.</math>  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
:<math>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}
</math>
Thus, for any value of ''n'' greater than one, there are regions within the box for which <math>P(x)=0</math>, 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
:<math>\langle x \rangle = \int_{-\infty}^{\infty} \psi^*(x) x \psi(x)\,\mathrm{d}x.</math>
For the particle in a box, it can be shown that the average position is always <math>\langle x \rangle = L/2</math>, regardless of the state of the particle.  In other words, the average position at which a particle in a box may be detected is exactly in the center of the quantum well; in agreement with a classical system.

==Higher-dimensional boxes==
[[Image:Particle2D.svg|thumb|320px|right|The wavefunction of a 2D well with n<sub>x</sub>=4 and n<sub>y</sub>=4]]
If a particle is trapped in a two-dimensional box, it may freely move in the <math>x</math> and <math>y</math>-directions, between barriers separated by lengths <math>L_x</math> and <math>L_y</math> 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
:<math>\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)</math>,
:<math>E_{n_x,n_y} = \frac{\hbar^2 k_{n_x,n_y}^2}{2m}</math>,
where the two-dimensional [[wavevector]] is given by
:<math>\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}}</math>.

For a three dimensional box, the solutions are
:<math>\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)</math>,
:<math>E_{n_x,n_y,n_z} = \frac{\hbar^2 k_{n_x,n_y,n_z}^2}{2m}</math>,
where the three-dimensional wavevector is given by
:<math>\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}}</math>.

An interesting feature of the above solutions is that when two or more of the lengths are the same (e.g. <math>L_x = L_y</math>), there are multiple wavefunctions corresponding to the same total energy. For example the wavefunction with <math>n_x = 2, n_y = 1</math> has the same energy as the wavefunction with <math>n_x = 1, n_y = 2</math>.  This situation is called ''[[Degenerate energy level|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.

==Related pages==
* [[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)

== References ==
{{reflist|2}}

== Bibliography ==
*{{cite book | last1=Bransden | first1=B. H. | last2 = Joachain | first2 = C. J. | title = Quantum mechanics | edition = 2nd | publisher = Pearson Education, Essex| year = 2000|isbn=0-582-35691-1}}
*{{cite book | author=Davies, John H.|title=The Physics of Low-Dimensional Semiconductors: An Introduction | publisher=Cambridge University Press|year=2006|isbn=0-521-48491-X|edition=6th reprint}}
*{{cite book | author=Griffiths, David J.|title=Introduction to Quantum Mechanics (2nd ed.) | publisher=Prentice Hall |year=2004 |isbn=0-13-111892-7}}

==Other websites==
* [http://scienceworld.wolfram.com/physics/InfiniteSquarePotentialWell.html Scienceworld] (Infinite 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]

{{DEFAULTSORT:Particle In A Box}}
[[Category:Quantum mechanics]]
[[Category:Physics]]
[[cs:Potenciálová jáma]]
[[de:Teilchen im Kasten]]
[[es:Partícula en una caja]]
[[fr:Particule dans une boîte]]
[[ko:상자속 입자]]
[[it:Particella in una scatola]]
[[he:בור פוטנציאל אינסופי]]
[[ja:井戸型ポテンシャル]]
[[pl:Cząstka w studni potencjału]]
[[pt:Partícula em uma caixa]]
[[sv:Partikel i låda]]
[[uk:Квантовий рух в прямокутній потенційній ямі]]
[[zh:無限深方形阱]]