# Flux Density through Guides with Microstructured Twisted Clad DB Medium

Authors: M. A. Baqir; P. K. Choudhury
Journal: Journal of Nanomaterials
 (2014)

---

## Abstract

The paper deals with the study of flux density through a newly proposed twisted clad guide containing DB medium. The inner core and the outer clad sections are usual dielectrics, and the introduced twisted windings at the core-clad interface are treated under DB boundary conditions. The pitch angle of twist is supposed to greatly contribute towards the control over the dispersion characteristics of the guide. The eigenvalue equation for the guiding structure is deduced, and the analytical investigations are made to explore the propagation patterns of flux densities corresponding to the sustained low-order hybrid modes under the situation of varying pitch angles. The emphasis has been put on the effects due to the DB twisted pitch on the propagation of energy flux density through the guide.


---

## Body

## 1. Introduction

Metamaterials are artificially designed mediums that owe unusual phenomena, such as reversal of Snell’s law and Doppler’s effect, negative reflection/refraction, and many others. These exotic features are basically due to the structures of materials rather than their compositions, and are used in many interesting applications, for example, cloaking, perfect lensing, and power confinement [1–4]. References [5, 6] describe negative index chiral metamaterial based on eight-crank molecule designed structure; homogeneous as well as isotropic behavior of chiral medium, based on the periodic inclusion of cranks, have been demonstrated.During the last couple of decades, complex structured guides have attracted the R&D community, primarily because of their varieties of potential applications that include optical sensing, integrated optics, and microwaves devices. The electromagnetic behavior of guides can be tailored by altering structural geometry, medium composition, and the nature of excited electromagnetic field. References [7–11] describe several forms of waveguide structures and their electromagnetic response. Reference [8] demonstrates wave propagation through chiral nihility metamaterial, a special class of chiral medium, in which the real part of permittivity and permeability simultaneously becomes zero, while its chirality remains nonzero. In the context of complex optical micro- and/or nanostructures, helical clad optical guides offer control over the dispersion behavior through suitable adjustments of the helix pitch angle [12–16]. Such helical forms or the twists can be written on fiber structures by the use of current advancements in nanotechnology.The usefulness of DB boundary conditions in the construction of spherical and cylindrical cloaks has been reported in the literature [3]. Within the context, the DB boundary conditions were introduced by Lindell et al. [17–19], in which the normal components of both electric and magnetic fields vanish at the interface [20–22]. DB medium can be easily realized in the case of anisotropic medium if the normal component of permittivity and permeability simultaneously becomes zero; that is,
(1)n^·B→=0,n^·D→=0.
In the case of isotropic medium, the boundary conditions can be written as
(2)n^·H→=0,n^·E→=0.In the present paper, we investigate the propagation through twisted clad guiding medium, wherein helical windings are introduced as DB medium at the dielectric core-clad interface. To the best of our knowledge, such guiding structures are not yet discussed in the literature, and the present work would be the first investigation of the fundamental electromagnetic behavior of microstructured guides with DB medium twisted clad. Rigorous analyses are made of the associated boundary-value problem implementing Maxwell’s equations, and the dispersion relations of the guide are deduced followed by the estimation of supported flux density patterns.


## 2. Theoretical Approach

We consider a cylindrical twisted clad guiding medium, as shown in Figure1, the core/clad sections of which are linear, isotropic, homogeneous, and lossless dielectrics. A DB based microstructured sheath helix is introduced at the core-clad interface, treated under the suitable boundary conditions. Although the sheath is tightly wound, each of the turns is insulated from the neighboring ones. In Figure 1, Ψ represents the angle that the conducting helical windings make with respect to the normal drawn on to the fiber surface in the longitudinal direction.Figure 1
Longitudinal view of the guiding structure.We assumencore and nclad as the refractive indices (RIs) of core and clad regions of the guide, respectively, and a as the core radius, the clad section being infinitely extended. The time t-harmonic and axis z-harmonic wave propagates down the guiding structure, considering the guide to be elongated along the z-direction. Such an assumption leads to the transverse electric/magnetic field components to be of the forms as(3a)Er1=-ju2[βuAJv′(ur)+jωμvrBJv(ur)]ej(ωt-βz+vφ),(3b)Hr1=ju2[jωμvrAJv(ur)-βuBJv′(ur)]ej(ωt-βz+vφ),(4a)
E
φ
1
=
-
j
u
2
[
j
β
v
r
A
J
v
(
u
r
)
-
ω
μ
u
B
J
v
′
(
u
r
)
]
e
j
(
ω
t
-
β
z
+
v
φ
)
,
(4b)
H
φ
1
=
-
j
u
2
[
ω
μ
u
A
J
v
′
(
u
r
)
+
j
β
v
r
B
J
v
(
u
r
)
]
e
j
(
ω
t
-
β
z
+
v
φ
)
,whereinu=(ncorek0)2-β2. Similarly, the transverse components of fields in the clad section can be written as(5a)Er2=-jw2[βwCKv′(wr)+jωμvwDKv(wr)]ej(ωt-βz+vφ),(5b)Hr2=jw2[jωμvrCKv(wr)-βwDKv′(wr)]ej(ωt-βz+vφ),(6a)
E
φ
2
=
-
j
w
2
[
j
β
v
r
C
K
v
(
w
r
)
-
ω
μ
w
D
K
v
′
(
w
r
)
]
e
j
(
ω
t
-
β
z
+
v
φ
)
,
(6b)
H
φ
2
=
-
j
w
2
[
ω
μ
w
C
K
v
′
(
w
r
)
+
j
β
v
r
D
K
v
(
w
r
)
]
e
j
(
ω
t
-
β
z
+
v
φ
)
,withw=β2-(ncladk0)2. In the above set of equations, μ and ε are, respectively, the permeability and the permittivity of the medium, ω is the angular frequency in the unbounded medium, β is the propagation constant in the medium, k0 is the free-space propagation constant, v is the azimuthal mode index, and primes represent differentiation with respect to the arguments of field functions. Also A, B, C, and D are unknown constants, to be estimated by the use of boundary conditions.Now, considering the existence of DB twisted clad, boundary conditions will essentially include vanishing normal components of electric/magnetic fields and continuous tangential components of the same at the medium interface. Implementation of this ultimately leads to the evolution of a set of four homogeneous equations, followed by the dispersion relation (as obtained by collecting the coefficients of unknown constants) given as(7)f(β)=14ω2μεw2{Kv-1(aw)+Kv+1(aw)}2sin2ψ×[v2ω2μεa2u4Jv2(au)(cosψ-sinψ)2×-β24{Jv-1(au)-Jv+1(au)}2(sinψ-cosψ)2]+{Kv(aw)cosψ+vβaw2Jv(au)sinψ}×{Kv(aw)cosψ-vβaw2Kv(aw)cosψ}×{v2μεω2a2u4Jv2(au)-14β2(Jv-1(au)-Jv+1(au))2}×(sinψ-cosψ)2=0.The unknown coefficients could be determined solving the aforesaid four equations obtained after implementing the boundary conditions. After some rigorous mathematical steps, the coefficientsA, B, and C can be determined in terms of the coefficient D as follows:
(8)A=-ju2aJv′(ua)Θ(u2aγΔ-Ξξ)+w2a(Θξ-ωεKv(wa))ωεaw2γ(Θξ-ωεKv(wa)δ)D,B=uJv(ua)Θ(u2aγΔ-Ξξ)+w2a(Θξ-ωεKv(wa)δ)w2aγ(Θξ-ωεKv(aw)δ)D,C=ju2aγΔ-ω2μεuw3Kv′(wa)ξωεKv(aw)(w2acosψ-βvsinψ)ξ-δD,
where
(9)γ=ω2μεu2Jv(ua)Jv′(ua)sinψ-βJv′(ua)(u2aJv(ua)cosψ+βvJv(ua)sinψ),ξ=Jv2(ua)(u2acosψ+βvsinψ)-βu2a2Jv′(ua)sinψ,Θ=ωϵKv(wa)(w2acosψ-βvsinψ),Δ=w2Kv(wa)cosψ+βvKv(wa)sinψ,Ξ=ω2μεw3Kv(wa)sinψ,δ=ωεuwa2Kv(wa)sinψ.Now, the flux density in the guide along the direction of propagation can be estimated, and, following the procedure in [11], it will finally assume the forms corresponding to the core/clad sections, as follows:
(10)(Sz)Core=β2u4(ε|A|2+μ|B|2){(vr)2Jv2(ur)+(u2)2σ},(Sz)Clad=β2w4(ε|C|2+μ|D|2){(vr)2Kv2(wr)+(w2)2Λ},
with
(11)σ=Jv-12(ur)+Jv+12(ur)-2Jv-1(ur)Jv+1(ur),Λ=Kv-12(wr)+Kv+12(wr)-2Kv-1(wr)Kv+1(wr).
Equations (10) can be written in terms of only one unknown constant, for example, in terms of D, as determined by the use of (8).


## 3. Results and Discussion

We now investigate the propagation features of flux density through DB twisted clad microstructured guide. However, for this purpose, it remains essential to evaluate the allowed values of propagation constants corresponding to different sustained modes under the conditions when the DB twists are typically oriented perpendicular as well as parallel to the direction of wave propagation. For the computational purpose, we consider the RI values of the core and the clad sections as 1.485 and 1.470, respectively, and the operating wavelength is taken to be 1.55μm. The chosen parameters are, indeed, the widely accepted ones in the case of modeling and simulations related to the electromagnetic properties of optical fibers. Also, we consider two typical values of core radius as 10 μm and 20 μm, the clad radius being fixed at 100 μm. The behavior of normalized flux density through core/clad regions is studied considering three low-order hybrid EH01, EH11, and EH21 modes having propagation constants, as obtained by the use of eigenvalue (7), as shown in Table 1.Table 1
Hybrid mode propagation constants.
Pitch angleψ
Core radiusa
Mode
Propagation constantβ (m−1)
0°
10 µm
EH01
5.970 × 106
EH11
5.975 × 106
EH21
5.968 × 106
20 µm
EH01
5.960 × 106
EH11
5.965 × 106
EH21
5.970 × 106
90°
10 µm
EH01
5.850 × 106
EH11
5.890 × 106
EH21
5.870 × 106
20 µm
EH01
5.950 × 106
EH11
5.980 × 106
EH21
5.965 × 106Figures2(a) and 2(b), respectively, demonstrate the normalized flux density patterns through core and clad regions of guide having core radius 10 μm under the situation when the DB twist is oriented perpendicular to the direction of propagation; that is, Ψ=0°. In all the figures, the results corresponding to the EH01, EH11, and EH21 modes are represented by black dash-dot, solid red, and dashed blue lines. We observe from Figure 2(a) that the flux generally exhibits an oscillatory trend in the core region of the guiding medium. Also, in the central region of the core, EH11 hybrid mode has higher energy than the other two EH01 and EH21 modes. Further, flux densities due to the EH01 and EH11 modes decrease with the increase in radial distance, whereas that due to the EH21 mode exhibits a minor increase only.Flux density through the core (a) and clad (b) corresponding to 10μm core radius and 0° helix pitch.
(a)
(b)In the clad section, as Figure2(b) shows, the aforesaid three hybrid modes exhibit similar flux characteristics, and it remains extremely small as compared to that in the core region. At the core-clad interface, it shows relatively large value, and then step-falls with a little increase in clad radius, and become almost vanishing as the radius is further increased.Figures3(a) and 3(b) correspond to the situations when the DB helical windings are put parallel to the optical axis, that is, the helix pitch Ψ=90°. We observe in this case that the confinements of flux in the core section (Figure 3(a)) exhibit almost similar patterns, as observed in the case when DB twists are oriented in the direction perpendicular to the wave propagation (Figure 2(a)), but its values are drastically reduced corresponding to all the sustained modes. Otherwise, the EH01 and EH11 modes correspond to higher power than the EH21 mode in the central region of core section. In the clad section too, the flux is further reduced (Figure 3(b)) as compared to the case of Figure 2(b), exhibiting similar density patterns, as observed before, for the situation of Ψ=0°. Thus, comparing the results of Figures 2 and 3, we notice the effect of introducing DB sheath helix structure at the dielectric core-clad interface. We find that making the angle of pitch (of DB helix) as 90° greatly reduces the propagation of power, and therefore, a change in pitch would make the guide to amplify or attenuate the signal depending upon lowering or increasing its value. This clearly demonstrates the importance of such guiding structures.Flux density through the core (a) and clad (b) corresponding to 10μm core radius and 90° helix pitch.
(a)
(b)Figures4 and 5 correspond to the situation when the core radius is doubled, that is, increased to 20 μm. We observe that the trends of flux density patterns are almost similar to the situation as observed before for 10 μm, but the oscillations of flux are increased in this case. This is very much obvious due to more room for the fields to be sustained, thereby causing dispersion and so forth. However, in this as well, we observe that the increase of pitch of DB twist essentially reduces the flow of flux density, as evidenced by Figures 4 and 5.Flux density through the core (a) and clad (b) corresponding to 20μm core radius and 0° helix pitch.
(a)
(b)Flux density through the core (a) and clad (b) corresponding to 20μm core radius and 90° helix pitch.
(a)
(b)


## 4. Conclusion

The eigenvalue equation is deduced for the aforesaid structure, and the propagation features of the low-order hybrid modes are estimated. Investigations are made of the energy flux density patterns corresponding to the sustained modes, and the effects are observed due to the alterations of pitch angles of the microstructured DB twists introduced at the dielectric core-clad interface. The flux patterns reveal profound effect of helix pitch, and making the DB windings parallel to the axis of guide greatly reduces power confinements in the guiding region. As such, these guides may be used for the purpose of amplifying and/or attenuating signals in optical systems.



---

*Source: 629651-2014-01-02.xml*