Jones calculus

In optics, polarized light can be described using the Jones calculus, discovered by R. C. Jones in 1941. Polarized light is represented by a Jones vector, and linear optical elements are represented by Jones matrices. When light crosses an optical element the resulting polarization of the emerging light is found by taking the product of the Jones matrix of the optical element and the Jones vector of the incident light. Note that Jones calculus is only applicable to light that is already fully polarized. Light which is randomly polarized, partially polarized, or incoherent must be treated using Mueller calculus.

The Jones vector

The Jones vector describes the polarization of light in free space or another homogeneous isotropic non-attenuating medium, where the light can be properly described as transverse waves. Suppose that a monochromatic plane wave of light is travelling in the positive z-direction, with angular frequency ω and wavevector k = (0,0,k), where the wavenumber k = ω/c. Then the electric and magnetic fields E and H are orthogonal to k at each point; they both lie in the plane "transverse" to the direction of motion. Furthermore, H is determined from E by 90-degree rotation and a fixed multiplier depending on the wave impedance of the medium. So the polarization of the light can be determined by studying E. The complex amplitude of E is written

${\begin{pmatrix}E_{x}(t)\\E_{y}(t)\\0\end{pmatrix}}={\begin{pmatrix}E_{{0x}}e^{{i(kz-\omega t+\phi _{x})}}\\E_{{0y}}e^{{i(kz-\omega t+\phi _{y})}}\\0\end{pmatrix}}={\begin{pmatrix}E_{{0x}}e^{{i\phi _{x}}}\\E_{{0y}}e^{{i\phi _{y}}}\\0\end{pmatrix}}e^{{i(kz-\omega t)}}$.

Note that the physical E field is the real part of this vector; the complex multiplier serves up the phase information. Here $i$ is the imaginary unit with $i^{2}=-1$.

The Jones vector is then

${\begin{pmatrix}E_{{0x}}e^{{i\phi _{x}}}\\E_{{0y}}e^{{i\phi _{y}}}\end{pmatrix}}\;.$

Thus, the Jones vector represents the amplitude and phase of the electric field in the x and y directions.

The sum of the squares of the absolute values of the two components of Jones vectors is proportional to the intensity of light. It is common to normalize it to 1 at the starting point of calculation for simplification. It is also common to constrain the first component of the Jones vectors to be a real number. This discards the overall phase information that would be needed for calculation of interference with other beams.

Note that all Jones vectors and matrices on this article employ the convention that the phase of the light wave is given by $\phi =kz-\omega t$, a convention used by Hecht. Under this convention, increase in $\phi _{x}$ (or $\phi _{y}$) indicates retardation (delay) in phase, while decrease indicates advance in phase. For example, a Jones vectors component of $i$ ($=e^{{i\pi /2}}$) indicates retardation by $\pi /2$ (or 90 degree) compared to 1 ($=e^{{0}}$). Collett uses the opposite definition for the phase ($\phi =\omega t-kz$). The reader should be wary of the choice of convention when consulting references on the Jones calculus.

The following table gives the 6 common examples of normalized Jones vectors.

Polarization	Corresponding Jones vector	Typical ket Notation
Linear polarized in the x-direction Typically called 'Horizontal'	${\begin{pmatrix}1\\0\end{pmatrix}}$	$|H\rangle$
Linear polarized in the y-direction Typically called 'Vertical'	${\begin{pmatrix}0\\1\end{pmatrix}}$	$|V\rangle$
Linear polarized at 45° from the x-axis Typically called 'Diagonal' L+45	${\frac {1}{{\sqrt 2}}}{\begin{pmatrix}1\\1\end{pmatrix}}$	$|D\rangle ={\frac {1}{{\sqrt 2}}}(|H\rangle +|V\rangle )$
Linear polarized at −45° from the x-axis Typically called 'Anti-Diagonal' L-45	${\frac {1}{{\sqrt 2}}}{\begin{pmatrix}1\\-1\end{pmatrix}}$	$|A\rangle ={\frac {1}{{\sqrt 2}}}(|H\rangle -|V\rangle )$
Right Hand Circular Polarized Typically called RCP or RHCP	${\frac {1}{{\sqrt 2}}}{\begin{pmatrix}1\\-i\end{pmatrix}}$	$|R\rangle ={\frac {1}{{\sqrt 2}}}(|H\rangle -i|V\rangle )$
Left Hand Circular Polarized Typically called LCP or LHCP	${\frac {1}{{\sqrt 2}}}{\begin{pmatrix}1\\+i\end{pmatrix}}$	$|L\rangle ={\frac {1}{{\sqrt 2}}}(|H\rangle +i|V\rangle )$

A general vector that points to any place on the surface is written as a ket $|\psi \rangle$. When employing the Poincaré sphere (also known as the Bloch sphere), the basis kets ($|0\rangle$ and $|1\rangle$) must be assigned to opposing (antipodal) pairs of the kets listed above. For example, one might assign $|0\rangle$ = $|H\rangle$ and $|1\rangle$ = $|V\rangle$. These assignments are arbitrary. Opposing pairs are

• $|H\rangle$ and $|V\rangle$
• $|D\rangle$ and $|A\rangle$
• $|R\rangle$ and $|L\rangle$

The polarization of any point not equal to $|R\rangle$ or $|L\rangle$ and not on the circle that passes through $|H\rangle ,|D\rangle ,|V\rangle ,|A\rangle$ is known as elliptical polarization.

Jones matrices

The Jones matrices are operators that act on the Jones vectors defined above. These matrices are implemented by various optical elements such as lenses, beam splitters, mirrors, etc. Each matrix represents projection onto a one-dimensional complex subspace of the Jones vectors. The following table gives examples of Jones matrices for polarizers:

Optical element	Corresponding Jones matrix
Linear polarizer with axis of transmission horizontal[1]	${\begin{pmatrix}1&0\\0&0\end{pmatrix}}$
Linear polarizer with axis of transmission vertical[1]	${\begin{pmatrix}0&0\\0&1\end{pmatrix}}$
Linear polarizer with axis of transmission at ±45° with the horizontal[1]	${\frac {1}{2}}{\begin{pmatrix}1&\pm 1\\\pm 1&1\end{pmatrix}}$
Right circular polarizer[1]	${\frac {1}{2}}{\begin{pmatrix}1&-i\\i&1\end{pmatrix}}$
Left circular polarizer[1]	${\frac {1}{2}}{\begin{pmatrix}1&i\\-i&1\end{pmatrix}}$

Phase retarders

Phase retarders introduce a phase shift between the vertical and horizontal component of the field and thus change the polarization of the beam. Phase retarders are usually made out of birefringent uniaxial crystals such as calcite, MgF₂ or quartz. Uniaxial crystals have one crystal axis that is different from the other two crystal axes (i.e., n_i ≠ n_j = n_k). This unique axis is called the extraordinary axis and is also referred to as the optic axis. An optic axis can be the fast or the slow axis for the crystal depending on the crystal at hand. Light travels with a higher phase velocity through an axis that has the smallest refractive index and this axis is called the fast axis. Similarly, an axis which has the highest refractive index is called a slow axis since the phase velocity of light is the lowest along this axis. "Negative" uniaxial crystals (e.g., calcite CaCO₃, sapphire Al₂O₃) have n_e < n_o so for these crystals, the extraordinary axis (optic axis) is the fast axis, whereas for "positive" uniaxial crystals (e.g., quartz SiO₂, magnesium fluoride MgF₂, rutile TiO₂), n_e > n _o and thus the extraordinary axis (optic axis) is the slow axis.

Any phase retarder with fast axis equal to the x- or y-axis has zero off-diagonal terms and thus can be conveniently expressed as

${\begin{pmatrix}e^{{i\phi _{x}}}&0\\0&e^{{i\phi _{y}}}\end{pmatrix}}$

where $\phi _{x}$ and $\phi _{y}$ are the phase offsets of the electric fields in $x$ and $y$ directions respectively. In the phase convention $\phi =kz-\omega t$, define the relative phase between the two waves as $\epsilon =\phi _{y}-\phi _{x}$. Then a positive $\epsilon$ (i.e. $\phi _{y}$ > $\phi _{x}$) means that $E_{y}$ doesn't attain the same value as $E_x$ until a later time, i.e. $E_x$ leads $E_{y}$. Similarly, if $\epsilon <0$, then $E_{y}$ leads $E_x$. For example, if the fast axis of a quarter wave plate is horizontal, then the phase velocity along the horizontal direction is ahead of the vertical direction i.e., $E_x$ leads $E_{y}$. Thus, $\phi _{x}<\phi _{y}$ which for a quarter wave plate yields $\phi _{y}=\phi _{x}+\pi /2$.

In the opposite convention $\phi =\omega t-kz$, define the relative phase as $\epsilon =\phi _{x}-\phi _{y}$. Then $\epsilon >0$ means that $E_{y}$ doesn't attain the same value as $E_x$ until a later time, i.e. $E_{x}$ leads $E_{y}$.

Phase retarders	Corresponding Jones matrix
Quarter-wave plate with fast axis vertical[2][note 1]	$e^{{i\pi /4}}{\begin{pmatrix}1&0\\0&-i\end{pmatrix}}$
Quarter-wave plate with fast axis horizontal[2]	$e^{{i\pi /4}}{\begin{pmatrix}1&0\\0&i\end{pmatrix}}$
Quarter-wave plate with fast axis at angle $\theta$ w.r.t the horizontal axis	${\displaystyle e^{i\pi /4}{\begin{pmatrix}\cos ^{2}\theta +i\sin ^{2}\theta &(1-i)\sin \theta \cos \theta \\(1-i)\sin \theta \cos \theta &\sin ^{2}\theta +i\cos ^{2}\theta \end{pmatrix}}}$
Half-wave plate with fast axis at angle $\theta$ w.r.t the horizontal axis[3]	${\begin{pmatrix}\cos 2\theta &\sin 2\theta \\\sin 2\theta &-\cos 2\theta \end{pmatrix}}$
Arbitrary birefringent material (as phase retarder)[4]	$\begin{pmatrix} e^{i\eta/2} \cos^2\theta+e^{-i\eta/2} \sin^2\theta & (e^{i\eta/2}-e^{-i\eta/2}) e^{-i\phi} \cos\theta \sin\theta \\ (e^{i\eta/2}-e^{-i\eta/2}) e^{i\phi} \cos\theta \sin\theta & e^{i\eta/2} \sin^2\theta+e^{-i\eta/2} \cos^2\theta \end{pmatrix}$

The special expressions for the phase retarders can be obtained by taking suitable parameter values in the general expression for a birefringent material. The general expression can be obtained In the general expression:

• The relative phase retardation induced between the fast axis and the slow axis is given by $\eta = \phi_y - \phi_x$
• $\theta$ is the orientation of the fast axis with respect to the x-axis.
• $\phi$ is the circularity.

Note that for linear retarders, $\phi$ = 0 and for circular retarders, $\phi$ = ± $\pi$/2, $\theta$ = $\pi$/4. In general for elliptical retarders, $\phi$ takes on values between - $\pi$/2 and $\pi$/2.

Axially rotated elements

Assume an optical element has its optic axis perpendicular to the surface vector for the plane of incidence and is rotated about this surface vector by angle θ/2 (i.e., the principal plane, through which the optic axis passes, makes angle θ/2 with respect to the plane of polarization of the electric field of the incident TE wave). Recall that a half-wave plate rotates polarization as twice the angle between incident polarization and optic axis (principal plane). Therefore, the Jones matrix for the rotated polarization state, M(θ), is

${\displaystyle M(\theta )=R(-\theta )\,M\,R(\theta ),}$

where ${\displaystyle R(\theta )={\begin{pmatrix}\cos \theta &\sin \theta \\-\sin \theta &\cos \theta \end{pmatrix}}.}$

This agrees with the expression for a half-wave plate in the table above. These rotations are identical to beam unitary splitter transformation in optical physics given by

$R(\theta )={\begin{pmatrix}r&t'\\t&r'\end{pmatrix}}$

where the primed and unprimed coefficients represent beams incident from opposite sides of the beam splitter. The reflected and transmitted components acquire a phase θ_r and θ_t, respectively. The requirements for a valid representation of the element are ^[5]

$\theta _{{\text{t}}}-\theta _{{\text{r}}}+\theta _{{\text{t'}}}-\theta _{{\text{r'}}}=\pm \pi$

and $r^{*}t'+t^{*}r'=0.$

Both of these representations are unitary matrices fitting these requirements; and as such, are both valid.

Arbitrarily rotated elements

This would involve a three-dimensional rotation matrix. See Garam Yun for work done on this.^[6][7]

See also

• Mueller calculus
• Stokes parameters
• Polarization

Notes

1. The prefactor $e^{{i\pi /4}}$ appears only if one defines the phase delays in a symmetric fashion; that is, $\phi _{x}=-\phi _{y}=\pi /4$. This is done in Hecht^[2] but not in Fowles.^[1] In the latter reference the Jones matrices for a quarter-wave plate have no prefactor.

References

1. Fowles, G. (1989). Introduction to Modern Optics (2nd ed.). Dover. p. 35.
2. Hecht, E. (2001). Optics (4th ed.). p. 378. ISBN 0805385665.
3. Gerald, A.; Burch, J.M. (1975). Introduction to Matrix Methods in Optics (1st ed.). John Wiley & Sons. ISBN 0471296856.
4. Obtainment of the polarizing and retardation parameters of a non-depolarizing optical system from the polar decomposition of its Mueller matrix, Optik, Jose Jorge Gill and Eusebio Bernabeu,76, 67-71 (1987).
5. Am. J. Phys. 57 (1), 66 (1988).
6. Three-dimensional polarization ray-tracing calculus I: definition and diattenuation, Applied Optics, Garam Yun, Karlton Crabtree, and Russell A. Chipman,50, 2855-2865 (2011).
7. Three-dimensional polarization ray-tracing calculus II: retardance, Applied Optics, Garam Yun, Stephen C. McClain, and Russell A. Chipman,50, 2866-2874 (2011).

Further reading

• E. Collett, Field Guide to Polarization, SPIE Field Guides vol. FG05, SPIE (2005). ISBN 0-8194-5868-6.
• D. Goldstein and E. Collett, Polarized Light, 2nd ed., CRC Press (2003). ISBN 0-8247-4053-X.
• E. Hecht, Optics, 2nd ed., Addison-Wesley (1987). ISBN 0-201-11609-X.
• Frank L. Pedrotti, S.J. Leno S. Pedrotti, Introduction to Optics, 2nd ed., Prentice Hall (1993). ISBN 0-13-501545-6
• A. Gerald and J.M. Burch, Introduction to Matrix Methods in Optics,1st ed., John Wiley & Sons(1975). ISBN 0-471-29685-6
• Jones, R. Clark (1941). "A new calculus for the treatment of optical systems, I. Description and Discussion of the Calculus". Journal of the Optical Society of America. 31 (7): 488–493. doi:10.1364/JOSA.31.000488.
• Hurwitz, Henry; Jones, R. Clark (1941). "A new calculus for the treatment of optical systems, II. Proof of three general equivalence theorems". Journal of the Optical Society of America. 31 (7): 493–499. doi:10.1364/JOSA.31.000493.
• Jones, R. Clark (1941). "A new calculus for the treatment of optical systems, III The Sohncke Theory of optical activity". Journal of the Optical Society of America. 31 (7): 500–503. doi:10.1364/JOSA.31.000500.
• Jones, R. Clark (1942). "A new calculus for the treatment of optical systems, IV". Journal of the Optical Society of America. 32 (8): 486–493. doi:10.1364/JOSA.32.000486.
• Fymat, A. L. (1971). "Jones's Matrix Representation of Optical Instruments. I: Beam Splitters". Applied Optics. 10 (11): 2499–2505. Bibcode:1971ApOpt..10.2499F. doi:10.1364/AO.10.002499. PMID 20111363.
• Fymat, A. L. (1971). "Jones's Matrix Representation of Optical Instruments. 2: Fourier Interferometers (Spectrometers and Spectropolarimeters)". Applied Optics. 10 (12): 2711–2716. Bibcode:1971ApOpt..10.2711F. doi:10.1364/AO.10.002711.
• Fymat, A. L. (1972). "Polarization Effects in Fourier Spectroscopy. I: Coherency Matrix Representation". Applied Optics. 11 (1): 160–173. Bibcode:1972ApOpt..11..160F. doi:10.1364/AO.11.000160. PMID 20111472.
• Gill, Jose Jorge; Bernabeu, Eusebio (1987). "Obtainment of the polarizing and retardation parameters of a non-depolarizing optical system from the polar decomposition of its Mueller matrix,". Optik. 76: 67–71.
• Brosseau, Christian; Givens, Clark R.; Kostinksi, Alexander B. (1993). "Generalized trace condition on the Mueller-Jones polarization matrix". Journal of the Optical Society of America A. 10 (10): 2248–2251. Bibcode:1993JOSAA..10.2248B. doi:10.1364/JOSAA.10.002248.
• McGuire, James P.; Chipman, Russel A. (1994). "Polarization aberrations. 1. Rotationally symmetric optical systems". Applied Optics. 33 (22): 5080–5100. doi:10.1364/AO.33.005080. PMID 20935891.
• Pistoni, Natale C. (1995). "Simplified approach to the Jones calculus in retracing optical circuits". Applied Optics. 34 (34): 7870–7876. Bibcode:1995ApOpt..34.7870P. doi:10.1364/AO.34.007870. PMID 21068881.
• Moreno, Ignacio; Yzuel, Maria J.; Campos, Juan; Vargas, Asticio (2004). "Jones matrix treatment for polarization Fourier optics". Journal of Modern Optics. 51 (14): 2031–2038. Bibcode:2000JMOp...51.2031M. doi:10.1080/09500340408232511.
• Moreno, Ivan (2004). "Jones matrix for image-rotation prisms". Applied Optics. 43 (17): 3373–3381. Bibcode:2004ApOpt..43.3373M. doi:10.1364/AO.43.003373. PMID 15219016.

External links

• Jones Calculus written by E. Collett on Optipedia