advanced30 min studyLesson 217 in Physics

Stokes parameters and partial polarization

Replace a single coherent polarization vector by a positive coherency matrix and connect its Stokes coordinates to measurable analyzer intensities.

On this page

What you’ll understand

  • Construct the coherency matrix from averaged field products.
  • Calculate Stokes parameters and degree of polarization with a declared circular sign.
  • Distinguish an incoherent polarization mixture from coherent amplitude addition.

A light beam may fluctuate among different polarization states during a detector exposure. Assigning one Jones vector to its mean electric field fails, because that mean can vanish while the measured intensity is large. The coherency matrix instead averages products of complex field components. Its four real degrees of freedom can be written as Stokes parameters, which connect polarization statistics to a practical set of analyzer measurements.

Average quadratic quantities before describing the state

Use a common transverse basis and the same negative-sign time factor as in Jones calculus. Define Jab=⟨EaEb∗⟩J_{ab}=\langle E_aE_b^*\rangle, including the common impedance conversion so the entries have irradiance units. Diagonal entries are nonnegative component intensities; off-diagonal entries retain their mutual coherence. Hermitian symmetry follows directly from complex conjugation. Positivity follows because every projected intensity is the average squared magnitude of a projected field.

For one fixed coherent polarization, JJ is an outer product and has rank one. For an incoherent mixture, add the coherency matrices weighted by their intensities or probabilities; do not add Jones amplitudes as though relative phase were stable. An equal incoherent mixture of orthogonal linear polarizations has zero off-diagonal entries. A coherent equal-amplitude sum of those same directions is instead linearly polarized at forty-five degrees and has a nonzero off-diagonal entry.

J=(⟨∣Ex∣2⟩⟨ExEy∗⟩⟨EyEx∗⟩⟨∣Ey∣2⟩)=12(S0+S1S2−iS3S2+iS3S0−S1)J=\begin{pmatrix}\langle |E_x|^2\rangle&\langle E_xE_y^*\rangle\\\langle E_yE_x^*\rangle&\langle |E_y|^2\rangle\end{pmatrix}=\frac12\begin{pmatrix}S_0+S_1&S_2-iS_3\\S_2+iS_3&S_0-S_1\end{pmatrix}
All J entries and Stokes parameters use the same irradiance units after a common impedance factor. S1=Jxx−Jyy, S2=2 Re Jxy, and S3=−2 Im Jxy. With exp(−i omega t), the normalized vector (1,i)/sqrt(2) has positive S3.↔ Scroll wide equations horizontally.

The interior of the sphere is physically meaningful

The total irradiance is S0S_0, while the length of the remaining Stokes vector measures its polarized part. Positivity bounds that length by S0S_0. A nonzero beam on the boundary is fully polarized, and a beam at the centre has no preferred polarization at this order of description. Interior points represent partial polarization. Normalizing the three polarization coordinates by their own length would incorrectly move every partially polarized state onto the fully polarized boundary.

The two eigenvalues of the coherency matrix are half the total intensity plus or minus half the polarization-vector length. They describe an incoherent decomposition into two orthogonal polarization states. An unpolarized contribution can also be separated from a fully polarized contribution, but neither decomposition means that individual photons carry a classical hidden polarization label. These are equivalent descriptions of measured second-order statistics, not unique microscopic histories.

P=S12+S22+S32S0,0≤P≤1,j±=S0±S12+S22+S322\mathcal P=\frac{\sqrt{S_1^2+S_2^2+S_3^2}}{S_0},\quad 0\leq\mathcal P\leq1,\qquad j_\pm=\frac{S_0\pm\sqrt{S_1^2+S_2^2+S_3^2}}{2}
The degree of polarization P is defined for S0>0 and is dimensionless; j± have irradiance units. The inequalities assume a physical positive semidefinite two-transverse-component coherency matrix. Polarization degree is undefined for a beam of exactly zero intensity.↔ Scroll wide equations horizontally.

Analyze intensities and keep coherence questions separate

An ideal linear analyzer at angle α\alpha transmits half the sum of total intensity and the appropriate double-angle projection of the linear Stokes components. Measuring orthogonal pairs at zero and ninety degrees, then forty-five and one hundred thirty-five degrees, determines S0,S1,S2S_0,S_1,S_2. Circular analyzers determine S3S_3 after their transmitted circular basis states have been defined. Circular-sign conventions vary, so a table without a basis definition is incomplete.

These Stokes parameters describe polarization at one location and within a chosen temporal or spectral averaging window. They do not determine mutual coherence between separated points or separate beams. A fully polarized source can still have limited temporal coherence; a partially polarized field can contain a coherent polarized component. The four-parameter transverse description also presumes one propagation direction or a suitable paraxial collection mode. Strongly nonparaxial, angularly unresolved fields may require a larger correlation description.

Read a coherency matrix

In irradiance units of W/m2\mathrm{W/m^2}, a beam has coherency matrix J=(31−i1+i2)J=\begin{pmatrix}3&1-i\\1+i&2\end{pmatrix}. Find its four Stokes parameters, degree of polarization, and eigenvalues. Check that it is physically allowed.

  1. The diagonal sum gives S0=3+2=5S_0=3+2=5 and the difference gives S1=3−2=1S_1=3-2=1. The off-diagonal entry has real part one and imaginary part minus one, so S2=2S_2=2 and S3=2S_3=2.

  2. The polarization-vector length is 12+22+22=3 W/m2\sqrt{1^2+2^2+2^2}=3\,\mathrm{W/m^2}. Divide by total intensity to obtain degree P=3/5=0.60\mathcal P=3/5=0.60, which is less than unity.

  3. The eigenvalues are (5+3)/2=4(5+3)/2=4 and (5−3)/2=1 W/m2(5-3)/2=1\,\mathrm{W/m^2}. Both are positive, confirming a physical matrix; equivalently its determinant is 6−∣1−i∣2=4>06-|1-i|^2=4>0.

The result

The Stokes vector is (5,1,2,2) W/m², its degree of polarization is 0.60, and its eigenvalues are 4 and 1 W/m².

Explore the model

Interactive lab

Interactive lab · quantum

Only the fully polarized boundary

Pause the instrument. For each run, Reset and set BOTH Mixing angle θ and Relative phase φ. Compare θ=0°, φ=0° with θ=90°, φ=90°, predicting a move from positive Bloch z to positive Bloch y. Identify Bloch coordinates with normalized (S2,S3,S1). Every displayed state is pure; there is no mixture or degree-of-polarization control.

Only the fully polarized boundaryA normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Current readings: Probability |0⟩: 71.1%; Probability |1⟩: 28.9%; Bloch (x, y, z): 0.78, 0.45, 0.42. All parameters are adjustable with the labeled controls below.xyz|0⟩|1⟩Basis probabilitiesP(0)71.1%P(1)28.9%
Bloch vector / P(0)P(1)Phase latitude
Probability |0⟩
71.1%
Probability |1⟩
28.9%
Bloch (x, y, z)
0.78, 0.45, 0.42
t = 0 s

A normalized two-level pure state, not an interacting field simulation. Playback adds one radian per second to the relative phase; basis probabilities stay fixed. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThe lab illustrates points on the fully polarized sphere only. It cannot display the interior coherency matrix in the worked example or infer temporal coherence.

Check your understanding

An equal incoherent mixture of x and y polarized beams has the same total intensity as their coherent in-phase equal-amplitude superposition. Which distinction is correct?

Try it yourself

A beam is an incoherent mixture of 6 W/m26\,\mathrm{W/m^2} in x polarization and 2 W/m22\,\mathrm{W/m^2} in y polarization. Find its Stokes vector, degree of polarization, and intensity transmitted by an ideal 45∘45^\circ linear analyzer.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Add component coherency matrices, with no stable cross term. The analyzer transmits half of each component at forty-five degrees.

Compare with the solution

The matrix is diagonal with entries six and two, so (S0,S1,S2,S3)=(8,4,0,0) W/m2(S_0,S_1,S_2,S_3)=(8,4,0,0)\,\mathrm{W/m^2}. The polarization degree is 4/8=0.504/8=0.50. A forty-five-degree analyzer passes 6/2+2/2=4 W/m26/2+2/2=4\,\mathrm{W/m^2}. In the general analyzer formula, cos⁡2α=0\cos2\alpha=0 and sin⁡2α=1\sin2\alpha=1, giving (S0+S2)/2=4(S_0+S_2)/2=4 as well. Treating the fields as coherent in-phase amplitudes would introduce a positive cross term and give a different answer; that operation contradicts the specified incoherent mixture.

Take it with you

  • Coherency matrices average field products and remain positive semidefinite.
  • Partial polarization occupies the interior of the normalized Stokes sphere.
  • Local polarization statistics do not specify spatial or temporal mutual coherence.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.