Wave function and probability interpretation

A tossed coin has a definite face-up outcome the instant it lands, but before that you can only speak of odds. Quantum mechanics says a particle is stranger still: even in principle there is no hidden trajectory waiting to be revealed, only a wave function ψ(x,t)\psi(x,t) whose square tells you the odds of finding the particle here or there. This single idea — Born's probability interpretation — replaces Newton's sharp orbits with a cloud of possibility, and it is the object every later quantum calculation manipulates. This node builds the toolkit around ψ\psi: what it means, how to normalize it so probabilities add to one, how to extract measurable averages like ⟨x⟩\langle x\rangle and ⟨p⟩\langle p\rangle from it, and what "well-behaved" requires of a physically acceptable wave function. Everything here stays within one dimension and uses only integrals you already know from calculus — no Schrödinger equation yet, that is the next chapter's job.

The core insight

The wave function ψ(x,t)\psi(x,t) is not itself observable, but ∣ψ(x,t)∣2 dx|\psi(x,t)|^2\,dx is the probability of finding the particle in [x,x+dx][x,x+dx], so a normalized ψ\psi (with ∫∣ψ∣2dx=1\int|\psi|^2dx=1) lets every measurable average be computed by sandwiching the right operator between ψ∗\psi^* and ψ\psi.

The core idea

From trajectory to probability cloud. Classical mechanics tracks a particle along one definite path x(t)x(t). Quantum mechanics abandons that picture entirely: all the physical information about a particle in one dimension is packed into a single complex-valued function ψ(x,t)\psi(x,t), and ψ\psi itself carries no direct physical meaning — you cannot measure it on a dial. What you can measure are probabilities, and Born's rule is the bridge between the mathematics and the lab.

Born's rule and probability density. The rule states that the squared magnitude of ψ\psi gives a probability density in space:

∣ψ(x,t)∣2 dx=probability of finding the particle in the interval [x, x+dx] at time t.|\psi(x,t)|^2\,dx = \text{probability of finding the particle in the interval } [x,\,x+dx] \text{ at time } t.

where: ψ(x,t)\psi(x,t) is the (generally complex) wave function, ∣ψ∣2=ψ∗ψ|\psi|^2=\psi^*\psi is real and never negative, and dxdx is an infinitesimal length so that ∣ψ∣2|\psi|^2 itself has units of probability per unit length — a density, not a probability. To get an actual probability you must integrate ∣ψ∣2|\psi|^2 over a finite stretch of the line.

Normalization. Because the particle is somewhere on the line with certainty, the probabilities over all space must sum to exactly 1:

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

where: the integral runs over the whole region the particle can occupy. A ψ\psi obeying this is called normalized. If a solution of the wave equation instead integrates to some finite N≠1N\ne 1, rescaling ψnorm=ψ/N\psi_{\text{norm}}=\psi/\sqrt N restores normalization without changing the physics, since ∣ψnorm∣2=∣ψ∣2/N|\psi_{\text{norm}}|^2 = |\psi|^2/N integrates to exactly 1.

Expectation values. Once ψ\psi is normalized, every measurable quantity is obtained by an "operator sandwich" — the appropriate operator placed between ψ∗\psi^* and ψ\psi, then integrated over all space. Position is represented by plain multiplication by xx, so

⟨x⟩=∫−∞∞ψ∗(x,t) x ψ(x,t) dx,\langle x\rangle = \int_{-\infty}^{\infty}\psi^*(x,t)\,x\,\psi(x,t)\,dx,

the mean position weighted by the probability density ∣ψ∣2|\psi|^2. Momentum is represented instead by a derivative operator, p^=−iℏ ddx\hat p = -i\hbar\,\dfrac{d}{dx}, so that

⟨p⟩=∫−∞∞ψ∗(x,t)(−iℏ∂∂x)ψ(x,t) dx.\langle p\rangle = \int_{-\infty}^{\infty}\psi^*(x,t)\left(-i\hbar\frac{\partial}{\partial x}\right)\psi(x,t)\,dx.

where: ℏ=h/2π\hbar = h/2\pi is the reduced Planck constant, ψ∗\psi^* is the complex conjugate of ψ\psi, and the integral is taken over all xx for which ψ\psi is defined. This is the general pattern behind every expectation value in quantum mechanics: pick the operator for the quantity you want, sandwich it, integrate.

Good behaviour of ψ\psi. Not every function is an acceptable wave function. Physically reasonable ψ\psi must be single-valued, and ψ\psi together with dψ/dxd\psi/dx must stay continuous wherever the potential V(x)V(x) is finite (the derivative can jump only across an infinite discontinuity in VV). ψ\psi must also vanish, or stay bounded, as x→±∞x\to\pm\infty, otherwise the normalization integral could not converge to a finite number. These conditions are not decoration — the next chapter uses exactly this continuity to match solutions of the Schrödinger equation across a potential step or well.

Figure
A bell shaped probability density curve over position with the region under the curve between two marked points shaded to represent the probability of finding the particle there
Probability density |psi(x)|^2 and probability as area position x |psi(x)|^2 x1 x2 shaded area = P(x1 to x2) Total area under the whole curve equals 1 — the normalization condition.
Derivation ANormalizing a wave function to fix an unknown constant

The goal of this derivation is to show, with nothing beyond a single definite integral, how the normalization condition ∫−∞∞∣ψ∣2 dx=1\int_{-\infty}^{\infty}|\psi|^2\,dx = 1 is turned into an equation that pins down an otherwise-unknown multiplicative constant in front of a wave function, which is the single most common calculation this topic asks for on a problem set.

Step 1. Start from a wave function known up to an overall real constant AA, written ψ(x)=A f(x)\psi(x) = A\,f(x) on some region, and zero outside it, where f(x)f(x) is a given shape (for example a sine, an exponential, or a Gaussian). Since AA is real, ∣ψ(x)∣2=A2∣f(x)∣2|\psi(x)|^2 = A^2 |f(x)|^2, and the constant factors straight out of the normalization integral.

Step 2. Write the normalization condition explicitly over the region where ψ\psi is nonzero, say a≤x≤ba\le x\le b: A2∫ab∣f(x)∣2 dx=1A^2\displaystyle\int_a^b |f(x)|^2\,dx = 1. Everything about the shape of the wave function is now buried inside the single number I=∫ab∣f(x)∣2 dxI = \int_a^b|f(x)|^2\,dx, which you evaluate using ordinary calculus (trigonometric identities, exponential integrals, and so on, depending on ff).

Step 3. Solve for the constant: A2=1/IA^2 = 1/I, so A=1/IA = 1/\sqrt{I} (taking the positive root by convention, since an overall phase or sign of AA has no physical effect on ∣ψ∣2|\psi|^2). This is the normalization constant.

Step 4. Sanity-check the result: AA must have units of (length)−1/2^{-1/2} so that A2A^2 times a length (from the integral) is dimensionless, matching the requirement that ∫∣ψ∣2dx=1\int|\psi|^2dx=1 is a pure number. If your computed AA carries the wrong units, an arithmetic slip has occurred upstream.

Worked example

A particle confined to 0≤x≤L0\le x\le L has, at some instant, the real wave function ψ(x)=Asin⁡(πx/L)\psi(x) = A\sin(\pi x/L) on that interval and ψ(x)=0\psi(x)=0 outside it. Find AA, and then find the probability that the particle is found in the left half of the box, 0≤x≤L/20\le x\le L/2.

Step 1. Normalize. We need ∫0L∣ψ(x)∣2dx=1\displaystyle\int_0^L|\psi(x)|^2dx=1, i.e. A2∫0Lsin⁡2(πx/L) dx=1A^2\displaystyle\int_0^L\sin^2(\pi x/L)\,dx=1. Using sin⁡2θ=12(1−cos⁡2θ)\sin^2\theta=\tfrac12(1-\cos2\theta), the integral evaluates to L/2L/2. So A2(L/2)=1A^2(L/2)=1, giving A=2/LA=\sqrt{2/L}.

Step 2. Compute the left-half probability. P(0≤x≤L/2)=∫0L/22Lsin⁡2(πx/L) dx=2L[x2−L4πsin⁡ ⁣(2πxL)]0L/2=2L(L4−0)=12P(0\le x\le L/2) = \displaystyle\int_0^{L/2}\frac{2}{L}\sin^2(\pi x/L)\,dx = \dfrac{2}{L}\left[\dfrac{x}{2}-\dfrac{L}{4\pi}\sin\!\left(\dfrac{2\pi x}{L}\right)\right]_0^{L/2} = \dfrac{2}{L}\left(\dfrac{L}{4}-0\right) = \dfrac12.

Step 3. Extend to the expectation value. With this same normalized ψ\psi, ⟨x⟩=∫0Lψ∗xψ dx=2L∫0Lxsin⁡2(πx/L) dx\langle x\rangle = \displaystyle\int_0^L\psi^*x\psi\,dx = \dfrac{2}{L}\int_0^Lx\sin^2(\pi x/L)\,dx. Splitting via the same identity, the first term gives 12⋅L22=L24\tfrac12\cdot\tfrac{L^2}{2}=\tfrac{L^2}{4} and the cosine term integrates to zero over the interval (a standard result for ∫0Lxcos⁡(2πx/L)dx\int_0^Lx\cos(2\pi x/L)dx), so ⟨x⟩=2L⋅L24=L2\langle x\rangle = \dfrac{2}{L}\cdot\dfrac{L^2}{4} = \dfrac{L}{2}.

Step 4. Interpret. Both results are exactly what symmetry demands: ∣ψ∣2∝sin⁡2(πx/L)|\psi|^2\propto\sin^2(\pi x/L) is symmetric about the midpoint x=L/2x=L/2, so half the probability sits in each half of the box, and the mean position is precisely the centre. Final answers: A=2/LA=\sqrt{2/L}, P(left half)=1/2P(\text{left half})=1/2, and ⟨x⟩=L/2\langle x\rangle = L/2.

Common traps

  • Treating ψ\psi itself as a probability — it is ∣ψ∣2|\psi|^2, a real non-negative quantity, that plays that role; ψ\psi can be complex or negative.
  • Skipping normalization before computing an expectation value: if ∫∣ψ∣2dx≠1\int|\psi|^2dx\ne 1, the raw formula ⟨x⟩=∫ψ∗xψ dx\langle x\rangle=\int\psi^*x\psi\,dx gives the wrong number.
  • Confusing "probability" with "probability density" — ∣ψ(x)∣2|\psi(x)|^2 has units of probability per unit length; you must integrate over an interval to get a number between 0 and 1.
  • Believing ⟨x⟩\langle x\rangle tells you where the particle actually is on any single measurement — it is a statistical mean over many identically prepared particles, not an outcome.

Ask the tutor: want to work through normalizing a Gaussian wave packet, or see how the continuity of ψ\psi and dψ/dxd\psi/dx is actually applied at a potential step?

Why can psi be complex if only real numbers come out of experiments?

Because psi itself is never the thing you measure — only ∣ψ∣2=ψ∗ψ|\psi|^2 = \psi^*\psi is, and that combination is automatically real for any complex ψ\psi, exactly the way multiplying a complex number by its own conjugate always yields a real, non-negative result. Think of psi as bookkeeping machinery: it carries phase information essential for interference effects later in the course, but every prediction that reaches the lab comes out through ∣ψ∣2|\psi|^2 or through an operator sandwich, both of which are guaranteed real.

Do I always integrate from minus infinity to plus infinity to normalize?

Only if the wave function is genuinely nonzero everywhere. In practice psi is very often confined to some finite region, like the particle-in-a-box example, and zero outside it, so the normalization integral only needs to run over the region where psi is nonzero — extending it further just adds zero. Always check where psi is actually defined before setting up the integral limits.

Is there a quick way to know whether a wave function is acceptable without doing the full integral?

A fast first check is convergence: does ∣ψ∣2|\psi|^2 die off (or stay bounded in a finite box) as x goes to plus or minus infinity, so that the area under it can possibly be finite? If it blows up, no rescaling can ever normalize it, and it is not physically acceptable. Beyond that, check that psi and its derivative are continuous wherever the potential is finite — a psi with a kink at a point where V is finite is a red flag before you even attempt the integral.

Before this

Where this leads