STEMviz
Physics · Quantum Mechanics and Modern Physics

The Schrodinger Equation (Time-Dependent and Time-Independent)

The time-dependent Schrödinger equation iℏ ∂tΨ=H^Ψi\hbar\,\partial_t\Psi=\hat H\Psi is the postulated dynamical law: first order in time (deterministic), linear (superposition), norm-preserving. For time-independent VV, separation of variables gives stationary states Ψ=ψ(x)e−iEt/ℏ\Psi=\psi(x)e^{-iEt/\hbar} obeying the time-independent equation −ℏ22mψ′′+Vψ=Eψ-\tfrac{\hbar^2}{2m}\psi''+V\psi=E\psi; the general state is Ψ=∑ncnψne−iEnt/ℏ\Psi=\sum_n c_n\psi_n e^{-iE_n t/\hbar}. Bound-state energies are quantised; free states are non-normalisable plane waves that must be built into spreading wave packets with vg=ℏk0/mv_g=\hbar k_0/m.

The Schrödinger Equation (Time-Dependent and Time-Independent)

Chapter 2 said that a particle's state is a wave function and that measurements are statistical. It did not say how the state changes with time. That is the job of the Schrödinger equation, the dynamical law of non-relativistic quantum mechanics. Like Newton's second law it is a postulate — not derived from anything more basic — justified because its predictions match experiment to extraordinary precision. This chapter states it, extracts the time-independent Schrödinger equation by separation of variables, explains why its solutions ("stationary states") are the backbone of every bound-state calculation in Unit D, and treats the one unbound case that is fully solvable by hand: the free particle and its wave packets.

The organising idea is that the Hamiltonian operator H^\hat H does two things at once. It generates time evolution, through iℏ ∂tΨ=H^Ψi\hbar\,\partial_t\Psi = \hat H\Psi; and its eigenvalue problem H^ψ=Eψ\hat H\psi = E\psi picks out the special states of definite energy, in terms of which every other state is expanded. Solving a quantum system means finding the spectrum and eigenfunctions of its H^\hat H; time evolution is then automatic, each eigenstate merely acquiring a phase e−iEt/ℏe^{-iEt/\hbar}.

How to read this chapter. Two equations do all the work. The TDSE iℏ ∂tΨ=−ℏ22m∂x2Ψ+VΨi\hbar\,\partial_t\Psi = -\dfrac{\hbar^2}{2m}\partial_x^2\Psi + V\Psi governs how any state evolves. Separation of variables reduces it, for time-independent VV, to the TISE −ℏ22mψ′′+Vψ=Eψ-\dfrac{\hbar^2}{2m}\psi'' + V\psi = E\psi, an eigenvalue problem whose bound-state solutions exist only for discrete EE. Chapters 6–8 are three worked instances (well, oscillator, angular momentum); this chapter is the machinery they share.

1. The Time-Dependent Schrödinger Equation

For a particle of mass mm in a potential V(x,t)V(x,t), the wave function obeys

iℏ ∂Ψ∂t=−ℏ22m ∂2Ψ∂x2+V(x,t) Ψ  ≡  H^Ψi\hbar\,\frac{\partial\Psi}{\partial t} = -\frac{\hbar^2}{2m}\,\frac{\partial^2\Psi}{\partial x^2} + V(x,t)\,\Psi \;\equiv\; \hat H\Psi

The right-hand side defines the Hamiltonian operator H^=p^2/2m+V(x^)\hat H = \hat p^2/2m + V(\hat x), the quantum version of the classical energy p2/2m+Vp^2/2m + V. Key features:

  • It is first order in time: Ψ(x,0)\Psi(x,0) everywhere determines Ψ(x,t)\Psi(x,t) for all later tt. The evolution is deterministic — the randomness in quantum mechanics is entirely in measurement, not in the dynamics.
  • It is linear: if Ψ1,Ψ2\Psi_1,\Psi_2 are solutions so is c1Ψ1+c2Ψ2c_1\Psi_1 + c_2\Psi_2. This is what makes superposition consistent.
  • The explicit ii forces Ψ\Psi to be complex and makes the evolution unitary (norm-preserving), which is why ∫∣Ψ∣2 dx\int|\Psi|^2\,dx stays equal to 11.

The equation cannot be derived from classical mechanics. It can be motivated: for a free particle, demanding that ei(kx−ωt)e^{i(kx-\omega t)} be a solution with E=ℏωE = \hbar\omega and p=ℏkp = \hbar k related by E=p2/2mE = p^2/2m forces exactly the free (V=0V=0) form. Adding VV is the natural guess, and it works.

2. Separation of Variables and the Time-Independent Equation

When V=V(x)V = V(x) does not depend on time, look for product solutions Ψ(x,t)=ψ(x) φ(t)\Psi(x,t) = \psi(x)\,\varphi(t). Substituting and dividing by ψφ\psi\varphi:

iℏ φ′(t)φ(t)=−ℏ22m ψ′′(x)ψ(x)+V(x)i\hbar\,\frac{\varphi'(t)}{\varphi(t)} = -\frac{\hbar^2}{2m}\,\frac{\psi''(x)}{\psi(x)} + V(x)

The left side depends only on tt, the right only on xx; both must equal a constant, which has the dimensions of energy and is called EE. This splits the problem into

iℏ φ′(t)=E φ(t)  ⇒  φ(t)=e−iEt/ℏi\hbar\,\varphi'(t) = E\,\varphi(t) \;\Rightarrow\; \varphi(t) = e^{-iEt/\hbar}   −ℏ22m d2ψdx2+V(x) ψ=E ψ  (time-independent Schro¨dinger equation)\boxed{\;-\frac{\hbar^2}{2m}\,\frac{d^2\psi}{dx^2} + V(x)\,\psi = E\,\psi\;}\qquad\text{(time-independent Schrödinger equation)}

The TISE is an eigenvalue equation H^ψ=Eψ\hat H\psi = E\psi: find the functions ψ\psi and numbers EE for which applying H^\hat H just rescales ψ\psi. The full separated solution is

ΨE(x,t)=ψE(x) e−iEt/ℏ\Psi_E(x,t) = \psi_E(x)\,e^{-iEt/\hbar}

3. Stationary States and Why They Are Called That

A separated solution ΨE=ψE(x)e−iEt/ℏ\Psi_E = \psi_E(x)e^{-iEt/\hbar} is a stationary state. Despite the oscillating phase, everything measurable about it is time-independent:

  • ∣ΨE∣2=∣ψE∣2|\Psi_E|^2 = |\psi_E|^2 — the probability density does not move.
  • ⟨Q⟩=∫ψE∗ Q^ ψE dx\langle Q\rangle = \int\psi_E^*\,\hat Q\,\psi_E\,dx for any operator Q^\hat Q with no explicit time dependence — the phases e±iEt/ℏe^{\pm iEt/\hbar} cancel.
  • The energy is sharp: ⟨H⟩=E\langle H\rangle = E and ⟨H2⟩=E2\langle H^2\rangle = E^2, so ΔE=0\Delta E = 0. A measurement of energy in a stationary state returns EE with certainty.
  • The probability current is constant in time (and, in 1D, constant in xx).

"Stationary" refers to the observables, not the wave function. These states are the quantum analogue of the normal modes of a vibrating string: each oscillates at its own frequency ω=E/ℏ\omega = E/\hbar and nothing about its shape changes.

4. The General Solution as a Superposition

The TISE typically has many solutions ψn(x)\psi_n(x) with energies EnE_n (discrete for bound states; see Section 5). Because the TDSE is linear, the general solution is a superposition of stationary states:

Ψ(x,t)=∑ncn ψn(x) e−iEnt/ℏ\Psi(x,t) = \sum_n c_n\,\psi_n(x)\,e^{-iE_n t/\hbar}

The constants cnc_n are fixed by the initial condition: at t=0t = 0, Ψ(x,0)=∑ncnψn(x)\Psi(x,0) = \sum_n c_n\psi_n(x), so by orthonormality of the ψn\psi_n,

cn=∫ψn∗(x) Ψ(x,0) dxc_n = \int \psi_n^*(x)\,\Psi(x,0)\,dx

This is the entire method. To evolve any state: (1) expand the initial Ψ(x,0)\Psi(x,0) in energy eigenstates to get the cnc_n; (2) attach the phase factor e−iEnt/ℏe^{-iE_n t/\hbar} to each term; (3) resum. The hard part is always step (1) — finding the spectrum and the cnc_n. Once done, time evolution is trivial.

Worked example — a two-term superposition in a box. Take Ψ(x,0)=12[ψ1(x)+ψ2(x)]\Psi(x,0) = \tfrac{1}{\sqrt2}[\psi_1(x) + \psi_2(x)] in an infinite well. Then Ψ(x,t)=12[ψ1e−iE1t/ℏ+ψ2e−iE2t/ℏ]\Psi(x,t) = \tfrac{1}{\sqrt2}[\psi_1 e^{-iE_1 t/\hbar} + \psi_2 e^{-iE_2 t/\hbar}], and

∣Ψ(x,t)∣2=12[ψ12+ψ22+2ψ1ψ2cos⁡ ⁣(E2−E1ℏ t)]|\Psi(x,t)|^2 = \tfrac12\left[\psi_1^2 + \psi_2^2 + 2\psi_1\psi_2\cos\!\left(\frac{E_2 - E_1}{\hbar}\,t\right)\right]

The probability density sloshes from side to side at the angular frequency ω21=(E2−E1)/ℏ\omega_{21} = (E_2 - E_1)/\hbar — a genuinely time-dependent state, unlike either eigenstate alone. Its energy is not sharp: a measurement gives E1E_1 or E2E_2, each with probability 12\tfrac12.

5. Energy Quantisation for Bound States

For a potential that confines the particle (e.g. V→∞V\to\infty as ∣x∣→∞|x|\to\infty, or a well deep enough to trap it), the requirement that ψ\psi be normalisable — ψ→0\psi\to 0 at infinity, ψ\psi and ψ′\psi' continuous — cannot be met for arbitrary EE. Generic EE produces a solution that diverges at large ∣x∣|x|; only a discrete set E1<E2<E3<⋯E_1 < E_2 < E_3 < \cdots gives a decaying, square-integrable ψn\psi_n. Bound-state energies are quantised. Two universal features:

  • There is a lowest energy E1>Vmin⁡E_1 > V_{\min} — the ground state. Confinement plus Δx Δp≳ℏ/2\Delta x\,\Delta p \gtrsim \hbar/2 forbids a state of zero kinetic energy (zero-point energy).
  • The nnth state has n−1n-1 nodes (interior zeros), and the ψn\psi_n can be chosen real and are mutually orthogonal, ∫ψm∗ψn dx=δmn\int\psi_m^*\psi_n\,dx = \delta_{mn}.

For energies above the top of the potential the particle is unbound: the spectrum is continuous and the eigenfunctions are non-normalisable scattering states (Chapter 6).

6. The Free Particle and Wave Packets

With V=0V = 0 the TISE is −ℏ22mψ′′=Eψ-\dfrac{\hbar^2}{2m}\psi'' = E\psi, solved by ψk(x)=eikx\psi_k(x) = e^{ikx} with E=ℏ2k22mE = \dfrac{\hbar^2 k^2}{2m}, so the stationary solutions are

Ψk(x,t)=ei(kx−ωt),ω=ℏk22m\Psi_k(x,t) = e^{i(kx - \omega t)}, \qquad \omega = \frac{\hbar k^2}{2m}

These are not normalisable (∣Ψk∣2=1|\Psi_k|^2 = 1 everywhere), so no free particle is in a single plane-wave state. Physical states are wave packets — continuous superpositions,

Ψ(x,t)=12π∫−∞∞ϕ(k) ei(kx−ω(k)t) dk\Psi(x,t) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^{\infty} \phi(k)\,e^{i(kx - \omega(k)t)}\,dk

with ϕ(k)\phi(k) the Fourier transform of Ψ(x,0)\Psi(x,0). A packet peaked around k0k_0 with spread Δk\Delta k describes a particle localised to Δx∼1/Δk\Delta x \sim 1/\Delta k around a moving centre, with Δx Δp≳ℏ/2\Delta x\,\Delta p \gtrsim \hbar/2.

7. Phase Velocity, Group Velocity, and Spreading

For the free particle the dispersion relation is ω(k)=ℏk2/2m\omega(k) = \hbar k^2/2m, which is nonlinear in kk. Two velocities follow:

vphase=ωk=ℏk2m,vgroup=dωdk∣k0=ℏk0mv_{\text{phase}} = \frac{\omega}{k} = \frac{\hbar k}{2m}, \qquad v_{\text{group}} = \frac{d\omega}{dk}\bigg|_{k_0} = \frac{\hbar k_0}{m}

The group velocity is the speed of the packet's envelope, and it equals the classical particle speed p0/mp_0/m — as it must, matching d⟨x⟩/dt=⟨p⟩/md\langle x\rangle/dt = \langle p\rangle/m from Chapter 2. The phase velocity is half of it and has no direct physical meaning for a particle.

Because ω′′(k)=ℏ/m≠0\omega''(k) = \hbar/m \ne 0, different Fourier components travel at different group speeds, so the packet spreads as it moves: an initial width Δx0\Delta x_0 grows, for large tt, as Δx(t)≈ℏt2m Δx0\Delta x(t) \approx \dfrac{\hbar t}{2m\,\Delta x_0}. A tightly localised particle disperses fastest. This spreading is a real, observed effect and is absent only for special potentials (the harmonic oscillator, whose linear dispersion in nn keeps packets coherent).

8. Hermiticity, Conservation, and the Role of the Energy Basis

Because H^\hat H is Hermitian:

  • its eigenvalues EnE_n are real (energies are measurable);
  • its eigenfunctions form a complete orthonormal set, so the expansion Ψ=∑ncnψne−iEnt/ℏ\Psi = \sum_n c_n\psi_n e^{-iE_nt/\hbar} is always possible;
  • time evolution is unitary: ∑n∣cn∣2=1\sum_n |c_n|^2 = 1 is preserved, since each ∣cn∣2|c_n|^2 is separately constant (the phases have unit modulus). Probabilities of energy outcomes never change in time — energy is conserved in the statistical sense.

More generally, for any operator Q^\hat Q with no explicit time dependence,

d⟨Q⟩dt=iℏ ⟨[H^,Q^]⟩\frac{d\langle Q\rangle}{dt} = \frac{i}{\hbar}\,\big\langle[\hat H,\hat Q]\big\rangle

so ⟨Q⟩\langle Q\rangle is a constant of the motion precisely when [H^,Q^]=0[\hat H,\hat Q] = 0. Energy itself (Q^=H^\hat Q = \hat H) is always conserved for a time-independent Hamiltonian.

9. Common Pitfalls

  • Calling ψE(x)e−iEt/ℏ\psi_E(x)e^{-iEt/\hbar} "time-dependent physics". The phase is unobservable; every measurable quantity in a stationary state is constant. It is a state of definite energy, nothing about it evolves.
  • Using the TISE when VV depends on time. Separation of variables fails if V=V(x,t)V = V(x,t); you must solve the full TDSE (or use time-dependent perturbation theory, Chapter 12).
  • Forgetting the phase factors when evolving a superposition. Ψ(x,t)=∑cnψne−iEnt/ℏ\Psi(x,t) = \sum c_n\psi_n e^{-iE_n t/\hbar} — each term carries its own EnE_n. Dropping them (or using a common phase) loses all the dynamics.
  • Normalising a plane wave. eikxe^{ikx} has infinite norm; build a wave packet. Quantities like ⟨x⟩\langle x\rangle are undefined for a bare plane wave.
  • Confusing phase and group velocity. The particle moves at vg=ℏk0/m=p0/mv_g = \hbar k_0/m = p_0/m; vp=ℏk/2mv_p = \hbar k/2m is not the particle speed.
  • Assuming the energy is sharp in a general state. Only stationary states have ΔE=0\Delta E = 0. A superposition of eigenstates has a spread of energies with probabilities ∣cn∣2|c_n|^2.
  • Expecting a free wave packet to hold its shape. Free packets spread, because ω(k)\omega(k) is nonlinear. Only the harmonic oscillator (and other special systems) avoid dispersion.
  • Thinking quantisation is imposed by hand. It emerges: the normalisability/boundary conditions on the TISE admit solutions only for discrete EE.

Summary

  • TDSE: iℏ ∂tΨ=H^Ψi\hbar\,\partial_t\Psi = \hat H\Psi with H^=−ℏ22m∂x2+V\hat H = -\dfrac{\hbar^2}{2m}\partial_x^2 + V. First-order in time (deterministic), linear (superposition), complex/unitary (norm-preserving).
  • Separation (for V=V(x)V = V(x)): Ψ=ψ(x)e−iEt/ℏ\Psi = \psi(x)e^{-iEt/\hbar} with the TISE −ℏ22mψ′′+Vψ=Eψ-\dfrac{\hbar^2}{2m}\psi'' + V\psi = E\psi, an eigenvalue problem H^ψ=Eψ\hat H\psi = E\psi.
  • Stationary states: ∣Ψ∣2|\Psi|^2, expectation values and jj are time-independent; energy is sharp, ΔE=0\Delta E = 0.
  • General solution: Ψ(x,t)=∑ncnψn(x)e−iEnt/ℏ\Psi(x,t) = \sum_n c_n\psi_n(x)e^{-iE_n t/\hbar}, cn=∫ψn∗Ψ(x,0) dxc_n = \int\psi_n^*\Psi(x,0)\,dx. Evolve by expanding, phasing, resumming.
  • Bound states: normalisability quantises the energy, E1<E2<⋯E_1 < E_2 < \cdots; ground-state energy >Vmin⁡> V_{\min} (zero-point energy); ψn\psi_n has n−1n-1 nodes and ∫ψm∗ψn=δmn\int\psi_m^*\psi_n = \delta_{mn}.
  • Free particle: Ψk=ei(kx−ωt)\Psi_k = e^{i(kx-\omega t)}, ω=ℏk2/2m\omega = \hbar k^2/2m, non-normalisable; physical states are wave packets Ψ=12π∫ϕ(k)ei(kx−ωt) dk\Psi = \tfrac{1}{\sqrt{2\pi}}\int\phi(k)e^{i(kx-\omega t)}\,dk.
  • Velocities: vp=ℏk/2mv_p = \hbar k/2m, vg=ℏk0/m=p0/mv_g = \hbar k_0/m = p_0/m (the particle speed). Nonlinear ω(k)\omega(k) makes packets spread as Δx(t)∼ℏt/(2mΔx0)\Delta x(t) \sim \hbar t/(2m\Delta x_0).
  • Conservation: H^\hat H Hermitian ⇒\Rightarrow real EnE_n, complete basis, ∣cn∣2|c_n|^2 constant in time; d⟨Q⟩dt=iℏ⟨[H^,Q^]⟩\dfrac{d\langle Q\rangle}{dt} = \dfrac{i}{\hbar}\langle[\hat H,\hat Q]\rangle, so [H^,Q^]=0[\hat H,\hat Q] = 0 means ⟨Q⟩\langle Q\rangle is conserved.