# Evanescent Field Enhancement in Liquid Crystal Optical Fibers: A Field Characteristics Based Analysis

Authors: P. K. Choudhury
Journal: Advances in Condensed Matter Physics
 (2013)

---

## Abstract

The paper presents the analysis of the electromagnetic wave propagation through liquid crystal optical fibers (LCOFs) of two different types—conventional guides loaded with liquid crystals (addressed as LCOFs) and those with additional twists due to conducting helical windings (addressed as HCLCOFs). More precisely, the three-layer optical waveguide structures are considered along with its outermost region being loaded with radially anisotropic liquid crystal material and the inner regions being made of usual silica, as used in conventional optical fibers. In addition to that, LCOF with twists introduced in the form of conducting helical windings at the interface of the silica core and the liquid crystal clad is also taken into account. Emphasis has been put on the power confinements by the lower-order TE modes sustained in the different sections of the LCOF structure. The results demonstrate useful applications of these guides in integrated optics as the power sustained in the liquid crystal section by the excited TE modes remains very high. In the case of twisted clad liquid crystal guides, descriptions are limited to the nature of dispersion relation only under the TE mode excitation, and corresponding to the cases of helix orientations being parallel and perpendicular to the optical axis.


---

## Body

## 1. Introduction

Optical fibers are of immense use in communication technology, which may be either for the long distance use, for example telecommunication needs, or the short distance ones, for example sensing and so forth. In either case, optical fibers can be designed to act as active and/or passive device. The propagation characteristics of various types of optical fibers have been discussed in the literature. In this respect, the fiber itself may either have wide varieties of geometries ranging from conventional to the unconventional ones, or the guide may possess different forms of material compositions. The characteristics of electromagnetic wave propagation through optical fibers are greatly dependent on the type of material which the guide is composed of.Among the others, fibers made of liquid crystals [1–3] fall into the category of complex mediums and proved to be greatly attractive for multifarious applications [4, 5] because these materials present the distinct characteristic of anisotropic nature in their physical properties [6, 7]. Investigators have reported liquid crystal optical fibers (LCOFs) of different forms of constructions in respect of geometry as well as material distributions. LCOFs may consist of liquid crystals either in the core [8–10] or in the clad section [9–11] of the guide. High optical anisotropy of liquid crystals makes them as demanding smart materials for their usage in a wide range of applications in the area of integrated optics [7, 12]. These can be used in electric field sensing as the macroscopic optical properties of liquid crystals can be altered under the influence of external electrical field [13]. Yoshino et al. [14] demonstrated the use of LCOFs in electric field sensing based on the inherent property of liquid crystals that these exhibit the largest electrooptic effect among known materials. In all the aforesaid applications, the clad region of the guide consists of liquid crystal. For the fabrication of directional couplers, LCOFs are used with the liquid crystal section sandwiched between two side polished fibers, and coaxial fibers are embedded in liquid crystals.Applications of liquid crystals make it clear that these exhibit distinct feature of polarization anisotropy which may have two different directions of orientation-radial and azimuthal anisotropies; the latter one is less discussed in the literature. In the present paper, we consider the nematic form of liquid crystals along with its radially anisotropic nature to constitute the fiber clad, the core region being made of standard silica. In forming such LCOFs, the radial anisotropy of liquid crystal can be obtained by capillary action after inserting the liquid crystal section (of the fiber) into a capillary tube coated with N,N-dimethyl-N-octadecyl-3-aminopropyltrimethoxysilyl chloride [11].Considering the importance of the usage of liquid crystals in optical fibers, such guides with tapered structures have been discussed before [15]. Since the propagation of power through a guide plays the determining role for specific applications of the guide, the scope of the present paper lies in the analysis of LCOFs in respect of power confinement factor. We consider a three-layer LCOF, of which the outermost clad is made of radially anisotropic liquid crystal material, and determine the power propagation characteristics corresponding to two lower-order TE modes in the guides of different dimensions. It has been observed that the TE modes transport relatively high amount of power through the liquid crystal clad. These interesting features of radially anisotropic LCOFs indicate their usefulness in the fabrication of different coupling devices used in the area of integrated optics, and also, in optical sensing where the method of evanescent field absorption is primarily implemented.The investigations are extended to the case of LCOFs with conducting helical windings (i.e., HCLCOFs—the helical clad liquid crystal optical fibers) introduced at the core-clad interface of the guide. In such fibers, the angle of pitch (of the helix in use) exhibits dominant characteristics to alter the wave propagation, and the working principles are based on the theory of low- and medium-power travelling wave tubes (TWTs) [16]. Implementing the concept in the case of dielectric optical fibers, studies have appeared in the literature for fibers having circular and elliptical cross-sections [17–20]. However, in the present paper, as stated before, the existence of the liquid crystal layer will modify the propagation characteristics. The dispersion relation for such types of guides (i.e. HCLCOFs) is deduced, and the features are reported corresponding to the cases of helical turns being parallel as well as perpendicular to the optical axis.


## 2. The Analytical Approach

### 2.1. The Case of LCOFs

Figure1 illustrates the cross-sectional view of the three-layer LCOF, of which the (infinitely extended) outermost clad is made of radially anisotropic liquid crystal material; the core and the inner clad regions being constructed of linear, homogeneous, and isotropic dielectrics having the refractive index (RI) values as n1 and n2 (n1>n2). Also, no and ne, respectively, are the ordinary and the extraordinary RIs of the liquid crystal section. Further, ρ1 and ρ2  are, respectively, the radii of the core and the inner clad regions. The dashed lines in the radial direction (in Figure 1) indicate the radial orientation of liquid crystal molecules. We consider the time harmonic waves with angular frequency ω (in the unbounded medium) propagate along the z-direction with the propagation constant β, and the wave is harmonic in the direction of propagation too. Of the liquid crystal section, the principal axes coincide with the z-axis, and the extraordinary principal axis is radially oriented. Therefore, the liquid crystal region will have the RI distribution as
(1)nρ=ne,nϕ=nzwithne>n1>n2>no,
where nρ, nϕ, and nz are, respectively, the RI values along the ρ-, ϕ-, and z-directions. Now, following the formulations presented in [21], the coupled wave equations for the transverse field components can be written as [22](2a)(∇t2+k02nρ2-β2-1ρ2)eρ=-2ρ2∂eϕ∂ϕ+(1-nρ2nz2)∂∂ρ{1ρ∂∂ρ(ρeρ)}+(1-nϕ2nz2)∂∂ρ(1ρ∂eϕ∂ϕ)(2b)(∇t2+k02nϕ2-β2-1ρ2)eϕ=2ρ2∂eρ∂ϕ+1ρ2{(1-nρ2nz2)∂∂ρ(ρ∂eρ∂ϕ)+(1-nϕ2nz2)∂2eϕ∂ϕ2},where ∇t2 is the Laplacian operator in the cylindrical coordinate system and k0 is the free-space propagation constant.Figure 1
Cross-sectional view of a three-layer LCOF with the outermost section filled with liquid crystal (the dashed lines in the radial direction indicate the orientation of the liquid crystal molecules).Liquid crystal is anisotropic in nature, and therefore, the sustained modes in the LCOF under consideration will essentially contain the electric field componentsEρ,Eϕ, and Ez, each facing different RI values in their respective directions. However, in the present work, we consider the lower-order TE modes (namely, TE01). In this case, there is only one transverse electrical field component eϕ, which is independent of the coordinate ϕ, making thereby one to have eρ=0 and ∂eϕ/∂ϕ=0 corresponding to this mode. Implementing these conditions and using (2a) and (2b), after some mathematical steps, it can be shown that the electric/magnetic field components in this case will be ultimately given as(3a)Hρ)I=-C1βωμ0J1(γ1ρ)exp{j(ωt-βz)}(3b)Hz)I=C1jωμ0{γ1J1′(γ1ρ)+1ρJ1(γ1ρ)}exp{j(ωt-βz)}(3c)Hρ)II=-βωμ0{C2K1(γ2ρ)+C3I1(γ2ρ)}exp{j(ωt-βz)}(3d)Hz)II=jωμ0[C2{γ2K1′(γ2ρ)+1ρK1(γ2ρ)}+C3{wI1′(γ2ρ)+1ρI1(γ2ρ)}exp{j(ωt-βz)}](3e)Hρ)III=-C4βωμ0K1(γ3ρ)exp{j(ωt-βz)}(3f)Hz)III=C4jωμ0{γ3K1′(γ3ρ)+1ρK1(γ3ρ)}×exp{j(ωt-βz)}.In (3a), (3b), (3c), (3d), (3e), and (3f), the superscripts I,II, and III, respectively, represent the situations in the LCOF core, inner clad and the outer liquid crystal clad sections of the guide. Also, μ0 is the free-space permeability, C1,C2,C3, and C4 are the arbitrary constants to be determined by the boundary conditions, and J(·),K(·), and I(·) represent Bessel and the modified Bessel functions. Further, the prime indicates the differentiation with respect to the argument, and the quantities γ1,γ2, and γ3 are defined as(4a)γ1=n12k02-β2(4b)γ2=β2-n22k02(4c)γ3=β2-no2k02Making use of the above-mentioned field components, as stated in (3a), (3b), (3c), (3d), (3e), and (3f) and implementing the electromagnetic boundary conditions at the layer interfaces of the LCOF, values of the arbitrary constants can be determined. By using those, the expressions corresponding to power [23] transmitted through the different sections of the LCOF can be evaluated as
(5)Pcore=C12πωμ0{γ1∫0ρ1ρJ1(γ1ρ)J1′(γ1ρ)dρ+∫0ρ2{J1(γ1ρ)}2dρ}(6)Pinnerclad=C12πωμ0[Ψ2{γ2∫ρ1ρ2ρI1(γ2ρ)I1′(γ2ρ)dρ+∫ρ1ρ2(I1(γ2ρ))2dρ}+Ω2{γ2∫ρ1ρ2ρK1(γ2ρ)K1′(γ2ρ)dρ+∫ρ1ρ2(K1(γ2ρ))2dρ}+ΦΨ{2∫ρ1ρ2K1(γ2ρ)I1(γ2ρ)dρ+γ2∫ρ1ρ2ρK1(γ2ρ)I1′(γ2ρ)dρ+γ2∫ρ1ρ2ρI1(γ2ρ)K1′(γ2ρ)dρ}](7)Pouterclad=C12πωμ0(ΦK1(γ2ρ2)+ΨI1(γ2ρ2)K1(γ3ρ2))2×[∫ρ2∞{K1(γ3ρ)}2dρ+γ3∫ρ2∞ρK1(γ3ρ)K1′(γ3ρ)dρ],
where
(8)Φ=(J1(γ1ρ1){γ2I1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-I1(γ2ρ1){γ1J1′(γ1ρ1)+1ρ1J1(γ1ρ1)})×(K1(γ2ρ1){γ2I1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-I1(γ2ρ1){γ2K1′(γ2ρ1)+1ρ1K1(γ2ρ1)})-1,(9)Ψ=(J1(γ1ρ1){γ2K1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-K1(γ2ρ1){γ1J1′(γ1ρ1)+1ρ1J1(γ1ρ1)})×(I1(γ2ρ1){γ2K1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-K1(γ2ρ1){γ2I1′(γ2ρ1)+1ρ1K1(γ2ρ1)})-1.
Above derived (5), (6), and (7), respectively, represent the expressions of the power propagating through the LCOF core, inner dielectric clad, and the outer liquid crystal clad regions. It is to be noted at this point that the evaluation of power requires the values of all the arbitrary constants C1,C2,C3, and C4 to be written in terms of only one constant (say C1), which can be determined by a normalization condition taking into consideration the input power. If Ptotal represents the total power transmitted through the LCOF by the TE01 modes, that is,
(10)Ptotal=Pcore+Pinnerclad+Pouterclad,
then Pcore/Ptotal(≡Γcore), Pinnerclad/Ptotal (≡Γinnerclad), and Pouterclad/Ptotal (≡Γouterclad) will, respectively, determine the relative power (or the power confinement factor) in the LCOF core, inner clad, and the outer clad sections.


### 2.2. The Case of HCLCOFs

Figure2 shows the longitudinal view of the HCLCOF under consideration, which incorporates some modifications in the previous guide structure illustrated in Figure 1 in the form of additional conducting helical windings loaded at the interface of the core and the inner dielectric clad sections. As the fiber core is wrapped with a conducting sheath helix at a constant angle (called as the helix pitch angle) ψ around the core-inner clad boundary, the structure essentially has a circular cylindrical surface with high conductivity in a preferential direction. This way, helical clad fibers offer the additional parameter to suitably control the modal behaviour of the guide, that is, by varying the pitch angle of helix—an aspect of great technological importance [17–20].Figure 2
Longitudinal view of the HCLCOF under consideration.Coming to the analytical steps, the field components in the different sections of the HCLCOFs will be represented by (3a), (3b), (3c), (3d), (3e), and (3f). However, the boundary conditions will be modified in this case [24], the implementation of which will finally provide an 8×8 determinant formed by the coefficients of the unknown constants. A nontrivial solution to the set of equations can be obtained if the determinant is equated to zero, that is,
(11)(x1x200000000x3x4x5x600x7-x8-x9-x10x11x1200-x13x14x15x16-x17-x180000y1y2y3y400000000y5y600y7y8-y9-y10-y11y1200-y13-y14y15y16y17-y18)×(ABCDEFGH)=(00000000).
In  (11), the used symbols have their meanings as
(12)x1=(sinψ+νβγ12ρacosψ)Jν(γ1ρa),x2=iωμ0γ1cosψJν′(γ1ρa),x3=(sinψ+νβγ22ρacosψ)Kν(γ2ρa),x4=(sinψ+νβγ22ρacosψ)Iν(γ2ρa),x5=iωμ0γ2cosψKν′(γ2ρa),x6=iωμ0γ2cosψIν′(γ2ρa),x7=(cosψ-νβγ12ρasinψ)Jν(γ1ρa),x8=iωμ0γ1Jν′(γ1ρa),x9=(cosψ-νβγ22ρasinψ)Kν(γ2ρa),x10=(cosψ-νβγ22ρasinψ)Iν(γ2ρa),x11=iωμ0γ2sinψKν′(γ2ρa),x12=iωμ0γ2sinψIν′(γ2ρa),x13=iωε1γ1cosψJν′(γ1ρa),x14=(sinψ+νβγ12ρacosψ)Jν(γ1ρa),x15=iωε2γ2cosψKν′(γ2ρa),x16=iωε2γ2cosψIν′(γ2ρa),x17=(sinψ+νβγ22ρacosψ)Kν(γ2ρa),x18=(sinψ+νβγ22ρacosψ)Iν(γ2ρa),y1=(sinψ+νβγ22ρbcosψ)Kν(γ2ρb),y2=(sinψ+νβγ22ρbcosψ)Iν(γ2ρb),y3=iωμ0γ2cosψKν′(γ2ρb),y4=iωμ0γ2cosψIν′(γ2ρb),y5=(sinψ+νβγ32ρbcosψ)Kν(γ3ρb),y6=iωμ0γ3cosψKν′(γ3ρb),y7=(cosψ-νβγ22ρbsinψ)Kν(γ2ρb),y8=(cosψ-νβγ22ρbsinψ)Iν(γ2ρb),y9=iωμ0γ2Kν′(γ2ρb),y10=iωμ0γ2Iν′(γ2ρb),y11=(cosψ-νβγ32ρbsinψ)Kν(γ3ρb),y12=iωμ0γ3sinψKν′(γ3ρb),y13=iωε2γ2cosψKν′(γ2ρb),y14=iωε2γ2cosψKν′(γ2ρb),y15=(sinψ+νβγ22ρbcosψ)Kν(γ2ρb),y16=(sinψ+νβγ22ρbcosψ)Iν(γ2ρb),y17=iωε3γ3cosψKν′(γ3ρb),y18=(sinψ+νβγ32ρbscosψ)Kν(γ3ρb).
Equation (11) is the eigenvalue equation corresponding to the HCLCOF structure, solutions to which will provide the modes sustained in the guide.


## 2.1. The Case of LCOFs

Figure1 illustrates the cross-sectional view of the three-layer LCOF, of which the (infinitely extended) outermost clad is made of radially anisotropic liquid crystal material; the core and the inner clad regions being constructed of linear, homogeneous, and isotropic dielectrics having the refractive index (RI) values as n1 and n2 (n1>n2). Also, no and ne, respectively, are the ordinary and the extraordinary RIs of the liquid crystal section. Further, ρ1 and ρ2  are, respectively, the radii of the core and the inner clad regions. The dashed lines in the radial direction (in Figure 1) indicate the radial orientation of liquid crystal molecules. We consider the time harmonic waves with angular frequency ω (in the unbounded medium) propagate along the z-direction with the propagation constant β, and the wave is harmonic in the direction of propagation too. Of the liquid crystal section, the principal axes coincide with the z-axis, and the extraordinary principal axis is radially oriented. Therefore, the liquid crystal region will have the RI distribution as
(1)nρ=ne,nϕ=nzwithne>n1>n2>no,
where nρ, nϕ, and nz are, respectively, the RI values along the ρ-, ϕ-, and z-directions. Now, following the formulations presented in [21], the coupled wave equations for the transverse field components can be written as [22](2a)(∇t2+k02nρ2-β2-1ρ2)eρ=-2ρ2∂eϕ∂ϕ+(1-nρ2nz2)∂∂ρ{1ρ∂∂ρ(ρeρ)}+(1-nϕ2nz2)∂∂ρ(1ρ∂eϕ∂ϕ)(2b)(∇t2+k02nϕ2-β2-1ρ2)eϕ=2ρ2∂eρ∂ϕ+1ρ2{(1-nρ2nz2)∂∂ρ(ρ∂eρ∂ϕ)+(1-nϕ2nz2)∂2eϕ∂ϕ2},where ∇t2 is the Laplacian operator in the cylindrical coordinate system and k0 is the free-space propagation constant.Figure 1
Cross-sectional view of a three-layer LCOF with the outermost section filled with liquid crystal (the dashed lines in the radial direction indicate the orientation of the liquid crystal molecules).Liquid crystal is anisotropic in nature, and therefore, the sustained modes in the LCOF under consideration will essentially contain the electric field componentsEρ,Eϕ, and Ez, each facing different RI values in their respective directions. However, in the present work, we consider the lower-order TE modes (namely, TE01). In this case, there is only one transverse electrical field component eϕ, which is independent of the coordinate ϕ, making thereby one to have eρ=0 and ∂eϕ/∂ϕ=0 corresponding to this mode. Implementing these conditions and using (2a) and (2b), after some mathematical steps, it can be shown that the electric/magnetic field components in this case will be ultimately given as(3a)Hρ)I=-C1βωμ0J1(γ1ρ)exp{j(ωt-βz)}(3b)Hz)I=C1jωμ0{γ1J1′(γ1ρ)+1ρJ1(γ1ρ)}exp{j(ωt-βz)}(3c)Hρ)II=-βωμ0{C2K1(γ2ρ)+C3I1(γ2ρ)}exp{j(ωt-βz)}(3d)Hz)II=jωμ0[C2{γ2K1′(γ2ρ)+1ρK1(γ2ρ)}+C3{wI1′(γ2ρ)+1ρI1(γ2ρ)}exp{j(ωt-βz)}](3e)Hρ)III=-C4βωμ0K1(γ3ρ)exp{j(ωt-βz)}(3f)Hz)III=C4jωμ0{γ3K1′(γ3ρ)+1ρK1(γ3ρ)}×exp{j(ωt-βz)}.In (3a), (3b), (3c), (3d), (3e), and (3f), the superscripts I,II, and III, respectively, represent the situations in the LCOF core, inner clad and the outer liquid crystal clad sections of the guide. Also, μ0 is the free-space permeability, C1,C2,C3, and C4 are the arbitrary constants to be determined by the boundary conditions, and J(·),K(·), and I(·) represent Bessel and the modified Bessel functions. Further, the prime indicates the differentiation with respect to the argument, and the quantities γ1,γ2, and γ3 are defined as(4a)γ1=n12k02-β2(4b)γ2=β2-n22k02(4c)γ3=β2-no2k02Making use of the above-mentioned field components, as stated in (3a), (3b), (3c), (3d), (3e), and (3f) and implementing the electromagnetic boundary conditions at the layer interfaces of the LCOF, values of the arbitrary constants can be determined. By using those, the expressions corresponding to power [23] transmitted through the different sections of the LCOF can be evaluated as
(5)Pcore=C12πωμ0{γ1∫0ρ1ρJ1(γ1ρ)J1′(γ1ρ)dρ+∫0ρ2{J1(γ1ρ)}2dρ}(6)Pinnerclad=C12πωμ0[Ψ2{γ2∫ρ1ρ2ρI1(γ2ρ)I1′(γ2ρ)dρ+∫ρ1ρ2(I1(γ2ρ))2dρ}+Ω2{γ2∫ρ1ρ2ρK1(γ2ρ)K1′(γ2ρ)dρ+∫ρ1ρ2(K1(γ2ρ))2dρ}+ΦΨ{2∫ρ1ρ2K1(γ2ρ)I1(γ2ρ)dρ+γ2∫ρ1ρ2ρK1(γ2ρ)I1′(γ2ρ)dρ+γ2∫ρ1ρ2ρI1(γ2ρ)K1′(γ2ρ)dρ}](7)Pouterclad=C12πωμ0(ΦK1(γ2ρ2)+ΨI1(γ2ρ2)K1(γ3ρ2))2×[∫ρ2∞{K1(γ3ρ)}2dρ+γ3∫ρ2∞ρK1(γ3ρ)K1′(γ3ρ)dρ],
where
(8)Φ=(J1(γ1ρ1){γ2I1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-I1(γ2ρ1){γ1J1′(γ1ρ1)+1ρ1J1(γ1ρ1)})×(K1(γ2ρ1){γ2I1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-I1(γ2ρ1){γ2K1′(γ2ρ1)+1ρ1K1(γ2ρ1)})-1,(9)Ψ=(J1(γ1ρ1){γ2K1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-K1(γ2ρ1){γ1J1′(γ1ρ1)+1ρ1J1(γ1ρ1)})×(I1(γ2ρ1){γ2K1′(γ2ρ1)+1ρ1K1(γ2ρ1)}-K1(γ2ρ1){γ2I1′(γ2ρ1)+1ρ1K1(γ2ρ1)})-1.
Above derived (5), (6), and (7), respectively, represent the expressions of the power propagating through the LCOF core, inner dielectric clad, and the outer liquid crystal clad regions. It is to be noted at this point that the evaluation of power requires the values of all the arbitrary constants C1,C2,C3, and C4 to be written in terms of only one constant (say C1), which can be determined by a normalization condition taking into consideration the input power. If Ptotal represents the total power transmitted through the LCOF by the TE01 modes, that is,
(10)Ptotal=Pcore+Pinnerclad+Pouterclad,
then Pcore/Ptotal(≡Γcore), Pinnerclad/Ptotal (≡Γinnerclad), and Pouterclad/Ptotal (≡Γouterclad) will, respectively, determine the relative power (or the power confinement factor) in the LCOF core, inner clad, and the outer clad sections.


## 2.2. The Case of HCLCOFs

Figure2 shows the longitudinal view of the HCLCOF under consideration, which incorporates some modifications in the previous guide structure illustrated in Figure 1 in the form of additional conducting helical windings loaded at the interface of the core and the inner dielectric clad sections. As the fiber core is wrapped with a conducting sheath helix at a constant angle (called as the helix pitch angle) ψ around the core-inner clad boundary, the structure essentially has a circular cylindrical surface with high conductivity in a preferential direction. This way, helical clad fibers offer the additional parameter to suitably control the modal behaviour of the guide, that is, by varying the pitch angle of helix—an aspect of great technological importance [17–20].Figure 2
Longitudinal view of the HCLCOF under consideration.Coming to the analytical steps, the field components in the different sections of the HCLCOFs will be represented by (3a), (3b), (3c), (3d), (3e), and (3f). However, the boundary conditions will be modified in this case [24], the implementation of which will finally provide an 8×8 determinant formed by the coefficients of the unknown constants. A nontrivial solution to the set of equations can be obtained if the determinant is equated to zero, that is,
(11)(x1x200000000x3x4x5x600x7-x8-x9-x10x11x1200-x13x14x15x16-x17-x180000y1y2y3y400000000y5y600y7y8-y9-y10-y11y1200-y13-y14y15y16y17-y18)×(ABCDEFGH)=(00000000).
In  (11), the used symbols have their meanings as
(12)x1=(sinψ+νβγ12ρacosψ)Jν(γ1ρa),x2=iωμ0γ1cosψJν′(γ1ρa),x3=(sinψ+νβγ22ρacosψ)Kν(γ2ρa),x4=(sinψ+νβγ22ρacosψ)Iν(γ2ρa),x5=iωμ0γ2cosψKν′(γ2ρa),x6=iωμ0γ2cosψIν′(γ2ρa),x7=(cosψ-νβγ12ρasinψ)Jν(γ1ρa),x8=iωμ0γ1Jν′(γ1ρa),x9=(cosψ-νβγ22ρasinψ)Kν(γ2ρa),x10=(cosψ-νβγ22ρasinψ)Iν(γ2ρa),x11=iωμ0γ2sinψKν′(γ2ρa),x12=iωμ0γ2sinψIν′(γ2ρa),x13=iωε1γ1cosψJν′(γ1ρa),x14=(sinψ+νβγ12ρacosψ)Jν(γ1ρa),x15=iωε2γ2cosψKν′(γ2ρa),x16=iωε2γ2cosψIν′(γ2ρa),x17=(sinψ+νβγ22ρacosψ)Kν(γ2ρa),x18=(sinψ+νβγ22ρacosψ)Iν(γ2ρa),y1=(sinψ+νβγ22ρbcosψ)Kν(γ2ρb),y2=(sinψ+νβγ22ρbcosψ)Iν(γ2ρb),y3=iωμ0γ2cosψKν′(γ2ρb),y4=iωμ0γ2cosψIν′(γ2ρb),y5=(sinψ+νβγ32ρbcosψ)Kν(γ3ρb),y6=iωμ0γ3cosψKν′(γ3ρb),y7=(cosψ-νβγ22ρbsinψ)Kν(γ2ρb),y8=(cosψ-νβγ22ρbsinψ)Iν(γ2ρb),y9=iωμ0γ2Kν′(γ2ρb),y10=iωμ0γ2Iν′(γ2ρb),y11=(cosψ-νβγ32ρbsinψ)Kν(γ3ρb),y12=iωμ0γ3sinψKν′(γ3ρb),y13=iωε2γ2cosψKν′(γ2ρb),y14=iωε2γ2cosψKν′(γ2ρb),y15=(sinψ+νβγ22ρbcosψ)Kν(γ2ρb),y16=(sinψ+νβγ22ρbcosψ)Iν(γ2ρb),y17=iωε3γ3cosψKν′(γ3ρb),y18=(sinψ+νβγ32ρbscosψ)Kν(γ3ρb).
Equation (11) is the eigenvalue equation corresponding to the HCLCOF structure, solutions to which will provide the modes sustained in the guide.


## 3. Results and Discussion

The two types of fiber structures, namely, the LCOF and the HCLCOF, can now be analyzed. In the case of LCOF, the relative distributions of power in the three different regions of the guide are considered, whereas the analyses corresponding to HCLCOF are limited to the modal dispersion relations and the field cutoffs. Figures3(a), 3(b), and 3(c), respectively, illustrate the nature of power propagation in the LCOF core, inner clad, and the outer clad regions, as plotted against the propagation constant β corresponding to the situation of TE mode excitation in the guide. In our illustrative cases, we consider four different values of core radius for the computational purpose, namely, 5 μm, 15 μm, 30 μm, and 50 μm, and the clad radius is kept fixed at 125 μm. Further, the RI values of core and the inner clad regions are taken to be n1=1.462 and n2=1.458, respectively, and the operating wavelength λ0 is kept fixed as 633 nm. The outermost region of the guides (LCOF and HCLCOF) is considered to be the nematic liquid crystal named as BDH mixture 14616, which has the ordinary and the extraordinary RI values as no=1.457 and ne=1.5037, respectively.(a) Power confinement by the TE modes in the LCOF core. (b)   Power confinement by the TE modes in the LCOF inner clad. (c)  Power confinement by the TE modes in the LCOF outer clad.
(a)
(b)
(c)Now, focussing on the computations of the relative distribution of power in LCOF, we observe from Figure3(a) that the confinement of power shows an initial increase with the propagation constant β corresponding to almost all the fiber core dimensions considered here. However, the increase is more prominent in the case of LCOFs of larger core radius. We further find that, near the higher β-tail, power confinements are observed to be increasing with that of the core dimensions, which is very much expected as the fibers of larger dimensions essentially possess the capability of confining more amount of power in the fiber core. In the lower β-tail regime, we observe some fluctuations in power, and corresponding to 50 μm core radius, almost 59% of power is confined into the core section. In the medium range of the allowed propagation constants, the confinement patterns exhibit a kind of uniformity corresponding to lower dimensions of the LCOF. The largest core radius chosen here is 50 μm, and the confinement shows more fluctuations corresponding to this case. This is attributed to the fact that the LCOF becomes extremely multimode in this case, resulting thereby some sort of modal interference, which affects the uniformity in the distribution of power. The feature related to the kind of uniform power distribution over propagating modes is more required in telecommunication purposes, and LCOFs of smaller dimensions would possibly meet the requirements. However, more promising needs of LCOFs are expected to be in optical sensing or other integrated optic applications, which would generally require more amount of power to be confined in the outermost clad. Figure 3(a) shows the confinement to be almost about 15%–25% distributed over major portion of the allowed values of propagation constants.Figures3(b) and 3(c), respectively, correspond to the relative distribution of power in the LCOF inner dielectric and the outer liquid crystal clad sections. We observe that the power confinement patterns are almost similar in both of these sections, and also, the trends of power patterns are just opposite to what we observed in the LCOF core, which is quite expected as the way confinement increases or decreases (with the propagation constant) in the core section, it will decrease or increase, respectively, in the two clads. The noticeable fact remains that the average value of power in the two clad sections is appreciably higher (within the longer β-range) than what is observed in the LCOF core. More precisely, corresponding to lower LCOF core dimensions, within the uniform range of the allowed values of propagation constants, approximately 45% power is shared by each of the clad regions, whereas around only 10% power transmission is observed in the LCOF core at the same time. This essentially indicates a substantial increase of power in the two clad sections, which is attributed to the presence of nematic radially anisotropic liquid crystal material in the outermost clad region. The distribution of power, however, does not seem to be that uniform in the lower β-tail regime, and the uniformity is more affected with the increase in core dimension owing to the reason of possibilities related to modal interference in this case. Finally, the combined observations of the uniformity of power patterns corresponding to the LCOFs of smaller dimensions and higher power confinements in the outermost liquid crystal section (Figures 3(a) and 3(b)) essentially indicate the possibility of useful integrated optic applications of these guides, for example, in field coupling and/or optical sensing needs.Now, turning the attention towards the case of HCLCOFs, we make an attempt to study the dispersion relation expressed by (11). In this stream, we consider the values of the fiber core radius as 60.38 μm and the inner clad radius as 120.07 μm. The operating wavelength is chosen as 1.55 μm, and the RI values of the core and the inner clad sections are taken to be 1.5 and 1.46, respectively. The physical properties of the outermost liquid crystal section are considered to be the same as in the case of LCOF discussed above. Equation (11) is plotted corresponding to four different values of the azimuthal index v, namely, v=1, v=2, v=3, and v=4, and the introduced conducting helical turns may have perpendicular or parallel orientations with respect to the optical axis, making thereby the respective helical pitch angles as ψ=0° and ψ=90°.Figure4(a) illustrates the plots of the left hand side ℱ of (6) in the case when the helix pitch angles are perpendicular to the optical axis, that is, ψ = 0°. In this case, we observe that, corresponding to all values of the azimuthal indices, the plots intersect the ℱ=0 axis; the points of intersection essentially determine the propagation constant of the corresponding mode. We observe that the modal propagation constants corresponding to v=1 and v=3 are very close, though with the situation β)v=1≤β)v=2. Further, the modes with v=2 and v=4 attain almost similar β-values with β)v=2≤β)v=4.(a)   Plot of the dispersion relation for the TE modes corresponding toψ=0°.  (b) Plot of the cutoff characteristics for the TE modes corresponding to ψ=0°.
(a)
(b)Figure4(b) describes the situation when (11) is solved under the limit γ22→0, indicating thereby the respective plots (of ℱC) to correspond the case of field cutoffs. We observe that the cutoffs for the modes with v=1 and v=3 are almost the same, and that with v=2 and v=4 too are very much similar. However, we observe from Figures 4(a) and 4(b) that the cutoff β-values, as seen in Figure 4(b), are less than that obtained from the knowledge of β-values noticed from the intersections in Figure 4(a), which essentially concludes that the results illustrated in Figures 4(a) and 4(b)  are consistent.It is noteworthy that the helix pitch angleψ introduced in HCLCOF has the ability to alter the dispersion characteristics of the guide. In order to demonstrate this, plots are made corresponding to 90° pitch, and those corresponding to the dispersion equation and the field cutoffs are illustrated in Figures 5(a) and 5(b). In both these figures, profound effect due to helix pitch can clearly be observed. It is found that the β-values are roughly increased in this case. Making the pitch ψ=90° leaves the impression of eliminating the helical turns as those are now oriented just parallel to the optical axis. Thus, it can be conceptually understood that the presence of conducting helical turns (i.e., ψ=0°) introduces the modal propagation constants to be lower. It has been reported before that helix pitch makes the propagating waves to attenuate up to some extent—the property which can clearly be seen in the case of HCLCOF as the waves propagate with lower propagation constants under the perpendicular orientation (with respect to the optical axis) of helical turns.(a) Plot of the dispersion relation for the TE modes corresponding toψ=90°. (b)   Plot of the cutoff characteristics for the TE modes corresponding to ψ=90°.
(a)
(b)


## 4. Conclusion

From the foregoing discussions, inference can be drawn that, under TE mode excitation, LCOFs present higher amount of power confinement in the inner dielectric and the outer liquid crystal clads—the feature that can be effectively used in several integrated optic applications including electromagnetic field coupling and optical sensing. It has been found that, with the increase in propagation constants, the confinement patterns become almost uniform, which is not observed corresponding to LCOFs of higher core diameters. Higher confinement of power in the cladding regions can be interpreted as if the power is leaking off the LCOF core, and this feature can be attributed to the presence of radially anisotropic liquid crystal material in the outermost section of the guide.Furthermore, the fundamental investigations of HCLCOFs reveal that the introduction of conducting helical windings at the interface of the core-inner dielectric clad sections makes the guide to attain the attenuation property. The possibility of suitable adjustments of the helix pitch leaves the option to alter the attenuation characteristics of the guide, the feature quite indispensable in several optical applications.



---

*Source: 504868-2013-03-14.xml*