Simulations of birefringent gratings as polarizing color separator in backlight
for flat-panel displays

doi:10.1364/OE.15.005789

« Show journal navigation

Simulations of birefringent gratings as polarizing color separator in backlight for flat-panel displays

Man Xu, H. Paul Urbach, and Dick K. G. de Boer »View Author Affiliations

Optics Express, Vol. 15, Issue 9, pp. 5789-5800 (2007)
http://dx.doi.org/10.1364/OE.15.005789

View Full Text Article

Acrobat PDF (416 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. Configuration
3. Theory
4. Simulations and discussions
5. Conclusions
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (14)
Cited By
Figures (8)
Metrics

Abstract

A color and polarization separating backlight can be obtained by using a surface-relief grating made of birefringent material as an outcoupling structure on top of the lightguide. A rigorous finite element diffraction model was applied to study the polarization effect of such a grating. The diffraction of plane waves by the anisotropic grating was studied for general conical incidence.

1. Introduction

Liquid crystal displays (LCDs) have been adapted to a broad range of applications as convenient flat panel displays. But the light efficiency of the current LCD systems is below 10%. The low efficiency is mainly because of the use of the absorption-based polarizers and color filters. The absorption by the polarizers and color filter is approximately 50% and 70%, respectively. Therefore, alternative approaches for polarization- and color-separation have to be considered. Ideally one would like to use one device which can separate both color and polarization and which is not based on absorption. By applying a diffraction grating on the lightguide in the LCDs, colors can be separated[1

Y. Taira, D. Nakano, H. Numata, A. Nishikai, S. Ono, F. Yamada, M. Suzuki, M. Noguchi, R. Singh, and E.G. Colgan, “Low-power LCD using a novel optical system,” SID 02 Digest , 1313–1315 (2002). [CrossRef]

][2

F. Yamada, S. Ono, and Y. Taira, “Dual layered very thin flat surface micro prism array directly molded in an LCD cell,” Eurodisplay 2002 , 339–342 (2002).

]. If the grating is made of anisotropic material, one can achieve that light of one polarization is diffracted and light of the orthogonal polarization is not[3

Dick K. G. de Boer, Roberto Caputo, Hugo J. Cornelissen, Chris M. van Heesch, Eefje J. Hornix, and Martin J.J. Jak, “Diffractive grating structures for colour-separating backlights,” Photonics in Multimedia, Proc. SPIE 6196, (2006).

]. The configuration of such a backlight for LCDs is shown in Fig. 1. This structure was placed in a side-lit configuration with a cold-cathode fluorescence lamp (CCFL). The grooves of the grating are filled with a birefringent material having its optical axis parallel to the grooves. The ordinary refractive index of the birefringent material is matched to the grating material, whereas its extraordinary refractive index is significantly higher. As a result, light with polarization parallel to the grooves is diffracted, whereas for the relevant angles of incidence, the light with polarization perpendicular to the grooves is totally reflected back into the lightguide.

Fig. 1. Configuration of the polarized color-separating backlight.

The principle and performance of both the non-polarized and polarized version of color-separating backlights have been discussed in a previous paper [3

]. In a color-separating backlight the light distribution, which contains a multitude of angles and mainly three wavelengths, is diffracted into three slightly overlapping angular distributions for the main colors. We showed experimental results of the (polarized) color-separated angular distribution of the light extracted from a structure like that of Fig. 1. We also showed simulations using a commercially available computer program [4

www.gsolver.com

], based on rigorous coupled wave theory for non-birefringent materials. We were able to make a good comparison between measured and simulated data in the plane perpendicular to the grooves. However, it was not possible to do simulations for other (conical) directions. Note that in general a backlight contains light of several directions and all these should be taken into account to understand its efficiency.

The aim of this paper is to discuss a method for simulating diffraction by birefringent gratings at conical angles, applied to the problem of polarized color-separating backlights. To study the performance of the diffraction gratings containing birefringent material, a numerical program is required. We applied a numerical code based on finite element method (FEM) in the simulations [5

X. Wei, H.P. Urbach, and A.J.H. Wachters, “Finite Element Model for Three-Dimensional Optical Scattering Problems,” J. Opt. Soc. Am. A , 24, 866 (2007). [CrossRef]

]. With our program any arbitrary conical angle of incidence can be calculated rigorously. Results of simulations are shown in this paper for when the incident wave vectors are in the plane perpendicular to the grooves (normal incidence) and for wave vectors with conical incidence.

2. Configuration

The example and numerical values that we use are similar to those of Ref.[3

]. The geometry and configuration of the grating in our study is depicted in Fig. 2. A surface-relief grating with a period of 400 nm and duty cycle of 0.5 is applied onto a lightguide made of polycarbonate. The grooves of the grating are filled with a birefringent material which is liquid crystal (LC) having its optical axis along the grooves. The wavelength dependent refractive indices are shown in Table.1. The grating is rectangular and the depth of the grooves is 140 nm.

Fig. 2. Configuration of the grating with the computational domain Ω.

Table 1. The refractive indices of polycarbonate and of liquid crystal for three colors.

λ	Polycarbonate n	LC n_o	LC n_e
450	1.61	1.54	1.82
535	1.58	1.53	1.77
632	1.58	1.52	1.76

3. Theory

3.1. Anisotropic media that are translation invariant with respect to one direction

For a configuration that consists of isotropic materials and that is translational invariant with respect to one coordinate, say the y-coordinate, it is well-known that Maxwell’s equations are equivalent to a system of two coupled second order partial differential equations for only the components E_y and H_y [6

Max Born and Emil Wolf, “Rigorous diffraction theory,” in Principles of Optics, (The University Press, Cambridge, 2005), pp.633–673.

]. This is also true when the sources and fields depend on the y-coordinate by a factor e^ik_yy , for some k_y . When k_y = 0, the two partial differential equations for E_y and H_y are uncoupled and therefore one can separate s- and p-polarization. For s-polarization, H_y vanishes whereas for p-polarization E_y is zero. When k_y ≠ 0, the two partial differential equations are coupled and the polarizations mix.

These observation apply in particular to periodic diffraction gratings that consist of isotropic materials. In this section we want to investigate to what extent the separation in two polarizations is still possible in the case of gratings made of anisotropic materials. We will fist show that in an anisotropic medium of which the dielectric permittivity ε͇ is independent of the y-direction, all field components can be expressed in terms of E_y and H_y .

The time-harmonic source-free Maxwell equations in non-magnetic anisotropic materials are,

\nabla \times E = iω μ_{0} H,

(1)

\nabla \times H = - iω ε_{0} ε͇ E,

(2)

where ε ₀ and μ ₀ are the dielectric permittivity and the magnetic permeability of vacuum respectively, and ε͇ is the relative dielectric permittivity of the anisotropic crystal which is a complex symmetric tensor of rank 2 given by:

ε͇ = (\begin{matrix} ε_{xx} & ε_{xy} & ε_{xz} \\ ε_{yx} & ε_{yy} & ε_{yz} \\ ε_{zx} & ε_{zy} & ε_{zz} \end{matrix}),

(3)

with respect to the cartesian coordinate system (x, y, z). We write the tensor ε͇ as a sum of its real and imaginary parts,

ε_{ij} = ε_{ij}^{'} + {iε}_{ij}^{″}, (i, j = x, y, z)

(4)

The real and imaginary parts of ε͇ commute, which means that for any vector v we have,

ε′ ε″ v = ε″ ε′v .

(5)

This implies that the tensors ε′ and ε″ have the same eigenvectors and that therefore they can be diagonalized on the same (orthogonal) basis.

The configuration and the material properties are assumed to be invariant with respect to translations in the y-direction, but may vary in a cross-sectional plane y = constant. Hence ε͇ may be a function of (x, z): ε͇(x, z).

We consider an electromagnetic field which depends harmonically on y, i.e.

\frac{\partial}{\partial y} E = {ik}_{y} E,

(6)

\frac{\partial}{\partial y} H = {ik}_{y} H,

(7)

for some real k_y . It is easy to see that by using Maxwell’s equations, the x- and z-components of the electromagnetic field can be expressed in terms of the y-components:

(\binom{E_{x}}{E_{z}}) = - \frac{ω^{2} ε_{0} μ_{0}}{D} 𝒩 (\binom{ε_{xy}}{ε_{zy}}) E_{y} + \frac{i}{D} 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{{ωμ}_{0} H_{y}}) + \frac{i}{D} 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{{ωμ}_{0} H_{y}})

(8)

(\binom{- H_{z}}{H_{x}}) = \frac{k_{y}}{{ωμ}_{0}} (\binom{E_{x}}{E_{z}}) + \frac{i}{{ωμ}_{0}} (\binom{\frac{\partial}{\partial x}}{\frac{\partial}{\partial z}}) E_{y}

= - \frac{ω^{2} ε_{0} μ_{0}}{D} \frac{k_{y}}{{ωμ}_{0}} 𝒩 (\binom{ε_{xy}}{ε_{zy}}) E_{y} + \frac{i}{D} \frac{k_{y}}{ω μ_{0}} 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{{ωμ}_{0} H_{y}}) + \frac{i}{D} \frac{k_{y}}{ω μ_{0}} 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}}) + \frac{i}{ω μ_{0}} (\binom{\frac{\partial}{\partial x}}{\frac{\partial}{\partial y}}) E_{y},

(9)

where

D = (ω^{2} ε_{0} μ_{0} ε_{xx} - k_{y}^{2}) (ω^{2} ε_{0} μ_{0} ε_{zz} - k_{y}^{2}) - {(ω^{2} ε_{0} μ_{0})}^{2} ε_{xz} ε_{zx}

= {(ω^{2} ε_{0} μ_{0})}^{2} (ε_{xx} ε_{zz} - ε_{xz} ε_{zx}) - ω^{2} ε_{0} μ_{0} k_{y}^{2} (ε_{xx} + ε_{zz}) + k_{y}^{4} .

(10)

𝒩 = (\begin{matrix} ω^{2} ε_{0} μ_{0} ε_{zz} - k_{y}^{2} & - ω^{2} ε_{0} μ_{0} ε_{xz} \\ - ω^{2} ε_{0} μ_{0} ε_{zx} & ω^{2} ε_{0} μ_{0} ε_{xx} - k_{y}^{2} \end{matrix}),

(11)

𝒬 = 𝒩 (\begin{matrix} 0 & - 1 \\ 1 & 0 \end{matrix}) = (\begin{matrix} {- ω}^{2} ε_{0} μ_{0} ε_{xz} & - (ω^{2} ε_{0} μ_{0} ε_{zz} - k_{y}^{2}) \\ ω^{2} ε_{0} μ_{0} ε_{xx} - k_{y}^{2} & ω^{2} ε_{0} μ_{0} ε_{zx} \end{matrix}),

(12)

Note that D is the determinant of matrix Q;. It is assumed here that D ≠ 0. The y-components of Maxwell’s equations (1) and (2) imply:

- iω μ_{0} H_{y} = \overset{͂}{\nabla} ∙ (\begin{matrix} 0 & 1 \\ - 1 & 0 \end{matrix}) (\binom{E_{x}}{E_{z}}),

(13)

iω ε_{0} ε_{yy} E_{y} = - iω ε_{0} (ε_{yx} ε_{yz}) (\binom{E_{x}}{E_{z}}) - \overset{͂}{\nabla} ∙ (\binom{H_{z}}{{- H}_{x}}),

(14)

where ∇͂ is defined by,

\overset{͂}{\nabla} = (\binom{\frac{\partial}{\partial x}}{\frac{\partial}{\partial z}}) .

(15)

By substitution of Eqs. (8), (9) into the right-hand sides of Eqs. (13), (14) we get a coupled system for E_y and H_y ,

- iω μ_{0} H_{y} = - ω^{2} ε_{0} μ_{o} \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒬^{T} (\binom{ε_{xy}}{ε_{zy}}) E_{y}] + i \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒬^{T} \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω_{0} μ_{0} H_{y}})] + i \overset{͂}{\nabla} ∙ [\frac{1}{D} ℳ \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})],

(16)

i ω ε_{0} ε_{yy} E_{y} = i \frac{ω^{2} ε_{0} μ_{0}}{D} ω ε_{0} (ε_{yx} ε_{yz}) 𝒩 (\binom{ε_{xy}}{ε_{zy}}) E_{y} + \frac{ω ε_{0}}{D} (ε_{yx} ε_{yz}) 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})

+ \frac{ω ε_{0}}{D} (ε_{yx} ε_{yz}) 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}}) + \frac{k_{y}}{ω μ_{0}} ω^{2} ε_{0} μ_{0} \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒩 (\binom{ε_{xy}}{ε_{zy}}) E_{y}]

- i \frac{k_{y}}{ω μ_{0}} \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})] - i \frac{k_{y}}{ω μ_{0}} \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})] - \frac{i}{ω μ_{0}} Δ E_{y},

(17)

where

ℳ = (\begin{matrix} ω^{2} ε_{0} μ_{0} ε_{xx} - k_{y}^{2} & ω^{2} ε_{0} μ_{0} ε_{xz} \\ ω^{2} ε_{0} μ_{0} ε_{zx} & ω^{2} ε_{0} μ_{0} ε_{zz} - k_{y}^{2} \end{matrix}) .

(18)

Note that M = DN^-1. Even when k_y = 0, both Eqs. (16) and (17) contain E_y and H_y . Consider the case of a grating made of anisotropic materials which is periodic with respect to the x and invariant with respect to translations in the y-direction. If the material above the grating is isotropic and the incident wave is in the (x, z)-plane perpendicular to the grooves, i.e. when k_y = 0, then when the incident field is s- or p-polarized, the total field will be a mixture of both polarizations. In this respect the anisotropic and isotropic gratings totally differ.

Now we consider an anisotropic medium for which the y-axis (i.e. the axis of translational invariance) is principal axis of ε͇ for all points (x, y, z). Then the dielectric tensor is of the form,

ε͇ = (\begin{matrix} ε_{xx} & 0 & ε_{xz} \\ 0 & ε_{yy} & 0 \\ ε_{zx} & 0 & ε_{zz} \end{matrix}),

(19)

and the coupled system Eqs. (16) and (17) becomes:

- i ω μ_{0} H_{y} = i \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒬^{T} \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω_{0} μ_{0} H_{y}})] + i \overset{͂}{\nabla} ∙ [\frac{1}{D} ℳ \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})],

(20)

i ω ε_{0} ε_{yy} E_{y} = - i \frac{k_{y}}{ω μ_{0}} \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})] - i \frac{k_{y}}{ω μ_{0}} \overset{͂}{\nabla} ∙ [\frac{1}{D} 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})] - \frac{i}{ω μ_{0}} Δ E_{y},

(21)

Furthermore Eqs. (8) and (9) become

(\binom{E_{x}}{E_{z}}) = \frac{i}{D} 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}}) + \frac{i}{D} 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}})

(22)

(\binom{{- H}_{z}}{H_{x}}) = \frac{i}{D} \frac{k_{y}}{ω μ_{0}} 𝒩 \frac{\partial}{\partial x} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}}) + \frac{i}{D} \frac{k_{y}}{ω μ_{0}} 𝒬 \frac{\partial}{\partial z} (\binom{k_{y} E_{y}}{ω μ_{0} H_{y}}) + \frac{i}{ω μ_{0}} (\binom{\frac{\partial}{\partial x}}{\frac{\partial}{\partial z}}) E_{y},

(23)

Now assume again that k_y = 0, then Eq. (20) contains only H_y and Eq. (21) contains only E_y :

ω^{2} μ_{0} ε_{0} H_{y} + \frac{\partial}{\partial x} (\frac{ε_{xx}}{ε_{xx} ε_{zz} - ε_{xz} ε_{zx}} \frac{\partial H_{y}}{\partial x} + \frac{ε_{zx}}{ε_{xx} ε_{zz} - ε_{xz} ε_{zx}} \frac{\partial H_{y}}{\partial z}) + \frac{\partial}{\partial z} (\frac{ε_{xz}}{ε_{xx} ε_{zz} - ε_{xz} ε_{zx}} \frac{\partial H_{y}}{\partial x} + \frac{ε_{zz}}{ε_{xx} ε_{zz} - ε_{xz} ε_{zx}} \frac{\partial H_{y}}{\partial z}) = 0,

(24)

ω^{2} ε_{0} μ_{0} ε_{yy} E_{y} + \frac{\partial^{2} E_{y}}{{\partial z}^{2}} + \frac{\partial^{2} E_{y}}{{\partial x}^{2}} = 0 .

(25)

Then the system is uncoupled for E_y and H_y . Eqs. (24) and (25) imply that when k_y = 0 we can distinguish two types of polarization, namely:

S-polarization H_y = 0, H_x ≠ 0, H_z ≠ 0, E_x = E_z = 0, E_y ≠ 0;

P-polarization E_y = 0, E_x ≠ 0, E_z ≠ 0, H_x = H_z = 0, H_y ≠ 0.

Hence when the y-axis is the principal axis, the splitting in uncoupled s- and p-polarization for k_y = 0 is analogous to the isotropic case. But, E_x and E_z have other values than in the isotropic case because Eq. (22) depends on off-diagonal elements of the tensor.

In our specific case, the LC is uniaxial material. For such uniaxial materials, of which the y-axis coincides with the optical axis, the dielectric permittivity tensor is diagonal:

ε͇ = (\begin{matrix} ε_{xx} & 0 & 0 \\ 0 & ε_{yy} & 0 \\ 0 & 0 & ε_{zz} \end{matrix}),

(26)

with ε_xx = ε_zz = n ² _o and ε_yy = n ² _e. When k_y = 0 the s- and p-polarization can be separated. We can choose n = n _o for p-polarization and n = n_e for s-polarization for analysis.

In general when k_y ≠ 0, the system is coupled as can be easily seen from Eqs. (16) and (17).

3.2. Numerical Method

In a general three-dimensional anisotropic grating problem all the electric and magnetic field components are coupled when k_y ≠ 0. Several authors have studied the scattering of an incident plane wave by an anisotropic grating. The Coupled Wave Method for isotropic gratings has been extended to anisotropic gratings by Rokushima et al. [7

K. Rokushima and J. Yamakita, “Analysis of anisotropic dielectric gratings,” J. Opt. Soc. Am. A , 73, 901 (1983). [CrossRef]

], Glytsis et al. [8

E.N. Glytsis and T.K. Gaylord, “Rigorous three-dimensional coupled-wave diffraction analysis of single and cascaded anisotropic gratings,” J. Opt. Soc. Am. A , 4, 2061 (1987). [CrossRef]

] and Mori et al. [9

S. Mori, K. Mukai, J. Yamakita, and K. Rokushima, “Analysis of dielectric lamellar gratings coated with anisotropic layers,” J. Opt. Soc. Am. A , 7, 1661 (1990). [CrossRef]

]. The Chandezon method (or C-method) has been applied to anisotropic gratings by Harris et al. [10

J.B. Harris, T.W. Preist, E.L. Wood, and J.R. Sambles, “Conical diffraction from multicoated gratings containing uniaxial materials,” J. Opt. Soc. Am. A , 13, 803 (1996). [CrossRef]

]. Li [11

L. Li, “Reformulation of the Fourier modal method for surface-relief gratings made with anisotropic materials,” J. Mod. Opt. , 45, 1313 (1998). [CrossRef]

] has extended his Fourier modal method (FMM) to anisotropic gratings. The FMM is a modification of the coupled wave method to improve convergence for the case of TM polarization. In Li’s paper conical incident angles are allowed.

In this paper we apply the Finite Element Method (FEM) [12

Xiuhong Wei , Three Dimensional Rigorous Model for Optical Scattering Problems , PhD thesis, Optics Research Group, Delft University of Technology, August 2006.

], [5

X. Wei, H.P. Urbach, and A.J.H. Wachters, “Finite Element Model for Three-Dimensional Optical Scattering Problems,” J. Opt. Soc. Am. A , 24, 866 (2007). [CrossRef]

] to model anisotropic gratings. The FEM is a general numerical method for solving boundary value problems in mathematical physics. When applied to Maxwell’s equations for 2D gratings, a 2D computational domain Ω (Fig.2) in a plane perpendicular to the grooves is defined which has width in the x-direction of one period and which contains all nonplanar and anisotropic materials in its interior. The domain is truncated in the vertical direction by a so-called Perfect Matched Layer (PML)[13

J.P. Berenger, “Perfectly matched layer for the absorption of electromagnetic waves,” Journal of Computational Physics , 114(2), 185–200 (1994). [CrossRef]

]. On this computational domain, a boundary value problem for the vector Helmholtz equation for either the electric or the magnetic field is formulated. (Note that we do not solve the coupled system for the E_y - and H_y -components that was derived in the previous section). The computational domain is meshed using triangles or quadrilaterals and the electric of magnetic field is approximated by edge elements.

The FEM is relatively difficult to implement, especially in 3D. However, it is very flexible because both periodic and non-periodic configurations and all kinds of materials can be modeled. In particular, anisotropic and inhomogeneous materials with arbitrary principal axes can be treated without problems. In this respect the FEM differs substantially from the Couple Wave Method and the FMM for which the extension to anisotropic materials is highly non-trivial.

To ensure that the FEM code gives accurate results when applied to the anisotropic gratings considered in this paper, we have compared the results obtained for a flat multi-layer consisting of anisotropic media with an analytic model based on the expansion in terms of plane waves. Such a multi-layer is a very good test case for the FEM code because, in contrast with for example the FMM, the FEM does not simplify when applied to an anisotropic multi-layer. We found that when the edge elements of lowest order are used whith 50 mesh points per wavelength, the computed field is quite accurate for all conical angles of incidence.

3.3. Birefringent gratings

The periodic grating which is applied onto the lightguide, was designed in such a way that for all angles of incidence of interest, only the -1^st order is transmitted and propagates. If we consider a plane wave incident from the lightguide on to the grating with wave vector kⁱ and incident angles θⁱ and ϕⁱ, the tangential components of the wave vector are given by,

k_{x}^{i} = k_{0} n^{i} \cos ϕ^{i} \sin θ^{i},

(27)

k_{y}^{i} = k_{0} n^{i} \sin ϕ^{i} \sin θ^{i},

(28)

where k ₀ = 2π/λ is the wave number in vacuum and nⁱ is the refractive index in the medium of the incident wave. Polar angle θ and azimuthal angle ϕ are defined and depicted in Fig. 2. The wave vector of a diffracted transmitted wave can be expressed as,

k_{x}^{d} = k_{0} n^{d} \cos ϕ^{d} \sin θ^{d},

(29)

k_{y}^{d} = k_{0} n^{d} \sin ϕ^{d} \sin θ^{d},

(30)

where n^d is the refractive index of the medium into which the light is partially transmitted and ϕ^d and θ^d are diffraction angles. The components of the wave vector of the m^th transmitted order and the components of the incident wave vector are related by:

k_{x}^{d} = k_{x}^{i} + \frac{2 π ∙ m}{p},

(31)

k_{y}^{d} = k_{y}^{i},

(32)

where p is the period or pitch of the grating. The angular directions of the m^th diffracted order can be calculated for varying incident angles θⁱ and ϕⁱ from the following equations:

n^{d} \cos ϕ^{d} \sin θ^{d} = n^{i} \cos ϕ^{i} \sin θ^{i} + \frac{mλ}{p},

(33)

n^{d} {\sin ϕ}^{d} \sin θ^{d} = n^{i} \sin ϕ^{i} \sin θ^{i} .

(34)

In our grating the pitch p = 400 nm and therefore for most wavelengths and for the angles of incidence 48° ≤ θⁱ ≤ 90° and 0° ≤ ϕⁱ ≤ 90° only the - 1^st order is propagating. For blue light with short wavelength the -2^nd transmitted order is also propagating, but only for very large incident angles.

4. Simulations and discussions

4.1. Normal incidence

As we discussed in Section 3.1, when the principal axis of the anisotropic material is along the grooves and the plane of incidence is perpendicular to the grooves, we can distinguish s- and p-polarized waves. In this case, the polarization is conserved. When the incident field is s-polarized, then one transmitted diffracted order propagates which is s-polarized as well. According to our definition of angles, ϕ^d = ϕⁱ = 0 in this case. We will call this case normal incidence.

Fig. 3. Relative intensity of the -1^st diffracted transmitted order as function of θⁱ for ϕⁱ = 0 and for three wavelengths namely 450, 535 and 632 nm for (a) s-polarization and (b) p-polarization.

The dependence of the intensity of the -1^st diffracted transmitted order as function of incident angle θⁱ is shown in Fig. 3 for both s- and p-polarization. The transmission of s- polarization is much larger than that of p-polarization. The small nonzero intensity of the transmitted order for p-polarization is due to the small mismatch between the refractive index of the polycarbonate and the ordinary refractive index of the LC. It can be seen that the first transmitted order appears above a certain incident angle. The 0^th transmitted orders become evanescent above 40°, and the intensity of the -1^st transmitted order then changes abruptly. The field increases as the angle becomes larger. For a side-lit lightguide system, the main rays have incident angles θ; larger than 48°. For the angles of interest, the transmission of s-polarization is high comparing to that of p-polarization. But the absolute transmission for s-polarization is still fairly small. This is not a problem in our application. Because the lightguide is very long and therefore the light that reflected back into the lightguide is recycled many times.

4.2. General incidence

In a more general case, when the plane of incidence is not perpendicular to the grooves, the diffracted transmitted field of the s- and p-polarized incident fields couple. Hence the transmitted field does not have the same polarization as the incident field anymore. We did simulations for both s- and p-polarized incident fields. The diffracted transmitted field is decomposed into two linear polarized components. One is along the direction of the grooves, which we call s-polarization. The other is perpendicular to the grooves, which is called p-polarization. Please notice that the definition of the transmitted s- and p-polarizations differs from that of the incident polarizations. The definition of the transmitted polarization agrees with what is customary [3

Fig. 4. Relative intensity of the -1^st diffracted transmitted order as function of incident angle ϕⁱ for incident angle ϕⁱ = 67° and for the three wavelengths 450, 535 and 632 nm. RI in the titles is an abbreviation of relative intensity, and the first subscript denotes the polarization of the diffracted transmitted order whereas the second subscript indicates the polarization of the incident field.

In a backlight configuration as that of Fig. 1 light is present with a multitude of angles (Lambertian). Typically, the polar angles θⁱ range from 48° to 90°. The distribution in diffracted angles θ^d has a typical width of 20°. The maximum of the distribution for diffracted green light is close to θⁱ = 67°. With fixed polar angle θⁱ = 67°, the azimuthal angle ϕⁱ was varied from 0° to 90°. The calculated intensities of the -1^st transmitted order are shown in Fig. 4 for the three wavelengths 450, 535 and 632 nm.

From those figures we can see that for θⁱ = 67° the diffracted order becomes evanescent when the incident angle ϕⁱ is larger than approximately 40°. The transmitted field strongly depends on the angle. The transmission of s-polarization is relatively high at small incident angle ϕⁱ, whereas the transmission of p-polarization is low.

The angular distribution of the wave components of the diffracted transmitted field was calculated according to the Eqs. (31) and (32), and is illustrated in Fig. 5. To give insight in the angular distribution of the diffracted light, the relative intensity of the -1^st diffracted transmitted order from Fig. 4 is plotted in Fig. 6 as function of the diffraction angle θ^d (instead of angle ϕⁱ as is done in Fig. 4).

Fig. 5. Angular distribution (θ^d , ϕ^d) of the -1^st diffracted order for different colors, blue (450 nm), green (535 nm) and red (632 nm) for θⁱ = 67° and for varying ϕⁱ (-90° ≤ϕⁱ ≤90°).

The transmission of s-polarization is much higher than that of p-polarization. For s-polarized incident light, especially at small diffraction angle θ^d , the large intensity difference between transmitted s-component and p-component results in a high contrast ratio, as is shown in Fig. 7(a). The contrast ratio for unpolarized light is somewhat lower, as shown in Fig. 7(b), due to the diffracted orders of p-polarization. By using materials with better matching refractive indices for the p-polarization, a higher contrast ratio can be achieved. Our measurements of the angular distribution of extracted s- and p-polarized light [3

], repeated in Fig.8, are in agreement with these simulated results. In practice the contrast ratio is not higher than 23, presumably due to non-perfect alignment of the liquid crystal in the grating.

Fig. 6. Relative intensity of the -1st diffracted transmitted order as function of diffraction angle θ^d , for incident angle θⁱ = 67°, and for varying 0° ≤ ϕⁱ ≤ 90°, for the three wavelengths 450, 535 and 632 nm. RI in the titles is an abbreviation for relative intensity, and the first subscript denotes the polarization of the transmitted field whereas the second subscript indicates the polarization of the incident field.

5. Conclusions

Diffractive outcoupling can be used to make efficient backlights. A diffraction grating containing birefringent material was studied with a rigorous model. The simulation results showed that the diffracted transmitted intensity strongly depends on incident angles and has a high contrast ratio for all angles of interest. The simulation results are in good agreement with experiments. In the case presented here, the efficiencies are not very high (around 2%). Note that these efficiencies refer to one interaction with the grating. For a long lightguide the total extraction efficiency will be much higher. Depending on the application, with a long or short lightguide, a higher extraction efficiency may be desired. This can be achieved by the use of a higher index contrast or by adapting the grating depth[3

]. If not all light is extracted before the other end of the lightguide is reached, it is desirable to use a diffusing mirror at that side to recycle the light. To avoid mixing of different colors, we proposed recently[14

R. Caputo, L. De Sio, M.J.J. Jak, E.J. Hornix, D.K.G. de Boer, H.J. Cornelissen, and M.P.C. Krijn, “New system concept for colour separating backlights,” Asia Display 2007, in press.

] a new concept using a grating with a smaller period resulting in the blue light diffracted along the normal and green and red around it. If combined with a suitable pixel lay-out and a proper lens array that directs the light to the desired pixels, an efficient display system can be achieved.

The rigorous model can be used as a tool to optimize gratings for certain applications. In general it is suitable for simulations of subwavelength structure consisting of anisotropic materials for variant applications, such as, switchable holographic grating made of liquid-crystal films

Fig. 7. Contrast ratio between calculated intensities of the s- and p-polarized components of the -1^st transmitted order as function of diffraction angle θ^d for ϕⁱ = 67°, and for varying 0° ≤ ϕⁱ ≤ 90°, (a) for s-polarization incidence, and (b) for unpolarized incident field.

Fig. 8. Measured angular distribution of color-separated luminance for s-polarized (left) and p-polarized (center) light and angular distribution of s/p contrast ratio (right) for color-separating polarized backlight structure with TL 213 (refractive indices of Table 1) as birefringent material. (Note that the azimuthal angles are shifted by 90° with respect to those used before.)

Acknowledgments

The research is part of STW - Dutch Technology Foundation project DOE.6823.

References and links

1.	Y. Taira, D. Nakano, H. Numata, A. Nishikai, S. Ono, F. Yamada, M. Suzuki, M. Noguchi, R. Singh, and E.G. Colgan, “Low-power LCD using a novel optical system,” SID 02 Digest , 1313–1315 (2002). [CrossRef]
2.	F. Yamada, S. Ono, and Y. Taira, “Dual layered very thin flat surface micro prism array directly molded in an LCD cell,” Eurodisplay 2002 , 339–342 (2002).
3.	Dick K. G. de Boer, Roberto Caputo, Hugo J. Cornelissen, Chris M. van Heesch, Eefje J. Hornix, and Martin J.J. Jak, “Diffractive grating structures for colour-separating backlights,” Photonics in Multimedia, Proc. SPIE 6196, (2006).
4.	www.gsolver.com
5.	X. Wei, H.P. Urbach, and A.J.H. Wachters, “Finite Element Model for Three-Dimensional Optical Scattering Problems,” J. Opt. Soc. Am. A , 24, 866 (2007). [CrossRef]
6.	Max Born and Emil Wolf, “Rigorous diffraction theory,” in Principles of Optics, (The University Press, Cambridge, 2005), pp.633–673.
7.	K. Rokushima and J. Yamakita, “Analysis of anisotropic dielectric gratings,” J. Opt. Soc. Am. A , 73, 901 (1983). [CrossRef]
8.	E.N. Glytsis and T.K. Gaylord, “Rigorous three-dimensional coupled-wave diffraction analysis of single and cascaded anisotropic gratings,” J. Opt. Soc. Am. A , 4, 2061 (1987). [CrossRef]
9.	S. Mori, K. Mukai, J. Yamakita, and K. Rokushima, “Analysis of dielectric lamellar gratings coated with anisotropic layers,” J. Opt. Soc. Am. A , 7, 1661 (1990). [CrossRef]
10.	J.B. Harris, T.W. Preist, E.L. Wood, and J.R. Sambles, “Conical diffraction from multicoated gratings containing uniaxial materials,” J. Opt. Soc. Am. A , 13, 803 (1996). [CrossRef]
11.	L. Li, “Reformulation of the Fourier modal method for surface-relief gratings made with anisotropic materials,” J. Mod. Opt. , 45, 1313 (1998). [CrossRef]
12.	Xiuhong Wei , Three Dimensional Rigorous Model for Optical Scattering Problems , PhD thesis, Optics Research Group, Delft University of Technology, August 2006.
13.	J.P. Berenger, “Perfectly matched layer for the absorption of electromagnetic waves,” Journal of Computational Physics , 114(2), 185–200 (1994). [CrossRef]
14.	R. Caputo, L. De Sio, M.J.J. Jak, E.J. Hornix, D.K.G. de Boer, H.J. Cornelissen, and M.P.C. Krijn, “New system concept for colour separating backlights,” Asia Display 2007, in press.

OCIS Codes
(050.1950) Diffraction and gratings : Diffraction gratings
(160.1190) Materials : Anisotropic optical materials
(160.3710) Materials : Liquid crystals

ToC Category:
Diffraction and Gratings

History
Original Manuscript: January 24, 2007
Revised Manuscript: February 28, 2007
Manuscript Accepted: April 14, 2007
Published: April 26, 2007

Citation
Man Xu, H. Paul Urbach, and Dick K. G. de Boer, "Simulations of birefringent gratings as polarizing color separator in backlight for flat-panel displays," Opt. Express 15, 5789-5800 (2007)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-15-9-5789

Sort: Year | Journal | Reset

References

Y. Taira, D. Nakano, H. Numata, A. Nishikai, S. Ono, F. Yamada, M. Suzuki, M. Noguchi, R. Singh, and E. G. Colgan, "Low-power LCD using a novel optical system," SID 02 Digest, 1313-1315 (2002). [CrossRef]
F. Yamada, S. Ono, and Y. Taira, "Dual layered very thin flat surface micro prism array directly molded in an LCD cell," Euro display 2002, 339-342 (2002).
D. K. G. de Boer, R. Caputo, H. J. Cornelissen, C. M. van Heesch, E. J. Hornix, and M. J. J. Jak, "Diffractive grating structures for colour-separating backlights," Proc. SPIE 6196, 61960R (2006).
www.gsolver.com
X. Wei, H. P. Urbach, and A. J. H. Wachters, "Finite element model for three-dimensional optical scattering problems," J. Opt. Soc. Am. A 24, 866 (2007). [CrossRef]
M. Born and E. Wolf, "Rigorous diffraction theory," in Principles of Optics, (The University Press, Cambridge, 2005), pp. 633-673.
K. Rokushima and J. Yamakita, "Analysis of anisotropic dielectric gratings," J. Opt. Soc. Am. A 73, 901 (1983). [CrossRef]
E. N. Glytsis and T. K. Gaylord, "Rigorous three-dimensional coupled-wave diffraction analysis of single and cascaded anisotropic gratings," J. Opt. Soc. Am. A 4, 2061 (1987). [CrossRef]
S. Mori, K. Mukai, J. Yamakita, and K. Rokushima, "Analysis of dielectric lamellar gratings coated with anisotropic layers," J. Opt. Soc. Am. A 7, 1661 (1990). [CrossRef]
J. B. Harris, T. W. Preist, E. L. Wood, and J. R. Sambles, "Conical diffraction from multicoated gratings containing uniaxial materials," J. Opt. Soc. Am. A 13, 803 (1996). [CrossRef]
L. Li, "Reformulation of the Fourier modal method for surface-relief gratings made with anisotropic materials," J. Mod. Opt. 45, 1313 (1998). [CrossRef]
X. Wei, Three dimensional rigorous model for optical scattering problems, PhD thesis, Optics Research Group, Delft University of Technology, August 2006.
J. P. Berenger, "Perfectly matched layer for the absorption of electromagnetic waves," Journal of Computational Physics 114, 185-200 (1994). [CrossRef]
R. Caputo, L. De Sio, M. J. J. Jak, E. J. Hornix, D. K. G. de Boer, H. J. Cornelissen, and M. P. C. Krijn, "New system concept for colour separating backlights," Asia Display 2007, in press.

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper