## 1. Introduction

Magnetic materials have attracted considerable attention for their applications in high density magneto-optical storage devices, which is not only appealing scientifically, but also makes the magnetic materials promising for a wide range of applications [1-11]. Due to their unique optical properties, magnetic materials have been introduced into the photonic crystals, forming magnetophotonic crystals (MPCs), i.e., the photonic crystals with at least one magnetic material component [2,3]. Now the optical properties can be mediated and controlled by the external electric or magnetic field due to the existence of magnetic materials. The special importance of MPCs can be ascribed to the existence of magneto-optical effects, for example, Kerr effect and Faraday rotation [12-28]. They were discovered by Kerr and Faraday, respectively, and are now widely used in integrated optics and magneto-optical devices for magnetic domain imaging, mapping of hysteresis loops and high density recording [6].

Magnetic materials with large magneto-optical responses are always the attractive ones used in magnetophotonic crystals. In contrast to the corresponding three-dimensional magnetic structures, magnetic materials exhibit a larger magneto-optical effects due to the light's confinement in the MPCs, offering a genuine chance to put their optical responses into applications. The optimized MPCs have been shown to behave like mixed systems, with a coexistence of high transmittance and large magneto-optical effects [12-28]. All possible configurations are proposed to achieve strong magneto-opitcal effects, including the ordinary cavity-based, multilayered periodic and aperiodic structures. Furthermore, the diffracted magneto-optical enhancement is also demonstrated in the grating structures theoretically and experimentally, which greatly reduces the thickness of the device contrary to multilayered structures and miniaturizes magneto-optical devices in integrated optics [29,30].

The ability to tune the otpical properties by an external stimulus is a key issue of modern optoelectronics. Although most attention has so far been focused on the magneto-optical properties of the given structures, tunable magneto-optical devices have important applications and will be respected in optical switches and displays. There are few reports concerning with the tunable magnetophotonic crystals [31-33]. For example, it is possible to manipulate the magnetic order of magnetic conducting spheres using the magnetic field, thus forming the tunable magnetophotonic crystals [31]. Semiconductor quantum well has been discussed as possible candidates for achieving artificial tunable MPCs [32]. From the application-oriented perspective it would be very desirable if tunable magneto-optical effects could be achieved by controlling applied electric fields, thus making such effect a potential proposition. Therefore, to search the alternate scheme is of great importance in tunable magneto-optical electronics. It was known that nematics liquid crystals (NLCs) are good choices for tunable photonic crystals due to their unique sensitivity to temperature, the electric, magnetic field, or lights itself [34]. However, seldom reports concerning tunable magneto-optical effects based on the NLCs have been presented in the literatures. Therefore, it would be interesting to investigate a possible way of creating electrically controlled magneto-optical effects in the MPCs with the NLCs. The application of liquid crystals in the MPCs offers new opportunities for the tunable optoelectronic devices.

Liquid crystals are materials that display a phase of matter whose properties lie between those of a conventional liquid and a solid crystal. They are a class of materials particularly attractive for liquid crystal displays and optical electronic applications due to their high sensitivity to the external stimulus and have been studied experimentally extensively [35-37]. There are three basis kinds of liquid crystals: NLCs, cholesteric liquid crystals and smectic liquid crystals [38-47]. Here we only focus attention on the NLCs, which is characterized by molecules that have no positional order but tend to align along the same direction. Due to thermal random motion, friction and collision between molecules, not all molecules align along a certain direction and their directions vary around the average direction randomly. This average direction is referred to as the orientation of the liquid crystal, which stands for the average direction of most molecules. This parameter, i.e., the director, is an important factor to denote the liquid crystal's properties. Normally, the light waves with electric fields perpendicular or parallel to the director of a NLC have ordinary _{
o
} or extraordinary _{
e
} refractive indices, respectively. Therefore, the refractive index of a NLC may be changed between _{
o
} and _{
e
} by controlling the orientation of the directors using the applied electric field or adjusting the temperature [34]. Specifically, the refractive index of NLC is especially easy to control using the external electric field in one-dimensional cases [38-40]. Previous studies have demonstrated, experimentally as well as theoretically, how tunability is brought to electro-optical systems by employing NLCs, such as the realization of tunable band gaps in one-, two-, and three-dimensional photonic crystals based on NLCs [34, 38-40]. Tunable negative refraction is also achieved at the interface between air and the NLC, and that the refraction can be adjusted by an applied electric field or by the temperature [48-50]. Multistable all-optical switching have been demonstrated by utilizing an unique nonlinear coupling between light and the NLC in a periodic dielectric structure [51]. The absorption peak of nanoparticles doped in the NLC can be controlled by the external electric field [52]. The light-induced reorientational effects have been observed in a one-dimensional photonic crystal with the NLCs [53]. Especially, the orientation of the magntic anisotropy from ferromagnetic nanorods have been manipulated by means of electric fields in the composite of liquid crystals and ferromagnetic nanorods [54]. While the investigations of the magnetic properties in a electrically controlled way is still in infancy, the current experimental work suggests a possible chance of realizing the tunable MPCs with liquid crystals for fundamental research and practical applications.

Moreover, the particular nonreciprocity of magnetic materials, i.e., the polarized state (right/left handed polarized light) will switch into the opposite one (left/right handed polarized light) when the direction of the incident light is reversed, makes the MPCs remarkable workbenches for the unidirectional electromagnetic transmission [55-61]. Such a character shown in the cavity-based, periodic and aperiodic multilayered structures has made them good candidates to be isolators, unidirectional transmitted MPCs and one way waveguides with MPCs [57-59]. It has been unraveled that both ternary and microcavity configurations combined with waveguides can generate the isolators and circulators [55]. The asymmetric isolators are shown to be much efficient compared with symmetric ones [62]. The performance of small antennas embedded within MPCs constructed from periodic arrangements of homogenous and anisotropic material layers was demonstrated that the extraordinary high gain and enhanced power reception can be achieved [63]. Furthermore, it has been demonstrated that the optical Tamm states (OTS) can be observed in one dimensional MPCs in recent theoretical and experimental works [64-70]. There will be considerable interest to exploit and identify the OTS in the MPCs. The ability of creating and manipulating the surface states is central to the development of subwavelength microscopy. However, few reports concerning the tunable OTS in MPCs have been presented [66]. Therefore, it is interesting to find whether the NLCs can be used to realize the OTS in MPCs, or more importantly, whether such a configuration will support the controllable OTS. Controlling the OTS by electrical means is particularly interesting, as that would allow the integration of magnetic materials and the NLC with conventional optical electronics. Due to these unique optical properties, the MPCs with the NLC will hold substantial promise for possible optoelectronic devices.

In this chapter, we present a theoretical investigation on the NLC-based MPCs by performing a 4×4 transfer matrix method [71-73]. The optical properties of the MPCs with the NLCs will be reminiscent of both components. The coupling between the incident waves and the magnetic materials creates the magneto-optical effects, while the component of the NLCs is now responsible for its tunability occurring in the present configurations. The chapter is organized as follows: In Sec. II, we give a general description of the 4×4 transfer matrix method, which is a consequence of the combination of Maxwell equations and constitutive relations. The relations between the reflected, transmitted field vectors and the incident field vectors are listed to elucidate the magnitude of magneto-optical effects. In Section III, the case of periodic structure is considered, i.e., one dimensional MPC whose unit cell is composed of alternating NLCs and magnetic materials. In this part, the NLC is treated as a simple isotropic dielectric material approximately [74]. In Sec. IV, one-dimensional MPC infiltrated with the NLC is investigated, where the intrinsic anisotropic properties of the NLC are taken into account. Combined with transfer matrix method and a piecewise homogeneity approximation method for the NLC [75], the magneto-optical effects of MPCs with the NLC is definitely achieved. In Sec. V, the tunable OTS is reported theoretically in the MPCs with the NLC. Finally, we give a summary about this work.

Our investigations reveal a variety of interesting results which are the main characteristics of one-dimensional MPCs with the NLC. The magneto-optical effects appear the peaks due to the abnormal dispersion relations at the edge of the band gaps. It is shown the shift of the peaks and the enhanced values in the magneto-optical spectrum as the permittivity of NLCs increases, which demonstrates that the magneto-optical effects can be altered by changing the permittivity in the NLCs. When the NLC-based cavity structure is subject to the applied electric field, the magneto-optical effects are observed to exhibit the similar trends as those of the periodic one treating the NLC as an approximate isotropic dielectric material. A significant shift of the peaks in the magneto-optical spectrum is observed with the external voltage due to the directors' reorientation. In addition, the existence of intra-Brillouin-zone band gaps is demonstrated and they are predicted to be dependent on the applied voltages in the MPC with the NLC. If the present MPC is attached to the other photonic crystals with two dielectric materials, the OTS is also observed in the transmission spectrum and can be controlled by the applied voltages. Our results provide an approach to modify a interface state in a controllable manner, which is significant to its potential applications. Moreover, it is possible to design a larger tunable magneto-optical effects by simply selecting the appropriate liquid crystals in the MPCs with the NLC. The ability of creating the electrically-controlled magneto-optical effects can provide an extra space to be explored in the future. Such tunability of the magneto-optical effect may be useful for future application in electro-optical devices.

## 2. Theoretical treatment

A 4×4 transfer matrix method is an effective way to explore the optical properties of one dimensional multilayered configurations, which has been proved to be accurate enough to tackle the band structures of the uniaxial, bi-axially anisotropic dielectric materials and liquid crystals [71-72]. For the sake of convenience, we describe the theoretical method and detailed treatment in this section.

Fig. 1 shows a typical multilayered structure, which is sandwiched to the incident medium with the refractive index _{
i
} and transmitted medium with the refractive index _{
t
}, respectively. The incident wave vector is expressed as
* ω* is the angular frequency of the incident wave, and

*stands for the velocity of light in the vacuum. Here an electromagnetic wave impinges on the multilayered structure with an angle*c

θ

_{ i }, which is confined in the

*–*x

*plane indicating*z

k

_{ y }= 0,

v

_{ x }=

n

_{ i }sin

*and*θ

v

_{ z }=

n

_{ i }cos

*. The dielectric permittivity and magnetic permeability of any component are characterized by the tensors:*θ

For a monochrome light propagating along the * z* - axis with a given frequency

*, based on Maxwell's equations, we have [71]*ω

where Ψ(* z*) is the field vector provided as

*) is the Berrmeman matrix which is the function of the permittivity, permeability of the material and the incident wave vector. The subscript*z

*denotes the transpose operator. The components of the Berreman matrix are provided with [71]*t

Eq. (2) has a solution of the form:

Within the frame of transfer matrix method, the field vector can be expressed as,

where * T*(

*,*z

z

_{0}) is the transfer matrix provided with the Cauchy and Hamilton theorem [71]:

These coefficients _{
i
} are:

where

_{
i
} are the eigenvalues of the Berreman matrix. In contrast with the isotropic materials, it is more complicate due to the existence of the non-diagonal elements in the dielectric permittivity or magnetic permeability tensors. Even the incident plane wave is a linear polarization, it will form the non-linear polarized waves. Therefore, it typically yields the four solutions for the anisotropic materials, two of which represent left- and right- handed polarized forward propagating waves and another two are the corresponding backward propagating waves. If the individual layer occupies the spatial region of _{0}< * z*<

z

_{0}+

*, the field vectors Ψ(*d

z

_{0}+

*) and Ψ(*d

z

_{0}) can be connected with

where _{
i
} is the single-layer transfer matrix and * d* is the thickness of this layer. For the finite photonic crystals, the total transfer matrix

T

_{ total }of the system can be achieved by multiplying the individual transfer matrix together. The NLC is of spatially inhomogeneity, whose treatment needs to adapt the piecewise homogeneity approximation method [75].

*(*p

*) denotes electromagnetic wave parallel (perpendicular) to the incident plane*s

*–*x

*, respectively. The incident, reflected and transmitted field vectors are expressed as*z

Due to the continuous boundary conditions of electronic and magnetic fields, the equation Ψ_{
t
} = * T* (Ψ

_{ i }+ Ψ

_{ r }) should be satisfied. Accordingly, the relations between the reflected, transmitted and incident waves are connected with [79]

where * R* (

*) is the Jones reflectance (transmitted) matrices and*T

r

_{ ij }(

t

_{ ij }) is the ratio of the incident

*polarized electric field and the reflected (transmitted)*j

*polarized electric field. Then, the magneto-optical effects, i.e., Kerr effect and Faraday effect can be determined by [79]*i

where
* p*(

*) polarized wave, respectively. Based on these equations, the magneto-optical effects and the corresponding transmission/reflectance coefficients can be obtained directly.*s

## 3. Electrically-controlled magneto-optical effects in magnetophotonic crystals consisting of magnetic materials and NLC

In this part, we investigate the magneto-optical effects of a one dimensional MPC whose unit cell is composed of a magnetic material and NLC and give a detailed explanation to the obtained results.

Generally, the permittivity of NLC is provided by
_{
o
} and _{
e
} are the respective dielectric permittivities for light polarized parallel and perpendicular to the director axis
* ϴ* with respect to the optical wave vector. Normally, losses are negligible for typical micron thick NLC in the optical frequency region. In the following numerical simulation, we treat the NLC as a homogeneous isotropic dielectric layer with the permittivity of

*of the aligned NLC with respect to the wave vector, which results in the variation of the extraordinary refractive index from*ϴ

Fig. 2 shows a schematic illustration of a one-dimensional finite MPC, i.e., (NLC/magnetic materials). In the calculations, the total periodic number * N* is taken as 8 and the thicknesses of NLC and magnetic layers are

d

_{1}= 500

*and*nm

d

_{2}= 20

*, respectively. The input and output media are supposed to be vacuum and glass, i.e.,*nm

n

_{ i }= 1 and

n

_{ t }= 1

*5, respectively. The dielectric tensor of magnetic layers has an antisymmetric form and is denoted as [80]*.

where the magnetization vector makes an angle Θ with the * z* axis and its projection on the

*–*x

*plane makes the azimuth Φ with the*y

*axis. As an illustrative example, the permittivity parameters of magnetic material are taken as*x

ε

_{ xx }=

ε

_{ yy }=

ε

_{ zz }= –4

*8984 + 19*.

*415*.

*,*i

ε

_{ xy }= –

ε

_{ yx }= 0.4322 + 0.0058

*[81] while the permeability of magnetic material is nearly 1 at the optical frequency. Note that there are three cases according to the different values of Θ and Φ: the polar Kerr effect (Θ = Φ = 0), the longitudinal Kerr effect (Θ=*i

ε

_{ LC }are calculated when the light is incident normally to MPC. Since normal incident light leads to decoupling between

*- and*p

*-polarized waves and the two polarized lights are degenerate, only results for the*s

*-wave are presented in the following figures.*p

The solid, dash and dotted lines in Fig. 3 (a) describe the behaviors of _{
p
} for _{
LC
} = 2* ,* 2

*5*.

*3, respectively. Though the choice of the value of*,

ε

_{ LC }seems arbitrary, the general tendency can still be seen clearly in the figures. Compared with the wavelength dependence of

Ө

_{ p }of a single magnetic layer as shown in the inset of Fig. 3(a), the Kerr rotation angles exhibited in the NLC-based MPC structures are enhanced dramatically, especially at certain wavelengths. The maximum value of

Ө

_{ p }increases from

Ө

_{ pmax }= –0

*86° for*.

ε

_{ LC }= 2 up to –8° for

ε

_{ LC }= 3. The latter is nearly 40 times larger than that of the single magnetic layer(

Ө

_{ pmax }= –0

*19°). The enhancement of*.

Ө

_{ p }compared with the single magnetic layer originates from the weak localization of light caused by the multiple interference due to the nonreciprocal properties of the Kerr effect [2]. Furthermore, the emergence of the

Ө

_{ p }-peaks is also an important feature of the multilayered structure in our MPC model. It is well-known that the Kerr effect is strengthened dramatically due to the abnormal dispersion relation at the edges of the photonic band gaps [28]. To make it more clear, Fig. 3 (b) shows the dependence of the reflectance of MPC on the wavelength. We clearly observe that the positions of the remarkable enhancement of

Ө

_{ p }are coincident with the positions of photonic band gap exactly.

Note that all the peaks in Fig. 3(a) shift towards longer wavelengths when the NLC's permittivity increases, which is particularly interesting since it reflects the good tunability of the Kerr effect in our model. A similar shift is also observed for the photonic band gaps as shown in Fig. 3(b). This result is reasonable since the positions of photonic band gaps depend on the wave impedance ratio of two components [82, 83], while the latter changes with the NLC's permittivity in the present case. For a rough estimation, the central wavelength of the NLC-based MPC structure can be evaluated by _{
c
} = 2(_{
LC
}
_{1}+_{
m
}
_{2}) approximately [84], where _{
m
} is the complex refractive index of magnetic layers and invariable with temperature or applied field. So the increase of the refractive index of NLC will directly lead to the enhancement of the central wavelength of the MPC structure. Since the refractive index of liquid crystal may be tuned between the ordinary and extraordinary refractive indices by controlling the orientation of the molecules under the influence of the applied electric field, the upper (lower) limit of the red shift is determined by the extraordinary (ordinary) permittivity, i.e., _{
LC
} = 3 (_{
LC
} = 2). Therefore, the NLC with larger optical anisotropy may behave as a better candidate for realizing larger tunability.

Fig. 3(a) also shows the significant increase of peak height with the permittivity of NLC in the spectrums of the Kerr rotation angle, which can be explained according to the circular birefringence. Different refractive indices due to different localization conditions, i.e., the localization wavelengths for left- and right- circular polarized lights, leads to the rotation of the reflected light [2]. In photonic crystals, the effective parameters of the unit cell can be expressed by the effective homogeneous anisotropic medium method [85]. The effective permittivity of MPC consisting of magnetic materials and NLC is given by

where
_{
LC
}, the difference between the dielectric permittivity of the left- and right- circular polarized lights is enlarged, which leads to a relatively larger Kerr rotation angle at a longer wavelength.

Generally, refractive indices of liquid crystal are fundamentally interesting and practically useful parameters in the calculation of optical properties. By changing the permittivity of the liquid crystals, the magneto-optical effects varied in a controlled fashion, indicating that the engineering of magneto-optical devices relies on the control of the applied electric field. Here we focus attention on the change of the orientation of the molecules under the influence of the applied electric field at a fixed temperature. Actually, besides the external electric field, temperature is also an important factor affecting the liquid crystal refractive indices, which offers an alternative approach to control the magneto-optical effects, i.e., a temperature tunable MPC.

## 4. Voltage-controlled magneto-optical effects in cavity- based magnetophotonic crystal with the NLC

We have predicted that the opportunity of creating electrically controlled Kerr effect in magnetic multilayered structures [76]. It was suggested that a one-dimensional MPC composed of alternating NLC and magnetic materials can create tunable Kerr effect by considering the properties of liquid crystals and thus provide for tunable MPCs. However, we just employed an approximate isotropic treatment of NLC to analyze magneto-optical effects, which is a rough theoretical evaluation. It is generally known that the directors of NLCs exhibit inhomogeneous distribution under the influence of the applied electrical field [40]. Therefore, a rigorous anisotropic treatment of NLCs is employed to consider the NLC director’s spatially inhomogeneous property upon an applied external voltage, which is based on the Newton method and continuous elastic theories. Although we expect such the tunable magneto-optical effect to be seen with a variety of patterns containing NLCs, we investigated a multilayered structure infiltrated with NLC, as shown schematically in Fig. 4. The defect of the NLC ensures that the present structure is sensitive to the external electric field and therefore shows the tunability of magneto-optical effects for observations at normal incidence.

The gray and white regions are magnetic materials, yttrium-iron-garnet(YIG) and non-magnetic materials gadolinium-gallium garnet (GGG) with permittivity _{2}, respectively. The magnetization vector is oriented along the z axis in the YIG layers, which are characterized by antisymmetrical permittivity and have the following nonzero components, _{
xx
} = _{
yy
} = _{
zz
} = _{1}, _{
xy
} = –_{
yx
} = * ig*, where

*is the gyrotropy of the magnetic layer. The thicknesses of magnetic and nonmagnetic materials are*g

d

_{1}and

d

_{2}, respectively. The NLC is characterized by a long-range uniaxial orientational order of rod-like anisotropic molecules along a common director, whose dielectric tensor is given by [86]

where * α* is the tilt angle between the principal axis of NLC molecule and the x axis. Δ

*is the optical dielectric anisotropy given as*ε

n

_{0}and

n

_{ e }are ordinary and extraordinary refractive indices, respectively. Upon an applied external voltage, the distribution of NLC director exhibits spatially inhomogeneous property. The distribution of

*is given in the NLC slab by the Ossen Frank elastic theory when an electric field is applied along the z axis. The relations between*α

*and the voltage*α

*applied to NLC slab can be achieved in the light of the following equations [87]*V

with * η* = sin

α

_{ m }and

*= (*k

k

_{3}–

k

_{1})/

k

_{1}. The parameters

k

_{1}and

k

_{3}are the Frank elastic constants for splay and bend modes of the nematic director distortion, respectively.

α

_{ m }is the maximum tilt angle in the NLC slab, Δ

ε

_{0}is the dielectric anisotropy at zero frequency. Note that the direction of NLC molecules is affected by the external voltage only when the applied voltage

*exceeds a threshold value*V

For the one-dimensional multilayered structure (GGG|YIG)^{n}(YIG|GGG)^{m}NLC(YIG|GGG)^{n}, at the near infrared wavelength * λ* = 1.55

*m, YIG and GGG are characterized by*μ

ε

_{1}= 5.5 +

*0.0025,*i

*= (1 –*g

*0*i

*15) ×10*.

^{–2}, and

ε

_{2}= 3.709, respectively [80]. The thicknesses are

d

_{1}= 0.466

*and*D

d

_{2}= 0.534

*, respectively.*D

k

_{1}= 7.4pN,

k

_{3}= 10.2pN, Δε

_{0}= 1.43 ×10

^{–10}F

*m,*/

n

_{ e }= 1.75,

n

_{0}= 1

*54 for typical NLC (5CB) [72],*.

d

_{ LC }=

*. The repetition numbers are*D

*= 16 and*n

*= 29, respectively.*m

Fig. 5(a) corresponds to the distribution of * α* with different voltage

*= 0*V

*2*,

V

_{ c }, 5

V

_{ c }for

*-polarized wave. The anisotropic character of NLC can be considered exactly by utilizing the obtained distribution of*p

*. The calculated results for*α

ϴ

_{ p }as functions of normalized frequency

*/(2*ωD

*) with different voltages are plotted in Fig. 5(b). The solid, dashed, and dot-dashed lines correspond to*πc

ϴ

_{ p }with

*= 0, 2*V

V

_{ c }, 5

V

_{ c }, respectively. It can be seen clearly that the spectrum of

ϴ

_{ p }moves at different voltages, which is a direct manifestation of voltage-controlled Kerr effect. The frequency dependence of

ϴ

_{ p }in the present case exhibits the following features. First, the positions of

ϴ

_{ p }-peaks are dependent on the applied voltage due to the reorientation of the director. It is seen that the positions of peaks move toward high frequencies with the applied voltages. This is because of NLC's high sensitivity to the applied voltage. When the applied voltage is very small, the distribution of the director aligns

*-direction due to the surface anchoring. With the increase of applied voltage, the energy of electric field adds gradually. When the applied voltage is bigger than the threshold value, the function of electric field plays a leading role. In such a case, the alignment of the director along the electric field can be observed. The reorientation of the NLC's director will lead to the change of its permittivity, which directly affects the Kerr rotation angles. Second, the amplitudes of*x

ϴ

_{ p }depend strongly on the applied voltage, i.e., the heights of

ϴ

_{ p }-peaks increase considerably with the increase of

*. For example, the maximum of*V

ϴ

_{ p }at

*= 5*V

V

_{ c }(

ϴ

_{ p max}= –30°) is approximately 3 times larger than that of

ϴ

_{ p }(

ϴ

_{ p max}= –10°) at

*= 0, which describes the large difference in the indices of refraction between the right and left circularly polarized light. This can be understood according to the circular birefringence, in which different refractive indices due to different localization wavelengths for left- and right- circularly polarized lights, leads to the large rotation of the reflected light [2]. In addition, the positions of*V

ϴ

_{ p }-peaks agree with those of the dips in the reflection spectrum exactly. Take

*= 5*V

V

_{ c }as an example, we calculated

ϴ

_{ p }and the corresponding reflection spectrum, which are exhibited in Fig. 5(c) and (d). Such a phenomenon is directly related to the edges of photonic band gaps, where the abnormal dispersion relation leads to strong photonic localization. As a result,

ϴ

_{ p }can be enhanced at these specific frequencies in the present structure [28].

Thus, the controllable magneto-optical effects have been predicted in the MPCs with the NLC, which will be a critical problem in the optoelectronic applications.

## 5. Voltage-controlled tamm state in periodic magnetophotonic crystal with the NLC

The Tamm state proposed by Tamm is one of the most fundamental physical properties of the interface, which means electron states can occur in the energy band gap at a crystal surface [68]. In analogy with the electronic case, the OTS can be realized at the interface between two different photonic crystals with all isotropic dielectric materials. Conventional surface waves with a wave vector exceeding that of light in an incident medium decay exponentially away from the surface. In comparison with the conventional surface waves, the OTS can be formed for both the * s*– and

*–polarized waves and occurs even at normal incidence [70]. In this part, we consider a one-dimensional MPC composed of the NLC and magnetic material, as illustrated in Fig. 6 (a). A linearly*p

*-polarized monochromatic light impinges normally onto the structure. The thicknesses of the NLC and magnetic materials are*p

d

_{1}and

d

_{2}, respectively. Due to the periodicity of the structure the eigenwaves satisfy Bloch's theorem, which leads to the equation det[

*–*P

e

^{ ikd }] = 0. Here

*is the total propagation matrix of unit cell. Based on these equations, we may illustrate the characteristics of the normalized angular frequency and*P

*in the present configuration.*kd

The parameters used in our calculations are taken as _{1} = 3 and * g* = 0.5, where a larger value of

*is used to make the band gaps easily observable.*g

k

_{1}= 7.4pN,

k

_{3}= 10.2pN, Δ

ε

_{0}= 1.43 × 10

^{-10}F

*m,*/

n

_{ e }= 1.75,

n

_{0}= 1.54 for the NLC (5CB)[72], with the thicknesses

d

_{1}= 0.5

*,*D

d

_{2}= 0.5

*, respectively. The spatially inhomogeneous property of the directors is also taken into account in this calculation. In Fig. 7 (a), we present the*D

*-*ωd / c

*characteristics under the applied voltages*kd

*= 0*V

*2*,

V

_{ c }

*5*,

V

_{ c }, respectively. From the figure it is clear that the intra-Brillouin-zone band gaps do occur and obviously depend on the applied voltages. This is a remarkable result in itself, as it indicates the existence and tunability of the OTS with the NLC. To understand qualitatively the origin of the appearance of intra-Brillouin-zone band gaps, we start with the simple case, i.e., we set the applied voltage as zero, where the distribution of the directors aligns along

*-direction due to the surface anchoring. Now the NLC is a general anisotropic dielectric layer, where the eigenmodes are the ordinary and extraordinary waves when the incident wave enters into the NLC [88], then they become left- and right-circularly polarized waves in the magnetic materials. It appears an approximate standing wave due to the strong interference between these forward and backward waves, which leads to the formation of the intra-Brillouin-zone band gaps [65]. More interesting, the*x

*dependence on*kd

*in the present case exhibits the following feature, i.e., the intra-Brillouin-zone band gap moves with the applied voltages. The driving mechanism for manipulating the OTS arise from the intrinsic properties of the NLC. It is mainly caused by the sensitivity of the NLC on the applied voltage, which has significantly influence on the distribution of the NLC's directors. As widely known, the directors of the NLC molecules are affected by the external electric fleld only when the voltage*ωd/c

*exceeds a threshold value*V

V

_{ c }, where

V

_{ c }is given by

*/*V

V

_{ c }> 1. The reorientation of the NLC's directors leads to the variation of the effective refractive index, which change the optical path directly. Subsequently, the positions of standing waves will be finally changed with the applied voltages.

Furthermore, two stacked photonic crystals are generally used to illuminate the existence of the OTS in experiments [64], which can be observed by studying the transmission/reflection spectra with a sharp narrow peak/dip in the finite photonic crystal. Thus, it is necessary for us to provide the transmission spectrum. We design the two stacked photonic crystals, shown in Fig. 6(b), one of which consists of two dielectric materials and the other is composed of the NLC and magnetic material. Now we use the realistic material parameters at the wavelength * λ* = 1.55

*,*μm

*07,*.

ε

_{1}= 5.58 and

*= –0.00198 for Bi:DyIG. The thicknesses are*g

*696*.

*,*μm

*37*.

*,*μm

d

_{ NLC }= 1

*and*μm

*, respectively [67]. Fig. 7 (b) depicts the transmission of the two photonic crystals as a function of wavelength at normal incidence. It is observed that a sharp narrow peak indeed appears in the transmission spectrum. In particular, we notice two photonic crystals generate photonic gaps nearly from 1.12 – 1.17*μm

*and 1.13 – 1.16*μm

*, respectively. However, the peak appears in the band gaps of two photonic crystals, which significantly demonstrates the sharp peak comes from the resonant tunneling of the electromagnetic wave through the OTS at the interface. In addition, the peak undergoes slightly blueshift with the applied voltages, which manifest itself a good candidate to be a controlled OTS structure. Finally, it is noted that the existence of OTS survives for*μm

*-polarized wave but the tunability disappears because*s

ε

_{ yy }is not dependent on the external voltage.

## 6. Conclusion

The tunable magneto-optical effects in magnetophotonic crystals with the NLC have been reviewed in this chapter. Two types of MPCs are studied, i.e., MPCs consisting of alternate magnetic materials and the NLC; the cavity-based MPCs with the NLC. Both of them exhibit the tunable magneto-optical properties qualitatively different from those of the unvariable MPCs without the NLC. It is predicted theoretically that the magneto-optical effects can be manipulated by the applied voltages in the theoretical treatments, i.e., the 4×4 transfer matrix method in combination with a piecewise homogeneity approximation for liquid crystals. A significant shift of the peaks' positions in magneto-optical spectrum is observed toward high frequencies with the applied voltages due to the high sensitivity of the directors on the external electric fields. By applying a tunable electric field, the voltage-induced reorientation of the directors of the NLCs alters the values of the dielectric permittivity in the NLC, thus leading to alternating magneto-optical effects. In addition, the tunable intra-Brillouin-zone band gaps can be realized in the MPCs with the NLC, which implies the existence of the OTS at the interface between two photonic crystals. The methods described here may possible to be implemented on many other periodic multilayered structures, cavity-based periodic structure, quasi-periodic and aperiodic structures, thus, opening the possibility to systematically investigate the magneto-optical effects in a controlled way. The structures we propose here is not only interesting in itself but also allows easy access to fabricate it in the experiments due to the mature techniques for one-dimensional systems. In fact, apart from the external applied electric field, the variation of temperature can also change the permittivity in the NLC component, and then alter the magneto-optical effects. It is obvious that there are a great space for the potential application of the MPCs with NLC. These results highlight an intriguing avenue for future investigations in the development of tunable liquid crystal-based magnetophotonic crystals optoelectronic devices.

## Acknowledgments

This project was supported by the National Natural Science Foundation of China under Grant No.10774107 and Jiangsu provincial innovation project under Grant No. ZY320607.

## References

- 1.
Mansuripur M. 1995 - 2.
Inoue M. Fujikawa R. Baryshev A. Khanikaev A. Lim P. B. Uchida H. Aktsipetrov O. Fedyanin A. Murzina T. Granovsky A. 2006 "Magnetophotonic crystals", J. Phys. D: Appl. Phys. 39, R151. - 3.
Mitsuteru Inoue. Alexander V. Baryshev Alexander B. Khanikaev Maxim E. Dokukin Kwanghyun Chung. Jinheo Hiroyuki. Takagi Hironaga Uchida. Pang Boey. Lim Jooyoung Kim. 2008 "Magnetophotonic Materials and Their Applications", IEICE TRANS. ELECTRON, 91, 1630. - 4.
Fujikawa R. Tanizaki K. Baryshev A. V. Lim P. B. Shin K. H. Uchida H. Inoue M. 2006 "Magnetic field sensors using magnetophotonic crystals", Proc. SPIE, 6369, 63690,. - 5.
Park J. H. Takagi H. Nishimura K. et al. 2003 "Magneto-optic spatial light modulators driven by an electric field", J. Appl. Phys., 93, 8525,. - 6.
Park H. J. Cho J. K. Nishimura K. Inoue M. 2002 "Magneto-optic spatial light modulator for volumetric digital recording system", Jpn. J. Appl. Phys., 41, 1813,. - 7.
Wang Z. Fan S. 2005 "Optical circulators in two-dimensional magneto-optical photonic crystals", Opt. Lett., 30, 1989,. - 8.
Merzlikin A. M. Vinogradov A. P. Inoue M. et al. 2006 "The Faraday effect in two-dimensional magneto-photonic crystals", Journal of Magnetism and Magnetic Materials 300, 108,. - 9.
Merzlikin A. M. Vinogradov A. P. 2006 "Superprism effect in 1D photonic crystal", Opt. Commun. 259, 700,. - 10.
Belotelov VI Kotov VA Zvezdin AK 2006 "New magnetooptical materials on a nanoscale", Phase Transitions, 79, 1135. - 11.
Belotelov VI Zvezdin AK 2005 "Magneto-optical properties of photonic crystals", J. Opt. Soc. Am. B 22, 286 - 12.
Inoue M. Arai K. Fujii T. Abe M. 1998 "Magneto-optical properties of one-dimensional photonic crystals composed of magnetic and dielectric layers", J. Appl. Phys. 83, 6768. - 13.
Inoue M. Fujii T. 1997 "A theoretical analysis of magneto-optical Faraday effect of YIG films with random multilayer structures", J. Appl. Phys. 81, 5659, - 14.
Kato H. Inoue M. 2002 "Reflection-mode operation of oned imensional magnetophotonic crystals for use in film-based magneto-optical isolator devices", J. Appl. Phys. 91, 7017,. - 15.
Kato H. Matsushita T. Takayama A. et al. 2002 "Properties of one dimensional magnetophotonic crystals for use in optical isolator devices", IEEE Trans. Magn. 38, 3246,. - 16.
Kato H. Matsushita T. Takayama A. Egawa M. Nishimura K. Inoue M. 2003 "Theoretical analysis of optical and magneto-optical properties of one-dimensional magnetophotonic crystals", J. Appl. Phys. 93, 3906 - 17.
Miguel Levy 2006 "Normal modes and birefringent magnetophotonic crystals", J. Appl. Phys. 99, 073104. - 18.
Miguel Levy Rong Li 2006 "Polarization rotation enhancement and scattering mechanisms in waveguide magnetophotonic crystals", Appl. Phys. Lett. 89, 121113. - 19.
Levy M. Jalali A. A. 2007 "Band structure and Bloch states in birefringent one-dimensional magnetophotonic crystals: An analytical approach, J. Opt. Soc. Am. B 24, 1603,. - 20.
Khartsev SI; Grishin AM; 2005 "[Bi3Fe5O12/Gd3Ga5O12](m) magneto-optical photonic crystals", Appl. Phys. Lett. 87, 122504. - 21.
Khartsev S. I. Grishin A. M. 2007 "High performance [Bi3Fe5O12/ Sm3Ga5O12]m magneto-optical photonic crystals", J. Appl. Phys. 101, 053906,. - 22.
Khanikaev A. B. Baryshev A. V. Inoue M. et al. 2005 "Two-dimensional magnetophotonic crystal: Exactly solvable model", Phys. Rev. B 72, 035123, - 23.
Sakaguchi S. Sugimoto N. 1999 "Multilayer films composed of periodic magneto-optical and dielectric layers for use as Faraday rotators", Opt. Commun. 162, 64 - 24.
Merzlikin A. M. Vinogradov A. P. Dorofeenko A. V. et al. 2007 "Controllable Tamm states in magnetophotonic crystal", Physica B 394, 277,. - 25.
Vinogradov A. P. Dorofeenko A. V. Erokhin S. G. et al. 2006 "Surface state peculiarities in one-dimensional photonic crystal interfaces", Phys. Rev. B 74, 045128, - 26.
Lyubchanskii I. L. Dadoenkova N. N. Lyubchanskii M. I. et al. 2006 "Response of two-defect magnetic photonic crystals to oblique incidence of light: Effect of defect layer variation", J. Appl. Phys. 100, 096110,. - 27.
Kahl S. Grishin A. M. 2004 "Enhanced Faraday rotation in all-garnet magneto-optical photonic crystal", Appl. Phys. Lett. 84, 1438 - 28.
Zhdanov A. G. Fedyanin A. A. Aktsipetrov O. A. Kobayashi D. Uchida H. Inoue M. 2006 "Enhancement of Faraday rotation at photonic-band-gap edge in garnet-based magnetophotonic crystals", Journal of Magnetism and Magnetic Materials 300, 253 - 29.
Yuehui Lu Min Hyung Cho JinBae Kim YoungPak Lee Jooyull Rhee Jae-Hwang Lee 2008 "Control of Diffracted Magneto-Optical Enhancement in Ni Gratings", IEEE TRANSACTIONS ON MAGNETICS 44, 3300, - 30.
Lu Y. H. Cho M. H. Kim J. B. Lee G. J. Lee Y. P. Rhee J. Y. 2008 "Magneto-optical enhancement through gyrotropic gratings", Optics Express 16, 5378 - 31.
Golosovsky M. Neve-Oz Y. Davidov D. 2005 "Magnetic-field-tunable photonic stop band in a three-dimensional array of conducting spheres", Phys. Rev. B 71, 195105 - 32.
Jiang-Tao Liu Kai Chang 2007 "Tunable giant Faraday rotation of exciton in semiconductor quantum wells embedded in a microcavity", Appl. Phys. Lett. 90, 061114 - 33.
Chernovtsev S. V. Belozorov D. P. Tarapov S. I. 2007 "Magnetically controllable 1D magnetophotonic crystal in millimetre wavelength band", J. Phys. D: Appl. Phys. 40, 295C299 - 34.
Graugnard E. King J. S. Jain S. Summers C. J. Zhang-Williams Y. Khoo I. C. 2005 "Electricfield tuning of the Bragg peak in large-pore TiO2 inverse shell opals", Phys. Rev. B 72, 233105, - 35.
Mach P. Wiltzius P. Megens M. Weitz D. A. K. Lin H. Lubensky T. C. Yodh A. G. 2002 "Electro-optic response and switchable Bragg diffraction for liquid crystals in colloid-templated materials", Phys. Rev. E 65, 65, - 36.
Larsen T. T. Bjarklev A. Hermann D. S. Broeng J. 2003 "Optical devices based on liquid crystal photonic bandgap fibers", Optics Express 11, 2589, - 37.
Hu W. Dickie R. Cahill R. Gamble H. Ismail Y. Fusco V. Linton D. Grant N. Rea S. 2007 "Liquid crystal tunable mm wave frequency selective surface", IEEE Microw. Wireless Compon. Lett. 17, 667, - 38.
Cos J. Ferre-Borrull J. Pallares J. Marsal L. F. 2009 "Tunable FabryCProt filter based on one-dimensional photonic crystals with liquid crystal components", Optics Communications 282, 1220, - 39.
Arkhipkin V. G. Gunyakov V. A. Myslivets S. A. Ya V. Zyryanov Shabanov V. F. 2007 "Angular tuning of defect modes spectrum in the one-dimensional photonic crystal with liquid-crystal layer", Eur. Phys. J. E 24, 297 - 40.
Ryotaro Ozakia Hiroshi Moritake Katsumi Yoshino Masanori Ozaki 2007 "Analysis of defect mode switching response in one-dimensional photonic crystal with a nematic liquid crystal defect layer", J. Appl. Phys. 101, 033503 - 41.
Jiun-Yeu Chen Lien-Wen Chen 2005 "Polarization-dependent filters based on chiral photonic structures with defects", J. Opt. A: Pure Appl. Opt. 7 558, - 42.
Hiroyuki Yoshida Chee Heng Lee Akihiko Fujii Masanori Ozaki 2007 "Tunable Chiral Photonic Defect Modes in Locally Polymerized Cholesteric Liquid Crystals", Mol. Cryst. Liq. Cryst. 477, 255, - 43.
Yuhua Huang Ying Zhou Charlie Doyle Shin-Tson Wu 2006 "Tuning the photonic band gap in cholesteric liquid crystals by temperature-dependent dopant solubility", Optics Express, 14, 1236 - 44.
Scaramuzza N. Ferrero C. Carbone B. V. Versace C. 1995 "Dynamics of selective reflections of cholesteric liquid crystals subject to electric fields", J. Appl. Phys. 77, 572 - 45.
Kuniaki Konishi. Benfeng Bai. Xiangfeng Meng. Petri Karvinen. Jari Turunen. Yuri P. Svirko Makoto-Gonokami. Kuwata. 2008 "Observation of extraordinary optical activity in planar chiral photonic crystals", Optics Express, 16, 7189 - 46.
Gorkunov M V Osipov M A 2008 "Molecular theory of layer contraction in smectic liquid crystals, J. Phys. Condens. Mater 20, 465101 - 47.
Barbero G. Komitov L. 2009 "Temperature-induced tilt transition in the nematic phase of liquid crystal possessing smectic C-nematic phase sequence", J. Appl. Phys. 105, 064516 - 48.
Qian Zhao Lei Kang Bo Li Ji Zhou Hong Tang Baizhe Zhang 2006 "Tunable negative refraction in nematic liquid crystals", Appl. Phys. Lett. 89, 221918 - 49.
Lei Kang Qian Zhao Bo Li Ji Zhoua Hao Zhu 2007 "Experimental verification of a tunable optical negative refraction in nematic liquid crystals", Appl. Phys. Lett. 90, 181931 - 50.
Jeremy A. Bossard Xiaotao Liang. Ling Li. Seokho Yun. Douglas H. Werner Brian E. Weiner Theresa S. Mayer Paul F. Cristman Andres Diaz. Khoo I. C. 2008 "Tunable Frequency Selective Surfaces and Negative-Zero-Positive Index Metamaterials Based on Liquid Crystals", IEEE Transactions on Antennas and Propagation, 56, 1308, - 51.
Miroshnichenko A. E. Brasselet E. Kivshar Y. S. 2008 "All-optical switching and multistability in photonic structures with liquid crystal defects", Appl. Phys. Lett. 92, 253306 - 52.
Li-Hsuan Hsu Kuang-Yao Lo Shih-An Huang Chi-Yen Huang Chung-Sung Yang 2008 "Irreversible redshift of transmission spectrum of gold nanoparticles doped in liquid crystals", Appl. Phys. Lett. 92, 181112 - 53.
Laudyn U. A. Miroshnichenko A. E. Krolikowski W. et al. 2008 "Observation of light-induced reorientational effects in periodic structures with planar nematic-liquid-crystal defects", Appl. Phys. Lett. 92, 203304 - 54.
Lin T. J. Chen C. C. Lee W. et al. 2008 "Electrical manipulation of magnetic anisotropy in the composite of liquid crystals and ferromagnetic nanorods", Appl. Phys. Lett. 93, 013108 - 55.
Alexander B. Khanikaev Steel M. J. 2009 "Low-symmetry magnetic photonic crystals for non-reciprocal and unidirectional devices", Optics Express 17, 5265 - 56.
Wang Z. Chong Y. D. Joannopoulos J. D. Soljacic M. 2008 "Reflection-Free One-Way Edge Modes in a Gyromagnetic Photonic Crystal", Phys. Rev. Lett. 100, 01390501 - 57.
Yu Z. Veronis G. Wang Z. Fan S. 2008 "One-Way Electromagnetic Waveguide Formed at the Interface between a Plasmonic Metal under a Static Magnetic Field and a Photonic Crystal", Phys. Rev. Lett. 100, 02390201 - 58.
Yu Z. Wang Z. Fan S. 2007 "One-way total reflection with one-dimensional magneto-optical photonic crystals", Appl. Phys. Lett. 90, 121133 - 59.
Figotin A. Vitebskiy I. 2003 "Electromagnetic unidirectionality in magnetic photonic crystals", Phys. Rev. B 67, 165210 - 60.
Figotin A. Vitebskiy I. 2001 "Nonreciprocal magnetic photonic crystals", Phys. Rev. E 63, 066609 - 61.
Jung K. Y. Donderici B. Teixeira F. L. 2006 "Transient analysis of spectrally asymmetric magnetic photonic crystals with ferromagnetic losses", Phys. Rev. B 74, 165207 - 62.
Ruiyi Chen Dongjie Tao Haifeng Zhou Yinlei Hao Jianyi Yang Minghua Wang Xiaoqing Jiang 2009 "Asymmetric multimode interference isolator based on nonreciprocal phase shift", Optics Communications 282 862. - 63.
Mumcu G. Sertel K. Volakis J. L. 2006 "Miniature Antennas and Arrays Embedded Within Magnetic Photonic Crystals", IEEE 5, 168 - 64.
Goto T. Dorofeenko A. V. Merzlikin A. M. Baryshev A. V. Vinogradov A. P. Inoue M. Lisyansky A. A. Granovsky A. B. 2008 "Optical Tamm States in One-Dimensional Magnetophotonic Structures", Phys. Rev. Lett. 101, 113902 - 65.
Fei Wang Akhlesh Lakhtakiab 2008 "Intra-Brillouin-zone bandgaps due to periodic misalignment in one-dimensional magnetophotonic crystals", Appl. Phys. Lett. 92, 011115 - 66.
Merzlikina A. M. Vinogradova A. P. Dorofeenkoa A. V. Inoueb M. Levyc M. Granovskyd A. B. 2007 "Controllable Tamm states in magnetophotonic crystal", Phys. B 394, 277 - 67.
Vinogradov A. P. Dorofeenko A. V. Erokhin S. G. Inoue M. Lisyansky A. A. Merzlikin A. M. Granovsky A. B. 2006 "Surface state peculiarities in one-dimensional photonic crystal interfaces", Phys. Rev. B 74, 045128 - 68.
Tamm I. E. 1932 "On the possible bound states of electrons on a crystal surface " Phys. Z. Sowjetunion, 1, 733 - 69.
Kaliteevski M. Iorsh I. Brand S. Abram R. A. Chamberlain J. M. Kavokin A. V. Shelykh I. A. 2007 "Tamm plasmon-polaritons: Possible electromagnetic states at the interface of a metal and a dielectric Bragg mirror", Phys. Rev. B 76, 165415 - 70.
Kavokin A. V. Shelykh I. A. Malpuech G. 2005 "Lossless interface modes at the boundary between two periodic dielectric structures", Phys. Rev. B 72, 233102 - 71.
Wohler H. Haas G. Fritsch M. Mlynski D. A. 1988 "Faster 4 x 4 matrix method for uniaxial inhomogeneous media", J. Opt. Soc. Am. A 5, 1554 - 72.
Berreman D. W. 1973 "Optics in smoothly varying anisotropic planar structures: Application to liquid-crystal twist cells", J. Opt. Soc. Am. 63, 1374 - 73.
Da H. X. Xu C. Li Z. Y. 2005 "Magneto-optical effect of left-handed material", Eur. Phys. J. B 45, 347 - 74.
Yao-Yu Wang Lien-Wen Chen 2006 "Tunable negative refraction photonic crystals achieved by liquid crystals ", Optics Express 14, 10580 - 75.
Lakhtakia A. Messier R. 2005 "Sculptured Thin Films: Nanoengineered Morphology and Optics", SPIE Press, Bellingham, WA, USA,. - 76.
Da H. X. Xu P. Wu J. C. Li Z. Y. 2008 "Electrically controlled Kerr effect in magnetophotonic crystals based on nematic liquid crystals", J. Appl. Phys. 104, 033911 - 77.
Hai-xia Da. Zi-qiang Huang. Li Z. Y. 2009 "Voltage-controlled Kerr effect in magnetophotonic crystal", Opt. Lett. 34, 356 - 78.
Hai-xia Da Zi-qiang Huang Li Z. Y. - 79.
Chun-Yeol You Sung-Chul Shin 1997 "Novel method to determine the off-diagonal element of the dielectric tensor in a magnetic medium", Appl. Phys. Lett. 70, 2595 - 80.
Zvezdin A. Kotov V. 1997 Modern Magneto-optics and Magneto-optical Materials (IOP, Bristol,). - 81.
Balasubramanian K. Marathay A. Macleod H. A. 1988 "Modeling Magneto-optical thin-film media for optical-data storage", Thin Solid Films 164, 391 - 82.
Kee C. S. Park J. Y. Kim S. J. Song H. C. Kwon Y. S. Myung N. H. Shin S. Y. Lim H. 1999 "Essential parameter in the formation of photonic band gaps", Phys. Rev. E 59, 4695 - 83.
Kee C.S. Kim J.E. Park H. Y. Lim H. 1999 "Roles of wave impedance and refractive index in photonic crystals with magnetic and dielectric properties", IEEE Trans. Microwave Theory Tech. 47, 2148 - 84.
Saib A. Vanhoenacker-Janvier D. I. Huynen Laboratoire Encinas A. Piraux L. Ferain E. Legras R. 2003 "Magnetic photonic band-gap material at microwave frequencies based on ferromagnetic nanowires", Appl. Phys. Lett. 83, 2378 - 85.
Tarkanyan R. H. Niarchos D. G. 2006 "Effective negative refractive index in ferromagnet-semiconductor superlattices", Optics Express 14, 5433 - 86.
Khoo I.C. 1995 - 87.
Liebert L. 1978 - 88.
Yeh P. 1979 "Electromagnetic propagation in birefringent layered media", J. Opt. Soc. Am. 69, 742