Skip to content

These exercise sets give a structured practice path through canonical wave mechanics. They are not meant to replace the exercises on individual model pages; they gather recurring techniques and make the progression explicit.

Before starting a set, read Worked Examples Index, skim Problem-Solving Patterns, and keep How to Solve a Wave-Mechanics Problem nearby. Each problem should be solved with stated assumptions, boundary conditions, normalization convention, and at least one check.

  1. For the infinite square well on 0<x<L0\lt x\lt L, verify that ψn(x)=2/Lsin⁡(nπx/L)\psi_n(x)=\sqrt{2/L}\sin(n\pi x/L) is normalized.
Solution

Use

∫0Lsin⁡2nπxL dx=L2.\int_0^L \sin^2\frac{n\pi x}{L}\,dx=\frac{L}{2}.

Then

∫0L∣ψn(x)∣2 dx=2L∫0Lsin⁡2nπxL dx=1.\int_0^L \lvert\psi_n(x)\rvert^2\,dx =\frac{2}{L}\int_0^L \sin^2\frac{n\pi x}{L}\,dx =1.
  1. A wavefunction on the line has probability density ρ(x)=Ae−2κ∣x∣\rho(x)=Ae^{-2\kappa\lvert x\rvert} with κ>0\kappa\gt0. Find AA.

  2. Explain why ∣ψ(x)∣2\lvert\psi(x)\rvert^2 is a probability density, not a probability, and state what must be integrated to get a probability on an interval.

  3. In spherical coordinates, explain why a radial probability density is not just ∣R(r)∣2\lvert R(r)\rvert^2.

  4. Starting from ∫dx ∣x⟩⟨x∣=I^\int dx\,\lvert x\rangle\langle x\rvert=\hat I, derive the position-space inner product ⟨ϕ∣ψ⟩=∫dx ϕ∗(x)ψ(x)\langle\phi\vert\psi\rangle=\int dx\,\phi^*(x)\psi(x).

Checks for the other problems
  1. 1=2A∫0∞e−2κxdx=A/κ1=2A\int_0^\infty e^{-2\kappa x}dx=A/\kappa, so A=κA=\kappa.
  2. P(a<x<b)=∫ab∣ψ(x)∣2dxP(a<x<b)=\int_a^b|\psi(x)|^2dx; the density has units of inverse length.
  3. With a normalized angular function, the shell probability is r2∣R(r)∣2drr^2|R(r)|^2dr because d3r=r2dr dΩd^3r=r^2dr\,d\Omega.
  4. Insert the identity: ⟨ϕ∣I∣ψ⟩=∫dx ⟨ϕ∣x⟩⟨x∣ψ⟩=∫dx ϕ∗ψ\langle\phi|I|\psi\rangle=\int dx\,\langle\phi|x\rangle\langle x|\psi\rangle=\int dx\,\phi^*\psi.

Good companion pages: Wavefunctions and Probability Density, Normalization Conventions, Coordinate Representation.

Set 2: Boundary Conditions and Bound States

Section titled “Set 2: Boundary Conditions and Bound States”
  1. For an infinite square well, show that the boundary conditions ψ(0)=ψ(L)=0\psi(0)=\psi(L)=0 force k=nπ/Lk=n\pi/L.

  2. Explain why the finite square well does not use ψ=0\psi=0 at the edges.

  3. Derive the derivative jump condition for V(x)=λδ(x)V(x)=\lambda\delta(x) by integrating the stationary Schrödinger equation across a small interval around the origin.

Solution

The stationary equation is

−ℏ22mψ′′(x)+λδ(x)ψ(x)=Eψ(x).-\frac{\hbar^2}{2m}\psi''(x) +\lambda\delta(x)\psi(x) =E\psi(x).

Integrate from −ϵ-\epsilon to ϵ\epsilon:

−ℏ22m[ψ′(ϵ)−ψ′(−ϵ)]+λψ(0)=E∫−ϵϵψ(x) dx.-\frac{\hbar^2}{2m} \left[\psi'(\epsilon)-\psi'(-\epsilon)\right] +\lambda\psi(0) =E\int_{-\epsilon}^{\epsilon}\psi(x)\,dx.

The right side vanishes as ϵ→0\epsilon\to0 for a finite wavefunction. Therefore

ψ′(0+)−ψ′(0−)=2mλℏ2ψ(0).\psi'(0^+)-\psi'(0^-) =\frac{2m\lambda}{\hbar^2}\psi(0).
  1. For a symmetric finite well, use parity to reduce the matching problem to even and odd sectors. State which matching condition is applied at the center and which is applied at the well edge.

  2. Build a finite-difference Hamiltonian for a deep finite well. Vary the outside barrier height and verify that the lowest energies approach the infinite-well values.

Checks for the other problems
  1. The first wall removes the cosine solution; the second requires sin⁡(kL)=0\sin(kL)=0. Nonzero states have k=nπ/Lk=n\pi/L, n≥1n\ge1.
  2. Finite walls admit decaying tails, matched continuously with their first derivatives to the interior solution.
  3. At the regular center, use ψ′(0)=0\psi'(0)=0 for even states and ψ(0)=0\psi(0)=0 for odd states. At the right edge, match the logarithmic derivative to −κ-\kappa from the exterior decaying solution.
  4. Hold the interior potential at zero while raising the exterior barrier. Check grid refinement and domain enlargement separately from the finite-barrier limit; otherwise grid error can conceal the approach to the hard-wall energies.

Good companion pages: Boundary Conditions, Boundary Conditions Table, Finite Square Well, Delta Function Potential.

  1. For a plane wave AeikxAe^{ikx} with k>0k\gt0, compute the one-dimensional probability current.
Solution

Use

j=ℏmIm⁡(ψ∗dψdx).j=\frac{\hbar}{m} \operatorname{Im}\left(\psi^*\frac{d\psi}{dx}\right).

For ψ=Aeikx\psi=Ae^{ikx},

ψ∗dψdx=A∗e−ikx ikAeikx=ik∣A∣2.\psi^*\frac{d\psi}{dx} =A^*e^{-ikx}\,ikAe^{ikx} =ik\lvert A\rvert^2.

Thus

j=ℏkm∣A∣2.j=\frac{\hbar k}{m}\lvert A\rvert^2.
  1. Explain why reflection and transmission coefficients are current ratios rather than merely amplitude-squared ratios.

  2. For a potential step with E>V0E\gt V_0, derive r=(k−q)/(k+q)r=(k-q)/(k+q) and identify why TT includes the factor q/kq/k.

  3. In the opaque-barrier approximation, TT is proportional to e−2κae^{-2\kappa a}. What happens to TT when the barrier width aa is increased by Δa\Delta a?

  4. Numerically propagate a Gaussian packet toward a rectangular barrier. Check that the total norm is conserved and estimate reflected plus transmitted probability after the packets separate.

Checks for the other problems
  1. A plane wave’s current is its density times velocity; different velocities require different factors even at equal amplitude.
  2. Continuity gives 1+r=t1+r=t and k(1−r)=qtk(1-r)=qt, with k=2mE/ℏk=\sqrt{2mE}/\hbar and q=2m(E−V0)/ℏq=\sqrt{2m(E-V_0)}/\hbar. Thus r=(k−q)/(k+q)r=(k-q)/(k+q), t=2k/(k+q)t=2k/(k+q), and T=4kq/(k+q)2T=4kq/(k+q)^2.
  3. At fixed energy and height, T(a+Δa)/T(a)≃e−2κΔaT(a+\Delta a)/T(a)\simeq e^{-2\kappa\Delta a} in the opaque regime.
  4. Before boundary reflections occur, integrate density in the left, barrier, and right regions. Their sum must stay near one; only after central probability becomes negligible may the left and right integrals alone be identified with final reflection and transmission probabilities.

Good companion pages: Probability Current, Potential Step, Reflection and Transmission Coefficients, Rectangular Barrier Tunneling.

Set 4: Oscillators and Finite-Dimensional Models

Section titled “Set 4: Oscillators and Finite-Dimensional Models”
  1. If ω\omega doubles in the harmonic oscillator, how do the level spacing ℏω\hbar\omega and oscillator length ℓ=ℏ/(mω)\ell=\sqrt{\hbar/(m\omega)} change?

  2. Explain why the harmonic-oscillator ground-state energy is not zero even though the potential minimum is zero.

  3. Starting from a^∣0⟩=0\hat a\lvert0\rangle=0, derive the Gaussian ground-state wavefunction in position space.

  4. For ∣ψ⟩=(∣0⟩+3 ∣1⟩)/2\lvert\psi\rangle=(\lvert0\rangle+\sqrt3\,\lvert1\rangle)/2, list the possible oscillator energies and their probabilities.

Solution

The coefficients are c0=1/2c_0=1/2 and c1=3/2c_1=\sqrt3/2. Therefore the possible energies are

E0=12ℏω,E1=32ℏω,E_0=\frac{1}{2}\hbar\omega, \qquad E_1=\frac{3}{2}\hbar\omega,

with probabilities

P(E0)=∣c0∣2=14,P(E1)=∣c1∣2=34.P(E_0)=\lvert c_0\rvert^2=\frac{1}{4}, \qquad P(E_1)=\lvert c_1\rvert^2=\frac{3}{4}.
  1. Explain why the oscillator is a bridge to field modes, while a two-level system is a bridge to qubits. Keep the answer at the level of Hilbert-space structure and spectra.
Checks for the other problems
  1. At fixed mass, the spacing doubles and the length becomes ℓ/2\ell/\sqrt2.
  2. Localizing at the potential minimum increases momentum uncertainty. Minimizing (Δp)2/(2m)+mω2(Δx)2/2(\Delta p)^2/(2m)+m\omega^2(\Delta x)^2/2 subject to ΔxΔp≥ℏ/2\Delta x\Delta p\ge\hbar/2 gives ℏω/2\hbar\omega/2.
  3. The equation is (x/ℓ+ℓ d/dx)ψ0=0(x/\ell+\ell\,d/dx)\psi_0=0, so ψ0=π−1/4ℓ−1/2e−x2/(2ℓ2)\psi_0=\pi^{-1/4}\ell^{-1/2}e^{-x^2/(2\ell^2)} after normalization.
  4. An oscillator has an unbounded occupation-number ladder; an ideal two-level system has a two-dimensional Hilbert space. A physical qubit device needs an isolated, controllable pair of levels, even when additional device levels exist.

Good companion pages: Quantum Harmonic Oscillator, Ladder-Operator Solution, Number States, Two-Level Systems.

Set 5: Three Dimensions, Angular Structure, and Magnetic Fields

Section titled “Set 5: Three Dimensions, Angular Structure, and Magnetic Fields”
  1. For a three-dimensional box with equal sides LL, find all positive integer triples (nx,ny,nz)(n_x,n_y,n_z) with nx2+ny2+nz2=6n_x^2+n_y^2+n_z^2=6 and state the degeneracy.

  2. Explain why mm degeneracy is expected in a rotationally invariant central potential, while hydrogen’s extra ℓ\ell degeneracy is special to the Coulomb potential.

  3. Starting from u(r)=rR(r)u(r)=rR(r), show how the radial normalization ∫0∞∣u(r)∣2 dr=1\int_0^\infty\lvert u(r)\rvert^2\,dr=1 follows from the three-dimensional normalization when the angular part is normalized.

  4. For the rigid rotor, show that the energy difference between adjacent levels is

EJ+1−EJ=2B(J+1),B=ℏ22I.E_{J+1}-E_J=2B(J+1), \qquad B=\frac{\hbar^2}{2I}.
  1. In Landau gauge, use guiding-center spacing to estimate the degeneracy of one Landau level in a rectangle of area A=LxLyA=L_xL_y.
Solution

For a particle with charge magnitude ∣q∣\lvert q\rvert in a uniform magnetic field, the magnetic length is

ℓB=ℏ∣q∣B.\ell_B=\sqrt{\frac{\hbar}{\lvert q\rvert B}}.

Periodic boundary conditions in the yy direction give spacing Δky=2π/Ly\Delta k_y=2\pi/L_y. The guiding center spacing is therefore

Δx0=ℓB2Δky=2πℓB2Ly.\Delta x_0=\ell_B^2\Delta k_y =\frac{2\pi\ell_B^2}{L_y}.

The number of allowed centers across width LxL_x is approximately

N=LxΔx0=A2πℓB2=∣q∣BAh.N=\frac{L_x}{\Delta x_0} =\frac{A}{2\pi\ell_B^2} =\frac{\lvert q\rvert BA}{h}.

Boundary corrections are ignored in this bulk estimate.

Checks for the other problems
  1. The triples are (1,1,2)(1,1,2), (1,2,1)(1,2,1), and (2,1,1)(2,1,1), giving degeneracy three.
  2. Rotations relate the 2ℓ+12\ell+1 states within a fixed angular-momentum multiplet. They do not require equal energies for different ℓ\ell; the Coulomb potential has an additional dynamical symmetry.
  3. Substitute ∣u∣2=r2∣R∣2|u|^2=r^2|R|^2 in the radial integral after the angular integral has given one.
  4. Subtraction gives B[(J+1)(J+2)−J(J+1)]=2B(J+1)B[(J+1)(J+2)-J(J+1)]=2B(J+1).

Good companion pages: Three-Dimensional Box, Radial Schrödinger Equation, Degeneracy of the Hydrogen Atom, Rigid Rotor, Landau Levels.

For each set, solve one direct calculation, one interpretation question, and one derivation before moving on. If a derivation fails, write down the exact condition that blocked it: a boundary condition, a normalization convention, a missing identity, or a physical assumption.

For computational problems, record the grid, time step, boundary conditions, norm error, and comparison target. A numerical answer without these checks is not yet reproducible.

  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • E. Merzbacher, Quantum Mechanics, 3rd ed., Wiley, 1998.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.