# 7.2: 3-Dimensional Models

- Page ID
- 95746

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\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]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\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]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)## Particle in a 3D Infinite Square Well

The infinite square well in three dimensions has the same property as the one-dimensional box – the potential is zero everywhere inside, and instantly becomes infinity at the boundaries. The one-dimensional case had a specified length, but we will not saddle this infinite well with the same width in all three directions, meaning we will confine the particle to a rectangular prism, not a cube. We will define our coordinate system so that the walls are parallel to the three planes, and we will place the origin at one of the box corners. The lengths of the walls along the \(x\), \(y\), and \(z\) axes we will call \(L_x\), \(L_y\), \(L_z\), respectively.

**Figure 7.2.1 Three-Dimensional Infinite Square Well**

Mathematically, the potential is written:

\[ V\left(x,y,z\right) = \left\{ \begin{array}{l} 0 & 0<x<L_x\;\;and\;\;0<y<L_y\;\;and\;\;0<z<L_z \\ \infty & \text{elsewhere} \end{array} \right. \]

As we saw for the one-dimensional box, we can use a combination of two oppositely-moving plane waves (for each of the three axes) to construct a wave function of definite energy that vanishes at the walls (i.e. sinusoidal functions). The individual solutions to the differential equations in \(x\), \(y\), and \(z\) are the same as before, with two exceptions: Each dimension involves a separate harmonic number \(n\), so like the 1-dimensional well, all the wave functions are sines:

\[ X\left(x\right) = A_x\sin\frac{n_x\pi \;x}{L_x} \;\;\;\;\; Y\left(y\right) = A_y\sin\frac{n_y\pi \;y}{L_y} \;\;\;\;\; Z\left(z\right)= A_z\sin\frac{n_z\pi \;z}{L_z},\;\;\;\;\; n_x,\;n_y,\;n_z=1,\;2,\;\dots\]

Before we move on to the energy spectrum, let's construct the spatial full wave function by multiplying the partial wave functions. We also have to deal with normalizing the full wave function. Normalization does not tell us anything about the values of \(A_x\), \(A_y\), and \(A_z\), but their product must equal the normalization constant for the full wave function. The normalization integral is over all three dimensions and integrals over \(x\) will only affect \(X\left(x\right)\), and similarly for \(y\) and \(z\), so the normalization constant for the full wave function turns out to be the same as the product of the normalization constants for the three separate one-dimensional wave functions.

\[ \psi_{n_xn_yn_z}\left(x,y,z\right) = \sqrt{\frac{8}{L_xL_yL_z}}\sin\frac{n_x\pi \;x}{L_x}\sin\frac{n_y\pi \;y}{L_y}\sin\frac{n_z\pi \;z}{L_z},\;\;\;\;\; n_x,\;n_y,\;n_z=1,\;2,\;\dots \]

In keeping with our notation of labeling the wave function with the quantum numbers, we have labeled the energy eigenstate wave function accordingly.

Plugging the wave function back into the stationary-state Schrödinger equation, we get the following energy spectrum:

\[E_{n_xn_yn_z} = \left(\frac{n_x^2}{L_x^2}+\frac{n_y^2}{L_y^2}+\frac{n_z^2}{L_z^2}\right)\frac{h^2}{8m}~,~~~n_x,\;n_y,\;n_z=1,\;2,\;\dots\]

Alert

*One might be tempted to think that the ground state energy of this particle occurs when the \(n\) along the longest dimension is 1, while the others are zero, but **It should be emphasized that the minimum value of all three \(n\) values is 1. None of the three modes can provide a zero contribution to the energy.*

Another way that the three-dimensional case differs from the one-dimensional case is apparent if we consider the hierarchy of the energy spectrum. Suppose for example, we wanted to draw an energy-level diagram for this spectrum. We know that \(\psi_{111}\) is the ground state, but which quantum state would be the first excited state? The answer depends upon the dimensions of the well. If the \(L_x\) is greater than the other two box dimensions, then the smallest increase in the total energy will come from incrementing \(n_x\) from 1 to 2, and the first excited state would be \(\psi_{211}\). What about the second excited state? Well, now we need even more information. If \(L_x\) is only slightly greater than \(L_y\) dimension (and both are longer than \(L_z\)), then the second excited state would be \(\psi_{221}\). But if \(L_x\) is significantly longer than the other dimensions, then the second excited state would be \(\psi_{311}\). In other words, with three quantum numbers, we have lost the ability (at least for this case) to express the energy levels with a single integer.

## The 3D Harmonic Oscillator

As our final example of a potential that allows for separation of variables in Cartesian coordinates, we consider the three dimensional harmonic oscillator, which has a potential that is a sum of functions purely of \(x\), \(y\), and \(z\). In general, the spring constants are different for each direction, so:

\[ V\left(x,y,z\right) = \frac{1}{2}\kappa_x x^2 + \frac{1}{2}\kappa_y y^2 + \frac{1}{2}\kappa_z z^2 \]

This sort of model might be useful for crystal lattices where the bonds differ along the different dimensions. Plugging this into Equation 7.1.10 results in an equation with three separated functions again:

\[ \left[\frac{1}{X}\frac{\partial^2X}{\partial x^2} + \frac{\kappa_x m}{\hbar^2} x^2 \right] + \left[\frac{1}{Y}\frac{\partial^2Y}{\partial y^2} + \frac{\kappa_y m}{\hbar^2} y^2 \right] + \left[\frac{1}{Z}\frac{\partial^2Z}{\partial z^2} + \frac{\kappa_z m}{\hbar^2} z^2 \right] = -\frac{2mE}{\hbar^2} \]

Following the same procedure as before gives us three separate differential equations. This decoupling maneuver once again leaves us with three wave function pieces, which are multiplied together to get the full wave function. As with the case of the square well, the energy contributions of the partial wave functions are added together to give the total energy of the state:

\[ E_{n_xn_yn_z} = \left(n_x + \frac{1}{2}\right)\hbar\sqrt{\frac{\kappa_x}{m}} + \left(n_y + \frac{1}{2}\right)\hbar\sqrt{\frac{\kappa_y}{m}} + \left(n_z + \frac{1}{2}\right)\hbar\sqrt{\frac{\kappa_z}{m}} \]

This gives us a nice way to describe bonds along different axes in general, but of particular interest is the *isotropic harmonic oscillator*, which involves equal spring constants (\(\kappa\)) in all three directions. In this case, the energy spectrum reduces to:

\[ E_{n_xn_yn_z} = \left(n_x + n_y + n_z + \frac{3}{2}\right)\hbar\omega_c~,~~~ \omega_c = \sqrt{\frac{\kappa}{m}} \]

Notice that even with the simplification of isotropy, three quantum numbers are required to define the state.

## Degeneracy

Notice that unlike one-dimensional potentials, in these cases a single quantum number does not define the energy. But there is even more to it than that. Looking at the case of the three-dimensional box again, suppose it has three equal sides: \(L_x=L_y=L_z=L\). In this case, there exist three distinct quantum states that possess the same total energy, namely \(\psi_{211}\), \(\psi_{121}\), and \(\psi_{112}\). These states clearly possess equal energies, and yet they are distinct because, for example, the states \(\psi_{211}\) and \(\psi_{121}\) yield different uncertainties in the \(x\)-component of the particle's position (\(\psi_{211}\) has two antinodes along the \(x\) axis, while \(\psi_{121}\) has only one). A similar thing occurs (only more dramatically) for the isotropic harmonic oscillator, as any combination of \(n_x\), \(n_y\), and \(n_z\) that gives the same sum will result in the same energy.

When multiple quantum states yield the same energy, they are said to be *degenerate*, and if there are a total of \(j\) distinct states for the same energy, that energy level is said to be \(j\)*-fold degenerate*. Typically degeneracy comes about due to obvious symmetries, such as in the cases mentioned above. All we need to do is rename our axes, and the states morph into each other, so naturally the energies are the same. But occasionally degeneracies arise unexpectedly, through what can only really be described as a coincidence. These are called *accidental degeneracies*. An example of one of these for the three-dimensional square well with all three sides of equal length arises for the states \(\psi_{511}\), \(\psi_{151}\), \(\psi_{115}\), and \(\psi_{333}\) – this energy level is 4-fold degenerate, rather than the "expected" 3-fold degeneracy. Naturally the first three states are not unexpectedly degenerate, but the fourth seems to come from out of nowhere.

Symmetric quantum systems are common in physics, and degeneracy follows them everywhere. This can cause difficulty in developing theory, as some internal structure can be obscured when different configurations result in the same energy spectrum. The trick then is to introduce an external perturbation that breaks the symmetry, thereby separating otherwise degenerate states. The analogous case for the cubical box would be squeezing or stretching one of the dimensions slightly. This also can work in the other direction – we might see unexpected additional spectral lines that indicate that there is additional structure present that breaks the symmetry we thought existed. So the analogous case for this is an infinite well that we think should be cubical, but provides a spectrum with energy levels landing between those that we compute.