Schrödinger Equation

The Schrödinger equation is the linear partial differential equation that governs the spatial distribution and temporal evolution of the quantum wavefunction of a non-relativistic physical system. Formulated by the Austrian physicist Erwin Schrödinger in January 1926, the equation plays a role in quantum mechanics analogous to Newton's second law of motion or Hamilton's equations in classical mechanics. Rather than predicting a deterministic classical trajectory with simultaneously defined position and linear momentum, the Schrödinger equation acts on a complex-valued state function, denoted Ψ(r, t), from which all measurable physical observables—including discrete energy levels, angular momenta, and spatial electron densities—are extracted through linear Hermitian operators.

In atomic physics and quantum chemistry, solving the time-independent form of the Schrödinger equation for an electron moving within an electrostatic nuclear potential yields the quantized spatial wavefunctions known as atomic orbitals and naturally generates the principal, azimuthal, and magnetic quantum numbers. Combined with the relativistic spin requirement and the Pauli exclusion principle, these solutions establish the ground-state electron configuration of the elements and provide the mathematical foundation for molecular orbital theory and chemical bonding.

The Time-Dependent Wave Equation and Hamiltonian Operator

For a single particle of mass m moving through a potential energy field V(r, t), the general time-dependent Schrödinger equation (TDSE) takes the form:

*i*ℏ (∂Ψ(r, t) / ∂*t*) = Ĥ Ψ(r, t)

Here, i is the imaginary unit (i² = −1), ℏ = h / (2π) ≈ 1.05457 × 10⁻³⁴ J·s is the reduced Planck constant, and Ĥ is the Hamiltonian operator representing the total energy of the system. In three-dimensional Cartesian coordinates, the non-relativistic Hamiltonian is the sum of the translational kinetic energy operator T̂ (derived from particle velocity and momentum p = −*i*ℏ∇) and the potential energy operator V̂:

Ĥ = −(ℏ² / 2*m*) ∇² + V(r, t)

where ∇² = ∂²/∂*x*² + ∂²/∂*y*² + ∂²/∂*z*² is the Laplacian differential operator. Because the time-dependent equation is first-order in time and contains the factor i, its general solutions exhibit oscillatory complex phase behavior proportional to exp(−*iEt*/ℏ), ensuring that total probability is conserved over time via unitary time evolution.

Separation of Variables and the Eigenvalue Problem

When the potential energy V(r) is independent of time—as in an isolated atom or stationary molecule in the absence of external electromagnetic fields—the total wavefunction can be factored into spatial and temporal components via separation of variables:

Ψ(r, t) = ψ(r) · exp(−*iEt* / ℏ)

Substituting this product into the time-dependent equation and dividing out the exponential time phase exp(−*iEt* / ℏ) yields the time-independent Schrödinger equation (TISE):

Ĥ ψ(r) = E ψ(r)

Mathematically, the TISE is a linear Sturm–Liouville eigenvalue equation in which the Hamiltonian operator Ĥ acts on a stationary spatial eigenfunction ψ(r) to return that same function multiplied by a real scalar eigenvalue E, which corresponds to the quantized total energy of the stationary state.

Born's Probabilistic Interpretation and Normalization

In 1926, Max Born proposed the statistical interpretation of the wavefunction: while ψ(r) itself is complex-valued and not directly measurable, the real, non-negative quantity |Ψ(r, t)|² = Ψ*(r, t)Ψ(r, t) represents the probability density per unit volume of detecting the particle at position r at time t. For a wavefunction to be physically acceptable, it must be continuous, single-valued, possess a continuous first spatial derivative (except at pathological potential singularities), and be square-integrable over all space so that the normalization integral equals unity:

∫_−∞^+∞ Ψ*(r, t) Ψ(r, t) *d*τ = 1

The average measurable value of any physical observable A—represented by a Hermitian operator —is calculated as the expectation value ⟨*A*⟩ = ∫ ψ*  ψ *d*τ.

Only a small class of idealized low-dimensional or highly symmetric potential energy functions V(r) admit closed-form analytical solutions to the time-independent Schrödinger equation. Nevertheless, these canonical models form the conceptual and mathematical basis for spectroscopy and atomic structure.

Particle in a Box and Harmonic Oscillator

Confining a particle of mass m to a one-dimensional potential well of length L bounded by impenetrable walls (V = ∞ at x ≤ 0 and x ≥ L) forces the wavefunction ψ(x) to vanish at both boundaries. This boundary condition restricts the de Broglie wavelength to λ_n = 2*L*/n, quantizing the translational kinetic energy of the particle moving along the box axis into discrete levels E_n = n²h²/(8*mL*²) and establishing a non-zero ground-state zero-point energy (E₁ > 0) required by the Heisenberg uncertainty principle. Similarly, solving the Schrödinger equation for a parabolic restoring potential V(x) = (1/2)*k*_fx² yields evenly spaced vibrational energy levels E_v = ℏω(v + 1/2) with zero-point vibrational energy (1/2)ℏω, where Gaussian-weighted Hermite polynomials penetrate into the classically forbidden region where E_v < V(x) via quantum tunneling.

The Hydrogen Atom and Spherical Harmonics

For a one-electron hydrogenic ion of nuclear charge +*Ze* and reduced electron–nucleus mass μ = m_em_N/(m_e + m_N) interacting via the spherically symmetric Coulomb potential V(r) = −*Ze*²/(4πε₀r), transforming the Laplacian into spherical polar coordinates (r, θ, φ) separates the spatial wavefunction into a radial component and an angular component:

ψ_n,*l*,*m*(r, θ, φ) = R_n,*l*(r) · Y_l,*m*(θ, φ)

Enforcing single-valuedness and square-integrability on the three separated differential equations naturally introduces three integer quantum numbers: the principal quantum number n = 1, 2, 3... (governing the radial Laguerre polynomials and energy E_n = −13.6057 eV · Z²/n²), the azimuthal quantum number l = 0, 1, ..., n − 1 (quantizing orbital angular momentum magnitude L = ℏ[l(l + 1)]^1/2), and the magnetic quantum number m_l = −*l*, ..., +*l* (quantizing the z-projection L_z = m_lℏ).

The Three-Body Problem and Born–Oppenheimer Approximation

In any atom or molecule containing two or more electrons (such as helium, H₂, or H₂O), the interelectronic Coulomb repulsion terms e²/(4πε₀r_ij)—which depend inversely on the instantaneous distance r_ij = |*r*_i − r_j| between electrons i and j—couple the spatial coordinates and prevent exact analytical separation of variables. To make molecular systems tractable, Max Born and J. Robert Oppenheimer introduced the Born–Oppenheimer approximation in 1927. Because atomic nuclei are 1,836 to >10⁵ times more massive than electrons, nuclear velocities and translational kinetic energies are negligible on the attosecond timescale of electronic motion. Consequently, the nuclear kinetic energy operator can be neglected when solving the electronic Schrödinger equation at fixed internuclear distances R:

Ĥ_elec ψ_elec(r; R) = E_elec(R) ψ_elec(r; R)

Repeating this calculation across a grid of nuclear geometries generates the adiabatic potential energy surface (PES) on which nuclear vibrations and chemical reactions occur.

Variational Principle and Hartree–Fock Theory

To approximate the coupled many-electron wavefunction ψ_elec, quantum chemistry relies on the variational theorem, which states that the expectation value of the Hamiltonian calculated with any normalized trial wavefunction Φ_trial is an upper bound to the true ground-state energy E₀:

⟨Φ_trial | Ĥ | Φ_trial⟩ ≥ E₀

In the Hartree–Fock self-consistent field (SCF) method, the N-electron wavefunction is approximated as a single antisymmetrized Slater determinant of one-electron spin-orbitals, satisfying the Pauli exclusion principle. Each electron is treated as moving within the static average electrostatic and exchange field generated by the nuclei and the remaining N − 1 electrons, reducing the coupled N-body problem to N coupled one-electron integro-differential equations solved iteratively until convergence. Post-Hartree–Fock methods (such as Møller–Plesset perturbation theory and Coupled Cluster) and Density Functional Theory (DFT) incorporate the remaining instantaneous dynamic electron correlation required for quantitative thermochemical accuracy.

Why does the Schrödinger equation fail for relativistic electrons in heavy elements?

The standard Schrödinger equation is non-relativistic because it combines a second-order spatial derivative (derived from classical translational kinetic energy p²/2*m*) with a first-order time derivative, violating Lorentz invariance. In heavy elements (such as gold, Z = 79, or mercury, Z = 80), inner 1*s* electrons reach velocities exceeding 50% of the speed of light c. Accurate modeling requires Paul Dirac's 1928 relativistic wave equation, which naturally predicts electron spin (s = 1/2), spin-orbit coupling, and the relativistic contraction of s orbitals.

What is the physical meaning of nodes in a Schrödinger wavefunction?

A node is a point, plane, or spherical surface where the spatial wavefunction ψ(r) passes through zero and changes algebraic phase sign (from + to −), making the probability density |ψ(r)|² zero. In any hydrogenic atomic orbital, the total number of nodes equals n − 1, partitioned into l angular (planar or conical) nodes and n − l − 1 radial (spherical) nodes. Higher node counts correspond to steeper wavefunction curvature (∇²ψ) and higher electronic kinetic energy.

How does the Schrödinger equation explain quantum tunneling in chemical reactions?

When a particle of energy E encounters a finite potential energy barrier V₀ > E of thickness a, classical mechanics forbids penetration. However, solving the Schrödinger equation inside the barrier (E − V₀ < 0) yields a real decaying exponential solution ψ(x) ∝ exp(−κ*x*), where κ = [2*m*(V₀ − E)]^1/2 / ℏ. If the barrier is narrow and particle mass m is small (such as proton or hydride transfer in enzymatic catalysis), ψ(x) retains non-zero amplitude on the product side.

1. Schrödinger, E. (1926). "An Undulatory Theory of the Mechanics of Atoms and Molecules." Physical Review, 28(6), 1049–1070. https://doi.org/10.1103/PhysRev.28.1049
2. Born, M. (1926). "Zur Quantenmechanik der Stoßvorgänge." Zeitschrift für Physik, 37(12), 863–867. https://doi.org/10.1007/BF01397477
3. Born, M., & Oppenheimer, R. (1927). "Zur Quantentheorie der Molekeln." Annalen der Physik, 389(20), 457–484. https://doi.org/10.1002/andp.19273892002
4. Szabo, A., & Ostlund, N. S. (1996). Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. Dover Publications. https://doi.org/10.1021/ed067pA280.2