# Waveguide Propagation of Light in Polymer Porous Films Filled with Nematic Liquid Crystals

Authors: A. D. Kiselev; S. V. Pasechnik; D. V. Shmeliova; A. P. Chopik; D. A. Semerenko; A. V. Dubtsov
Journal: Advances in Condensed Matter Physics
 (2019)

---

## Abstract

We theoretically analyze the waveguide regime of light propagation in a cylindrical pore of a polymer matrix filled with liquid crystals assuming that the effective radial optical anisotropy is biaxial. From numerical analysis of the dispersion relations, the waveguide modes are found to be sensitive to the field-induced changes of the anisotropy. The electrooptic properties of the polymer porous polyethylene terephthalate (PET) films filled with the nematic liquid crystal 5CB are studied experimentally and the experimental results are compared with the results of the theoretical investigation.


---

## Body

## 1. Introduction

Liquid crystals (LCs) are known to be of primary importance in the modern display industry. A unique combination of physical properties of LCs provides a means to control light propagation through micron-sized LC layers using electric fields at low operating voltages and small energy consumption. The key feature of liquid crystal materials that underlies excellent electrooptical characteristics of LCs and makes them promising for photonic applications is the high sensitivity of LC orientational structures to both external fields and conditions at bounding surfaces. In particular, the use of LC in microstructured waveguides [1, 2] and waveguides based on photonic crystals [3, 4] is mainly determined by this feature providing an efficient tool to adjust optical characteristics of waveguide photonic devices by means of the external (electric) field.Utilization of LCs in nondisplay applications such as photonic devices [5] and a variety of sensors [6, 7] faces a number of challenges. These include a slow response time of orientational transformations for the devices of THz photonics based on thick layers of LC [8], high operating voltages needed to tune propagation of light in photonics fibers filled with LC [9], and relatively low sensitivity of planar LC layers to impurities in biosensensing applications [7]. Some of these problems can be solved by usage of composite materials, like PDLC, LC emulsions, and porous films filled with LC.In this paper the electrooptical properties of the porous polymer films filled with nematic liquid crystals will be our primary concern. Recently such type of composite material—the polymer polyethylene terephthalate (PET) film—with normally oriented open pores of submicron and micron sizes, filled with LC, was experimentally studied in [10, 11]. It was shown that such systems are very promising for applications in photonic waveguide and terahertz (THz) devices [11, 12].In particular, we have observed pronounced changes of light intensity transmitted through the film induced by applying low frequency ac electric field. Similarly, modulation of the light intensity can also be caused by heating of the sample via absorption of blue light by the azo-dye layer adsorbed on the internal surfaces of the pores. The hypothesis put forward in [11] to explain the experimental results assumes that breaking of the waveguide regime of light propagation is responsible for the light modulation.These results motivate theoretical considerations of this paper that aimed to examine the waveguide regime of light propagation in cylindrical LC pores. In our calculations, we shall treat the LC material as an effective medium with variable anisotropy embedded in an isotropic environment. Though similar approach was previously used in [2], it was not applied to analyze the waveguide regime. Rigorous analysis of the problem based on numerical solution of Maxwell’s equations [13] showed that waveguide propagation of light in a cylindrical cavity surrounded by isotropic media is crucially influenced by the orientational structure formed inside the pore [14–16]. On the other hand, the type of the orientational configuration is governed by a number of parameters such as the Frank elastic constants, the surface anchoring strengths and orientation of the easy axes, and the diameter of LC pores [17–21]. The paper is organized as follows.In Section2 we present the details of our theoretical analysis and systematically study the effects of effective LC anisotropy on the waveguide regime of light propagation. In Section 3, we describe the experimental results obtained for the porous PET films with micron and submicron pores filled with LC 5CB. In our discussion of these results we, in particular, demonstrate a qualitative agreement between the experimental data and the theoretical predictions. Concluding remarks are given in Section 4.


## 2. Theory

### 2.1. Basic Equations

We begin with the Maxwell equations for a monochromatic electromagnetic wave(1)∇×E=ik0μH,∇×H=-ik0D,D=ϵE,where k0=ω/c is the free space wavenumber and ϵ is the dielectric tensor, which propagate along the axis of a cylindrical waveguide (the z axis). Assuming that the electromagnetic field can be written in the factorized form {E,H}ei(kzz-ωt), we shall use the cylindrical coordinate system {ρ,ϕ,z} and separate out the longitudinal and transverse parts of both the electromagnetic field and the gradient operator as follows:(2)E=ET+Ezz^,∇=∇T+z^∂z=∇T+ikzz^,(3)∇T=x^∂x+y^∂y=e^ρ∂ρ+e^ϕρ-1∂ϕ,where e^ρ=cos⁡ϕx^+sin⁡ϕy^ and e^ϕ=-sin⁡ϕx^+cos⁡ϕy^. With the help of the relation z^×(∇×E)=∇TEz-ikzET, we obtain the transverse part of Maxwell’s equations (1) in the following form:(4)∇TEz-ikzET=ik0μz^×HT,z^×∇THz-ikzz^×HT=ik0DT.In what follows, we, following the results of the paper [2], restrict ourselves to the case of effective radial biaxial anisotropy characterized by the constitutive relations of the following form:(5)DT=ϵTET=ϵρ00ϵϕEρEϕ,Dz=ϵzzEz.Note that the case of radial anisotropy was also considered in the theory of light scattering by optically anisotropic cylindrical particles developed in [22, 23].From equations (4) it is not difficult to express the transverse components of the electromagnetic field {ET,HT}, in terms of the longitudinal components {Ez,Hz}, as follows:(6)ikz2-k02μϵTET=kz∇TEz-k0μz^×∇THz,(7)ikz2-k02μϵTz^×HT=kzz^×∇THz-k0ϵT∇TEz.For convenience, we shall introduce the dimensionless parameters(8)r=k0ρ,qz=kzk0,qρ2=μϵρ-qz2=nρ2-qz2,qϕ2=μϵϕ-qz2=nϕ2-qz2and write down the expressions for the transverse components in the explicit form:(9)Eρ=iqρ2rqzr∂rEz+μ∂ϕHz,(10)Eϕ=iqϕ2rqz∂ϕEz-μr∂rHz,(11)Hρ=iqϕ2rqzr∂rHz-ϵϕ∂ϕEz,(12)Hϕ=iqρ2rqz∂ϕHz+ϵρr∂rEz.Equations for the longitudinal components can be obtained using the divergence relations:(∇·D)=(∇·H)=0 (these follow from Maxwell’s equations (1)), which can be conveniently recast into the following form:(13)∇T·DT=-ikzDz,∇T·HT=z^×∇T·z^×HT=-ikzHz.After eliminating the transverse components from the relations (6)–(7) with the help of Eq. (13), we have(14)∇TϵTkz2-k02μϵT-1∇T-ϵzzEz=qz∇Tkz2-k02μϵT-1z^×∇THz(15)z^×∇Tμkz2-k02μϵT-1z^×∇T-μHz=qzz^×∇Tkz2-k02μϵT-1∇TEz,where the operators are given by(16)-r2∇TϵTkz2-k02μϵT-1∇T=ϵρqρ2r∂r2+ϵϕqϕ2∂ϕ2,(17)-r2z^×∇Tkz2-k02μϵT-1z^×∇T=1qϕ2r∂r2+1qρ2∂ϕ2,(18)-r2∇Tkz2-k02μϵT-1z^×∇T=r2z^×∇Tkz2-k02μϵT-1∇T=1qϕ2-1qρ2r∂r∂ϕ=μϵρ-ϵϕqρ2qϕ2r∂r∂ϕ.For the Fourier harmonics of the longitudinal components(19)Ez=Emreimϕ,Hz=-iHmreimϕ,∂ϕ→im,where m∈Z is the azimuthal number enumerating the harmonics, we derive the following system of equations:(20)ϵρqϕ2r∂r2-ϵϕqρ2m2+ϵzzq2r2mκzr∂rmκzr∂rμqρ2r∂r2-μqϕ2m2+μq2r2EmHm=0,where q2=qρ2qϕ2, κz=qzμΔϵ and Δϵ=ϵρ-ϵϕ.


### 2.2. Boundary Conditions and Dispersion Relations

For the special case of a zero mode withm=0, the solution of the system (20) can be readily found and is given by(21)Ez=E0r=A0J0rzρqρr,iHz=H0r=B0J0qϕr,rzρ2=ϵzzϵρ.In general case, we introduce the matrix notations(22)Aα=ϵρqϕ2α2-ϵϕqρ2m2mκzαmκzαμqρ2α2-μqϕ2m2,B=q2ϵzz00μand will search for the solutions of the system (20) that are regular at r=0 (on the waveguide axis) in the form of a power series:(23)Em,αrHm,αr=rα∑n=0∞cnrn.Substitution of the expansion (23) into Eq. (20) gives the equation for the coefficient c0 and the recurrence relations for cn:(24)A0c0=0,cn=-An-1Bcn-2,An≡Aα+n.The values of the parameter α can be found from the singularity condition for the matrix A0:(25)det⁡Aα=nρ2q2α4-α2m2nϕ2qρ4+nρ2qϕ4+κz2+nϕ2q2m4=0,giving two positive values of the parameter α∈{α+,α-} that define two solutions which are regular at r=0. Note that, similar to the case of Bessel functions of the first kind Jν(x), we can use the expression for the determinant of the matrix An (at α=α±):(26)det⁡An=nρ2q2α+n4-α4-mκz2α+n2-α2to prove convergence of the power series (23), representing an entire function.The obtained solution can be used not only for analyzing the waveguide regime, but also for describing the light scattering by radially anisotropic scatterers without recourse to the simplifying assumptions such asEz=0 used in the literature [22, 23].Now, following [2], we concentrate on the important special case with ϵρ=ϵϕ≡ϵt. In this case we have(27)Ez=Emr=AmJmrztqtr,iHz=Hmr=BmJmqtr,rzt2=ϵzzϵt,where qt2=nt2-qz2.Outside the waveguide atr>k0R≡x, the solution that exponentially decreases at infinity, Ez,Hz∝e-r/r, is given by(28)Ez=Emr=amKmqcr,iHz=Hmr=bmKmqcr,qc2=qz2-nc2where Km(x) is the modified Bessel function [24] and nc=μcϵc is the refractive index of the ambient medium.From the continuity conditions for the tangential components of the electromagnetic field at the waveguide boundary(29)Ezx+0=Ezx-0,Hzx+0=Hzx-0,(30)Eϕx+0=Eϕx-0,Hϕx+0=Hϕx-0,x≡k0Rwe have(31)Cmαmβm=0,αm=AmJmrztxt=amKmxc,βm=BmJmxt=bmKmxc,where xt=qtx, xc=qcx and the matrix Cm is given by(32)Cm=qzmnt2-nc2μqc2RJmxt+μcqt2RKmxcϵtqc2RJmrztxt+ϵcqt2RKmxcqzmnt2-nc2,(33)RJmx=xJm′xJmx=m-xJm+1xJmx,RKmx=xKm′xKmx=m-xKm+1xKmx.Dispersion relations for the harmonics of the electromagnetic field, linking kz and k0=ω/c, can be deduced from the singularity condition for the matrix (32):(34)det⁡Cmqz,x=0,nc2≤qz2=kzk02≤nt2,x=k0R.


### 2.3. Numerical Results

We begin our analysis with the case of waveguide modes with the zero azimuthal number:m=0. In this case, the dispersion relations are simplified and reduced to the relations for two types of the waves: TE waves with Ez=0(35)μqc2RJ0xt=-μcqt2RK0xc⇒μJ1xtxtJ0xt=-μcK1xcxcK0xc and TH waves with Hz=0(36)ϵtqc2RJ0rztxt=-ϵcqt2RK0xc⇒ϵzJ1rztxtrztxtJ0rztxt=-ϵcK1xcxcK0xc.The roots of Eqs. (35) and (36) are located between the corresponding zeros of the Bessel function J0 and will be numbered in ascending order by the index l. In our subsequent calculations, we shall use the estimates taken from the paper [11]: the pores of radius R that ranged from 200 nm to 600 nm are filled with the nematic liquid crystal 5CB characterized by the refractive indices ne≈1.72 and no≈1.5 (μ=1), the effective refractive index of the polymer matrix is nc≈1.57 (μc=1), and the wavelength of the laser radiation is λ=632 nm. Clearly, in this case, we have 2≈1.98≤k0R≤5.96≈6.Figure1 shows how liquid crystal anisotropy (the value of the longitudinal refractive index nz) affects the dispersion curves for the waveguide modes TE0l and TH0l. It can be seen that the curves for the TE modes, which are computed for the planar director configuration with the effective refractive index nt=ne=1.72, are independent of nz (see Eq. (35)), whereas the difference between the dispersion curves for the TH and TE modes becomes more pronounced as the difference between the refractive indices nt and nz increases. The radial distributions of nonvanishing components of the electromagnetic fields for the waveguide modes TH01 and TE01 are presented in Figures 2 and 3, respectively.Figure 1
Dependence of the dimensionless parameterk0R≡ωR/c on qz=kz/k0 for TEml (Ez=0) and THml (Hz=0) waves with m=0: (a) l=1 and (b) l=2 (index l enumerates the branches of roots in ascending order) at different values of the refractive index nz.
(a)
m = 0 and l=1 (b)
m = 0 and l=2Figure 2
Electromagnetic field versus the dimensionless radial parameterr=k0ρ for the mode TH01 (Hz=0) at (a) nz=1.4 and (b) nz=1.5.
(a)
m = 0 and nz=1.4 (b)
m = 0 and nz=1.5Figure 3
Electromagnetic field versus the dimensionless radial parameterr=k0ρ for the mode TE01 (Ez=0) at (a) nz=1.4 and (b) nz=1.5.
(a)
m = 0 and nz=1.4 (b)
m = 0 and nz=1.5The results illustrating effect ofnz on the hybrid modes (Ez≠0 and Hz≠0) with m=1 and m=2 are plotted in Figure 4. It is seen that as compared to the TH modes the curves for the hybrid modes are much less sensitive to variations in nz. Interestingly, for such modes, similar to the case of TH01 and TE01 modes, the waveguide regime may occur in the porous films, where 2≤k0R≤6.Figure 4
Dependence of the dimensionless parameterk0R≡ωR/c on qz=kz/k0 for the waves with (a) m=1 and (b) m=2 at different values of the refractive index nz.
(a)
m = 1 (b)
m = 2When the electric field is applied, the effective refractive indices of the LC (in a typical situationnt decreases, whereas nz increases) and the waveguide regime in the porous film can be suppressed. This is obvious if the refractive index nt, which approaches no<nc as the field increases, becomes less than nc, so that the condition of total internal reflection appears to be violated. Figure 5 shows the dispersion curves for the case when the configuration distorted by the field is close to axial with nz=ne, but nt is still higher than nc. It is interesting that under these conditions, for a film with 2≤k0R≤6, the conditions for the waveguide propagation regime are satisfied only for the hybrid mode with m=1.Figure 5
Dimensionless parameterk0R≡ωR/c as a function of qz=kz/k0 for the waves with 0≤m≤2 at nz=1.72>nt=1.6≈nc.Radial distributions for the components of the electric and magnetic fields are shown in Figures6 and 7, respectively. These figures, in particular, demonstrate that the degree of the field localization inside the pore decreases as the director configuration approaches the axial structure.Figure 6
Components of electromagnetic field as a function of the dimensionless radial parameterr=k0ρ for the mode with m=1 at (a) nz=1.5 and (b) nz=1.72.
(a)
m = 1 and nz=1.5 (b)
m = 1 and nz=1.72Figure 7
Components of electromagnetic field as a function of the dimensionless radial parameterr=k0ρ for the mode with m=1 at (a) nz=1.5 and (b) nz=1.72.
(a)
m = 1 and nz=1.5 (b)
m = 1 and nz=1.72


## 2.1. Basic Equations

We begin with the Maxwell equations for a monochromatic electromagnetic wave(1)∇×E=ik0μH,∇×H=-ik0D,D=ϵE,where k0=ω/c is the free space wavenumber and ϵ is the dielectric tensor, which propagate along the axis of a cylindrical waveguide (the z axis). Assuming that the electromagnetic field can be written in the factorized form {E,H}ei(kzz-ωt), we shall use the cylindrical coordinate system {ρ,ϕ,z} and separate out the longitudinal and transverse parts of both the electromagnetic field and the gradient operator as follows:(2)E=ET+Ezz^,∇=∇T+z^∂z=∇T+ikzz^,(3)∇T=x^∂x+y^∂y=e^ρ∂ρ+e^ϕρ-1∂ϕ,where e^ρ=cos⁡ϕx^+sin⁡ϕy^ and e^ϕ=-sin⁡ϕx^+cos⁡ϕy^. With the help of the relation z^×(∇×E)=∇TEz-ikzET, we obtain the transverse part of Maxwell’s equations (1) in the following form:(4)∇TEz-ikzET=ik0μz^×HT,z^×∇THz-ikzz^×HT=ik0DT.In what follows, we, following the results of the paper [2], restrict ourselves to the case of effective radial biaxial anisotropy characterized by the constitutive relations of the following form:(5)DT=ϵTET=ϵρ00ϵϕEρEϕ,Dz=ϵzzEz.Note that the case of radial anisotropy was also considered in the theory of light scattering by optically anisotropic cylindrical particles developed in [22, 23].From equations (4) it is not difficult to express the transverse components of the electromagnetic field {ET,HT}, in terms of the longitudinal components {Ez,Hz}, as follows:(6)ikz2-k02μϵTET=kz∇TEz-k0μz^×∇THz,(7)ikz2-k02μϵTz^×HT=kzz^×∇THz-k0ϵT∇TEz.For convenience, we shall introduce the dimensionless parameters(8)r=k0ρ,qz=kzk0,qρ2=μϵρ-qz2=nρ2-qz2,qϕ2=μϵϕ-qz2=nϕ2-qz2and write down the expressions for the transverse components in the explicit form:(9)Eρ=iqρ2rqzr∂rEz+μ∂ϕHz,(10)Eϕ=iqϕ2rqz∂ϕEz-μr∂rHz,(11)Hρ=iqϕ2rqzr∂rHz-ϵϕ∂ϕEz,(12)Hϕ=iqρ2rqz∂ϕHz+ϵρr∂rEz.Equations for the longitudinal components can be obtained using the divergence relations:(∇·D)=(∇·H)=0 (these follow from Maxwell’s equations (1)), which can be conveniently recast into the following form:(13)∇T·DT=-ikzDz,∇T·HT=z^×∇T·z^×HT=-ikzHz.After eliminating the transverse components from the relations (6)–(7) with the help of Eq. (13), we have(14)∇TϵTkz2-k02μϵT-1∇T-ϵzzEz=qz∇Tkz2-k02μϵT-1z^×∇THz(15)z^×∇Tμkz2-k02μϵT-1z^×∇T-μHz=qzz^×∇Tkz2-k02μϵT-1∇TEz,where the operators are given by(16)-r2∇TϵTkz2-k02μϵT-1∇T=ϵρqρ2r∂r2+ϵϕqϕ2∂ϕ2,(17)-r2z^×∇Tkz2-k02μϵT-1z^×∇T=1qϕ2r∂r2+1qρ2∂ϕ2,(18)-r2∇Tkz2-k02μϵT-1z^×∇T=r2z^×∇Tkz2-k02μϵT-1∇T=1qϕ2-1qρ2r∂r∂ϕ=μϵρ-ϵϕqρ2qϕ2r∂r∂ϕ.For the Fourier harmonics of the longitudinal components(19)Ez=Emreimϕ,Hz=-iHmreimϕ,∂ϕ→im,where m∈Z is the azimuthal number enumerating the harmonics, we derive the following system of equations:(20)ϵρqϕ2r∂r2-ϵϕqρ2m2+ϵzzq2r2mκzr∂rmκzr∂rμqρ2r∂r2-μqϕ2m2+μq2r2EmHm=0,where q2=qρ2qϕ2, κz=qzμΔϵ and Δϵ=ϵρ-ϵϕ.


## 2.2. Boundary Conditions and Dispersion Relations

For the special case of a zero mode withm=0, the solution of the system (20) can be readily found and is given by(21)Ez=E0r=A0J0rzρqρr,iHz=H0r=B0J0qϕr,rzρ2=ϵzzϵρ.In general case, we introduce the matrix notations(22)Aα=ϵρqϕ2α2-ϵϕqρ2m2mκzαmκzαμqρ2α2-μqϕ2m2,B=q2ϵzz00μand will search for the solutions of the system (20) that are regular at r=0 (on the waveguide axis) in the form of a power series:(23)Em,αrHm,αr=rα∑n=0∞cnrn.Substitution of the expansion (23) into Eq. (20) gives the equation for the coefficient c0 and the recurrence relations for cn:(24)A0c0=0,cn=-An-1Bcn-2,An≡Aα+n.The values of the parameter α can be found from the singularity condition for the matrix A0:(25)det⁡Aα=nρ2q2α4-α2m2nϕ2qρ4+nρ2qϕ4+κz2+nϕ2q2m4=0,giving two positive values of the parameter α∈{α+,α-} that define two solutions which are regular at r=0. Note that, similar to the case of Bessel functions of the first kind Jν(x), we can use the expression for the determinant of the matrix An (at α=α±):(26)det⁡An=nρ2q2α+n4-α4-mκz2α+n2-α2to prove convergence of the power series (23), representing an entire function.The obtained solution can be used not only for analyzing the waveguide regime, but also for describing the light scattering by radially anisotropic scatterers without recourse to the simplifying assumptions such asEz=0 used in the literature [22, 23].Now, following [2], we concentrate on the important special case with ϵρ=ϵϕ≡ϵt. In this case we have(27)Ez=Emr=AmJmrztqtr,iHz=Hmr=BmJmqtr,rzt2=ϵzzϵt,where qt2=nt2-qz2.Outside the waveguide atr>k0R≡x, the solution that exponentially decreases at infinity, Ez,Hz∝e-r/r, is given by(28)Ez=Emr=amKmqcr,iHz=Hmr=bmKmqcr,qc2=qz2-nc2where Km(x) is the modified Bessel function [24] and nc=μcϵc is the refractive index of the ambient medium.From the continuity conditions for the tangential components of the electromagnetic field at the waveguide boundary(29)Ezx+0=Ezx-0,Hzx+0=Hzx-0,(30)Eϕx+0=Eϕx-0,Hϕx+0=Hϕx-0,x≡k0Rwe have(31)Cmαmβm=0,αm=AmJmrztxt=amKmxc,βm=BmJmxt=bmKmxc,where xt=qtx, xc=qcx and the matrix Cm is given by(32)Cm=qzmnt2-nc2μqc2RJmxt+μcqt2RKmxcϵtqc2RJmrztxt+ϵcqt2RKmxcqzmnt2-nc2,(33)RJmx=xJm′xJmx=m-xJm+1xJmx,RKmx=xKm′xKmx=m-xKm+1xKmx.Dispersion relations for the harmonics of the electromagnetic field, linking kz and k0=ω/c, can be deduced from the singularity condition for the matrix (32):(34)det⁡Cmqz,x=0,nc2≤qz2=kzk02≤nt2,x=k0R.


## 2.3. Numerical Results

We begin our analysis with the case of waveguide modes with the zero azimuthal number:m=0. In this case, the dispersion relations are simplified and reduced to the relations for two types of the waves: TE waves with Ez=0(35)μqc2RJ0xt=-μcqt2RK0xc⇒μJ1xtxtJ0xt=-μcK1xcxcK0xc and TH waves with Hz=0(36)ϵtqc2RJ0rztxt=-ϵcqt2RK0xc⇒ϵzJ1rztxtrztxtJ0rztxt=-ϵcK1xcxcK0xc.The roots of Eqs. (35) and (36) are located between the corresponding zeros of the Bessel function J0 and will be numbered in ascending order by the index l. In our subsequent calculations, we shall use the estimates taken from the paper [11]: the pores of radius R that ranged from 200 nm to 600 nm are filled with the nematic liquid crystal 5CB characterized by the refractive indices ne≈1.72 and no≈1.5 (μ=1), the effective refractive index of the polymer matrix is nc≈1.57 (μc=1), and the wavelength of the laser radiation is λ=632 nm. Clearly, in this case, we have 2≈1.98≤k0R≤5.96≈6.Figure1 shows how liquid crystal anisotropy (the value of the longitudinal refractive index nz) affects the dispersion curves for the waveguide modes TE0l and TH0l. It can be seen that the curves for the TE modes, which are computed for the planar director configuration with the effective refractive index nt=ne=1.72, are independent of nz (see Eq. (35)), whereas the difference between the dispersion curves for the TH and TE modes becomes more pronounced as the difference between the refractive indices nt and nz increases. The radial distributions of nonvanishing components of the electromagnetic fields for the waveguide modes TH01 and TE01 are presented in Figures 2 and 3, respectively.Figure 1
Dependence of the dimensionless parameterk0R≡ωR/c on qz=kz/k0 for TEml (Ez=0) and THml (Hz=0) waves with m=0: (a) l=1 and (b) l=2 (index l enumerates the branches of roots in ascending order) at different values of the refractive index nz.
(a)
m = 0 and l=1 (b)
m = 0 and l=2Figure 2
Electromagnetic field versus the dimensionless radial parameterr=k0ρ for the mode TH01 (Hz=0) at (a) nz=1.4 and (b) nz=1.5.
(a)
m = 0 and nz=1.4 (b)
m = 0 and nz=1.5Figure 3
Electromagnetic field versus the dimensionless radial parameterr=k0ρ for the mode TE01 (Ez=0) at (a) nz=1.4 and (b) nz=1.5.
(a)
m = 0 and nz=1.4 (b)
m = 0 and nz=1.5The results illustrating effect ofnz on the hybrid modes (Ez≠0 and Hz≠0) with m=1 and m=2 are plotted in Figure 4. It is seen that as compared to the TH modes the curves for the hybrid modes are much less sensitive to variations in nz. Interestingly, for such modes, similar to the case of TH01 and TE01 modes, the waveguide regime may occur in the porous films, where 2≤k0R≤6.Figure 4
Dependence of the dimensionless parameterk0R≡ωR/c on qz=kz/k0 for the waves with (a) m=1 and (b) m=2 at different values of the refractive index nz.
(a)
m = 1 (b)
m = 2When the electric field is applied, the effective refractive indices of the LC (in a typical situationnt decreases, whereas nz increases) and the waveguide regime in the porous film can be suppressed. This is obvious if the refractive index nt, which approaches no<nc as the field increases, becomes less than nc, so that the condition of total internal reflection appears to be violated. Figure 5 shows the dispersion curves for the case when the configuration distorted by the field is close to axial with nz=ne, but nt is still higher than nc. It is interesting that under these conditions, for a film with 2≤k0R≤6, the conditions for the waveguide propagation regime are satisfied only for the hybrid mode with m=1.Figure 5
Dimensionless parameterk0R≡ωR/c as a function of qz=kz/k0 for the waves with 0≤m≤2 at nz=1.72>nt=1.6≈nc.Radial distributions for the components of the electric and magnetic fields are shown in Figures6 and 7, respectively. These figures, in particular, demonstrate that the degree of the field localization inside the pore decreases as the director configuration approaches the axial structure.Figure 6
Components of electromagnetic field as a function of the dimensionless radial parameterr=k0ρ for the mode with m=1 at (a) nz=1.5 and (b) nz=1.72.
(a)
m = 1 and nz=1.5 (b)
m = 1 and nz=1.72Figure 7
Components of electromagnetic field as a function of the dimensionless radial parameterr=k0ρ for the mode with m=1 at (a) nz=1.5 and (b) nz=1.72.
(a)
m = 1 and nz=1.5 (b)
m = 1 and nz=1.72


## 3. Experiment

### 3.1. Experimental Procedure

In our experiments, we have used two types of PET films that are characterized by both micron-sized (the diameter isd=5μm) and submicron-sized (d=0.6μm) pores without any surfactant treatment. The PET films of thickness h=23μm were obtained by heavy ion bombarding with subsequent chemical etching. Axes of open-end pores were normal to the surface of the films. The PET films were rinsed two times with ethanol, blown, and then vacuumed to remove a solvent. The films were filled with nematic liquid crystal (4-n-pentyl-4’-cyanobiphenyl, 5CB) in isotropic phase and then slowly cooled down to room temperature. Afterwards, the films were wiped to remove parasitic LC layers on their surfaces.Two different cells and experimental set-ups were used for the optical study of light propagation through the PET films filled with 5CB. In the case of the film with submicron pores, we measured the integral intensity of He-Ne laser (the wavelength isλ=0.633μm) beam passed through the PET film. The linearly polarized beam was normally incident on the PET film sandwiched between two glass substrates with ITO electrodes and placed on an optical bench. The intensity of light was measured at fixed angles β between the optical axis of the film and the polarization plane of the beam using the scheme without an analyzer to avoid the birefringence effects.In the films with micron pores, we studied the optical textures of 5CB within the individual micron-sized pores of the PET film. In this case, the film was stacked between two glass plates. The lower glass plate contained ITO layer with inter-electrode gap of the width g=50μm needed to arrange in-plane electric field, applied to the sample. Polarized optical microscopy (POM) with 1000x magnification was used to study optical textures of 5CB in both crossed and parallel polarizers for the LC cell differently oriented with respect to them.


### 3.2. Experimental Results and Discussion

The electrooptical response of the sample at different frequencies of the applied field is shown in Figure8. At high frequency (f=5 kHz), it looks like a typical response of planar layers of LC with the positive dielectric permittivity anisotropy in crossed polarizers, oriented at 45∘ with respect to the optical axis. So, the dielectric mechanism of action of electric field on LC, giving rise to overall rotation of the director towards the pore’s axis, is responsible for the observed effect.Figure 8
The dependence of a light intensityI passed through the PET-5CB film as a function of the dimensionless time parameter t/T at different frequencies f=0.3…5 kHz of the applied electric field: T=1/2f; the voltage U=170 V.In particular, the field induced distortions of the initial configuration results in the decrease of the effective anisotropy of LC, introduced by the theoretical model, and disappearance of the waveguide modes, sensitive to this parameter. It leads to decreasing of the light intensity registered in the experiment. In the case of strong electric fields, one can expect the breaking of surface anchoring [2] and formation of the axial LC configuration. This configuration is characterized by the ordinary refractive index which is lower than the indexes of the polymer matrix, so that the waveguide regime appears to be suppressed. The effective value of the electric field strength (of order 1 V/μm) is comparable with the critical values for such orientational transformation obtained via the computer simulation [2] for a nematic mixture E7 with material parameters close to those for 5CB.So, the breaking of the surface anchoring with corresponding light scattering can be responsible for strong decrease of the light intensity observed in our experiments. It is also important to notice that such type of response is drastically different from the low frequency response shown in Figure8. This response can be explained by the additional influence of ion’s motion [12]. Nevertheless, the strong decrease of the instant values of light intensity, registered in this case, can be also explained by the abovementioned mechanism.We also directly studied the propagation of light through the single micrometer-sized pores (d=5μm) of the PET film filled with LC by analyzing the optical textures. Figure 9(a) shows the microscopic images of the PET-LC film differently oriented with respect to the crossed (parallel) polarizers. For the crossed polarizers, the majority of pores reveal the cross-like texture indicating orientation of the director field either along or perpendicular to the dark lines. Such behavior may take place as a result of different director configurations within pores.Figure 9
(a) Microscopic images of a single pore of the PET film filled with LC and (b) the intensity profiles measured along the dashed lines for the film rotated by different angles with respect to the polarizers.
(a) (b)In contrast to the case of crossed polarizers, a bright state appears for the film rotated about the normal by 45 degrees. This points to the existence of the symmetry axis for the director configuration, which is typical for both planar (the planar polar configuration) and homeotropic (the axial polar configuration) anchoring. At the same time, the dark-to-bright transition as the film is rotated with respect to the crossed polarizers results in the formation of the ring in the central region of the pores. It should also be noted that the bright spot surrounded by the dark area within the pores was observed in the parallel polarizers (see the lower panel of Figure9(a)). Usually observation of the dark areas in the parallel polarizers points to the twisting orientational LC structure. Observation of the bright spot surrounded by the dark ring can be explained as a result of waveguide propagation of light in the pore. Since the defectless cross-like figure in crossed polarizers corresponds to the uniform LC structure within the pores, observation of the dark area surrounding the bright spot indicates violation of the propagation of polarization light. In this case, the distribution of brightness (the light intensity as we assume in this study) strongly depends on the polarization state of light (see Figure 9(b)).The dark-to-bright transition was also caused by the in-plane ac electric field applied to the PET-LC film (see Figure10(a)). In this case, the electric field induced the twist torque that reorients the bulk director of LC within the pores in the plane of the PET film. It produces the rotation of the symmetry axis of the LC structure along the field lines. Therefore, the symmetry axis of LC configuration twisted from the side of the cell with inter-electrode gap to the glass substrate without ITO that mimics the helicoid LC structure within the pores.Figure 10
(a) Microscopic images of a single pore filled with LC for different values of the applied voltage (the voltageU is ranged from 0 to 120 V and f=3 kHz) and (b) the intensity profiles measured along the dashed line.
(a) (b)Since the electric coherence length (the length that characterizes the thickness of a boundary layer, with the orientation stabilized by surface in the presence of electric field) decreases with the voltage, it explains the increase in the intensity of light passed through the single pore (see Figure10(b)).From the above discussion, we arrive at the conclusion that the waveguide light propagation strongly depends on the orientational LC configuration formed inside the pores and can be effectively controlled by applying either in-plane or out-of-plane electric fields to the PET-LC film.


## 3.1. Experimental Procedure

In our experiments, we have used two types of PET films that are characterized by both micron-sized (the diameter isd=5μm) and submicron-sized (d=0.6μm) pores without any surfactant treatment. The PET films of thickness h=23μm were obtained by heavy ion bombarding with subsequent chemical etching. Axes of open-end pores were normal to the surface of the films. The PET films were rinsed two times with ethanol, blown, and then vacuumed to remove a solvent. The films were filled with nematic liquid crystal (4-n-pentyl-4’-cyanobiphenyl, 5CB) in isotropic phase and then slowly cooled down to room temperature. Afterwards, the films were wiped to remove parasitic LC layers on their surfaces.Two different cells and experimental set-ups were used for the optical study of light propagation through the PET films filled with 5CB. In the case of the film with submicron pores, we measured the integral intensity of He-Ne laser (the wavelength isλ=0.633μm) beam passed through the PET film. The linearly polarized beam was normally incident on the PET film sandwiched between two glass substrates with ITO electrodes and placed on an optical bench. The intensity of light was measured at fixed angles β between the optical axis of the film and the polarization plane of the beam using the scheme without an analyzer to avoid the birefringence effects.In the films with micron pores, we studied the optical textures of 5CB within the individual micron-sized pores of the PET film. In this case, the film was stacked between two glass plates. The lower glass plate contained ITO layer with inter-electrode gap of the width g=50μm needed to arrange in-plane electric field, applied to the sample. Polarized optical microscopy (POM) with 1000x magnification was used to study optical textures of 5CB in both crossed and parallel polarizers for the LC cell differently oriented with respect to them.


## 3.2. Experimental Results and Discussion

The electrooptical response of the sample at different frequencies of the applied field is shown in Figure8. At high frequency (f=5 kHz), it looks like a typical response of planar layers of LC with the positive dielectric permittivity anisotropy in crossed polarizers, oriented at 45∘ with respect to the optical axis. So, the dielectric mechanism of action of electric field on LC, giving rise to overall rotation of the director towards the pore’s axis, is responsible for the observed effect.Figure 8
The dependence of a light intensityI passed through the PET-5CB film as a function of the dimensionless time parameter t/T at different frequencies f=0.3…5 kHz of the applied electric field: T=1/2f; the voltage U=170 V.In particular, the field induced distortions of the initial configuration results in the decrease of the effective anisotropy of LC, introduced by the theoretical model, and disappearance of the waveguide modes, sensitive to this parameter. It leads to decreasing of the light intensity registered in the experiment. In the case of strong electric fields, one can expect the breaking of surface anchoring [2] and formation of the axial LC configuration. This configuration is characterized by the ordinary refractive index which is lower than the indexes of the polymer matrix, so that the waveguide regime appears to be suppressed. The effective value of the electric field strength (of order 1 V/μm) is comparable with the critical values for such orientational transformation obtained via the computer simulation [2] for a nematic mixture E7 with material parameters close to those for 5CB.So, the breaking of the surface anchoring with corresponding light scattering can be responsible for strong decrease of the light intensity observed in our experiments. It is also important to notice that such type of response is drastically different from the low frequency response shown in Figure8. This response can be explained by the additional influence of ion’s motion [12]. Nevertheless, the strong decrease of the instant values of light intensity, registered in this case, can be also explained by the abovementioned mechanism.We also directly studied the propagation of light through the single micrometer-sized pores (d=5μm) of the PET film filled with LC by analyzing the optical textures. Figure 9(a) shows the microscopic images of the PET-LC film differently oriented with respect to the crossed (parallel) polarizers. For the crossed polarizers, the majority of pores reveal the cross-like texture indicating orientation of the director field either along or perpendicular to the dark lines. Such behavior may take place as a result of different director configurations within pores.Figure 9
(a) Microscopic images of a single pore of the PET film filled with LC and (b) the intensity profiles measured along the dashed lines for the film rotated by different angles with respect to the polarizers.
(a) (b)In contrast to the case of crossed polarizers, a bright state appears for the film rotated about the normal by 45 degrees. This points to the existence of the symmetry axis for the director configuration, which is typical for both planar (the planar polar configuration) and homeotropic (the axial polar configuration) anchoring. At the same time, the dark-to-bright transition as the film is rotated with respect to the crossed polarizers results in the formation of the ring in the central region of the pores. It should also be noted that the bright spot surrounded by the dark area within the pores was observed in the parallel polarizers (see the lower panel of Figure9(a)). Usually observation of the dark areas in the parallel polarizers points to the twisting orientational LC structure. Observation of the bright spot surrounded by the dark ring can be explained as a result of waveguide propagation of light in the pore. Since the defectless cross-like figure in crossed polarizers corresponds to the uniform LC structure within the pores, observation of the dark area surrounding the bright spot indicates violation of the propagation of polarization light. In this case, the distribution of brightness (the light intensity as we assume in this study) strongly depends on the polarization state of light (see Figure 9(b)).The dark-to-bright transition was also caused by the in-plane ac electric field applied to the PET-LC film (see Figure10(a)). In this case, the electric field induced the twist torque that reorients the bulk director of LC within the pores in the plane of the PET film. It produces the rotation of the symmetry axis of the LC structure along the field lines. Therefore, the symmetry axis of LC configuration twisted from the side of the cell with inter-electrode gap to the glass substrate without ITO that mimics the helicoid LC structure within the pores.Figure 10
(a) Microscopic images of a single pore filled with LC for different values of the applied voltage (the voltageU is ranged from 0 to 120 V and f=3 kHz) and (b) the intensity profiles measured along the dashed line.
(a) (b)Since the electric coherence length (the length that characterizes the thickness of a boundary layer, with the orientation stabilized by surface in the presence of electric field) decreases with the voltage, it explains the increase in the intensity of light passed through the single pore (see Figure10(b)).From the above discussion, we arrive at the conclusion that the waveguide light propagation strongly depends on the orientational LC configuration formed inside the pores and can be effectively controlled by applying either in-plane or out-of-plane electric fields to the PET-LC film.


## 4. Conclusion

In this paper we have theoretically studied the light propagation in a cylindrical waveguide filled with the material characterized by the biaxial radial optical anisotropy. The analysis of the obtained dispersion relations for the waveguide modes showed that, when the electric field is applied to the porous film, the waveguide propagation regime is possible for the hybrid mode with the azimuthal numberm=1, until the total internal reflection condition is violated. We have also presented some experimental results on the waveguide regime of light propagation in polymer PET porous films with micro-sized and submicro-sized pores filled with the nematic liquid crystal 5CB in the presence of electric fields. In particular, we have found that the pronounced decrease of the transmitted light intensity induced by the strong electric field applied along the submicro-sized pore’s axis can be explained by the suppression of the propagating waveguide modes resulting from the surface anchoring breaking and the formation of the axial configuration. At the same time, in-plane electric field, applied to the micro-sized pores, results in essential changes of the inhomogeneous light intensity distribution across the pore’s cross-section, which can be assigned to the electrically induced twist deformation of the initial director configuration.The part of our analysis was concerned with the special case whereϵρ=ϵϕ≡ϵt. It relies on the assumption that the optical properties of the pores can be described by the effective dielectric tensor averaged over the LC director distribution: ϵij(eff)=ϵoδij+(ϵe-ϵo)〈didj〉. In [2], this approximation was applied to the escaped radial configuration in 5CB cells with the director field of the form d^=cos⁡χ(ρ)e^ρ+sin⁡χ(ρ)z^ characterized by the angle χ(ρ). In this case, the principle values of the effective dielectric tensor are ϵt=ϵo+(ϵe-ϵo)〈cos2⁡χ〉/2 and ϵzz=ϵo+(ϵe-ϵo)〈sin2χ〉, where 〈cos2χ〉=2R-2∫0R‍cos2χ(ρ)ρdρ.Heuristically, validity of the above approximation can be justified using the condition of weakly inhomogeneous anisotropy that can be written in the form(37)k0Rδn<1,where δn=max⁡ϵeff-ϵ is the refractive index variation characterizing the difference between the effective and local values of the refractive indices. For the escaped radial configuration in the 5CB cell without applied electric field, it is not difficult to estimate the refractive index variation at about δn≈0.1. In deriving this estimate, we have assumed strong homeotropic anchoring conditions so that cos⁡χ(R) is close to unity. Electric field applied across the film and weak anchoring conditions are the factors that may considerably reduce δn. Our results, however, indicate that, for micron-sized pores, the approximation based on the effective dielectric tensor breaks down.Our concluding remark is that, in addition to the description of the waveguide regime in optically anisotropic cylindrical waveguides, our analytical results can be used to extend the class of exactly solvable problems of the theory of light scattering. Further investigations suggest solving the problem of light transmission for optically anisotropic porous films that, in particular, will require a more sophisticated analysis using the methods of the effective medium theory.



---

*Source: 1539865-2019-02-13.xml*