Skip to content

Normalization

Normalization fixes the total weight of a quantum state so that the Born rule defines a probability law. For an ordinary pure-state vector,

⟨ψ∣ψ⟩=1.\langle\psi|\psi\rangle=1.

This condition is not cosmetic. If a complete measurement is applied to a unit vector, its outcome probabilities sum to one. The same abstract condition appears as a coefficient sum in a discrete basis, an integral for a wavefunction, or a trace condition for a density operator.

Several conventions that share the word “normalization” must be kept distinct:

ObjectNormalization statementMeaning
ordinary pure state⟨ψ∣ψ⟩=1\langle\psi\rvert\psi\rangle=1total probability one
density operatorTr⁡ρ=1\operatorname{Tr}\rho=1total probability one
conditional branchTr⁡ρ~=p≤1\operatorname{Tr}\widetilde\rho=p\le1branch probability retained
continuum basis ket⟨a∣a′⟩=δ(a−a′)\langle a\rvert a'\rangle=\delta(a-a')distributional basis convention
box mode⟨n∣m⟩=δnm\langle n\rvert m\rangle=\delta_{nm}finite-volume basis convention
scattering modeprescribed incident fluxcomparison of rates or cross sections

Only the first two describe normalized physical states directly. Delta, box, and flux normalization organize idealized basis states or scattering solutions. A subnormalized operator retains the probability that a selected branch occurred.

Let a projective measurement have mutually orthogonal projectors {Pa}\{P_a\} with

∑aPa=I.\sum_a P_a=I.

For a vector ∣ψ⟩|\psi\rangle, the Born weights are

w(a)=⟨ψ∣Pa∣ψ⟩.w(a)=\langle\psi|P_a|\psi\rangle.

Summing over all outcomes gives

∑aw(a)=⟨ψ∣(∑aPa)∣ψ⟩=⟨ψ∣ψ⟩.\begin{aligned} \sum_a w(a) &= \langle\psi| \left(\sum_a P_a\right) |\psi\rangle \\ &= \langle\psi|\psi\rangle. \end{aligned}

Thus the weights form a normalized probability distribution precisely when ∥ψ∥=1\|\psi\|=1. For a normalized POVM {Ea}\{E_a\}, the same argument uses ∑aEa=I\sum_aE_a=I.

If ∣χ⟩|\chi\rangle is nonzero but not normalized, probabilities can be written as ratios:

p(a)=⟨χ∣Pa∣χ⟩⟨χ∣χ⟩.p(a) = \frac{\langle\chi|P_a|\chi\rangle} {\langle\chi|\chi\rangle}.

This ratio is invariant under any nonzero complex rescaling ∣χ⟩↦λ∣χ⟩|\chi\rangle\mapsto\lambda|\chi\rangle. Unit normalization packages the denominator once so subsequent probability formulas are simpler.

The Born Rule owns the probability postulate itself. This page owns the scale condition that makes the postulate produce total probability one.

For any nonzero vector ∣χ⟩|\chi\rangle with finite norm,

∥χ∥=⟨χ∣χ⟩>0.\|\chi\| = \sqrt{\langle\chi|\chi\rangle}>0.

A normalized representative is

∣ψ⟩=∣χ⟩∥χ∥.|\psi\rangle = \frac{|\chi\rangle}{\|\chi\|}.

Indeed,

⟨ψ∣ψ⟩=⟨χ∣χ⟩∥χ∥2=1.\langle\psi|\psi\rangle = \frac{\langle\chi|\chi\rangle} {\|\chi\|^2} =1.

The zero vector cannot be normalized because its norm vanishes. It represents no physical pure state: every Born weight computed from it is zero, and it does not define a ray.

Normalization also fails when ∥χ∥=∞\|\chi\|=\infty. Such an expression may be a generalized eigenfunction, a useful asymptotic solution, or simply an inadmissible candidate state. Dividing it by an informal “infinite norm” does not create a Hilbert-space vector.

A physical pure state is a ray rather than one preferred vector. If

∣χ′⟩=λ∣χ⟩,λ≠0,|\chi'\rangle = \lambda|\chi\rangle, \qquad \lambda\ne0,

then the two vectors define the same ray. Normalizing each gives

∣χ′⟩∥χ′∥=λ∣λ∣∣χ⟩∥χ∥.\frac{|\chi'\rangle}{\|\chi'\|} = \frac{\lambda}{|\lambda|} \frac{|\chi\rangle}{\|\chi\|}.

The magnitude ∣λ∣|\lambda| disappears, while its phase remains as a global phase. Therefore normalization selects the unit sphere inside the nonzero vectors but does not select one unique vector on each ray.

The sequence of quotients is conceptually useful:

nonzero vectors⟶unit vectors⟶physical rays,\begin{aligned} \text{nonzero vectors} &\longrightarrow \text{unit vectors} \\ &\longrightarrow \text{physical rays}, \end{aligned}

where the first step removes positive scale and the second identifies constant phases. The geometry of the final quotient belongs to Projective Hilbert Space.

In an orthonormal basis {∣n⟩}n=1d\{|n\rangle\}_{n=1}^d, write

∣χ⟩=∑n=1dan∣n⟩.|\chi\rangle = \sum_{n=1}^d a_n|n\rangle.

Orthonormality gives

∥χ∥2=∑n=1d∣an∣2.\|\chi\|^2 = \sum_{n=1}^d|a_n|^2.

If the sum is nonzero, normalized coefficients are

cn=an∑m∣am∣2.c_n = \frac{a_n} {\sqrt{\sum_m|a_m|^2}}.

For example,

∣χ⟩=2∣0⟩−i∣1⟩|\chi\rangle = 2|0\rangle-i|1\rangle

has squared norm 55, so

∣ψ⟩=25∣0⟩−i5∣1⟩.|\psi\rangle = \frac{2}{\sqrt5}|0\rangle - \frac{i}{\sqrt5}|1\rangle.

The basis probabilities are 4/54/5 and 1/51/5. The relative phase −i-i does not affect this particular measurement but can affect measurements in another basis.

Nonorthogonal Coordinates and the Gram Matrix

Section titled “Nonorthogonal Coordinates and the Gram Matrix”

The coefficient-sum rule is not valid in an arbitrary basis. Let {∣bi⟩}\{|b_i\rangle\} be a nonorthogonal basis and

∣χ⟩=∑ici∣bi⟩.|\chi\rangle = \sum_i c_i|b_i\rangle.

Define the Gram matrix

Gij=⟨bi∣bj⟩.G_{ij}=\langle b_i|b_j\rangle.

Then

∥χ∥2=∑i,jci∗Gijcj=c†Gc.\|\chi\|^2 = \sum_{i,j}c_i^*G_{ij}c_j = c^\dagger Gc.

The normalized coordinate column is therefore

cnorm=cc†Gc.c_{\mathrm{norm}} = \frac{c}{\sqrt{c^\dagger Gc}}.

For linearly independent basis vectors, GG is positive definite, so c†Gc>0c^\dagger Gc>0 for every nonzero column. If the spanning set is linearly dependent, GG is singular and coordinate descriptions are not unique; one must first identify the actual vector or remove redundancy.

This example emphasizes the invariant rule: normalize using the Hilbert-space inner product, not by applying a remembered component formula outside its assumptions.

For a countable orthonormal basis,

∣χ⟩=∑n=1∞an∣n⟩,|\chi\rangle = \sum_{n=1}^{\infty}a_n|n\rangle,

the vector belongs to the Hilbert space only if

∑n=1∞∣an∣2<∞.\sum_{n=1}^{\infty}|a_n|^2<\infty.

It can be normalized when this sum is positive and finite. The sequence an=1/na_n=1/\sqrt n, for example, does not define an element of ℓ2\ell^2 because

∑n=1∞1n\sum_{n=1}^{\infty}\frac1n

diverges. An expression may be a formally meaningful series without being a normalizable state.

For a normalized state, truncating after NN terms gives the subnormalized vector

∣ψN⟩=∑n=1Ncn∣n⟩|\psi_N\rangle = \sum_{n=1}^{N}c_n|n\rangle

with retained weight

wN=⟨ψN∣ψN⟩=∑n=1N∣cn∣2.w_N = \langle\psi_N|\psi_N\rangle = \sum_{n=1}^{N}|c_n|^2.

If wN>0w_N>0, the normalized truncation is

∣ψ^N⟩=∣ψN⟩wN.|\widehat\psi_N\rangle = \frac{|\psi_N\rangle}{\sqrt{w_N}}.

Renormalizing makes a valid state inside the truncated subspace, but it does not erase truncation error. The discarded probability weight is 1−wN1-w_N, and the fidelity with the original pure state is wNw_N.

For a particle on the line, normalization is

∫−∞∞dx ∣ψ(x)∣2=1.\int_{-\infty}^{\infty} dx\,|\psi(x)|^2=1.

If

ψ(x)=Af(x),\psi(x)=A f(x),

then

∣A∣=(∫dx ∣f(x)∣2)−1/2,|A| = \left( \int dx\,|f(x)|^2 \right)^{-1/2},

provided the integral is finite and nonzero. Normalization fixes ∣A∣|A| but not the constant phase of AA.

For example, with α>0\alpha>0,

ψ(x)=Ae−αx2/2\psi(x)=Ae^{-\alpha x^2/2}

has

∫dx e−αx2=πα,\int dx\,e^{-\alpha x^2} = \sqrt{\frac\pi\alpha},

so

∣A∣=(απ)1/4.|A| = \left(\frac\alpha\pi\right)^{1/4}.

The function e+αx2/2e^{+\alpha x^2/2} cannot be normalized on the line because its squared modulus is not integrable. No finite constant A≠0A\ne0 repairs that failure.

Square Integrability and Almost-Everywhere Equality

Section titled “Square Integrability and Almost-Everywhere Equality”

Normalizable position wavefunctions live in an L2L^2 space. The condition

∫dx ∣ψ(x)∣2<∞\int dx\,|\psi(x)|^2<\infty

states square integrability, while equality to one selects a unit-norm representative.

An L2L^2 vector is an equivalence class of functions that agree almost everywhere. Changing a wavefunction at finitely many points does not alter its norm or any interval probability. Conversely, pointwise decay is not by itself enough to guarantee square integrability; the rate and measure matter.

For a power-law tail

∣ψ(x)∣∼∣x∣−s|\psi(x)|\sim |x|^{-s}

on the line, the tail contribution behaves as

∫R∞dx x−2s,\int_R^\infty dx\,x^{-2s},

which converges only for s>1/2s>1/2. Local singularities require a separate integrability check. The mathematical structure is canonical in L2 Spaces.

Normalization belongs to a function and a measure together. In coordinates q=(q1,…,qd)q=(q^1,\ldots,q^d), one may have

1=∫dq J(q)∣ψ(q)∣2,1 = \int dq\,J(q)|\psi(q)|^2,

where J(q)J(q) is a Jacobian density. In spherical coordinates,

1=∫0∞dr∫0πdθ×∫02πdφ r2sin⁡θ∣ψ(r,θ,φ)∣2.\begin{aligned} 1 &= \int_0^\infty dr \int_0^\pi d\theta \\ &\quad\times \int_0^{2\pi}d\varphi\, r^2\sin\theta |\psi(r,\theta,\varphi)|^2. \end{aligned}

The radial factor and angular Jacobian are part of the probability measure. Dropping them changes the norm.

If ψ\psi is normalized with respect to ddqd^dq, then

[ψ]=[q]−d/2[\psi]=[q]^{-d/2}

when all coordinates carry the same dimension. For a one-particle position wavefunction in dd spatial dimensions,

[ψ]=L−d/2.[\psi]=L^{-d/2}.

The abstract ket has no position-space units of this kind; the units belong to its coordinate representative and chosen continuum normalization.

One can absorb the Jacobian into a redefined function:

ψ~(q)=J(q)1/2ψ(q),\widetilde\psi(q) = J(q)^{1/2}\psi(q),

so that

1=∫dq ∣ψ~(q)∣2.1=\int dq\,|\widetilde\psi(q)|^2.

The represented operators must be transformed at the same time.

With the site’s symmetric Fourier convention,

ϕ(p)=12πℏ∫dx e−ipx/ℏψ(x),\phi(p) = \frac{1}{\sqrt{2\pi\hbar}} \int dx\,e^{-ipx/\hbar}\psi(x),

and a normalized state satisfies

∫dp ∣ϕ(p)∣2=1.\int dp\,|\phi(p)|^2=1.

The equality

∫dx ∣ψ(x)∣2=∫dp ∣ϕ(p)∣2\int dx\,|\psi(x)|^2 = \int dp\,|\phi(p)|^2

is norm preservation under a unitary Fourier transform. It is not an independent normalization condition on a second state.

The momentum wavefunction has units

[ϕ(p)]=[p]−1/2[\phi(p)]=[p]^{-1/2}

in one dimension. Relabelling momentum by p=ℏkp=\hbar k requires

ψ~(k)=ℏ ϕ(ℏk)\widetilde\psi(k) = \sqrt\hbar\,\phi(\hbar k)

so that

∣ψ~(k)∣2dk=∣ϕ(p)∣2dp.|\widetilde\psi(k)|^2dk = |\phi(p)|^2dp.

The Jacobian is part of normalization. The full representation dictionary is in Momentum-Space Representation.

For a wavefunction with a discrete internal label ss, normalization sums over that label and integrates the continuous variables:

∑s∫dx ∣ψs(x)∣2=1.\sum_s \int dx\,|\psi_s(x)|^2=1.

For a spin-1/21/2 particle,

Ψ(x)=(ψ↑(x)ψ↓(x)),\Psi(x) = \begin{pmatrix} \psi_\uparrow(x)\\ \psi_\downarrow(x) \end{pmatrix},

so

∫dx Ψ(x)†Ψ(x)=∫dx ∣ψ↑(x)∣2+∫dx ∣ψ↓(x)∣2=1.\begin{aligned} \int dx\, \Psi(x)^\dagger\Psi(x) &= \int dx\, |\psi_\uparrow(x)|^2 \\ &\quad+ \int dx\, |\psi_\downarrow(x)|^2 \\ &=1. \end{aligned}

The individual components need not each have norm one. Their squared norms are the probabilities of the corresponding internal outcomes when the continuous coordinate is ignored.

For NN distinguishable particles in three dimensions,

1=∫∏j=1Nd3rj ∣Ψ(r1,…,rN)∣2.1 = \int \prod_{j=1}^{N}d^3r_j\, |\Psi(\mathbf r_1,\ldots,\mathbf r_N)|^2.

The wavefunction lives on 3N3N-dimensional configuration space and has units

[Ψ]=L−3N/2.[\Psi]=L^{-3N/2}.

If there are spin labels, normalization also sums over every spin configuration. Exchange symmetry for identical particles restricts the allowed functions but does not change the unit-norm condition.

For a product state

∣Ψ⟩=∣ψ⟩⊗∣ϕ⟩,|\Psi\rangle = |\psi\rangle\otimes|\phi\rangle,

the norm factors:

∥Ψ∥2=∥ψ∥2∥ϕ∥2.\|\Psi\|^2 = \|\psi\|^2\|\phi\|^2.

Normalized factors therefore give a normalized product. Entangled states are normalized by their full joint norm, not by trying to normalize nonexistent separate subsystem wavefunctions.

A general quantum state is represented by a positive trace-one operator:

ρ≥0,Tr⁡ρ=1.\rho\ge0, \qquad \operatorname{Tr}\rho=1.

For any complete POVM {Ea}\{E_a\},

∑ap(a)=∑aTr⁡(ρEa)=Tr⁡(ρ∑aEa)=Tr⁡ρ=1.\begin{aligned} \sum_a p(a) &= \sum_a\operatorname{Tr}(\rho E_a) \\ &= \operatorname{Tr} \left( \rho\sum_aE_a \right) \\ &= \operatorname{Tr}\rho =1. \end{aligned}

For a pure normalized vector,

ρψ=∣ψ⟩⟨ψ∣\rho_\psi = |\psi\rangle\langle\psi|

has

Tr⁡ρψ=⟨ψ∣ψ⟩=1.\operatorname{Tr}\rho_\psi = \langle\psi|\psi\rangle =1.

For an ensemble

ρ=∑jqj∣ψj⟩⟨ψj∣,\rho = \sum_j q_j |\psi_j\rangle\langle\psi_j|,

with normalized ∣ψj⟩|\psi_j\rangle, trace normalization requires

qj≥0,∑jqj=1.q_j\ge0, \qquad \sum_jq_j=1.

This decomposition is not unique, but Tr⁡ρ=1\operatorname{Tr}\rho=1 is representation independent. The state concept and positivity requirements are canonical in Density Operators.

Intermediate states need not always carry total weight one. Suppose a measurement outcome aa is represented by an operation that produces

ρ~a=MaρMa†.\widetilde\rho_a = M_a\rho M_a^\dagger.

Its trace is the probability of that branch:

p(a)=Tr⁡ρ~a≤1.p(a) = \operatorname{Tr}\widetilde\rho_a \le1.

The tilde signals that the operator is subnormalized. Conditional on the outcome occurring and p(a)>0p(a)>0, the normalized state is

ρa=ρ~ap(a).\rho_a = \frac{\widetilde\rho_a}{p(a)}.

For a pure-state projection,

∣ψ~a⟩=Pa∣ψ⟩|\widetilde\psi_a\rangle = P_a|\psi\rangle

has squared norm

⟨ψ~a∣ψ~a⟩=⟨ψ∣Pa∣ψ⟩=p(a),\langle\widetilde\psi_a |\widetilde\psi_a\rangle = \langle\psi|P_a|\psi\rangle = p(a),

and the conditional unit vector is

∣ψa⟩=Pa∣ψ⟩p(a).|\psi_a\rangle = \frac{P_a|\psi\rangle}{\sqrt{p(a)}}.

One must not renormalize a collection of branches before recording their traces, because those traces are the outcome probabilities. The dynamical and interpretive statement belongs to State Update Rule.

Closed-system time evolution is unitary:

∣ψ(t)⟩=U(t,t0)∣ψ(t0)⟩,U†U=I.|\psi(t)\rangle = U(t,t_0)|\psi(t_0)\rangle, \qquad U^\dagger U=I.

Therefore

⟨ψ(t)∣ψ(t)⟩=⟨ψ(t0)∣U†U∣ψ(t0)⟩=⟨ψ(t0)∣ψ(t0)⟩.\begin{aligned} \langle\psi(t)|\psi(t)\rangle &= \langle\psi(t_0)| U^\dagger U |\psi(t_0)\rangle \\ &= \langle\psi(t_0)|\psi(t_0)\rangle. \end{aligned}

A normalized initial state remains normalized. In differential form, for a self-adjoint Hamiltonian and a state in the relevant domains,

ddt⟨ψ∣ψ⟩=⟨ψ˙∣ψ⟩+⟨ψ∣ψ˙⟩=iℏ⟨ψ∣H∣ψ⟩−iℏ⟨ψ∣H∣ψ⟩=0.\begin{aligned} \frac{d}{dt}\langle\psi|\psi\rangle &= \langle\dot\psi|\psi\rangle + \langle\psi|\dot\psi\rangle \\ &= \frac{i}{\hbar} \langle\psi|H|\psi\rangle - \frac{i}{\hbar} \langle\psi|H|\psi\rangle \\ &=0. \end{aligned}

Repeatedly renormalizing an exact closed-system solution should be unnecessary. In numerical work, norm drift is instead a diagnostic of discretization, solver, or implementation error unless the effective dynamics is intentionally nonunitary.

The structural reason for preservation is canonical in Unitary Time Evolution.

Continuous spectral bases are commonly normalized by a Dirac delta:

⟨a∣a′⟩=δ(a−a′).\langle a|a'\rangle = \delta(a-a').

Position and momentum examples are

⟨x∣x′⟩=δ(x−x′),⟨p∣p′⟩=δ(p−p′).\langle x|x'\rangle = \delta(x-x'), \qquad \langle p|p'\rangle = \delta(p-p').

This is not unit normalization. The generalized ket ∣a⟩|a\rangle is not normally an element of the physical Hilbert space, and δ(0)\delta(0) is not a large finite norm. The delta is a distribution specifying how continuum basis objects pair under integrals.

The identity resolution

I=∫da ∣a⟩⟨a∣I = \int da\,|a\rangle\langle a|

allows a normalizable state to be written

∣ψ⟩=∫da ψ(a)∣a⟩,|\psi\rangle = \int da\,\psi(a)|a\rangle,

with ordinary state normalization

∫da ∣ψ(a)∣2=1.\int da\,|\psi(a)|^2=1.

The generalized basis kets are delta normalized; the coefficient function of a physical state is square normalized. Confusing these two levels leads to illegal expressions such as ∣δ(a−a0)∣2|\delta(a-a_0)|^2 as an ordinary probability density.

Delta normalization depends on the spectral label. Let b=f(a)b=f(a) be monotonic. Since

δ(a−a′)=∣dbda∣δ(b−b′),\delta(a-a') = \left|\frac{db}{da}\right| \delta(b-b'),

normalized generalized bases are related by

∣b⟩=∣dadb∣1/2∣a(b)⟩.|b\rangle = \left|\frac{da}{db}\right|^{1/2} |a(b)\rangle.

Their coefficient functions obey

ψb(b)=∣dadb∣1/2ψa(a(b)).\psi_b(b) = \left|\frac{da}{db}\right|^{1/2} \psi_a(a(b)).

Consequently,

∣ψb(b)∣2db=∣ψa(a)∣2da.|\psi_b(b)|^2db = |\psi_a(a)|^2da.

This square-root Jacobian is the continuum analogue of a basis normalization factor. It appears when changing from momentum to wave number, from momentum to energy on a fixed branch, or between other spectral coordinates.

When ff is not one-to-one, each branch or degeneracy label must be retained. For a free particle, the energy E=p2/(2m)E=p^2/(2m) does not distinguish +p+p from −p-p, so an energy representation requires an additional channel label.

Box normalization replaces a continuum by a finite region with specified boundary conditions. On a periodic interval of length LL,

ψn(x)=1Leiknx,kn=2πnL.\psi_n(x) = \frac1{\sqrt L} e^{ik_nx}, \qquad k_n=\frac{2\pi n}{L}.

These modes satisfy

∫0Ldx ψm(x)∗ψn(x)=δmn.\int_0^Ldx\, \psi_m(x)^*\psi_n(x) = \delta_{mn}.

The finite-volume basis is genuinely orthonormal within the periodic Hilbert space. It is not the same object as a delta-normalized plane wave on the full line.

With

pn=ℏkn,Δp=2πℏL,p_n=\hbar k_n, \qquad \Delta p=\frac{2\pi\hbar}{L},

the large-box relation between a continuum momentum amplitude and discrete coefficients is

cn≃Δp ϕ(pn).c_n \simeq \sqrt{\Delta p}\,\phi(p_n).

Then

∑n∣cn∣2⟶∫dp ∣ϕ(p)∣2.\sum_n|c_n|^2 \longrightarrow \int dp\,|\phi(p)|^2.

At the basis level,

∣pn⟩box≃Δp ∣pn⟩δ.|p_n\rangle_{\mathrm{box}} \simeq \sqrt{\Delta p}\,|p_n\rangle_{\delta}.

The factor converts a Kronecker-normalized discrete ket to a Dirac-delta-normalized continuum ket. Densities of states arise because sums carry one state per spacing Δp\Delta p.

Boundary Conditions Are Part of the Convention

Section titled “Boundary Conditions Are Part of the Convention”

The phrase “normalize in a box” is incomplete without boundary conditions. Periodic boundary conditions yield plane-wave momentum modes. Hard-wall Dirichlet conditions yield standing waves instead. Other self-adjoint boundary conditions can shift or reorganize the spectrum.

The factor 1/L1/\sqrt L normalizes a periodic plane wave because its modulus is constant on an interval of length LL. It does not imply that the same plane wave is a physical eigenstate for every finite interval problem.

Box normalization is useful for:

  • turning integrals into sums;
  • defining finite-volume numerical bases;
  • counting modes and deriving densities of states;
  • regulating continuum expressions;
  • connecting Kronecker and Dirac deltas.

The artificial volume should disappear from physical continuum predictions after sums, amplitudes, and state densities are converted consistently.

Scattering calculations often use stationary states that are not square normalizable. One may choose their amplitude so the incident probability current has a prescribed value, commonly unit flux. This is flux normalization, not total-probability normalization.

For a one-dimensional plane wave

ψ(x)=Aeikx,p=ℏk,\psi(x)=Ae^{ikx}, \qquad p=\hbar k,

the probability current is

j=ℏkm∣A∣2=pm∣A∣2.j = \frac{\hbar k}{m}|A|^2 = \frac{p}{m}|A|^2.

Choosing unit incident flux would require

∣A∣2=m∣p∣|A|^2=\frac{m}{|p|}

for the magnitude convention ∣j∣=1|j|=1, up to any additional continuum or channel normalization factors. The dimensions of such a mode differ from those of a unit-normalized bound-state wavefunction.

Flux-normalized states are designed so reflected, transmitted, or scattered flux ratios yield probabilities or cross sections. The detailed convention depends on dimension, channels, relativistic versus nonrelativistic kinematics, and the definition of the scattering matrix. Those applications belong to the scattering volume and Normalization Conventions.

Standard non-normalizable objects include:

  • exact position or momentum eigenkets;
  • plane waves on all of space;
  • energy-normalized scattering eigenfunctions;
  • an unregularized Dirac delta used as a putative state;
  • formal solutions with nonintegrable growth or tails.

They can still be useful as generalized basis vectors, asymptotic modes, Green function sources, or limits of normalized wave packets. Their usefulness does not make them ordinary physical states.

For example, a family of normalized Gaussians can become increasingly narrow in position, but its limit is not a normalized L2L^2 vector. The sequence approaches a delta distribution only in a weak sense, while its momentum spread diverges. There is no unit-norm position eigenvector hidden at the endpoint of the Hilbert space.

Similarly, a sequence of broader normalized wave packets can approximate a plane wave locally while spreading its total probability over an ever larger region. The limiting plane wave is delta or flux normalized, not square normalized.

Suppose a normalized position wavefunction is known to lie in a region Δ\Delta after a successful post-selection. The unnormalized projected wavefunction is

ψ~Δ(x)=1Δ(x)ψ(x),\widetilde\psi_\Delta(x) = \mathbf 1_\Delta(x)\psi(x),

where 1Δ\mathbf 1_\Delta is the indicator function. Its squared norm is

pΔ=∫Δdx ∣ψ(x)∣2.p_\Delta = \int_\Delta dx\,|\psi(x)|^2.

If pΔ>0p_\Delta>0, the conditional normalized wavefunction is

ψΔ(x)=1Δ(x)ψ(x)pΔ.\psi_\Delta(x) = \frac{\mathbf 1_\Delta(x)\psi(x)} {\sqrt{p_\Delta}}.

This renormalization changes the state because it conditions on new information or a measurement event. It should not be confused with merely choosing a unit representative of an unchanged ray.

If pΔ=0p_\Delta=0, no conditional state for that event is defined by this formula. Division by zero is not repaired by assigning an arbitrary vector to an impossible branch.

On a grid xjx_j with quadrature weights wjw_j, the continuum norm is approximated by

∥ψ∥2≈∑jwj∣ψj∣2.\|\psi\|^2 \approx \sum_j w_j|\psi_j|^2.

For a uniform grid with spacing Δx\Delta x,

∥ψ∥2≈Δx∑j∣ψj∣2.\|\psi\|^2 \approx \Delta x\sum_j|\psi_j|^2.

The raw Euclidean norm ∑j∣ψj∣2=1\sum_j|\psi_j|^2=1 is therefore not generally the continuum normalization unless the stored array has already absorbed Δx\sqrt{\Delta x}.

For a nonorthogonal finite basis with overlap matrix SS, a coefficient vector cc is normalized by

c†Sc=1.c^\dagger Sc=1.

In generalized eigenvalue problems,

Hc=ESc,Hc=ESc,

the same overlap matrix defines the physical inner product. Normalizing with c†c=1c^\dagger c=1 instead can produce basis-dependent errors.

Useful numerical checks include:

  1. include quadrature or overlap weights;
  2. verify that the norm is real and nonnegative within tolerance;
  3. reject zero or nearly zero vectors before division;
  4. track discarded weight after truncation;
  5. monitor norm drift during nominally unitary evolution;
  6. distinguish deliberate subnormalization from numerical loss;
  7. avoid hiding large errors by renormalizing at every step.

Normalization Does Not Guarantee Physical Admissibility

Section titled “Normalization Does Not Guarantee Physical Admissibility”

Unit norm is necessary for an ordinary state vector, but it is not sufficient for every physical question. A normalized wavefunction may fail to lie in the domain of an unbounded observable. For example, it may have finite norm but infinite kinetic-energy expectation or undefined boundary derivatives.

Likewise, a trace-one operator is not a state unless it is also positive. The matrix

(200−1)\begin{pmatrix} 2&0\\ 0&-1 \end{pmatrix}

has trace one but is not positive and therefore is not a density operator.

Normalization controls total weight. Positivity, domains, symmetry constraints, boundary conditions, and finite expectation values impose additional requirements depending on the system and observable.

Before interpreting amplitudes, ask:

  1. What object is being normalized? An ordinary vector, density operator, generalized eigenfunction, box mode, flux mode, or conditional branch?
  2. What inner product or measure is used? Is there a Jacobian, overlap matrix, spin sum, or quadrature weight?
  3. Is the norm finite and nonzero? Otherwise unit normalization is impossible.
  4. What is the target condition? Unit norm, unit trace, delta function, Kronecker delta, prescribed flux, or retained branch probability?
  5. Are labels continuous or discrete? This determines integrals versus sums and Dirac versus Kronecker deltas.
  6. Has a variable been relabelled? Include the square-root Jacobian in amplitudes.
  7. Is the state conditional? Preserve the branch weight before dividing by it.
  8. Will a later limit be taken? Track all box-volume and density-of-states factors before removing the regulator.
  • Forgetting to normalize an ordinary state before applying the simplest Born formulas.
  • Trying to normalize the zero vector or a vector of infinite norm.
  • Thinking normalization removes global phase or chooses one unique vector on a ray.
  • Using ∑i∣ci∣2\sum_i|c_i|^2 in a nonorthogonal basis instead of c†Gcc^\dagger Gc.
  • Dropping coordinate Jacobians, spin sums, particle coordinates, quadrature weights, or overlap matrices.
  • Calling ∣ψ(x)∣2|\psi(x)|^2 a probability instead of a density with respect to a measure.
  • Treating delta normalization as ordinary unit normalization.
  • Squaring a Dirac delta as though it were an ordinary function.
  • Mixing box-normalized, delta-normalized, and flux-normalized amplitudes.
  • Renormalizing a conditional branch before recording its probability.
  • Assuming trace one is enough for a density operator without checking positivity.
  • Renormalizing every numerical time step and thereby hiding nonunitary error.
  • Assuming any unit-norm wavefunction lies in the domain of every observable.

This page owns the meaning and translation of normalization conventions across the core formalism. Detailed neighboring topics remain canonical elsewhere:

  • P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press, 1958, Chapters II and III.
  • J. von Neumann, Mathematical Foundations of Quantum Mechanics, Princeton University Press, 1955, Chapters II and III.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2021, Chapters 1 and 2.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994, Chapters 1, 4, and 5.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014, Chapters 2 and 3.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013, Chapters 3, 7, and 10.

Let

∣b1⟩=∣0⟩,∣b2⟩=∣0⟩+∣1⟩2,|b_1\rangle=|0\rangle, \qquad |b_2\rangle = \frac{|0\rangle+|1\rangle}{\sqrt2},

where ∣0⟩,∣1⟩|0\rangle,|1\rangle are orthonormal. Normalize ∣χ⟩=∣b1⟩+∣b2⟩|\chi\rangle=|b_1\rangle+|b_2\rangle using the Gram matrix.

Solution

The Gram matrix is

G=(11/21/21),G = \begin{pmatrix} 1&1/\sqrt2\\ 1/\sqrt2&1 \end{pmatrix},

and the coordinate column is

c=(11).c = \begin{pmatrix} 1\\ 1 \end{pmatrix}.

Therefore

∥χ∥2=c†Gc=2+2.\|\chi\|^2 = c^\dagger Gc = 2+\sqrt2.

A normalized vector is

∣ψ⟩=∣b1⟩+∣b2⟩2+2.|\psi\rangle = \frac{|b_1\rangle+|b_2\rangle} {\sqrt{2+\sqrt2}}.

The naive coefficient sum would give 22 and would miss the nonzero overlap.

Let

ψ(x)=A1+x2/a2,a>0.\psi(x) = \frac{A}{1+x^2/a^2}, \qquad a>0.

Find ∣A∣|A| on the real line.

Solution

Set u=x/au=x/a. Then

1=∣A∣2a∫−∞∞du(1+u2)2=∣A∣2aπ2.\begin{aligned} 1 &= |A|^2a \int_{-\infty}^{\infty} \frac{du}{(1+u^2)^2} \\ &= |A|^2a\frac\pi2. \end{aligned}

Thus

∣A∣=2πa.|A| = \sqrt{\frac{2}{\pi a}}.

Its phase remains arbitrary.

Suppose

ψ(r)=R(r)Yℓm(θ,φ),\psi(\mathbf r) = R(r)Y_{\ell m}(\theta,\varphi),

with ∫dΩ ∣Yℓm∣2=1\int d\Omega\,|Y_{\ell m}|^2=1. If u(r)=rR(r)u(r)=rR(r), show that three-dimensional normalization is equivalent to ∫0∞dr ∣u(r)∣2=1\int_0^\infty dr\,|u(r)|^2=1.

Solution

Using d3r=r2dr dΩd^3r=r^2dr\,d\Omega,

∫d3r ∣ψ(r)∣2=(∫0∞dr r2∣R(r)∣2)×(∫dΩ ∣Yℓm∣2)=∫0∞dr r2∣R(r)∣2=∫0∞dr ∣u(r)∣2.\begin{aligned} \int d^3r\,|\psi(\mathbf r)|^2 &= \left( \int_0^\infty dr\,r^2|R(r)|^2 \right) \\ &\quad\times \left( \int d\Omega\,|Y_{\ell m}|^2 \right) \\ &= \int_0^\infty dr\,r^2|R(r)|^2 \\ &= \int_0^\infty dr\,|u(r)|^2. \end{aligned}

The factor rr in u=rRu=rR absorbs the radial Jacobian into the represented function.

Let a normalized state have orthonormal-basis coefficients cnc_n. Define

∣ψN⟩=∑n=1Ncn∣n⟩|\psi_N\rangle = \sum_{n=1}^{N}c_n|n\rangle

and wN=∑n=1N∣cn∣2>0w_N=\sum_{n=1}^{N}|c_n|^2>0. Normalize the truncation and find its fidelity with the original state.

Solution

The normalized truncation is

∣ψ^N⟩=1wN∑n=1Ncn∣n⟩.|\widehat\psi_N\rangle = \frac1{\sqrt{w_N}} \sum_{n=1}^{N}c_n|n\rangle.

Its overlap with the full state is

⟨ψ∣ψ^N⟩=1wN∑n=1N∣cn∣2=wN.\begin{aligned} \langle\psi|\widehat\psi_N\rangle &= \frac1{\sqrt{w_N}} \sum_{n=1}^{N}|c_n|^2 \\ &= \sqrt{w_N}. \end{aligned}

Therefore the pure-state fidelity is

∣⟨ψ∣ψ^N⟩∣2=wN.|\langle\psi|\widehat\psi_N\rangle|^2 = w_N.

Renormalization does not restore the discarded weight 1−wN1-w_N.

A normalized qubit is

∣ψ⟩=α∣0⟩+β∣1⟩.|\psi\rangle = \alpha|0\rangle+\beta|1\rangle.

Project onto ∣0⟩|0\rangle. Find the unnormalized branch, its norm, and the conditional state when the outcome has nonzero probability.

Solution

With P0=∣0⟩⟨0∣P_0=|0\rangle\langle0|,

∣ψ~0⟩=P0∣ψ⟩=α∣0⟩.|\widetilde\psi_0\rangle = P_0|\psi\rangle = \alpha|0\rangle.

Its squared norm is

⟨ψ~0∣ψ~0⟩=∣α∣2,\langle\widetilde\psi_0 |\widetilde\psi_0\rangle = |\alpha|^2,

which is the outcome probability. If α≠0\alpha\ne0, the conditional state is

∣ψ0⟩=α∣α∣∣0⟩,|\psi_0\rangle = \frac{\alpha}{|\alpha|}|0\rangle,

which represents the same ray as ∣0⟩|0\rangle. If α=0\alpha=0, the branch is impossible and the division is undefined.

Suppose E=p2/(2m)E=p^2/(2m) with p>0p>0, and ϕ(p)\phi(p) is normalized on the positive momentum half-line. Find the energy amplitude χ(E)\chi(E) normalized with respect to dEdE.

Solution

On the positive branch,

p(E)=2mE,dpdE=mp(E).p(E)=\sqrt{2mE}, \qquad \frac{dp}{dE} = \frac{m}{p(E)}.

Probability invariance requires

∣χ(E)∣2dE=∣ϕ(p)∣2dp.|\chi(E)|^2dE = |\phi(p)|^2dp.

Therefore

χ(E)=mp(E)ϕ(p(E))\chi(E) = \sqrt{\frac{m}{p(E)}} \phi\bigl(p(E)\bigr)

up to a phase convention. If both signs of momentum are present, energy alone is not a complete label; one must retain a direction or channel index.

7. Norm conservation from the Schrödinger equation

Section titled “7. Norm conservation from the Schrödinger equation”

Assume H=H†H=H^\dagger and that all domain conditions needed below hold. Starting from

iℏ∣ψ˙⟩=H∣ψ⟩,i\hbar|\dot\psi\rangle=H|\psi\rangle,

show that d⟨ψ∣ψ⟩/dt=0d\langle\psi|\psi\rangle/dt=0.

Solution

The Schrödinger equation and its adjoint give

∣ψ˙⟩=−iℏH∣ψ⟩,⟨ψ˙∣=iℏ⟨ψ∣H.|\dot\psi\rangle = -\frac{i}{\hbar}H|\psi\rangle, \qquad \langle\dot\psi| = \frac{i}{\hbar}\langle\psi|H.

Hence

ddt⟨ψ∣ψ⟩=⟨ψ˙∣ψ⟩+⟨ψ∣ψ˙⟩=iℏ⟨ψ∣H∣ψ⟩−iℏ⟨ψ∣H∣ψ⟩=0.\begin{aligned} \frac{d}{dt}\langle\psi|\psi\rangle &= \langle\dot\psi|\psi\rangle + \langle\psi|\dot\psi\rangle \\ &= \frac{i}{\hbar}\langle\psi|H|\psi\rangle - \frac{i}{\hbar}\langle\psi|H|\psi\rangle \\ &=0. \end{aligned}

Self-adjointness makes the bra equation use the same HH.

Consider

A=(200−1).A = \begin{pmatrix} 2&0\\ 0&-1 \end{pmatrix}.

Show that AA has unit trace but cannot be a density operator. Give an effect for which the Born expression is negative.

Solution

The trace is

Tr⁡A=2−1=1.\operatorname{Tr}A = 2-1 =1.

However, with the positive rank-one effect

E=∣1⟩⟨1∣,E=|1\rangle\langle1|,

one obtains

Tr⁡(AE)=⟨1∣A∣1⟩=−1.\operatorname{Tr}(AE) = \langle1|A|1\rangle = -1.

A probability cannot be negative. The failure is that AA is not positive, so trace normalization alone does not make it a state.