Linear Wave Retarder

The intensity and polarization of light can be accurately controlled by transmitting it through various types of optical components. While it is possible to model complex changes in ray polarization by applying a customized Mueller matrix to rays at a single boundary, the same result can often be achieved by placing several simpler optical components in series. This tutorial model shows how the polarization of radiation can be manipulated using a combination of linear polarizers and linear wave retarders.

An ideal linear polarizer is a device that only transmits radiation for which the electric field lies within a plane. The intersection of this plane with the polarizer forms a line known as the transmission axis.

An ideal linear wave retarder is an optical device that applies a phase shift to radiation polarized in one direction with respect to radiation polarized in an orthogonal direction. The plane of linear polarization that corresponds to the direction of minimal phase retardation, when intersected with the surface of the wave retarder, yields a line known as the fast axis. The phase shift between radiation polarized parallel to and perpendicular to the fast axis is known as the retardance of the device.

By combining linear wave retarders and linear polarizers in series and varying their orientations and retardance values, it is possible to control the intensity of emitted radiation and to produce light in any state of linear, circular, or elliptical polarization.

Stokes Vectors and Optical Components

The intensity and polarization of a ray of light can be described by a 4-vector known as the Stokes vector, whose components are called the Stokes parameters. The first Stokes parameter is the ray intensity, the second and third parameters indicate linear polarization in various directions, and the fourth and final parameter corresponds to the degree of circular polarization. A more in-depth explanation of the physical meaning of the Stokes parameters can be found in the Ray Optics Module User’s Guide.

For example, a ray of natural (unpolarized) light with intensity I0 has Stokes vector

(1)

The Stokes vector of linearly polarized light varies depending on the direction of polarization. If the light is polarized in the direction of the x-axis, the Stokes vector is

(2)

For right-hand circularly polarized light the Stokes vector is

(3)

Varying states of elliptical polarization can be defined by specifying nonzero values of all of the Stokes parameters. In general, the light ray is fully polarized if the L2 norm of the second, third, and fourth parameters equals the first parameter or ray intensity.

In the Geometrical Optics interface, optical components such as linear polarizers and wave retarders are implemented as boundary conditions that do not affect the ray path but do change the Stokes parameters, which are stored as separate degrees of freedom for each ray. The effect of any optical component or system of optical components can be represented by a 4 × 4 Mueller matrix M. The Mueller matrix is multiplied by the Stokes vector of the incoming ray to produce a new Stokes vector for the transmitted ray,

where si is the Stokes vector of the incident ray. A rotation matrix R is included to account for the orientation of the optical component with respect to the coordinate system in which the ray’s Stokes parameters are defined.

Linear Polarizer

A linear polarizer transforms any incident light ray into a linearly polarized ray. Its Mueller matrix is

(4)

The orientation of the linear polarizer is specified by defining a direction called the transmission axis T. Light that is polarized in the direction parallel to the fast axis is completely transmitted, whereas light that is polarized in the orthogonal direction is completely blocked. The matrix R is then the rotation matrix from the coordinate system in which the Stokes parameters are defined to a coordinate system in which the transmission axis is parallel to the x-axis.

To see how the Mueller matrix in Equation 4 produces linearly polarized light, consider its effect on natural light (Equation 1):

(5)

The transmitted light has half the intensity of the incident light and is linearly polarized. Linear polarization in any other direction can be achieved by including the rotation matrix R in Equation 5. This is equivalent to changing the orientation of the linear polarizer; that is, rotating the transmission axis.

Linear Wave Retarder

A linear wave retarder applies a phase delay to linearly polarized light in one direction, relative to linearly polarized light in the orthogonal direction. Its Mueller matrix is

(6)

The orientation of the linear wave retarder is specified by defining a direction called the fast axis F. Light that is polarized in the direction parallel to the fast axis is subjected to the minimum phase delay, whereas light that is polarized in the orthogonal direction is subjected to the maximum phase delay. The matrix R is then the rotation matrix from the coordinate system in which the Stokes parameters are defined to a coordinate system in which the fast axis is parallel to the x-axis. The relative phase delay δ is called the retardance.

A linear wave retarder has no discernible effect on natural light; this can be seen by combining Equation 1 and Equation 6 for any value of δ:

A quarter-wave retarder (δ = π/2) can convert linearly polarized light into circularly polarized light if the incident light is linearly polarized at a 45° angle with respect to the fast axis:

A half-wave retarder (δ = π) can convert right-hand circularly polarized light into left-hand circularly polarized light and vice-versa:

Combining Optical Components

A combination of several optical components in series can be modeled by defining a Mueller matrix for each layer. Alternatively, it is possible to specify a single Mueller matrix, which is the product of the Mueller matrices for all of the optical components; if there are N optical components such that M1 is the Mueller matrix of the first component encountered and MN is the component of the last component encountered, the equivalent Mueller matrix Meq is

In this example, the individual polarizers and wave retarders are instead handled as separate entities, so that the effect of each optical component on the intensity and polarization of the transmitted light can be observed in greater detail.

Model Definition

The model geometry consists of three parallel surfaces, all of which are parallel to the xy-plane. The first and last boundaries are assigned the boundary condition, one with a transmission axis parallel to the x-axis and the other with a transmission axis parallel to the y-axis. The middle boundary is assigned the boundary condition with a fast axis parallel to the line y = x; it makes a 45° angle with the transmission axes of both polarizers.

Without the linear wave retarder, no radiation would propagate through the assembly of polarizers. However, by varying the retardance, the linearly polarized ray transmitted by the first polarizer can be converted to a ray of circular polarization, elliptical polarization, or linear polarization in a different direction before reaching the second polarizer, thus allowing some of the light to be transmitted.

Results and Discussion

When the linear wave retarder is given zero retardance, the series of optical devices consists of two linear polarizers with orthogonal transmission axes. As shown in Figure 1, the first linear polarizer reduces the ray intensity by half, as indicated by the color expression. The transmitted ray is also linearly polarized, as indicated by the deformation of the ray trajectory. The second linear polarizer reduces the ray intensity to zero.

When the linearly polarized ray is transmitted through a quarter-wave retarder, the intensity of the ray transmitted by the second linear polarizer is nonzero, as shown in Figure 2. The quarter-wave retarder converts the incident linearly polarized ray to a circularly polarized ray without changing its intensity. Then half of the intensity of the circularly polarized ray is transmitted through the second polarizer, which returns it to a state of linear polarization.

When the quarter-wave retarder is replaced with a half-wave retarder, the incident linearly polarized ray is instead converted to linearly polarized radiation with an orthogonal polarization. As a result, the second linear polarizer does not cause any noticeable decrease in the ray intensity, as shown in Figure 3.

Figure 1: A ray passing through two linear polarizers. The final intensity is zero.

Figure 2: Passage of a linearly polarized ray through a quarter-wave retarder between two linear polarizers. The final intensity is 1/4 of the initial intensity.

Figure 3: Passage of a linearly polarized ray through a half-wave retarder between two linear polarizers. The final intensity is 1/2 of the initial intensity.

Ray_Optics_Module/Tutorials/linear_wave_retarder

Modeling Instructions

From the menu, choose .

New

In the window, click

Model Wizard

In the window, click

In the tree, select .

Click .

Click

In the tree, select .

Click

Global Definitions

Parameters 1

In the window, under click .

In the window for , locate the section.

In the table, enter the following settings:


Name	Expression	Value	Description
delta	0	0	Retardance

Geometry 1

The geometry contains no domains and consists of three parallel surfaces, at which the boundary conditions will be applied.

Work Plane 1 (wp1)

In the toolbar, click

In the window for , locate the section.

In the text field, type 1.

Work Plane 1 (wp1)>Plane Geometry

In the window, click .

Work Plane 1 (wp1)>Square 1 (sq1)

In the toolbar, click

In the window for , locate the section.

From the list, choose .

Array 1 (arr1)

In the window, right-click and choose .

Select the object only.

In the window for , locate the section.

From the list, choose .

In the text field, type 3.

Locate the section. In the text field, type 1.

Click

Click the

button in the toolbar.

Geometrical Optics (gop)

Since no reflection at material discontinuities occurs in this model, set the number of secondary rays to zero.

In the window, under click .

In the window for , locate the section.

In the text field, type 0.

Locate the section. From the list, choose .

Release from Grid 1

In the toolbar, click

and choose .

In the window for , locate the section.

Specify the L0 vector as


0	x
0	y
1	z

Create a pair of linear polarizers with orthogonal transmission axes.

Linear Polarizer 1

In the toolbar, click

and choose .

Select Boundary 1 only.

Linear Polarizer 2

In the toolbar, click

and choose .

Select Boundary 3 only.

In the window for , locate the section.

Specify the T vector as


0	x
1	y
0	z

Next, add a quarter-wave retarder to enable a fraction of the energy of the ray to propagate through the second polarizer.

Linear Wave Retarder 1

In the toolbar, click

and choose .

Select Boundary 2 only.

In the window for , locate the section.

Specify the F vector as


1	x
1	y
0	z

In the δ text field, type delta.

Study 1

Parametric Sweep

In the toolbar, click

The study uses three different retardance values. For δ = 0 the boundary condition has no effect. For δ =π /2 a quarter-wave retarder is used. For δ =π a half-wave retarder is used.

In the window for , locate the section.

Click

In the table, enter the following settings:


Parameter name	Parameter value list	Parameter unit
delta (Retardance)	0 pi/2 pi	rad

Step 1: Ray Tracing

In the window, click .

In the window for , locate the section.

From the list, choose .

Click

In the dialog box, type 0.1 in the text field.

In the text field, type 4.

Click .

In the toolbar, click

Results

No Wave Retarder

In the window, expand the node, then click .

In the window for , locate the section.

From the list, choose .

In the text field, type No Wave Retarder.

Locate the section. Clear the check box.

Modify the default plot to indicate the intensity and polarization of the ray.

Ray Trajectories 1

In the window, click .

In the window for , locate the section.

Find the subsection. From the list, choose .

Color Expression 1

In the window, expand the node, then click .

In the window for , locate the section.

Click

In the dialog box, select in the tree.

Click .

Ray Trajectories 1

In the window, click .

In the window for , locate the section.

Find the subsection. From the list, choose .

In the text field, type 25.

Select the check box. In the associated text field, type 0.3.

Surface 1

In the window, right-click and choose .

In the window for , locate the section.

From the list, choose .

Locate the section. From the list, choose .

From the list, choose .

In the toolbar, click

Click the

button in the toolbar.

Compare the resulting image to Figure 1.

Quarter-Wave Retarder

Right-click and choose .

In the window for , type Quarter-Wave Retarder in the text field.

Locate the section. From the list, choose .

In the toolbar, click

Compare the resulting image to Figure 2. The quarter-wave retarder causes a phase shift between the electric field components parallel to and perpendicular to the fast axis. As a result, the linearly polarized ray becomes circularly polarized. Upon reaching the second linear polarizer, the ray intensity is reduced by half.

Half-Wave Retarder

Right-click and choose .

In the window for , type Half-Wave Retarder in the text field.

Locate the section. From the list, choose .

In the toolbar, click

Compare the resulting image to Figure 3. While propagating through the half-wave retarder, the ray remains linearly polarized, but the direction of polarization is rotated. The ray then propagates through the second linear polarizer without any loss of intensity.