Appendix A: Complex Numbers in Optics

Complex numbers provide one of the most elegant and powerful mathematical frameworks in all of optics. Far from being abstract mathematical curiosities, they offer profound insights into the nature of light itself. This appendix will guide you through the essential concepts, starting from the very basics and building toward sophisticated applications that illuminate the behavior of electromagnetic waves.

1A.1 What Are Complex Numbers and Why Do We Need Them?¶

1.1A.1.1 The Historical Motivation¶

Complex numbers emerged from a simple but frustrating problem: what is $\sqrt{-1}$? For centuries, mathematicians considered this question meaningless. But in the 16th century, while solving cubic equations, Italian mathematicians discovered that even when seeking real solutions, intermediate steps often required manipulating square roots of negative numbers.

Rather than dismiss these expressions as nonsensical, they began treating $\sqrt{-1}$ as a new kind of number, eventually denoted by $i$ (from “imaginary”). This bold step opened up an entirely new mathematical universe—one that turns out to perfectly describe oscillatory phenomena like light waves.

1.2A.1.2 Definition and Basic Structure¶

A complex number $z$ has the general form:

where $a$ and $b$ are ordinary real numbers, and $i$ is the imaginary unit defined by the fundamental property:

Let’s explore what this means through some examples:

1.3A.1.3 Arithmetic Operations¶

Complex numbers follow familiar arithmetic rules, with one key difference: we must remember that $i^2 = -1$.

Addition and Subtraction: Complex numbers add component-wise, just like vectors:

This makes intuitive sense: real parts combine with real parts, imaginary parts with imaginary parts.

Multiplication: Here’s where things get interesting. Using the distributive property and $i^2 = -1$:

Complex Conjugate: The complex conjugate of $z = a + bi$ is $z^* = a - bi$. This operation flips the sign of the imaginary part and has profound importance in physics—it’s how we extract physically meaningful quantities from complex expressions.

1.4A.1.4 Geometric Representation: The Complex Plane¶

One of the most illuminating ways to understand complex numbers is through their geometric representation in the complex plane. This is a coordinate system where:

• The horizontal axis represents the real part

• The vertical axis represents the imaginary part

• Each complex number $z = a + bi$ corresponds to the point $(a, b)$

This geometric view immediately reveals why complex numbers are so powerful for describing rotations and oscillations—operations that are fundamental to understanding waves.

1.5A.1.5 Polar Form: Magnitude and Phase¶

Every complex number can be expressed in polar form, which separates its magnitude from its direction:

where:

• $r = |z| = \sqrt{a^2 + b^2}$ is the magnitude or modulus

• $\theta = \arg(z) = \operatorname{atan2}(b,a)$ is the argument or phase when $z\ne0$; it is undefined at $z=0$

The polar form is particularly powerful because:

1. Multiplication becomes addition of phases: $r_1 e^{i\theta_1} \cdot r_2 e^{i\theta_2} = r_1 r_2 e^{i(\theta_1 + \theta_2)}$

2. Powers become simple: $(re^{i\theta})^n = r^n e^{in\theta}$

3. Physical meaning emerges: $r$ represents amplitude, $\theta$ represents phase

2A.2 Euler’s Formula: The Bridge Between Exponentials and Trigonometry¶

2.1A.2.1 The Most Beautiful Equation in Mathematics¶

Euler’s formula stands as one of the most remarkable discoveries in mathematics:

This equation creates a bridge between three seemingly unrelated mathematical concepts: exponential functions, trigonometry, and complex numbers. But why is this true, and why should we care?

2.2A.2.2 Understanding Why Euler’s Formula Works¶

2.3A.2.3 Physical Interpretation: Rotation in the Complex Plane¶

Euler’s formula has a beautiful geometric meaning: $e^{i\theta}$ represents a point on the unit circle at angle $\theta$ from the positive real axis. As $\theta$ increases, $e^{i\theta}$ traces out circular motion in the complex plane.

This insight is crucial for optics because light waves are oscillatory phenomena. The phase $\theta$ in Euler’s formula will correspond to the phase of an electromagnetic wave, while the magnitude can represent the amplitude.

2.4A.2.4 Special Cases and Identities¶

Several important special cases emerge from Euler’s formula:

2.5A.2.5 Properties of Complex Exponentials¶

Complex exponentials inherit wonderful properties from regular exponentials:

1. Multiplication Rule: $e^{i\alpha} \cdot e^{i\beta} = e^{i(\alpha + \beta)}$

2. Division Rule: $\frac{e^{i\alpha}}{e^{i\beta}} = e^{i(\alpha - \beta)}$

3. Power Rule: $(e^{i\theta})^n = e^{in\theta}$

4. Conjugate Rule: $(e^{i\theta})^* = e^{-i\theta}$

These properties make complex exponentials incredibly convenient for calculations involving waves and oscillations.

3A.3 Complex Numbers in Wave Physics¶

3.1A.3.1 Why Waves and Complex Numbers Are Perfect Partners¶

Before diving into optical applications, let’s understand why complex numbers are so naturally suited for describing waves. Consider a simple harmonic oscillator—perhaps a mass on a spring, or the electric field in a light wave. The position (or field strength) varies sinusoidally with time:

While this real function perfectly describes the physical motion, working with it mathematically can be cumbersome. Derivatives of cosines give sines, integrals mix sines and cosines, and keeping track of phase relationships becomes complex (pun intended).

3.2A.3.2 Electromagnetic Waves in Complex Notation¶

In optics, we deal with electromagnetic waves—oscillating electric and magnetic fields. A plane wave traveling in the positive $z$-direction can be written as:

Let’s unpack this expression:

The argument of the exponential, $(kz - \omega t + \phi)$, is called the phase of the wave. This single expression contains all the information about how the wave varies in space and time.

3.3A.3.3 Physical Fields from Complex Expressions¶

The actual, measurable electric field is always real. We extract it by taking the real part of our complex expression:

If $\vec{E}_0 = E_0 \hat{x}$ (linearly polarized in the $x$-direction), then:

This is the familiar sinusoidal wave we expect for light.

4A.4 Interference: Where Complex Numbers Shine¶

4.1A.4.1 The Superposition Principle¶

One of the most important principles in optics is superposition: when two or more waves overlap, the total field is simply the sum of the individual fields. With complex notation, this becomes beautifully simple.

Consider two waves with the same frequency but different amplitudes and phases:

The total field is:

4.2A.4.2 Intensity and the Complex Conjugate¶

The intensity of light is proportional to the time-averaged square of the real electric field. If $E_{\rm phys}(t)=\operatorname{Re}[\mathcal E e^{-i\omega t}]$, with complex peak amplitude $\mathcal E$, then

in a nonmagnetic, lossless medium of refractive index $n$. Thus relative intensities in the same medium are proportional to $|\mathcal E|^2$; the factor $1/2$ comes from averaging the real oscillation.

For our two-wave example, use the phasors $\mathcal E_1=E_{01}e^{i\phi_1}$ and $\mathcal E_2=E_{02}e^{i\phi_2}$. Their common space-time factor cancels in the modulus:

The last term is the interference term—it’s what creates the bright and dark fringes in interference patterns. Without complex numbers, deriving this result would involve tedious trigonometric identities.

5A.5 Complex Refractive Index: Absorption and Dispersion¶

5.1A.5.1 Beyond Simple Refraction¶

When light propagates through most materials, two important effects occur:

1. The phase velocity changes (refraction)

2. The wave weakens (absorption)

A complex refractive index describes phase propagation and attenuation at a given frequency:

where:

• $n$ is the familiar refractive index (affects phase velocity)

• $\kappa$ is the extinction coefficient (causes absorption); we reserve $k$ for wave number

5.2A.5.2 Wave Propagation in Absorbing Media¶

When a wave propagates through a medium with complex refractive index $\tilde{n}$, the wave number becomes complex:

The wave in the medium becomes:

5.3A.5.3 Physical Examples¶

Transparent glass and water have relatively small $\kappa$ in much of the visible spectrum, so a thin sample can transmit light. Metals such as gold and silver have substantial extinction and also reflect strongly at their surfaces. Numerical optical constants depend on wavelength, temperature, and the sample; they should be taken from a specified measured dataset when a quantitative calculation is needed.

6A.6 Advanced Applications in Optics¶

6.1A.6.1 Fresnel Coefficients and Reflection¶

When light hits an interface between two media, some reflects and some transmits. The Fresnel coefficients describe these processes and are naturally expressed using complex numbers.

For s-polarized light (electric field perpendicular to the plane of incidence):

where angles are related by Snell’s law: $n_1\sin\theta_1 = n_2\sin\theta_2$.

6.2A.6.2 Polarization States¶

Complex numbers provide an elegant description of polarization. A general elliptically polarized wave can be written as:

The polarization state is determined by:

• Amplitude ratio: $E_y/E_x$

• Phase difference: $\delta = \phi_y - \phi_x$

6.3A.6.3 Fourier Optics: Spatial Frequencies¶

Modern optics relies heavily on Fourier analysis, which is built on complex exponentials. The spatial Fourier transform of a function $f(x,y)$ is:

This mathematical framework enables:

• Diffraction analysis: Understanding how apertures affect light

• Image formation: How lenses create images

• Holography: Recording and reconstructing wavefronts

• Optical information processing: Manipulating spatial frequencies

7A.7 Computational Techniques¶

7.1A.7.1 Time Averaging with Complex Fields¶

A common calculation in optics involves time averaging oscillatory quantities. For the real field $E_{\rm phys}(t)=\operatorname{Re}[E_0e^{-i\omega t}]$:

The terms oscillating at $2\omega$ average to zero. By contrast, $|E_0e^{-i\omega t}|^2=|E_0|^2$ is the squared modulus of the complex representation and does not itself perform this real-field average.

7.2A.7.2 Differentiation in the Frequency Domain¶

Time and spatial derivatives become algebraic operations:

• $\frac{\partial}{\partial t}[e^{-i\omega t}] = -i\omega e^{-i\omega t}$

• $\frac{\partial}{\partial z}[e^{ikz}] = ik e^{ikz}$

This transforms differential equations (like Maxwell’s equations) into algebraic equations, greatly simplifying analysis.

7.3A.7.3 Phasor Diagrams¶

Complex numbers can be visualized as rotating vectors (phasors) in the complex plane. This visualization is particularly powerful for:

8A.8 Why Complex Numbers Work So Well in Optics¶

8.1A.8.1 Mathematical Elegance¶

Complex numbers work beautifully in optics because they match the mathematical structure of wave phenomena:

1. Linearity: Both wave equations and complex arithmetic are linear

2. Superposition: Adding waves ↔ Adding complex numbers

3. Harmonic motion: Natural connection to $e^{i\omega t}$

4. Phase relationships: Automatically preserved in calculations

8.2A.8.2 Physical Insight¶

Beyond mathematical convenience, complex numbers often reveal physical insights:

8.3A.8.3 Computational Power¶

Modern optics relies heavily on numerical computation. Complex numbers enable:

• FFT algorithms: Fast Fourier transforms for diffraction calculations

• Matrix methods: Describing optical systems using complex matrices

• Beam propagation: Numerical solution of wave equations

• Nonlinear optics: Analysis of intensity-dependent phenomena

9A.9 Common Pitfalls and How to Avoid Them¶

9.1A.9.1 Physical vs. Mathematical Fields¶

9.2A.9.2 Phase Conventions¶

9.3A.9.3 Complex Conjugation in Intensity Calculations¶

10A.10 Summary and Looking Forward¶

Complex numbers are not mere mathematical conveniences in optics—they are fundamental tools that reveal the deep structure of electromagnetic phenomena. They provide:

As you progress through this textbook, you’ll encounter complex numbers in virtually every advanced topic:

• Fourier optics: Spatial frequency analysis

• Laser physics: Gain and cavity analysis

• Nonlinear optics: Phase matching and frequency conversion

• Quantum optics: Photon statistics and coherence

• Metamaterials: Effective medium theory

11A.11 Worked Examples¶

11.1A.11.1 Young’s Double-Slit Interference¶

Solution:

Let’s place the two slits at $y = \pm d/2$ and consider a point $P$ on the screen at height $y$. The path difference between rays from the two slits is approximately:

$\Delta = \frac{yd}{L}$

for small angles.

The fields from the two slits at point $P$ are: $E_1 = E_0 e^{i(kL - \omega t)}$ $E_2 = E_0 e^{i(kL + k\Delta - \omega t)} = E_0 e^{i(kL - \omega t)} e^{ikyd/L}$

The total field is: $E_{\text{total}} = E_0 e^{i(kL - \omega t)}[1 + e^{ikyd/L}]$

To find the intensity, we calculate $|E_{\text{total}}|^2$: $|E_{\text{total}}|^2 = |E_0|^2 |1 + e^{ikyd/L}|^2$

Using $|1 + e^{i\phi}|^2 = (1 + e^{i\phi})(1 + e^{-i\phi}) = 2 + 2\cos\phi$:

If $I_s=(n\varepsilon_0c/2)|E_0|^2$ is the intensity from either source alone at the screen, then $I(y)=2I_s[1+\cos(kyd/L)] =4I_s\cos^2\left(\frac{kyd}{2L}\right)$.

11.2A.11.2 Transmission Through an Absorbing Slab¶

Solution:

This problem involves multiple physical processes:

1. Reflection at the first interface

2. Absorption within the slab

3. Reflection at the second interface

Step 1: First interface (air to glass) The Fresnel reflection coefficient for normal incidence is: $r_{12} = \frac{1 - \tilde{n}}{1 + \tilde{n}} = \frac{1 - (1.5 + 0.01i)}{1 + (1.5 + 0.01i)} = \frac{-0.5 - 0.01i}{2.5 + 0.01i}$

To compute this, multiply numerator and denominator by the complex conjugate of the denominator: $r_{12} = \frac{(-0.5 - 0.01i)(2.5 - 0.01i)}{(2.5 + 0.01i)(2.5 - 0.01i)} = \frac{-1.2501 - 0.0200i}{6.2501} \approx -0.200013 - 0.003200i$

The reflectance is $R_1=|r_{12}|^2\approx0.040015$.

The transmitted field coefficient at the front surface is $t_{12}=1+r_{12}\approx0.799987-0.003200i$.

Step 2: Propagation through the slab The wave number in the material is: $\tilde{k} = \frac{2\pi}{\lambda_0}\tilde{n} = \frac{2\pi}{\lambda_0}(1.5 + 0.01i)$

After propagating a distance $t$, the field is multiplied by: $e^{i\tilde{k}t} = e^{i(2\pi/\lambda_0)(1.5 + 0.01i)t} = e^{i(2\pi/\lambda_0)(1.5t)} e^{-(2\pi/\lambda_0)(0.01t)}$

The absorption factor is $e^{-(2\pi/\lambda_0)(0.01t)}$.

Step 3: Second interface (glass to air) The reverse field-reflection coefficient is $r_{21}=-r_{12}$. For this small extinction coefficient, approximate the exit-surface power loss by the same reflectance $R_2\approx R_1$.

Step 4: Total transmission In the stated single-pass approximation, the transmitted intensity fraction is $T\approx(1-R_1)(1-R_2)e^{-(4\pi\cdot0.01)t/\lambda_0} \approx0.9216e^{-0.1257t/\lambda_0}$.

11.3A.11.3 Circular Polarization Analysis¶

Solution:

Initial state: Unpolarized light can be represented as an incoherent mixture of all polarization states.

After linear polarizer at 45°: The transmitted field is linearly polarized. Let $E_0$ denote its peak amplitude after the polarizer: $\vec{E}_1 = E_0 \frac{(\hat{x} + \hat{y})}{\sqrt{2}}$

Quarter-wave plate action: A quarter-wave plate introduces a phase difference of $\pi/2$ between its fast and slow axes. With the fast axis along x, it multiplies the y-component by $e^{i\pi/2} = i$:

$\vec{E}_2 = E_0 \frac{(\hat{x} + i\hat{y})}{\sqrt{2}}$

Analysis: This is left-handed circular polarization! We can verify this by checking the time dependence:

$\vec{E}(t) = \text{Re}\left[E_0 \frac{(\hat{x} + i\hat{y})}{\sqrt{2}} e^{-i\omega t}\right]$ $= E_0 \frac{(\cos\omega t \hat{x} + \cos(\omega t - \pi/2)\hat{y})}{\sqrt{2}}$ $= E_0 \frac{(\cos\omega t \hat{x} + \sin\omega t \hat{y})}{\sqrt{2}}$

At $t = 0$: $\vec{E} = (E_0/\sqrt{2})\hat{x}$ At $t = \pi/(2\omega)$: $\vec{E} = (E_0/\sqrt{2})\hat{y}$ At $t = \pi/\omega$: $\vec{E} = -(E_0/\sqrt{2})\hat{x}$ At $t = 3\pi/(2\omega)$: $\vec{E} = -(E_0/\sqrt{2})\hat{y}$

The electric field vector rotates counterclockwise when viewed against the direction of propagation—this is left-handed circular polarization.

12A.12 Practice Problems¶

12.1A.12.1 Basic Complex Arithmetic¶

12.2A.12.2 Wave Interference¶

12.3A.12.3 Complex Refractive Index¶

12.4A.12.4 Polarization¶

13A.13 Historical Notes and Further Reading¶

13.1A.13.1 Historical Development¶

The development of complex numbers and their application to physics is a fascinating story of mathematical evolution:

Historical Timeline:

• 1545: Cardano’s “Ars Magna” first mentions square roots of negative numbers

• 1777: Euler introduces the notation $i$ for $\sqrt{-1}$

• 1799: Gauss proves the fundamental theorem of algebra using complex numbers

• 1806: Argand develops the geometric representation (complex plane)

• 1864: Maxwell’s equations are formulated, setting the stage for complex analysis in optics

• 1920s: Complex amplitudes become central to modern quantum mechanics

13.2A.13.2 Connections to Other Fields¶

Complex numbers in optics connect to many other areas of physics and mathematics:

13.3A.13.3 Recommended Further Reading¶

Further Reading:

Born, M., & Wolf, E. (1999). Principles of Optics (7th ed.). Cambridge University Press.

The definitive reference for classical optics, with extensive use of complex analysis.

Goodman, J. W. (2005). Introduction to Fourier Optics (3rd ed.). Roberts & Company.

Essential reading for understanding how complex analysis enables modern optical system design.

Hecht, E. (2017). Optics (5th ed.). Pearson.

Comprehensive undergraduate text with clear explanations of complex number applications.

Fowles, G. R. (1989). Introduction to Modern Optics (2nd ed.). Dover Publications.

Classic text emphasizing the wave nature of light and mathematical methods.

Jackson, J. D. (1999). Classical Electrodynamics (3rd ed.). Wiley.

Advanced treatment of electromagnetic theory with sophisticated use of complex analysis.

13.4A.13.4 Modern Computational Tools¶

Today’s optical research relies heavily on computational tools that leverage complex number arithmetic:

14A.14 Final Thoughts: The Beauty of Mathematical Physics¶

As you’ve seen throughout this appendix, complex numbers aren’t just computational tools—they’re windows into the fundamental nature of wave phenomena. The fact that $e^{i\pi} + 1 = 0$ connects the most important constants in mathematics is no accident. It reflects deep symmetries in the mathematical structures that govern our physical world.

In optics, these mathematical structures manifest as:

• Wave-particle duality: The complex wavefunction that bridges classical waves and quantum particles

• Symmetry principles: The invariances that lead to conservation laws

• Analytical properties: The mathematical constraints that ensure physical causality

As you continue your journey through optics, remember that the mathematics isn’t separate from the physics—it is the physics, expressed in its most elegant and universal form. Complex numbers are your passport to this deeper understanding.