Stokes parameters

The Stokes parameters are a set of values that describe the polarization state of electromagnetic radiation. They were defined by George Gabriel Stokes in 1852,[1][2] as a mathematically convenient alternative to the more common description of incoherent or partially polarized radiation in terms of its total intensity (I), (fractional) degree of polarization (p), and the shape parameters of the polarization ellipse. The effect of an optical system on the polarization of light can be determined by constructing the Stokes vector for the input light and applying Mueller calculus, to obtain the Stokes vector of the light leaving the system. The original Stokes paper was discovered independently by Francis Perrin in 1942[3] and by Subrahamanyan Chandrasekhar in 1947[4][5], who named it as the Stokes parameters.

Definitions

Polarisation ellipse, showing the relationship to the Poincaré sphere parameters ψ and χ.

The Poincaré sphere is the parametrisation of the last three Stokes' parameters in spherical coordinates.

Depiction of the polarization states on Poincaré sphere

The relationship of the Stokes parameters S₀, S₁, S₂, S₃ to intensity and polarization ellipse parameters is shown in the equations below and the figure at right.

${\displaystyle {\begin{aligned}S_{0}&=I\\S_{1}&=Ip\cos 2\psi \cos 2\chi \\S_{2}&=Ip\sin 2\psi \cos 2\chi \\S_{3}&=Ip\sin 2\chi \end{aligned}}}$

Here ${\displaystyle Ip}$, ${\displaystyle 2\psi }$ and ${\displaystyle 2\chi }$ are the spherical coordinates of the three-dimensional vector of cartesian coordinates ${\displaystyle (S_{1},S_{2},S_{3})}$. ${\displaystyle I}$ is the total intensity of the beam, and ${\displaystyle p}$ is the degree of polarization, constrained by ${\displaystyle 0\leq p\leq 1}$. The factor of two before ${\displaystyle \psi }$ represents the fact that any polarization ellipse is indistinguishable from one rotated by 180°, while the factor of two before ${\displaystyle \chi }$ indicates that an ellipse is indistinguishable from one with the semi-axis lengths swapped accompanied by a 90° rotation. The phase information of the polarized light is not recorded in the Stokes parameters. The four Stokes parameters are sometimes denoted I, Q, U and V, respectively.

If given the Stokes parameters one can solve for the spherical coordinates with the following equations:

${\displaystyle {\begin{aligned}I&=S_{0}\\p&={\frac {\sqrt {S_{1}^{2}+S_{2}^{2}+S_{3}^{2}}}{S_{0}}}\\2\psi &=\mathrm {arctan} {\frac {S_{2}}{S_{1}}}\\2\chi &=\mathrm {arctan} {\frac {S_{3}}{\sqrt {S_{1}^{2}+S_{2}^{2}}}}\\\end{aligned}}}$

Stokes vectors

The Stokes parameters are often combined into a vector, known as the Stokes vector:

${\displaystyle {\vec {S}}\ ={\begin{pmatrix}S_{0}\\S_{1}\\S_{2}\\S_{3}\end{pmatrix}}={\begin{pmatrix}I\\Q\\U\\V\end{pmatrix}}}$

The Stokes vector spans the space of unpolarized, partially polarized, and fully polarized light. For comparison, the Jones vector only spans the space of fully polarized light, but is more useful for problems involving coherent light. The four Stokes parameters are not a preferred coordinate system of the space, but rather were chosen because they can be easily measured or calculated.

Note that there is an ambiguous sign for the ${\displaystyle V}$ component depending on the physical convention used. In practice, there are two separate conventions used, either defining the Stokes parameters when looking down the beam towards the source (opposite the direction of light propagation) or looking down the beam away from the source (coincident with the direction of light propagation). These two conventions result in different signs for ${\displaystyle V}$, and a convention must be chosen and adhered to.

Examples

Below are shown some Stokes vectors for common states of polarization of light.

${\displaystyle {\begin{pmatrix}1\\1\\0\\0\end{pmatrix}}}$	Linearly polarized (horizontal)
${\displaystyle {\begin{pmatrix}1\\-1\\0\\0\end{pmatrix}}}$	Linearly polarized (vertical)
${\displaystyle {\begin{pmatrix}1\\0\\1\\0\end{pmatrix}}}$	Linearly polarized (+45°)
${\displaystyle {\begin{pmatrix}1\\0\\-1\\0\end{pmatrix}}}$	Linearly polarized (−45°)
${\displaystyle {\begin{pmatrix}1\\0\\0\\1\end{pmatrix}}}$	Right-hand circularly polarized
${\displaystyle {\begin{pmatrix}1\\0\\0\\-1\end{pmatrix}}}$	Left-hand circularly polarized
${\displaystyle {\begin{pmatrix}1\\0\\0\\0\end{pmatrix}}}$	Unpolarized

Alternate explanation

A monochromatic plane wave is specified by its propagation vector, ${\displaystyle {\vec {k}}}$, and the complex amplitudes of the electric field, ${\displaystyle E_{1}}$ and ${\displaystyle E_{2}}$, in a basis ${\displaystyle ({\hat {\epsilon }}_{1},{\hat {\epsilon }}_{2})}$. The pair ${\displaystyle (E_{1},E_{2})}$ is called a Jones vector. Alternatively, one may specify the propagation vector, the phase, ${\displaystyle \phi }$, and the polarization state, ${\displaystyle \Psi }$, where ${\displaystyle \Psi }$ is the curve traced out by the electric field as a function of time in a fixed plane. The most familiar polarization states are linear and circular, which are degenerate cases of the most general state, an ellipse.

One way to describe polarization is by giving the semi-major and semi-minor axes of the polarization ellipse, its orientation, and the direction of rotation (See the above figure). The Stokes parameters ${\displaystyle I}$, ${\displaystyle Q}$, ${\displaystyle U}$, and ${\displaystyle V}$, provide an alternative description of the polarization state which is experimentally convenient because each parameter corresponds to a sum or difference of measurable intensities. The next figure shows examples of the Stokes parameters in degenerate states.

Definitions

The Stokes parameters are defined by

${\displaystyle {\begin{aligned}I&\equiv \langle E_{x}^{2}\rangle +\langle E_{y}^{2}\rangle \\&=\langle E_{a}^{2}\rangle +\langle E_{b}^{2}\rangle \\&=\langle E_{r}^{2}\rangle +\langle E_{l}^{2}\rangle ,\\Q&\equiv \langle E_{x}^{2}\rangle -\langle E_{y}^{2}\rangle ,\\U&\equiv \langle E_{a}^{2}\rangle -\langle E_{b}^{2}\rangle ,\\V&\equiv \langle E_{r}^{2}\rangle -\langle E_{l}^{2}\rangle .\end{aligned}}}$

where the subscripts refer to three different bases of the space of Jones vectors: the standard Cartesian basis (${\displaystyle {\hat {x}},{\hat {y}}}$), a Cartesian basis rotated by 45° (${\displaystyle {\hat {a}},{\hat {b}}}$), and a circular basis (${\displaystyle {\hat {l}},{\hat {r}}}$). The circular basis is defined so that ${\displaystyle {\hat {l}}=({\hat {x}}+i{\hat {y}})/{\sqrt {2}}}$.

The symbols ⟨⋅⟩ represent expectation values. The light can be viewed as a random variable taking values in the space C² of Jones vectors ${\displaystyle (E_{1},E_{2})}$. Any given measurement yields a specific wave (with a specific phase, polarization ellipse, and magnitude), but it keeps flickering and wobbling between different outcomes. The expectation values are various averages of these outcomes. Intense, but unpolarized light will have I > 0 but Q = U = V = 0, reflecting that no polarization type predominates. A convincing waveform is depicted at the article on coherence.

The opposite would be perfectly polarized light which, in addition, has a fixed, nonvarying amplitude -- a pure sine curve. This is represented by a random variable with only a single possible value, say ${\displaystyle (E_{1},E_{2})}$. In this case one may replace the brackets by absolute value bars, obtaining a well-defined quadratic map

${\displaystyle {\begin{matrix}I\equiv |E_{x}|^{2}+|E_{y}|^{2}=|E_{a}|^{2}+|E_{b}|^{2}=|E_{r}|^{2}+|E_{l}|^{2}\\Q\equiv |E_{x}|^{2}-|E_{y}|^{2},\\U\equiv |E_{a}|^{2}-|E_{b}|^{2},\\V\equiv |E_{r}|^{2}-|E_{l}|^{2}.\end{matrix}}}$

from the Jones vectors to the corresponding Stokes vectors; more convenient forms are given below. The map takes its image in the cone defined by |I |² = |Q |² + |U |² + |V |², where the purity of the state satisfies p = 1 (see below).

The next figure shows how the signs of the Stokes parameters are determined by the helicity and the orientation of the semi-major axis of the polarization ellipse.

Representations in fixed bases

In a fixed (${\displaystyle {\hat {x}},{\hat {y}}}$) basis, the Stokes parameters when using an increasing phase convention are

${\displaystyle {\begin{aligned}I&=|E_{x}|^{2}+|E_{y}|^{2},\\Q&=|E_{x}|^{2}-|E_{y}|^{2},\\U&=2\mathrm {Re} (E_{x}E_{y}^{*}),\\V&=-2\mathrm {Im} (E_{x}E_{y}^{*}),\\\end{aligned}}}$

while for ${\displaystyle ({\hat {a}},{\hat {b}})}$, they are

${\displaystyle {\begin{aligned}I&=|E_{a}|^{2}+|E_{b}|^{2},\\Q&=-2\mathrm {Re} (E_{a}^{*}E_{b}),\\U&=|E_{a}|^{2}-|E_{b}|^{2},\\V&=2\mathrm {Im} (E_{a}^{*}E_{b}).\\\end{aligned}}}$

and for ${\displaystyle ({\hat {l}},{\hat {r}})}$, they are

${\displaystyle {\begin{aligned}I&=|E_{l}|^{2}+|E_{r}|^{2},\\Q&=2\mathrm {Re} (E_{l}^{*}E_{r}),\\U&=-2\mathrm {Im} (E_{l}^{*}E_{r}),\\V&=|E_{r}|^{2}-|E_{l}|^{2}.\\\end{aligned}}}$

Properties

For purely monochromatic coherent radiation, it follows from the above equations that

${\displaystyle Q^{2}+U^{2}+V^{2}=I^{2},}$

whereas for the whole (non-coherent) beam radiation, the Stokes parameters are defined as averaged quantities, and the previous equation becomes an inequality:[6]

${\displaystyle Q^{2}+U^{2}+V^{2}\leq I^{2}.}$

However, we can define a total polarization intensity ${\displaystyle I_{p}}$, so that

${\displaystyle Q^{2}+U^{2}+V^{2}=I_{p}^{2},}$

where ${\displaystyle I_{p}/I}$ is the total polarization fraction.

Let us define the complex intensity of linear polarization to be

${\displaystyle {\begin{aligned}L&\equiv |L|e^{i2\theta }\\&\equiv Q+iU.\\\end{aligned}}}$

Under a rotation ${\displaystyle \theta \rightarrow \theta +\theta '}$ of the polarization ellipse, it can be shown that ${\displaystyle I}$ and ${\displaystyle V}$ are invariant, but

${\displaystyle {\begin{aligned}L&\rightarrow e^{i2\theta '}L,\\Q&\rightarrow {\mbox{Re}}\left(e^{i2\theta '}L\right),\\U&\rightarrow {\mbox{Im}}\left(e^{i2\theta '}L\right).\\\end{aligned}}}$

With these properties, the Stokes parameters may be thought of as constituting three generalized intensities:

${\displaystyle {\begin{aligned}I&\geq 0,\\V&\in \mathbb {R} ,\\L&\in \mathbb {C} ,\\\end{aligned}}}$

where ${\displaystyle I}$ is the total intensity, ${\displaystyle |V|}$ is the intensity of circular polarization, and ${\displaystyle |L|}$ is the intensity of linear polarization. The total intensity of polarization is ${\displaystyle I_{p}={\sqrt {|L|^{2}+|V|^{2}}}}$, and the orientation and sense of rotation are given by

${\displaystyle {\begin{aligned}\theta &={\frac {1}{2}}\arg(L),\\h&=\operatorname {sgn}(V).\\\end{aligned}}}$

Since ${\displaystyle Q={\mbox{Re}}(L)}$ and ${\displaystyle U={\mbox{Im}}(L)}$, we have

${\displaystyle {\begin{aligned}|L|&={\sqrt {Q^{2}+U^{2}}},\\\theta &={\frac {1}{2}}\tan ^{-1}(U/Q).\\\end{aligned}}}$

Relation to the polarization ellipse

In terms of the parameters of the polarization ellipse, the Stokes parameters are

${\displaystyle {\begin{aligned}I_{p}&=A^{2}+B^{2},\\Q&=(A^{2}-B^{2})\cos(2\theta ),\\U&=(A^{2}-B^{2})\sin(2\theta ),\\V&=2ABh.\\\end{aligned}}}$

Inverting the previous equation gives

${\displaystyle {\begin{aligned}A&={\sqrt {{\frac {1}{2}}(I_{p}+|L|)}}\\B&={\sqrt {{\frac {1}{2}}(I_{p}-|L|)}}\\\theta =&{\frac {1}{2}}\arg(L)\\h&=\operatorname {sgn}(V).\\\end{aligned}}}$

Relationship to Hermitian operators and quantum mixed states

From a geometric and algebraic point of view, the Stokes parameters stand in one-to-one correspondence with the closed, convex, 4-real-dimensional cone of nonnegative Hermitian operators on the Hilbert space C². The parameter I serves as the trace of the operator, whereas the entries of the matrix of the operator are simple linear functions of the four parameters I, Q, U, V, serving as coefficients in a linear combination of the Stokes operators. The eigenvalues and eigenvectors of the operator can be calculated from the polarization ellipse parameters I, p, ψ, χ.

The Stokes parameters with I set equal to 1 (i.e. the trace 1 operators) are in one-to-one correspondence with the closed unit 3-dimensional ball of mixed states (or density operators) of the quantum space C², whose boundary is the Bloch sphere. The Jones vectors correspond to the underlying space C², that is, the (unnormalized) pure states of the same system. Note that phase information is lost when passing from a pure state |φ⟩ to the corresponding mixed state |φ⟩⟨φ|, just as it is lost when passing from a Jones vector to the corresponding Stokes vector.

See also

• Mueller calculus
• Jones calculus
• Polarization (waves)
• Rayleigh Sky Model
• Stokes operators
• Polarization mixing

Notes

1. Stokes, G. G. (1852). On the composition and resolution of streams of polarized light from different sources. Transactions of the Cambridge Philosophical Society, 9, 399.
2. S. Chandrasekhar 'Radiative Transfer, Dover Publications, New York, 1960, ISBN 0-486-60590-6, page 25
3. Perrin, F. (1942). Polarization of light scattered by isotropic opalescent media. The Journal of Chemical Physics, 10(7), 415-427.
4. "S. Chandrasekhar - Session II". Oral History Interviews. AIP. 18 May 1977.
5. Chandrasekhar, S. (1947). The transfer of radiation in stellar atmospheres. Bulletin of the American Mathematical Society, 53(7), 641-711.
6. H. C. van de Hulst Light scattering by small particles, Dover Publications, New York, 1981, ISBN 0-486-64228-3, page 42

References

• E. Collett, Field Guide to Polarization, SPIE Field Guides vol. FG05, SPIE (2005). ISBN 0-8194-5868-6.
• E. Hecht, Optics, 2nd ed., Addison-Wesley (1987). ISBN 0-201-11609-X.
• William H. McMaster (1954). "Polarization and the Stokes Parameters". Am. J. Phys. 22: 351. Bibcode:1954AmJPh..22..351M. doi:10.1119/1.1933744.
• William H. McMaster (1961). "Matrix representation of polarization". Rev. Mod. Phys. 33: 8. Bibcode:1961RvMP...33....8M. doi:10.1103/RevModPhys.33.8.

External links