Normalization of the Wavefunction

Every electronic structure calculation run today — Hartree-Fock, density functional theory, coupled-cluster, configuration interaction — begins with a basis set of normalized atomic orbitals (Gaussian-type or Slater-type functions), and every quantum chemistry software package (Gaussian, ORCA, Psi4, NWChem) internally renormalizes basis functions and computes overlap matrices $S_{ij} = \int \phi_i^*\phi_j\,d\tau$ whose diagonal elements must equal 1 for correctly normalized basis functions — a routine sanity check built into every self-consistent-field calculation. In spectroscopy, normalized wavefunctions are what make transition dipole moments and absorption intensities quantitatively meaningful: the intensity of a UV-Vis or IR transition is proportional to $|\langle\psi_f|\hat\mu|\psi_i\rangle|^2$, and this number is only physically interpretable as a transition probability if both $\psi_i$ and $\psi_f$ are properly normalized beforehand. NMR chemical shift and coupling constant calculations, similarly, rest on normalized spin and spatial wavefunctions throughout. In materials science and nanotechnology, normalization underlies the design of quantum dots and quantum wells, where particle-in-a-box-like confinement determines emission color; the normalization constant $\sqrt{2/L}$ derived above appears directly in calculating the radiative recombination rates and absorption cross-sections of semiconductor nanocrystals used in QLED displays and biological fluorescent tags. In quantum computing, normalized state vectors ($\sum_i|c_i|^2=1$ for qubit superpositions) are the literal software objects manipulated by quantum circuits, making the discrete form of the normalization condition an everyday computational constraint rather than a textbook abstraction.

Ask the tutor: Once I've normalized a wavefunction, how do I use it to compute an expectation value like the average energy or average position?

Think of $|\psi(x)|^2\,dx$ as answering a concrete question: if you searched for the particle in the sliver of space between $x$ and $x+dx$, how likely are you to find it there? Since the particle exists somewhere with certainty, adding up (integrating) this probability over every sliver of space — the entire domain accessible to the particle — must give exactly 1. If instead the raw output of the Schrödinger equation gives you $\int|\psi|^2\,dx = 7$, that is not a contradiction of physics; it just means your $\psi$ has the right shape but the wrong scale. Multiply it by $1/\sqrt7$ and the new function $N\psi$ has all the same nodes, curvature, and relative probabilities, but now integrates to exactly 1. This rescaling freedom exists because the time-independent Schrödinger equation,

is linear: if $\psi$ solves it, so does $c\psi$ for any constant $c$, real or complex, with the same eigenvalue $E$. The differential equation fixes the functional form — the nodes, the decay rate, the oscillation pattern — but the amplitude is a free parameter left for the physical boundary condition (total probability = 1) to pin down. This is precisely analogous to solving $y'' = -y$ and getting $y = A\sin x$: the equation doesn't care what $A$ is, but a physical constraint (say, a fixed initial displacement) will. Normalization is a global statement — one number, one condition — that constrains the entire function at once, in contrast to the local conditions (continuity of $\psi$, continuity of $d\psi/dx$ except at infinite discontinuities of $V$) that were already imposed to make $\psi$ an "acceptable" wavefunction in the first place. A function can be perfectly continuous, single-valued, and finite everywhere and still fail to be a valid quantum state if $\int|\psi|^2 d\tau$ diverges — such functions (plane waves $e^{ikx}$ over infinite space, for instance) are not normalizable in the strict sense, which is why free-particle problems require special treatment (wave packets, box normalization, or the Dirac delta normalization used for continuum states). The picture generalizes immediately to three dimensions and to many particles. For a single particle in three dimensions,

and in spherical coordinates, appropriate for atoms, the volume element $d\tau = r^2\sin\theta\,dr\,d\theta\,d\phi$ carries an explicit $r^2\sin\theta$ weight that must not be forgotten — a classic source of factor-of-something errors when normalizing atomic orbitals. For an $N$-electron wavefunction $\Psi(\mathbf{r}_1,\ldots,\mathbf{r}_N)$, normalization means integrating $|\Psi|^2$ over all electron coordinates simultaneously, and it underlies the interpretation of $|\Psi|^2 d\tau_1 \cdots d\tau_N$ as the joint probability of electron 1 in $d\tau_1$, electron 2 in $d\tau_2$, and so on.

The normalization condition, in its most general form for a single particle, is

where: $\psi^*$ is the complex conjugate of the wavefunction, $|\psi|^2 = \psi^*\psi$ is the probability density, $d\tau$ is the volume element appropriate to the coordinate system, and the integral runs over the entire domain accessible to the particle (all of space, or the length of a box, or all angles on a sphere).

Reading: this equation says nothing about the shape of $\psi$ — it is one scalar constraint that fixes only the overall magnitude of the function. If $\phi(\mathbf r)$ is a solution of the Schrödinger equation that is not normalized (i.e. $\int|\phi|^2 d\tau = S \ne 1$ for some finite positive number $S$), then the physically correct, normalized wavefunction is

where: $N$ is the normalization constant, real and positive by convention, and $S$ is the "raw" normalization integral computed directly from the un-normalized function.

Reading: you compute $S$ once by brute-force integration, and $N = S^{-1/2}$ rescales $\phi$ so that $\int|N\phi|^2 d\tau = N^2 S = 1$ automatically. For a superposition (linear combination) of orthonormal basis functions $\phi_i$ with coefficients $c_i$,

where: orthonormal means $\int \phi_i^*\phi_j\,d\tau = \delta_{ij}$ (the Kronecker delta), and $|c_i|^2$ is interpreted as the probability that a measurement of the corresponding observable yields the eigenvalue associated with $\phi_i$.

Reading: this is the discrete analogue of $\int|\psi|^2 d\tau=1$ — normalization becomes a statement that probabilities of distinct outcomes sum to one, foreshadowing the full machinery of Expectation Values built from expansions in a complete orthonormal basis. Finally, normalization is preserved automatically by time evolution under the time-dependent Schrödinger equation,

where: $\Psi(\mathbf r, t)$ is a solution of $i\hbar\,\partial \Psi/\partial t = \hat H \Psi$, and the zero follows from $\hat H$ being Hermitian.

Reading: normalize once at $t=0$ and, so long as $\hat H$ is Hermitian (true for any physically sensible Hamiltonian), the wavefunction remains normalized for all later times — probability is conserved, never created or destroyed.

Example 1 — Foundational: normalizing an asymmetric linear combination in a box. Consider a particle in a one-dimensional box of length $L = 1.20\ \text{nm}$, prepared not in a single stationary state but in the superposition

where $\psi_1,\psi_2$ are the (already normalized) $n=1,2$ box eigenfunctions from the derivation above, and both $c_1,c_2$ are taken real for simplicity. Because $\psi_1$ and $\psi_2$ are orthonormal ($\int\psi_1\psi_2\,dx = 0$, $\int\psi_i^2 dx = 1$), the cross term vanishes upon integrating $|\Phi|^2$:

Normalization requires $c_1^2 + c_2^2 = 1$, so

The normalized superposition is $\Phi(x) = 0.60\,\psi_1(x) + 0.80\,\psi_2(x)$, and $|c_1|^2 = 0.36$, $|c_2|^2 = 0.64$ are the probabilities that an energy measurement returns $E_1$ or $E_2$ respectively — they sum to 1.00, as required. Note the asymmetry deliberately chosen (0.60/0.80, a 3-4-5 triple) rather than the equal-mixing case $c_1=c_2=1/\sqrt2$, to confirm the Pythagorean-style relation is being applied correctly rather than masked by symmetry.

Example 2 — Applied: normalizing a hydrogen 2p$_z$ orbital used in bonding descriptions. The (un-normalized, radial-times-angular) hydrogen 2p$_z$ orbital has the form

a shape used routinely when constructing $\sigma$ and $\pi$ molecular orbitals in valence-bond and Hückel-adjacent pictures. Find the normalization constant $N$ such that $\psi_{2p_z} = N\phi_{2p_z}$ satisfies $\int|\psi_{2p_z}|^2 d\tau = 1$. Separate the integral into radial and angular parts, since $\phi_{2p_z} = r e^{-r/2a_0} \times \cos\theta$:

The radial factor is $r^4 e^{-r/a_0}$ (careful: $r^2$ from $d\tau$ times $r^2$ from $\phi^2$), so using $\int_0^\infty r^n e^{-kr}dr = n!/k^{n+1}$ with $n=4$, $k=1/a_0$:

The polar-angle integral uses the substitution $u=\cos\theta$, $du = -\sin\theta\,d\theta$:

The azimuthal integral gives $\int_0^{2\pi}d\phi = 2\pi$. Multiplying the three pieces,

Solving,

The normalized orbital is

matching the standard hydrogenic orbital tables once the conventional radial prefactor $1/\sqrt{32\pi a_0^3}$ is written with $r$ measured in units of $a_0$ folded in — the key chemical point being that the same normalization procedure (separate radial and angular integrals, multiply, invert) applies to every hydrogenic orbital regardless of $\ell$ and $m$.

Example 3 — Exam/competition level: normalizing a Gaussian wave packet and checking the free-particle limit. A localized wave packet used to model an electron just after a fast ionization event is modeled (at $t=0$) by the un-normalized Gaussian

with $k_0$ an arbitrary real wavevector. Because $|\phi(x)|^2 = e^{-\alpha(x-x_0)^2}$ (the phase factor $e^{ik_0x}$ has modulus 1 and drops out of $|\phi|^2$), normalization requires

Shifting variables $u = x - x_0$ (the shift doesn't change the limits, since they are $\pm\infty$), and using the standard Gaussian integral $\int_{-\infty}^\infty e^{-\alpha u^2}\,du = \sqrt{\pi/\alpha}$,

Plugging in $\alpha = 4.50\ \text{nm}^{-2}$:

So the normalized packet is $\psi(x) = 1.0938\ \text{nm}^{-1/2}\, e^{-2.25(x-0.75\,\text{nm})^2/\text{nm}^2}\, e^{ik_0 x}$, independent of $x_0$ and $k_0$ (only the width $\alpha$ sets the scale — shifting the packet's center or boosting it to a different momentum never changes how much normalization constant it needs, since both operations only multiply $\phi$ by a phase or translate the peak, neither of which affects $\int|\phi|^2 dx$).

Limiting case check. As the packet is allowed to spread out — physically, $\alpha \to 0^+$, meaning the Gaussian becomes extremely broad and approaches an (unnormalizable) plane wave — the normalization constant behaves as

This is exactly the expected pathology: a genuine plane wave $e^{ik_0x}$ (the $\alpha \to 0$ limit) cannot be normalized over infinite space at all ($\int_{-\infty}^\infty 1\,dx$ diverges), and the vanishing of $N$ signals that the true normalized state degenerates into something requiring delta-function (continuum) normalization rather than the ordinary $\int|\psi|^2=1$ condition used here. Checking the opposite limit, $\alpha \to \infty$ (packet collapsing to a point), gives $N \to \infty$, consistent with squeezing all the probability into a vanishingly narrow spike, which needs an ever-larger amplitude to keep the area under $|\psi|^2$ fixed at 1.

Normalization as a formal requirement traces back to Max Born's 1926 proposal that $|\psi|^2$ be interpreted as a probability density — a genuinely radical reading of a quantity that Schrödinger himself had originally hoped to interpret as a literal charge or mass density smeared through space. Born's insight, offered almost as a brief remark in a paper on collision theory, supplied the missing physical postulate that turns Schrödinger's differential equation into a predictive theory of measurement outcomes, and normalization is the direct, unavoidable mathematical consequence of taking that postulate seriously.