26

Birefringence Measurement

Yukitoshi Otani

CONTENTS

26.1  Introduction

26.2  Applications of Retardance and Birefringence Measurement

26.3  About Polarized Light

26.4  Polarization State

26.5  What Is Birefringence?

26.6  Notation of Polarization States

26.6.1  Jones Vector

26.6.2  Stokes Parameter

26.6.3  Poincaré Sphere

26.7  Jones Method and Mueller Method

26.8  Birefringence Measurement

26.8.1  Linear Polariscope

26.8.2  Circular Polariscope

26.8.3  Senarmont Method

26.8.4  Photoelastic Modulator Method

26.8.5  Optical Heterodyne Method

26.8.6  Two-Dimensional Birefringence Measurement by Phase-Shifting Method

References

26.1    INTRODUCTION

A number of polarization and birefringence users, also known as a polarization society, have shown remarkable development in the last few decades. Most of the recent advanced technologies are based on the polarization technology including birefringence [1,2,3,4]. Many new interesting topics such as liquid crystal (LC) displays, LC projectors, optical disks, and optical communications are proposed using the polarization effect. A birefringence measurement is one of the visualizing methods for internal conditions of materials. We can understand information on stress distribution and macromolecular orientation distribution inside materials from birefringence mapping.

This chapter describes the polarization situation mathematically and explains how birefringence distribution is measured. As we have to handle the two parameters, birefringence perplexes us, because birefringence has a magnitude and an azimuthal direction and is not a scalar.

Birefringence measurement has a long history that originated from photoelasticity [5]. A polariscope for the photoelastic phenomenon was first discovered by Sir David Brewster in 1816. We can know the stress distribution from the fringe intensity. Photoelasticity has been expanded widely to the field of experimental mechanics. However, its activities have been broken down after computer simulation as the famous finite element method. Moreover, fringe counting is difficult because precise measurements are required. We have to measure birefringence that is less than one fringe. Polarimetric measurements, especially Mueller matrix polarimetry, are powerful tools for mapping birefringence [6]. It is also possible to analyze birefringence parameters from elements of a partial matrix. Moreover, if we can also neglect the diattenuation, which means no absorption dependences, we could propose many unique methods for mapping birefringence.

26.2    APPLICATIONS OF RETARDANCE AND BIREFRINGENCE MEASUREMENT

Recently, many optical instruments have used polarization effects. Table 26.1 shows remarkable applications of polarization and birefringence effects. We can divide features of optical applications into two groups as information and energy. Most of the applications are used for information such as manufacture, optical communication, biotechnology, nanotechnology, astronomy, and remote sensing, but some reports for energy application become increasingly unique applications of polarization. Polarization has contributed to the technological advances of recent years.

Representative examples are an LC display, an optical pickup for data storage, an optical isolator, and an optical multiplexer for optical communications. Birefringence mapping is an important method to visualize and analyze crystal orientation, stress distribution, polymer molecule orientation, and thin films for any kind of material. In optically isotropic materials, we have a problem with birefringence when the materials are applied to a stress field. This can possibly be explained by the fact that the speed of light changes depending on the stress fields caused by the force inside the material because of mass density. External forces make the atomic distance narrower along the force and broaden the distance perpendicular to it. In recent years, many optical instruments have become small and precise. Deterioration in quality causes birefringence, even if it is small. For example, birefringence in pickup lens is caused by astigmatism. Some experimental results are shown later.

TABLE 26.1
Remarkable Applications of Polarization and Birefringence

Feature

Area

Applications

Characteristics

Information

Manufacturing technology (LC)

Phase film, polarization film

Development, orientation

Substrate

Polarization

Birefringence

Birefringence

Manufacturing technology (pickup)

Optical elements

Birefringence

Manufacturing technology (injection process)

Orientation of polymer molecule

Birefringence

Manufacturing technology (semiconductor photolithography)

Optical elements

Birefringence

Manufacturing technology (metrology)

White-light interferometry

Geometric phase

Nanotechnology

Polarization elements

Structural birefringence

Optical communications

Polarization analysis

Optical isolator

Polarization control

Polarization

Optical rotation

Polarization

Biotechnology

Glucose sensor

Cancer diagnosis

Cell observation

Visualization

Orientation of polymer molecule

Optical rotation

Polarization

Birefringence

Polarization

Fluorescent polarization

Astronomy

Coronagraph

Observation of radiation mechanism

Geometric phase

Zeeman effect

Remote sensing

Visualization

Earth observation (surface, aerosol, cloud, sea)

Polarization

Polarization

Energy

Manufacturing process

Laser processing

Polarization

Manipulation

Laser trapping

Polarization

Major optical devices using polarized light are LC displays in optical displays, pickups in optical memories, and multiplexing in optical communications or isolators. In addition, the acquisition of internal information is possible in the field of measurement. Nowadays, the acquisition of an object’s internal information is highly required in the field of optoelectronics performing in various fields such as material development including glasses and transparent plastics, LC displays, optical parts such as polymers and crystals, or optical disks and thin film products and also in the development of new products. In addition, in practical use, because the shape of optoelectronics is ever more complex, specialized, and large-sized, the material’s original double refraction and the double refraction due to residual stresses created inside while processing cause various problems. In addition, orientation information, important in the field of polymers, can be taken as double refraction.

In recent years, polarized light has been used in a wide variety of fields. Especially, its importance is increasing due to the recent discovery of high-resolution optical parts. The term polarized light has many meanings such as double refraction, optical rotating power, dichroism, and depolarization. Here, we set our goal to the active use of polarizers. Thus, we begin with the questions “What is polarized light? and What is double refraction?” and then discuss the basis to actively use polarized light. It would be a pleasure if this would help in the understanding of research papers and catalogs on polarized light or the introduction of polarimeters.

26.3    ABOUT POLARIZED LIGHT

Although light is said to have wave nature and particle nature, only particle nature is considered when dealing with polarized light. That is, light can be treated as a lateral electromagnetic wave. Consequently, light consists of electric vibration and magnetic vibration, but generally, only electric vibration is considered. In this case, we deal with electric vibration by decomposing it in two perpendicular vibrating directions. Suppose the z-axis to be in the direction of light in the right-hand coordinate system, x-axis in the horizontal direction, and y-axis in the vertical direction. Just for reference, the wavelength of visible light is 400–700 nm.

Since light is a lateral wave, an electric vibrating plane is formed on the xz surface as shown in Figure 26.1. This plane is called plane of polarization. The H-axis represents the magnetic vibrating plane. If a polarizer made with small enough gratings compared to light wave or with absorptiveness in a particular plane of polarization is installed, light wave penetrates or is intercepted depending on the vibrating direction. That is, light that has penetrated a polarizer will have only one vibrating direction. Here, since the vibrating direction is linear, it is called linear polarization, and the plane in the vibrating direction is called plane of polarization. Generally, p-polarized light (parallel in German) and s-polarized light (senkrecht) are used to express the polarization state. When light incidents upon a medium, the angle between the incident light’s direction and normal to the surface is called the angle of incidence. The vibrating components inside the surface are called p-polarized light, and those normal to the surface, s-polarized light.

26.4    POLARIZATION STATE

Polarization states are linear polarization, circular polarization, and elliptical polarization. Let us consider these polarization states. As mentioned earlier, a light wave can be decomposed into xy-axis electric vibration vector components. Figure 26.2 shows the relationship between the electric vibration vectors Ex and Ey before and after passing through a sample with the retardance δ(–2π(none)d/λ). Suppose A to be the amplitude, t the time, ω the angular frequency, and λ the wavelength, the electric vibration components Ex and Ey can be expressed as

Image

FIGURE 26.1  Electric vector of light. (a) Polarized light and (b) electric vector of light.

Image

FIGURE 26.2  Electric field before and after passing through a sample.

Ey=Ayexpi(ωt2πλz+φy)Ex=Axexpi(ωt2πλz+φx)Ey=Ayexpi(ωt2πλz+φy)Ex=Axexpi(ωt2πλz+φx)

(26.1)

provided that the phase difference δ is

δ =φyφxδ=φyφx

(26.2)

Image

FIGURE 26.3  Composition of the electric vibration vectors Ex and Ey. (a) Linear polarization in case of Ay = 0, (b) linear polarization in case of Ax = 0, and (c) linear polarization Ax = Ay, δ = 0.

It is possible to determine the statuses of linear polarization, circular polarization, and elliptical polarization from the amplitude A and the phase difference δ.

As shown in Figure 26.3, linear polarization takes place when the phase difference is 0 or π. As shown in Figure 26.3a, horizontal linear polarizations vibrating in the x-axis direction are obtained when the amplitude is Ax = 0, vertical linear polarization when Ay = 0 as shown is Figure 26.3b, and linear polarization with azimuth 45° when Ax = Ay = 1, as shown in Figure 26.3, because it is the composition of the amplitude vectors of x and y directions.

Next, circular polarization takes place when δ is 90° as shown in Figure 26.4. In this case, the relationship between the phase differences of Ex and Ey is given as a composition of vectors as shown in Figure 26.4a and traces the locus shown in Figure 26.4c. When observed from the direction of light, the vector’s rotation is counterclockwise as shown in Figure 26.4b, and so the locus will be as shown in Figure 26.4c. Thus, this kind of polarization state is called left-hand circular polarization. In addition, this is when Ey has a phase delay of 90°(–90°) when compared with Ex. On the other hand, if Ey has a phase lead of 90°(+90°), the locus is clockwise as opposed to Figure 26.4. This kind of polarization state is called right-hand circular polarization.

With these facts considered, Equation 26.1 can be expressed as

(Eyay)2+(Exax)22(Eyay)(Exax)cosδ =sin2δ(Eyay)2+(Exax)22(Eyay)(Exax)cosδ=sin2δ

(26.3)

Linear polarization is obtained when the phase difference is 0° or 180°. When δ is 90° and Ax = Ay, Equation 26.3 is expressed in the equation of a circle. When the phase difference is anything other than these, the polarization state will have its vectors directed elliptically, and its locus will be as shown in Figure 26.4c. This kind of polarization state is called elliptical polarization.

Image

FIGURE 26.4  Left circular polarization light. (a) Polarization state, (b) direction of polarization vector, and (c) trajectory.

26.5    WHAT IS BIREFRINGENCE?

Birefringence has, as it is written, two indexes of refraction. That is, the propagation velocity differs depending on the polarization direction. Figure 26.5 shows how an optical wave transmits double refractive objects with indexes of refraction n// and n, respectively. Figure 26.5b separates and aligns vertically each plane of polarization of Figure 26.5a. Phase difference Δ depending on the double refraction occurs between orthogonal polarizations. In Figure 26.5, since light Ex vibrating in the x direction is transmitted through the sample slowly, birefringent phase difference (retardation) depending on the transmitting velocity arises in the light emitted from the sample. In other words, the transmitting velocity in the x direction differs from that in the y direction. That is, their indexes of refraction are different. Since the refractive index is the ratio between the velocity of light c′ transmitted through a substance and the velocity of light c in vacuum,

n=ccn=cc

(26.4)

the double refraction Δn is given as

Δn=nn//Δn=nn//

(26.5)

Generally, birefringent phase difference Δ,

Δ=2πl/λΔnΔ=2πl/λΔn

(26.6)

Image

FIGURE 26.5  An example that the light passes though a birefringence sample. (a) Birefringence sample and (b) difference of light speed along the polarization difference refractive index.

TABLE 26.2
Determination of the Retardation Unit

Waves

1 Wave (λ)

Degrees

360°

Radians

2π rad

Nanometers

λ nm

λ: wavelength

since phase display (2π = 360°) can be expressed as follows:

Δ=360Δnl/λΔ=360Δnl/λ

(26.7)

and can be standardized by wavelength (nm) in Table 26.2.

The direction in which light is fastly transmitted (phase lead) is called the fast axis. On the other hand, the direction in which light is slowly transmitted (phase delay) is called the slow axis. The general term for the fast axis and the slow axis is the principal axis. In addition, samples with double refraction are called anisotropic and those without double refraction, such as glasses, are called isotropic in optics.

Thus, determination of the size (birefringent phase difference) and the principal axis direction (the fast axis and the slow axis) of a double refraction is needed to measure inner information of a sample.

The s-polarization is called (transverse electric) TE-mode, and p-polarization is called (transverse magnetic) TE-mode.

26.6    NOTATION OF POLARIZATION STATES

The notation of polarization states can be classified into those using visual techniques and those using mathematical techniques. The Poincaré sphere is a visual technique, and the Mueller method and the Jones method are mathematical methods.

26.6.1    JONES VECTOR

To express the polarization state mathematically, let us refer to Equation 26.1, which expresses optical waves decomposed in x and y components. When this equation is expressed in a matrix, it becomes

E=[AxAx2+Ax2exp(iδ2)AxAx2+Ay2exp(iδ2)]E=AxAx2+Ax2exp(iδ2)AxAx2+Ay2exp(iδ2)

(26.8)

For, only relative phases should be considered when dealing with general polarization states:

E=[ExEy]=expi(ωt2πλ+φx)[AxAyexpiδ]E=[ExEy]=expi(ωt2πλ+φx)[AxAyexpiδ]

(26.9)

In addition, if the polarization status is the only one to be considered, the amplitude must be standardized. When the phase is also transformed,

E=[AxAy]=[AxAyexpi(δ)]E=[AxAy]=[AxAyexpi(δ)]

(26.10)

For example, since in horizontal linear polarizations Ax = 1, Ay = 0, and δ = 0, Jones vector J is obtained as follows by using the above equation.

J=[10]J=[10]

(26.11)

Also, for instance, in clockwise circular polarizations, Ax = 1, Ay = 1, and δ = 90°. This is when Ey has a lead of 90°(+90°), and Jones vector J is expressed as follows by using the complex number i, the same as the notation of an electric current.

J=12[1i]J=12[1i]

(26.12)

All polarization states can similarly be expressed with Jones vectors. The Jones vectors of representative polarization states are shown in Table 26.3.

TABLE 26.3
Various Forms of Standard Normalized Jones Vectors

Image

26.6.2  STOKES PARAMETER

Here, unlike the Jones vector, which uses two elements to express a polarization state, we use the Stokes parameter made up of four elements, s0 to s3.

Let us refer to Equation 26.1, which expresses optical waves decomposed in x and y components, just like the Jones vector. When this equation is expressed in matrices, it becomes

s0=ExE*x+EyE*y=A2x+A2ys1=ExE*xEyE*y=A2xA2ys2=E*xEy+ExE*y=2AxAycosΔs3=iE*xEyExE*y=2AxAysinΔwherei2=1s0=ExEx+EyEy=A2x+A2ys1=ExExEyEy=A2xA2ys2=ExEy+ExEy=2AxAycosΔs3=iExEyExEy=2AxAysinΔwherei2=1

(26.13)

The four elements could also be expressed with I, M, C, S or I, Q, U, V. When these are put together in the matrix S,

S=[s0s1s2s3]=[IMCS]=[IQUV]S=s0s1s2s3=IMCS=IQUV

(26.14)

In the Stokes parameter, s0 stands for the light intensity, stands for the horizontal linear polarization component, s2 stands for the 45° straight line, and s3 stands for the right-hand circular polarization. Generally, light intensity is treated as the unit intensity.

Based on these, we deal with representative polarization states. For example, since Ax = 1, Ay = 0, and δ = 0 in horizontal linear polarization,

S=[1100]S=1100

(26.15)

Likewise, in right-handed circular polarization, Ax = 1, Ay = 1, and δ = 90°.

S=[1001]S=1001

(26.16)

All these forms are shown in Table 26.4.

TABLE 26.4
Various Forms of Standard Normalized Stokes Vectors

Image

The relationship between the Stokes parameter’s four elements s0 to s3 depends on the polarization state.

In totally polarized lights,

s20=s21+s22+s23s20=s21+s22+s23

(26.17)

In partially polarized lights,

s20>s21+s22+s23s20>s21+s22+s23

(26.18)

Thus, the degree of polarization V is defined as an index of how much light is polarized.

V=s21+s22+s23s0V=s21+s22+s23s0

(26.19)

26.6.3    POINCARÉ SPHERE

The possibility of expressing the polarization states mathematically was characteristic to the Jones vector and the Stokes parameter. However, the polarization states, hard enough to perceive, are even more difficult to understand when expressed in matrices. The Poincaré sphere is a method to visually understand the polarization states.

As mentioned in the section dealing with the Stokes parameter, since Equation 26.17 is satisfied in totally polarized lights, when each axis of spatial Cartesian coordinates is expressed as s1, s2, and s3, a sphere with a radius of s0 is obtained. In addition, if the coordinates of point S, located on the sphere, are expressed with latitude 2β and longitude 2θ, the following equations are obtained.

s1=s0cos2βcoss2=s0cos2βsins3=s0sin2βs1=s0cos2βcos2θs2=s0cos2βsin2θs3=s0sin2β

(26.20)

We will not go into details, but the correspondence to the azimuth θ of the long axis and to the ellipticity angle β in elliptical polarizations can be considered. This situation is illustrated in Figure 26.6. It is easier to understand when the Poincaré sphere is treated as the Earth. A right-hand polarization and a left-hand polarization are created in the upper and lower hemispheres respectively. On the equator, a linear polarization is obtained by the ellipticity of 0°. A right-hand and a left-hand circular polarization are acquired in the upper and the lower poles, respectively, and elliptical polarizations are distributed between each pole and the equator. Polarization states symmetrical about the center of the sphere have a difference of 90° in the azimuth angle of the major axis.

Image

FIGURE 26.6  Poincaré sphere representation of polarization state.

26.7    JONES METHOD AND MUELLER METHOD

All polarization states could be described using the Jones vector and the Stokes parameter. As shown in Figure 26.7, once the incident light and the output light’s polarization states are identified, the polarization elements of an optical device, which is a black box, can be determined subsequently.

The Jones calculus is written as

JE=E[JxxJxyJyxJyy][ExEy]=[ExEy]JE=E[JxxJyxJxyJyy][ExEy]=[ExEy]

(26.21)

Image

FIGURE 26.7  Jones matrices and Mueller calculus.

The Mueller calculus is indicated as

MS=S[M00M01M02M03M10M11M12M13M20M21M22M23M30M31M32M33][s0s1s2s3]=[s0s1s2s3]MS=SM00M10M20M30M01M11M21M31M02M12M22M32M03M13M23M33s0s1s2s3=s0s1s2s3

(26.22)

The polarization elements that can be obtained from these methods are put together in Table 26.5.

Only multiplication of the matrices is needed if the polarization elements are lined together. Additionally, the light intensity obtained at the end in the case of Jones matrix is

I=Ex+Ey2I=Ex+Ey2

(26.23)

In the case of Mueller matrix, the light intensity is directly given by s0.

To conclude, the polarization degrees of freedom of the Jones and the Mueller matrices are as shown in Figure 26.7. Both the Jones method and the Mueller method give double refraction, double absorption (dichroism), circular double refraction (optical rotating power), circular double absorption (circular double dichroism), and others.

Their major difference is that as opposed to the Jones vector treating only totally polarized lights, the Mueller matrix can deal with all polarization states including depolarization.

26.8    BIREFRINGENCE MEASUREMENT

26.8.1    LINEAR POLARISCOPE

In this section, we progress to observe and measure stress distribution. Birefringence distribution can be observed from an intensity distribution inserted in a sample between two polarizers as in Figure 26.8. If you analyze this fringe, you can determine the quantity of stress.

In this chapter, three nonmodulated methods and four modulated methods are introduced in the following pages.

Figure 26.9 shows two polarizers and a light source. The optical axes of the polarizer are set perpendicular to each other. A sample is inserted between the polarizers.

The relation among the polarization elements is given using a Mueller matrix and Stokes vector by

S=P(0°)R(Δ,θ)P(90°)SS=P(0°)R(Δ,θ)P(90°)S

(26.24)

where the Stokes vector S′ and S are the output and input of the light, P(0°), R(Δ, θ), and P(90°) are Mueller matrix of polarizer with 0° of azimuthal direction, sample with Δ retardance and θ azimuthal direction and polarizer with 90° of azimuthal direction, respectively.

The retardance Δ indicates stress-induced birefringence by σ1 and σ2.

[s0s1s2s3]=12[1100110000000000][10000cos22θ+sin22θcosΔ(1cosΔ)sin2θcossin2θsinΔ0(1cosΔ)sin2θcossin22θ+cos22θcosΔcos2θsinΔ0sin2θsinΔcos2θsinΔcosΔ]×12[1100110000000000][1000]=14[sin22θ(1cosΔ)sin22θ(1cosΔ)00]=12[sin22θsin2(Δ/2)sin22θsin2(Δ/2)00]s0s1s2s3=12110011000000000010000cos22θ+sin22θcosΔ(1cosΔ)sin2θcos2θsin2θsinΔ0(1cosΔ)sin2θcos2θsin22θ+cos22θcosΔcos2θsinΔ0sin2θsinΔcos2θsinΔcosΔ×1211001100000000001000=14sin22θ(1cosΔ)sin22θ(1cosΔ)00=12sin22θsin2(Δ/2)sin22θsin2(Δ/2)00

(26.25)

TABLE 26.5
Polarization Elements of Jones Matrix and Mueller Matrix

Image

Image

FIGURE 26.8  Visualization of birefringence.

Image

FIGURE 26.9  Plane polariscope.

The intensity I observed at the eye is given from the first term of s0

I=s0=12sin22θsin2(Δ/2)=14sin22θ(1cosΔ)I=s0=12sin22θsin2(Δ/2)=14sin22θ(1cosΔ)

(26.26)

This means the intensity changes sinusoidally along the retardance Δ and the azimuth angle makes the intensity of the periodic circle 90° independently.

Figure 26.10 shows stress mapping by the plane polariscope [16] (http://www.luceo.co.jp). A plastic disk substrate and square plate are used as samples.

In Figure 26.10a, the intensity distribution of a disk substrate changes bright to dark radially because of the birefringence distribution. The intensity distribution is modulated sinusoidally along the circumferential direction and this period is 90°. This means the azimuth direction of the birefringence indicates concentricity.

The three pictures in Figure 26.10b through d mean merely rotating the sample at the same stress condition. The stress around the surrounding part is larger than that at the center section. This intensity distribution is also modulated sinusoidally in respect of the sample rotation.

Image

FIGURE 26.10  Stress mapping by the plane polariscope. (a) Plastic disk substrate, (b) stressed plastic plate stress direction is 0°, (c) rotation angle of 45°, and (d) rotation angle of 90°.

It is difficult to determine both the retardance and azimuth direction simultaneously due to the image of the plane polariscope.

26.8.2    CIRCULAR POLARISCOPE

Figure 26.11 shows a circular polariscope. It consists of a circular polarizer and a circular analyzer. A sample is set between quarter wave plates. The relation among the polarization elements is given using a Mueller matrix and Stokes vector by

S=P(0°)Q(45°)R(Δ,θ)Q(45°)P(90°)SS=P(0°)Q(45°)R(Δ,θ)Q(45°)P(90°)S

(26.27)

where

the Stokes vector S′ and S are the output and input of the light

P and Q are Mueller matrix of polarizer and quarter wave plate with azimuthal direction angle between parentheses, respectively

R(Δ, θ) is a Mueller matrix of the sample with Δ retardance and θ azimuthal direction

The retardance Δ indicates stress-induced birefringence by σ1 and σ2. All Stokes parameters and Mueller matrices are shown as

[s0s1s2s3]=12[1100110000000000][1000000100100100]×[10000cos22θ+sin22θcosΔ(1cosΔ)sin2θcossin2θsinΔ0(1cosΔ)sin2θcossin22θ+cos22θcosΔcos2θsinΔ0sin2θsinΔcos2θsinΔcosΔ]×[1000000100100100]12[1100110000000000][1000]=14[1cosΔ1cosΔ00]s0s1s2s3=1211001100000000001000000100100100×10000cos22θ+sin22θcosΔ(1cosΔ)sin2θcos2θsin2θsinΔ0(1cosΔ)sin2θcos2θsin22θ+cos22θcosΔcos2θsinΔ0sin2θsinΔcos2θsinΔcosΔ×10000001001001001211001100000000001000=141cosΔ1cosΔ00

(26.28)

Image

FIGURE 26.11  Circular polariscope.

Image

FIGURE 26.12  Stress mapping by the circular polariscope. (a) Plastic disk substrate, (b) stressed plastic plate stress direction is 0°, (c) rotation angle of 45°, and (d) rotation angle of 90° (http://www.luceo.co.jp).

Finally, the intensity I observed at the eye is given from the first term of s0

I=s0=14(1cosΔ)=12sin2(Δ/2)I=s0=14(1cosΔ)=12sin2(Δ/2)

(26.29)

The intensity distribution means the retardation is changing sinusoidally along the retardance Δ but the azimuth angle is independent of the intensity change.

Figure 26.12 shows stress mapping by the circular polariscope (http://www.luceo.co.jp), the same as that in Figure 26.10.

26.8.3    SENARMONT METHOD

The requirement of birefringence measurement includes stress mapping to measure the state of polarization following the sample, to measure the state of the polarization of incident beam to the sample, or to determine the state of polarization in advance. In addition, it measures the alignment of direction of the sample to incident light or measures both retardation and azimuthal direction. Below are some measurement methods including the Senarmont method, the optical heterodyne method, and the two-dimensional birefringence measurement, using the phase-shifting technique.

The Senarmont method is a classical method for stress mapping. Because it is economical, it is often used for industrial inspection. It has features such as easy alignment, an equal method of ellipsometry, and high accuracy even when using a low-quality wave plate. An optical arrangement of the Senarmont method is shown in Figure 26.13. It consists of a polarizer, a quarter wave plate with 0°, and a rotating analyzer. A sample is set between quarter wave plates with the fixed value of azimuthal direction as 45°. The relation among the polarization elements is given using a Mueller matrix and Stokes vector by

Image

FIGURE 26.13  Optical arrangement of Senarmont method.

S=P(θ)Q(0°)R(Δ,45°)P(90°)SS=P(θ)Q(0°)R(Δ,45°)P(90°)S

(26.30)

where

the Stokes vector S’ and S are the output and input of the light

P and Q are Mueller matrices of polarizer and quarter wave plate with azimuthal direction angle between parentheses, respectively

θ is the azimuthal direction of the analyzer

R(Δ,45°) is a Mueller matrix of the sample with Δ retardance and 45° azimuthal direction

The retardance Δ indicates stress-induced birefringence by σ1 and σ2 as shown in Equation 26.11. All Stokes parameters and Mueller matrices are shown as

[s0s1s2s3]=12[1cos2θsin2θ0cos2θcos22θsin2θcos0sin2θsin2θcossin22θ00000][1000010000010010][10000cosΔ0sinΔ00100sinΔ0cosΔ]×12[1100110000000000][1000]=14[1cosΔcos2θsinΔcoscos2θcos22θcosΔsin2θcossinΔsin2θsincos2θcosΔsin22θsinΔ0]s0s1s2s3=121cos2θsin2θ0cos2θcos22θsin2θcos2θ0sin2θsin2θcos2θsin22θ00000100001000001001010000cosΔ0sinΔ00100sinΔ0cosΔ×1211001100000000001000=141cosΔcos2θsinΔcos2θcos2θcos22θcosΔsin2θcos2θsinΔsin2θsin2θcos2θcosΔsin22θsinΔ0

(26.31)

Finally, the intensity I observed at the eye is given from the first term of s0

I=s0=1cos(Δ2θ)I=s0=1cos(Δ2θ)

(26.32)

The intensity modulates sinusoidally along the retardance Δ and the azimuthal direction θ of the analyzer. This means we can determine the retardation Δ from the azimuthal angle θ of the analyzer if we find the dark fringe by the rotating analyzer. The retardation is

Δ=2θΔ=2θ

(26.33)

26.8.4    PHOTOELASTIC MODULATOR METHOD

One of the popular methods used for many kinds of industries by Hinds is a birefringence measurement system using a photoelastic modulator (PEM) [7,8,9,10]. The optical arrangement consists of a polarizer, a PEM with 0° azimuthal direction, and an analyzer with 0° and 45° azimuthal angle (Figure 26.14). The PEM modulates the retardance δ with temporal modulation. A sample is set between PEM and the analyzer. The relation among the polarization elements is given using Mueller matrices and Stokes vectors by

Image

FIGURE 26.14  Optical arrangement of birefringence measurement.

S=P(45°)X(Δ,θ)PEM(δ,0°)P(45°)SS=P(45°)X(Δ,θ)PEM(δ,0°)P(45°)S

(26.34)

where

the azimuthal angle of the analyzer is 45°

the Stokes vector S′ and S are the output and input of the light, respectively

P is the Mueller matrix of the polarizer

PEM(δ, 0) is the Mueller matrix of the polarizer with the temporal phase modulation δ

X(Δ, θ) is the Mueller matrix of the sample with Δ retardance and θ azimuthal direction

The retardance Δ also indicates stress-induced birefringence by σ1 and σ2 as shown in Equation 26.11. All Stokes parameters and Mueller matrices are shown as

[s0s1s2s3]=12[1010000010100000][10000cos22θ+sin22θcosΔ(1cosΔ)sin2θcossin2θsinΔ0(1cosΔ)sin2θcossin22θ+cos22θcosΔcos2θsinΔ0sin2θsinΔcos2θsinΔcosΔ]×[1000010000cosδsinδ00sinδcosδ][1010]s0s1s2s3=12101000001010000010000cos22θ+sin22θcosΔ(1cosΔ)sin2θcos2θsin2θsinΔ0(1cosΔ)sin2θcos2θsin22θ+cos22θcosΔcos2θsinΔ0sin2θsinΔcos2θsinΔcosΔ×1000010000cosδsinδ00sinδcosδ1010

(26.35)

Finally, the detected intensity I is given from the first term of s0

I=s0=I0[1+sinδ(cos2θsinΔ)cosδ(sin22θ+cos22θcosΔ)]=I0[1+sin(δ0sinωt)cos2θsinΔcos(δ0sinωt)(sin22θ+cos22θcosΔ)]=I0{1+cos2θsinΔ[2J1(δ0)sinωt](sin22θ+cos22θcosΔ)[J0(δ0)+2J2(δ0)cos2ωt]}I=s0=I0[1+sinδ(cos2θsinΔ)cosδ(sin22θ+cos22θcosΔ)]=I0[1+sin(δ0sinωt)cos2θsinΔcos(δ0sinωt)(sin22θ+cos22θcosΔ)]=I0{1+cos2θsinΔ[2J1(δ0)sinωt](sin22θ+cos22θcosΔ)[J0(δ0)+2J2(δ0)cos2ωt]}

(26.36)

where Jn is the nth order’s Bessel function.

The phase modulation δ is

δ =δ0sinωtδ=δ0sinωt

(26.37)

where

t is the time

δ0 is the maximum amplitude of retardance

ω is the angular frequency

Using appropriate modulation δ0 we can eliminate the J0 and J1 terms. If the azimuthal angle θ is 0°, the retardation is obtained using a lock-in amplifier:

Δ=sin1[kI(ω)J1(δ0)]Δ=sin1[kI(ω)J1(δ0)]

(26.38)

In this case, the calibration process is to determine the parameter k.

One more equation is required to determine both retardation and azimuthal direction, simultaneously. The different angles of the analyzer such as 0° are written as

S=P(0°)X(Δ,θ)PEM(δ,0°)P(45°)SS=P(0°)X(Δ,θ)PEM(δ,0°)P(45°)S

(26.39)

The detected intensity I is given from the first term of s0,

I45°=I0{1+sin(δ0sinωt)sin2θsinΔ+cos(δ0sinωt)((1cosΔ)sin2θcossinsinΔ)}=I0{1+sin2θsinΔ[2J1(δ0)sinωt]((1cosΔ)sin2θcossinsinΔ)[J0(δ0)+2J2(δ0)cos2ωt]}I45°=I0{1+sin(δ0sinωt)sin2θsinΔ+cos(δ0sinωt)((1cosΔ)sin2θcos2θsin2θsinΔ)}=I0{1+sin2θsinΔ[2J1(δ0)sinωt]((1cosΔ)sin2θcos2θsin2θsinΔ)[J0(δ0)+2J2(δ0)cos2ωt]}

(26.40)

Using appropriate modulation δ0 we can eliminate the J0 and J1 terms. The amplitude signal of a lock-in amplifier is Δ at the 0° of the analyzer and Δ45° at the 45° of the analyzer. Therefore, retardation is observed by

θ =12tan1Δ45°Δ0°θ=12tan1Δ45°Δ0°

(26.41)

and the azimuthal direction is

Δ=Δ20°+Δ245°Δ=Δ20°+Δ245°

(26.42)

Figure 26.15 is a measured result of birefringence distribution of CaF2 [111]. CaF2 is only a lens material for UV lithography. A He–Ne laser with 632.8 nm wavelength is used for a light source. The size of material is 40 mm diameter and 7 mm thickness. Retardation distribution is 0.5–7 nm. This distribution is caused by residual stress. If you change a light source to the vacuum ultraviolet such as 157 nm, intrinsic birefringence could be observed.

Image

FIGURE 26.15  Optical arrangement of birefringence measurement by optical heterodyne method.

26.8.5  OPTICAL HETERODYNE METHOD

The third method is a birefringence measurement using optical heterodyne detection [10]. Figure 26.15 is an optical arrangement. It consists of a stabilized transverse Zeeman laser (STZL), a half wave plate, and the linear polarizer. The laser STZL has two orthogonal linearly polarized modes and their angular frequencies are ω1 and ω2, and ωb = |ω1 – ω2|. The frequency difference is ωb ≈ 105 rad. The half wave plate and linear polarizer are mounted on a rotating stage. The rotating angle of HWP is 2θ and the polarizer is always 45° between two beams so that they interfere with the sample; when two linear polarizations are aligned with the sample fast axis, the phase shifts one way. About 90° later, the phase shifts in the opposite direction.

The relation among the polarization elements is given using Mueller matrices and Stokes vectors by

S=LP2θ+45XΔ,φHWθSS=LP2θ+45XΔ,φHWθS

(26.43)

The Stokes vector is expressed by

S=[ax2+ay2ax2ay22axaycosωbt2axaysinωbt]S=ax2+ay2ax2ay22axaycosωbt2axaysinωbt

(26.44)

All Stokes parameters and Mueller matrices are shown as

[s0s1s2s3]=12[1cos(4θ90°)sin(4θ+90°)0cos(4θ+90°)cos2(4θ+90°)sin(4θ+90°)cos(4θ+90°)0sin(4θ+90°)sin(4θ+90°)cos(4θ+90°)sin2(4θ+90°)00000][10000cos22ϕ+sin22ϕcosΔ(1cosΔ)sin2ϕcos2ϕsin2ϕsinΔ0(1cosΔ)sin2ϕcos2ϕsin22ϕ+cos22ϕcosΔcos2ϕsinΔ0sin2ϕsinΔcos2ϕsinΔcosΔ]×[10000cos4θsin4θ00sin4θcos4θ00001][ax2+ay2ax2ay22axaycosωbt2axaysinωbt]s0s1s2s3=121cos(4θ+90°)sin(4θ+90°)0cos(4θ90°)cos2(4θ+90°)sin(4θ+90°)cos(4θ+90°)0sin(4θ+90°)sin(4θ+90°)cos(4θ+90°)sin2(4θ+90°)0000010000cos22ϕ+sin22ϕcosΔ(1cosΔ)sin2ϕcos2ϕsin2ϕsinΔ0(1cosΔ)sin2ϕcos2ϕsin22ϕ+cos22ϕcosΔcos2ϕsinΔ0sin2ϕsinΔcos2ϕsinΔcosΔ×10000cos4θsin4θ00sin4θcos4θ00001ax2+ay2ax2ay22axaycosωbt2axaysinωbt

(26.45)

Finally, the detected intensity I is given from the first term of s0

I=S0=ax2+ay22axaycos[ωbtΔcos(4θ2ϕ)]I=S0=ax2+ay22axaycos[ωbtΔcos(4θ2ϕ)]

(26.46)

where ωb is called the beat signal of the optical heterodyne method and Δ ≪ 1.

The optical heterodyne technique is one of the famous interferometers, as shown in Figure 26.16. The electric fields with ω1 and ω2 are shown in Figure 26.16a. The interference signal at the reference detector and measurement detector are shown in Figure 26.16b and c, respectively. The phase difference in Figure 26.16 means the phase Δ · cos(4θ – 2φ) of Equation 26.46. This retardance Δ and azimuth angle φ are easy to determine by using a discrete Fourier transform. The resolution of the birefringence measurement is 0.01 nm of the retardation at 0.2° of the azimuthal direction.

Image

FIGURE 26.16  Optical heterodyne detection. (a) Reference detector measures signal coupled from back of He–Ne laser through LP (45°), (b) from reference detector, and (c) from measurement detector.

Image

FIGURE 26.17  Stress mapping of CaF2. (a) CaF2, (b) retardation, and (c) azimuth direction.

Figure 26.17 is birefringence distribution of CaF2. The size of this sample is 70 mm in diameter. The distribution of the retardance is less than 25°.

26.8.6    TWO-DIMENSIONAL BIREFRINGENCE MEASUREMENT BY PHASE-SHIFTING METHOD

Figure 26.18 represents a birefringence distribution measurement by the phase-shifting technique [11,12]. Figure 26.18 is an optical arrangement. It consists of a quarter wave plate, a retarder that uses a Babinet-Soleil compensator, and a linear polarizer. A sample with Δ retardation and φ azimuthal direction is set after the retarder. The retarder can work to modulate both retardation δ and azimuthal direction θ.

The relation among the polarization elements is given using Mueller matrices and Stokes vectors by

S=P(θ+45°)X(Δ,φ)R(δ,θ)SS=P(θ+45°)X(Δ,φ)R(δ,θ)S

(26.47)

where

the Stokes vector S′ and S are the output and input of the light, respectively

P is Mueller matrix of the polarizer with azimuthal direction angle θ

X(Δ, φ) is a Mueller matrix of the sample with Δ, retardance and φ azimuthal direction

R(δ, θ) is a Muller matrix of the retarder

Image

FIGURE 26.18  Birefringence measurement by phase-shifting technique.

All Stokes parameters and Mueller matrices are calculated as

[s0s1s2s3]=12[1sin2θcos2θ0sin2θsin22θsin2θcos0cos2θsin2θcoscos22θ00000][1000010Δsin2ϕ001Δcos2ϕ0Δsin2ϕΔcos2ϕ1]×[10000cos22θ+sin22θcosδ(1cosδ)sin2θcossin2θsinδ0(1cosδ)sin2θcossin22θ+cos22θcosδcos2θsinδ0sin2θsinδcos2θsinδcosδ][1001]s0s1s2s3=121sin2θcos2θ0sin2θsin22θsin2θcos2θ0cos2θsin2θcos2θcos22θ000001000010Δsin2ϕ001Δcos2ϕ0Δsin2ϕΔcos2ϕ1×10000cos22θ+sin22θcosδ(1cosδ)sin2θcos2θsin2θsinδ0(1cosδ)sin2θcos2θsin22θ+cos22θcosδcos2θsinδ0sin2θsinδcos2θsinδcosδ1001

(26.48)

Finally, the detected intensity I is determined by the first term of s0

I=s0=12I0[1cos(tan1[Δcos(2θ2ϕ)]+δ)]I=s0=12I0[1cos(tan1[Δcos(2θ2ϕ)]+δ)]

(26.49)

where retardation Δ ≪ 1.

The intensity modulates sinusoidally by changing the phase δ of the retarder in Equation 26.49. If we apply the phase-shifting technique with π/2 of each phase shift shown in Figure 26.19 [13,14], we can determine the initial phase by

Δcos(2θM2ϕ)=I3I1I0I2=ΦMΔcos(2θM2ϕ)=I3I1I0I2=ΦM

(26.50)

The retardation Δ and azimuthal direction φ analytically substitute the phase ΦM corresponding to the azimuthal angle θ.

Image

FIGURE 26.19  Four steps of phase-shifting technique.

The birefringence phase difference is given by

Δ=12(Φ2Φ0)2+(Φ3Φ1)2Δ=12(Φ2Φ0)2+(Φ3Φ1)2

(26.51)

The azimuthal direction of

ϕ=12tan1Φ3Φ1Φ0Φ2ϕ=12tan1Φ3Φ1Φ0Φ2

(26.52)

This means 16 images are enough to determine birefringence mapping. The resolution of retardance is ±1°.

Figure 26.20 shows a birefringence distribution measurement using two LC retarders for highspeed measurement. Three LC retarders are required to determine all polarization states, generally. However, two retarders are enough to determine the limited retardation in this experiment. Two LC cells whose azimuthal directions are set at 0° and 45°, respectively, are used. If the incident light to the LC retarder is circular polarization, phase φ (ellipticity) and azimuthal direction θ by LC retarders are expressed as

δ1=sin1(sinϕsin2θ)δ2=tan1(tanϕcos2θ)δ1δ2=sin1(sinϕsin2θ)=tan1(tanϕcos2θ)

(26.53)

where δ1 and δ2 are retardation of the LC modulator. Figure 26.21 represents a two-dimensional accuracy check using a 15° retarder at three different orientations with 30°, 0°, and −30°. The dots indicate the measurement points, and the directions of the line indicate the azimuthal direction.

Now, a process to determine birefringence distribution of a plastic disk is shown in Figure 26.22. The original image in Figure 26.22a is the same as that of a polariscope. After analyzing the birefringence distribution using 16 images, the birefringence distribution such as retardation and azimuthal angle are shown in Figure 26.22b and c. A vector map shown in Figure 26.22d indicates both magnitude of retardance and azimuthal direction as a whole. A cross section of AA′ in Figure 26.22b is displayed in Figure 26.22e.

Image

FIGURE 26.20  Birefringence distribution measurement using two LC retarders.

Image

FIGURE 26.21  Two-dimensional accuracy check; same 15° retarder at three orientations. (a) 30° of azimuthal direction, (b) 0° of azimuthal direction, and (c) −30° of azimuthal direction.

Image

FIGURE 26.22  Process to determine birefringence distribution of plastic disk (3.5 in.). (a) Original image, (b) birefringence distribution, (c) azimuthal angle distribution, (d) vector map, and (e) cross section AA.

This birefringence distribution is mainly caused by polymer orientation. The residual stress mapping of the plastic disk is shown in Figure 26.23. We can compare the affective mapping in Figure 26.23a to the defective one in Figure 26.23b. Both cross sections are shown at the bottom. The birefringence distribution of the defective sample indicates the birefringence change in concentric circles and two parts of the radial direction. The part corresponds to the water channel for cooling water, indicating the residual stress caused by thermal stress.

Figure 26.24 demonstrates stress-induced birefringence by uniaxial stretching. A sample is a gelatin seat measuring 10 × 10 mm. An experimental setup for uniaxial stretching is shown in Figure 26.24a. The stress is indicated by a moving micrometer. Birefringence and azimuthal direction distributions along the length of extension are shown in Figure 26.24b and c, respectively. This relation is shown in Figure 26.25. There are not many changes until 50 ∝m because of the bio-sample [15].

Image

FIGURE 26.23  Plastic disk inspection. (a) Good product and (b) defective product.

Image

FIGURE 26.24  Change of birefringence under uniaxial stretching. (a) Experimental setup, (b) birefringence distribution, and (c) azimuth distribution.

Image

FIGURE 26.25  Change of birefringence under uniaxial stretching.

REFERENCES

1.   P.S. Theocaris and E.E. Gdoutos: Matrix Theory of Photoelasticity, Springer, Berlin, Germany, pp. 105–131, 1979.

2.  R.M.A. Azzam and N.M. Bashara: Ellipsometry and Polarized Light, North-Holland, New York, pp. 153–268, 1987.

3.  D. Goldstein: Polarized Light, Marcel Dekker, New York, pp. 553–623, 2003.

4.  W.A. Shurclif: Polarization Light, Harvard University Press, Cambridge, U.K., pp. 124–148, 1962.

5.  H. Aben and C. Guillement: Photoelasticity of Glass, Springer, Berlin, Germany, pp. 51–77, 1993.

6.  R.A. Chipman: Polarimetry, Handbook of Optics II, Michael Bass (Editor in Chief) McGraw-Hill, New York, Chapter 22, pp. 22.1–22.37, 1995.

7.  Y. Shindo and H. Hanabusa: An improved highly sensitive instrument for measuring optical birefringence, Polym. Commun., 25, 378–382, 1984.

8.  B. Wang: Linear birefringence measurement at 157 nm, Rev. Sci. Instrum., 74, 1386, 2003.

9.  B. Wang and T.C. Oakberg: A new instrument for measuring both the magnitude and angle of low level linear birefringence, Rev. Sci. Instrum., 70, 3847, 1999.

10.  N. Umeda, S. Wakayama, S. Arakawa, A. Takayanagi, and H. Kohwa: Fast birefringence measurement using right and left hand circularly polarized laser, Proc. SPIE, 2873, 119–122, 1996.

11.  Y. Otani, T. Shimada, T. Yoshizawa, and N. Umeda: Two-dimensional birefringence measurement using the phase shifting technique, Opt. Eng., 33(5), 1604–1609, 1994.

12.  Y. Otani, T. Shimada, and T. Yoshizawa: The local-sampling phase shifting technique for precise two-dimensional birefringence measurement, Opt. Rev., 1(1), 103–106, 1994.

13.  J.E. Greivenkamp and J.H. Bruning: Phase shifting interferometry, Optical Shop Testing, Daniel Malacara, (Ed.), John Wiley, New York, Chapter 14, pp. 501–598, 1992.

14.  K. Creath: Temporal phase measurement methods, Interferogram Analysis, David W. Robinson and T. Reid (Eds.), Institute of Physics Publishing, London, U.K., Chapter 4, pp. 94–140, 1993.

15.  M. Ebisawa, Y. Otani, and N. Umeda: Microscopic birefringence imaging by the local-sampling phase shifting technique, Proc. SPIE, 5462, pp. 45–50, (Biophotonics Micro and Nano Imaging, Strasbourg, France), 2004.

16.  Luceo Co., Ltd, Circular polarimetry, http://www.luceo.co.jp/technical/pdf/enhenko.pdf (accessed on January 5, 2015).

..................Content has been hidden....................

You can't read the all page of ebook, please click here login for view all page.
Reset
18.118.16.81