Achromatic Retarders in Polarization Optics: Michael Kraemer Tom Baur
Achromatic Retarders in Polarization Optics: Michael Kraemer Tom Baur
Achromatic Retarders in Polarization Optics: Michael Kraemer Tom Baur
optics
Michael Kraemer
Tom Baur
Michael Kraemer, Tom Baur, “Achromatic retarders in polarization optics,” Opt. Eng. 58(8), 082406 (2019),
doi: 10.1117/1.OE.58.8.082406.
Abstract. Retarders or waveplates are tools for polarization modification in bulk optical systems. These devices
usually have a strong wavelength dependence in their performance, making them suitable for use over a wave-
length band on the order of a few percent of the center wavelength for which they are made. Display and tunable
laser applications are examples that can require consistent polarization modification over a much broader wave-
length range. We discuss new methods and designs for dramatically increasing range of performance and
review older methods as well. We show examples of achievable performance using modern polymer and liquid
crystal materials. © 2019 Society of Photo-Optical Instrumentation Engineers (SPIE) [DOI: 10.1117/1.OE.58.8.082406]
Keywords: ferroelectric liquid crystals; achromatic; retarder; waveplate; liquid crystal; polarization; superachromat.
Paper 181372SS received Sep. 24, 2018; accepted for publication Dec. 26, 2018; published online Jan. 24, 2019.
Fig. 1 Two separate TIRs in the Fresnel rhomb each contribute a 45-
deg phase change in the p and s components of an incoming light ray,
thereby producing a quarter-wave retarder. Fresnel invented this
prism in 1817 to convert incoming linearly polarized light into outgoing Fig. 3 For internal reflection angle α ¼ 55° 5 0 , the retardance versus
circularly polarized light. wavelength curve for BK7 glass yields the smallest deviation from
90 deg, namely a deviation of 1°31′ over the 400- to 1600-nm wave-
length range, taking dispersion of refractive index according to
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi Schott’s data sheets8 into account.
δðλÞ cos α nðλÞ2 · sin2 α − 1
tan ¼− ; (1)
nðλÞ · sin2 α
EQ-TARGET;temp:intralink-;e001;63;548
2
• influence of stress birefringence in the rhomb
for all incidence angles α that obey the TIR condition: material,4
• influence of surface nonuniformities and of surface
1 films,10
arcsin ≤ α ≤ 90 deg : (2)
• lateral displacement of the beam during passage
EQ-TARGET;temp:intralink-;e002;63;499
nðλÞ
through the device (for quarter wave rhomb).11
In the Fresnel rhomb, two phase shifts of 45 deg each are
introduced by two separate TIRs, totaling up to the desired Additional surface films can be deposited on the rhomb
90-deg phase shift of a quarter wave retarder. In Fig. 2, two surface to lower the phase retardation maximum for a single
quarter wave Fresnel rhombs are connected to form a half TIR in Fig. 4 to 45 deg, or to increase retardance achroma-
wave Fresnel rhomb, that is one that shifts the phase of ticity or to decrease the retardance dependence on the angle
an incoming ray by 180 deg. For BK7 glass with a refractive of incidence, see Refs. 4 and 12 and references therein.
index of 1.5151 at 632.8 nm, for example, α ¼ 55° 5 0 produ-
ces the necessary 45-deg phase shift. Figure 3 shows retard-
ance versus wavelength dependence for a 400- to 1600-nm 3.2 N-Crystals Combination Plates: Achromatic
Fresnel rhomb made from Schott N-BK 7® glass, using the Retardance through Dispersion of Birefringence
dispersion formula from Schott.8 The disadvantages of the In birefringent crystals such as quartz or magnesium fluo-
Fresnel rhomb are ride, the birefringence is not constant as a function of wave-
length (which would give rise to a retardance inversely
• size, as compared with fractions of a mm for, e.g., proportional to wavelength), but depends on the wavelength,
polymer stack retarders, too. This phenomenon, which is called dispersion of birefrin-
• sensitive dependence of the resulting phase shift on the gence, can be utilized to construct achromatic retarders by
angle of incidence, which restricts the angular aperture stacking two or more waveplates of different birefringent
considerably,9
• rotation of retarder fast axis with the azimuth angle for
non-normal incidence,
Fig. 2 Doubling the quarter wave Fresnel rhomb from Fig. 1 yields a
half wave retarder with four TIRs, each contributing 45 deg to the Fig. 4 Dependence of rhomb retardance on internal reflection angle
retardance. Note that for the half wave rhomb the outgoing beam for BK7 glass at 632.8 nm (n ¼ 1.5151). Best for increased angular
is collinear with the incoming one (no displacement). For further aperture is TIR retardance graph being tangent to dashed 45-deg
rhomb geometries, see Ref. 7. retardance line.4
materials with fast axes either parallel or perpendicular to from 0.25 wavelengths (i.e., quarter wave retardance) in the
each other. The advantages of such a crystal combination range from 400 to 1000 nm. An even better fit for the same
plate are wavelength range, yielding a quarter wave achromat with
λ∕100 tolerance, is a combination of a 104.6-μm-thick mag-
• The stack of crystal plates is typically a few millimeters nesium fluoride MgF2 ) waveplate and a 24.1-μm-thick
thick, less bulky than a prism retarder. ammonium dihydrogen phosphate (ADP) waveplate with
• The direction of the stack optic axis does not vary with parallel optic axes, as shown in the same Fig. 5.
wavelength, as is the case for the Pancharatnam-type However, crystal plates of a thickness of a few hundred
plate combinations described in Sec. 3.3.2. (There microns or less are very fragile and therefore difficult to pol-
are other types of crystal waveplate combinations ish to a given retardance goal. Instead, it is necessary to sub-
(see Ref. 13) where fast or slow axes are not aligned; stitute the 104.6-μm-thick MgF2 waveplate by a pair of two
these show achromaticity over a shorter wavelength thicker MgF2 waveplates, e.g., with thicknesses 2104.6 and
range, however, and experience a variation of the stack 2000 μm, in a manner that their optic axes are perpendicular
fast axis with wavelength.) to each other so that their retardances are subtracted and yield
• Crystal waveplates have a higher threshold to with- a net retardance of the desired 2104.6 − 2000 ¼ 104.6 μm.
stand laser or solar flux than polymer retarders. Similarly, the unrealistically thin 24.1-μm ADP waveplate
• Achromats of arbitrary target retardances can be real- can be replaced by a combination of two thicker ADP plates
ized by increasing or decreasing the thicknesses of all in subtraction mode (i.e., with perpendicular optic axes), e.g.,
plates in a given achromatic crystal stack by the same by a 2024.1 μm and a 2000-μm ADP plate with net retard-
scaling factor. ance 2024.1 − 2000 ¼ 24.1 μm. This method is also appli-
cable to realize small net retardances in the Pancharatnam-
The retardance (in units of wavelength) of a stack of N type combination retarder stacks that will be discussed in
crystal plates at wavelength λ (in microns), assuming normal Sec. 3.3.6 and in the application example in Sec. 4.1.
incidence of the ray on the plates, is given as According to Beckers,14 the substitution of crystal wave-
plates of thickness D by two waveplates of the same material
X
N
ϵ i − ωi in subtraction mode (i.e., with fast axes crossed with thick-
RðλÞ ¼ di ; (3) nesses D1 and D2 such that D ¼ D1 − D2 ) not only allows
λ
EQ-TARGET;temp:intralink-;e003;63;484
X
N
εi − ωi ϕ2 cos2 α sin2 α
RðλÞ ¼ di 1 − −
λ ωi εi
EQ-TARGET;temp:intralink-;e005;326;268
i¼0
2ωi
þ Oðϕ4 Þ: (5)
Fig. 9 Measured retardance versus wavelength dependence for a quarter wave five-layer stretched
polymer superachromat for 380 to 1100 nm (measured data points as circles) that was manufactured
in October 2016 at Meadowlark Optics, compared with the theoretical curve obtained by modeling with
Mueller matrices, taking dispersion of birefringence into account.
where the first and third waveplates are of equal retardance sin f · π2
2δ1 with fast axes parallel, whereas the center plate has sinðf · 2δ1 Þ ¼ sin 2δ1 ; (8)
cos Δ
EQ-TARGET;temp:intralink-;e008;326;507
EQ-TARGET;temp:intralink-;e006;63;353
QW achromat design
400 to 800 nm
(1 octave) const. No. of Design Retardance FA
birefringence plates parameters deviation variation
2δ1
such as polyvinyl alcohol20 possess a nearly constant bire-
fringence for visible light, in general the dispersion of
By comparison of Eq. (8) with Eq. (10) and of Eq. (9) birefringence may lead to less achievable achromaticity in
with Eq. (11), we see that in the limit f → 0 Pancharatnam’s a given wavelength range and must be taken into account
equations go over into Title’s equations, so that Title’s ach- when designing achromats. With the following mathematical
romat can be considered a limit case of the Pancharatnam transformation, all achromat models may be generalized to
achromat for small wavelength ranges. Figure 10 shows
the stack retardance and stack fast axis angle graphs for
Table 2 Retardance deviation (in degrees) and fast axis variation for
both approaches for the fit with least retardance deviation the three-half wave achromat designs shown in Fig. 11 for a 400- to
from 0.25 waves (quarter wave achromat) over the wave- 800-nm wavelength range (i.e., one octave).
length range. As a third solution, we show the one obtained
by numerical optimization for a five-plate-stack designed by
the authors. HW achromat design
400 to 800 nm
(1 octave) const. No. of Design Retardance FA
birefringence plates parameters deviation variation
3.3.4 N-layer designs
Title, conditions at λ0 3 λ0 ¼ 533.3 nm 14°22′ 5°27′
By stacking an odd number of waveplates symmetric about
the center plate (as for plate retardances and plate fast axis Pancharatnam, 3 λ0 ¼ 533.3 nm, 4°25′ 7°8′
orientations), a tighter tolerance can be realized for quarter conditions at ð1 f Þλ0 f ¼ 0.2851
wave, half wave, and other achromats. The authors have Meadowlark, global 5 proprietary 0°16′ 1°41′
min. of retardance
designed and manufactured five-plate polymer film achro- deviation
mats for quarter wave, third wave, and half wave require-
ments. In Ref. 20, Title proposed and provided error Title, 3 × 3 9 λ0 ¼ 533.3 nm, 0°13′26″ 0°23′5″
analysis for a half wave achromat obtained by stacking combination retardance
three of his three-plate achromats with parallel fast axes 92°13′
for the outer three-stack elements and with fast axis of the
Pancharatnam, 9 λ0 ¼ 533.3 nm, 0°0′23″ 0°2′12″
center three-stack rotated by 60 deg with respect to the 3 × 3 combination f ¼ 0.2851
fast axis of the outer three-stack elements. This nine-stack
Table 3 Retardance deviation (in degrees) and fast axis variation for
the three half wave achromat designs shown in Fig. 12 for a 400- to
1600-nm wavelength range (i.e., two octaves). The nine-stack
consisting of three Pancharatnam three-stacks provides the least
retardance deviation over the two octave range, while the stacking
of three title three-stacks yields better achromaticity near the design
wavelength.
HW achromat design
400 to 1600 nm
(two octaves) const. No. of Design Retardance FA
birefringence plates parameters deviation variation
λ0 bðλÞ
δðλÞ ¼ δðλ0 Þ · · : (12)
λ bðλ0 Þ
EQ-TARGET;temp:intralink-;e012;326;518
λ bðλ0 Þ
λ0 ¼
EQ-TARGET;temp:intralink-;e013;326;365
bðλÞ
¼λ· ; (13)
bðλÞ
bðλ0 Þ
λ0
δ 0 ðλ 0 Þ ¼ δðλ0 Þ · ; (14)
λ0
EQ-TARGET;temp:intralink-;e014;326;293
have a birefringence increasing with wavelength; an example Here, c is a compensation parameter to be chosen so that
is apophyllite,7,23,24 a mineral with a “natural” retardance the retardance versus wavelength function of the pair is
achromaticity in the visible, which is, however, “rare in usa- close to constant, thus making the pair as retardance achro-
ble optical quality specimens.”23 Certain stretched polymer matic as possible. If we then symmetrically stack three
films have increasing birefringence, also known as reverse such achromatic waveplate pairs of suitable retardances
or abnormal wavelength dispersion and can in some cases to form a 3 × 2 Pancharatnam combination waveplate, the
be used as single-layer achromatic retarder.25,26 achromaticity is far better than for the simple Pancharatnam
As an example, we will use Pancharatnam’s design for a three-stack of material A. As an example, we want to build
quarter wave retarder using three single magnesium fluoride a mixed dispersion 3 × 2 Pancharatnam-type achromat for
(MgF2 ) plates in the wavelength interval from 2 to 5 μm. the wavelength range 400 to 1600 nm (two octaves). With
According to Ref. 27, a best fit to measured data for MgF2 design wavelength λ0 ¼ 1000 nm, we transform the wave-
birefringence as function of wavelength λ (in microns) is length λ to the wavelength λ 0 in which the birefringence
obtained in the model equation: dispersion appears constant:
Iλ2 Jλ2 λ
bðλÞ ¼ H þ þ ; (16) λ0 ¼ : (17)
λ2 − G λ2 − L bðλÞ
c þ 1 − c bðλ
EQ-TARGET;temp:intralink-;e017;326;609
EQ-TARGET;temp:intralink-;e016;63;609
with coefficients H ¼ −19.364 · 10−3 , I ¼ 30.992 · 10−3 , For a particular choice of polymers for materials A and B,
G ¼ 2.3253 · 10−3 , J ¼ 40.060 · 10−3 , and L ¼ 388.37. 0 ∕λ 0
the wavelength interval ratio λmax min obtains its minimum
With λmin ¼ 2, λmax ¼ 5, and the design wavelength at approximately c ¼ 2.76, yielding the wavelength trans-
λ0 ¼ 3.5, we transform to zero-dispersion wavelength formation depicted in Fig. 15 with values between λmin 0
¼
bðλ0 Þ
λ 0 in the range from λmin
0
¼ λmin · bðλ min Þ
0
¼ 1.84 to λmax ¼ 0
755.5 nm and λmax ¼ 1461.6 nm. So we see that the origi-
bðλ0 Þ nal two octaves of 400 to 1600 nm have shrunk down to
λmax · bðλ max Þ
¼ 5.82 (see Fig. 13), and, choosing f ¼ 0.29
less than one octave in terms of zero-dispersion wavelength,
in Pancharatnam’s method, we obtain the retardance and
see Fig. 16. The “doubling” of the shape is a result of the fact
fast axis graphs shown in Fig. 14 for zero-dispersion
that the wavelength transformation is not one-to-one, see
wavelength λ 0 and for actual wavelength λ.
Fig. 15.
3.3.6 Mixed dispersion Pancharatnam-type six-plate
combinations 3.3.7 Effect of oblique incidence on retardance and
Instead of a single waveplate of material A with birefringence fast axis direction of a Pancharatnam-type
function bðλÞ and retardance (in waves) 2δ at the design combination plate
wavelength λ0 , one can consider a pair (of same resulting Consider a stack of N parallel waveplates of the same
retardance 2δ) consisting of a waveplate of a zero-dispersion uniaxially birefringent material with positive birefringence
material B and retardance ðc þ 1Þ · 2δ in subtraction mode βðλÞ that is assumed to be very small compared with the
with a waveplate of material A with retardance c · 2δ. ordinary refractive index (small-birefringence approxima-
tion).28 Let the normal-incidence plate retardances be
r1 ðλÞ; r2 ðλÞ; : : : ; rN ðλÞ at a certain wavelength, and let
the fast axes be in the surface plane oriented at angles
α1 ; α2 ; : : : ; αN relative to a reference axis in the surface
plane. Since the birefringence is much smaller than the ordi-
nary refractive index, the ordinary and extraordinary ray are
parallel in zero order approximation, so that a linearly polar-
ized beam that is obliquely incident on the waveplate stack
will transverse the stack in a straight line with an internal
incidence angle ϕ determined from the external incidence
angle with Snell’s law and with an azimuth angle α to
the reference axis in the surface plane that is unaffected
by refraction. According to Hale and Day,28 the phase
delay of the extraordinary wave relative to the ordinary
wave in plate no. i in the small-birefringence aproximation
is given as
Fig. 14 In terms of the zero-dispersion wavelength λ 0 , the “classic” Pancharatnam three-stack solution
(i.e., assuming constant birefringence) in the range from 1.84 to 5.82 μm is equivalent to the original
problem for λ between 2 and 5 μm.
EQ-TARGET;temp:intralink-;e021;326;391 cos θi ¼ k^ · Ω
^ i ¼ sin ϕ sinðα − αi Þ: (21)
1 − sin2 ϕ sin2 ðα − αi Þ
ri ¼ ri · : (22)
cos ϕ
EQ-TARGET;temp:intralink-;e022;326;338
cos ϕ
tion in the plane normal to the beam. We suggest to define ωi
as the polar angle in the s–p coordinate system in the beam
the optic axis of plate no. i is given by the unit vector: normal plane, where the unit vector s^ ¼ ð− sin α; cos α; 0ÞT
! is both normal to the plane of incidence and normal to
− sin αi the beam, and the unit vector p^ ≔ k^ × s^ is in the plane of
^i ¼
Ω
EQ-TARGET;temp:intralink-;e020;63;112 cos αi ; (20) incidence. With this convention, the waveplate no. i under
0 oblique incidence is equivalent to a waveplate normal to
Fig. 16 In the transformed wavelengths λ 0 , the original two-octave (400 to 1600 nm) Pancharatnam prob-
lem is reduced to one with less than an octave of wavelength (from 755.5 to 1461.6 nm) if the wavelength
transformation to zero-birefringence wavelength is not one-to-one, see Fig. 15.
the beam with the normal-incidence retardance ri in Eq. (22) variation of the stack fast axis over a given wavelength
and the fast axis orientation βi given as range, depend on incidence and azimuth angles of the
beam relative to the compound retarder. For compound
half wave retarders with three layers (by Pancharatnam19),
EQ-TARGET;temp:intralink-;e024;63;359
Fig. 18 For ideal (modeled) Pancharatnam-type half wave achromats in a 400- to 800-nm wavelength
ranges of 3, 5, 7, and 9 layers, the dependence of retardance deviation and fast axis variation (i.e., differ-
ence of maximum and minimum compound fast axis angle in the wavelength range) are compared, for
each incidence angle taking the maximum (worst case) over all wavelengths and all azimuth angles.
Table 4 Design retardance and fast axis orientation for the three
pairs of Pancharatnam elements that form the SAR (in the 2 to
5 μm wavelength range) for calibration of the Cryo-NIRSP of the
DKIST, currently being built in Hawai’i.
Fig. 22 Modeled retardance graph (left) and fast axis graph (right) for γ 1 R1 þ γ 2 R2 ¼ 0;
EQ-TARGET;temp:intralink-;e027;326;422 (27)
open and closed state of half wave retarder between crossed
polarizers. we ensure that the change of retardance of material 1 due to
a temperature gradient is compensated by an equal, but of
opposite sign, change of retardance of material 2. The con-
Following Hariharan,31 this issue can be relieved by
dition that the two retardances R1 and R2 must add up to a
making the FLC half wave retarder achromatic by placing
given target retardance yields a second equation to calculate
it between two nonswitchable half wave polymer retarders the two unknown plate retardances R1 and R2 . Achromaticity
with parallel fast axes to form a Pancharatnam-type half (independence of retardance from wavelength) and athermal-
wave achromat, see Fig. 21 for operation between crossed ity (independence of retardance from temperature) cannot, in
Fig. 23 Comparison of modeled closed state (a) and open state (b) polarized transmission intensities for
achromatic Pancharatnam-stack (solid lines) and nonachromatic single (dashed lines) FLC light valves
between crossed polarizers. Since for the closed state, the fast axis of the single-layer FLC retarder is
aligned with the axes of the crossed polarizers, the closed state transmission is zero over the entire
wavelength range (for perfect polarizers).
Fig. 24 Measured open and closed state transmissions and contrast ratio for an achromatic FLC device
built by Meadowlark Optics compared with the measured contrast ratio of a single layer FLC device.
The fringing oscillations seen in the open state transmission graph are due to reflections from the
glass windows that encase the FLC center waveplate.
22. F. Pockels, Lehrbuch der Kristalloptik, Teubner, Leipzig and Berlin 31. P. Hariharan and P. Ciddor, “Achromatic switchable polarization
(1906). rotators,” Opt. Eng. 36, 952–956 (1997).
23. A. M. Title and W. J. Rosenberg, “Achromatic retardation plates,” Proc. 32. P. Hariharan, P. Ciddor, and M. Roy, “Improved switchable achromatic
SPIE 307, 120–125 (1981). phase shifters, part 2,” Opt. Eng. 44, 105603 (2005).
24. M. Françon, S. Mallick, and J. Vulmière, “Use of apophyllite in polari- 33. S. Guimond and D. Elmore, “Polarizing views,” oemagazine (2004),
zation interferometers—achromatic half-wave plate,” J. Opt. Soc. Am. http://spie.org/newsroom/polarizing-views?SSO=1 (03 January 2017).
55, 1553–1553 (1965).
25. A. Uchiyama and T. Yatabe, “P-5: characteristics and applications of
new wide-band retardation films,” in SID 01 Digest, pp. 566–569 Michael Kraemer received his PhD in mathematics from the
(2001). University of Montana–Missoula in 2001 and his PhD in physics from
26. C. D. West and A. S. Makas, “The spectral dispersion of birefringence, Colorado State University in 2006. His current research interests
especially of birefringent plastic sheets,” J. Opt. Soc. Am. 39, 791–794 include numerical computations in polarization optics, optical meas-
(1949). uring techniques, and achromatic devices.
27. S. R. Sueoka, “Polarization optical components of the Daniel K. Inouye
Solar Telescope,” PhD thesis, University of Arizona (2016). Tom Baur is the founder and CTO of Meadowlark Optics. He has
28. P. D. Hale and G. W. Day, “Stability of birefringent linear retarders been developing polarization-related products for the company for
(waveplates),” Appl. Opt. 27, 5146–5153 (1988).
29. D. F. Elmore et al., “The Daniel K. Inouye Solar Telescope first light almost 40 years. He is the 2017 winner of the SPIE G.G. Stokes
instruments and critical science plan,” Proc. SPIE 9147, 9147071 Award and has a dual career in precision optics and cattle ranching
(2014). on the eastern plains of Colorado.
30. W. H. Schubert, E. Petrak, and T. G. Baur, “Measurement of polariza-
tion assemblies for the Daniel K. Inouye Solar Telescope,” Proc. SPIE
9369, 93690N (2015).