$$\require{cancel}$$

# 13.7: Particle in a 2-dimensional box

$$\newcommand{\vecs}{\overset { \rightharpoonup} {\mathbf{#1}} }$$ $$\newcommand{\vecd}{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}{\| #1 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}{\| #1 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$

We learned from solving Schrödinger’s equation for a particle in a one-dimensional box that there is a set of solutions, the stationary states, for which the time dependence is just an overall rotating phase factor, and these solutions correspond to definite values of the energy. An alternative way of finding that set of solutions is separation of variables. The basic strategy is to assume that the solution to the wave equation can be factored into a product of two functions, one depending only on time, the other on the spatial variable,

$\Psi(x,t) = \psi(x)\varphi(t)$

If this solution is substituted in the Schrödinger equation, and the result divided by $$\Psi(x,t)$$, we find

$i \hbar \dfrac{\dfrac{\partial \phi(t)}{\partial t}}{\phi(t)} = \dfrac{- \dfrac{\hbar^2}{2m} \dfrac{\partial^2 \psi(x)}{\partial x^2} + V(x)\psi(x)}{\psi(x)}$

On writing the equation in this form, it is clear that the left hand side is only a function of $$t$$, not of $$x$$, and the right hand side is only a function of $$x$$! This can only make sense if in fact both sides are the same constant. If we denote this constant by $$E$$, we can write two equations:

$i \hbar \dfrac{\partial \phi(t)}{\partial t} = E \phi(t)$

$- \dfrac{\hbar^2}{2m} \dfrac{\partial^2 \psi(x)}{\partial x^2} + V(x)\psi(x) = E \psi(x)$

The solution to the first equation just gives the phase time dependence,

$\varphi(t) = A e^{-iEt/k}$

and the second is the time independent Schrödinger equation as before. The solutions to this equation are determined by the boundary conditions on $$ψ$$, in general there is a sequence of such eigenstates labeled by a quantum number n = 0,1,2,3,…, with corresponding values $$E_0$$, $$E_1$$, …, which are put in the corresponding $$\varphi(t)$$.

## A Two Dimensional Box

Let us now consider the Schrödinger equation for an electron confined to a two dimensional box, $$0 < x < a$$ and $$0 < y < b$$. That is to say, within this rectangle the electron wavefunction behaves as a free particle ($$V(x,y) = 0$$), but the walls are impenetrable so the wavefunction $$\Psi(x,y,t)=0$$ at the walls. What do we expect the wavefunction to look like?

First notice that the separation of variables trick given above for one dimension works equally well here, writing

$\Psi(x,y,t) = \psi(x,y) \varphi(t)$

gives

$i \hbar \dfrac{\partial \phi(t)}{\partial t} = E \phi(t)$

$- \dfrac{\hbar^2}{2m} \left (\dfrac{\partial^2 \psi(x)}{\partial x^2} + \dfrac{\partial ^2 \psi(x)}{\partial y^2} \right) = E \psi(x)$

The surprising thing at this point is that we can do the separation of variables trick again—we can write

$\psi(x,y) = f(x)g(t)$

and substitute in the above equation to find

$- \dfrac{\hbar^2}{2m} \left (\dfrac{\dfrac{\partial^2 f(x)}{\partial x^2}}{f(x)} + \dfrac{\dfrac{\partial ^2 g(x)}{\partial y^2}}{g(y)} \right) = E$

Again, we have an equation in which only one term is x-dependent, so it must be a constant (which we take to be negative for future convenience),

$\dfrac{\dfrac{\partial^2 f(x)}{\partial x^2}}{f(x)} = -C$

so

$\dfrac{\partial^2 f(x)}{\partial x^2} = -C f(x)$

This is exactly the same equation we dealt with in the one dimensional case, so we know

$f(x) = A \sin \dfrac{n\pi x}{a}$

with $$n$$ an integer, and the constant $$C$$ is equal to $$n^2π^2/a^2$$. Hence the energy levels in this rectangular well are given by

$E_n = \dfrac{\hbar^2}{2m} \left( \dfrac{n_x^2\pi^2}{a^2} +\dfrac{n_y^2\pi^2}{b^2} \right)$

with $$n_x$$ and $$n_y$$ are the two quantum numbers needed to label each state.

The first nice probability distributions for a particle in a 2D box are viewable below.

wiki.page("Visualizations_and_Simulations/CalcPlot3D/Probability_Wave_Function")

Figure $$\PageIndex{1}$$: Visualization of the probability for a particle in a 2D box. The quantum numbers ($$n_x$$ and $$n_y$$) can be varied in the upper left.

## Degeneracy

Two distinct wavefunctions are said to be degenerate if they correspond to the same energy. If the sides a, b of the rectangle are such that a/b is irrational (the general case), there will be no degeneracies. The most degenerate case is the square, a = b, for which clearly $$E_{m,n} = E_{n,m}$$. Degeneracies in quantum physics are most often associated with symmetries in this way.

Figure $$\PageIndex{2}$$: Contours of the $$n_x=2$$ and $$n_y=3$$ wavefunction (left) and $$n_x=3$$ and $$n_y=2$$ wavefunction (right).

We give here examples of wavefunctions (3,2) and (2,3) for a rectangle. These are contour maps for the time-independent solution, with white being the highest point. These two wavefunctions do not correspond to the same energy, although they would, of course, for a square.