Skip to content

One-Dimensional Problems

Attempt these from memory before reading.

  • What are the energy levels of the infinite square well, and what sets the spacing between them?
  • What is the lowest energy of the well, and why is it not zero?
  • What makes the harmonic oscillator soluble in closed form?
  • What does the parity of a state mean, and when is it a good quantum number?
  • What is the tunnelling probability’s dependence on barrier width?

Check your answers against the text. Prompt 2 is the one to sit with: a particle in a box cannot sit still, and that is a statement about the size of the box rather than about the particle.

Motivation: boundary conditions are the physics

Section titled “Motivation: boundary conditions are the physics”

Nothing in the Schrodinger equation knows that a well has walls. The equation inside the well is the same free-particle equation, and its solutions are sines and cosines with any wavelength. Every number that comes out of the infinite square well — the energies En=n2π2ℏ2/2mL2E_n = n^2\pi^2\hbar^2/2mL^2, the orthogonality of the states, the nodes — comes from the boundary conditions and nothing else.

That is worth stating plainly, because the standard presentation makes it look as though solving the equation produces the quantisation. It does not. Quantisation is a statement that only certain wavelengths fit, and “fit” is a boundary-condition word.

The zero-point energy is the same argument in disguise. Localising a particle to width LL costs momentum of order ℏ/L\hbar/L, so the kinetic energy is at least ℏ2/2mL2\hbar^2/2mL^2. The exact answer E1=π2ℏ2/2mL2E_1 = \pi^2\hbar^2/2mL^2 is that estimate with π2\pi^2 in the right place. Shrink the box and the ground state energy rises, which is why an electron confined to an atom has tens of electron-volts and a proton confined to a nucleus has millions.

Tunnelling is the same boundary-condition story in the other direction. Inside a barrier with E<V0E < V_0, the solutions are exponentials rather than oscillations — not because the equation changed, but because the sign of E−VE - V flipped and turned the differential equation from oscillatory to exponential. The decay constant κ\kappa is the only new object, and the transmission probability goes as e−2κae^{-2\kappa a}. Everything else is bookkeeping.

A particle of mass mm in a potential V(x)=0V(x) = 0 for 0<x<L0 \lt x \lt L and V(x)=∞V(x) = \infty otherwise.

Derivation. Inside the well, the time-independent Schrodinger equation is:

−ℏ22md2ϕdx2=Eϕ  ⟹  d2ϕdx2+k2ϕ=0-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} = E\phi \implies \frac{d^2\phi}{dx^2} + k^2\phi = 0

Where k=2mE/ℏk = \sqrt{2mE}/\hbar. The general solution is:

ϕ(x)=Asin⁡(kx)+Bcos⁡(kx)\phi(x) = A\sin(kx) + B\cos(kx)

Boundary conditions: ϕ(0)=ϕ(L)=0\phi(0) = \phi(L) = 0.

From ϕ(0)=0\phi(0) = 0: B=0B = 0, so ϕ(x)=Asin⁡(kx)\phi(x) = A\sin(kx).

From ϕ(L)=0\phi(L) = 0: sin⁡(kL)=0\sin(kL) = 0, which requires kL=nπkL = n\pi for n=1,2,3,…n = 1, 2, 3, \ldots

Therefore kn=nπ/Lk_n = n\pi/L and:

En=ℏ2kn22m=n2π2ℏ22mL2E_n = \frac{\hbar^2 k_n^2}{2m} = \frac{n^2\pi^2\hbar^2}{2mL^2}

Normalisation. ∫0L∣A∣2sin⁡2(nπx/L) dx=∣A∣2L/2=1\int_0^L |A|^2\sin^2(n\pi x/L)\,dx = |A|^2 L/2 = 1 Giving A=2/LA = \sqrt{2/L}.

Solutions:

ϕn(x)=2Lsin⁡(nπxL),En=n2π2ℏ22mL2,n=1,2,3,…\phi_n(x) = \sqrt{\frac{2}{L}}\sin\left(\frac{n\pi x}{L}\right), \quad E_n = \frac{n^2 \pi^2 \hbar^2}{2mL^2}, \quad n = 1, 2, 3, \ldots

Properties:

  • The ground state (n=1n = 1) has the lowest energy E1>0E_1 > 0 (zero-point energy).
  • Energy levels are not equally spaced; En∝n2E_n \propto n^2.
  • There are (n−1)(n - 1) nodes in the nn-th eigenstate.

V(x)=12mω2x2V(x) = \frac{1}{2}m\omega^2 x^2.

Define the ladder operators (creation and annihilation operators):

a^=mω2ℏ(x^+ip^mω),a^†=mω2ℏ(x^−ip^mω)\hat{a} = \sqrt{\frac{m\omega}{2\hbar}}\left(\hat{x} + \frac{i\hat{p}}{m\omega}\right), \quad \hat{a}^\dagger = \sqrt{\frac{m\omega}{2\hbar}}\left(\hat{x} - \frac{i\hat{p}}{m\omega}\right)

Commutation relation. Using [x^,p^]=iℏ[\hat{x}, \hat{p}] = i\hbar:

[a^,a^†]=mω2ℏ ⁣[x^+ip^mω, x^−ip^mω]=12ℏ(−i)(iℏ)+12ℏ(i)(−iℏ)=1[\hat{a}, \hat{a}^\dagger] = \frac{m\omega}{2\hbar}\!\left[\hat{x} + \frac{i\hat{p}}{m\omega},\, \hat{x} - \frac{i\hat{p}}{m\omega}\right] = \frac{1}{2\hbar}(-i)(i\hbar) + \frac{1}{2\hbar}(i)(-i\hbar) = 1

Inversion. From the definitions:

x^=ℏ2mω(a^+a^†),p^=−imωℏ2(a^−a^†)\hat{x} = \sqrt{\frac{\hbar}{2m\omega}}(\hat{a} + \hat{a}^\dagger), \quad \hat{p} = -i\sqrt{\frac{m\omega\hbar}{2}}(\hat{a} - \hat{a}^\dagger)

Hamiltonian in terms of ladder operators. Substituting into H^=p^2/(2m)+mω2x^2/2\hat{H} = \hat{p}^2/(2m) + m\omega^2\hat{x}^2/2:

H^=ℏω ⁣(a^†a^+12)\hat{H} = \hbar\omega\!\left(\hat{a}^\dagger\hat{a} + \frac{1}{2}\right)

Where we used a^a^†=[a^,a^†]+a^†a^=1+a^†a^\hat{a}\hat{a}^\dagger = [\hat{a}, \hat{a}^\dagger] + \hat{a}^\dagger\hat{a} = 1 + \hat{a}^\dagger\hat{a}.

Number operator. N^=a^†a^\hat{N} = \hat{a}^\dagger\hat{a} So H^=ℏω(N^+1/2)\hat{H} = \hbar\omega(\hat{N} + 1/2).

Proof that a^\hat{a} and a^†\hat{a}^\dagger lower and raise the energy. Compute [H^,a^][\hat{H}, \hat{a}]:

[H^,a^]=ℏω[a^†a^,a^]=ℏω(a^†[a^,a^]+[a^†,a^]a^)=−ℏω a^[\hat{H}, \hat{a}] = \hbar\omega[\hat{a}^\dagger\hat{a}, \hat{a}] = \hbar\omega(\hat{a}^\dagger[\hat{a}, \hat{a}] + [\hat{a}^\dagger, \hat{a}]\hat{a}) = -\hbar\omega\,\hat{a}

So H^a^∣n⟩=(En−ℏω)a^∣n⟩\hat{H}\hat{a}|n\rangle = (E_n - \hbar\omega)\hat{a}|n\rangle: a^\hat{a} lowers energy by ℏω\hbar\omega. Similarly, [H^,a^†]=+ℏω a^†[\hat{H}, \hat{a}^\dagger] = +\hbar\omega\,\hat{a}^\dagger.

Let ∣n⟩|n\rangle be an eigenstate with H^∣n⟩=En∣n⟩\hat{H}|n\rangle = E_n|n\rangle. Then:

a^∣n⟩=cn∣n−1⟩,a^†∣n⟩=cn+1∣n+1⟩\hat{a}|n\rangle = c_n|n-1\rangle, \quad \hat{a}^\dagger|n\rangle = c_{n+1}|n+1\rangle

The constants follow from normalisation. Since a^†a^∣n⟩=n∣n⟩\hat{a}^\dagger\hat{a}|n\rangle = n|n\rangle:

∥cn∣n−1⟩∥2=⟨n∣a^†a^∣n⟩=n\|c_n|n-1\rangle\|^2 = \langle n|\hat{a}^\dagger\hat{a}|n\rangle = n

Therefore:

a^∣n⟩=n ∣n−1⟩,a^†∣n⟩=n+1 ∣n+1⟩\hat{a}|n\rangle = \sqrt{n}\,|n-1\rangle, \quad \hat{a}^\dagger|n\rangle = \sqrt{n+1}\,|n+1\rangle

Ground state. The lowering process must terminate: a^∣0⟩=0\hat{a}|0\rangle = 0. This gives the Differential equation:

(x+ℏmωddx)ϕ0(x)=0  ⟹  ϕ0(x)=(mωπℏ)1/4exp⁡ ⁣(−mωx22ℏ)\left(x + \frac{\hbar}{m\omega}\frac{d}{dx}\right)\phi_0(x) = 0 \implies \phi_0(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\exp\!\left(-\frac{m\omega x^2}{2\hbar}\right)

Energy spectrum. En=ℏω(n+1/2)E_n = \hbar\omega(n + 1/2) for n=0,1,2,…n = 0, 1, 2, \ldots The zero-point energy E0=ℏω/2>0E_0 = \hbar\omega/2 \gt 0 is a direct consequence of [x^,p^]=iℏ[\hat{x}, \hat{p}] = i\hbar.

The eigenfunctions involve Hermite polynomials HnH_n:

ϕn(x)=(mωπℏ)1/412nn!Hn ⁣(mωℏ x)e−mωx2/(2ℏ)\phi_n(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4} \frac{1}{\sqrt{2^n n!}} H_n\!\left(\sqrt{\frac{m\omega}{\hbar}}\,x\right) e^{-m\omega x^2/(2\hbar)}

The first few Hermite polynomials are H0(ξ)=1H_0(\xi) = 1, H1(ξ)=2ξH_1(\xi) = 2\xi, H2(ξ)=4ξ2−2H_2(\xi) = 4\xi^2 - 2.

Example 5.1. Using the ladder operators, find ϕ1(x)\phi_1(x) from ϕ0(x)\phi_0(x).

Solutionϕ1(x)∝a^†ϕ0(x)=mω2ℏ(x−ℏmωddx)ϕ0(x)\phi_1(x) \propto \hat{a}^\dagger\phi_0(x) = \sqrt{\frac{m\omega}{2\hbar}}\left(x - \frac{\hbar}{m\omega}\frac{d}{dx}\right)\phi_0(x) =mω2ℏ ⁣(x+ℏmω⋅mωxℏ)ϕ0(x)=mω2ℏ⋅2x ϕ0(x)= \sqrt{\frac{m\omega}{2\hbar}}\!\left(x + \frac{\hbar}{m\omega}\cdot\frac{m\omega x}{\hbar}\right)\phi_0(x) = \sqrt{\frac{m\omega}{2\hbar}}\cdot 2x\,\phi_0(x)

Normalising gives ϕ1(x)=(mωπℏ)1/42mωℏ x e−mωx2/(2ℏ)\phi_1(x) = \left(\frac{m\omega}{\pi\hbar}\right)^{1/4}\sqrt{\frac{2m\omega}{\hbar}}\,x\,e^{-m\omega x^2/(2\hbar)}.

V(x)=0V(x) = 0 everywhere. The Schrodinger equation:

−ℏ22md2ϕdx2=Eϕ-\frac{\hbar^2}{2m}\frac{d^2\phi}{dx^2} = E\phi

Solutions: ϕk(x)=12πeikx\phi_k(x) = \frac{1}{\sqrt{2\pi}} e^{ikx} with E=ℏ2k22mE = \frac{\hbar^2 k^2}{2m}.

The energy spectrum is continuous (all E≥0E \geq 0). The eigenfunctions are not normalisable (plane Waves); physical states are wave packets constructed by superposition.

The parity operator Π^\hat{\Pi} reflects the coordinate: Π^ψ(x)=ψ(−x)\hat{\Pi}\psi(x) = \psi(-x).

Properties:

  • Π^2=I^\hat{\Pi}^2 = \hat{I} So eigenvalues are ±1\pm 1.
  • Even functions (ψ(−x)=ψ(x)\psi(-x) = \psi(x)) have parity +1+1.
  • Odd functions (ψ(−x)=−ψ(x)\psi(-x) = -\psi(x)) have parity −1-1.
  • If V(x)=V(−x)V(x) = V(-x) (symmetric potential), then [H^,Π^]=0[\hat{H}, \hat{\Pi}] = 0 So energy eigenstates can be chosen to have definite parity.

Theorem 5.1. For a symmetric potential V(x)=V(−x)V(x) = V(-x) The energy eigenstates are either even Or odd.

Proof. Since [H^,Π^]=0[\hat{H}, \hat{\Pi}] = 0 There exists a simultaneous eigenbasis. Let H^ϕ=Eϕ\hat{H}\phi = E\phi and Π^ϕ=πϕ\hat{\Pi}\phi = \pi\phi where π=±1\pi = \pm 1. Then ϕ(−x)=πϕ(x)\phi(-x) = \pi\phi(x) So ϕ\phi is either even (π=+1\pi = +1) or odd (π=−1\pi = -1). ■\blacksquare

This theorem explains why the infinite square well, harmonic oscillator, and finite square well Eigenstates all have definite parity: their potentials are all symmetric about the origin.

Theorem 5.2 (Virial Theorem). For a stationary state of a Hamiltonian H^=p^2/(2m)+V(x^)\hat{H} = \hat{p}^2/(2m) + V(\hat{x}):

2⟨T⟩=⟨x V"(x)⟩2\langle T \rangle = \langle x\,V"(x) \rangle

Where TT is the kinetic energy.

Proof. Using Ehrenfest’s theorem for the operator G^=x^p^\hat{G} = \hat{x}\hat{p}:

ddt⟨x^p^⟩=iℏ⟨[H^,x^p^]⟩=0\frac{d}{dt}\langle\hat{x}\hat{p}\rangle = \frac{i}{\hbar}\langle[\hat{H}, \hat{x}\hat{p}]\rangle = 0

For a stationary state. Computing the commutator:

[H^,x^p^]=[p^22m+V,x^p^]=12m[p^2,x^]p^+[x^p^,V]+x^[V,p^][\hat{H}, \hat{x}\hat{p}] = \left[\frac{\hat{p}^2}{2m} + V, \hat{x}\hat{p}\right] = \frac{1}{2m}[\hat{p}^2, \hat{x}]\hat{p} + [\hat{x}\hat{p}, V] + \hat{x}[V, \hat{p}] =−iℏmp^p^+x^[V,p^]+x^[V,p^]=−iℏp^2m+2iℏx^ V′(x)= \frac{-i\hbar}{m}\hat{p}\hat{p} + \hat{x}[V, \hat{p}] + \hat{x}[V, \hat{p}] = \frac{-i\hbar\hat{p}^2}{m} + 2i\hbar\hat{x}\,V'(x)

Setting d⟨x^p^⟩/dt=0d\langle\hat{x}\hat{p}\rangle/dt = 0 and dividing by iℏi\hbar:

−⟨p^2⟩m+2⟨x^ V′(x^)⟩=0-\frac{\langle\hat{p}^2\rangle}{m} + 2\langle\hat{x}\,V'(\hat{x})\rangle = 0 −2⟨T⟩+⟨x V′(x)⟩=0  ⟹  2⟨T⟩=⟨x V′(x)⟩■-2\langle T \rangle + \langle x\,V'(x) \rangle = 0 \implies 2\langle T \rangle = \langle x\,V'(x) \rangle \qquad \blacksquare

Applications. For the harmonic oscillator (V∝x2V \propto x^2): 2⟨T⟩=2⟨V⟩2\langle T \rangle = 2\langle V \rangle So ⟨T⟩=⟨V⟩=E/2\langle T \rangle = \langle V \rangle = E/2. For the hydrogen atom (V∝−1/rV \propto -1/r): 2⟨T⟩=−⟨V⟩2\langle T \rangle = -\langle V \rangle So ⟨T⟩=−E\langle T \rangle = -E and ⟨V⟩=2E\langle V \rangle = 2E.

Consider V(x)=−V0V(x) = -V_0 for ∣x∣<a|x| \lt a and V(x)=0V(x) = 0 for ∣x∣>a|x| \gt a Where V0>0V_0 \gt 0.

Define k=2m(E+V0)/ℏk = \sqrt{2m(E + V_0)}/\hbar (inside) and κ=−2mE/ℏ\kappa = \sqrt{-2mE}/\hbar (outside). Note that k2+κ2=2mV0/ℏ2k^2 + \kappa^2 = 2mV_0/\hbar^2.

Even parity solutions. Inside: ϕ(x)=Acos⁡(kx)\phi(x) = A\cos(kx). Outside: ϕ(x)=Be−κ∣x∣\phi(x) = Be^{-\kappa|x|}.

Matching ϕ\phi and ϕ′\phi' at x=ax = a and dividing the two conditions:

ktan⁡(ka)=κk\tan(ka) = \kappa

Odd parity solutions. Inside: ϕ(x)=Asin⁡(kx)\phi(x) = A\sin(kx). Outside: ϕ(x)=Be−κ∣x∣\phi(x) = Be^{-\kappa|x|} (with sign For x<0x \lt 0). Matching gives:

−kcot⁡(ka)=κ-k\cot(ka) = \kappa

These are transcendental equations solved graphically. Define z=kaz = ka and z0=a2mV0/ℏ2z_0 = a\sqrt{2mV_0/\hbar^2}.

The even condition becomes tan⁡z=z02/z2−1\tan z = \sqrt{z_0^2/z^2 - 1} and the odd condition becomes −cot⁡z=z02/z2−1-\cot z = \sqrt{z_0^2/z^2 - 1}. The number of bound states is N=⌊2z0/π⌋+1N = \lfloor 2z_0/\pi \rfloor + 1. There is always at least one bound state (the even ground state).

For E>0E \gt 0 The particle has enough energy to escape. Define k1=2mE/ℏk_1 = \sqrt{2mE}/\hbar (outside) And k2=2m(E+V0)/ℏk_2 = \sqrt{2m(E + V_0)}/\hbar (inside). The solutions are oscillatory everywhere. The Transmission coefficient is:

T=11+V024E(E+V0)sin⁡2(2k2a)T = \frac{1}{1 + \dfrac{V_0^2}{4E(E + V_0)}\sin^2(2k_2 a)}

Resonances occur when 2k2a=nπ2k_2 a = n\pi (integer multiples of π\pi), giving T=1T = 1: the well Becomes perfectly transparent.

Example 5.3. A finite square well has V0=5 eVV_0 = 5\,\mathrm{eV} and 2a=1 nm2a = 1\,\mathrm{nm}. Estimate the Number of bound states for an electron.

Solution

Compute z0=a2meV0/ℏz_0 = a\sqrt{2m_e V_0}/\hbar:

z0=(0.5×10−9)2(9.109×10−31)(5)(1.602×10−19)1.055×10−34z_0 = (0.5 \times 10^{-9})\frac{\sqrt{2(9.109 \times 10^{-31})(5)(1.602 \times 10^{-19})}}{1.055 \times 10^{-34}} =(5×10−10)1.460×10−481.055×10−34=(5×10−10)3.821×10−241.055×10−34= (5 \times 10^{-10})\frac{\sqrt{1.460 \times 10^{-48}}}{1.055 \times 10^{-34}} = (5 \times 10^{-10})\frac{3.821 \times 10^{-24}}{1.055 \times 10^{-34}} =(5×10−10)(3.622×1010)=18.11= (5 \times 10^{-10})(3.622 \times 10^{10}) = 18.11

The number of bound states is N=⌊2z0/π⌋+1=⌊36.22/π⌋+1=⌊11.53⌋+1=12N = \lfloor 2z_0/\pi \rfloor + 1 = \lfloor 36.22/\pi \rfloor + 1 = \lfloor 11.53 \rfloor + 1 = 12.

(Actually, the formula is N=⌊z0/(π/2)⌋+1N = \lfloor z_0/(\pi/2) \rfloor + 1 only when counting the number of Intersections. With z0/(π/2)=18.11/1.571=11.53z_0/(\pi/2) = 18.11/1.571 = 11.53 There are 11 full intersections plus one Partial, giving about 11 or 12 bound states.)

Consider V(x)=−αδ(x)V(x) = -\alpha\delta(x) where α>0\alpha \gt 0.

The wave function is ψ(x)=Aeκx\psi(x) = Ae^{\kappa x} for x<0x \lt 0 and ψ(x)=Be−κx\psi(x) = Be^{-\kappa x} for x>0x \gt 0 Where κ=−2mE/ℏ\kappa = \sqrt{-2mE}/\hbar.

Matching conditions.

  1. Continuity: A=BA = B at x=0x = 0.

  2. Discontinuity in derivative (integrating the Schrodinger equation across x=0x = 0):

ψ′(0+)−ψ′(0−)=−2mαℏ2ψ(0)\psi'(0^+) - \psi'(0^-) = -\frac{2m\alpha}{\hbar^2}\psi(0)

This gives −κB−κA=−2mαA/ℏ2-\kappa B - \kappa A = -2m\alpha A/\hbar^2 And since A=BA = B:

κ=mαℏ2\kappa = \frac{m\alpha}{\hbar^2}

The bound state energy is:

E=−ℏ2κ22m=−mα22ℏ2E = -\frac{\hbar^2\kappa^2}{2m} = -\frac{m\alpha^2}{2\hbar^2}

The normalised wave function is ψ(x)=κ e−κ∣x∣\psi(x) = \sqrt{\kappa}\,e^{-\kappa|x|}. There is exactly one bound state.

For a particle of energy E=ℏ2k2/(2m)E = \hbar^2 k^2/(2m) incident from the left:

ψ(x)={eikx+Re−ikxx<0Teikxx>0\psi(x) = \begin{cases} e^{ikx} + Re^{-ikx} & x \lt 0 \\ Te^{ikx} & x \gt 0 \end{cases}

Applying the matching conditions at x=0x = 0:

1+R=T,ik(T−1−R)=−2mαℏ2T1 + R = T, \quad ik(T - 1 - R) = -\frac{2m\alpha}{\hbar^2}T

Solving:

T=ikik−mα/ℏ2=11+imα/(ℏ2k)T = \frac{ik}{ik - m\alpha/\hbar^2} = \frac{1}{1 + im\alpha/(\hbar^2 k)} R=−mα/ℏ2ik−mα/ℏ2=−imα/ℏ2ik+mα/ℏ2R = \frac{-m\alpha/\hbar^2}{ik - m\alpha/\hbar^2} = \frac{-im\alpha/\hbar^2}{ik + m\alpha/\hbar^2}

The transmission and reflection coefficients:

∣T∣2=11+(mα)2/(ℏ4k2)=11+mα2/(2ℏ2E),∣R∣2=1−∣T∣2|T|^2 = \frac{1}{1 + (m\alpha)^2/(\hbar^4 k^2)} = \frac{1}{1 + m\alpha^2/(2\hbar^2 E)}, \quad |R|^2 = 1 - |T|^2

Note that even for very high energies (E→∞E \to \infty), ∣R∣2→(mα)2/(ℏ4k2)≠0|R|^2 \to (m\alpha)^2/(\hbar^4 k^2) \neq 0: The delta function always reflects some probability, unlike a smooth potential which becomes Transparent at high energies. This is because the delta function has an infinitely sharp feature At x=0x = 0 that scatters waves of all wavelengths.

Consider a rectangular barrier V(x)=V0V(x) = V_0 for 0<x<a0 \lt x \lt a and V(x)=0V(x) = 0 otherwise, with E<V0E \lt V_0.

Inside the barrier, the Schrodinger equation gives exponentially decaying and growing solutions:

ψ(x)=Ceκx+De−κx,κ=2m(V0−E)ℏ2\psi(x) = Ce^{\kappa x} + De^{-\kappa x}, \quad \kappa = \sqrt{\frac{2m(V_0 - E)}{\hbar^2}}

For a thick barrier (κa≫1\kappa a \gg 1), the growing solution CeκxCe^{\kappa x} is negligible at the Far edge, and the transmission coefficient simplifies to:

T≈16E(V0−E)V02 e−2κaT \approx \frac{16E(V_0 - E)}{V_0^2}\,e^{-2\kappa a}

The exponential factor e−2κae^{-2\kappa a} is the hallmark of quantum tunnelling: the probability of Penetration decreases exponentially with barrier width and height.

Example 5.2. An electron with E=5E = 5 eV approaches a barrier of height V0=10V_0 = 10 eV and Width a=0.5a = 0.5 nm. Calculate TT.

Solutionκ=2(9.109×10−31)(10−5)(1.602×10−19)(1.055×10−34)2=1.302×1020=1.141×1010  m−1{\kappa = \sqrt{\frac{2(9.109 \times 10^{-31})(10 - 5)(1.602 \times 10^{-19})}{(1.055 \times 10^{-34})^2}} = \sqrt{1.302 \times 10^{20}} = 1.141 \times 10^{10}\;\mathrm{m}^{-1}} 2κa=2(1.141×1010)(5×10−10)=11.412\kappa a = 2(1.141 \times 10^{10})(5 \times 10^{-10}) = 11.41 T≈16(5)(5)100 e−11.41=4.0×e−11.41=4.0×1.097×10−5=4.4×10−5T \approx \frac{16(5)(5)}{100}\,e^{-11.41} = 4.0 \times e^{-11.41} = 4.0 \times 1.097 \times 10^{-5} = 4.4 \times 10^{-5}

The electron has roughly a 0.004%0.004\% chance of tunnelling through this barrier.

Application: alpha decay. Alpha decay can be understood as quantum tunnelling through the Coulomb Barrier. The Geiger-Nuttall law, which relates the decay constant to the alpha particle energy, Follows directly from the exponential dependence of TT on the barrier width.

Application: scanning tunnelling microscope (STM). In an STM, a small voltage is applied between A sharp tip and a conducting surface. Electrons tunnel across the gap, producing a current that Depends exponentially on the tip-surface distance: I∝e−2κdI \propto e^{-2\kappa d}. This allows atomic- Resolution imaging of surfaces, as a change in distance of 0.10.1 nm changes the current by a factor of about 10.

One-dimensional quantum problems are the laboratory where quantum weirdness becomes visible. The infinite square well is like a ball bouncing in a perfectly elastic box: the standing waves that fit inside determine the allowed energies, and the zero-point energy means the ball can never be perfectly still. The harmonic oscillator is the quantum version of a pendulum: the ladder operators let you climb up and down the energy ladder one rung at a time, and each rung costs exactly one quantum of energy. The delta function potential is an infinitely sharp spike that still manages to bind a particle, showing that even an infinitely narrow potential can trap a quantum state. Quantum tunneling is the most dramatic departure from classical physics: a particle can pass through a barrier it classically cannot climb over, like a ball rolling through a wall. The thinner and lower the barrier, the more likely the tunnel, which is how nuclear decay and scanning tunneling microscopes work.

ClaimTrue?Counterexample
A particle in a box can have zero energynoE1=π2ℏ2/2mL2>0E_1 = \pi^2\hbar^2/2mL^2 > 0; zero energy would need zero momentum and a position known to width LL, violating the uncertainty principle
The infinite well is the limit of the finite wellyes, and it is a useful checkas V0→∞V_0 \to \infty the finite-well energies approach the infinite-well ones from below
Tunnelling probability is e−2κae^{-2\kappa a}only for a thick barrierthe exact TT has a prefactor 16E(V0−E)/V0216E(V_0-E)/V_0^2, which matters when EE is close to V0V_0
Parity is always a good quantum numbernoit requires the potential to be even about the origin; shift the well to [0,L][0,L] and parity is lost
A delta potential gives a continuous spectrumnothe attractive delta has exactly one bound state, and it is the only 1D potential with a single one
The harmonic oscillator states are evenly spaced by ℏω\hbar\omegayesthis is special and is why it is exactly soluble; anharmonic terms split the spacing
Higher-energy states have more nodesyesand this is not a coincidence: n−1n-1 nodes is what makes the states orthogonal
Wave functions must be continuousnoψ\psi is continuous only where VV is finite; at an infinite wall it is zero, and at a delta potential ψ\psi is continuous but ψ′\psi' is not

The first row is the one that generalises furthest. The zero-point energy is not a curiosity of the square well — it is why atoms do not collapse, why helium does not freeze at absolute zero, and why the cosmological constant is not obviously zero.

The last row is the one people trip on when moving from the infinite to the finite well: the matching conditions are different, and using the wrong ones gives eigenvalues that are close enough to look plausible.

  • Tunnelling in real devices. The e−2κae^{-2\kappa a} law is the whole of scanning tunnelling microscopy, flash memory and resonant tunnelling diodes; the open problems are inelastic tunnelling and interactions between the tunnelling particle and the barrier.
  • Quantum well engineering. Superlattices and quantum wells are engineered by choosing layer widths and compositions, so the level structure becomes a design variable rather than a given.
  • The WKB approximation beyond 1D. WKB gives the exponential factor and the pre-factor, and its failure near turning points is exactly what the connection formulae patch; how well it carries into higher dimensions is still studied.
  • Stochastic simulation. Propagating the Schrodinger equation numerically for a particle in a well is the standard test problem for time integrators, because the exact propagator is unitary and finite steps are not.
  • Many-particle wells. Adding a second particle to a well makes the problem interacting, and the level structure changes qualitatively; that is where the field’s hard problems begin.

Mix these with the postulates, wave functions, operators, Linear Algebra and Fourier methods.

  1. (Operators.) Compute ⟨x⟩\langle x\rangle and ⟨p⟩\langle p\rangle for the nn-th state of the infinite square well without integrating, using parity: ϕn(x)\phi_n(x) is symmetric or antisymmetric about L/2L/2, so the integrands are odd.
  2. (Linear algebra.) Show that the infinite-well states 2/Lsin⁡(nπx/L)\sqrt{2/L}\sin(n\pi x/L) are orthonormal on [0,L][0,L]. Where does the orthogonality come from — the equation, or the boundary conditions?
  3. (Wave functions.) Show that the virial theorem 2⟨T⟩=n⟨V⟩2\langle T\rangle = n\langle V\rangle holds for the harmonic oscillator and that ⟨T⟩=⟨V⟩\langle T \rangle = \langle V \rangle for the infinite well. Why do the two potentials give different results?
  4. (Classical mechanics.) Compare the classical period of a particle in the infinite well with the quantum beating period between nn and n+1n+1. Show that they agree as n→∞n \to \infty, which is the correspondence principle working rather than being assumed.
  5. (Operators.) Show that the tunnelling exponent κ\kappa is the imaginary part of the wave number inside the barrier, and hence that e−2κae^{-2\kappa a} is ∣eika∣2|e^{ik a}|^2 with kk analytically continued. That continuation is why the same equation gives both.
  6. (Research.) Derive the WKB connection formulae at a single turning point, and state the condition under which they are valid. Then explain why a double turning point needs an extra phase of π/4\pi/4.

Mistake 1: Assuming energy quantisation always occurs Energy is quantised only for bound states (particles trapped in a potential well). For scattering states (E>0E > 0 for a finite potential), the energy spectrum is continuous. Students often assume all quantum systems have discrete energy levels, which is only true for confined particles.

Mistake 2: Forgetting the zero-point energy The ground state energy of a quantum system is never zero for a confining potential. For the infinite square well, E1=π2ℏ2/(2mL2)>0E_1 = \pi^2\hbar^2/(2mL^2) > 0. This is a consequence of the uncertainty principle: confining a particle to a region of size LL requires momentum ∼ℏ/L\sim \hbar/L, giving kinetic energy ∼ℏ2/(2mL2)\sim \hbar^2/(2mL^2). Students sometimes set the lowest energy to zero as in classical mechanics.

Mistake 3: Misunderstanding quantum tunnelling Tunnelling does not require the particle to “have” energy greater than the barrier height inside the barrier. The wave function extends into the classically forbidden region with exponentially decreasing amplitude. The particle’s energy is E<V0E < V_0 throughout, and the tunnelling probability decreases exponentially with barrier width and height.