Wavefunctions and Nodes in the 1D Box

The sine eigenfunctions and their node-counting rule are not merely a classroom exercise; they are the working model behind several genuinely quantitative applications. In quantum dot photophysics, semiconductor nanocrystals confine electrons and holes in a roughly box-like potential, and the particle-in-a-box formula for $E_n$ (extended to three dimensions) correctly predicts that smaller dots have larger HOMO–LUMO-type gaps and hence blue-shifted emission — the basis of size-tunable quantum dot displays and biological fluorescent labels. In conjugated polymer and dye chemistry, the free-electron model used in Example 2 gives surprisingly good zeroth-order estimates of absorption wavelengths for linear polyenes and cyanine dyes, precisely because the $\pi$-electron system behaves like electrons confined along the conjugation length, with nodal patterns in the frontier orbitals dictating which transitions are allowed. In spectroscopy and photochemistry, the parity/node-based selection rule derived in Example 3 is the direct ancestor of the Laporte rule governing d–d and $\pi$–$\pi^*$ transitions in centrosymmetric molecules, and of vibrational overtone selection rules in the harmonic-oscillator approximation. In scanning tunneling microscopy and surface science, standing electron waves confined by step edges or corrals on metal surfaces show exactly the sine-like nodal patterns predicted by 1D and 2D box models, directly imaged as ripples in STM topographs. Finally, in teaching and model-building, the node-counting rule here is the conceptual seed for understanding nodal surfaces in atomic orbitals (radial and angular nodes of 1s, 2s, 2p, 3d, etc.) and in molecular orbital theory, where the number of nodes in a linear combination of atomic orbitals tracks its bonding/antibonding character and hence its relative energy — the same curvature-costs-energy logic derived rigorously here in one dimension.

Ask the tutor: How would the node-counting and symmetry rules change if the 1D box had walls of finite (not infinite) potential height?

Consider a particle of mass $m$ confined to $0 \le x \le L$ by walls of infinite potential energy outside that interval, and zero potential energy inside. Inside the box the time-independent Schrödinger equation reduces to a free-particle equation,

whose general solution is a combination of sine and cosine. The infinite walls impose boundary conditions: the wavefunction must vanish at $x=0$ and at $x=L$, because the particle has zero probability of being found where the potential is infinite, and $\psi$ must be continuous. Cosine does not vanish at the origin unless its coefficient is zero, so the solution collapses to pure sine, and the condition at $x=L$ then quantizes the wavelength: only sine functions completing a half-integer number of full periods between the walls survive. This is the entirety of the quantization mechanism — no probabilistic argument, just geometry meeting a differential equation. The resulting eigenfunctions are indexed by $n=1,2,3,\dots$, and each one is a single sine arch (or a chain of arches) fit exactly between the two walls. The lowest state, $n=1$, is one smooth hump with no zero crossings between the walls — the ground state. The first excited state, $n=2$, is shaped like a full sine wave: it has one point inside the box, at the exact center, where $\psi$ crosses zero. This interior zero crossing is called a node. The state $n=3$ has two interior nodes; in general $\psi_n$ has $n-1$ of them, evenly spaced at $x = kL/n$ for $k=1,\dots,n-1$. Nodes are not a bookkeeping curiosity — they are the visible fingerprint of kinetic energy. The kinetic energy operator involves a second spatial derivative, which measures curvature. A wavefunction that must cross zero more often in the same fixed length $L$ must bend more sharply between crossings, and sharper bending means a larger (in magnitude) second derivative, hence higher kinetic energy. Since the potential energy inside the box is zero, all of the energy is kinetic, and the energy ladder $E_n = n^2 h^2/(8mL^2)$ is, in a completely literal sense, an accounting of curvature — equivalently, of nodes. A second structural feature worth isolating is symmetry. If the origin is shifted to the box's midpoint (so the walls sit at $x = -L/2$ and $x=+L/2$), every odd-$n$ wavefunction becomes an even function ($\psi(-x)=\psi(x)$, symmetric about the center), and every even-$n$ wavefunction becomes an odd function ($\psi(-x) = -\psi(x)$, antisymmetric about the center). This alternation is a direct consequence of the reflection symmetry of the box itself and will reappear, in modified form, in every centrosymmetric potential a chemist ever draws — including the harmonic oscillator and the nodal patterns of atomic $s$, $p$, $d$ orbitals.

The normalized stationary-state wavefunctions of the 1D box of length $L$ are

where: $n$ is the quantum number, $L$ is the box length, and $\sqrt{2/L}$ is the normalization constant fixing $\int_0^L |\psi_n|^2\,dx = 1$. Reading it: each $n$ selects a sine wave that completes exactly $n$ half-wavelengths across the box, automatically satisfying $\psi_n(0)=\psi_n(L)=0$. The number of interior nodes — zero crossings strictly between the walls, not counting the walls themselves — is

where: $n$ is the same quantum number as above. Reading it: the ground state ($n=1$) is nodeless; each successive excited state adds exactly one more interior node. The interior nodes sit at evenly spaced positions

where: $k$ indexes the nodes in order from the left wall. Reading it: for $\psi_n$, nodes divide the box into $n$ equal segments, each containing one arch (lobe) of the wavefunction. The associated energies are

where: $h$ is Planck's constant, $m$ is the particle mass, $L$ is the box length. Reading it: this is derived and discussed at length in the prerequisite treatment of quantization; here the point is that $E_n \propto n^2$ is the direct arithmetic consequence of the node count, made explicit in the derivation below. Symmetry about the midpoint $x_0 = L/2$ is captured by writing $u = x - L/2$ and noting

where: $u$ is displacement from the center, and the sign $(-1)^{n+1}$ is $+1$ for odd $n$ and $-1$ for even $n$. Reading it: odd-$n$ states are mirror-symmetric (even functions) about the box center, even-$n$ states are mirror-antisymmetric (odd functions) about the center.

Example 1 — Foundational: node bookkeeping for $n=5$. For $n=5$, the number of interior nodes is $n-1=4$. Their positions, from $x_k = kL/n$ with $k=1,2,3,4$, are

Check by direct substitution: at $x_2 = 2L/5$, the sine argument is $n\pi x_2/L = 5\pi(2L/5)/L = 2\pi$, and $\sin(2\pi)=0$ — confirmed. The five lobes of $\psi_5$ occupy the five equal sub-intervals $(0,L/5)$, $(L/5,2L/5)$, $(2L/5,3L/5)$, $(3L/5,4L/5)$, $(4L/5,L)$, alternating in sign: $+,-,+,-,+$. This alternation is required because $\psi_5$ must cross zero at each node, flipping sign each time it does. As a limiting-case check: setting $n=1$ in $n-1$ gives zero nodes, consistent with the nodeless ground-state hump shown in the first figure — the formula degrades correctly to the simplest case.

Example 2 — Applied: a conjugated "molecular wire" and its HOMO node count. Model a linear polyene chain — a common device for estimating $\pi$-electron transition energies — as a 1D box of length $L = 1.10\ \text{nm}$ containing $N=6$ delocalized $\pi$ electrons (three conjugated double bonds, e.g., an idealized hexatriene fragment). Filling levels two electrons at a time (Pauli exclusion, ignoring spin details beyond occupancy), the electrons occupy $n=1,2,3$, so the HOMO is $n=3$ and the LUMO is $n=4$. The HOMO wavefunction $\psi_3$ has $n-1=2$ interior nodes, located at

The LUMO $\psi_4$ has $n-1=3$ interior nodes at $L/4=0.275\ \text{nm}$, $L/2=0.550\ \text{nm}$, $3L/4=0.825\ \text{nm}$. The HOMO$\to$LUMO transition therefore adds exactly one more node to the occupied-orbital pattern, consistent with the general rule that the lowest-energy allowed excitation always adds one node. The transition energy, using electron mass $m_e = 9.109\times10^{-31}\,\text{kg}$ and $L=1.10\times10^{-9}\,\text{m}$:

Compute the denominator first:

With $h=6.626\times10^{-34}\,\text{J}\,\cdot\,\text{s}$, $h^2 = 4.394\times10^{-67}\ \text{J}^2\text{s}^2$, so

Converting to a wavelength via $\lambda = hc/\Delta E$ with $c=2.998\times10^8\,\text{m/s}$:

This falls in the visible region (blue-green edge), the right order of magnitude for conjugated polyene absorption — the crude particle-in-a-box model, built entirely on the node-counting machinery above, already predicts a physically sensible color.

Example 3 — Exam/competition level: symmetry selection rule from node parity. A dipole transition between two box states $\psi_n$ and $\psi_{n'}$ is allowed only if the transition dipole integral

is nonzero (the operator is written relative to the box center to expose the symmetry cleanly; a constant shift of origin does not change whether the integral vanishes). Determine, using only the symmetry rule $\psi_n(L/2+u) = (-1)^{n+1}\psi_n(L/2-u)$, whether the $n=2 \to n=4$ transition is allowed, and contrast with $n=2\to n=3$. Substitute $u = x - L/2$ so the integral becomes $\int_{-L/2}^{L/2}\psi_{n'}(L/2+u)\,u\,\psi_n(L/2+u)\,du$. Under $u\to -u$, each wavefunction picks up a sign $(-1)^{n+1}$ or $(-1)^{n'+1}$, and the factor $u$ itself is odd, contributing $(-1)$. The full integrand transforms under $u\to-u$ as

For the integral over the symmetric interval $[-L/2,L/2]$ to be nonzero, the integrand must be an even function of $u$, i.e., the exponent $n+n'+1$ must be even, i.e., $n+n'$ must be odd — one of $n,n'$ even and the other odd. For $n=2\to n=4$: $n+n' = 6$, which is even, so $n+n'+1=7$ is odd — the integrand is odd in $u$ and the integral vanishes over the symmetric interval. This transition is forbidden by symmetry (it is the box analogue of a "g→g" or "u→u" forbidden transition in a centrosymmetric molecule). For $n=2\to n=3$: $n+n'=5$, odd, so $n+n'+1=6$ is even — the integrand is even and the integral is generally nonzero. This transition is allowed. The general selection rule that falls out is: transitions between two box states are dipole-allowed only when $n$ and $n'$ have opposite parity (one even, one odd), exactly mirroring the Laporte-type rules for centrosymmetric molecules and atomic orbitals, where allowed transitions always connect states of opposite inversion symmetry (g$\leftrightarrow$u).

The particle-in-a-box problem does not trace to a single named discovery in the way, say, the Balmer formula does — it is the direct spatial-boundary-value analogue of standing-wave problems studied throughout 19th-century mathematical physics (vibrating strings, organ pipes, membranes). Its quantum-mechanical form crystallized in the earliest years of wave mechanics, immediately following Schrödinger's 1926 papers, as physicists and chemists sought the simplest exactly solvable potential to illustrate quantization. Erwin Schrödinger himself, in constructing wave mechanics, drew explicitly on the mathematics of vibrating strings and resonant cavities — the same boundary-value machinery that gives the box its sine eigenfunctions.