Complex Numbers & Euler's Formula — Complex numbers and...

SECTION 1

Historical Context & Motivation

The story of complex numbers begins not in quantum mechanics but in Renaissance algebra, where mathematicians first encountered the square root of negative one while solving cubic equations. Initially dismissed as "impossible" or "imaginary," these numbers gradually proved indispensable across mathematics and physics. The pivotal moment for physical chemistry arrived in the early twentieth century, when Schrödinger, Heisenberg, and Dirac realized that the fundamental equation governing atomic and molecular behavior required complex-valued solutions—not as a mathematical convenience, but as an intrinsic feature of nature.

Understanding why complex numbers are necessary in quantum mechanics, rather than merely tolerated, is essential for any physical chemist working with spectroscopy, molecular orbital theory, or statistical thermodynamics. Euler's formula provides the bridge between exponential and trigonometric representations of oscillatory phenomena, making it the single most important identity in quantum chemistry's mathematical toolkit.

1545

Cardano's Ars Magna

Gerolamo Cardano publishes solutions to cubic equations involving √(−1), calling them "sophistic" quantities. These are the earliest recorded manipulations of imaginary numbers in formal mathematics.

1748

Euler's Formula Published

Leonhard Euler establishes the identity e^iθ = cos θ + i sin θ in his Introductio in analysin infinitorum, unifying exponential growth with circular motion.

1926

Schrödinger's Wave Equation

Erwin Schrödinger formulates the time-dependent wave equation with an explicit factor of i, making complex numbers an irreducible part of quantum theory rather than a computational shortcut.

1930

Dirac's Quantum Mechanics Textbook

Paul Dirac's Principles of Quantum Mechanics formalizes the bra-ket notation in complex Hilbert space, cementing the role of complex inner products in computing expectation values and transition probabilities.

2022

Experimental Proof of Necessity

Experiments by Chen et al. (Nature, 2022) demonstrate that no real-valued formulation of quantum mechanics can reproduce all observed correlations, confirming that complex numbers are physically essential, not merely convenient.

The central question this lesson addresses is both practical and conceptual: How do complex numbers and Euler's formula enable us to represent, manipulate, and extract physical predictions from quantum mechanical wavefunctions? By the end, you will be able to move fluently between Cartesian, polar, and exponential forms of complex numbers—and apply these representations to problems in spectroscopy, tunneling, and molecular orbital theory.

SECTION 2

Core Principles & Definitions

Before diving into quantum applications, we must establish the algebraic and geometric foundations of complex arithmetic. A complex number z is written z = a + bi, where a is the real part Re(z), b is the imaginary part Im(z), and i = √(−1) satisfies i² = −1. This Cartesian form maps directly onto a two-dimensional plane—the Argand plane—where the horizontal axis represents the real component and the vertical axis the imaginary component.

Complex Conjugate (z*)

Given z = a + bi, the complex conjugate is z* = a − bi. Multiplying z × z* yields |z|² = a² + b², always a real, non-negative number. In QM, ψ*ψ gives the probability density.

Modulus and Phase

The modulus |z| = √(a² + b²) gives the distance from the origin; the phase θ = arctan(b/a) gives the angle from the positive real axis. Together they define the polar form z = |z|e^iθ.

Euler's Formula

e^iθ = cos θ + i sin θ connects the complex exponential to trigonometric functions. Setting θ = π yields Euler's identity: e^iπ + 1 = 0.

Born Interpretation

The wavefunction ψ(x,t) is complex-valued, but the Born rule states that the physical observable—probability density—is |ψ|² = ψ*ψ, guaranteeing a real, non-negative quantity that integrates to unity over all space.

Time Evolution Factor

Stationary states carry a time-dependent phase factor e^−iEt/ℏ. This complex exponential encodes oscillation without changing the probability density, since |e^−iEt/ℏ|² = 1.

✦ KEY TAKEAWAY

Think of a complex number as a vector with both a length (modulus) and a direction (phase). In quantum mechanics, the length squared tells you how much probability is present, while the phase tells you how wavefunctions interfere—constructively or destructively—much like two ocean waves that can reinforce or cancel each other depending on their relative crest positions. Euler's formula is the Rosetta Stone translating between the oscillatory (trigonometric) and exponential languages used to describe these phases.

SECTION 3

Visualizing Complex Numbers in the Argand Plane

The Argand plane (also called the complex plane) provides geometric intuition for every algebraic operation on complex numbers. In quantum chemistry, this visualization is indispensable: multiplication by e^iθ corresponds to a rotation of the complex vector by angle θ, while multiplication by a real scalar r stretches or compresses the vector's length. Time evolution of a quantum state is literally a rotation in the complex plane at angular frequency ω = E/ℏ.

The complex number z = a + bi is plotted as a point in the Argand plane. The cyan vector represents z with modulus |z| and phase angle θ (shown in amber). The green dashed vector shows the complex conjugate z* = a − bi, reflected across the real axis. Multiplying z by z* gives |z|², which in QM corresponds to the probability density.

Several features of this diagram are worth noting for quantum applications. First, the modulus |z| is invariant under phase rotations: if you multiply z by e^iφ, the vector rotates but its length remains unchanged. This is precisely why the probability density |ψ|² does not change under time evolution for stationary states—the time-dependent phase factor e^−iEt/ℏ is a pure rotation. Second, complex conjugation reflects the vector across the real axis, flipping the sign of the imaginary part. When two complex amplitudes are multiplied—as in the overlap integral ⟨ψ₁|ψ₂⟩—their phases interact, producing constructive or destructive interference depending on whether the net phase difference is near 0 or π.

SECTION 4

Mathematical Framework

The mathematical machinery connecting complex numbers to quantum mechanics centers on three key results: Euler's formula itself, the time-dependent Schrödinger equation, and the inner product in complex Hilbert space. We present each with its physical interpretation relevant to physical chemistry applications.

EULER'S FORMULA

e^(iθ) = cos θ + i sin θ

Where θ is a real number (the phase angle in radians), e is Euler's number (≈ 2.71828), and i = √(−1). This identity is derived by comparing the Taylor series of e^iθ with those of cos θ and sin θ. It allows any sinusoidal oscillation to be written as a complex exponential, dramatically simplifying differential equations in quantum mechanics.

POLAR FORM OF A COMPLEX NUMBER

z = |z| e^(iθ) = r(cos θ + i sin θ)

Where r = |z| = √(a² + b²) is the modulus and θ = arg(z) = arctan(b/a) is the argument. This representation is powerful because multiplication becomes addition of phases: z₁z₂ = r₁r₂ e^{i(θ₁+θ₂)}.

TIME-DEPENDENT SCHRÖDINGER EQUATION

iℏ ∂ψ/∂t = Ĥψ

The explicit factor of i on the left side is not optional—it guarantees that the time evolution operator Û(t) = e^−iĤt/ℏ is unitary (probability-preserving). For a stationary state with energy E, the solution takes the form ψ(x,t) = φ(x) e^−iEt/ℏ, where φ(x) satisfies the time-independent Schrödinger equation.

COMPLEX INNER PRODUCT (BORN RULE)

⟨ψ₁|ψ₂⟩ = ∫ ψ₁*(x) ψ₂(x) dx

The complex conjugate ψ₁* in the integrand ensures that the norm ⟨ψ|ψ⟩ = ∫|ψ|² dx is real and non-negative. This is the Born rule: |ψ(x)|² dx gives the probability of finding the particle in [x, x+dx]. Overlap integrals ⟨ψ₁|ψ₂⟩ determine transition amplitudes and selection rules in spectroscopy.

💡 Why i in the Schrödinger Equation?

If we replaced i with 1 in iℏ ∂ψ/∂t = Ĥψ, we would get a real diffusion equation (like the heat equation), where ψ would decay exponentially to zero. The factor of i transforms exponential decay into oscillation (via Euler's formula), ensuring that probability is conserved for all time. This is the mathematical reason quantum mechanics requires complex numbers—not real ones.

SECTION 5

Complex Numbers in Quantum Chemistry Applications

In physical chemistry, complex numbers appear in several crucial contexts beyond the abstract Schrödinger equation. The following diagram illustrates how Euler's formula connects to the particle-in-a-ring model—a prototype for aromatic π-electron systems in organic molecules such as benzene. In this system, the wavefunctions are pure complex exponentials, and the quantization of angular momentum arises directly from the requirement that the wavefunction be single-valued after a full rotation.

Left: A particle constrained to a ring of radius R. The wavefunction ψ_m(φ) = (2π)^−1/2 e^imφ must be single-valued: ψ(φ+2π) = ψ(φ), which forces m to be an integer (0, ±1, ±2, …). Right: The resulting energy levels E_m = m²ℏ²/2I, each doubly degenerate for |m| > 0, modeling the Hückel π-electron levels of aromatic rings.

The single-valuedness requirement is a direct consequence of Euler's formula. After a full revolution, e^im(φ+2π) = e^imφ × e^i2πm. For this to equal e^imφ, we need e^i2πm = 1, which by Euler's formula requires cos(2πm) + i sin(2πm) = 1, and this holds if and only if m is an integer. This is a beautiful example of how quantization emerges from the topology of the complex exponential rather than from any ad hoc postulate.

Common QM contexts where complex numbers play an essential role

QM Context	Where Complex Numbers Appear	Physical Consequence
Particle on a Ring	ψ_m = (2π)^−1/2 e^imφ	Angular momentum quantization; Hückel MOs for aromatic molecules
Hydrogen Atom (Φ part)	Φ_m(φ) = (2π)^−1/2 e^imφ	Magnetic quantum number m_l; orbital angular momentum projection
Time Evolution	ψ(x,t) = φ(x) e^−iEt/ℏ	Stationary states oscillate in phase; probability density \|ψ\|² is time-independent
Perturbation Theory	Transition amplitudes c_fi = ⟨f\|V\|i⟩/(E_i−E_f)	Complex matrix elements determine transition probabilities and selection rules in spectroscopy
Fourier Transforms	f̃(k) = ∫ f(x) e^−ikx dx	Momentum-space representation of wavefunctions; uncertainty principle

SECTION 6

Worked Example: Probability Density from a Superposition State

Consider a particle on a ring in the superposition state Ψ(φ,t) = (1/√2)[ψ₁(φ)e^{−iE₁t/ℏ} + ψ₋₁(φ)e^{−iE₋₁t/ℏ}], where ψ_m = (2π)^−1/2 e^imφ. Since E₁ = E₋₁ = ℏ²/2I, compute the probability density |Ψ|² and show it is time-independent.

Probability Density of a Ring Superposition State

Step 1 — Write the Full WavefunctionSubstitute the explicit forms of ψ₁ and ψ₋₁ into the superposition. Since E₁ = E₋₁ ≡ E, the time factor is the same for both terms: Ψ(φ,t) = (1/√2) × (2π)^−1/2 [e^iφ + e^−iφ] × e^−iEt/ℏ

Ψ(φ,t) = (1/√π) cos φ × e^−iEt/ℏ

Step 2 — Apply Euler's FormulaUse Euler's formula to simplify the sum of complex exponentials: e^iφ + e^−iφ = (cos φ + i sin φ) + (cos φ − i sin φ) = 2 cos φ. Thus the spatial part simplifies to a purely real standing wave.

e^iφ + e^−iφ = 2 cos φ

Step 3 — Compute |Ψ|²Take Ψ*Ψ. The complex conjugate flips the sign of i in the time factor: |Ψ|² = Ψ* × Ψ = (1/√π) cos φ × e^+iEt/ℏ × (1/√π) cos φ × e^−iEt/ℏ The time-dependent exponentials cancel because e^+iEt/ℏ × e^−iEt/ℏ = e⁰ = 1.

|Ψ(φ,t)|² = (1/π) cos²φ — time-independent

Step 4 — Verify NormalizationIntegrate over the full ring (0 to 2π): ∫₀^2π (1/π) cos²φ dφ = (1/π) × π = 1. The wavefunction is properly normalized, confirming total probability equals unity.

∫₀²π |Ψ|² dφ = 1 ✓

Step 5 — Physical InterpretationThe superposition of two counter-rotating complex exponentials (m = +1 and m = −1) produces a real standing wave proportional to cos²φ. The probability density is peaked at φ = 0 and φ = π and vanishes at φ = π/2 and 3π/2. This is directly analogous to how the real p_x orbital of the hydrogen atom is formed by combining the m = +1 and m = −1 spherical harmonics. The time-independence results from the degeneracy E₁ = E₋₁, so both components rotate at the same rate in the complex plane and the interference pattern is stationary.

SECTION 7

Comparing Complex Number Representations

Choosing the right representation of a complex number—Cartesian, polar, or exponential—can dramatically simplify a calculation. Each form has strengths suited to particular operations common in quantum chemistry. The table below compares these representations and identifies when each is most efficient.

Three equivalent representations of a complex number and their operational strengths

Representation	Form	Best For	Limitations
Cartesian	z = a + bi	Addition/subtraction; separating real and imaginary parts; computing expectation values	Multiplication and division are algebraically tedious; obscures phase information
Polar / Trigonometric	z = r(cos θ + i sin θ)	Visualizing phase and amplitude; de Moivre's theorem for integer powers	Notation is verbose; less compact for differential equations
Exponential (Euler)	z = re^iθ	Multiplication, division, powers, roots; solving ODEs; time evolution in QM; Fourier transforms	Addition requires converting back to Cartesian form; can obscure individual components

✦ KEY TAKEAWAY

Think of the three representations as different coordinate systems for the same map. Just as you might use Cartesian coordinates for grid navigation and polar coordinates for radar sweeps, in quantum chemistry you should use the Cartesian form for addition (e.g., superposing wavefunctions), the exponential form for multiplication (e.g., applying time evolution operators), and the polar form for visualization (e.g., sketching orbital phase diagrams). Fluency in switching between them is a core skill for physical chemistry problem-solving.

SECTION 8

Connection to Advanced Quantum Theory

The algebra of complex numbers introduced in this lesson is the foundation upon which more sophisticated mathematical structures in quantum chemistry are built. As you progress into graduate-level topics, you will encounter complex Hilbert spaces, unitary operators, and Hermitian matrices—all of which rely on the properties of complex conjugation, modulus, and phase that we have developed here.

How complex number concepts extend into advanced quantum chemistry

Concept in This Lesson	Advanced Extension	Where It Appears
Complex conjugate z*	Hermitian adjoint A† (transpose + complex conjugate of a matrix)	Observable operators must satisfy A = A† to guarantee real eigenvalues
Modulus \|z\|² = zz*	Trace of density matrix ρ = \|ψ⟩⟨ψ\|	Mixed states in statistical mechanics; decoherence theory
Phase factor e^iθ	Unitary time evolution Û = e^−iĤt/ℏ	Quantum dynamics; propagator methods; path integrals
Inner product ⟨ψ₁\|ψ₂⟩	Complex-valued matrix elements ⟨n\|Â\|m⟩	Transition dipole moments; Fermi's golden rule; selection rules in spectroscopy
e^imφ on a ring	Spherical harmonics Y_l^m(θ,φ)	Atomic orbitals; angular momentum coupling; term symbols

One particularly important forward-looking connection is the relationship between complex exponentials and group theory. The set of complex numbers with unit modulus, {e^iθ : θ ∈ [0, 2π)}, forms the unitary group U(1), which is the symmetry group underlying electromagnetic gauge invariance. In molecular symmetry, the irreducible representations of cyclic groups C_n are built from the n-th roots of unity—complex numbers of the form e^i2πk/n. Thus, the character tables you use to determine allowed IR and Raman transitions are deeply rooted in the complex arithmetic developed in this lesson.

SECTION 9

Practice Problems

PROBLEM 1 — CONCEPTUAL

The time-dependent Schrödinger equation contains an explicit factor of i. If we replaced i with 1 (making the equation real), what type of physical behavior would the solutions describe instead of quantum oscillation? Explain why this change would violate conservation of probability.

PROBLEM 2 — BASIC CALCULATION

Convert the complex number z = −3 + 3i to polar (exponential) form re^iθ. Express θ in radians as an exact fraction of π.

PROBLEM 3 — INTERMEDIATE

A particle on a ring has wavefunction ψ₂(φ) = (2π)^−1/2 e^2iφ. (a) Show that ψ₂ is an eigenfunction of the angular momentum operator L̂_z = −iℏ d/dφ and find the eigenvalue. (b) Compute ψ₂*ψ₂ and interpret the result physically.

PROBLEM 4 — APPLIED

In the Hückel model of benzene (6 carbon atoms in a ring), the MO energies are given by E_k = α + 2β cos(2πk/6) for k = 0, 1, …, 5, where α and β are the Coulomb and resonance integrals (β < 0). Using Euler's formula, show that the six wavefunctions can be written as φ_k(j) = (1/√6) e^i2πkj/6, and identify which values of k correspond to degenerate energy levels.

PROBLEM 5 — CRITICAL THINKING

A recent paper (Chen et al., Nature, 2022) demonstrated that no real-valued formulation of quantum mechanics can reproduce all experimentally observed correlations. Consider a two-state system with basis {|0⟩, |1⟩} and Hamiltonian H with complex off-diagonal elements H₀₁ = Δe^iδ. (a) Show that Hermiticity requires H₁₀ = Δe^−iδ. (b) Argue why the phase δ cannot, in general, be removed by a change of basis, and discuss what physical observable might depend on this irreducible complex phase.

SUMMARY

Summary & Key Concepts

Complex numbers z = a + bi are not an abstract curiosity but the native language of quantum mechanics. Every wavefunction is complex-valued, and Euler's formula e^iθ = cos θ + i sin θ provides the essential bridge between oscillatory (trigonometric) and exponential representations. The modulus |z| determines amplitudes and probabilities via the Born rule (|ψ|² = ψ*ψ), while the phase θ governs interference between quantum states.

In physical chemistry, these tools appear everywhere: from the time evolution factor e^−iEt/ℏ that makes stationary states oscillate without changing probability, to the particle-on-a-ring model where single-valuedness of e^imφ enforces angular momentum quantization, to Hückel MO theory where discrete Fourier modes on a ring give the π-orbital energies of aromatic molecules. Mastery of the three representations—Cartesian for addition, exponential for multiplication and time evolution, and polar for visualization—equips you to tackle the full spectrum of quantum chemistry problems with confidence and efficiency.