Chapter 2: The Wave Function and the Schrödinger Equation — The Rules of the Game

"Where did we get that [equation] from? Nowhere. It is not possible to derive it from anything you know. It came out of the mind of Schrödinger." — Richard Feynman, The Feynman Lectures on Physics

In Chapter 1, we watched classical physics break down. The ultraviolet catastrophe, the photoelectric effect, atomic stability, wave-particle duality — these are not minor glitches in an otherwise perfect theory. They are evidence that nature operates by fundamentally different rules at the atomic scale. We ended that chapter with a question: If classical physics is wrong, what replaces it?

This chapter answers that question. We are going to build — from physical reasoning and mathematical necessity — the complete framework that governs quantum mechanics at its most basic level. By the end, you will have three things:

1. A mathematical object — the wave function $\psi(x,t)$ — that encodes everything knowable about a quantum system.
2. An equation of motion — the Schrödinger equation — that tells us how $\psi$ evolves in time.
3. A measurement rule — the Born rule — that connects the mathematics to what we actually observe in the laboratory.

These three elements constitute the "rules of the game." Everything else in this textbook — from the hydrogen atom to quantum computing — is the application of these rules to increasingly interesting and complex situations.

Let us begin.

2.1 From Experiments to a Theory: What Do We Need?

Before writing any equations, let us think carefully about what the experiments from Chapter 1 demand of any replacement theory. This is not merely a historical exercise — it is the physical reasoning that motivates the mathematical structure we are about to build.

The Experimental Constraints

Recall the key results:

Quantization. Energy comes in discrete packets. Planck's constant $\hbar$ (or equivalently $h = 2\pi\hbar$) sets the scale. The hydrogen atom has energy levels $E_n = -13.6\,\text{eV}/n^2$, not a continuum.

Wave-particle duality. Electrons, photons, and all quantum objects exhibit both wave-like behavior (interference, diffraction) and particle-like behavior (localized detection events). De Broglie's relation $\lambda = h/p$ connects the two descriptions.

Probabilistic outcomes. When we send identically prepared electrons through a double slit, each electron lands at a definite spot — but where it lands is not predictable. Only the statistical pattern is reproducible. The distribution of many events matches an interference pattern.

Non-classical statistics. The double-slit pattern is not the sum of the two single-slit patterns. Something interferes, but it is not a classical wave (because we detect individual particles).

What the Theory Must Provide

These constraints force our hand:

1. We need a wave-like mathematical object. The interference patterns demand something that can exhibit constructive and destructive interference — something with phase and amplitude. Call it $\psi$.

2. We need $\psi$ to give probabilities, not certainties. The individual outcomes are unpredictable; only the statistical distribution is deterministic.

3. We need an equation governing how $\psi$ changes in time. Just as Newton's second law governs how a classical particle evolves, we need a dynamical equation for $\psi$.

4. The equation must be linear. Interference requires the superposition principle: if $\psi_1$ and $\psi_2$ are valid states, then $\psi_1 + \psi_2$ must also be a valid state. Only a linear equation guarantees this.

5. The equation must reduce to classical physics in the appropriate limit. Quantum mechanics is not a replacement for classical mechanics — it is a generalization. When objects are large and heavy (when $\hbar$ is negligible compared to the relevant action), we must recover Newton's laws.

💡 Key Insight. Notice how the experimental facts constrain the mathematical structure. This is not a case of physicists choosing an equation out of thin air. The mathematics is forced upon us by what nature does. Feynman's quote at the chapter opening is slightly misleading — the Schrödinger equation was not arbitrary. It is the simplest equation consistent with everything we know.

Spaced Review from Chapter 1

Before we proceed, make sure you can answer these quickly:

🔄 Check Your Understanding. 1. What is the de Broglie wavelength of an electron with kinetic energy 10 eV? (Answer: $\lambda = h/p = h/\sqrt{2mE} \approx 0.39\,\text{nm}$.) 2. Why does the photoelectric effect require $E = h\nu$? What fails if energy is spread continuously across the wavefront? 3. In the double-slit experiment with single electrons, what happens if you close one slit? Why is this significant?

2.2 The Wave Function $\psi(x,t)$

Definition

The wave function $\psi(x,t)$ is a complex-valued function that completely describes the quantum state of a particle. For a single particle moving in one dimension, $\psi$ assigns a complex number to every point $x$ in space and every time $t$:

$$\psi: \mathbb{R} \times \mathbb{R} \to \mathbb{C}, \quad (x,t) \mapsto \psi(x,t).$$

The wave function is not directly observable. You will never measure $\psi$ in a laboratory. What you can measure is the probability of finding the particle in a given region of space, and this is where Max Born's brilliant insight enters.

The Born Rule

🔵 Historical Note. In 1926, Max Born proposed a radical interpretation of Schrödinger's wave function: $|\psi(x,t)|^2$ gives the probability density for finding the particle at position $x$ at time $t$. This was controversial — Schrödinger himself preferred to interpret $|\psi|^2$ as a literal charge density. Born's statistical interpretation won out, and he received the Nobel Prize for it in 1954, nearly thirty years later.

The Born rule states:

$$\boxed{P(a \leq x \leq b, t) = \int_a^b |\psi(x,t)|^2\, dx}$$

The quantity $|\psi(x,t)|^2\,dx$ is the probability of finding the particle between $x$ and $x + dx$ at time $t$. We call $|\psi(x,t)|^2$ the probability density.

Several points require emphasis:

$\psi$ itself is not a probability. The wave function is complex-valued; probabilities are real and non-negative. You must take the modulus squared $|\psi|^2 = \psi^*\psi$ to extract physical predictions. The wave function is sometimes called a probability amplitude to distinguish it from probability itself.

$|\psi|^2$ is a density, not a probability. You cannot ask "what is the probability of finding the particle at exactly $x = 3$?" — the answer to that is zero for any continuous distribution. You must ask about a range of positions.

$\psi$ has units. In one dimension, $|\psi|^2$ is a probability per unit length, so $[\psi] = \text{m}^{-1/2}$. In three dimensions, $[\psi] = \text{m}^{-3/2}$.

Normalization

Since the particle must be found somewhere, the total probability must equal one:

$$\boxed{\int_{-\infty}^{\infty} |\psi(x,t)|^2\, dx = 1}$$

This is the normalization condition. A wave function satisfying this condition is said to be normalized. If a wave function is not normalized, we can (usually) fix this by multiplying by an appropriate constant:

$$\psi_{\text{normalized}} = A\,\psi, \quad A = \left(\int_{-\infty}^{\infty} |\psi(x,t)|^2\, dx\right)^{-1/2}.$$

⚠️ Common Misconception. Not every mathematical function can be normalized. The function $\psi(x) = e^{ikx}$ (a pure plane wave) has $|\psi|^2 = 1$ everywhere, so $\int_{-\infty}^{\infty} |\psi|^2\,dx = \infty$. This function is not normalizable in the usual sense. We will deal with such functions using wave packets and the Dirac delta normalization later (Ch 6, Ch 9). For now, we require $\psi$ to be square-integrable: $\int |\psi|^2\,dx < \infty$.

What $\psi$ Is NOT

Let us be very precise about what the wave function is not:

1. $\psi$ is not a physical wave like a water wave or a sound wave. Those are real-valued disturbances in a physical medium. $\psi$ is complex-valued, has no medium, and exists in configuration space (which for $N$ particles is $3N$-dimensional, not 3-dimensional).

2. $\psi$ does not describe the trajectory of a particle. In classical mechanics, we track $x(t)$ — the position as a function of time. In quantum mechanics, there is no trajectory. There is only $\psi(x,t)$, from which we can calculate the probability of various outcomes.

3. $|\psi|^2$ is not a measure of our ignorance. This is a subtle and important point. When we say "the particle has a 30% probability of being found in region $A$," we do not mean that the particle is in some definite place and we just do not know where. The standard interpretation (which we will examine critically in Ch 24 and Ch 28) says that the particle does not have a definite position prior to measurement.

Worked Example 2.1: Normalizing a Wave Function

Consider the wave function $\psi(x) = Ae^{-\alpha x^2}$ where $\alpha > 0$ and $A$ is a normalization constant.

Task: Find $A$ such that $\psi$ is normalized.

Solution: We require

$$\int_{-\infty}^{\infty} |A|^2 e^{-2\alpha x^2}\,dx = 1.$$

The Gaussian integral gives us $\int_{-\infty}^{\infty} e^{-2\alpha x^2}\,dx = \sqrt{\pi/(2\alpha)}$. Therefore:

$$|A|^2 \sqrt{\frac{\pi}{2\alpha}} = 1 \implies A = \left(\frac{2\alpha}{\pi}\right)^{1/4}.$$

The normalized wave function is $\psi(x) = \left(\frac{2\alpha}{\pi}\right)^{1/4} e^{-\alpha x^2}$.

This Gaussian wave function will appear throughout quantum mechanics. It is the ground state of the quantum harmonic oscillator (Ch 4) and plays a central role in coherent states (Ch 27).

🔄 Check Your Understanding. 1. What are the units of $\alpha$ in the example above? 2. Verify that $|\psi(x)|^2$ has the shape of a bell curve centered at $x = 0$. What is its width (standard deviation)? 3. If a particle's wave function is $\psi(x) = B\,x\,e^{-\alpha x^2}$, find $B$. (Hint: the integral $\int_{-\infty}^{\infty} x^2 e^{-2\alpha x^2}\,dx = \frac{1}{2}\sqrt{\pi/(2\alpha)^3}$.)

2.3 The Schrödinger Equation: Time-Dependent Form

We have our mathematical object ($\psi$) and our measurement rule (Born rule). Now we need the equation of motion. How does $\psi$ evolve in time?

Motivation from de Broglie Waves

Consider a free particle (no forces) with definite momentum $p$ and energy $E$. From Chapter 1, we know:

• De Broglie: $p = \hbar k$ where $k = 2\pi/\lambda$ is the wave number.
• Planck-Einstein: $E = \hbar\omega$ where $\omega$ is the angular frequency.

The simplest wave with these properties is:

$$\psi(x,t) = Ae^{i(kx - \omega t)}.$$

This is a plane wave propagating in the $+x$ direction. Let us see what mathematical operations extract the physical quantities:

Energy. Take the time derivative:

$$i\hbar \frac{\partial \psi}{\partial t} = i\hbar \cdot (-i\omega) \psi = \hbar\omega\,\psi = E\,\psi.$$

Momentum. Take the spatial derivative:

$$-i\hbar \frac{\partial \psi}{\partial x} = -i\hbar \cdot (ik)\psi = \hbar k\,\psi = p\,\psi.$$

Kinetic energy. Take the second spatial derivative:

$$-\frac{\hbar^2}{2m}\frac{\partial^2 \psi}{\partial x^2} = -\frac{\hbar^2}{2m}(ik)^2\psi = \frac{\hbar^2 k^2}{2m}\psi = \frac{p^2}{2m}\psi = T\,\psi,$$

where $T = p^2/(2m)$ is the kinetic energy.

For a free particle, the total energy is just the kinetic energy: $E = p^2/(2m)$. Combining the expressions above:

$$i\hbar \frac{\partial\psi}{\partial t} = -\frac{\hbar^2}{2m}\frac{\partial^2\psi}{\partial x^2}.$$

This is the Schrödinger equation for a free particle.

Adding a Potential

Now suppose the particle moves in a potential $V(x)$. The total energy is $E = T + V = p^2/(2m) + V(x)$. We simply add the potential energy:

$$\boxed{i\hbar \frac{\partial \psi(x,t)}{\partial t} = \left[-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2} + V(x)\right]\psi(x,t)}$$

This is the time-dependent Schrödinger equation (TDSE) in one dimension. It is the fundamental dynamical equation of (nonrelativistic) quantum mechanics.

The Hamiltonian Operator

The quantity in square brackets has a name. Define the Hamiltonian operator:

$$\hat{H} = -\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2} + V(x) = \frac{\hat{p}^2}{2m} + V(\hat{x}),$$

where $\hat{p} = -i\hbar\,\partial/\partial x$ is the momentum operator (we will say more about this in Section 2.5). Then the TDSE takes the compact form:

$$\boxed{i\hbar \frac{\partial \psi}{\partial t} = \hat{H}\psi}$$

This is beautiful in its simplicity: the rate of change of the quantum state is determined by the Hamiltonian, just as the rate of change of a classical state is determined by the Hamiltonian in classical mechanics.

🔵 Historical Note. Erwin Schrödinger published this equation in 1926, in a series of four papers that constitute one of the most remarkable achievements in the history of physics. He was 38 years old and working at the University of Zurich. Legend has it that the key insight came during a Christmas holiday in the Swiss Alps — a story we explore in Case Study 2 of this chapter.

Key Properties of the TDSE

1. First order in time. Unlike Newton's second law ($\ddot{x} = F/m$, second order in time), the TDSE is first order in $t$. This means that given $\psi(x,0)$ — the wave function at one instant — the future evolution is completely determined. There is no need to specify $\partial\psi/\partial t$ separately (as you would need $\dot{x}(0)$ for Newton's law).

2. Linear. If $\psi_1$ and $\psi_2$ are both solutions, then $c_1\psi_1 + c_2\psi_2$ is also a solution (for any complex constants $c_1, c_2$). This is the superposition principle — the single most important property of quantum mechanics.

3. Deterministic. Given $\psi(x,0)$, the Schrödinger equation determines $\psi(x,t)$ for all future $t$ with absolute certainty. The randomness in quantum mechanics enters only through the measurement process (the Born rule), not through the evolution of the wave function.

4. Preserves normalization. If $\psi$ is normalized at $t = 0$, it remains normalized for all $t > 0$. We will prove this in Section 2.6.

💡 Key Insight. Point 3 deserves emphasis. The Schrödinger equation is completely deterministic. There is no randomness in how $\psi$ evolves — the randomness is in what happens when you measure. This is a persistent source of confusion. The quantum state evolves smoothly and predictably until a measurement collapses it. The tension between these two types of evolution (smooth Schrödinger evolution vs. abrupt measurement collapse) is the measurement problem, one of the deepest open questions in physics. We will confront this directly in Chapter 28.

🔄 Check Your Understanding. 1. Why must the Schrödinger equation be first order in time but second order in space? 2. If $\psi_1(x,t)$ and $\psi_2(x,t)$ are both solutions of the TDSE for the same potential $V(x)$, show that $3\psi_1 - 2i\psi_2$ is also a solution. 3. The classical wave equation $\partial^2 u/\partial t^2 = v^2 \partial^2 u/\partial x^2$ is second order in time and real. What physical consequences follow from the TDSE being first order in time and involving complex numbers?

2.4 The Time-Independent Schrödinger Equation

The time-dependent Schrödinger equation is the fundamental law, but solving it directly is often difficult. A powerful technique — separation of variables — allows us to reduce it to a simpler equation in many important cases.

Separation of Variables

Assume the potential $V$ does not depend on time (this covers an enormous range of physical problems). We look for solutions of the form:

$$\psi(x,t) = \phi(x)\,f(t),$$

where $\phi$ depends only on $x$ and $f$ depends only on $t$. Substituting into the TDSE:

$$i\hbar\,\phi(x)\frac{df}{dt} = f(t)\left[-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} + V(x)\phi(x)\right].$$

Dividing both sides by $\phi(x)f(t)$:

$$i\hbar\frac{1}{f}\frac{df}{dt} = \frac{1}{\phi}\left[-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} + V(x)\phi\right].$$

The left side depends only on $t$; the right side depends only on $x$. Since $x$ and $t$ are independent variables, both sides must equal a constant. Call it $E$:

$$i\hbar\frac{df}{dt} = Ef, \qquad \hat{H}\phi(x) = E\phi(x).$$

The Time Part

The first equation has the immediate solution:

$$f(t) = e^{-iEt/\hbar}.$$

This is a complex exponential — a pure oscillation in the complex plane with angular frequency $\omega = E/\hbar$, confirming the Planck-Einstein relation $E = \hbar\omega$.

The Spatial Part: The Time-Independent Schrödinger Equation

The second equation is the time-independent Schrödinger equation (TISE):

$$\boxed{\hat{H}\phi(x) = E\phi(x) \qquad \Longleftrightarrow \qquad -\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} + V(x)\phi(x) = E\phi(x)}$$

This is an eigenvalue equation: the operator $\hat{H}$ acting on $\phi$ gives back $\phi$ multiplied by a constant $E$. The function $\phi$ is an eigenfunction of $\hat{H}$, and $E$ is the corresponding eigenvalue (the energy).

Stationary States

The full separable solution is:

$$\psi(x,t) = \phi(x)\,e^{-iEt/\hbar}.$$

These are called stationary states. The name comes from a remarkable property: the probability density is time-independent:

$$|\psi(x,t)|^2 = |\phi(x)|^2 |e^{-iEt/\hbar}|^2 = |\phi(x)|^2.$$

The probability of finding the particle at any location does not change with time. The particle is not "sitting still" — but its probability distribution is static. Stationary states are the quantum analogue of classical equilibrium configurations.

Why are stationary states important? Two reasons:

1. They are the building blocks of general solutions. Because the Schrödinger equation is linear, the general solution is a superposition of stationary states:

$$\psi(x,t) = \sum_n c_n \phi_n(x) e^{-iE_n t/\hbar},$$

where the sum runs over all energy eigenstates. (For continuous spectra, the sum becomes an integral.) Determining the coefficients $c_n$ from the initial condition $\psi(x,0)$ is a Fourier-type problem that we will tackle in Chapter 3.

2. They give the allowed energies. Not every value of $E$ produces a normalizable $\phi$. The requirement that $\phi$ be physically acceptable (normalizable, continuous, etc. — see Section 2.7) restricts $E$ to a discrete set of values for bound states. This is the origin of energy quantization — the very first quantum mystery from Chapter 1, now derived rather than postulated.

Worked Example 2.2: Verifying a Stationary State

Suppose a particle in a potential $V(x)$ has the wave function $\psi(x,t) = \phi(x)e^{-iEt/\hbar}$ where $\phi(x)$ satisfies the TISE. Show that $\langle E\rangle = E$ (the expectation value of energy equals the eigenvalue) and that the variance in energy is zero.

Solution:

The expectation value of $\hat{H}$ is:

$$\langle E \rangle = \int_{-\infty}^{\infty} \psi^* \hat{H}\psi\,dx = \int_{-\infty}^{\infty} \phi^* e^{iEt/\hbar}\cdot E\phi\, e^{-iEt/\hbar}\,dx = E\int_{-\infty}^{\infty} |\phi|^2\,dx = E.$$

For the variance, we need $\langle E^2\rangle$:

$$\langle E^2 \rangle = \int_{-\infty}^{\infty} \psi^* \hat{H}^2\psi\,dx = E^2 \int_{-\infty}^{\infty} |\phi|^2\,dx = E^2.$$

Therefore $\sigma_E^2 = \langle E^2\rangle - \langle E\rangle^2 = E^2 - E^2 = 0$.

A stationary state has a definite energy — no uncertainty whatsoever. This is one of the hallmarks of an eigenstate.

🔗 Connection. The fact that eigenstates have zero variance in their eigenvalue will be generalized in Chapter 6 to all observables: if a system is in an eigenstate of operator $\hat{A}$, then a measurement of $A$ yields the eigenvalue with certainty. This leads directly to the uncertainty principle (Ch 6) — if you cannot be in a simultaneous eigenstate of both $\hat{x}$ and $\hat{p}$, then you cannot have definite values of both position and momentum.

2.5 Operators and Observables: A First Look

We have already seen the Hamiltonian operator $\hat{H}$ and the momentum operator $\hat{p}$. Let us now discuss these more systematically. This is a preview of the full operator formalism developed in Chapter 6 — here we introduce only what we need to compute measurable quantities.

The Central Idea

In quantum mechanics, every measurable quantity (called an observable) is associated with a mathematical operator. An operator is a rule that takes a function and produces another function. The key operators for a single particle in one dimension are:

The prescription is: replace $p$ by $\hat{p} = -i\hbar\,\partial/\partial x$ and $x$ by $\hat{x} = x$. This is called the correspondence principle for operators.

Where Does $\hat{p} = -i\hbar\,\partial/\partial x$ Come From?

We already saw the answer in Section 2.3. For a plane wave $\psi = Ae^{i(kx - \omega t)}$:

$$-i\hbar\frac{\partial\psi}{\partial x} = \hbar k\,\psi = p\,\psi.$$

The operator $-i\hbar\,\partial/\partial x$ extracts the momentum from a state of definite momentum. We define this to be the momentum operator for general wave functions.

Why the factor of $-i$? Because we need a Hermitian operator (one whose eigenvalues are real), since momentum is a real-valued observable. We will prove Hermiticity in Chapter 6.

Expectation Values

Given a normalized wave function $\psi(x,t)$, the expectation value of an observable $A$ is:

$$\boxed{\langle A \rangle = \int_{-\infty}^{\infty} \psi^*(x,t)\,\hat{A}\,\psi(x,t)\,dx}$$

This is the average value you would obtain if you prepared many identical systems in the state $\psi$ and measured $A$ on each one.

Position:

$$\langle x \rangle = \int_{-\infty}^{\infty} \psi^*\,x\,\psi\,dx = \int_{-\infty}^{\infty} x\,|\psi|^2\,dx.$$

This is simply the weighted average of $x$ with the probability density $|\psi|^2$ — exactly what you would expect.

Momentum:

$$\langle p \rangle = \int_{-\infty}^{\infty} \psi^*\left(-i\hbar\frac{\partial}{\partial x}\right)\psi\,dx.$$

This is not simply the weighted average of some function — the operator acts on $\psi$ by differentiation. The order matters: $\hat{p}$ acts to the right on $\psi$, and then we multiply by $\psi^*$ from the left and integrate.

Energy:

$$\langle E \rangle = \langle H \rangle = \int_{-\infty}^{\infty} \psi^*\,\hat{H}\,\psi\,dx.$$

Worked Example 2.3: Expectation Values for a Gaussian

For the normalized Gaussian $\psi(x) = \left(\frac{2\alpha}{\pi}\right)^{1/4}e^{-\alpha x^2}$, compute $\langle x \rangle$, $\langle x^2 \rangle$, $\langle p \rangle$, and $\langle p^2 \rangle$.

$\langle x \rangle$:

$$\langle x\rangle = \sqrt{\frac{2\alpha}{\pi}}\int_{-\infty}^{\infty} x\,e^{-2\alpha x^2}\,dx = 0,$$

because the integrand is an odd function of $x$ (and the Gaussian is symmetric about $x=0$).

$\langle x^2\rangle$:

$$\langle x^2\rangle = \sqrt{\frac{2\alpha}{\pi}}\int_{-\infty}^{\infty} x^2\,e^{-2\alpha x^2}\,dx = \sqrt{\frac{2\alpha}{\pi}}\cdot\frac{1}{2}\sqrt{\frac{\pi}{(2\alpha)^3}} = \frac{1}{4\alpha}.$$

$\langle p\rangle$:

$$\langle p\rangle = \int_{-\infty}^{\infty}\psi^*\left(-i\hbar\frac{d}{dx}\right)\psi\,dx = -i\hbar\sqrt{\frac{2\alpha}{\pi}}\int_{-\infty}^{\infty}e^{-\alpha x^2}(-2\alpha x)e^{-\alpha x^2}\,dx = 0,$$

again by odd symmetry.

$\langle p^2\rangle$:

$$\frac{d^2\psi}{dx^2} = \left(\frac{2\alpha}{\pi}\right)^{1/4}(4\alpha^2 x^2 - 2\alpha)e^{-\alpha x^2}.$$

$$\langle p^2\rangle = -\hbar^2\int_{-\infty}^{\infty}\psi^*\frac{d^2\psi}{dx^2}\,dx = \hbar^2\alpha.$$

Note the uncertainty product:

$$\sigma_x\sigma_p = \sqrt{\langle x^2\rangle - \langle x\rangle^2}\cdot\sqrt{\langle p^2\rangle - \langle p\rangle^2} = \sqrt{\frac{1}{4\alpha}}\cdot\sqrt{\hbar^2\alpha} = \frac{\hbar}{2}.$$

This is exactly $\hbar/2$ — the minimum allowed by the Heisenberg uncertainty principle (which we will derive in Chapter 6). The Gaussian wave function is a minimum-uncertainty state, saturating the uncertainty bound. This is not a coincidence; it is a deep property of Gaussians.

🔗 Connection. The uncertainty principle $\sigma_x \sigma_p \geq \hbar/2$ is one of the most famous results in quantum mechanics. Here we have verified it for a specific case. The general proof requires the Cauchy-Schwarz inequality and the commutator $[\hat{x},\hat{p}] = i\hbar$, which we will develop in Chapter 6.

🔄 Check Your Understanding. 1. For the Gaussian wave function above, compute $\langle T \rangle = \langle p^2\rangle/(2m)$. What is the kinetic energy of a particle in this state? 2. If the particle is in a potential $V(x) = \frac{1}{2}m\omega^2 x^2$ (harmonic oscillator), compute $\langle V \rangle$ and $\langle E \rangle = \langle T \rangle + \langle V \rangle$. What value of $\alpha$ minimizes $\langle E \rangle$? 3. Can $\langle p \rangle$ be nonzero for a real-valued wave function? Prove your answer.

2.6 Probability Current and Conservation

One of the key properties of the Schrödinger equation is that it conserves probability. If a wave function is normalized at $t = 0$, it stays normalized forever. This is not a trivial statement — it is a consequence of the specific mathematical structure of the TDSE. Let us prove it.

Proof of Probability Conservation

We want to show that $\frac{d}{dt}\int_{-\infty}^{\infty}|\psi|^2\,dx = 0$.

Compute:

$$\frac{\partial}{\partial t}|\psi|^2 = \frac{\partial}{\partial t}(\psi^*\psi) = \psi^*\frac{\partial\psi}{\partial t} + \frac{\partial\psi^*}{\partial t}\psi.$$

From the TDSE: $\frac{\partial\psi}{\partial t} = \frac{1}{i\hbar}\hat{H}\psi = \frac{1}{i\hbar}\left(-\frac{\hbar^2}{2m}\frac{\partial^2\psi}{\partial x^2} + V\psi\right)$.

Taking the complex conjugate (and noting that $V$ is real):

$$\frac{\partial\psi^*}{\partial t} = -\frac{1}{i\hbar}\left(-\frac{\hbar^2}{2m}\frac{\partial^2\psi^*}{\partial x^2} + V\psi^*\right).$$

Substituting:

$$\frac{\partial |\psi|^2}{\partial t} = \frac{1}{i\hbar}\left[-\frac{\hbar^2}{2m}\left(\psi^*\frac{\partial^2\psi}{\partial x^2} - \frac{\partial^2\psi^*}{\partial x^2}\psi\right)\right] = -\frac{\partial}{\partial x}\left[\frac{\hbar}{2mi}\left(\psi^*\frac{\partial\psi}{\partial x} - \frac{\partial\psi^*}{\partial x}\psi\right)\right].$$

(The potential energy terms cancel — this is crucial and is why $V$ must be real.)

The Continuity Equation

Define the probability current density:

$$\boxed{j(x,t) = \frac{\hbar}{2mi}\left(\psi^*\frac{\partial\psi}{\partial x} - \frac{\partial\psi^*}{\partial x}\psi\right) = \frac{\hbar}{m}\,\text{Im}\left(\psi^*\frac{\partial\psi}{\partial x}\right)}$$

Then we have the continuity equation:

$$\boxed{\frac{\partial \rho}{\partial t} + \frac{\partial j}{\partial x} = 0}$$

where $\rho = |\psi|^2$ is the probability density and $j$ is the probability current density.

This has exactly the same form as the continuity equation in fluid dynamics ($\partial\rho/\partial t + \nabla\cdot\mathbf{j} = 0$) or charge conservation in electromagnetism. It says: probability is neither created nor destroyed — it can only flow from one region to another.

📊 By the Numbers. The continuity equation guarantees that $\frac{d}{dt}\int_{-\infty}^{\infty}|\psi|^2\,dx = -[j]_{-\infty}^{+\infty} = 0$ for normalizable wave functions (since $\psi \to 0$ as $x \to \pm\infty$). Probability is globally conserved.

Physical Interpretation

The probability current $j(x,t)$ tells you the rate at which probability flows past the point $x$. If $j > 0$, probability flows to the right; if $j < 0$, to the left.

Example: For a plane wave $\psi = Ae^{i(kx-\omega t)}$:

$$j = \frac{\hbar k}{m}|A|^2 = \frac{p}{m}|A|^2 = v|A|^2.$$

The probability current is proportional to the velocity — exactly what you would expect classically for a beam of particles. This is reassuring: the quantum formalism reduces to sensible results in cases where we have classical intuition.

🧪 Experiment. The probability current plays a central role in tunneling (Ch 3). When a particle encounters a potential barrier, the probability current on the far side of the barrier is smaller than on the near side — the "missing" current is reflected. The ratio of transmitted to incident probability current is the tunneling probability. This is how tunnel diodes and scanning tunneling microscopes work.

2.7 Normalization, Boundary Conditions, and Physical Solutions

Not every mathematical solution to the Schrödinger equation is a physical wave function. We need to impose conditions that ensure $\psi$ represents a realizable quantum state.

Requirements for Physical Wave Functions

A physically acceptable wave function must satisfy:

1. Square-integrability (normalizability). $\int_{-\infty}^{\infty}|\psi|^2\,dx < \infty$. The particle must be findable somewhere with total probability 1. Functions that blow up at infinity (like $e^{+\alpha x^2}$) or do not decay (like $e^{ikx}$, with an important caveat) are not normalizable in the usual sense.

2. Continuity. $\psi(x)$ must be continuous everywhere. A discontinuous wave function would imply an infinite probability current at the discontinuity (since $j$ involves $\partial\psi/\partial x$).

3. Continuous first derivative. $d\psi/dx$ must be continuous wherever $V(x)$ is finite. If $V$ has an infinite discontinuity (like an infinite potential wall), then $d\psi/dx$ may be discontinuous there, but $\psi$ itself must still be continuous.

4. Single-valuedness. $\psi(x)$ must have a unique value at each point. (This becomes important in three dimensions with angular coordinates, where the requirement $\psi(\phi + 2\pi) = \psi(\phi)$ leads to quantization of angular momentum — see Chapter 5.)

Why These Conditions Lead to Quantization

Here is the conceptual payoff. For a particle in a bound state (confined to a finite region by a potential), the boundary conditions at $x \to \pm\infty$ require $\psi \to 0$. Combined with the TISE, these boundary conditions can only be satisfied for specific values of $E$.

Consider, schematically, the TISE:

$$-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} = [E - V(x)]\phi.$$

In a region where $E > V(x)$ (the classically allowed region), $\phi$ oscillates. In a region where $E < V(x)$ (the classically forbidden region), $\phi$ is exponential. For $\phi$ to decay to zero at both $x \to +\infty$ and $x \to -\infty$, the oscillating and exponential solutions must match up correctly at the boundaries. This "matching" works only for special values of $E$.

This is exactly how quantization arises — not from a postulate, but from the mathematics of boundary conditions. Energy quantization is as natural and inevitable as the quantization of standing wave frequencies on a guitar string.

⚠️ Common Misconception. Students sometimes think quantization is a separate postulate of quantum mechanics. It is not. Quantization emerges automatically from the Schrödinger equation plus boundary conditions. The Schrödinger equation is the postulate; quantization is a consequence.

Worked Example 2.4: Why $E$ Is Quantized (Qualitative)

Consider a particle in a potential well — say, $V(x) = 0$ for $0 < x < L$ and $V = \infty$ outside (the infinite square well, which we solve in detail in Chapter 3).

Inside the well, the TISE is $-\frac{\hbar^2}{2m}\phi'' = E\phi$, with solutions $\phi = A\sin(kx) + B\cos(kx)$ where $k = \sqrt{2mE}/\hbar$.

The boundary conditions require $\phi(0) = 0$ and $\phi(L) = 0$ (since $\psi$ must vanish where the potential is infinite).

$\phi(0) = 0$ forces $B = 0$, so $\phi = A\sin(kx)$.

$\phi(L) = 0$ forces $\sin(kL) = 0$, so $kL = n\pi$ for $n = 1, 2, 3, \ldots$

Therefore:

$$E_n = \frac{\hbar^2 k_n^2}{2m} = \frac{n^2\pi^2\hbar^2}{2mL^2}, \quad n = 1, 2, 3, \ldots$$

The energy is quantized. Only specific values are allowed. The ground state energy $E_1 = \pi^2\hbar^2/(2mL^2)$ is not zero — there is a minimum energy. This zero-point energy is a purely quantum effect with no classical analogue.

🔗 Connection. We will solve the infinite square well completely in Chapter 3, including the full set of normalized eigenfunctions, orthogonality, and the expansion of arbitrary initial conditions. The quantum harmonic oscillator (Ch 4) provides a more physically realistic example where the same logic of boundary conditions leading to quantized energies plays out with Hermite polynomials instead of sines.

🔄 Check Your Understanding. 1. For the infinite square well, why is $n = 0$ not allowed? (What would the wave function look like?) 2. If you double the width $L$ of the well, what happens to the energy levels? Is this consistent with the classical limit? 3. Can the wave function itself be discontinuous at a point where $V$ has a finite discontinuity? What about at an infinite discontinuity?

2.8 Superposition: The First Truly Quantum Idea

We have saved the most important concept for the end of the main development. Everything we have discussed so far — wave functions, the Schrödinger equation, the Born rule — is surprising and counterintuitive. But superposition is the feature that makes quantum mechanics fundamentally, irreducibly different from any classical theory.

The Superposition Principle

Because the Schrödinger equation is linear, if $\psi_1(x,t)$ and $\psi_2(x,t)$ are both solutions, then

$$\psi(x,t) = c_1\psi_1(x,t) + c_2\psi_2(x,t)$$

is also a solution, for any complex constants $c_1$ and $c_2$. This is the superposition principle.

Applied to stationary states, the general solution is:

$$\psi(x,t) = \sum_n c_n \phi_n(x)\,e^{-iE_n t/\hbar}.$$

This looks innocent. After all, classical wave equations also have superposition. But in quantum mechanics, superposition has a radically different meaning because of the Born rule.

Why Quantum Superposition Is Fundamentally Different

🧩 Productive Struggle. Before reading further, consider this: if $\phi_1(x)$ is a state where the particle is localized near $x = a$, and $\phi_2(x)$ is a state where the particle is localized near $x = b$, what is the state $\psi = \frac{1}{\sqrt{2}}(\phi_1 + \phi_2)$? Is the particle at $a$, at $b$, or somewhere else? What does $|\psi|^2$ look like? Take a moment to think about this before reading on.

Consider a particle in the superposition $\psi = \frac{1}{\sqrt{2}}(\phi_1 + \phi_2)$ where $\phi_1$ and $\phi_2$ are states localized at positions $a$ and $b$ respectively. The probability density is:

$$|\psi|^2 = \frac{1}{2}|\phi_1|^2 + \frac{1}{2}|\phi_2|^2 + \frac{1}{2}(\phi_1^*\phi_2 + \phi_2^*\phi_1) = \frac{1}{2}|\phi_1|^2 + \frac{1}{2}|\phi_2|^2 + \text{Re}(\phi_1^*\phi_2).$$

The first two terms are what you would get from a classical mixture — half the time at $a$, half the time at $b$. But the third term is the interference term $\text{Re}(\phi_1^*\phi_2)$. This term can be positive or negative, leading to constructive or destructive interference. It has no classical analogue.

In classical probability theory, if a ball is in box $A$ with probability 1/2 and in box $B$ with probability 1/2, then the probability of finding it in any region is just $\frac{1}{2}P_A + \frac{1}{2}P_B$. There is no interference term. The ball is in one box or the other; we just do not know which.

In quantum mechanics, the interference term is real and measurable. The double-slit experiment (Ch 1) is precisely a demonstration of this interference term — and the fringe pattern cannot be explained by saying "the electron goes through one slit or the other."

🚪 THRESHOLD CONCEPT: Superposition Is Fundamentally New

The idea: A quantum superposition $\psi = c_1\psi_1 + c_2\psi_2$ is not a statement of ignorance. It does not mean "the system is in state $\psi_1$ or state $\psi_2$, and we don't know which." It means the system is in a genuinely new state $\psi$ that is neither $\psi_1$ nor $\psi_2$ but has properties of both — and properties that neither has alone (interference).

Why it matters: This single idea underlies everything in quantum mechanics: interference (Ch 1, 7), the uncertainty principle (Ch 6), entanglement (Ch 11, 24), quantum computing (Ch 25, 40), and the measurement problem (Ch 28). If you understand superposition correctly, the rest of quantum mechanics is a series of increasingly sophisticated applications of this principle. If you misunderstand it as classical ignorance, you will be confused by every subsequent chapter.

The test: You understand superposition when you can explain why the double-slit interference pattern disappears when you observe which slit the electron goes through — not because the observation physically disturbs the electron (a common but incorrect explanation), but because the observation eliminates the superposition by correlating the electron's path with a macroscopic detector.

Analogy (and its limits): Superposition is sometimes compared to a coin spinning in the air — it is "both heads and tails" until it lands. This analogy captures the indeterminacy but misses the interference. A better analogy: superposition is like a musical chord. A chord is not "either C or E or G" — it is a genuinely new sound that contains all three notes simultaneously and has properties (harmony, beating, timbre) that no single note has alone.

Superposition and Measurement

If the particle is in a superposition of energy eigenstates $\psi = \sum_n c_n \phi_n e^{-iE_n t/\hbar}$, and we measure the energy, what happens?

The Born rule tells us: we obtain the eigenvalue $E_n$ with probability $|c_n|^2$. After the measurement, the system is in the state $\phi_n$ — the superposition has collapsed to a single eigenstate.

The normalization condition requires $\sum_n |c_n|^2 = 1$ — the probabilities of all possible outcomes sum to one.

⚖️ Interpretation. The "collapse" of the wave function upon measurement is one of the most debated aspects of quantum mechanics. Does the wave function really collapse, or does the universe split into branches (many-worlds interpretation)? Is the wave function even real, or is it just a calculational tool (epistemic interpretation)? We flag these questions here and will explore them seriously in Chapter 24 (Bell's theorem and entanglement) and Chapter 28 (the measurement problem). For now, we use "collapse" as a practical computational rule: after measuring energy and obtaining $E_n$, compute subsequent evolution using $\psi = \phi_n$.

Worked Example 2.5: Superposition in the Infinite Well

A particle in an infinite square well of width $L$ is in the state:

$$\psi(x,0) = A\left[\phi_1(x) + 2\phi_2(x)\right],$$

where $\phi_n(x) = \sqrt{2/L}\sin(n\pi x/L)$ are the normalized energy eigenstates.

(a) Find $A$.

Since $\phi_1$ and $\phi_2$ are orthonormal ($\int \phi_m^*\phi_n\,dx = \delta_{mn}$):

$$1 = |A|^2(1 + 4) = 5|A|^2 \implies A = \frac{1}{\sqrt{5}}.$$

(b) What is $\psi(x,t)$?

$$\psi(x,t) = \frac{1}{\sqrt{5}}\left[\phi_1(x)e^{-iE_1 t/\hbar} + 2\phi_2(x)e^{-iE_2 t/\hbar}\right].$$

(c) If you measure the energy, what values might you get, and with what probabilities?

$E_1 = \pi^2\hbar^2/(2mL^2)$ with probability $|c_1|^2 = 1/5 = 20\%$.

$E_2 = 4\pi^2\hbar^2/(2mL^2)$ with probability $|c_2|^2 = 4/5 = 80\%$.

(d) Compute $\langle E \rangle$.

$$\langle E\rangle = \frac{1}{5}E_1 + \frac{4}{5}E_2 = \frac{1}{5}\cdot\frac{\pi^2\hbar^2}{2mL^2} + \frac{4}{5}\cdot\frac{4\pi^2\hbar^2}{2mL^2} = \frac{17\pi^2\hbar^2}{10mL^2}.$$

Note that $\langle E \rangle$ is not an eigenvalue — it is between $E_1$ and $E_2$. You will never measure this value in a single experiment. Each measurement gives either $E_1$ or $E_2$; the average over many experiments is $\langle E \rangle$.

⚠️ Common Misconception. Students often confuse the expectation value with a single measurement result. The expectation value is the statistical average over many identically prepared experiments. A single measurement of energy will always yield an eigenvalue — never anything in between.

🔄 Check Your Understanding. 1. Is $|\psi(x,t)|^2$ time-dependent for the superposition in Example 2.5? (Yes! Show this explicitly.) 2. What is the period of the time-dependent probability oscillation? 3. If, after measuring energy and finding $E_2$, you immediately measure energy again, what do you get? Why?

2.9 What the Wave Function Is (and Isn't): Interpretive Honesty

We end this chapter with the hardest question of all: What IS the wave function?

This is not a question with a universally agreed-upon answer. After nearly a century of debate, physicists remain divided. We owe you honesty about this, and we owe you a clear picture of what we know for certain versus what remains contested.

What We Know for Certain

These facts are not interpretation-dependent. They are true regardless of which interpretation of quantum mechanics you prefer:

1. The wave function is a complete description. Given $\psi(x,t)$, you can calculate the probability of any measurement outcome. There is no hidden variable that $\psi$ fails to capture — at least, not in standard quantum mechanics. (Hidden-variable theories are a separate topic; see Ch 24.)

2. The wave function evolves deterministically. Between measurements, $\psi$ obeys the Schrödinger equation, which is completely deterministic.

3. The Born rule works. $|\psi|^2$ gives the correct probability distribution for position measurements, and the analogous rules for other observables (Ch 6) give correct predictions. This has been verified by every experiment in the history of quantum mechanics, without a single exception.

4. Interference is real. The superposition of wave functions produces measurable interference effects. This is not a matter of interpretation — it is an experimental fact.

5. The wave function is complex-valued, and the complex phase matters. The relative phase between components of a superposition affects interference and hence affects measurable probabilities. This is not a mathematical convenience — it is physical.

What Remains Debated

⚖️ Interpretation. The following questions are genuinely open:

Is $\psi$ real or epistemic? The "ontic" view says $\psi$ is a real feature of the physical world — as real as an electric field. The "epistemic" view says $\psi$ represents our knowledge of the system, not the system itself. Both views have sophisticated defenders. The PBR theorem (2012) places strong constraints on epistemic interpretations, but the debate continues.

What happens during measurement? The standard ("Copenhagen") interpretation says the wave function collapses instantaneously upon measurement, but does not explain how or why. The many-worlds interpretation says there is no collapse — instead, the universe splits into branches. The Bohmian interpretation says particles have definite positions at all times, guided by $\psi$. The decoherence program (Ch 33) explains why measurement outcomes look classical, but does not fully resolve the measurement problem.

Does the wave function exist for a single system, or only for an ensemble? The ensemble interpretation says $\psi$ describes a statistical ensemble of identically prepared systems, not a single electron. Modern experiments with single trapped ions and single photons have made this position harder to maintain, but it still has adherents.

Our Pragmatic Position

In this textbook, we adopt the following stance:

1. We treat $\psi$ as a real, objective description of the quantum state. This is the most natural interpretation for working physicists and the one that makes the mathematics most transparent.

2. We use "collapse" as a computational tool without claiming it is a physical process. When we say "after measurement, the state collapses to $\phi_n$," we mean: use $\phi_n$ for all subsequent calculations.

3. We flag interpretive questions honestly whenever they arise, and we devote full chapters to them (Ch 24, Ch 28) when we have the mathematical tools to discuss them properly.

4. We encourage you to form your own view — but only after you have mastered the formalism. It is premature to argue about what $\psi$ "means" until you can calculate what it predicts.

💡 Key Insight. The remarkable thing about quantum mechanics is that physicists who disagree completely about what $\psi$ is agree completely about what $\psi$ predicts. Every interpretation gives the same experimental predictions (if they did not, we could settle the debate by experiment). This means you can learn to use quantum mechanics perfectly well without resolving the interpretive questions. But you should know the questions exist, and you should find the lack of consensus intellectually honest rather than troubling.

The Measurement Problem: A Preview

Here is the core puzzle, stated as sharply as we can at this stage:

The Schrödinger equation is linear and deterministic. It does not contain the word "measurement" or "observer." It does not explain why a single definite outcome occurs when you look at a detector. If you model the detector as a quantum system (which it is), and the observer as a quantum system (which they are), and apply the Schrödinger equation to the combined system, you get a superposition of "detector reads UP, observer sees UP" and "detector reads DOWN, observer sees DOWN." You do not get a single outcome.

Yet we always experience a single outcome.

This is the measurement problem. We will return to it in Chapter 6 (where we formalize measurement as a postulate), Chapter 23 (where density matrices let us describe incomplete knowledge), Chapter 24 (where Bell's theorem constrains the possibilities), and Chapter 28 (where we give it the full treatment it deserves).

For now, the practical rule is: compute with the Schrödinger equation, and extract predictions with the Born rule. This works. It has always worked. It will always work. The question of why it works is one of the deepest in all of science.

🔗 Connection. Chapter 3 begins applying these rules to real problems. We solve the infinite square well, the finite square well, and the potential barrier (tunneling). For the first time, you will see quantization and tunneling emerge directly from the mathematics — not as postulates, but as consequences of the Schrödinger equation plus boundary conditions.

Chapter Summary

Let us collect the key results:

The wave function $\psi(x,t)$ is a complex-valued function that completely describes the quantum state of a particle.

The Born rule states that $|\psi(x,t)|^2\,dx$ is the probability of finding the particle between $x$ and $x+dx$.

Normalization: $\int_{-\infty}^{\infty}|\psi|^2\,dx = 1$.

The time-dependent Schrödinger equation:

$$i\hbar\frac{\partial\psi}{\partial t} = \hat{H}\psi = \left[-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial x^2} + V(x)\right]\psi.$$

The time-independent Schrödinger equation (for time-independent potentials):

$$\hat{H}\phi = E\phi.$$

Stationary states: $\psi(x,t) = \phi(x)e^{-iEt/\hbar}$, with time-independent $|\psi|^2$.

Operators: Observables correspond to operators. The momentum operator is $\hat{p} = -i\hbar\,\partial/\partial x$.

Expectation values: $\langle A\rangle = \int\psi^*\hat{A}\psi\,dx$.

Probability current: $j = \frac{\hbar}{m}\text{Im}(\psi^*\partial\psi/\partial x)$, satisfying $\partial\rho/\partial t + \partial j/\partial x = 0$.

Superposition: If $\psi_1$ and $\psi_2$ are solutions, so is $c_1\psi_1 + c_2\psi_2$. This is the fundamental quantum principle.

Physical constraints on $\psi$: continuous, single-valued, square-integrable, with continuous first derivative (where $V$ is finite).

Project Checkpoint: Quantum Simulation Toolkit v0.2

In Chapter 1, you set up the toolkit repository and created a constants module. Now we add the core Wavefunction class.

New module: wavefunction.py

Your Wavefunction class should:

1. Accept a function $\psi(x)$ and a spatial grid as input.
2. Compute the normalization integral numerically (trapezoidal rule).
3. Normalize the wave function.
4. Compute $|\psi|^2$ (probability density).
5. Compute expectation values $\langle x \rangle$, $\langle x^2 \rangle$, and $\langle p \rangle$ (the last using a finite-difference approximation to $d\psi/dx$).
6. Plot $\psi(x)$ (real and imaginary parts) and $|\psi(x)|^2$.

See code/project-checkpoint.py for the implementation.

Test your class on: - The Gaussian $\psi(x) = Ce^{-x^2}$ — verify that $\langle x\rangle = 0$, $\sigma_x = 1/2$. - The function $\psi(x) = C\,x\,e^{-x^2/2}$ — verify that $\langle x\rangle = 0$ (by symmetry) but $\langle x^2\rangle \neq 0$. - A superposition $\psi(x) = C[e^{-(x-2)^2} + e^{-(x+2)^2}]$ — verify the interference term in $|\psi|^2$.

🔗 Connection. In Chapter 3, we will add a finite-difference Schrödinger solver to the toolkit, allowing you to compute energy eigenvalues and eigenfunctions numerically for arbitrary potentials.

Looking Ahead

With the rules of the game established, we are ready to play. Chapter 3 applies the Schrödinger equation to the simplest possible systems: the infinite square well, the finite square well, and the potential barrier. These "toy models" are not trivial — they teach us the essential physics of quantization, tunneling, and scattering that recur throughout quantum mechanics.

The infinite square well will give us our first complete set of energy eigenstates and demonstrate how to expand an arbitrary initial condition as a superposition. The finite well will show us that quantum particles "leak" into classically forbidden regions. And the barrier will reveal tunneling — a phenomenon with no classical analogue that powers technologies from tunnel diodes to scanning tunneling microscopes.

But first, test your understanding with the exercises that follow.

Next: Chapter 3 — Solving the Schrödinger Equation: Wells, Barriers, and the Art of Quantum Problem-Solving

Previous: Chapter 1 — The Quantum Revolution: When Classical Physics Broke Down