In 1926, Erwin Schrödinger constructed arrays from differential operators and complete orthogonal function systems, showing how their products reproduce the array calculus. His paper connected wave mechanics with matrix mechanics; the translation's sections 2–3 develop the product and commutation rules.[ref]
The modern proof is a change of orthonormal basis. If form a complete orthonormal basis of the physical Hilbert space, map a wavefunction to its coefficients:
Completeness makes this map unitary: it preserves the inner product and reconstructs from its coefficients. Inserting the complete orthonormal basis proves the matrix multiplication rule. The Schrödinger equation becomes
For a time-independent Hamiltonian operator, the Heisenberg picture transfers the same evolution to the observables, preserving measurement predictions. Unitary equivalence requires matching operator domains as well as formal entries. A finite selection of eigenfunctions gives a compression, not an exact replacement of the full Hilbert space; a continuous spectrum needs its spectral integrals. This equivalence concerns representations of specified operators, and does not remove the ordering question when quantizing an arbitrary classical expression, as Weyl ordering illustrates.
Matrix mechanics 2026-10-18
Matrix mechanics represents an observable by an array of matrix elements between chosen quantum states. In an orthonormal basis of energy eigenstates, each a stationary state, the diagonal entries of the Hamiltonian operator are the energy levels; off-diagonal entries of other observables describe couplings between those states. The same array can therefore answer two different questions: which values can a measurement return, and which pairs of energy eigenstates can a specified interaction connect? The eigenvalues answer the first question; the relevant matrix elements, together with the energy differences, answer the second.
This makes linear algebra a practical language for atomic physics. One can calculate energy levels, electric-dipole selection rules and probability amplitudes for transitions without assigning an Electron an unobserved classical orbit. The three-state example below separates a missing optical coupling from a missing energy level; the finite oscillator commutator boundary defect explains why a finite array requires care when approximating a quantum harmonic oscillator.
Consider an original ideal quantum system with positive scales and
Interpret as one component of an electric dipole moment, so its coupling to the corresponding electric field is the electric-dipole interaction . The gaps are , but : this interaction has no direct first-order coupling. The nonzero squared matrix elements are and . These are coupling weights; actual spectral line intensities also depend on populations and radiation factors. In particular, there is no promise that an experiment populates the upper state of every allowed transition.
Figure 1.
Matrix mechanics: three energy levels, an observable array and two nonzero coupling weights
. Original ideal three-state model. The gap exists even though its electric dipole matrix element vanishes. The plotted weights are squared couplings, not measured emission intensities. Details.
A measurement of itself has possible outcomes , since
These are eigenvalues of , not the energy eigenvalues . Moreover, although : two factors connect through . This algebraic statement does not turn into a nonzero direct coupling. Higher-order probability amplitudes for transitions require the appropriate dynamical calculation.