Calculation of reflectance and transmittance of optical birefringent networks based on cholesteric liquid crystals

In this paper we calculated reflectance and transmittance of different optical birefringent networks sandwiched between two isotropic media. To model optical phenomena in the considered systems, we applied the 4×4 matrix method. Some interesting reflectance and transmittance spectra and polar plots for different system parameters and arbitrary incident monochromatic light were reported. The illustrated and discussed results can be useful for understanding optical phenomena in numerous birefringent media. Especially, the reflectance and transmittance spectra computed in this article can be suitable for analysis of more advanced contemporary optical systems, i.e. photonic band gap materials based on cholesteric liquid crystals.


INTRODUCTION
Birefringence is the optical property of materials having refractive indices, which depend both on the polarization and propagation direction of the incident light.The best characterized materials possessing these properties are crystals, because their refractive indices are well defined due to their specific structures.Among the mentioned crystalline structures, special importance have the so-called liquid crystals (LCs).Liquid crystals are the interesting fourth state of matter that lies between the crystalline solid and amorphous liquid states.In some temperature ranges, these materials exhibit both the properties of liquids and the properties of crystalline solid state.Most of them are organic compounds with elongated molecules that influence their physical properties.
The classification of LCs includes three types of liquid crystalline phases, i.e. nematic, smectic, and cholesteric liquid crystals (CLCs).The abovementioned groups differ in physical properties, especially in optical ones.However, due to the strongly periodic helical structure, cholesteric liquid crystals are characterized by unique optical properties.First of them is the abovementioned optical birefringence, i.e. this substance have two different refractive indices.It means that in such an optical media, the propagating wave of the monochromatic light decomposes into two waves, which propagate with different velocities.CLCs also exhibit optical activity and circular dichroism.However, the most popular engineering applications of them are caused by selective reflection of the incident light.The color (wavelength) of the reflected light depends on mechanical stress, temperature, external electric and magnetic fields, etc.
Due to the abovementioned properties, LCs have been used in practical applications for many years, for instance in liquid crystal displays (LCDs).Nowadays, most of phones, computers, television sets and other electronic devices implement LCDs as a one of the most important optical display systems for high-performance flat panel display of information.This is mainly arises from such advantages as flat panel, large area, light weight, high definition, low driving voltages and low power consumption (Yeh and Gu, 1999).However, CLCs seem to be also the promising candidates for numerous different photonic applications.This is why, currently, a lot of attention is paid to the photonic crystals (PCs) having a three-dimensional (3D) ordered structure with a periodicity of the optical wavelength, whose optical properties, due to their periodic dielectric structures, are similar to the properties of LCs.In such materials, novel physical concepts have been predicted theoretically and different engineering applications of PCs have been proposed (Ozaki et al., 2007;Cardador et al., 2017;Escobar et al., 2018).
Due to the abovementioned advantages, PCs have attracted much attention from both fundamental studies and practical applications in optical circuits, devices and interconnects (Guo et al., 2010;Notomi, 2010;Prasmusinto et al., 2017).Among the most known new applications of PCs, one can mention light-emitting diodes (Wierer et al., 2009;Wiesmann et al., 2009), optical switches (Belotti et al., 2008;Nozaki et al., 2014), colour filters (Huang and Zhang, 2011;Lin et al., 2015), gas sensors (Hodgkinson and Tatam, 2013;Xu et al., 2013;Cardador et al., 2017) and lasers (Noda, 2010;Hirose et al., 2014;Ortega et al., 2018).Recently, photonic crystals find applications in building all-optical computers and communication systems, using photonic crystal logic gates for faster computers and mobile devices with high both speed processing of data and transmission rates.An interesting literature review in this area can be found in the review paper (Hussein et al., 2018).Also, the most recent articles pay much attention to the theoretical study and new applications of PCs.For instance, Cardador et al. (Cardador et al., 2017) designed and optimized theoretically a macroporous silicon photonic crystal for its use in gas sensing applications and IR optical filters, confirming them as promising candidates for precise optical gas sensors.Tavousi and Heidarzadeh (Tavousi and Heidarzadeh, 2018) realised a high efficiency multichannel drop filter by using a square-lattice-type photonic crystal and a circular ring resonator.Mandal (Mandal, 2018) demonstrated a visible frequency plasmonic perfect absorber made of a thin metal layer with an array of cylindrical grooves used to achieve perfect absorption at tunable wavelengths.Escobar et al. (Escobar et al., 2018) fabricated an active 1D photonic waveguide of porous silicon capable of concentrating light emission, which have potential applications in different integrated optoelectronic devices.In other paper, Zhang et al. (Zhang et al., 2018) proposed a novel compact photonic crystal ring resonator configuration to realise high-efficiency waveguided add-drop filtering.
In general, photonic crystals contain a periodic distribution of both refractive indices in one, two, or three dimensions, and can be used to prohibit, confine and control light propagation in a specific wavelength bandwidth (Lin et al., 2013;Cardador et al., 2017).One of the most important forms of PCs are cholesteric liquid crystals (CLCs) that spontaneously form periodic helical structures (molecules are self-assembled).If the pitch of the CLC helix is of the order of magnitude of the wavelength of the visible light, the material can exhibit photonic properties.In this case, CLCs present a photonic band gap.It means that circularly polarized light with the same handedness as that of the helix cannot propagate in a certain range of frequencies.For instance, the photonic character of CLCs is the basis for the application of these materials as low-threshold distributed feedback lasers (Ortega et al., 2018).
The mentioned helical structure of PCs leads to selective reflection for circularly polarized light in different optical systems (Morris et al., 2008;Choi et al., 2009;Balamurugan and Liu, 2016;Nemati et al., 2018;Ortega et al., 2018).The birefringence is typically limited to 0.3-0.4 for colourless organic compounds, and therefore the bandwidth is usually less than 100 nm in the visible spectrum (Balamurugan and Liu, 2016).However, there are also materials with different bandwidths.For instance, Fleming et al. (Fleming et al., 2002) reported measurements and simulations of a 3D tungsten crystal that has a large photonic bandgap at infrared wavelengths (from about 8 to 20 µm).An interesting review of the methods of fabrication of various photonic band gap materials based on cholesteric liquid crystals are presented in one of the recent review papers (Balamurugan and Liu, 2016).Especially, that review paper covers the methodology used recently, when the fabrication of photonic band gap materials arose.
It should be noticed that optical properties of LCs strongly depend on the orientation of the liquid crystal molecules.By applying voltage across the transparent electrodes, an electric field obtained inside the liquid crystal can be used to control the orientation of these molecules.This leads to a change in the property of the considered optical systems (Yeh and Gu, 1999).In turn, the performance of photonic crystal devices depends strongly on their geometry, mechanical stress and defects inside (Hart et al., 2012).An interesting literature review on the important defect types, their origins as well as interactions during the bulk crystal growth from the melt and selected epitaxial processes is given in one of the recent review papers (Rudolph, 2016).Some aspects of the effect of anisotropy on defect mode features in CLCs was investigated by Gevorgyan (Gevorgyan, 2018).Yang et al. (Yang et al., 2015) investigated the effect of stress-induced anisotropy on localised mode in photonic crystal and they found that the change in the crystal thickness, caused by stress, leads to localised modes to move linearly in the transmission spectrum.Also recently, Avendanno and Martinez (Avendanno and Martinez, 2018) investigated the transmission spectra in a 1-dimensional dielectric multilayer photonic structure, which contains a cholesteric liquid crystal elastomer layer as a defect.It should be emphasized that although the issues mentioned here are really interesting, they are not considered in detail in this study.crystals Dariusz Grzelczyk et al.
Latin American Journal of Solids and Structures, 2019, 16(1 Thematic Section), e105 3/20 In this paper, we focused on a theoretical analysis of optical phenomena in different birefringent media, with particular reference to the structure of the CLCs.In order to calculate the reflectance and the transmittance of the investigeted optical system, we implemented a computer algorithm which realizes the 4×4 matrix method.This method gives the exact solutions, since it takes into account both the effect of refraction and multiple reflections at the interfaces of the dielectric discontinuities.And although we investigated only the chosen birefringent networks, the implemented algorithm can be adopted to model optical phenomena of other optical systems, for instance other types of photonic crystals with a known distribution of the dielectric structure.
The rest of the paper is organized as follows.Section 2 summarizes the matrix methods applied to birefringent networks, especially the conventional Jones calculus, the extended Jones 2×2 matrix method and the Berreman 4×4 matrix method.The investigated optical system and crucial stages of the implemented matrix method are presented in Section 3. Section 4 contains calculations and discussion of the reflectance and the transmittance obtained for different optical systems.The article is concluded in the last Section 5 with a brief summary of the main obtained results.

METHODS -A BRIEF LITERATURE REVIEW
Optical phenomena and the transmission of light through birefringent optical media (networks) have been treated by different methods.In the conventional Jones calculus, each optical element (wave plate, liquid crystal layer, etc.) is represented by a 2×2 matrix, and reflection/refraction of the incident light at the plate surfaces (dielectric discontinuities) are neglected.Moreover, this method is limited to normally incident light and does not explain the leakage of off-axis light through a pair of crossed ideal polarizers (Gu and Yeh, 1999).For treating the transmission of off-axis light in a general birefringent optical media, the extended Jones matrix method should be used, which takes into consideration the single reflection at the interfaces.This method is adequate for numerous practical applications, and has been widely used in the analysis of numerous optical birefringent systems.Gu and Yeh (Gu and Yeh, 1999) applied the extended Jones matrix method in the analysis of compensators for liquid crystal displays.MacGregor (MacGregor, 1990) presented a model for computing the conoscopic figures produced by a uniformly oriented liquid crystal in a thin cell observed with a polarizing microscope.Lien (Lien, 1990a;1990b) used Jones matrix methods for the high pretilt twisted nematic cell and for the twisted nematic liquid-crystal display at oblique incidence.Ong (Ong, 1991a;1991b) analysed electro-optics of electrically controlled birefringence liquid-crystal and twisted nematic liquid-crystal displays at oblique angle of the incident light.In 2000, Li (Li, 2000) applied the Jones matrix method to ellipsometry, i.e. for investigating the dielectric properties of thin film.In 2010, Chen and co-investigators (Chen et al., 2010) used a new 2×2 matrix method and discussed the polarization state of transmitted waves through the cholesteric liquid crystals.Most recently, with the help of the concept of 2×2 transfer matrix method, Kumar et al. (Kumar et al., 2018) investigated comparison between dielectric-dielectric photonic crystal and birefringent-dielectric structure with and without the gradation in the thickness of the layers of the proposed structure.
On the contrary to the aforementioned methods dedicated to the birefringent networks, the exact solutions can be obtained by using the 4×4 matrix method, which takes into account both the effect of refraction and multiple reflections between plate interfaces (Gu and Yeh, 1999).This approach has been applied by many researches for different optical systems.For instance, Schwelb (Schwelb, 1986) analysed lossy gyroelectromagnetic layers in polar and longitudinal orientation.Chen et al. (Chen et al., 1996) employed the 4×4 matrix method to liquid crystal displays.Using the 4×4 matrix method, Ivanov and Sementsov (Ivanov and Sementsov, 2000) modelled the propagation of electromagnetic waves in stratified bianisotropic chiral structures.Yang et al. (Yang et al., 2015) also used this method to investigate the effect of stress-induced anisotropy on localized mode in photonic crystal.Most recently, the 4×4 matrix method were also used by Ortega and his co-investigators (Ortega et al., 2018) to study different kinds of cells of cholesteric liquid crystal lasers.Due to abovementioned anisotropic properties of CLCs and PCs, also in this paper we used the exact 4×4 matrix method, which is shortly presented in Section 3.

The 4×4 matrix method
In this study, we adopted the 4×4 matrix method on the basis of mathematical formalism presented in monograph (Yeh and Gu, 1999).Therefore, only the main and crucial stages of this method are presented below.
In this method, the dielectric tensors ˆ() n ε for the individual liquid crystal layers ( 1,2,..., n N  ) and bounding substrates ( 0 n  and 1 nN  ) are determined first.For both substrates where is the vacuum permittivity.In case of uniaxial liquid crystal medium, in each -th n layer, the orientation of the optical axes () n c of the molecules in the laboratory coordinate system can be described by two angles, i.e.
where x , y and z are the unit vectors of the laboratory coordinate system.For an arbitrary orientation of the optical axis In the second step, the wave vectors ) of all four elementary waves in all layers ( 0,1,2,..., 1 nN  ) are calculated.Let , the z components are determined directly from the wave equation in momentum space is the vector of the electric field, whereas  is the magnetic permeability (here the considered optical medium is a diamagnetic substance).The above presented wave equation ( 7) has the following scalar form In order to obtain non-trivial plane-wave solutions, the determinant of the 3×3 matrix must vanish in equation ( 8).
This condition gives an equation in  that yields four roots ) propagate backward.Two of these roots have always the opposite sign, but the same absolute values.These roots correspond to the ordinary waves, i.e. the wave with positive root propagates forward, whereas the wave with a negative root propagates backward.Two other roots (with different sign and, in general, different absolute values) correspond to the extraordinary waves.Here also the wave with a positive root propagates forward, whereas the wave with negative root propagates backward.
The third stage of the implemented method is devoted to computing the optical polarization vectors ) representing the directions of the electric field as well as the corresponding vectors representing the directions of the magnetic field.For each isotropic media, these vectors, for s and p waves, can be obtained as follow and In each -th n layer of the birefringent medium, these vectors, for o and e waves, have the following forms The vectors , which correspond to the directions of the magnetic field H of all elementary waves ( 1, 2,3, 4

 
) propagating in the layered anisotropic optical medium can be expressed as In the fourth stage of the method, transition matrices between individual layers are calculated, based on the assumption of continuity of the tangential components of the vectors at the dielectric interfaces.In the applied approach, we expressed the amplitudes  of all elementray waves in the upper substrate.This relation is given by the transmittion matrix 44 M  .With this assumption, the electric field E in the z coordinate can be expressed as the sum of the electric fields of all elementary waves, in the form where denote the amplitudes of all elementary waves.Similarly, the magnetic field H can be expressed as a sum of the magnetic fields of all elementary waves, in the form At the interface , both x and y components of electric and and magnetic fields ( x E , y E , x H and y H ) must be continuous, and it gives the set of the following four equations where ) is the thickness of -th n layer and equals dN , while 1 0 N t   .Equations ( 15) have the following scalar form where

x p x p x p x p y p y p y p y p D
x q x q x q x q y q y q y q y q , (17) and where The fraction of the electromagnetic energy reflected is while the fraction of the electromagnetic energy transmitted is which depend on the matrix elements ij M and on the polarization state of the incident s and p waves (i.e.
while the fraction of the energy transmitted is given by

RESULTS
In what follows we illustrated and discussed some calculations of the reflectance R (and transmittance T) spectra of the normally incident light, the reflectance R vs. the angle  , as well as colour polar plots in the full range of the angles  and  of the incident light.In the presented reflectance spectra, the colours visible under the chart correspond to the colours of the light with wavelength presented on the x-axis.In turn, in all polar plots presented in this paper, the radial position corresponds to the angle  (from 0, in the middle, to the 90°, at the edges of the plotted polar distributions), whereas the circumferential position corresponds to the angle  (from 0 to 360°, as it is shown in plotted distributions).In all cases, we used a monochromatic non-polarized light and different values of the parameters of the investigated system.To test the implemented method, we considered the propagation of light within the homogenous isotropic medium and the propagation of light in the homogenous isotropic medium sandwiched between the other isotropic medium first.Next, we investigated the reflection properties of liquid crystals with different textures, i.e. nematic, homeotropic and helical structures (cholesteric liquid crystal).All calculations were carried out for N =1000.

Propagation of light in a homogenous isotropic medium
To test the developed algorithm, we considered the case when the light propagates through the homogenous and lossless isotropic optical medium first, i.e. we assumed that  , in each case we obtained R = 0 and T = 1.This is why there are no dielectric discontinuities at the interfaces of the imaginary layers in such optical systems.Therefore, there is no electromagnetic energy reflected, i.e. the incident light is completely transmitted through the investigated optical system.Moreover, we tested the propagation of not normally incident light (i.e. for arbitrary angles  and  , which describe the orientation of the wave wector of the incident beam of the light).In these cases, we also obtained R = 0 and T = 1.

Propagation of light in a homogenous isotropic medium sandwiched between other isotropic medium
In this subsection we considered the propagation of light within an isotropic medium sandwiched between other isotropic medium (see parameters of the system in Tab. 1). Figure 2 presents reflectance spectra of the investigated system for d = 10 μm and different values of the refractive index s n of the substrates.It can be clearly seen that we obtained R + T = 1 in each case, and therefore usually only reflectance spectra are presented in the further analysis.The obtained results show that the reflection coefficient R is relatively low, and oscillates locally between the minimum and maximum values.This amplitude is greater when the difference between refraction indices is greater.In the presented cases, we have the situation similar to the interferences which occur in Fabry-Perot resonator, whereas the role of the optical mirrors is played by the dielectric interfaces at the surfaces of contacting optical media with different refraction indices.
Figure 3 presents reflectance spectra of the system for different values of the thickness d of the isoptropic medium.In all cases, the maximum value of the amplitude of oscillations is the same (R ≈ 0.03), since the difference between the refractive indices of both isotropic media are the same in all presented cases.However, it is clearly visible that for larger values of the thickness d, the reflection coefficient R oscillates with a higher frequency.

 
and depicted in Fig. 7.It should be noted that in this case, due to inhomegenous structure of the considered medium, the reflection coefficient R, visible in polar plots depicted in Fig. 7, strongly depends both on the angle  and the angle  .Due to the symmetry of the tensor of the electrical permeability along the x and y axes of the laboratory coordinate system, also the presented distributions of the reflectance (polar plots) are symmetrical with respect to these axes of the coordinate system.

 
and visible in Fig. 9 are symmetrical with respect to the z-axis, i.e. the reflectance strongly depends on the angle  , while it does not depend on the angle  .The obtained curve of the reflectance is characterized by the reflection bandwitch, i.e. the range of the wavelength of the incident light where R = 0.5.The natural non-polarized incident light is the equal mixture of both right-and leftcircularly polarized light.Within the bandwidth, right-circularly polarized light is reflected by a right-handed helix, while left-circularly polarized light is transmitted by the cholesteric liquid crystal, and therefore values of both R and T coefficients are equal 0.5.Concluding, we confirmed the occurrence of the selective reflection of the incident radiation within the bandwidth and transmission of both polarization states outside this bandwidth.  -i.e. the middle of the selective reflection bandwidth presented in Fig. 10).As can be seen, the obtained plot does not depend on the angle  .However, we observed an interesting dependence on the angle  presented in Figure 12(b), i.e. along the radius of the polar plot presented in Fig. 12(a).For small values of angle  , we obtained T = 0.5 (selective reflection of the incident non-polarized light).For larger values of angle  , the reflectance R decreases (due to mutual interferences in CLC both right-and left-circularly polarized light can be transmitted).For large values of angle  (i.e. 90    ), the reflectance reaches 1. Below we also presented interesting dependences of the reflection coefficient R on the angle  of the incident light.
For all cases presented in Figs.13-14, we took 0   , however, for thick layers of the CLC, this angle does not have a significant impact on the obtained dependences.

CONCLUSIONS
The paper attempts to model optical phenomena in different birefringent optical networks sandwiched between two isotropic media.To test the implemented algorithm, we considered the propagation of the light within homogenous isotropic medium and the propagation of the light in homogenous isotropic medium sandwiched within other isotropic medium first.Next, we calculated the reflection coefficient the liquid crystals with different structures, especially cholesteric liquid crystals.Optical phenomena were modeled using the 4×4 matrix method, which takes into account the effect of refraction and multiple reflections between plate interferences, both ordinary and extraordinary waves that occur in birefringent networks.As a result, some interesting reflection and transmission spectra as well as polar plots of the reflection or transmission coefficient have been obtained, reported and discussed.The presented computer simulations allow to understand optical properties of the optical system analysed in this paper for different system parameters.And although the presented results have been obtained for birefringent networks with the positive dielectric anositropy (i.e. ).In addition, assuming different distributions of the optical axis of the investigated liquid crystal, we are able to examine both reflection and transmittance spectra for the deformated LCs caused by the applied external electric field, external mechanical stresses or defects inside the crystal.These issues can be an interesting area for further research.Nevertheless, the results obtained in this paper can be adopted for understanding of optical phenomena in photonic crystals, which have attracted much attention in recent years.This is important from the engineering point of view, because they offer unique optical properties (including ability to control the propagation of light) and can be used as basic elements of all-optical integrated circuits in the future.

Figure 1
Figure1(a) presents the model of the cholesteric liquid crystal embedded in the Cartesian coordinate system xyz and sandwiched between two (i.e.lower and upper) homogenous isotropic media.The entire considered CLC with thickness d is divided into N equal multiple layers parallel to the xy plane.Therefore, each layer can be considered as an anisotropic homegeneous medium.The incident light propagates in the lower medium (n = 0), next partially propagates in the liquid crystal (n = 1,2, ..., N) and the upper isotropic medium (n = N+1).The c-axis (optical axis) of the CLC lies in xy

Figure 1
Figure1(b) presents orientation of the wave vector k of the incident non-polarized monochromatic light with wavelength  , in the lower isotropic medium.The angle  is measured from the z-axis to the direction of the wave wector k , while the angle  is measured from x-axis to the projection of the k director on the xy plane.

Figure 1 :
Figure 1: Model of the investigated optical network sandwiched within two isotropic media (a) and orientation of the wave vector k of the arbitrarily incident light (b).
) is measured from the z-axis to the direction of the vector ) is measured for x-axis to the projection of the () n c director on the xy plane.As a result,


s and p wave vectors of the incoming waves propagating forward, respectively, whereas 2 s and p wave wectors of the reflected waves propagating backward, respectively, in the lower substrate.Moreover, s and p wave vectors transmitted to the upper substrate, respectively, whereas s and p wave vectors reflected in the upper substrate (their amplitudes are equals zero).In an anisotropic medium ( e wave vectors of the transmitted waves, respectively, whereas e wave vectors of the reflected waves, respectively.Since the whole layered birefringent medium is homogeneous in the xy plane,  and  remain the same throughout the whole layered medium.The components of the wave vectors in both x-and y-direction of all elementary waves can be obtained as follow is the wavelength of the monochromatic light in vacuum).The z components of the vectors solved analitycally.The roots can be either real or complex.Because all the coefficients of this equation are real, complex roots are always in conjugate pairs.As a result, there exist four modes of propagation within each -th n homogeneous layer.Two of them (with positive roots all elementary waves in the lower substrate as a function of the Calculation of reflectance and transmittance of optical birefringent networks based on cholesteric liquid crystals Dariusz Grzelczyk et al.Latin American Journal of Solids and Structures, 2019, 16(1 Thematic Section)

A
denote the amplitudes of the electric field of s and p waves that incident the anisotropic medium from the lower substrate, respectively.Moreover, let s B and p B denote the amplitudes of the reflected s and p waves at the interface 0 z  , respectively, whereas s C and p C are the amplitudes of the s and p waves transmitted to the upper substrate at the interface zd  .As it has been already mentioned, the amplitudes of the reflected waves in the upper isotropic medium equal zero.The relation between all mentioned amplitudes can be presented in the vectorcomponents of the matrix 44  M .This relation gives the values of the unknown amplitudes s B , p B of the reflected waves and s C , p C of the transmitted waves, as a function of the known amplitudes s A and p A of the incident waves.Moreover, based on the relation (20), we calculated both the reflection and transmission (20) can be also given in the following form of two matrix equations

E
deal with an incident beam of non-polarized light, therefore the incident light can be written is the intensity of the electric field of the non-polarized light, and simultaneously.Then, the fraction of the energy reflected is given by Calculation of reflectance and transmittance of optical birefringent networks based cholesteric liquid crystals Dariusz Grzelczyk et al.Latin American Journal of Solids and Structures, 2019, 16(1 Thematic Section), e105 10/20 this case, for arbitrary values of the parameters d ,  ,  ,

Figure 2 :
Figure 2: Reflectance spectra for 10 μm d  and different values of the refractive index

Figure 3 :
Figure 3: Reflectance spectra for different values of the thickness d.

Figure 4
Figure 4 presents the reflectance R vs. the angle  for different values of the thickness d and 560 nm   Fig.5also shows that the reflectance R strongly depends on the angle  , while it is not dependent on the angle  (due to isotropy of the considered optical medium).

Figure 4 :
Figure 4: Reflectance vs. the angle  for different values of the thickness d.

Figure 6
Figure 6 presents reflectance spectra of the optical system for different values of the thickness d of the LC cell.The reflection coefficient oscillates locally between the minimum and the maximum, with different amplitude, depending on the wavelength of the incident light.These oscillations are larger for larger values of the thickness d.This regularity is also visible on the polar plots obtained for 560 nm

Figure 6 :
Figure 6: Reflectance spectra for different values of the thickness d.

Figure 7 :
Figure 7: Polar plots of reflectance for different values of d : (a)

Figure 8 :
Figure 8: Reflectance spectra for different values of the thickness d.
Calculation of reflectance and transmittance of optical birefringent networks based on cholesteric liquid crystals Dariusz Grzelczyk et al.Latin American Journal of Solids and Structures, 2019, 16(1 Thematic Section), e105 14/20

Figure 10 :
Figure 10: Reflectance and transmittance spectra of the normally incident non-polarized light.

Figure 11 :
Figure 11: Reflection spectra of the normally incident non-polarized light: (a) for different values of o n ; (b) for different values of e n ; (c) for different values of p ; (d) for different values of d.

Figure 11
Figure 11 presents reflection spectra for different parameters of the CLC.The calculations show that the bandwidth and its position in the wavelength spectrum strongly depends on the refractive indices o n and

Figure 12 :
Figure 12: Polar plot (a) and the plot profile along the radius of the plot (b) the reflectance.
Figure 13 shows the reflectance R vs. the angle  for different values of the wavelength  of the incidend light.As can be seen, does not occur, regardless of the angle  .In turn, the selective reflection for the not normally incident light.In addition, in Figure 14, similar dependences, i.e., the reflection coefficient R vs. the angle  , but for different values of the thickness d of the investigated CLC, are presented.Both in Figs. 13 and 14, strong oscillations of the reflectance R as a function of the angle  are visible.

Figure 15 :
Figure 15: Polar plots of reflectance for 520 nm   and different values of the thickness d :

Figure 16 :
Figure 16: Polar plots of reflectance for 540 nm   and different values of the thickness d :

Figure 17 :
Figure 17: Polar plots of reflectance for 545 nm   and different values of the thickness d :

Figure 18 :
Figure 18: Polar plots of reflectance for 560 nm   and different values of the thickness d : the implemented algorithm can be used for optical systems with negative dielectric anisotropy, i.e. (

Table 1 :
Parameters of the investigated two different isotropic optical media.

Table 2 :
Parameters of the investigated planar structure of LC and isotropic optical medium.

Table 3 :
Parameters of the investigated homeotropic structure of LC and isotropic optical medium.

Table 4 :
Parameters of the investigated helical structure of LC and isotropic optical medium.