Quantum Mechanics & Wave Functions — Schrödinger, Tunnelling, Double-Slit and Quantum Circuits

Why Quantum Mechanics Needs Simulation

Classical mechanics lets you picture a ball rolling down a hill. Quantum mechanics describes an electron as a complex-valued wave function ψ(x,t) whose squared magnitude |ψ|² gives a probability density. There is no intuitive mental image for this. The wave function doesn't live in ordinary space — it lives in a Hilbert space. When you "measure" the electron, the wave function collapses to a definite value in a way that no classical process can reproduce.

Simulation helps by rendering the otherwise invisible. An interactive Schrödinger wave packet lets you watch superposition and interference directly. Quantum tunnelling becomes visually obvious when you see the probability amplitude leak through a barrier that the particle classically couldn't cross. Grover's algorithm, which looks like pure abstract algebra on paper, snaps into focus when you watch amplitudes amplify iteratively toward the marked state.

This guide follows five progressive layers of quantum theory, each anchored to a simulation you can open and manipulate right now.

Part 1 — The Wave Function and Schrödinger's Equation

The Time-Dependent Schrödinger Equation

The central equation of quantum mechanics is the time-dependent Schrödinger equation (TDSE). It governs how a quantum state ψ(x,t) evolves over time in the presence of a potential V(x). Unlike Newton's F=ma — a differential equation for a trajectory — the TDSE is a differential equation for a probability amplitude. The wave function ψ is complex-valued; its modulus squared |ψ(x,t)|² gives the probability of finding the particle at position x at time t.

The Schrödinger simulation solves the 1D TDSE numerically using the Crank-Nicolson method, which is unconditionally stable and norm-conserving (it keeps ∫|ψ|² = 1 at every step). You can place a Gaussian wave packet — a localised probability bump with a mean momentum ℏk₀ — and watch it propagate, spread due to dispersion, and interact with barriers.

Part 2 — Quantum Tunnelling

Classically Forbidden Regions

Classical mechanics is absolute: a particle at position x with total energy E cannot enter a region where the potential energy V(x) > E. The kinetic energy K = E − V would be negative, which is impossible for a classical particle. Quantum mechanics has no such prohibition. The wave function can be non-zero inside a classically forbidden region — it decays exponentially, but doesn't vanish. If the barrier is thin enough, the wave function emerges on the other side with finite amplitude. This is quantum tunnelling.

Tunnelling is not a rare exotic effect. It is the mechanism behind alpha decay in radioactive nuclei, the operation of tunnel diodes, scanning tunnelling microscopy (STM), and enzymatic hydrogen transfer in biochemistry. The Sun itself depends on quantum tunnelling: without it, proton–proton fusion would be impossible at solar core temperatures (≈15 million K is far too cold to overcome the Coulomb barrier classically).

The scanning tunnelling microscope achieves sub-angstrom spatial resolution because tunnelling current falls by roughly an order of magnitude for every 1 Å increase in tip-to-sample distance. This extraordinary sensitivity — exponential dependence on gap width — turns each atom's topography directly into current fluctuations readable with a lock-in amplifier. The quantum tunnelling simulation captures this physics in one dimension with adjustable barrier width and height.

Part 3 — Atomic Structure: Hydrogen Orbitals

The Quantum Numbers of Hydrogen

The hydrogen atom is the only atom for which the Schrödinger equation has an exact analytic solution. Solving the 3D TDSE in spherical coordinates produces a set of wave functions ψ_nlm(r,θ,φ) labelled by three quantum numbers: n (principal, sets energy), l (orbital angular momentum, 0 ≤ l ≤ n−1), and m (magnetic quantum number, −l ≤ m ≤ l). Each unique (n,l,m) triple is an orbital — a spatial probability distribution for the electron.

The shapes of atomic orbitals determine chemical bonding. The 2p orbital has two lobes separated by a nodal plane — this geometry forces sigma bonds to form along the internuclear axis. The 3d orbitals' complex lobed geometries create the crystal field splitting responsible for the colours of transition metal compounds. Understanding the orbital shapes as probability distributions (rather than "shells") is the conceptual step from Bohr's model to modern quantum chemistry.

Part 4 — Quantum Computing: Qubits and Gates

From Classical Bits to Qubits

A classical bit is either 0 or 1. A qubit is a quantum system that can exist in a superposition of |0⟩ and |1⟩: |ψ⟩ = α|0⟩ + β|1⟩, where |α|² + |β|² = 1. The qubit's state is a unit vector in a 2D complex Hilbert space — visualised as a point on the Bloch sphere. Quantum gates are unitary operations that rotate this vector, and they are the quantum analogue of classical logic gates.

The power of quantum computing comes from entanglement: two entangled qubits cannot be described independently — their joint state has correlations that no classical probability distribution can replicate. This is what Bell's theorem proves: the CHSH inequality S ≤ 2 must hold for any local hidden variable theory, but quantum mechanics predicts — and experiment confirms — S = 2√2 ≈ 2.83. The quantum entanglement simulation runs a live CHSH test that you can verify by adjusting measurement angles.

Part 5 — Grover's Quantum Search Algorithm

Quadratic Speedup for Unstructured Search

Given an unsorted database of N items with one marked item, a classical computer requires O(N) queries on average to find it. Grover's algorithm solves the same problem in O(√N) quantum queries — a quadratic speedup achieved by exploiting quantum amplitude amplification. The algorithm requires only about π√N/4 iterations to reach near-certain success probability.

Each iteration of Grover's algorithm consists of two operations: (1) an oracle that flips the sign (phase) of the marked state's amplitude, and (2) a diffusion operator that reflects all amplitudes about their mean. The net effect is to increase the marked state's amplitude by approximately 2/√N per step. After ~π√N/4 steps, the marked state has near-unit amplitude and measurement finds it with high probability.

Grover's algorithm demonstrates a key principle of quantum advantage: it doesn't require the database to have any structure — it works on any oracle that can identify the marked item. The amplitude amplification technique it introduces is also a subroutine in more powerful quantum algorithms. In the simulation you can vary the database size (N), watch the sinusoidal probability evolution, and see that stopping at exactly the optimal iteration count matters — run too many steps and the probability starts falling again.

Quantum Connections Across the Collection

Quantum mechanics underpins simulations across many other categories on the platform:

• Blackbody Radiation — Planck's quantum hypothesis (E=hν) resolved the ultraviolet catastrophe of classical physics and laid the foundation for quantum mechanics.
• Mass–Energy (E=mc²) — nuclear fission and fusion are quantum processes; the binding energy curve is calculated from the Bethe-Weizsäcker semi-empirical mass formula.
• Nuclear Fusion — the proton–proton chain in the Sun requires quantum tunnelling through the Coulomb barrier.
• Brownian Motion — at nanoscale, thermal fluctuations become comparable to quantum zero-point energy; De Broglie wavelength λ = h/mv starts to matter.
• Crystal Structures — the metallic, ionic, and covalent bonding types in the crystal structures simulation are all quantum mechanical in origin.

Algorithms & Methods