Skip to content

Wave-Mechanics Postulates

Wave mechanics is the position-representation form of nonrelativistic quantum mechanics. A pure state is represented by a wavefunction, observables act as multiplication or differential operators, position probabilities are integrals of a probability density, and closed-system dynamics is governed by the Schrödinger equation together with an operator domain and boundary conditions.

This representation is adapted to particles moving through space, so it is the practical language of wells, barriers, tunneling, wave packets, atoms, and scattering. It is not a theory separate from the Minimal Postulates. A wavefunction is a coordinate representation of a state, just as a column is a basis representation of a finite-dimensional state.

For one spinless particle in a configuration region Ω⊆Rd\Omega\subseteq\mathbb R^d:

  1. State: a pure state is a ray represented by ψ∈L2(Ω,dμ)\psi\in L^2(\Omega,d\mu).

  2. Normalization: choose a representative satisfying ∫Ω∣ψ∣2dμ=1\int_\Omega\lvert\psi\rvert^2d\mu=1.

  3. Observables: suitable self-adjoint operators act on wavefunctions, often through multiplication or differentiation on declared domains.

  4. Position Born rule: for a measurable region R⊆ΩR\subseteq\Omega,

    Pr⁡(r∈R)=∫R∣ψ(r)∣2dμ(r).\Pr(\mathbf r\in R) = \int_R \lvert\psi(\mathbf r)\rvert^2 d\mu(\mathbf r).
  5. Dynamics: a closed-system wavefunction satisfies

    iℏ∂ψ∂t=H^ψ.i\hbar \frac{\partial\psi}{\partial t} = \hat H\psi.
  6. Boundary conditions: the Hilbert space, operator domain, and boundary conditions complete the physical specification of differential operators.

  7. Composition: a multiparticle wavefunction lives on configuration space; product functions represent product states, while nonfactorizable functions represent entanglement.

The first six items are the traditional wave-mechanics package. The seventh makes explicit how the abstract composition postulate appears in coordinates.

The simplest setting is a single spinless particle with Cartesian coordinate r∈R3\mathbf r\in\mathbb R^3. Its position representation is

H≅L2(R3,d3r).\mathcal H \cong L^2(\mathbb R^3,d^3r).

If the particle is confined to a region Ω\Omega, one instead uses an appropriate L2(Ω,dμ)L^2(\Omega,d\mu) space and an operator domain encoding the boundary conditions. The measure dμd\mu is part of the representation.

This page assumes:

  • a nonrelativistic particle with an external time parameter;
  • a scalar wavefunction unless spin or internal components are stated;
  • standard complex Hilbert-space quantum mechanics;
  • sufficiently regular potentials for the displayed differential formulas;
  • normalized states when probabilities are quoted.

Spin introduces multicomponent wavefunctions, mixed states require density operators or integral kernels, identical particles require permutation symmetry, and variable particle number leads toward Fock space. Relativistic quantum field theory is not obtained by merely replacing the Schrödinger equation with a relativistic-looking one-particle equation.

Let ∣ψ⟩\lvert\psi\rangle be an abstract pure state. Its position-space wavefunction is the amplitude

ψ(r)=⟨r∣ψ⟩.\psi(\mathbf r) = \langle\mathbf r\mid\psi\rangle.

The symbols ∣r⟩\lvert\mathbf r\rangle are generalized position eigenkets, not normalizable vectors in L2L^2. They are distributional objects used to express the position representation. The physical state is the normalizable vector reconstructed from its amplitudes, not an exact-position ket.

Multiplying every amplitude by one constant phase does not change the pure state:

ψ(r)∼eiαψ(r).\psi(\mathbf r) \sim e^{i\alpha}\psi(\mathbf r).

Relative phase does matter. If

ψ=c1ψ1+c2ψ2,\psi = c_1\psi_1+c_2\psi_2,

then

∣ψ∣2=∣c1ψ1∣2+∣c2ψ2∣2+2Re⁡(c1c2∗ψ1ψ2∗).\begin{aligned} \lvert\psi\rvert^2 &= \lvert c_1\psi_1\rvert^2 + \lvert c_2\psi_2\rvert^2\\ &\quad+ 2\operatorname{Re} \left( c_1c_2^*\psi_1\psi_2^* \right). \end{aligned}

The last term is the interference term. A common overall phase cancels, whereas changing the phase between the two components changes this term.

An element of L2(Ω,dμ)L^2(\Omega,d\mu) is technically an equivalence class of functions. Two representatives define the same Hilbert-space vector if they differ only on a set of measure zero. Consequently, assigning a different value to a wavefunction at one isolated point does not create a new state or change a position probability.

This fact does not make pointwise regularity irrelevant. Differential operators require domains whose representatives possess suitable differentiability and boundary behavior. The Hilbert space tells us which states have finite norm; the operator domain tells us which of those states an unbounded operator may act on.

Not every state is represented by one wavefunction. A density operator can be represented by an integral kernel

ρ(r,r′)=⟨r∣ρ^∣r′⟩.\rho(\mathbf r,\mathbf r') = \langle\mathbf r\rvert \hat\rho \lvert\mathbf r'\rangle.

For a pure state,

ρψ(r,r′)=ψ(r)ψ(r′)∗.\rho_\psi(\mathbf r,\mathbf r') = \psi(\mathbf r)\psi(\mathbf r')^*.

The diagonal gives the position density,

ρψ(r,r)=∣ψ(r)∣2,\rho_\psi(\mathbf r,\mathbf r) = \lvert\psi(\mathbf r)\rvert^2,

but off-diagonal entries carry coherence information. The Density-Matrix Formulation is the canonical general-state presentation.

Normalization, Inner Products, and Measure

Section titled “Normalization, Inner Products, and Measure”

A normalized scalar wavefunction obeys

∫Ω∣ψ(r)∣2dμ(r)=1.\int_\Omega \lvert\psi(\mathbf r)\rvert^2 d\mu(\mathbf r) =1.

The inner product is

⟨ϕ∣ψ⟩=∫Ωϕ(r)∗ψ(r)dμ(r).\langle\phi\mid\psi\rangle = \int_\Omega \phi(\mathbf r)^* \psi(\mathbf r) d\mu(\mathbf r).

Orthogonality means this integral vanishes. Completeness of an orthonormal discrete basis {φn}\{\varphi_n\} means that a suitable state can be expanded as

ψ(r)=∑ncnφn(r),\psi(\mathbf r) = \sum_n c_n\varphi_n(\mathbf r),

with

cn=∫Ωφn(r)∗ψ(r)dμ(r).c_n = \int_\Omega \varphi_n(\mathbf r)^* \psi(\mathbf r) d\mu(\mathbf r).

Parseval’s identity then gives

∑n∣cn∣2=∫Ω∣ψ(r)∣2dμ(r).\sum_n\lvert c_n\rvert^2 = \int_\Omega \lvert\psi(\mathbf r)\rvert^2 d\mu(\mathbf r).

In Cartesian coordinates,

dμ(r)=d3r=dx dy dz.d\mu(\mathbf r) = d^3r = dx\,dy\,dz.

In spherical coordinates,

d3r=r2sin⁡θ dr dθ dϕ.d^3r = r^2\sin\theta \,dr\,d\theta\,d\phi.

More generally, a coordinate change r=r(q)\mathbf r=\mathbf r(\mathbf q) introduces the Jacobian:

ddr=∣det⁡∂r∂q∣ddq.d^dr = \left\lvert \det\frac{\partial\mathbf r}{\partial\mathbf q} \right\rvert d^dq.

The wavefunction and measure must be interpreted together. Omitting the Jacobian changes normalization, expectation values, and probabilities.

Because a probability is dimensionless, a Cartesian dd-dimensional wavefunction has units

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

Thus ∣ψ(r)∣2\lvert\psi(\mathbf r)\rvert^2 is a density with units L−d\mathsf L^{-d}, not a dimensionless probability. A probability appears only after integration over a region.

The site-wide normalization convention is summarized in Wavefunction Normalization.

Postulate 2: Operators in Position Representation

Section titled “Postulate 2: Operators in Position Representation”

An abstract operator becomes a rule acting on wavefunctions. In Cartesian position representation,

(r^ψ)(r)=r ψ(r),(\hat{\mathbf r}\psi)(\mathbf r) = \mathbf r\,\psi(\mathbf r),

while momentum acts formally as

(p^ψ)(r)=−iℏ∇ψ(r).(\hat{\mathbf p}\psi)(\mathbf r) = -i\hbar\nabla\psi(\mathbf r).

For a spinless particle in a scalar potential,

H^=p^ 22m+V(r^,t)\hat H = \frac{\hat{\mathbf p}^{\,2}}{2m} + V(\hat{\mathbf r},t)

becomes

(H^ψ)(r)=−ℏ22m∇2ψ(r)+V(r,t)ψ(r).(\hat H\psi)(\mathbf r) = -\frac{\hbar^2}{2m} \nabla^2\psi(\mathbf r) + V(\mathbf r,t)\psi(\mathbf r).

An expectation value is

⟨A⟩ψ=∫Ωψ(r)∗(A^ψ)(r)dμ(r),\langle A\rangle_\psi = \int_\Omega \psi(\mathbf r)^* (\hat A\psi)(\mathbf r) d\mu(\mathbf r),

provided ψ\psi belongs to the required domain and the integral exists.

The symbols −iℏ d/dx-i\hbar\,d/dx and −ℏ2d2/(2m dx2)-\hbar^2d^2/(2m\,dx^2) are differential expressions. A physical operator also requires:

  1. a Hilbert space;
  2. a domain of functions on which it acts;
  3. boundary or asymptotic conditions;
  4. a specification of how singular points are handled.

The same differential expression can define different operators with different spectra. On a line, a box, a ring, and a half-line, the kinetic-energy expression is accompanied by different spaces and domains.

In finite dimensions, Hermitian matrices are automatically self-adjoint. For unbounded differential operators, symmetry of matrix elements on a convenient set is not enough; the operator and its adjoint must have the same domain. See Hermitian vs Self-Adjoint Operators for the canonical distinction.

Postulate 3: The Born Rule in Position Space

Section titled “Postulate 3: The Born Rule in Position Space”

For a normalized pure state, the probability of detecting the particle in a measurable region R⊆ΩR\subseteq\Omega is

Pr⁡(r∈R)=∫R∣ψ(r)∣2dμ(r).\Pr(\mathbf r\in R) = \int_R \lvert\psi(\mathbf r)\rvert^2 d\mu(\mathbf r).

The associated spectral projector is formally

P^R=∫R∣r⟩⟨r∣dμ(r),\hat P_R = \int_R \lvert\mathbf r\rangle \langle\mathbf r\rvert d\mu(\mathbf r),

so the abstract Born rule gives

⟨ψ∣P^R∣ψ⟩=∫R∣ψ(r)∣2dμ(r).\langle\psi\rvert \hat P_R \lvert\psi\rangle = \int_R \lvert\psi(\mathbf r)\rvert^2 d\mu(\mathbf r).

The density at a point is not the probability of that exact point. For a continuous distribution,

Pr⁡(r=r0)=0\Pr(\mathbf r=\mathbf r_0) =0

under ordinary conditions, while

Pr⁡(r∈ddr)≈∣ψ(r)∣2ddr\Pr(\mathbf r\in d^dr) \approx \lvert\psi(\mathbf r)\rvert^2d^dr

describes a sufficiently small detector cell. Real detectors have finite spatial resolution and are more accurately represented by effects associated with finite regions or response functions.

The position-space wavefunction does not privilege position as the only observable. If A^\hat A has a normalized discrete eigenfunction φa(r)\varphi_a(\mathbf r), then

Pr⁡(a)=∣∫Ωφa(r)∗ψ(r)dμ(r)∣2.\Pr(a) = \left\lvert \int_\Omega \varphi_a(\mathbf r)^* \psi(\mathbf r) d\mu(\mathbf r) \right\rvert^2.

For momentum in dd Cartesian dimensions, use the Fourier convention

ψ~(p)=1(2πℏ)d/2∫Rde−ip⋅r/ℏψ(r)ddr.\widetilde\psi(\mathbf p) = \frac{1} {(2\pi\hbar)^{d/2}} \int_{\mathbb R^d} e^{-i\mathbf p\cdot\mathbf r/\hbar} \psi(\mathbf r) d^dr.

Then

∫Rd∣ψ~(p)∣2ddp=1,\int_{\mathbb R^d} \lvert\widetilde\psi(\mathbf p)\rvert^2 d^dp =1,

and momentum probabilities are obtained from ∣ψ~(p)∣2\lvert\widetilde\psi(\mathbf p)\rvert^2. The functions ψ(r)\psi(\mathbf r) and ψ~(p)\widetilde\psi(\mathbf p) represent the same state in different generalized bases. See Momentum Representation and Born Rule for Continuous Spectra.

For a closed system, the abstract equation

iℏddt∣ψ(t)⟩=H^(t)∣ψ(t)⟩i\hbar \frac{d}{dt} \lvert\psi(t)\rangle = \hat H(t)\lvert\psi(t)\rangle

becomes

iℏ∂∂tψ(r,t)=H^(t)ψ(r,t).i\hbar \frac{\partial}{\partial t} \psi(\mathbf r,t) = \hat H(t)\psi(\mathbf r,t).

For the standard scalar Hamiltonian,

iℏ∂ψ∂t=−ℏ22m∇2ψ+V(r,t)ψ.i\hbar \frac{\partial\psi}{\partial t} = -\frac{\hbar^2}{2m} \nabla^2\psi + V(\mathbf r,t)\psi.

The initial wavefunction and Hamiltonian domain determine an initial-value problem:

ψ(r,t0)=ψ0(r).\psi(\mathbf r,t_0) = \psi_0(\mathbf r).

For a time-independent self-adjoint Hamiltonian with a discrete orthonormal eigenbasis,

H^φn=Enφn,\hat H\varphi_n = E_n\varphi_n,

an initial state expands as

ψ(r,0)=∑ncnφn(r),\psi(\mathbf r,0) = \sum_n c_n\varphi_n(\mathbf r),

and evolves as

ψ(r,t)=∑ncne−iEnt/ℏφn(r).\psi(\mathbf r,t) = \sum_n c_n e^{-iE_nt/\hbar} \varphi_n(\mathbf r).

Each energy component acquires a phase. Relative phases change in time and can alter interference, even though the coefficient magnitudes ∣cn∣\lvert c_n\rvert remain fixed.

The Schrödinger Equation page owns the abstract dynamics, while Time-Dependent Schrödinger Equation develops coordinate-space solution methods.

For a real scalar potential, define

ϱ(r,t)=∣ψ(r,t)∣2\varrho(\mathbf r,t) = \lvert\psi(\mathbf r,t)\rvert^2

and

j(r,t)=ℏmIm⁡(ψ∗∇ψ).\mathbf j(\mathbf r,t) = \frac{\hbar}{m} \operatorname{Im} \left( \psi^*\nabla\psi \right).

Combining the Schrödinger equation with its complex conjugate yields the local continuity equation

∂ϱ∂t+∇⋅j=0.\frac{\partial\varrho}{\partial t} + \nabla\cdot\mathbf j =0.

Integrating over a fixed region RR and applying the divergence theorem gives

ddt∫Rϱ ddr=−∫∂Rj⋅dS.\frac{d}{dt} \int_R\varrho\,d^dr = -\int_{\partial R} \mathbf j\cdot d\mathbf S.

Probability in a region changes only through flux across its boundary. For the whole configuration space, boundary behavior that makes the net outward flux vanish gives

ddt∫Ω∣ψ∣2dμ=0.\frac{d}{dt} \int_\Omega \lvert\psi\rvert^2d\mu =0.

This is the position-space expression of norm preservation under unitary evolution. A complex absorbing potential or an explicitly non-Hermitian effective Hamiltonian can add source or sink terms; such dynamics is not closed-system unitary evolution on the stated Hilbert space.

The full derivations and interpretation are in Continuity Equation and Probability Current.

Boundary Conditions Complete the Hamiltonian

Section titled “Boundary Conditions Complete the Hamiltonian”

Consider the one-dimensional kinetic-energy expression

T^=−ℏ22md2dx2.\hat T = -\frac{\hbar^2}{2m} \frac{d^2}{dx^2}.

Integration by parts on an interval [a,b][a,b] gives the boundary form

⟨ϕ∣T^ψ⟩−⟨T^ϕ∣ψ⟩=−ℏ22m[ϕ∗ψ′−(ϕ′)∗ψ]ab.\begin{aligned} &\langle\phi\mid\hat T\psi\rangle - \langle\hat T\phi\mid\psi\rangle\\ &\qquad= -\frac{\hbar^2}{2m} \Big[ \phi^*\psi' - (\phi')^*\psi \Big]_a^b. \end{aligned}

A self-adjoint realization requires a domain whose boundary data make the appropriate boundary form vanish and whose adjoint has the same domain. Different admissible conditions describe different physical systems.

Common examples include:

  • decay or scattering asymptotics on the full line;

  • Dirichlet conditions ψ(0)=ψ(L)=0\psi(0)=\psi(L)=0 for an ideal hard-wall box;

  • periodic conditions

    ψ(L)=ψ(0),ψ′(L)=ψ′(0)\psi(L)=\psi(0), \qquad \psi'(L)=\psi'(0)

    for a ring without a boundary twist;

  • Robin conditions, which relate a value and normal derivative at a boundary;

  • continuity of ψ\psi and ψ′\psi' across finite nonsingular potential steps;

  • a derivative jump fixed by the coupling at a delta-function potential.

Boundary conditions are not decorative constraints added after solving the differential equation. They help define the operator, determine its spectrum, and enforce the intended probability flux.

The canonical practical treatment is Boundary Conditions.

For two distinguishable particles,

HAB=L2(ΩA)⊗L2(ΩB)≅L2(ΩA×ΩB).\begin{aligned} \mathcal H_{AB} &= L^2(\Omega_A) \otimes L^2(\Omega_B)\\ &\cong L^2(\Omega_A\times\Omega_B). \end{aligned}

The joint wavefunction is

Ψ(rA,rB)=⟨rA,rB∣Ψ⟩.\Psi(\mathbf r_A,\mathbf r_B) = \langle \mathbf r_A,\mathbf r_B \mid\Psi\rangle.

A product state factorizes:

Ψ(rA,rB)=ψA(rA)ϕB(rB).\Psi(\mathbf r_A,\mathbf r_B) = \psi_A(\mathbf r_A) \phi_B(\mathbf r_B).

Most joint wavefunctions do not factorize. A nonfactorizable pure-state wavefunction is entangled relative to the declared subsystem split.

Writing dμAB=dμA dμBd\mu_{AB}=d\mu_A\,d\mu_B, the probability that AA lies in RAR_A and BB lies in RBR_B is

Pr⁡(RA,RB)=∫RA×RB∣Ψ∣2dμAB.\Pr(R_A,R_B) = \int_{R_A\times R_B} \lvert\Psi\rvert^2 d\mu_{AB}.

For identical particles, physical wavefunctions are symmetric or antisymmetric under exchange:

Ψ(r1,r2)=±Ψ(r2,r1),\Psi(\mathbf r_1,\mathbf r_2) = \pm \Psi(\mathbf r_2,\mathbf r_1),

with the plus sign for bosons and the minus sign for fermions when the displayed variables include all relevant one-particle degrees of freedom.

Spin adds discrete indices. A spin-1/21/2 position-space state, for example, is a spinor

ψ(r)=(ψ↑(r)ψ↓(r)),\psi(\mathbf r) = \begin{pmatrix} \psi_\uparrow(\mathbf r)\\ \psi_\downarrow(\mathbf r) \end{pmatrix},

normalized by summing component densities:

∫d3r(∣ψ↑∣2+∣ψ↓∣2)=1.\int d^3r \left( \lvert\psi_\uparrow\rvert^2 + \lvert\psi_\downarrow\rvert^2 \right) =1.

These are coordinate forms of tensor products, not exceptions to the abstract composition postulate.

The wave-mechanics statements translate the abstract formalism as follows:

  • A ray ∣ψ⟩\lvert\psi\rangle becomes a wavefunction ψ(r)=⟨r∣ψ⟩\psi(\mathbf r)=\langle\mathbf r\mid\psi\rangle.
  • The abstract inner product becomes an integral against the coordinate measure.
  • A self-adjoint observable becomes a multiplication, differential, or integral operator with a specified domain.
  • A position-region projector becomes multiplication by the indicator of that region, expressed formally through generalized position kets.
  • The Born trace or expectation rule becomes an integral of a probability density.
  • Unitary evolution becomes the time-dependent Schrödinger partial differential equation.
  • Tensor products become functions on product configuration spaces, possibly with spinor indices or exchange symmetry.

Changing representation does not change predictions. For example,

⟨ϕ∣ψ⟩=∫ddr ϕ(r)∗ψ(r)\langle\phi\mid\psi\rangle = \int d^dr\, \phi(\mathbf r)^* \psi(\mathbf r)

and, with a unitary Fourier convention,

⟨ϕ∣ψ⟩=∫ddp ϕ~(p)∗ψ~(p).\langle\phi\mid\psi\rangle = \int d^dp\, \widetilde\phi(\mathbf p)^* \widetilde\psi(\mathbf p).

The coordinate and momentum wavefunctions are different descriptions of the same Hilbert-space vectors.

Worked Example: The Infinite Square Well as a Package

Section titled “Worked Example: The Infinite Square Well as a Package”

Consider one particle confined to 0<x<L0<x<L by ideal hard walls. The Hilbert space is

H=L2(0,L),\mathcal H = L^2(0,L),

and the Hamiltonian acts as

H^=−ℏ22md2dx2\hat H = -\frac{\hbar^2}{2m} \frac{d^2}{dx^2}

on a domain with Dirichlet conditions

ψ(0)=ψ(L)=0.\psi(0)=\psi(L)=0.

The normalized stationary eigenfunctions and energies are

φn(x)=2Lsin⁡nπxL,\varphi_n(x) = \sqrt{\frac{2}{L}} \sin\frac{n\pi x}{L},

and

En=n2π2ℏ22mL2,n=1,2,….E_n = \frac{n^2\pi^2\hbar^2} {2mL^2}, \qquad n=1,2,\ldots.

An arbitrary normalized state in the discrete span has

ψ(x,0)=∑n=1∞cnφn(x),∑n=1∞∣cn∣2=1,\psi(x,0) = \sum_{n=1}^{\infty} c_n\varphi_n(x), \qquad \sum_{n=1}^{\infty} \lvert c_n\rvert^2 =1,

and evolves as

ψ(x,t)=∑n=1∞cne−iEnt/ℏφn(x).\psi(x,t) = \sum_{n=1}^{\infty} c_n e^{-iE_nt/\hbar} \varphi_n(x).

For the ground state,

ψ1(x,t)=2Lsin⁡πxLe−iE1t/ℏ.\psi_1(x,t) = \sqrt{\frac{2}{L}} \sin\frac{\pi x}{L} e^{-iE_1t/\hbar}.

Its position density is time independent because the phase has unit magnitude:

∣ψ1(x,t)∣2=2Lsin⁡2πxL.\lvert\psi_1(x,t)\rvert^2 = \frac{2}{L} \sin^2\frac{\pi x}{L}.

The Born probability for the left half is

Pr⁡(0<x<L/2)=2L∫0L/2sin⁡2πxL dx=12.\begin{aligned} \Pr(0<x<L/2) &= \frac{2}{L} \int_0^{L/2} \sin^2\frac{\pi x}{L}\,dx\\ &= \frac12. \end{aligned}

The current vanishes because the spatial eigenfunction can be chosen real:

j(x,t)=ℏmIm⁡(ψ1∗∂ψ1∂x)=0.j(x,t) = \frac{\hbar}{m} \operatorname{Im} \left( \psi_1^* \frac{\partial\psi_1}{\partial x} \right) =0.

This one model uses the Hilbert space, operator domain, boundary conditions, normalization, spectrum, unitary time dependence, position Born rule, and current check. The ideal walls are a model assumption; a finite physical barrier requires a different Hamiltonian and matching problem.

For a new problem:

  1. Declare configuration space. Is the particle on a line, interval, ring, plane, or three-dimensional region?
  2. Declare the measure and components. Include coordinate Jacobians and spinor indices.
  3. Specify the Hamiltonian as an operator. State its differential expression, domain, and boundary or asymptotic conditions.
  4. Choose and normalize the initial state.
  5. Solve the spectral or initial-value problem.
  6. Represent the measurement. For position, specify a detector region or response; for another observable, use its spectral data or effects.
  7. Compute probabilities and expectations.
  8. Check normalization, dimensions, boundary flux, and limiting cases.

The Hamiltonian is physical model input. The postulates explain how a stated Hamiltonian generates predictions; they do not derive its potential or boundary conditions from notation alone.

Verify

∫Ω∣ψ(r,t)∣2dμ=1\int_\Omega \lvert\psi(\mathbf r,t)\rvert^2 d\mu =1

at the initial time and after evolution. For exact self-adjoint evolution, later normalization follows theoretically; in a numerical calculation, drift diagnoses discretization or time-stepping error.

Check that ∣ψ∣2dμ\lvert\psi\rvert^2d\mu and every reported probability are dimensionless. Fourier conventions give position and momentum wavefunctions different physical units.

Confirm that the stated boundary conditions support the claimed self-adjoint operator and that the probability flux matches the physical boundary. A hard wall should not leak.

Check free-particle, weak-potential, large-box, short-time, or symmetry limits as appropriate. A formally normalized answer can still solve the wrong boundary-value problem.

When practical, compute a quantity in more than one representation. Position-space integration and expansion in an orthonormal eigenbasis should agree.

A single scalar ψ(r,t)\psi(\mathbf r,t) is insufficient for:

  • spin and other internal components;
  • mixed preparations and reduced subsystem states;
  • several particles without a configuration-space wavefunction;
  • identical-particle exchange symmetry unless it is imposed;
  • variable particle number;
  • relativistic quantum field theory;
  • general measurements whose detector response is not a sharp position region.

These cases extend the representation; they do not invalidate the abstract postulates. The danger is not using wave mechanics, but mistaking its simplest scalar form for the whole theory.

  1. Treating ψ\psi as a classical material wave. It is a complex probability-amplitude representation.
  2. Calling ∣ψ(r)∣2\lvert\psi(\mathbf r)\rvert^2 a probability. It is a density; integrate it over a region.
  3. Forgetting the coordinate measure. Jacobian factors are part of normalization and expectation values.
  4. Assigning physical meaning to one pointwise value of an L2L^2 representative. States are defined almost everywhere.
  5. Using generalized eigenfunctions as normalized states. Position kets and plane waves require distributional normalization or wave packets.
  6. Writing a differential expression without a domain. Boundary conditions can change the operator and spectrum.
  7. Confusing symmetric with self-adjoint. Domain equality matters for unbounded operators.
  8. Dropping relative phases. Global phase is irrelevant; relative phase controls interference and current.
  9. Assuming stationary means time-independent wavefunction. A stationary energy eigenstate carries a time-dependent phase while its density is time independent.
  10. Using one scalar wavefunction for every system. Spin, mixed states, many particles, and fields require richer structures.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018 — standard wave-mechanics development and examples.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994 — Hilbert-space foundations connected carefully to coordinate representations.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020 — abstract and wave-mechanics formulations with symmetries and spin.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Vol. 1, Wiley, 1977 — detailed wavefunction, operator, and measurement treatment.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013 — rigorous bridge between Hilbert-space quantum theory and differential operators.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics, Vol. I: Functional Analysis, revised and enlarged ed., Academic Press, 1980 — domains, self-adjointness, and spectral foundations.
  1. Normalize a Gaussian. In one dimension, let
ψ(x)=Aexp⁡ ⁣(−x24σ2),σ>0.\psi(x) = A \exp\!\left( -\frac{x^2}{4\sigma^2} \right), \qquad \sigma>0.

Find ∣A∣\lvert A\rvert and state its physical units.

Solution

The density is

∣ψ(x)∣2=∣A∣2exp⁡ ⁣(−x22σ2).\lvert\psi(x)\rvert^2 = \lvert A\rvert^2 \exp\!\left( -\frac{x^2}{2\sigma^2} \right).

Using the Gaussian integral,

∫−∞∞e−x2/(2σ2)dx=2π σ.\int_{-\infty}^{\infty} e^{-x^2/(2\sigma^2)}dx = \sqrt{2\pi}\,\sigma.

Normalization gives

∣A∣22π σ=1,\lvert A\rvert^2 \sqrt{2\pi}\,\sigma =1,

so

∣A∣=1(2πσ2)1/4.\lvert A\rvert = \frac{1} {(2\pi\sigma^2)^{1/4}}.

In one dimension, AA and ψ\psi have units L−1/2\mathsf L^{-1/2}. The phase of AA is a global phase and may be chosen conveniently.

  1. Use the spherical measure. A normalized spherically symmetric state has wavefunction ψ(r)=Ce−r/a\psi(\mathbf r)=Ce^{-r/a} with a>0a>0. Find ∣C∣\lvert C\rvert.
Solution

Spherical symmetry gives

1=4π∣C∣2∫0∞r2e−2r/a dr.1 = 4\pi\lvert C\rvert^2 \int_0^\infty r^2e^{-2r/a}\,dr.

Using

∫0∞r2e−βr dr=2β3\int_0^\infty r^2e^{-\beta r}\,dr = \frac{2}{\beta^3}

with β=2/a\beta=2/a gives

∫0∞r2e−2r/a dr=a34.\int_0^\infty r^2e^{-2r/a}\,dr = \frac{a^3}{4}.

Therefore,

πa3∣C∣2=1,\pi a^3\lvert C\rvert^2 =1,

and

∣C∣=1πa3.\lvert C\rvert = \frac{1}{\sqrt{\pi a^3}}.

Omitting the factor 4πr24\pi r^2 would give the wrong normalization and wrong units.

  1. Global and relative phase. Compare
ψα=eiαψ1+ψ22\psi_\alpha = e^{i\alpha} \frac{\psi_1+\psi_2}{\sqrt2}

with

ψβ=ψ1+eiβψ22,\psi_\beta = \frac{\psi_1+e^{i\beta}\psi_2}{\sqrt2},

where ψ1\psi_1 and ψ2\psi_2 overlap. Which phase can affect the position density?

Solution

The common factor in ψα\psi_\alpha cancels:

∣ψα∣2=12∣ψ1+ψ2∣2.\lvert\psi_\alpha\rvert^2 = \frac12 \lvert\psi_1+\psi_2\rvert^2.

It is a global phase. For ψβ\psi_\beta,

∣ψβ∣2=12(∣ψ1∣2+∣ψ2∣2)+Re⁡(e−iβψ1ψ2∗).\begin{aligned} \lvert\psi_\beta\rvert^2 &= \frac12 \Big( \lvert\psi_1\rvert^2 + \lvert\psi_2\rvert^2 \Big)\\ &\quad+ \operatorname{Re} \left( e^{-i\beta}\psi_1\psi_2^* \right). \end{aligned}

The relative phase β\beta changes the interference term wherever the components overlap.

  1. A boundary-form check. On [0,L][0,L], show that the kinetic-energy boundary form vanishes for two twice-differentiable functions satisfying Dirichlet conditions at both endpoints.
Solution

The boundary form is proportional to

[ϕ∗ψ′−(ϕ′)∗ψ]0L.\Big[ \phi^*\psi' - (\phi')^*\psi \Big]_0^L.

Dirichlet conditions give

ϕ(0)=ϕ(L)=0,ψ(0)=ψ(L)=0.\begin{aligned} \phi(0)&=\phi(L)=0,\\ \psi(0)&=\psi(L)=0. \end{aligned}

At each endpoint, the first term vanishes because ϕ=0\phi=0 and the second because ψ=0\psi=0. Hence the entire boundary form vanishes. This establishes symmetry on the stated domain; a full self-adjointness claim also requires comparing the adjoint domain.

  1. Periodic quantization. A free particle on a ring of circumference LL has ψ(x+L)=ψ(x)\psi(x+L)=\psi(x). Starting from ψk(x)=Aeikx\psi_k(x)=Ae^{ikx}, find the allowed kk, normalize the wavefunction, and give the energies.
Solution

Periodicity requires

eikL=1,e^{ikL}=1,

so

kn=2πnL,n∈Z.k_n = \frac{2\pi n}{L}, \qquad n\in\mathbb Z.

Normalization on one period gives

1=∣A∣2L,1 = \lvert A\rvert^2L,

so one may choose A=1/LA=1/\sqrt L. The energies are

En=ℏ2kn22m=2π2ℏ2n2mL2.E_n = \frac{\hbar^2k_n^2}{2m} = \frac{2\pi^2\hbar^2n^2} {mL^2}.

Unlike a plane wave on the full line, the periodic wavefunction is normalizable on its finite configuration space.

  1. Continuity equation in one dimension. For real V(x,t)V(x,t), start from the Schrödinger equation and its complex conjugate and show
∂∂t∣ψ∣2+∂j∂x=0,\frac{\partial}{\partial t} \lvert\psi\rvert^2 + \frac{\partial j}{\partial x} =0,

where

j=ℏmIm⁡(ψ∗∂ψ∂x).j = \frac{\hbar}{m} \operatorname{Im} \left( \psi^* \frac{\partial\psi}{\partial x} \right).
Solution

Multiply the Schrödinger equation by ψ∗\psi^* and its conjugate by ψ\psi, then subtract. The real potential terms cancel:

iℏ∂∂t∣ψ∣2=−ℏ22m(ψ∗ψ′′−ψ(ψ∗)′′).\begin{aligned} i\hbar \frac{\partial}{\partial t} \lvert\psi\rvert^2 &= -\frac{\hbar^2}{2m} \left( \psi^*\psi'' - \psi(\psi^*)'' \right). \end{aligned}

The spatial expression is a total derivative:

ψ∗ψ′′−ψ(ψ∗)′′=∂∂x(ψ∗ψ′−ψ(ψ∗)′).\psi^*\psi'' - \psi(\psi^*)'' = \frac{\partial}{\partial x} \left( \psi^*\psi' - \psi(\psi^*)' \right).

Dividing by iℏi\hbar and using

j=ℏ2mi(ψ∗ψ′−ψ(ψ∗)′)j = \frac{\hbar}{2mi} \left( \psi^*\psi' - \psi(\psi^*)' \right)

gives the continuity equation.

  1. A two-state box superposition. In the infinite square well, let
ψ(x,0)=φ1(x)+φ2(x)2.\psi(x,0) = \frac{ \varphi_1(x)+\varphi_2(x) }{\sqrt2}.

Write ψ(x,t)\psi(x,t) and identify the angular frequency of the interference term in ∣ψ(x,t)∣2\lvert\psi(x,t)\rvert^2.

Solution

Each eigenstate acquires its own phase:

ψ(x,t)=12[e−iE1t/ℏφ1(x)+e−iE2t/ℏφ2(x)].\begin{aligned} \psi(x,t) &= \frac1{\sqrt2} \Big[ e^{-iE_1t/\hbar}\varphi_1(x)\\ &\qquad+ e^{-iE_2t/\hbar}\varphi_2(x) \Big]. \end{aligned}

Because the box eigenfunctions can be chosen real,

∣ψ(x,t)∣2=12[φ1(x)2+φ2(x)2]+φ1(x)φ2(x)×cos⁡ ⁣(E2−E1ℏt).\begin{aligned} \lvert\psi(x,t)\rvert^2 &= \frac12 \left[ \varphi_1(x)^2+\varphi_2(x)^2 \right]\\ &\quad+ \varphi_1(x)\varphi_2(x)\\ &\qquad\times \cos\!\left( \frac{E_2-E_1}{\hbar}t \right). \end{aligned}

The interference oscillates with angular frequency

ω21=E2−E1ℏ=3π2ℏ2mL2.\omega_{21} = \frac{E_2-E_1}{\hbar} = \frac{3\pi^2\hbar} {2mL^2}.

The total norm remains one even though the spatial density changes.

  1. Factorization and entanglement. Let ϕ0\phi_0 and ϕ1\phi_1 be orthonormal one-particle wavefunctions. Show that
Ψ(xA,xB)=12[ϕ0(xA)ϕ0(xB)+ϕ1(xA)ϕ1(xB)].\begin{aligned} \Psi(x_A,x_B) &= \frac1{\sqrt2} \Big[ \phi_0(x_A)\phi_0(x_B)\\ &\qquad+ \phi_1(x_A)\phi_1(x_B) \Big]. \end{aligned}

cannot be written as f(xA)g(xB)f(x_A)g(x_B).

Solution

Suppose a product form existed. Expand

f=aϕ0+bϕ1+⋯ ,g=cϕ0+dϕ1+⋯ .\begin{aligned} f&=a\phi_0+b\phi_1+\cdots,\\ g&=c\phi_0+d\phi_1+\cdots. \end{aligned}

Matching the four coefficients in the displayed two-dimensional subspace would require

ac=12,bd=12,ac=\frac1{\sqrt2}, \qquad bd=\frac1{\sqrt2},

and

ad=0,bc=0.ad=0, \qquad bc=0.

The first two equations require a,b,c,da,b,c,d all to be nonzero, contradicting the last two. Hence the wavefunction is nonfactorizable and the two particles are entangled relative to the AA–BB split.