Polarization of Light

Introduction: The Nature of Polarization

Light is a transverse electromagnetic wave. The electric field oscillates perpendicular to the propagation direction. Optical effects depend almost entirely on the electric field oscillation. The way the electric field oscillates in space and time defines the polarization state of light.

Polarization is fundamentally important in photonics: - Optical communication: Polarization-division multiplexing (PDM) increases bandwidth - Display technology: Liquid crystal displays (LCDs) rely on polarization control - Stress analysis: Photoelasticity uses polarized light to visualize strain - Quantum optics: Single photons carry quantum information in polarization - Remote sensing: Depolarization indicates atmospheric properties

This lecture develops the mathematical formalism to describe and manipulate light polarization.

Electromagnetic Waves and Polarization

The Wave Equation from Maxwell’s Equations

Starting from Maxwell’s equations in vacuum:

\[\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}\] \[\nabla \times \mathbf{B} = \mu_0 \epsilon_0 \frac{\partial \mathbf{E}}{\partial t}\]

Taking the curl of the first equation yields:

\[\nabla^2 \mathbf{E} = \mu_0 \epsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2} \tag{1}\]

This is the wave equation, with wave speed \(c = 1/\sqrt{\mu_0 \epsilon_0}\).

Plane Wave Solutions

For a plane wave propagating in the \(+z\) direction:

\[\mathbf{E}(z,t) = \mathbf{E}_0 \cos(kz - \omega t + \phi)\]

where \(k = 2\pi/\lambda\) is the wavenumber and \(\omega = 2\pi f\) is the angular frequency.

The electric field has transverse components:

\[E_x(z,t) = E_{0x} \cos(kz - \omega t + \phi_x)\] \[E_y(z,t) = E_{0y} \cos(kz - \omega t + \phi_y) \tag{2}\]

At a fixed position (\(z = 0\)):

\[E_x(t) = E_{0x} \cos(\omega t - \phi_x)\] \[E_y(t) = E_{0y} \cos(\omega t - \phi_y)\]

The relative phase \(\delta = \phi_y - \phi_x\) determines the polarization state.

Polarization States

Linear Polarization

When \(\delta = 0\) or \(\delta = \pi\):

\[E_y = \frac{E_{0y}}{E_{0x}} E_x\]

The electric field vector oscillates along a straight line in the \(xy\) plane.

Circular Polarization

When \(E_{0x} = E_{0y} = E_0\) and \(\delta = \pm \pi/2\):

\[E_x(t) = E_0 \cos(\omega t)\] \[E_y(t) = E_0 \sin(\omega t) \quad \text{(right circular)}\]

The electric field vector rotates with constant magnitude \(E_0\).

Elliptical Polarization

For arbitrary amplitudes and phases, the E-field vector tip traces an ellipse.

Python Visualization: Polarization States

Electric field vector tip traces for different polarization states. (a) Linear: straight line at 45°. (b) Right circular: circle. (c) Elliptical: ellipse.
Polarization states: E-field vector traces in the transverse plane
Figure 1— Electric field vector tip traces for different polarization states. (a) Linear: straight line. (b) Right circular: circle. (c) Elliptical: ellipse.

Jones Vector Formalism

The Jones vector represents a plane wave:

\[\mathbf{J} = \begin{pmatrix} E_x \\ E_y \end{pmatrix} = \begin{pmatrix} E_{0x} e^{-i\phi_x} \\ E_{0y} e^{-i\phi_y} \end{pmatrix} \tag{3}\]

Common Jones Vectors

Horizontal: \(\mathbf{J}_H = \begin{pmatrix} 1 \\ 0 \end{pmatrix}\), Vertical: \(\mathbf{J}_V = \begin{pmatrix} 0 \\ 1 \end{pmatrix}\)

Right circular: \(\mathbf{J}_{RCP} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ -i \end{pmatrix}\), Left circular: \(\mathbf{J}_{LCP} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 \\ i \end{pmatrix}\)

Intensity: \(I = |\mathbf{J}|^2 = |E_x|^2 + |E_y|^2\)

Jones Matrices and Optical Elements

Polarizers

Horizontal: \(M_H = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}\), Vertical: \(M_V = \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix}\)

At angle \(\theta\): \(M(\theta) = \begin{pmatrix} \cos^2\theta & \cos\theta\sin\theta \\ \cos\theta\sin\theta & \sin^2\theta \end{pmatrix}\) {#eq-jones-polarizer}

Retardance Plates

Quarter-wave at 45°: \(M_{QWP} = \frac{1}{\sqrt{2}} \begin{pmatrix} 1 & i \\ i & 1 \end{pmatrix}\)

Half-wave at angle \(\alpha\): \(M_{HWP}(\alpha) = \begin{pmatrix} \cos 2\alpha & \sin 2\alpha \\ \sin 2\alpha & -\cos 2\alpha \end{pmatrix}\)

Rotation: \(R(\theta) = \begin{pmatrix} \cos\theta & -\sin\theta \\ \sin\theta & \cos\theta \end{pmatrix}\)

Jones Matrix Example

Code
# Example: H-polarized light → QWP at 45° → Right circular

import numpy as np

# Matrices
M_pol = np.array([[1, 0], [0, 0]], dtype=complex)
M_qwp = (1/np.sqrt(2)) * np.array([[1, 1j], [1j, 1]], dtype=complex)

# Initial: horizontal linear
J = np.array([1, 0], dtype=complex)
print("Initial: H-polarized")
print(f"  Intensity: {np.sum(np.abs(J)**2):.2f}")

# After polarizer
J = M_pol @ J
print("After H-polarizer: intensity =", np.sum(np.abs(J)**2).real)

# After QWP at 45°
J = M_qwp @ J
print("After QWP at 45°: Jones vector =", J)
print("Result: Right Circular Polarization!")
Initial: H-polarized
  Intensity: 1.00
After H-polarizer: intensity = 1.0
After QWP at 45°: Jones vector = [0.70710678+0.j         0.        +0.70710678j]
Result: Right Circular Polarization!

Malus’s Law

When linearly polarized light passes through a second polarizer at angle \(\theta\):

\[I_{\text{transmitted}} = I_0 \cos^2(\theta) \tag{4}\]

Malus’s Law visualization. (a) Cartesian plot showing transmitted intensity vs analyzer angle. (b) Polar plot showing the \(\cos^2(\theta)\) dependence.
Figure 2— Malus’s Law: Intensity through two polarizers as a function of angle between them.

Stokes Parameters

Stokes parameters describe any polarization state (fully or partially polarized):

\[S_0 = \langle E_x E_x^* \rangle + \langle E_y E_y^* \rangle \tag{5}\] \[S_1 = \langle E_x E_x^* \rangle - \langle E_y E_y^* \rangle\] \[S_2 = \langle E_x E_y^* + E_y E_x^* \rangle\] \[S_3 = i\langle E_x E_y^* - E_y E_x^* \rangle\]

Examples

Unpolarized: \(\mathbf{S} = I_0 \begin{pmatrix} 1 \\ 0 \\ 0 \\ 0 \end{pmatrix}\)

H-linear: \(\mathbf{S}_H = I_0 \begin{pmatrix} 1 \\ 1 \\ 0 \\ 0 \end{pmatrix}\)

Right circular: \(\mathbf{S}_{RCP} = I_0 \begin{pmatrix} 1 \\ 0 \\ 0 \\ 1 \end{pmatrix}\)

Degree of polarization: \(P = \frac{\sqrt{S_1^2 + S_2^2 + S_3^2}}{S_0}\)

The Poincaré Sphere

The Poincaré sphere geometrically represents all fully polarized states using normalized Stokes parameters:

\[x = S_1/S_0, \quad y = S_2/S_0, \quad z = S_3/S_0\]

Each point on the unit sphere represents a unique polarization state: - North pole: Right circular - South pole: Left circular
- Equator: All linear polarizations

Poincaré sphere: geometric representation of all fully polarized states. North pole: right circular (RCP), South pole: left circular (LCP), Equator: all linear polarizations (H, V, ±45°).
Poincaré sphere: geometric representation of polarization states
Figure 3— The Poincaré sphere represents all fully polarized states geometrically.

Müller Matrices

The Stokes vector transforms through optical elements via Müller matrices (\(4 \times 4\)):

Horizontal polarizer: \[M_H = \frac{1}{2} \begin{pmatrix} 1 & 1 & 0 & 0 \\ 1 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix}\]

Retarder (retardance \(\Gamma\)): \[M_{\text{ret}} = \begin{pmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & \cos\Gamma & \sin\Gamma \\ 0 & 0 & -\sin\Gamma & \cos\Gamma \end{pmatrix}\]

Müller matrices work for any polarization state.

Practical Polarizing Elements

Linear Polarizers

Wire-grid: Parallel metal wires transmit perpendicular polarization, reflect parallel.

Dichroic (Polaroid): Absorbing aligned molecules. Affordable, used in displays.

Glan-Thompson: Calcite crystals at Brewster angle. Extinction ratio \(> 10^{-8}\), highest quality.

Brewster angle: p-polarized light has zero reflection at \(\theta_B = \arctan(n)\).

Retardance Plates

Materials: Quartz, calcite, mica (birefringent crystals)

Retardance: \(\Gamma = \frac{2\pi \Delta n d}{\lambda}\) where \(\Delta n = |n_e - n_o|\)

Quarter-wave (QWP): \(\Gamma = \pi/2\) converts linear ↔︎ circular at 45°

Half-wave (HWP): \(\Gamma = \pi\) rotates linear polarization

Achromatic: Multiple plates cancel wavelength dependence

Optical Rotation

Optical rotation (chiral materials): \(\theta = \alpha \ell\) (rotatory power \(\alpha\))

Faraday rotation (magneto-optic): \(\theta = V B \ell\) (non-reciprocal, used in isolators)

Summary

  1. Polarization is the oscillation direction of the electric field perpendicular to propagation.

  2. Three states: linear (straight oscillation), circular (rotation), elliptical (combination).

  3. Jones vectors/matrices for fully polarized light; Stokes parameters/Müller matrices for any state.

  4. Poincaré sphere: Geometric representation of all polarization states.

  5. Optical elements: Polarizers (select axis), retarders (introduce phase), rotators (rotate axis).

  6. Key applications: Fiber optics, LCD displays, stress analysis, quantum optics.

The elegant mathematical structure—from Maxwell’s equations through Jones calculus to Poincaré geometry—demonstrates the remarkable unity underlying optical phenomena.

Experimental Connections

Polarization is one of the most accessible optical phenomena for hands-on experiments. The following activities connect the theory above to real laboratory practice:

Polarizer and analyzer (Malus’s law) Set up two linear polarizers on an optical rail with a white-light source or laser. Rotate the analyzer in 10° steps and measure the transmitted intensity with a photodetector. Plot \(I(\theta)\) and fit to \(I_0 \cos^2\theta\). Deviations from perfect extinction reveal the quality of your polarizers and any residual birefringence in optical components.

Quarter-wave plate: linear → circular Insert a quarter-wave plate between crossed polarizers. Rotate it until the fast axis is at 45° to the first polarizer — you should recover maximum transmission through the analyzer (since circular polarization has equal projections on any linear axis). Verify by rotating the analyzer: intensity should be nearly constant.

Stress birefringence (photoelasticity) Place a transparent plastic ruler or protractor between crossed polarizers and illuminate with white light. The colorful fringe patterns reveal internal stress from the manufacturing process. This directly visualizes how mechanical stress creates birefringence. Quantitative analysis connects stress to the retardance via the stress-optic coefficient.

Optical activity in sugar solutions Fill a glass tube with sugar solution of known concentration. Place between crossed polarizers and observe the rotation of the polarization plane. By varying concentration and tube length, verify the linear dependence \(\theta = [\alpha] \cdot c \cdot L\). This is the basis of saccharimetry (sugar concentration measurement in the food industry).

Brewster angle measurement Reflect a laser beam off a glass surface (microscope slide) and vary the angle of incidence. At Brewster’s angle \(\theta_B = \arctan(n)\), the reflected beam is perfectly linearly polarized. Verify with a polarizer. This gives an elegant measurement of the refractive index.

Further Reading

The following references are linked to the central Resources & Recommended Reading page:

  • Saleh & Teich, Ch. 6 — Polarization Optics. Systematic development of Jones and Müller calculus with many worked examples.
  • Hecht, Ch. 8 — Polarization. Excellent physical intuition, many photographs of real polarization effects.
  • Goldstein (2011)Polarized Light. Comprehensive reference for Stokes parameters, Müller matrices, and the Poincaré sphere.
  • Collett (2005)Field Guide to Polarization. Compact summary of all key formulae.
  • Born & Wolf, Ch. 1, 14–15 — Rigorous EM treatment of polarization and crystal optics.

Lecture 3: Polarization of Light — Introduction to Photonics