Open access peer-reviewed chapter

A Theoretical Model of the Holographic Formation of Controllable Waveguide Channels System in Photopolymer Liquid Crystalline Composition

By Artem Semkin and Sergey Sharangovich

Submitted: October 10th 2017Reviewed: February 4th 2018Published: August 1st 2018

DOI: 10.5772/intechopen.74838

Downloaded: 190


Rapid development of the integrated optics and photonics makes it necessary to create cheap and simple technology of optical waveguide systems formation. Photolithography methods, widely used for these tasks recently, require the production of a number of precision amplitude and phase masks. This fact makes this technology expensive and the formation process long. On another side there is a cheap and one-step holographic recording method in photopolymer compositions. Parameters of the waveguide system formed by this method are determined by recording geometry and material’s properties. Besides, compositions may contain liquid crystals that make it possible to create elements, controllable by external electric field. In this chapter, the theoretical model of the holographic formation of controllable waveguide channels system in photopolymer liquid crystalline composition is developed. Special attention is paid to localization of waveguides in the media caused by light field attenuation during the formation process.


  • photopolymer
  • liquid crystal
  • waveguide
  • holography

1. Introduction

The ability to form waveguide systems for optical and terahertz radiation in photopolymerizable compositions recently is of great interest among researchers: [1, 2, 3, 4]. Formed holographically or by photolithography methods, such waveguides are widely used in the integrated optics and photonics devices. Besides, it seems urgent to create the manageable light guides, in which the light propagation conditions can be controlled by external influences, such as an electric field. One of the possible solutions of this problem is a holographic recording of the waveguide channels in a photopolymer composition containing liquid crystals.

The aim of this chapter is to develop the theoretical model of holographic formation of controllable waveguide channels system in photopolymer liquid crystalline composition with dye sensitizer, also known as polymer-dispersed (PDLC) or polymer-stabilized (PSLC) liquid crystals.

2. Theoretical model

We consider the incidence of two plane monochromatic waves E0and E1with incidence angles θ0and θ1on the PDLC (PSLC) sample for two formation geometries: transmission (Figure 1a) and reflection (Figure 1b) ones.

Figure 1.

Waveguide channels holographic formation: (a) transmission recording geometry and (b) reflection recording geometry.

Thus, recording waves in the general case can be described as:


where e0,e1are the unit polarizations vectors of the beams; E0r,E1rare the spatial amplitude distributions; φ0r,φ1rare spatial phase distributions; k0,k1are the wave vectors; N0,N1are wave normals; αrtis the absorption coefficient; ω=2π/λ, λis the wavelength of the recording radiation.

Investigated material is characterized by optical anisotropy properties; thus, in the material, each of recording wave (Eq. (1)) will be divided into two mutually orthogonal ones called ordinary and extraordinary. So, in the sample, Eq. (1) should be rewritten as:


where m = o corresponds to ordinary waves and m = e corresponds to extraordinary waves, respectively.

For two recording geometries (see Figure 1), the spatial distributions of the forming field intensities are determined by the following expression:


where ITrt, IRrtare the spatial distributions of light field intensity for transmission and reflection geometries, respectively; I0r=E02r+E12r, mmr=2E0rE1rE02r+E12re0me1mare the local contrasts of interference patterns; Km=k0mk1m; dis the thickness of the material.

Under the influence of light field in photopolymer liquid crystalline composition with dye sensitizer, the dye molecule absorbs a light radiation quantum with the dye radical and the primary radical initiator formation. The radical of dye is not involved in further chemical reactions and turns to colorless leuco form.

Thus, during the waveguide channel system formation, dye concentration decreases. This fact causes the light-induced decreasing of light absorption. The absorption coefficient dependence of a mount of absorbed radiation can be written as [5]:


where αrt=α0Kdrtis the absorption coefficient with light-induced change taken into account; Kdrtis dye concentration; α0is the absorption of one molecule; βqis quantum yield of dye; and Iabsrtis the light intensity from Bouguer-Lambert law:


where IabsTrt, IabsRrtare defined for transmission and reflection geometries, respectively.

Then, the solution of Eq. (3) for light-induced absorption change for transmission and reflection geometries is obtained in [5] as follows:


where αsubis the substrate absorption coefficient and Kd0is the initial dye concentration.

The process of waveguide channels’ holographic formation is described by the kinetic equations system (KES), written for monomer concentration and refraction index [6, 7, 8]:


where Mnis the initial concentration of the monomer; Kg, Kbare parameters of the rate of growth and breakage of the polymer chain, respectively; βis the parameter of initiation reaction; τ0is lifetime of the excited state of the dye molecule; DMmrt, DLCmrtare the diffusion coefficients of the monomer and liquid crystal, respectively; δnp, δnlcare weight coefficients of the contribution of photopolymerization and diffusion processes; and Imrtis the intensity distribution (Eq. (3)).

Diffusion coefficients can be defined from the following equations:


where DMnand DLCnare the initial diffusion coefficients, respectively; sM, sLCare rates of reduction in time; Ln, Lmrtare the initial and current concentrations of liquid crystal.

Weight coefficients δnpand δnlcfrom Eq. (7) are found from the Lorentz-Lorentz formula [8]:


where αM, αP, αLCare the polarizability of monomer, polymer, and liquid crystal molecules, respectively; ρM, ρLCare the density of the monomer and liquid crystal, respectively; WM, WLCare molecular weights; lis the average length of polymeric chains; and nstis the refractive index of the composition prior to the start of the recording process, determined by the Lorentz-Lorentz formula from the refractive indices of the monomer and the liquid crystal [8]:


where nM, nLCare the monomer and liquid crystal refractive indices.

Because of periodical character of the forming fields’ intensities’ spatial distributions, solution of KES can be found as a sum of H spatial harmonics [8]:


where Mjmrτ=12πππMjmrτcosjKmrdKmr, njmrτ=12πππnjmrτcosjKmrdKmrare the monomer concentration and refractive index harmonics amplitudes, respectively; nstis the PDLC (PSLC) refractive index at τ=0; τ=t/TMmis the relative time; TMm=1DMnmKm2is the components’ diffusion characteristic time.

By substituting Eq. (13) to the KES (Eqs. (7 and 8)) and using the orthogonality of spatial harmonics, a system of coupled kinetic differential equations for the amplitudes of monomer concentration harmonics can be obtained for DMmrt=DLCmrt=DMnm(stable diffusion coefficients) and Kdrt=Kd0(stable absorption) [6]:


and also a system of differential equations for the amplitudes of refraction index harmonics:


In equation systems Eqs. (14) and (15), a coefficient matrix is introduced [6]:


where e1m=2bsm1+Lsm; e11m=2bsm1+3Lsm/2; e2m=2msm4bsm; e3m=2bsmLsm/2; msm=mmr; Lsm=Lmr=mmr2/16; bsm=bmr=TpmrTMmrare the parameters that characterize the ratio of polymerization and diffusion rates; and Tpmris the characteristic polymerization time:


TMmris the characteristic diffusion time:


Coefficients aj,lmrdescribe the contributions of photopolymerization and diffusion recording mechanisms. However, for analysis of equation systems Eqs. (14) and (15), it is convenient to introduce the coupling coefficients cj,lmr=aj,lmrj2δj,l(δj,lis the Kronecker symbol), which characterizes the coupling between j and l harmonics. The difference between coefficients cj,lmrand aj,lmrcharacterizes the contribution of monomer diffusion to the recording process and it is proportional to the second degree of the harmonic’s number. The increase of this contribution according to the harmonic’s number is due to grating period decrease and, respectively, the diffusion characteristic time decrease for this harmonic.

For solution of the coupled differential equations system (14), the initial conditions should be introduced:


The solution can be found using the operator method [6]. The general solution for the spatial amplitude profiles of monomer concentration harmonics will be:


where functional dependencies of coefficients λlmrare defined as the roots of the characteristic equation cj,lmrλlmr=0. Analysis shows that λlmris real, different, and negative. Coefficients Aj,lmrare defined as solutions of a system of linear algebraic equations:


where λlm=λlmr, Aj,lm=Aj,lmr, cj,lm=cj,lmr.

Then, by substituting (20) to (15) and by integrating the resulting differential equations with initial conditions, we get.


The general solution for the amplitude of j-harmonic of the refractive index can be found:


where npjmrτ=δnpl=0Haj,lmrq=0HAl,qmr1expλqmrτλqmr, nlcjmrτ=δnlcj2q=0HAj,qmr1expλqmrτλqmr, j = 0,…,H.

Thus, Eqs. (13), (20), and (23) are the general solutions of nonlinear photopolymerization diffusion holographic recording of waveguide channels system in PDLCs (PSLCs) in the case of stable diffusion coefficients and stable absorption. They define kinetics of spatial profiles of monomer concentration—Mmrτand refractive index nmrτ.

In case of high nonlinearity, it is possible to form the specific spatial profile of refractive index. The nonlinearity of recording is achieved by changing the ratio of the photopolymerization and diffusion mechanism contributions to the process of structure?s formation.

3. Numerical simulations

To investigate the formation processes, numerical simulations of first four refractive index harmonics kinetics for transmission and reflection geometries were made with the following parameters: λ=633nm and θ0=θ1=30°for transmission geometry (Figure 1a) and θ0=θ1=60°for reflection geometry (Figure 1b); E0r=E1r=1; φ0r,=φ1r=0; e0,e1are oriented with respect to extraordinary waves in material; d=10μm; δnp/δni=0.5; and for two values of parameter b=be=0.2and b=be=5in the absence of absorption. Values b=0.2and 5are chosen as an example to investigate two common cases: b<<1and b>>1. Simulations were made by Eq. (23); results are shown in Figure 2.

Figure 2.

Refractive index harmonics kinetics for b = 0.2 (a) and b = 5 (b).

As can be seen from Figure 2, in the case of predomination of polymerization (be=0.2, Figure 2a), the structure forms quickly, but amplitudes of higher harmonics are high, so the spatial profile of refractive index changes has an inharmonic character. In another case of diffusion predominance (be=5, Figure 2b), spatial profile is quasi-sinusoidal, but it is slower. These effects can be explained by the following. When polymerization is rapid, in the area of the maximum of intensity distribution (Eq. (3)), molecules of liquid crystal do not have time to diffuse, so the concentrations’ gradient is not high enough; monomer molecules do not diffuse from minimums of intensity distribution. Thus, profile of refractive index distribution becomes inharmonic. So, it can be concluded that ratio of polymerization and diffusion rates bmrτdefines the distribution nmrτ.

Corresponding distributions nerτ=50for two examined cases in the absence of absorption are shown in Figure 3.

Figure 3.

Refractive index distributions for b = 0.2 (a) and b = 5 (b).

In Figure 3, a new spatial coordinate is introduced: ξ=yfor transmission geometry (see Figure 1a) and ξ=zfor reflection geometry (see Figure 1b).

According to Eqs. (17) and (18), parameter bmrτdepends on formation field’s amplitude and phase distributions as well as material properties. So, by controlling these parameters one can create any distribution of refractive index nmrτ. As is shown in Figure 3, distributions can be quasi-rectangular (Figure 3a) and quasi-sinusoidal (Figure 3b), so they can be mentioned as waveguides systems.

It should be also noted that the amplitude of refractive index change when b=0.2(Figure 3a), due to high amplitudes of higher harmonics (see Figure 2a), is lower than in the case of b=5.

To investigate the impact of absorption on the spatial profile of structure along the thickness of the sample, numerical simulations were made with the following parameters: λ=633nm; θ0=θ1=30°for transmission geometry (Figure 1a) and θ0=θ1=60°for reflection geometry (Figure 1b); E0r=E1r=1; φ0r,=φ1r=0; e0,e1are oriented with respect to extraordinary waves in material; d=10μm; δnp/δni=0.5; b=be=5; αd=2Np. Simulations were made by Eq. (23) for the amplitude of the first harmonic of refractive index and for two formation geometries. Results are shown in Figure 4.

Figure 4.

Impact of light absorption: transmission geometry (a) and reflection geometry (b).

As shown in Figure 4, absorption of the formation field causes dependence of the first harmonic on z-coordinate (see Figure 1) during the recording process. Thus, the spatial profile of refractive index changes its character in time and space.

In Figure 4, some characteristic cases of localization of maximum of refractive index change in the thickness can be seen. In particular, in the case of transmission geometry (Figure 4a) at τ=2.5, maximum is localized between 1.5 and 3.5 μm of the material thickness. For reflection geometry (Figure 4b), one maximum is localized near 5 μm at τ=1.25, and there are two other maximums at τ>2.5. These effects take place because of dependence bmr. Absorption causes the dependence on z-coordinate and spatial distributions of amplitude, and the phase of recording field (Eqs. (1) and (2)) causes the dependence on x and y coordinates. In more general cases, there is also a temporal dependence bmrτcaused by photo-induced effects (Eq. (6)).

To illustrate the localization of waveguide channels into the material, in Figure 5, there are some more two-dimensional spatial distributions of refractive index, obtained by Eq. (23) for transmission and reflection geometries—b=be=5and for different values of τ.

Figure 5.

Two-dimensional refractive index distributions: transmission geometry (a–c) and reflection geometry (d–f).

It follows from the abovementioned that by controlling the spatial distribution of light field and taking the absorption effects into account, the waveguide systems localized into the material can be holographically formed. Localization is important because of the necessity to create a predetermined refractive index profile independent from properties of possible substrates (electrodes, glass, polymer layers, etc.). In particular, to create the waveguide system localized near the bottom substrate, formation geometry should be transmissive and formation should be stopped before τ=0.5, for localization into the center of the sample— 0.5<τ<2.5—for transmission and reflection geometries, and for τ>2.5, waveguides will be localized near the top substrate for transmission geometry and not localized for reflection geometry. The given values of τare valid for the formation parameters given earlier.

In our recent works [9, 10, 11, 12, 13, 14, 15, 16, 17, 18], we have seen that the impact of the external electric field on the holographic structure formed in PDLC (PSLC) leads to the refractive index change due to electro-optical orientation mechanisms, that is typical to liquid crystals, and the recorded structure can be “erased”. So, it can be supposed that the impact of external electric field on individual waveguides will lead to its “erasing”, and it is possible to “switch off” some waveguides or change the system’s period, and so on.

4. Conclusion

Thus, in this chapter, the theoretical model of holographic formation of controllable waveguide channels system in photopolymer liquid crystalline composition is developed. The most general cases are described by the developed model, and numerical simulations were made for plane recording waves and stable absorption cases. Special attention is paid to localization of waveguides in the media caused by light field attenuation during the formation process.

It is shown that parameters of the waveguide system formed by holography’s methods are determined by recording geometry and material’s properties. Also, by control of these parameters, waveguides can be localized into the sample that makes them independent from substrates. Also, introduced compositions contain liquid crystals that make it possible to create elements, controllable by external electric field.

Obtained results can be used for new photonics devices based on photopolymer liquid crystalline composition development.


The work is performed as a project that is part of Government Task of Russian Ministry of Education (project No. 3.1110.2017/4.6).

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Artem Semkin and Sergey Sharangovich (August 1st 2018). A Theoretical Model of the Holographic Formation of Controllable Waveguide Channels System in Photopolymer Liquid Crystalline Composition, Emerging Waveguide Technology, Kok Yeow You, IntechOpen, DOI: 10.5772/intechopen.74838. Available from:

chapter statistics

190total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Application of Numeric Routine for Simulating Transients in Power Line Communication (PLC) Systems

By Afonso José do Prado, Luis Henrique Jus, Melissa de Oliveira Santos, Elmer Mateus Gennaro, André Alves Ferreira, Thainá Guimarães Pereira, Aghatta Cioqueta Moreira, Juliana Semiramis Menzinger, Caio Vinícius Colozzo Grilo, Marinez Cargnin Stieler and José Pissolato Filho

Related Book

Emerging Microwave Technologies in Industrial, Agricultural, Medical and Food Processing

Edited by Kok Yeow You

First chapter

Introductory Chapter: RF/Microwave Applications

By Kok Yeow You

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us