Schrödinger Equation - Physical Chemistry 2

SECTION 1

Historical Context & Motivation

By the mid-1920s, physics had accumulated a wealth of experimental evidence—blackbody radiation, the photoelectric effect, atomic line spectra—that classical mechanics simply could not explain. The old quantum theory of Bohr and Sommerfeld had enjoyed partial success by imposing quantization conditions ad hoc on otherwise classical orbits, yet it failed for multi-electron atoms and could not predict transition intensities. Meanwhile, de Broglie's 1924 hypothesis that matter possesses wave-like properties demanded a proper dynamical equation—an equation that would govern the evolution of these matter waves just as Maxwell's equations govern electromagnetic waves. The stage was set for a unifying framework.

1900 • Planck's Quantum Hypothesis

Max Planck introduces energy quantization (E = nhν) to explain the blackbody radiation spectrum, planting the seed of quantum theory.

1913 • Bohr's Atomic Model

Niels Bohr proposes quantized electron orbits for hydrogen, successfully reproducing the Rydberg formula but lacking a deeper theoretical justification for the quantization rules.

1924 • De Broglie's Matter Waves

Louis de Broglie postulates that every particle has an associated wavelength λ = h/p, unifying the wave-particle duality first observed for photons to all matter.

1926 • Schrödinger's Wave Mechanics

Erwin Schrödinger publishes a series of landmark papers deriving a wave equation for quantum systems, demonstrating its equivalence to Heisenberg's matrix mechanics and solving the hydrogen atom exactly.

1927 • Born's Probabilistic Interpretation

Max Born interprets |Ψ|² as a probability density, completing the physical meaning of the wave function and establishing the statistical foundation of quantum mechanics.

The central question that Schrödinger confronted was deceptively simple: if particles behave as waves, what equation governs the propagation and evolution of those waves? Classical wave equations—such as the equation for vibrating strings or sound in air—are second-order in time, yet quantum mechanics would require something fundamentally different: an equation first-order in time and intimately linked to the energy-momentum relation through Planck's constant. The answer Schrödinger provided in 1926 not only resolved the question but also laid the mathematical cornerstone upon which virtually all of non-relativistic quantum mechanics rests.

SECTION 2

Core Principles & Definitions

Understanding the Schrödinger equation requires familiarity with several foundational concepts that together form the conceptual scaffolding of wave mechanics. The equation does not emerge in isolation; it rests on the postulates of quantum mechanics, the mathematical structure of Hilbert space, and the physical intuition provided by de Broglie's hypothesis. The following core ideas are essential before engaging with the equation itself.

The Wave Function Ψ

A complex-valued function Ψ(x, t) encodes all physically obtainable information about a quantum system. It is not directly observable; rather, its modulus squared |Ψ|² gives the probability density for finding the particle at position x and time t.

Operators & Observables

Every measurable quantity (position, momentum, energy) is represented by a linear, Hermitian operator acting on Ψ. The momentum operator is p̂ = −iℏ(∂/∂x) and the energy operator is Ê = iℏ(∂/∂t).

Superposition Principle

If Ψ₁ and Ψ₂ are valid solutions of the Schrödinger equation, then any linear combination c₁Ψ₁ + c₂Ψ₂ is also a valid solution. This linearity underpins quantum interference and entanglement.

Normalization Condition

Because |Ψ|² is a probability density, the total probability of finding the particle somewhere must equal unity: ∫|Ψ|²dx = 1. This constraint restricts the class of physically acceptable wave functions.

Hamiltonian Operator

The Hamiltonian Ĥ = T̂ + V̂ represents the total energy—kinetic plus potential. It generates time evolution and its eigenvalues are the allowed energies of the system.

✦

KEY TAKEAWAY

Think of the wave function Ψ as a detailed weather map for a quantum particle: just as a meteorological probability map tells you the likelihood of rain at each location without specifying where a particular raindrop will fall, |Ψ|² provides a probability landscape for where the particle is likely to be detected. The Schrödinger equation is the dynamical law that tells you how this probability landscape evolves over time—analogous to the atmospheric equations that evolve the weather forecast.

SECTION 3

Visualizing the Wave Function

The diagram below illustrates the wave function and its corresponding probability density for the first three energy eigenstates of a particle in a one-dimensional infinite square well—one of the most fundamental exactly solvable systems in quantum mechanics. The well confines the particle between x = 0 and x = L with infinitely high walls, forcing the wave function to vanish at the boundaries. Notice how the number of nodes (zero crossings) in Ψ_n equals n − 1, and how the probability density |Ψ_n|² reveals regions where the particle is most and least likely to be found.

Left column: the wave functions Ψ_n(x) for n = 1, 2, 3. The vertical walls represent the infinite potential barriers. Right column: the corresponding probability densities |Ψ_n|². Note how higher quantum numbers produce more nodes and a more uniform distribution—approaching the classical limit.

Several important features are visible in this diagram. First, the wave functions are sinusoidal: Ψ_n(x) = √(2/L) sin(nπx/L), which are the eigenfunctions of the Hamiltonian for this system. Second, the number of nodes directly encodes the quantum number and hence the energy: more nodes mean higher kinetic energy. Third, the probability density |Ψ_n|² is always non-negative and integrates to unity over the box. For the ground state (n = 1), the particle is most likely to be found near the center, whereas a classical particle bouncing back and forth would be equally likely to be found anywhere—a stark departure that highlights the non-classical nature of quantum confinement.

SECTION 4

Mathematical Framework

The Schrödinger equation exists in two primary forms. The time-dependent Schrödinger equation (TDSE) describes how a quantum state evolves in time, while the time-independent Schrödinger equation (TISE) arises upon separation of variables when the potential V does not depend on time. The TISE is an eigenvalue equation whose solutions yield the stationary states—the energy eigenfunctions and their corresponding eigenvalues.

TIME-DEPENDENT SCHRÖDINGER EQUATION

iℏ ∂Ψ(x,t)/∂t = ĤΨ(x,t) = [−(ℏ²/2m)(∂²/∂x²) + V(x,t)] Ψ(x,t)

Here i is the imaginary unit, ℏ = h/(2π) is the reduced Planck constant, m is the particle mass, V(x,t) is the potential energy, and Ĥ is the Hamiltonian operator. The equation is first-order in time, reflecting the deterministic evolution of probability amplitudes.

TIME-INDEPENDENT SCHRÖDINGER EQUATION

Ĥψ(x) = Eψ(x) → −(ℏ²/2m)(d²ψ/dx²) + V(x)ψ(x) = Eψ(x)

When V is independent of time, we write Ψ(x,t) = ψ(x)e^−iEt/ℏ. Substituting into the TDSE and dividing by the time factor yields this eigenvalue equation. E is the total energy eigenvalue, and ψ(x) is the spatial part of the wave function.

PARTICLE-IN-A-BOX ENERGY EIGENVALUES

Eₙ = n²π²ℏ² / (2mL²), n = 1, 2, 3, …

For a particle confined in a box of length L, the boundary conditions Ψ(0) = Ψ(L) = 0 select only discrete values of energy. Note that E₁ > 0: the particle can never have zero kinetic energy, a direct consequence of the Heisenberg uncertainty principle.

ℹ

Derivation Sketch

Starting from the classical energy expression E = p²/(2m) + V, one replaces p → −iℏ(∂/∂x) and E → iℏ(∂/∂t)—the canonical quantization prescription. Operating on Ψ(x,t) immediately yields the TDSE. This substitution is motivated by the de Broglie relations (p = ℏk, E = ℏω) applied to a plane wave Ψ = Ae^{i(kx − ωt)}. The requirement that Ψ satisfy a linear, first-order-in-time differential equation uniquely constrains the form of the Schrödinger equation.

SECTION 5

Solutions for Key Quantum Systems

The power of the Schrödinger equation is best appreciated through its application to model systems of increasing complexity. Each system introduces new physical features—degeneracy, angular momentum, tunneling—that recur throughout quantum chemistry. The diagram below compares the energy level structures of three canonical systems: the infinite well, the harmonic oscillator, and the hydrogen atom.

Energy level spacing comparison. The infinite well has quadratically increasing gaps (E ∝ n²). The harmonic oscillator has equally spaced levels (ΔE = ℏω). The hydrogen atom has levels that converge toward zero (E ∝ −1/n²), reflecting the Coulomb potential.

System	Potential V(x)	Energy Eigenvalues	Key Feature
Infinite Well	V = 0 inside, ∞ outside	E_n = n²π²ℏ²/(2mL²)	Discrete, quadratically growing spectrum; strict confinement
Harmonic Oscillator	V = ½kx² = ½mω²x²	E_n = (n + ½)ℏω	Equally spaced levels; zero-point energy E₀ = ½ℏω
Hydrogen Atom	V = −e²/(4πε₀r)	E_n = −13.6 eV / n²	Degeneracy in ℓ; converging levels; bound + continuum states

Comparison of three fundamental exactly solvable quantum systems

SECTION 6

Worked Example: Electron in a Quantum Dot

Consider an electron confined in a semiconductor quantum dot that can be modeled as a one-dimensional infinite square well of width L = 5.0 nm. We wish to find the ground-state energy, the wavelength of a photon emitted during the n = 2 → n = 1 transition, and the probability of finding the electron in the central third of the well in its ground state.

Electron in a 5.0 nm Quantum Dot

1
Step 1 — Identify Given Values & Relevant Equation

We use E_n = n²π²ℏ²/(2mL²). The known quantities are: m_e = 9.109 × 10⁻³¹ kg, L = 5.0 × 10⁻⁹ m, ℏ = 1.055 × 10⁻³⁴ J·s.
2
Step 2 — Compute Ground-State Energy (n = 1)

E₁ = (1)²π²(1.055 × 10⁻³⁴)² / [2 × (9.109 × 10⁻³¹) × (5.0 × 10⁻⁹)²]. Evaluating the numerator: π² × (1.113 × 10⁻⁶⁸) = 1.098 × 10⁻⁶⁷ J²·s². The denominator: 2 × 9.109 × 10⁻³¹ × 2.5 × 10⁻¹⁷ = 4.555 × 10⁻⁴⁷ J·m²·kg⁻¹ … wait, let's be precise. Denominator = 2mL² = 2 × (9.109 × 10⁻³¹) × (25.0 × 10⁻¹⁸) = 4.555 × 10⁻⁴⁷ kg·m².

E₁ = 1.098 × 10⁻⁶⁷ / (4.555 × 10⁻⁴⁷) = 2.41 × 10⁻²¹ J ≈ 0.0150 eV
3
Step 3 — Compute E₂ and the Transition Energy

Since E_n ∝ n², we have E₂ = 4E₁ = 4 × 0.0150 eV = 0.0601 eV. The photon energy for the 2 → 1 transition is ΔE = E₂ − E₁ = 3E₁ = 0.0451 eV.

ΔE = 0.0451 eV = 7.22 × 10⁻²¹ J
4
Step 4 — Photon Wavelength

Using λ = hc/ΔE: λ = (6.626 × 10⁻³⁴ × 2.998 × 10⁸) / (7.22 × 10⁻²¹) = 1.987 × 10⁻²⁵ / 7.22 × 10⁻²¹.

λ ≈ 27.5 μm (mid-infrared region)
5
Step 5 — Probability in the Central Third (Ground State)

P = ∫ from L/3 to 2L/3 of |Ψ₁|² dx = (2/L) ∫ from L/3 to 2L/3 sin²(πx/L) dx. Using the identity sin²(u) = [1 − cos(2u)]/2 and substituting u = πx/L: P = (2/L) × [x/2 − (L/(4π))sin(2πx/L)] evaluated from L/3 to 2L/3. This yields P = 1/3 + (1/(2π))sin(2π/3) − (−1/(2π))sin(2π/3) = 1/3 + (√3)/(2π).

P = 1/3 + √3/(2π) ≈ 0.609 (about 61%) — significantly larger than the classical 33%, confirming quantum localization near the center.

SECTION 7

Strengths, Limitations, and Scope

The Schrödinger equation is the workhorse of non-relativistic quantum mechanics, but it is important to recognize both its extraordinary reach and the boundaries where it must be supplemented or replaced by more advanced theories. The table below summarizes the principal strengths and limitations.

Strengths	Limitations
Exact solutions for hydrogen, harmonic oscillator, particle in a box, rigid rotor—foundational models in chemistry	Non-relativistic: fails for particles at speeds approaching c (need Dirac equation)
Naturally produces quantized energy levels, reproducing atomic spectra with extraordinary precision	No exact solutions for multi-electron atoms; requires approximate methods (HF, DFT, CI)
Linear equation: superposition principle enables powerful mathematical techniques (perturbation theory, variational method)	Does not naturally incorporate spin; spin must be added as a postulate or derived from Dirac theory
Generalizes seamlessly to 3D, many-body systems, and time-dependent problems	Cannot describe particle creation/annihilation (need quantum field theory)
Basis for all of computational quantum chemistry (molecular orbital theory, spectroscopy predictions)	Measurement problem and wave function collapse are not explained by the equation itself

Strengths and limitations of the non-relativistic Schrödinger equation

✦

KEY TAKEAWAY

The Schrödinger equation is to quantum chemistry what Newton's second law is to classical mechanics—it is the equation of motion for quantum systems. Just as F = ma breaks down at relativistic speeds and must be replaced by special relativity, the Schrödinger equation breaks down at relativistic energies and must be superseded by the Dirac equation. Within its domain of validity (v ≪ c, no particle creation), it remains astonishingly accurate and forms the backbone of modern computational chemistry.

SECTION 8

Connection to Advanced Quantum Theory

The Schrödinger equation is the first rung on a ladder of increasingly sophisticated quantum theories. Understanding where it sits in this hierarchy clarifies why physical chemists continue to use it for molecular problems while high-energy and condensed-matter physicists often invoke its generalizations. The table below contrasts the Schrödinger framework with its two most important extensions: the Dirac equation and quantum field theory (QFT).

Feature	Schrödinger Equation	Dirac Equation	QFT (Second Quantization)
Relativity	Non-relativistic (v ≪ c)	Fully relativistic	Fully relativistic
Spin	Added ad hoc (Pauli matrices)	Emerges naturally as a 4-component spinor	Encoded in field operators
Particle number	Fixed	Fixed (but predicts antiparticles)	Variable (creation/annihilation)
Typical use in chemistry	Almost all molecular calculations	Heavy-element relativistic corrections	Photochemistry, radiation–matter coupling
Mathematical form	PDE, first-order in t, second-order in x	PDE, first-order in t and x (4×4 matrix)	Operator equations on Fock space

Schrödinger equation in the hierarchy of quantum theories

For most problems in physical chemistry—molecular orbital theory, vibrational spectroscopy, reaction dynamics, NMR chemical shifts—the non-relativistic Schrödinger equation provides ample accuracy. Relativistic corrections become important primarily for elements with Z > 50, where inner-shell electrons reach appreciable fractions of c. Even in those cases, the standard approach is to solve the Schrödinger equation first and then apply perturbative relativistic corrections (scalar relativity, spin-orbit coupling) rather than solving the full Dirac equation. This pragmatic strategy underscores the centrality of Schrödinger's framework to the entire discipline.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the time-dependent Schrödinger equation is first-order in time rather than second-order, and discuss how this relates to the fact that the wave function is complex-valued. Why does this differ from the classical wave equation for a vibrating string?

Show answer

The classical wave equation ∂²u/∂t² = v²∂²u/∂x² is second-order in time and requires two initial conditions (u and ∂u/∂t), yielding real-valued solutions. The Schrödinger equation, derived from the energy-momentum relation E = p²/(2m) + V via the substitutions E → iℏ∂/∂t and p → −iℏ∂/∂x, involves only one time derivative because energy appears linearly in the non-relativistic relation. The factor of i ensures the equation is consistent with unitary time evolution (probability conservation). Because the equation is first-order in time, specifying Ψ(x, 0) alone determines Ψ(x, t) for all future times—but Ψ must be complex to carry both amplitude and phase information simultaneously.

Question 2: BASIC CALCULATION

An electron is confined in a one-dimensional infinite square well of width L = 1.0 nm. Calculate the energy (in eV) of the n = 3 state. Use m_e = 9.109 × 10⁻³¹ kg and ℏ = 1.055 × 10⁻³⁴ J·s.

Show answer

Using E_n = n²π²ℏ²/(2mL²), first compute E₁: E₁ = π²(1.055 × 10⁻³⁴)²/[2(9.109 × 10⁻³¹)(1.0 × 10⁻⁹)²] = (1.098 × 10⁻⁶⁷)/(1.822 × 10⁻⁴⁸) = 6.024 × 10⁻²⁰ J = 0.376 eV. Therefore E₃ = 9 × E₁ = 9 × 0.376 eV = 3.39 eV. This corresponds to a photon in the near-ultraviolet range, illustrating how nanometer-scale confinement produces electronvolt-scale energies.

Question 3: INTERMEDIATE

For the quantum harmonic oscillator with angular frequency ω, compute the expectation value ⟨x²⟩ for the ground state (n = 0) using the known ground-state wave function ψ₀(x) = (mω/πℏ)^1/4 exp(−mωx²/(2ℏ)). Express your answer in terms of ℏ, m, and ω.

Show answer

We need ⟨x²⟩ = ∫ from −∞ to ∞ x²|ψ₀|² dx = (mω/πℏ)^1/2 ∫ x² exp(−mωx²/ℏ) dx. Using the standard Gaussian integral ∫ x²

e^{−αx²}

dx = (1/2)√(π/α³) with α = mω/ℏ: ⟨x²⟩ = (mω/πℏ)^(1/2) × (1/2)√(πℏ³/(m³ω³)) = (1/2)(ℏ/(mω)). Thus ⟨x²⟩ = ℏ/(2mω). This is consistent with the zero-point energy E₀ = ½ℏω = ½mω²⟨x²⟩ + ½⟨p²⟩/m, where by the virial theorem for the harmonic oscillator ⟨T⟩ = ⟨V⟩ = E₀/2.

Question 4: APPLIED

A conjugated dye molecule has π electrons delocalized over a chain of length L = 1.2 nm. Using the free-electron (particle-in-a-box) model, estimate the wavelength of light absorbed when the molecule transitions from its highest occupied energy level to the lowest unoccupied level. Assume 6 π electrons and use the electron mass.

Show answer

With 6 π electrons, by the Pauli exclusion principle the levels n = 1, 2, 3 are doubly occupied (each holds 2 electrons with opposite spin). The HOMO is n = 3 and the LUMO is n = 4. The transition energy is ΔE = E₄ − E₃ = (4² − 3²)π²ℏ²/(2mL²) = 7π²ℏ²/(2mL²). Substituting: ΔE = 7 × (1.098 × 10⁻⁶⁷)/(2 × 9.109 × 10⁻³¹ × 1.44 × 10⁻¹⁸) = 7.686 × 10⁻⁶⁷/(2.624 × 10⁻⁴⁸) = 2.93 × 10⁻¹⁹ J. The wavelength is λ = hc/ΔE = (1.987 × 10⁻²⁵)/(2.93 × 10⁻¹⁹) ≈ 678 nm, which falls in the red region of the visible spectrum—consistent with the observed colors of cyanine dyes. This beautifully demonstrates how the particle-in-a-box model provides semi-quantitative predictions for conjugated molecular systems.

Question 5: CRITICAL THINKING

The time-independent Schrödinger equation Ĥψ = Eψ is an eigenvalue equation. Prove that if ψ₁ and ψ₂ are eigenfunctions of Ĥ with different eigenvalues E₁ ≠ E₂, then ψ₁ and ψ₂ are orthogonal: ∫ψ₁*ψ₂ dx = 0. State clearly where you use the Hermiticity of Ĥ.

Show answer

Start from Ĥψ₁ = E₁ψ₁ and Ĥψ₂ = E₂ψ₂. Take the inner product of the first equation with ψ₂: ∫ψ₂*Ĥψ₁ dx = E₁∫ψ₂*ψ₁ dx. Now use the Hermiticity of Ĥ, which states ∫ψ₂*Ĥψ₁ dx = ∫(Ĥψ₂)*ψ₁ dx = E₂*∫ψ₂*ψ₁ dx. Since Ĥ is Hermitian, its eigenvalues are real, so E₂* = E₂. Equating: E₂∫ψ₂*ψ₁ dx = E₁∫ψ₂*ψ₁ dx, which gives (E₁ − E₂)∫ψ₂*ψ₁ dx = 0. Since E₁ ≠ E₂, the integral must vanish: ∫ψ₂*ψ₁ dx = 0. This orthogonality of eigenfunctions belonging to distinct eigenvalues is fundamental—it ensures that different energy states are distinguishable and that the spectral decomposition of Ĥ is well-defined.

SUMMARY

Lesson Summary

Summary

The Schrödinger equation is the fundamental equation governing the behavior of non-relativistic quantum systems. In its time-dependent form (iℏ ∂Ψ/∂t = ĤΨ), it prescribes how the wave function Ψ(x, t) evolves deterministically under the influence of the Hamiltonian operator Ĥ = T̂ + V̂. The time-independent form (Ĥψ = Eψ) arises when the potential is time-independent, yielding stationary states and discrete energy eigenvalues that explain the quantized spectra of atoms and molecules.

Key exactly solvable systems—the particle in a box, the harmonic oscillator, and the hydrogen atom—illustrate how boundary conditions and potential shapes determine energy spectra and probability densities. The equation's linearity enables superposition and underpins all approximate methods (perturbation theory, variational method, Hartree–Fock). While limited to non-relativistic regimes and fixed particle number, the Schrödinger equation remains the indispensable foundation upon which virtually all of modern quantum chemistry is built.