Introduction to Stokes Parameters and Jones Matrices for Polarization States

In the study of wave optics, the polarization state of light describes the specific orientation and trajectory of the electric field vector as it propagates through space. While the fundamental concept—the vibration of light in a plane perpendicular to the direction of propagation—is straightforward, the practical application of this concept in engineering and physics requires rigorous mathematical frameworks.

In ideal theoretical models, we often deal with perfectly coherent, fully polarized light. However, real-world optical systems frequently encounter complex light states, including partially polarized light and unpolarized light. To accurately model, manipulate, and measure these states, two primary mathematical toolsets have been developed: Jones Calculus and the Stokes-Mueller Formalism. These two frameworks serve different purposes, catering to different levels of coherence and different experimental requirements.


Jones Calculus: The Coherent Approach

Developed by R. Clark Jones in 1941, Jones calculus is a powerful method for describing the polarization of fully polarized (coherent) light. It treats the electric field as a complex-valued vector, allowing for a precise description of both amplitude and phase.

The Jones Vector

For a monochromatic, fully polarized light wave, the electric field components in two orthogonal directions (typically $x$ and $y$) can be represented as a column vector known as the Jones vector:

$$
\mathbf{E} = \begin{pmatrix} E_x \ E_y \end{pmatrix}
$$

Here, $E_x$ and $E_y$ are complex numbers. The relative magnitudes of these components determine the degree of linear or circular polarization, while the relative phase between them determines the orientation of the polarization ellipse. This vector-based approach is exceptionally efficient for modeling laser beams and other highly coherent light sources.

The Jones Matrix and System Cascading

When light interacts with a non-depolarizing optical element—such as a linear polarizer, a quarter-wave plate, or a retarder—its polarization state is transformed. This transformation is mathematically represented by a $2 \times 2$ complex matrix called the Jones matrix ($\mathbf{J}$):

$$
\mathbf{E}{out} = \mathbf{J} \mathbf{E}{in}
$$

One of the most significant advantages of Jones calculus is its ability to handle cascaded optical systems. If a light beam passes through a series of optical components, the total transformation of the system can be found by simply multiplying the individual Jones matrices in the order they are encountered:

$$
\mathbf{J}_{total} = \mathbf{J}n \cdot \mathbf{J}{n-1} \dots \mathbf{J}_1
$$

This property makes Jones calculus the standard tool for designing precision laser optics, fiber-optic polarization controllers, and coherent communication systems.


Stokes Parameters and Mueller Matrices: The Intensity-Based Approach

While Jones calculus is elegant, it has a fundamental limitation: it cannot describe unpolarized light or the phenomenon of depolarization. In many practical scenarios—such as sunlight scattering in the atmosphere or light reflected from a rough surface—the light is only partially polarized. To address this, the Stokes-Mueller formalism was developed.

The Stokes Parameters

The Stokes parameters provide a way to describe the polarization state of a light beam using four real-valued components $(S_0, S_1, S_2, S_3)$. Unlike Jones vectors, which describe the field itself, Stokes parameters describe the time-averaged intensity of the light:

  • $S_0$: The total intensity of the light beam.
  • $S_1$: The difference between the intensity of horizontally polarized light and vertically polarized light.
  • $S_2$: The difference between the intensity of $+45^\circ$ polarized light and $-45^\circ$ polarized light.
  • $S_3$: The difference between the intensity of right-circularly polarized light and left-circularly polarized light.

Because these parameters are based on intensity, they can be directly measured in a laboratory using standard equipment like photodiodes and rotating polarizers. Furthermore, this framework seamlessly covers all states: for unpolarized light, $S_1 = S_2 = S_3 = 0$, while for fully polarized light, the relationship $S_0^2 = S_1^2 + S_2^2 + S_3^2$ holds.

The Mueller Matrix

To describe how an optical element affects these intensity-based parameters, we use a $4 \times 4$ real matrix known as the Mueller matrix ($\mathbf{M}$):

$$
\mathbf{S}{out} = \mathbf{M} \mathbf{S}{in}
$$

The Mueller matrix is more robust than the Jones matrix because it can account for depolarizing effects—processes where a coherent beam becomes partially unpolarized due to scattering or material imperfections. This makes it indispensable for analyzing complex, non-coherent, or stochastic optical environments.


Comparative Analysis: Choosing the Right Tool

The choice between Jones and Stokes-Mueller formalisms depends entirely on the nature of the light source and the physical phenomena being studied.

Feature Jones Calculus Stokes/Mueller Formalism
Mathematical Basis $2 \times 2$ Complex Matrices $4 \times 4$ Real Matrices
Light Type Fully polarized (Coherent) All states (Coherent, Partial, Unpolarized)
Information Content Includes absolute phase and amplitude Includes intensity and relative phase
Depolarization Cannot model depolarization Can model depolarization and scattering
Complexity Lower (2D complex space) Higher (4D real space)
Primary Use Case Laser systems, coherent optics, fiber sensing Remote sensing, imaging, atmospheric optics

Engineering Applications and Practical Implementation

Precision Laser Design (Jones Domain)

In a typical laser setup, an engineer might pass a horizontally polarized beam through a quarter-wave plate to create circular polarization. By representing the initial state as $\mathbf{E}_{in} = [1, 0]^T$ and applying the known Jones matrix for a quarter-wave plate, the output state can be predicted with near-perfect accuracy. This is critical in applications like optical trapping or interferometry, where phase precision is paramount.

Remote Sensing and Material Characterization (Stokes Domain)

In contrast, consider a satellite performing polarimetric remote sensing to detect oil spills on the ocean or to characterize forest canopies. The light reflected from these surfaces is highly scattered and partially polarized. A Jones matrix would fail here because the phase relationship between different parts of the beam is lost. Instead, engineers use Mueller matrices to characterize the "polarimetric signature" of the surface, allowing them to extract information about surface roughness, moisture content, and material composition.

Conclusion

In summary, Jones calculus and the Stokes-Mueller formalism are complementary pillars of polarization optics. Jones calculus offers a high-precision, computationally efficient method for analyzing coherent light and phase-sensitive systems. Meanwhile, the Stokes-Mueller formalism provides a comprehensive, experimentally accessible framework for describing the full spectrum of light, from perfectly ordered lasers to the chaotic, unpolarized light of the natural world. A deep understanding of both is essential for any advanced work in optical design, sensing, or imaging.