## Abstract

We describe the main results previously studied concerning the thermal and electrical tuning of photonic band gap structures and the temperature-dependent defect modes in multilayer photonic liquid crystals using nematic liquid crystal slabs in a twisted configuration. In addition to this, we present new results regarding the electrical control of defect modes in such multilayer structures. In order to achieve this goal, we establish and solve numerically the equation governing the twisted nematic configurations under the action of the external electric field by assuming arbitrary anchoring conditions at the boundaries. After this, we write Maxwell's equations in a 4 × 4 matrix representation and, by using the matrix transfer technique, we obtain the transmittance and reflectance for incident circularly polarized waves.

### Keywords

- photonic band bap
- electrical and thermal tuning
- nematic liquid crystal
- multilayer structure
- defect mode

## 1. Introduction

Photonic crystals (PCs) are artificial structures with spatially periodic dielectric permittivity whose interesting optical properties have attracted the attention of the scientific community since the seminal works made by Yablonovitch [1] and John [2]. The most attractive attribute of these periodic structures is the existence of photonic band gaps (PBGs) in which the propagation of electromagnetic waves is prohibited for a specific wavelength range. In one-dimensional PCs, this phenomenon is usually called Bragg reflection. Liquid crystals (LCs) are anisotropic intermediate phases between the solid and liquid states of matter that possess positional and orientational order just like those of the solid crystals, and they can flow as a conventional liquid [3]. For many decades, LCs have been used as optoelectronic substances on account of easy tunability of their properties under the stimuli of external agents as temperature, pressure and electromagnetic fields. This fact suggests the conception of new artificial structures by making a convenient combination of LCs with PCs, whose most prominent feature is the externally controlled PBG. First studies reported on this subject in three- and two-dimensional structures are attributed to Busch and John [4] and Leonard et al. [5], respectively. In [4], it was demonstrated the tunability of the PBG under the action of an external electric field meanwhile the temperature tuning when a nematic LC is infiltrated into the void regions of solid PCs is showed in [5].

Multilayer photonic liquid crystals (MPLCs) consisting of LCs alternated by transparent isotropic dielectric films have been previously studied. In Ref. [6], Ha et al. demonstrated experimentally simultaneous red, green and blue reflections (multiple PBGs) using the single-pitched polymeric cholesteric LC films. Later, Molina et al. [7] investigated the strong dependence of electric field on the PBG for incident waves of left- and right-circular polarization at arbitrary incidence angles using nematic liquid crystal (NLC) slabs in a twisted configuration. In a recent paper, Avendaño and Reyes [8] studied the optical band structure for reflectance and transmittance considering that the dielectric matrix of a similar one-dimensional photonic structure to that studied in Ref. [7] depends on temperature and wavelength. Twisted nematic LCs, where the molecular orientation exhibits a 90° twist, have proven technological advantages to control light flow. They have been used to switch effectively the pass of polarized light in nematic displays by means of a normally applied low-frequency electric field.

Surface anchoring plays an essential role in the science and technology of LCs. The structure of LCs in the bulk is different than that near the interface, and the boundary conditions established from this interface structure influence the behaviour of the LCs in the bulk. There are two cases of surface anchoring of particular interest. First, a strong anchoring case in which the molecules near the surface adopt a rigidly fixed orientation, and the anchoring energies are very large. Second, a weak anchoring case where the surface strengths are not strong enough to impose a well-defined molecular orientation at the interface, and the expression for the anchoring energy is some finite function that depends on the LC properties at the surface, the surface properties and the external fields (e.g., electric and magnetic fields) and temperature [9, 10]. Anchoring effects on the electrically controlled PBG in MPLCs were previously investigated by Avendaño [11]. They considered a generalization of the model studied in [7] for which arbitrary anchoring of the nematic at the boundaries is taken into account. They also found the nematic configuration versus the anchoring forces and the PBG under the action of a strong enough external field parallel to the periodicity axis, which is able to modify the configuration of the nematic-twisted LC in the whole material including at the boundaries of each nematic slab. Later, Avendaño and Martínez [12] theoretically exhibited that this system is able to produce an omnidirectional PBG that can be electrically controlled for circularly polarized incident waves. An omnidirectional PBG requires that there be no states in the given frequency range for propagation in any direction in the material for both polarizations, which implies the total reflectivity for all incident angles.

Resonant transmittance peaks in the PBG can be induced in PCs when defects are introduced in the periodic lattice. In this case, standing waves with a huge energy density are localized in the proximity of the defects (defect modes). Ozaki et al. [13] developed the first tunable PC/LC hybrid structure by using a planar aligned NLC as a defect layer sandwiched between two one-dimensional periodical multilayers (dielectric materials with high- and low-refractive-index layers stacked alternatively) and demonstrated the electrical tuning of the defect modes. Thermal tunability of one-dimensional PC/LC cells was demonstrated by Arkhipkin et al. [14]. Electrical-dependent defect mode in PC/LC hybrid structures using a twisted nematic LC as defect layer was studied by Lin et al. [15] and Timofeev et al. [16]. Thermal tuning of defect modes in MPLCs using twisted nematic LC was recently shown in Ref. [8].

In this chapter, we describe the main results previously studied concerning the thermal and electrical tuning of PBG and the temperature-dependent defect modes in MPLCs using NLC slabs in a twisted configuration. In addition to this, we present new results regarding the electrical control of defect modes in MPLCs. In order to achieve this goal, we establish and solve numerically the equation governing the twisted nematic configurations under the action of the external electric field by assuming arbitrary anchoring conditions at the boundaries. After this, we write Maxwell's equations in a 4 × 4 matrix representation and, by using the matrix transfer technique, we obtain the transmittance and reflectance for incident circularly polarized waves.

## 2. Nematic-twisted configuration

As said above, we are focused on a 1D structure consisting in *N* NLC slabs in a twisted configuration alternated by *N* transparent isotropic dielectric films as it is illustrated in Figure 1(a) and (b). For each of the NLC cells, the nematic is sandwiched between two dielectric layers in such a way that its director is aligned parallel in both frontiers. A twist is then imposed on the NLC by rotating an angle 2*φt*, one of the dielectric layers about its own normal direction. Because of the possibility of molecular reorientation under the influence of external stimuli, such as electromagnetic fields, the director **n** takes the general form

with *α*(*z*) and *φ*(*z*), the polar (zenithal) and azimuthal angles made by **n** with the *xy* plane and the *x*-axis, respectively. For the present physical system and assuming small distortions in the nematic [17], the expression that describes the elastic energy density of the nematic is

where the positive elastic moduli *K*_{1}, *K*_{2} and *K*_{3} refer to splay, twist and bend bulk deformations, respectively. At this point, it is important to mention that it has experimentally found that when the nematic temperature is increasing to its transition temperature *TNI* (where the NLC becomes isotropic), a reduction of the values of the elastic moduli is induced [3, 17]. On the other hand, if we take into account the presence of an external electric field **E**, the interaction of this field with the LC is described through the energy density **D** = *ε*_{0}**ε** ⋅ **E** characterized by the uniaxial dielectric tensor

with *εa* = *ε*_{∥} − *ε*_{⊥} the dielectric anisotropy of the medium and *ε*_{0} the permittivity of free space. Here, *δij* is the Kronecker delta, *ε*_{⊥} and *ε*_{∥} denote the relative dielectric permittivity perpendicular and parallel to the nematic axis, respectively, and they are related to the ordinary *no* and extraordinary *ne* refractive indices by

The study of confined nematic liquid crystals is strongly influenced by the physical properties of the boundary walls [18]. From the macroscopic-geometrical and microscopic interactions between the molecules of such surfaces and of the nematic, the alignment of the director **n** on the boundary surfaces, known as *anchoring*, can be completely determined. Once the anchoring conditions are established, the orientation of the NLC molecules at the substrate surface determined the director in the bulk.

Several methods and techniques for surface alignment have been developed [19, 20]. In the case of rubbed polymer films [21, 22], it has been observed that NLC molecules are *strongly anchored* at the surface, and the alignment is parallel to the grooves produced by the rubbing process. Also, the orientation of NLC molecules at the surfaces is preserved even if an external field (electric or magnetic) is applied and removed.

On the other hand, photoalignment [23] and nanostructuring polymer surfaces [24] are contact-free methods where it is induced a surface ordering that causes an anchoring of controllable strength, which corresponds to a *weak anchoring*. For this anchoring condition, alignment of the NLC molecules before and after the application of external fields is different.

Anchoring energy can be expressed in terms of the surface anchoring coefficients which are related to the interaction strength between the NLC and the wall substrate for the deviation of the easy axis along the correspondent directions. It is experimentally found that these coefficients are temperature dependent [18] and their values for specific NLCs can be obtained by using the dynamic light scattering [25]. Thus, if we write the director in terms of *α*(*z*) and *φ*(*z*), as in expression (1), the anchoring energy of each NLC slab can be expressed in terms of the surface anchoring coefficients *Wα* and *Wφ* [26] as follows:

which is an extension of the Rapini-Papoular model [9] and where *φt* the twist angle. These anchoring coefficients are measured in energy per area units.

Thus, strong anchoring conditions are achieved when the anchoring coefficients are sufficiently large and can be modelled by considering that *Wα*→∞ and *Wφ*→∞. In contrast, for weak anchoring conditions, it is taken that *Wα*→0 and *Wφ*→0. Another criterion to establish whether the anchoring is strong or weak is based on the *extrapolation length* [27].

Thus, the equations governing the equilibrium configuration of the system are obtained by considering specific anchoring conditions and by minimizing the total free energy

that can be achieved by considering strong or weak anchoring conditions.

### 2.1. Strong anchoring

Let us assume that the structure shown in Figure 1(a) is subjected to a DC electric field **E***dc* = (0, 0, *Edc*) parallel to *z*-axis, and we consider that the orientation of the director at the surfaces of each nematic cell are fixed and given by *α* = 0°, *m* = 1 , 2 , 3 , … , *N*. Here, *N* nematic layers, respectively. Under these circumstances and by using a standard variational calculus procedure, the minimum free-energy condition *δF* = 0 together with the restriction *δ***n** = 0 at the surface of each slab generate the equations [12]

where we have defined the dimensionless parameter *f*(*α*) and *g*(*α*) are defined as

In absence of the dc electric field, the polar angle *α*(*z*) = 0° for any value of *z* and, the solution of Eqs. (6) and (7) is simply

where *φm*(*z*) represents the configuration of the *m*th layer in the region (*m* − 1)(*d* + *h*) ≤ *z* ≤ (*m* − 1)(*d* + *h*) + *d* and the nematic director (1) is reduced to

### 2.2. Weak anchoring

In this case, we consider a free-end-point variation for which the director orientation is affected by the existence of finite anchoring coefficients [12]. This minimization procedure leads to the same set of coupled equations given by (6) and (7) subjected to boundary conditions at each layer:

With *Γα* = 1/*γa*,*Γ* = *γφ*/*γa*, *γa* = *Wαd*/*K*_{1} and *γφ* = *Wφd*/*K*_{1}.

## 3. Electromagnetic propagation in a layered medium

The interaction between electromagnetic fields and matter is governed by Maxwell's equations and their corresponding constitutive equations. Optical propagation in layered media can be studied by conveniently writing Maxwell's equations in a 4 × 4 matrix representation. In this matrix representation, the boundary conditions of waves impinging on material can be imposed in such a way that the transfer and scattering matrix formalism to obtain the transmittances and reflectances can be used in a natural way [28, 29].

### 3.1. 4 × 4 matrix representation

In systems where boundary conditions cannot be avoided, Maxwell's equations require the continuity of tangential components of electric **E** and magnetic **H** fields at the boundaries. In studying the optical properties of dielectric layers which are confined between parallel walls, it is useful to write the set of Maxwell's equations in a representation where only appears, at the same time, the transversal components of **E** and **H** (two components for **E** and two components for **H**). This formalism is frequently referred to as Marcuvitz-Schwinger representation [30]. If we consider that the optical properties of a multilayer structure depends only on spatial variable *z*, we define the time-harmonic transversal four-vector

with *ω* the angular frequency of the propagating wave and *kx* the transversal component of the wave vector. Maxwell's equations, inside a nonmagnetic medium, can be written in the following matrix form:

for which the 4 × 4 matrix **A**(*z*) is given by

where *εij* (*i* , *j* = *x* , *y* , *z*) represents the elements of dielectric matrix in the structure, *k*_{0} = 2*π*/*λ* = *ω*/*c* is the wavenumber in free space, *λ* is the wavelength and *c* denotes the speed of light in vacuum. Also the fields **e**(*z*) = (*ex*(*z*), *ey*(*z*), *ez*(*z*)) and **h**(*z*) = (*hx*(*z*), *hy*(*z*), *hz*(*z*)) are related to the electric **E**(*z*) and magnetic **H**(*z*) fields by the following expressions *ε*_{0} and *μ*_{0} the permittivity and permeability of free space, respectively.

For a homogeneous and isotropic dielectric medium, the matrix **ε**(*z*) is diagonal and independent of the position, whereas for a nematic slab, **ε**(*z*) depends on the local orientation of the principal axis of the liquid crystal molecules characterized by expression (3).

### 3.2. Boundary condition

Let us consider a multilayer structure where each of the layers is confined between two planes, and the whole structure is surrounded by air. An electromagnetic wave impinging from the left side of the multilayer structure will propagate through the sample, and it will be transmitted and reflected outside the medium (see Figure 1 (c)).

The general solution of the differential equation (16) for electromagnetic waves propagating in homogeneous media is the superposition of four plane waves: two left-going and two right-going waves. With this in mind, we state the procedure to find the amplitudes of the transmitted and reflected waves in terms of incident ones (at plane *z* = 0). This implies the definition of the following quantities [31]:

(i) The propagation matrix **U**(0, *z*) that is implicitly defined by the equations

where **U**(0, *z*) satisfies the same propagation equation (16) found for **ψ**:

This propagation matrix gives the right-side field amplitudes of the multilayer structure as function of the left-side ones.

(ii) For a specific value *z* = *z*_{0}, the transfer matrix is defined as **U**(0, *z*_{0}).

(iii) The scattering matrix **S** giving the output field as function of the incident one. The matrix **S** is defined through the relation **α***out* = **S** ⋅ **α***in*, where **α***in* and **α***out* are the amplitudes of the in-going and out-going waves.

To find out **S**, the field must expressed, in any one of the external media, as a superposition of planes waves by setting:

where

The relation **ψ** = **T** ⋅ **α** can be interpreted as a basis change in the four dimensional space of the *state vectors* **ψ**. The columns of **T** are the **ψ** vectors representing the four plane waves generated by the incident waves in the two external media (assumed as identical). The elements of vector **α** are the amplitudes of the four plane wave. The choice of the new basis could be different depending on the particular problem. By setting

the scattering matrix writes:

where the symbols + and *f* (− and *b*) mean forward (backward) propagating waves.

We point out that the methods of transfer and scattering matrices are very useful in studying the plane wave transmission and reflection from surfaces of multilayer structures.

Differential equation (16) can be formally integrated over a certain distance *z*_{0} of the medium

and by straight comparison of Eqs. (18) and (23), the transfer matrix **U**(0, *z*_{0}) is defined as:

where plane waves are incident and reflected in the half-space *z* < 0, and plane waves are transmitted on the half-space *z* > *z*_{0}.

It can be seen immediately that the problem of finding **U**(0, *z*_{0}) is reduced to find a method to integrate expression (24) on the whole multilayer structure. Because of the non-homogeneity of the medium proposed here, we consider it as broken up into many very thin parallel layers, each of them with homogeneous anisotropic optical parameters [32]. In this way, **U**(0, *z*_{0}) is obtained by multiplying iteratively the matrix for each sublayer from *z* = 0 to *z* = *z*_{0}.

### 3.3. Transmission and reflection by multilayer structures

As said above, the general solution of the differential equation (16) for electromagnetic waves propagating in homogeneous media is the superposition of forward and backward propagating waves. The obliquely incident and reflected electromagnetic fields in free half-space *z* ≤ 0 (Figure 1(c)), for an arbitrary polarization state which are solutions of equation (16), can be expressed as:

and

where **k**_{0} = (*k*_{0x}, *k*_{0y}, *k*_{0z}) = *k*_{0}(sin*θ*, 0, cos*θ*) is the wave vector of the incident wave making an angle *θ* with respect to the *z*-axis, *aL* and *aR* represent the amplitudes of left- and right-circularly polarized (LCP and RCP) components of incident wave, respectively, and *rL* and *rR* correspond to those of the reflected wave (see Figure 1(c)). The unit vectors **u** and **v** are defined as

with **u**_{x} , **u**_{y} , **u**_{z} the triad of Cartesian unit vectors. In the region *z* ≥ *z*_{0}, the transmitted electromagnetic field is

where *tL* and *tR* are the amplitudes of LCP and RCP components, respectively, of transmitted wave. As the tangential components of **e** and **h** must be continuous across the planes *z* = 0 and *z* = *z*_{0}, the boundary values **ψ**(0) and **ψ**(*z*_{0}) can be fixed as:

with

By using Eqs. (23), (24) and (29), the problem of reflection-transmission can be established as follows

where **M** = **P**^{−1} ⋅ **U**(0, *z*_{0}) ⋅ **P** and **U**(0, *z*_{0}) are defined in (24). Notice that the matrix equation (31) gives a set of coupled equations relating the amplitudes *aL*, *aR*, *rL* and *rR* (from *z* ≤ 0) to the transmitted amplitudes *tL* and *tR* (for *z* ≥ *z*_{0}).

The scattering matrix **S** relates the amplitudes *tL*, *tR*, *rL* and *rR* to the known incident amplitudes *aL* and *aR*. This relation can be expressed in terms of matrix **M** as [33]

where

and

Co-polarized coefficients have both subscripts identical meanwhile cross-polarized coefficients have different subscripts. The square of the amplitudes of *t* and *r* is the corresponding transmittance and reflectance, respectively; thus, *TRR* = |*tRR*|^{2} is the co-polarized transmittance corresponding to the transmission coefficient *tRR*, *TRL* = |*tRL*|^{2} is the cross-polarized transmittance corresponding to the transmission coefficient *tRL*, and so forth. In the absence of dissipation of energy inside the sample, the principle of conservation of energy must be satisfied from which we have that

Before ending this section, we mention that an alternative way to find the transmission and reflection coefficients is using the expressions given by (21) and (22). Also, the system of equations (31) can be solved numerically to find the scattering matrix.

## 4. Numerical results and discussion

In previous sections, we have presented in detail a general mathematical formalism to determine the reflectances and transmittances by multilayer structures. In this section, we apply this formalism to MPLCs using NLC slabs in a twisted configuration considering that circularly polarized light impinges on the structure in order to analyse the optical spectra and their dependence on external agents. In particular, we describe the main results previously studied concerning the thermal and electrical tuning of optical spectra and the temperature-dependent defect modes. In addition to this, we present new results regarding the electrical control of defect modes.

### 4.1. Electrical tuning of band structure and defect mode

In this section, we present the influence of the electric field on the optical band structure and defect mode by considering arbitrary anchoring conditions at the boundaries. To this aim, the equilibrium configuration of each NLC layer as a function of *σ* is obtained by solving the second order differential equations (6) and (7) for *α*(*z*) and *φ*(*z*) subjected to the conditions expressed in Eqs. (11)–(14). Then, this configuration is substituted into Eq. (23) in order to obtain the transfer matrix **M** as function of *σ* for circularly polarized incident waves.

Numerical calculations were performed by considering a NLC phase 5CB for which *K*_{1} = 0.62 × 10^{−11}*N*, *K*_{2} = 0.39 × 10^{−11}*N*, *K*_{3} = 0.82 × 10^{−11}*N* [17] and refractive indices at optical frequencies *φt* = 90°, and the homogeneous isotropic dielectric medium is zinc sulphide (ZnS) with refractive index *nd* = 2.35. The MPLC consists of *N* = 11 NLC layers alternating with *N* = 11 dielectric slabs with the same thickness. Finally, we report our results parameterizing all the spatial variables by the NLC thickness *d*. In this way, the dimensionless thickness of each NLC cell is *h*^{'} = *d*/*d* = 1, whereas for each ZnS slab is *h*^{'} = *h*/*d* = 1, and so forth.

Due to the competition between orientation produced by influence of the external electric field and by surface anchoring effects, we expect a deformation in the NLC only above a certain critical value *σc*. This critical electric field is expected to be maximum for the case of strong anchoring conditions, whereas for the weak anchoring case, *σc* will decrease as the surface forces get smaller [7, 11].

#### 4.1.1. Strong anchoring conditions

For strong anchoring conditions, the orientation of the nematic molecules at the walls of each NLC is specified in Section 2.1. The curves for *α*(*z*) and *φ*(*z*) are shown in Figure 2(a) and (b), respectively, as function of dimensionless variable *w* = *z*/*d* above the critical value *σc* = 3.26. As it can be noticed in Figure 2(a), an increment in the electric field involves the augmentation in the polar angle *α*. Owing to the influence of the external field, the nematic molecules tend to be aligned parallel to it (*z*-axis). As expected, for *σ* < *σc*, *α* = 0° for all values of *w*, which means that the director in this case is perpendicular to the *z*-axis. In Figure 2(b), we observe that for *σ* < *σc*, the curves for azimuthal angle *φ* are reduced to straight lines with slope equal to 2*φt* = 90°, that corresponds to the configuration of a pure twisted NLC. Above the critical value, the strong anchoring condition is really dominant on the parameter *φ* as the electric field increases. Indeed, most of molecules tend to spread far from *xz* plane.

Figure 3 exhibits the co-polarized and cross-polarized transmittances and reflectances for LCP and RCP waves impinging normally on the structure as function of the dimensionless parameter *d*/*λ* for continuous values of the electric field above the critical value *σc* (it is worth to mention that below this value, the not-shown curves are very similar to that of *σ* = *σc*). Note the strong influence of *σ* on the transmission and reflection spectra in enhancing and extinguishing bands. Indeed, Figure 3 clearly shows that for *σ* = *σc*, the curves for transmittances exhibit several stop bands of different widths in the plotted interval and, as *σ* increases, each stop band gets wider for co-polarized transmittances *TRR* and *TLL*. Also, cross-polarized transmittances *TLR* and *TRL* are totally absent for high enough values of the electric field. On the other hand, at the critical value, co-polarized reflectances *RRR* and *RLL* exhibit narrow-reflection bands with relatively high amplitudes and cross-polarized reflectances *RLR* and *RRL* show one-dominant high-reflection band. In this case, co-polarized reflectances reduce their band amplitudes practically to zero, and reflection bands of cross-polarized reflectances are highly enhanced for larger values of *σ*. These optical properties allow us to use this MPLC as an electrically shiftable universal rejection filter for incident RCP and LCP waves where, by increasing the electric field, one can highly enhance the cross-polarized reflection bands and supress the co-polarized ones.

In [7], it is shown that for a fixed value of *σ* the band structure of the reflectances and transmittances are shifted towards smaller wavelength regions as the incident angle *θ* increases. This behaviour results from the fact that for plane electromagnetic waves propagating obliquely with respect to the layer interfaces, only the normal component of the wave vector is involved in the photonic band formation. Hence, as the incident angle augments, the relative position of the bands is moved towards smaller wavelengths.

If one of the layers possesses a different size compared with the remaining ones, this layer can act as a defect, and an optical defect mode can be induced. Here, we specifically consider that the middle NLC-ZnS stack of the MPLC has a different size compared with the remaining ones. We choose specific values *dd* = 2*d*^{'} and *hd* = 2*h*^{'}, where *dd* and *hd* are the dimensionless thicknesses of the NLC and ZnS defect layers, respectively. Figure 4(a) and (b) displays the defect mode induced in the photonic band of the co-polarized transmittance *TRR* and cross-polarized reflectance *RLR*, respectively, by LCP waves impinging normally on the MPLC. We notice that as the parameter *σ* increases two important facts occur: (i) two defect modes with small amplitude are induced within the first stop band (see Figure 3) which gradually merge into only one; the position of the defect mode possessing the largest wavelength moves toward regions of smaller wavelengths, keeping fixed the position of the other one and (ii) the amplitude of the defect modes gets larger. Physically, the origin of the defect mode is the phase change due to the variation in the optical path length caused by the defective medium. Once the defect mode is created at specific position, it can be controlled by inducing reorientation in the nematic molecules by means of an external electric field [7]. Indeed, since the refractive index of the LC depends on the angle *β* between the wave vector **k** of the electromagnetic wave in the LC and the local orientation of the director **n**, the refractive index (and the optical path length) can be changed by varying *β*. At normal incidence and *σ* < *σc*, *β* = 90° for all positions *z*. Nevertheless, as *σ* increases, most of the molecules tend to be aligned parallel to *z*-axis (see Figure 2(a)) and *β*→0°. These results show that the amplitude of defect mode and its position can be tuned by a DC electric field.

#### 4.1.2. Weak anchoring conditions

It is experimentally found that for a LC phase 5CB, the polar anchoring *γα* is of the order of 10^{1}, and this value is one or two orders stronger than the azimuthal anchoring *γφ* [34]. Under these considerations, the values of the dimensionless anchoring parameters are taken as Γ = Γ_{α} = 0.1.

The curves for *α*(*z*) and *φ*(*z*) are shown in Figure 5(a) and **(b)**, respectively, as function of dimensionless variable *w* = *z*/*d* above the critical value *σc* = 2.86. In Figure 5(a), we can notice that, as *σ* augments, the values of *α* increase, getting a maximum at the middle of the cell. Because of the influence of external electric field, the polar angle at both borders enlarges by increasing *σ* highlighting the fact that even at the borders, the field is able to distort the configuration. Figure 5(b) shows two interesting phenomena: (i) for *σ* < *σc*, the curves are reduced to straight lines with slope equal to 2*φ*_{0c}, where *φ*_{0c} represents the azimuthal angle adopted by the MPLC at the walls of each NLC cell for values of electric field below the critical field; (ii) above the critical value, most of the molecules tend to acquire an angle *φt* = − 45° for 0 < *w* < 0.5 and *φt* = 45° for 0.5 < *w* < 1.

Figure 6 shows the co-polarized and cross-polarized transmittances and reflectances for LCP and RCP waves impinging normally on the structure as function of the dimensionless parameter *d*/*λ* for continuous values of the electric field above the critical value *σc* (below this value, the not-shown curves are very similar to those corresponding to *σ* = *σc*). Although, the optical properties shown in Figure 6 are qualitatively similar to those of Figure 3 where strong anchoring conditions were considered, we notice that in the case of weak anchoring conditions, the behaviour of transmittances and reflectances in Figure 6 is enhanced in comparison with Figure 3. Because of the strong influence of the electric field on the molecular orientation for all values of *z* (including the walls of each cell), the alignment of most of the nematic molecules parallel to *z*-axis occurs at smaller values of electric field unlike for strong anchoring. Hence, the phenomenon of extinguishing and enhancing bands is present at smaller values of *σ*.

Now, we induce a defect mode in the photonic band structure by generating a defect in the MPLC in the same way as explained in Section 4.1.1. Figure 7(a) and (b) displays the defect mode induced in the photonic band of the co-polarized transmittance *TRR* and cross-polarized reflectance *RLR*, respectively, for LCP waves impinging normally on the MPLC. Similar to the case of strong anchoring conditions, we can observe that when the parameter *σ* augments the amplitude of the defect modes gets larger, and the position of the defect mode possessing the largest wavelength moves toward regions of smaller wavelengths, while the position of the other defect mode remains fixed. These facts are enhanced in comparison to those of strong anchoring assumptions because of the strong influence of the electric field on the molecular orientation for all values of *z*, including the walls of each NLC slab in the MPLC. This implies that the defect-mode amplitude gets larger for smaller values of *σ* in comparison to that of the strong anchoring case.

### 4.2. Temperature-dependent band structure and defect mode

Here, we assume that the orientation of the director at the surfaces of each nematic cell is strongly anchored at the boundaries. In order to obtain the band structure, we apply the same mathematical procedure as depicted in Section 4.1, but in this case, we have to take into account that the director **n** is given by expressions (9) and (10) and the elements of dielectric tensor **ε**(*z*) depend on the wavelength and temperature [8]. By considering E7 LC mixture slabs and ZnS dielectric layers, it is found that in the interval of temperatures [15°C, 50°C] and for RCP waves impinging normally, the position of the photonic bands of the co-polarized transmittance *TRR* and cross-polarized reflectance *RLR* (analogous to those of Figure 3 with *σ* = 0) can be shifted from regions of small wavelengths toward regions of higher wavelengths by increasing the thickness *d* of the NLC layers. In addition to this, bandwidth increases for thicker layers, and new narrower transmission bands are created in regions of smaller wavelengths. In summary, the position, the width, and the number of bands augment as the thickness *d* is increased. Physically, when the magnitude of *d* gets larger, the optical path lengths increase, and hence, the wavelength zones of destructive or constructive interference are shifted towards higher wavelength regions. For constant thickness and temperature, it is observed that as the incident angle *θ* augments, the photonic bands undergo a shift towards smaller wavelengths and their widths get narrower. As said previously, this behaviour results from the fact that for plane electromagnetic waves propagating obliquely with respect to the layer interfaces, only the normal component of the wave vector is involved in the photonic band formation. Thus, as the incident angle augments the relative position of the bands moves towards smaller wavelengths, and the overall band is always closed up. On the other hand, for constant thickness and ~~a~~ fixed incident angle, as the temperature augments, the photonic bands move towards the short-wavelength region. Physically, since the average refractive index of the liquid crystal decreases as the temperature gets larger, the optical path length diminishes, and thus, the wavelength regions where the waves are able to undergo constructive or destructive interference shift towards smaller wavelengths zones.

In a similar way as demonstrated above, a defect mode can be induced by considering that the middle layer of the homogeneous and isotropic slabs (ZnS) has a different size compared with the remaining ones. If we consider normal incident RCP waves, for *hd* = 1.991*h*^{'} and temperature values in the interval [15°C, 50°C], the position of the defect mode induced in the photonic band of the co-polarized transmittance *TRR* and cross-polarized reflectance *RLR* (analogous to those of Figure 4 with *σ* = 0) shifts from larger wavelengths toward smaller ones as the temperature gets increasing. Because the origin of the defect mode is the phase change due to the variation in the optical path length caused by the defective medium, the defect wavelength can be shifted towards smaller wavelength regions as the temperature is increased by taking into account that the average refractive index of the NLC decreases as temperature increases.

## 5. Conclusion

We presented a series of results concerning the thermal and electrical tuning of photonic band gaps and defect modes in multilayer photonic liquid crystals consisting of liquid crystal layers alternated by transparent isotropic dielectric films using nematic liquid crystal slabs in a twisted configuration. We exhibited that the position and width of the band gaps can be electrically and thermally controlled. When one of the homogeneous and isotropic slabs has a different size compared with the remaining ones, a defect mode is induced in the band structure whose wavelength can be tuned. Tuning of the transmission and reflection bands and the defect mode investigated here could be useful in the implementation of tunable optical filters and waveguides.