CogitaVerse

Particle in a 3D Box | Full Derivation with Equation Numbers

Particle in a Three‑Dimensional Box

Complete derivation of the Schrödinger wave equation – with numbered equations and all steps

1. Introduction

Consider a particle of mass \(m\) confined inside a rectangular box of side lengths \(a, b, c\) along the \(x, y, z\) axes. The potential is zero inside the box and infinite outside:

\[ V(x,y,z) = \begin{cases} 0 & 0 < x < a,\; 0 < y < b,\; 0 < z < c,\\ \infty & \text{otherwise}. \end{cases} \] (1)

Figure 1: Particle in a three‑dimensional rectangular box.

2. Time‑Independent Schrödinger Equation

The Hamiltonian operator in three dimensions is:

\[ \hat{H} = -\frac{\hbar^{2}}{2m}\nabla^{2} + V(x,y,z),\qquad \nabla^{2} = \frac{\partial^{2}}{\partial x^{2}}+\frac{\partial^{2}}{\partial y^{2}}+\frac{\partial^{2}}{\partial z^{2}} \] (2)

The time‑independent Schrödinger equation \(\hat{H}\psi = E\psi\) gives:

\[ -\frac{\hbar^{2}}{2m}\left(\frac{\partial^{2}\psi}{\partial x^{2}}+\frac{\partial^{2}\psi}{\partial y^{2}}+\frac{\partial^{2}\psi}{\partial z^{2}}\right) + V\psi = E\psi \] (3)

Rearranging:

\[ \frac{\partial^{2}\psi}{\partial x^{2}}+\frac{\partial^{2}\psi}{\partial y^{2}}+\frac{\partial^{2}\psi}{\partial z^{2}} + \frac{2m}{\hbar^{2}}(E-V)\psi = 0 \] (4)

3. Solution Outside the Box (\(V=\infty\))

Outside the box, \(V=\infty\). Since \(E\) is finite, equation (4) becomes:

\[ \frac{\partial^{2}\psi}{\partial x^{2}}+\frac{\partial^{2}\psi}{\partial y^{2}}+\frac{\partial^{2}\psi}{\partial z^{2}} + \frac{2m}{\hbar^{2}}(E-\infty)\psi = 0 \] (5)

The term \(\frac{2m}{\hbar^{2}}(E-\infty)\to -\infty\). For the equation to hold, the wavefunction must be zero:

\[ \psi(x,y,z) = 0 \quad \text{outside the box}. \] (6)

Thus the particle cannot exist outside. The wavefunction vanishes at the boundaries and beyond.

4. Solution Inside the Box (\(V=0\))

Inside the box, \(V=0\). Equation (4) reduces to:

\[ \frac{\partial^{2}\psi}{\partial x^{2}}+\frac{\partial^{2}\psi}{\partial y^{2}}+\frac{\partial^{2}\psi}{\partial z^{2}} + \frac{2mE}{\hbar^{2}}\psi = 0 \] (7)

Define the constant:

\[ k^{2} = \frac{2mE}{\hbar^{2}} \] (8)

Then equation (7) becomes:

\[ \frac{\partial^{2}\psi}{\partial x^{2}}+\frac{\partial^{2}\psi}{\partial y^{2}}+\frac{\partial^{2}\psi}{\partial z^{2}} + k^{2}\psi = 0 \] (9)

We solve this using separation of variables. Assume:

\[ \psi(x,y,z) = X(x)Y(y)Z(z) \] (10)

Substituting (10) into (9) and dividing by \(XYZ\):

\[ \frac{1}{X}\frac{d^{2}X}{dx^{2}} + \frac{1}{Y}\frac{d^{2}Y}{dy^{2}} + \frac{1}{Z}\frac{d^{2}Z}{dz^{2}} + k^{2} = 0 \] (11)

Since each term depends on a different variable, each must equal a constant. Write:

\[ \frac{1}{X}\frac{d^{2}X}{dx^{2}} = -k_{x}^{2},\quad \frac{1}{Y}\frac{d^{2}Y}{dy^{2}} = -k_{y}^{2},\quad \frac{1}{Z}\frac{d^{2}Z}{dz^{2}} = -k_{z}^{2} \] (12)

with

\[ k_{x}^{2} + k_{y}^{2} + k_{z}^{2} = k^{2} = \frac{2mE}{\hbar^{2}} \] (13)

Each of the equations in (12) is a one‑dimensional harmonic oscillator equation. For \(X(x)\):

\[ \frac{d^{2}X}{dx^{2}} + k_{x}^{2}X = 0 \] (14)

The general solution is \(X(x) = A\sin(k_{x}x) + B\cos(k_{x}x)\). The boundary conditions require \(\psi=0\) at \(x=0\) and \(x=a\) (since the wavefunction must vanish at the infinite walls). At \(x=0\):

\[ X(0) = B = 0 \quad\Rightarrow\quad X(x) = A\sin(k_{x}x) \] (15)

At \(x=a\):

\[ X(a) = A\sin(k_{x}a) = 0 \quad\Rightarrow\quad k_{x}a = n_{x}\pi,\quad n_{x}=1,2,3,\dots \] (16)

Thus

\[ k_{x} = \frac{n_{x}\pi}{a} \] (17)

The normalized solution for \(X(x)\) is obtained by requiring \(\int_{0}^{a} |X|^{2}dx = 1\):

\[ \int_{0}^{a} A^{2}\sin^{2}\left(\frac{n_{x}\pi x}{a}\right) dx = A^{2}\frac{a}{2} = 1 \;\Rightarrow\; A = \sqrt{\frac{2}{a}} \] (18)

Hence

\[ X_{n_{x}}(x) = \sqrt{\frac{2}{a}}\sin\left(\frac{n_{x}\pi x}{a}\right) \] (19)

Analogously, for the \(y\) and \(z\) directions:

\[ Y_{n_{y}}(y) = \sqrt{\frac{2}{b}}\sin\left(\frac{n_{y}\pi y}{b}\right),\qquad Z_{n_{z}}(z) = \sqrt{\frac{2}{c}}\sin\left(\frac{n_{z}\pi z}{c}\right) \] (20)

with \(n_{y}, n_{z}=1,2,3,\dots\).

5. Total Wavefunction and Energy Quantization

Combining (10), (19) and (20), the complete normalized wavefunction is:

\[ \psi_{n_{x}n_{y}n_{z}}(x,y,z) = \sqrt{\frac{8}{abc}}\;\sin\left(\frac{n_{x}\pi x}{a}\right)\sin\left(\frac{n_{y}\pi y}{b}\right)\sin\left(\frac{n_{z}\pi z}{c}\right) \] (21)

From (13), (17) and the analogous expressions:

\[ k^{2} = \left(\frac{n_{x}\pi}{a}\right)^{2} + \left(\frac{n_{y}\pi}{b}\right)^{2} + \left(\frac{n_{z}\pi}{c}\right)^{2} \] (22)

Using (8) \(k^{2}=2mE/\hbar^{2}\), we get:

\[ \frac{2mE}{\hbar^{2}} = \pi^{2}\left(\frac{n_{x}^{2}}{a^{2}} + \frac{n_{y}^{2}}{b^{2}} + \frac{n_{z}^{2}}{c^{2}}\right) \] (23)

Solving for \(E\):

\[ E_{n_{x}n_{y}n_{z}} = \frac{\hbar^{2}\pi^{2}}{2m}\left(\frac{n_{x}^{2}}{a^{2}} + \frac{n_{y}^{2}}{b^{2}} + \frac{n_{z}^{2}}{c^{2}}\right) \] (24)

Since \(\hbar = h/(2\pi)\), \(\hbar^{2}\pi^{2} = \frac{h^{2}}{4}\). Substituting gives the familiar form:

\[ E_{n_{x}n_{y}n_{z}} = \frac{h^{2}}{8m}\left(\frac{n_{x}^{2}}{a^{2}} + \frac{n_{y}^{2}}{b^{2}} + \frac{n_{z}^{2}}{c^{2}}\right) \] (25)

Equation (25) is the energy quantization condition. The quantum numbers \(n_{x}, n_{y}, n_{z}\) are positive integers. The ground state is \((n_{x},n_{y},n_{z})=(1,1,1)\).

Example – Cubic box (\(a=b=c=L\)):
\[ E_{n_{x}n_{y}n_{z}} = \frac{h^{2}}{8mL^{2}}\left(n_{x}^{2}+n_{y}^{2}+n_{z}^{2}\right) \] (26)

Degeneracy occurs: e.g., (2,1,1), (1,2,1), (1,1,2) all have the same energy \(\frac{6h^{2}}{8mL^{2}}\).

6. Probability Density

The probability density is the square of the wavefunction:

\[ |\psi_{n_{x}n_{y}n_{z}}|^{2} = \frac{8}{abc}\;\sin^{2}\!\left(\frac{n_{x}\pi x}{a}\right)\sin^{2}\!\left(\frac{n_{y}\pi y}{b}\right)\sin^{2}\!\left(\frac{n_{z}\pi z}{c}\right) \] (27)

This is always non‑negative and integrates to 1 over the volume of the box. The maxima occur where all sine functions equal ±1.

Figure 2: Probability density contour for the ground state in the \(xy\)-plane (\(n_{x}=n_{y}=1\), \(z= c/2\)).

7. Discussion of Degeneracy

When the box is not cubic (\(a\neq b\neq c\)), the energy levels are generally non‑degenerate because the sum \(\frac{n_{x}^{2}}{a^{2}}+\frac{n_{y}^{2}}{b^{2}}+\frac{n_{z}^{2}}{c^{2}}\) takes different values for different triples. However, if two sides are equal (e.g., \(a=b\)), degeneracy appears. In a fully cubic box, degeneracy is high and is essential for understanding the electronic structure of atoms and quantum dots.

The complete set of wavefunctions (21) forms an orthonormal basis for the Hilbert space of square‑integrable functions on the box.

8. Video Lecture (Urdu/Hindi)

Watch Complete Lecture in Urdu/Hindi for Comprehensive Understanding

Step‑by‑step derivation, separation of variables, boundary conditions, normalization, and energy quantization – explained in Urdu/Hindi.

9. Summary of Key Results

  • Wavefunction: \(\displaystyle \psi_{n_{x}n_{y}n_{z}}(x,y,z) = \sqrt{\frac{8}{abc}}\sin\left(\frac{n_{x}\pi x}{a}\right)\sin\left(\frac{n_{y}\pi y}{b}\right)\sin\left(\frac{n_{z}\pi z}{c}\right)\)
  • Energy: \(\displaystyle E_{n_{x}n_{y}n_{z}} = \frac{h^{2}}{8m}\left(\frac{n_{x}^{2}}{a^{2}}+\frac{n_{y}^{2}}{b^{2}}+\frac{n_{z}^{2}}{c^{2}}\right)\)
  • Boundary conditions force the wavefunction to vanish at the walls, leading to three quantum numbers \(n_{x},n_{y},n_{z}\in\mathbb{N}\).
  • The ground state energy is non‑zero (zero‑point energy).
  • Degeneracy arises when the box has symmetric dimensions.
\[ \int_{0}^{a}\!\int_{0}^{b}\!\int_{0}^{c} |\psi|^{2}\,dxdydz = 1 \] (28)

These results form the basis for many advanced quantum systems, including quantum wells, wires, and dots.

Complete derivation – all equations numbered, no steps omitted. Interactive diagrams and video lecture included.

Download Complete Notes Below

Leave a Comment

Your email address will not be published. Required fields are marked *

Scroll to Top