A polarimetric experimental method was developed to determine the Jones matrix elements of transparent optical materials without sign ambiguity. A set of polarization dependent transmittance data of the samples was measured with polarizer – sample – analyzer system and another set of data was measured with polarizer – sample – quarterwave plate – analyzer. Two data sets were compared and mathematically analyzed to obtain the correct signs of the elements of the matrix. The Jones matrix elements of a quarterwave plate were determined to check the validity of the method. The experimentally obtained matrix elements of the quarterwave plate were consistent with the theoretical expectations. The same method was applied to obtain the Jones matrix elements of a twisted nematic liquid crystal panel.
I. INTRODUCTION
Liquid crystal panels (LCPs) are optoelectronic devices that can modulate the state of light. They are used for many practical applications including flat panel display, holography and optical information processing because of their inherent ability to spatially modulate a light beam in a programmable manner
[1

4]
. Twisted nematic (TN) type LCP spatial light modulators (SLMs) are good for intensity modulation and parallel aligned LCP SLMs are suitable for phase modulation. Some application areas of SLMs such as ultrashort laser pulse control require phase modulation
[5
,
6]
. However parallel aligned LCP SLMs are expensive, while TNLCP SLMs are relatively lowcost and easily available. This is why much effort has been made to use TNLCP SLMs for many applications
[7
,
8]
.
For programmable operation of a TNLCP spatial light modulator (SLM), the Jones matrix of the TNLCP needs to be determined accurately. However, manufacturers do not usually provide that information, and the user must determine the matrix. According to the literature, the TNLCP can be considered as a stack of thin slices acting as homogeneous uniaxial crystals
[9]
. The Jones matrix of a homogeneous uniaxial crystal slice can be written as a function of its ordinary and extraordinary indices of refraction as well as the orientation angle of the liquid crystal molecules
[9]
. The Jones matrix of the whole TNLCP can be described as the product of the individual Jones matrices of each homogeneous uniaxial crystal slice
[10]
. According to the approximation, the Jones matrix of the whole TNLCP can be written as
Here
c
_{1}
is a constant used to define the intensity loss,
i
is
and
Δϕ
is the phase change induced by the LCP at a given wavelength. The matrix describes the change in polarization state of the incoming beam. The goal of this study is to determine the signs as well as the absolute values of the matrix elements can be determined by polarimetric methods experimentally, but it has been thought that the signs of the elements can not be determined uniquely by polarimetric methods only. However it has been shown that the signs can be determined by polarimetric methods
[11]
.
In this paper, we demonstrated a polarimetric method for unambiguous determination of signs of the Jones matrix elements. By measuring the transmittance with an optical system of polarizerLCPanalyzer, we obtained the set of the unitary matrix elements. In addition, by measuring the transmittance after inserting a quarter wave plate in the system, we obtained another set of unitary matrix elements. By comparing two sets of unitary matrix elements, we could determine the unitary matrix elements accurately without the sign ambiguity.
Our method has some advantages with respect to other previous similar methods. First, our method improved the accuracy. In particular, polarizer and analyzer rotate in the same and opposite directions with 1degree angular increment from 0° to 180°, so Jones matrix elements were obtained from fitting a trigonometric function with many data points. Also, our method uses only two rotating elements (polarizer and analyzer) and except for the rotation of the polarizers the alignments of all the other optical elements were not disturbed during the measurement, which can reduce potential errors due to realignment.
II. THEORY
 2.1. Transmittance Intensity of the First Part of the Optical System (Quarter Wave Plate + Polarizer)
The polarimetric method to determine the Jones matrix is based on the transmittance measurement of the optical system. Let us first explain how to calculate the transmittance of the optical system. Our whole optical system consists of a HeNe laser light source, a quarter wave plate (QWP1), polarizer (P), TNLCP, analyzer (A) and a silicon photodiode detector as shown in
Fig. 1
.
The optical system used to measure the transmittance of a TNLCP. Azimuth angles of the polarizer and analyzer are ψ_{P} and ψ_{A}, respectively. BE: beam expander, QWP1 and QWP3: quarter wave plates, P: polarizer, LCP: liquid crystal panel, A: analyzer, L: lens, I: Iris, D: photodiode.
First we considered the transmittance intensity of the optical system consists of a light source, a quarterwave plate (QWP1) and a polarizer. The laser light source was linearly polarized, so the QWP1 was used to make the linearly polarized light into the circularly polarized light. The circularly polarized light was the input beam of the polarizer. The azimuth angle of the polarizer was defined as
ψ_{P}
. The relation between the Jones vector of the input beam (
a
,
b
‧
e
^{iδ}
) and that of the output beam (
x
_{out1}
,
y
_{out1}
) for the optical system can be written as Eq. (2)
Here
a
and
b
are the electric field components on the vertical (
x
,
z
) and horizontal (
y
,
z
) planes, respectively.
δ
is a phase difference between the electric field components.
Q
(
π
/4) is the Jones matrix of the QWP1 whose fast axis was aligned 45° with respect to the vertical (
x
) axis.
P
(
ψ_{P}
) is the Jones matrix of the polarizer whose azimuthal angle is
ψ_{P}
. The output beam intensity was determined by following equation,
The laser light was almost linearly polarized in the vertical (
x
,
z
) plane, so the electric field components satisfied a condition,
a
>>
b
.
 2.2. Transmittance Intensity of the Whole Optical System (Quarter Wave Plate + Polarizer + TNLCP + Analyzer)
We calculated the transmittance intensity of the whole optical system consists of a light source, a quarterwave plate (QWP1), a polarizer, a TNLCP and an analyzer. The TNLCP was placed between a polarizer and an analyzer whose azimuth angles were
ψ_{P}
and
ψ_{A}
, respectively. The output beam of the polarizer described in Eq. (2) was used as the input beam of the later optical system. The Jones vector of the output beam for the whole optical system can be written as Eq. (4).
 2.3. Definition of the Transmittance of the Optical System Composed of TNLCP and Analyzer
We defined the transmittance of the optical system composed of TNLCP and analyzer,
T_{ψP}
,
_{ψA}
, as the relative ratio between the output (I) beam intensity of Eq. (5) and the input beam intensity (I
_{1}
) of Eq. (3).
T_{ψP}
,
_{ψA}
can be expressed as Eq. (7),
Here
c
_{3}
is a constant.
III. EXPERIMENT
The optical system used in our study was shown in
Fig. 1
. The light source was a linearly polarized HeNe laser. The laser beam was expanded by using the beam expander (BE) and then was sent to the quarter wave plate (QWP1). The quarterwave plate made the linearly polarized light circularly polarized. The circularly polarized beam was directed to a rotating linear polarizer (P) and then was sent to the sample. The polarization state of the transmitted beam was investigated with an analyzer (A) and an additional quarter wave plate (QWP3). The intensity of the beam passed through the analyzer was measured with the photodiode (D). A quarterwave plate (QWP2) and a TNLCP SLM (Holoeye LC2002, SONY LCX016AL6) were used as samples.
The transmittance of the optical system (sample + analyzer) depends only on the sum (
ψ_{A}
+
ψ_{P}
) and the difference (
ψ_{A}

ψ_{P}
) of the angles of the polarizer and analyzer. After substituting
φ
=
ψ_{A}

ψ_{P}
and
ψ
=
ψ_{A}
+
ψ_{P}
into Eq. (7), the transmittance of the optical system,
T_{φ, ψ}
takes a simpler form.
Here C is a constant.
φ
_{0}
,
ψ
_{0}
are the initial phases of the trigonometric functions and A, B are the amplitudes of the functions. A set of transmittance
T_{φ, ψ}
data was collected by changing
φ
and
ψ
from 0° to 360°. The coefficients A and B, also the initial phase angles
φ
_{0}
and
ψ
_{0}
in Eq. (8) were determined by the NewtonRaphson fitting method.
IV. RESULTS AND DISCUSSION
 4.1. Determination of the Jones Matrix Elements of a Quarterwave Plate
To check the validity of our method, the Jones matrix elements of a quarter wave plate (QWP2) is experimentally determined and compared with the theoretical one. The fast axis of the QWP2 is set to make 45º with the
x
axis. If the quarterwave plate is ideal, then the Jones matrix should be
from Eq. (1), where
c
_{4}
is the constant representing the intensity loss. The matrix should be unitary, so the elements
f
,
h
,
g
and
j
should satisfy the condition,
f
^{2}
+
h
^{2}
+
g
^{2}
+
j
^{2}
= 1. To satisfy the constraints
f
,
j
,
h
and
g
should be
h
=
g
= 0. Eqs. (6), (7), and (8) show that
f
,
j
,
h
and
g
are connected to the coefficients A and B also the initial phase angles
φ
_{0}
and
ψ
_{0}
of the trigonometric functions representing the transmittance of the optical system, i.e.
The transmittance
T_{φ, ψ}
of the quarterwave plate (QWP2) and the analyzer was measured and fitted with Eq. (8) to determine A, B,
φ
_{0}
and
ψ
_{0}
. To make Eq. (8) simpler, some conditions are chosen to make the cosines of the sum and the difference of the angles constant. First the analyzer and polarizer rotated in the opposite direction with same angle increment, so that the ‘sum angle (
ψ_{A}
+
ψ_{P}
)’ is set to be zero, i.e.
ψ
=
ψ_{A}
+
ψ_{P}
=0. Then the transmittance as a function of the difference angle
φ
=
ψ_{A}

ψ_{P}
is measured as presented in
Fig. 2
(a). By fitting the data points with a
Transmittance (T) of the optical system composed of a quarterwave plate (QWP2), and an analyzer. Instead of the LCP shown in FIG. 1, a quarter wave plate (QWP2) was inserted. (a) The analyzer and the polarizer rotated in the opposite directions. (b) The analyzer and polarizer rotated together. The solid circles are the experimental result, and the line is the theoretical estimation.
cosine function, the information about
f
and
h
is obtained.
Similarly, the analyzer and the polarizer are rotated in the same direction with the same angle increment so that the ‘difference angle (
ψ_{A}

ψ_{P}
)’ is set to be zero, i.e.
φ
=
ψ_{A}

ψ_{P}
=0. Then the transmittance as a function of the sum angle
ψ
=
ψ_{A}
+
ψ_{P}
is measured as presented in
Fig. 2
(b). By fitting the data points with a cosine function, the information about
g
and
j
is obtained.
It should be noted however that, the relative signs of
g
and
j
elements to those of
f
or
h
were not uniquely defined because of the ambiguity of the phase of the trigonometric functions
φ
_{0}
and
ψ
_{0}
in Eq. (10).
To eliminate the sign ambiguity, an additional quarterwave plate (QWP3) was inserted between the quarterwave plate (QWP2) and the analyzer. The fast axis of QWP3 was aligned 45° with respect to the x axis. Including the contribution of the QWP3, the transmittance became a function of
φ
and
ψ
as shown in Eq. (11).
Here
c
_{5}
is a constant.
m
,
q
,
p
and
l
were connected to
f
,
h
,
g
and
j
by comparing Eq. (11) with Eq. (6). As a result,
Also the sum of the squares of the new elements should
Transmittance (T) of the optical system composed of a quarterwave plate (QWP2), additional quarter wave plate (QWP3), and an analyzer. Instead of LCP shown in FIG. 1, a quarter wave plate (QWP2) was inserted. (a) The analyzer and the polarizer rotated in the opposite directions. (b) The analyzer and the polarizer rotated together. The solid circles are the experimental result, and the line is theoretical estimation.
satisfy
m
^{2}
+
q
^{2}
+
p
^{2}
+
l
^{2}
= 2 because of the unitarity.
To determine
m
,
p
,
q
and
l
, first the analyzer and the polarizer are rotated in the opposite directions with the same angle increment so that the sum angle remains zero, i.e.
ψ
=
ψ_{A}
+
ψ_{P}
=0. The transmittance as a function of the difference angle
φ
=
ψ_{A}

ψ_{P}
is presented in
Fig. 3
(a). By data fitting, the information about
m
and
q
is obtained. Similarly, the analyzer and the polarizer are rotated together so that the difference angle remains zero, i.e.
φ
=
ψ_{A}

ψ_{P}
=0. Then the transmittance as a function of the sum angle
ψ
=
ψ_{A}
+
ψ_{P}
is presented in
Fig. 3
(b). By data fitting, the information about
p
and
l
is obtained.
The sign of the
f
,
g
,
j
, and
h
can be determined by comparing four element values of the first method (without QWP3) with four element values of the second method (with QWP3). The experimentally obtained Jones matrix elements of QWP2 are
f
= 0.706,
j
= 0.707,
h
= 0.017 and
g
= 0.026 consistent with the theory, which confirms the validity of our method.
 4.2. Determination of the Jones Matrix Elements of a TNLCP Spatial Light Modulator
The same method was used to determine the Jones matrix elements of a programmable TNLCP SLM. A programmable TNLCP SLM was installed between the polarizer and analyzer as shown in
Fig. 1
. The TNLCP can change the intensity and the phase of the input beam as the grey level is controlled by an external device. The corrected transmittances
T
_{0}
of the TNLCP and the analyzer for several different grey levels were measured. The meaning of the corrected transmittance
T
_{0}
will be discussed in subsection 4.3.
Fig. 4
shows the difference angle (
φ
=
ψ_{A}

ψ_{P}
) and the sum angle (
ψ
=
ψ_{A}
+
ψ_{P}
) dependence of the
T
_{0}
for
The corrected transmittance (T_{0}) of the optical system composed of a TNLCP, and an analyzer as shown in FIG. 1. (a) The analyzer and polarizer rotated in the opposite directions. (b) The analyzer and the polarizer rotated together. The solid circles are the corrected experimental result, and the line is the theoretical estimation.
The corrected transmittance (T_{0}) of the optical system composed of a TNLCP, a quarterwave plate (QWP3), and an analyzer as shown in FIG. 1. (a) The analyzer and the polarizer rotated in the opposite directions. (b) The analyzer and the polarizer rotated together. The solid circles are the corrected experimental result, and the line is the theoretical estimation.
one of the grey levels (grey level 192 of our LCP device). As demonstrated for the quarterwave plate (QWP2) case, the absolute values of the
f
,
g
,
j
, and
h
of the Jones matrix were determined by fitting the transmittance with Eq. (8).
To determine the signs of the Jones matrix elements, additional set of data were measured.
Fig. 5
shows the sum angle (
ψ
=
ψ_{A}
+
ψ_{P}
) and the difference angle (
φ
=
ψ_{A}

ψ_{P}
) dependence of
T
_{0}
. Note that additional quarterwave plate QWP3 was inserted between the TNLCP and the analyzer for obtaining the data set. The absolute values of the
m
,
q
,
p
, and
l
of the Jones matrix were determined by fitting the transmittance with the Eq. (11).
The sign of the
f
,
g
,
j
, and
h
can be determined by comparing four element values of the first method (without QWP3) with four element values of the second method (with QWP3). The experimentally obtained Jones matrix elements of TNLCP (at grey level 192 of our LCP device) were
f
= 0.64,
h
= 0.63,
g
= 0.44 and
j
= 0.02. The matrix elements of other grey levels can be obtained similarly.
 4.3. Corrections for the Transmittance Measurement
There is discrepancy between the theoretical transmittances described by Eqs. (8) and (11) and the measured transmittance. Such discrepancy is negligible for the quarterwave (QWP2), but is not negligible for the TNLCP.
Fig. 6
shows the measured transmittance (
T
, solid asterisks) of TNLCP, the errorcorrected transmittance (
T
_{0}
, solid circles), and the theoretical fitting function (solid line). They do not overlap perfectly and the discrepancy is especially large around the
π
/4 of the difference angle.
The intensity loss occurred at the TNLCP is the most
Uncorrected and corrected transmittance of the optical system composed of a TNLCP, and an analyzer as shown in FIG. 1. The solid asterisks are the transmittance without correction (T). The analyzer and polarizer rotated in the opposite directions. The solid circles are the transmittance (T_{0}) of the optical system with correction, and the line is the theoretical estimation.
Transmittance (T_{1}) of the TNLCP for (a) the clockwise polarizer rotation and (b) for the counterclockwise polarizer rotation as illustrated in the insets. ‘D’ represents the detector.
dominant reason of the discrepancy. The intensity loss can be caused by many factors including absorption, scattering or reflection. However, the intensity loss does not affect the fitting process significantly in our method.
Figure 7
shows the polarizer angle dependence of the transmittance of the TNLCP (
T
_{1}
). Note that
T
_{1}
was measured without the analyzer and the additional quarterwave plate QWP3. There are two polarizer rotation directions.
Fig. 7
(a) shows the
T
_{1}
for the clockwise polarizer rotation (See the inset of
Fig. 7
(a)). This corresponds to the polarize rotation for the difference angle (
φ
=
ψ_{A}

ψ_{P}
) measurement case shown in
Figs. 4
(a) and
5
(a).
Fig. 7
(b)
T
_{1}
for the counterclockwise polarizer rotation (See the inset of
Fig. 7
(b)). This corresponds to the polarize rotation for the sum angle (
ψ
=
ψ_{A}
+
ψ_{P}
) measurement case shown in
Figs. 4
(b) and
5
(b). The
T
_{1}
data show the oscillating behavior deviated from the constant value about 2%.
The final transmittance (
T
) of the whole optical system can be considered as the multiplication of the
T
_{1}
and the transmittance (
T
_{0}
) of the analyzer (or analyzer + the additional quarterwave plate QWP3, when QWP3 is inserted). To remove the polarizer dependence in the final transmittance of the whole optical system the transmittance
T
was divided with transmittance (
T
_{1}
) of TNLCP. This is the corrected transmittance (
T
_{0}
=
T
/
T
_{1}
) of the optical system also the transmittance (
T
_{0}
) of the analyzer (or analyzer + the additional quarterwave plate QWP3). Errors were negligible for the quarterwave plate QWP2 data, so the Jones matrix elements of QWP2 are obtained without the correction.
V. CONCLUSION
A polarimetric method was developed for determining the Jones matrix element without the sign ambiguity. The credibility of the method was checked by comparing the theoretical and experimental Jones matrix elements of a quarterwave plate. The method was used to determine the Jones matrix elements of a twisted nematic liquid crystal panel (TNLCP) successfully. The proposed method can be used to determine the Jones matrix elements of any other transparent optical materials.
Acknowledgements
This research was supported by the National Research Foundation of Korea (Grants No. 20100023535 and No 220 20111C00016).
Castellano J. A.
(2006)
“Modifying light”
American Scientist
94
438 
445
Barnes T. H.
,
Matsuda K.
,
Eiju T.
,
Matsumoto K.
,
Johnson F.
(1990)
“Joint transform correlator using a phase only spatial light modulator”
Jpn. J. Appl. Phys.
29
1293 
1296
Mok F.
,
Diep J.
,
Liu H.
,
Psaltis D.
(1986)
“Realtime computergenerated hologram by means of liquidcrystal television spatial light modulator”
Opt. Lett.
11
748 
750
Ostrovsky A. S.
,
MartínezNiconoff G.
,
Arrizón V.
,
MartínezVara P.
,
OlveraSantamaría M. A.
,
RickenstorffParrao C.
(2009)
“Modulation of coherence and polarization using liquid crystal spatial light modulators”
Opt. Express
17
5257 
5264
Wang Z.
,
Marks D. L.
,
Carney P. S.
,
Millet L. J.
,
Gillette M. U.
,
Mihi A.
,
Braun P. V.
,
Shen Z.
,
Prasanth S. G.
,
Popescu G.
(2011)
“Spatial light interference tomography”
Opt. Express
19
19907 
19918
Polachek L.
,
Oron D.
,
Silberberg Y.
(2006)
“Full control of the spectral polarization of ultrashort pulses”
Opt. Lett.
31
631 
633
Yocky D. A.
,
Bernes T. H.
,
Matsumoto K.
,
Ooyama N.
,
Matsuda K.
(1990)
“Simple measurement of the phase modulation capability of liquid crystal phaseonly light modulators”
Optik
84
140 
144
Ohkubo K.
,
Ohtsubo J.
(1993)
“Evaluation of LCTV as a spatial light modulator”
Opt. Commun.
102
116 
124
Lu K.
,
Saleh B. E. A.
(1990)
“Theory and design of the liquid crystal TV as an optical phase modulator”
Opt. Eng.
29
240 
246
Yariv A.
,
Yeh P.
1984
Optical Waves in Crystal
Wiley
New York, USA
Chapter 5.
Moreno I.
,
Velásquez P.
,
FernándezPousa C. R.
,
SánchezLópez M. M.
,
Mateos F.
(2003)
“Jones matrix method for predicting and optimizing the optical modulation properties of a liquidcrystal display”
J. Appl. Phys.
94
3697 