The Dispersion Characteristics of a One Dimensional Plasma Photonic Crystal Having Inhomogeneous Plasma Density Profiles

The dispersion characteristics of binary 1D-PPCs having inhomogeneous plasma in the unit cell are studied. Using the transfer matrix method the required dispersion relations are obtained. Here the linear and exponential plasma density profiles are considered and compared with the homogeneous plasma having uniform density profile. It is observed that the inhomogeneity in plasma layer highly affect dispersion curves. By comparing the dispersion curves obtained in all considered cases, it is found that the widths of band gaps and phase velocities are always larger for exponential density profile than the linear uniform density profiles in the considered frequency range.


Introduction
Photonic crystals (PCs) are periodically structured electromagnetic media in which certain range of electromagnetic (EM) waves are not allowed to propagate through the structure.This range of frequencies is called photonic band gaps (PBGs).The periodicity of the structure and the periodic variation of dielectric constant of different materials are the essential parameters for the formation of these PBGs [1].These PCs have achieved archival interest worldwide because for its numerous applications such as: Filters, optical switches, waveguides, cavities and design of more efficient layers, etc. [2].
Recently, plasma photonic crystals (PPCs) have attracted the attention of researchers because the properties of these PPCs are externally controllable and it possesses the characteristics of conventional PCs and plasma.The unit cell of PPCs consists of periodic arrangement of plasma and dielectric/air.The dispersion relation of binary one dimensional plasma photonic crystals 1D-PPCs by solving Maxwell equations using a method analogous to the Kronig-Penny's problem in quantum mechanics is derived by Hojo et al. [3] and discussed the effect of plasma density, thickness of plasma layer and dielectric constant of dielectric media on the dispersion characteristics.Subsequently many research groups analyzed the dispersion and optical properties of 1D-PPCs [3][4][5][6][7][8].The dispersion relation of binary 1D-PPCs considering the collision effect in plasma with obliquely incident EM waves is deduced by Guo Bin [6].Also the dispersion and propagation characteristics of ternary 1D-PPCs have been studied by Prasad et al. [7,8] and they found that the ternary 1D-PPCs provide additional degree of freedom to control the dispersion and propagation characteristics compared to binary 1D-PPCs.But most of the researchers have analyzed the dispersion characteristics, modal propagation characteristics, reflection and transmission coefficients of 1D-PPCs in which a unit cell consists of homogeneous plasma and homogeneous dielectric materials and they obtained that these properties are being controlled by plasma density, thickness of plasma layers, collisions in plasma layer, dielectric constant of dielectric materials and external magnetic fields [9,10].The PBGs of 1D-PPCs can also be improvise in desire range by considering inhomogeneous dielectric layer in the unit cell of PPCs [11,12].
Similarly, the dispersion characteristics and optical properties of PPCs can also be controlled by considering inhomogeneous plasma in the unit cell.This consideration is more practical also because homogeneous plasma having uniform density is rarely realized in the laboratory plasma [13].The inhomogeneities in the plasma layer may occur due to spatially dependent plasma density.Therefore, in present communication the propagation of EM waves in 1D-PPCs having inhomogeneous plasma in the unit cell is considered.By using the transfer matrix method (TMM) the dispersion relations of proposed structures are obtained.The dispersion characteristics of PPCs are compared by taking exponential density profile and linear density profile with uniform (constant) density in the plasma layer.The inhomogeneous plasma with linear density profile [14] and exponentially density profile [15] have been chosen because their closeed form solutions are available in terms of special mathematical functions such as Airy and Bessel functions.The paper is organized as follow: In Section 2 formulas for the dispersion relation of the proposed structure is given.The other necessary formulas used in this paper are also presented.Section 3 is devoted to result and discussion.A conclusion is drawn in Section 4.

Theoretical Modeling
A plane EM waves with angular frequency ω is assumed to obliquely incident on the 1D-PPCs structure.The unit cells of considered 1D-PPCs have periodic arrangement of inhomogeneous plasma and homogeneous dielectric material.The unit cell is shown in Figure 1.The dielectric material and inhomogeneous plasma are assumed to be non-magnetic.For simplicity, collisions in inhomogeneous plasma are neglected.The plasma density varies either linearly or exponentially in space and is given by n(x).The linear plasma density profile is given as: Similarly, the exponential plasma density profile is given as: where b is the width of plasma layer, p is gradation parameter for controlling variation of density in the plasma layer and n cr is the critical density [14] 2 2 0 cr n m e    .The permittivity profile for inhomogeneous plasma and dielectric media is given by: For linear plasma density profile, the permittivity is written as: ; 1 Here ε 1 is the dielectric constant of dielectric layer, , a and b are widths of dielectric and plasma layer respectively.The permittivity profile is linearly varying with space along x-direction.The electric field in the case of TE mode is in y-z plane.Along the z-direction there is no change in permittivity, so z-component of wave vector is conserved.The one dimensional wave equation for the spatial part of Typical field solution can be expressed as and using this in Equation ( 3) we can write the above equation both regions: inhomogeneous plasma layer and dielectric medium in the n th unit cell as: where β is the z-component of wave-vector and is given by   Here, it is assumed that EM waves incident from vacuum at certain angle θ.The solutions of above equations for electric fields in the n th unit cell are given by Ai and Bi are Airy functions.Now, using the continuity ition of electric field E(x) and its derivatives at boundaries and cond   and TMM, we obt atrix ents o plex a cell to corresponding amplitudes in nth cell.to rem, wave propagating in a e where m 11 , m , m and m 22 are elem f unit translation matrix [16], which relates com mplitudes of incident and reflected wave in   The constant K is known as the Bloch wave number.Here scr k is dropped for simplicity in notations and k 1x is K. Now the problem and is to determine K and ng Equation ( 6), odic condition (8) for Bloch wave is simply given: It follows from Equations ( 7) and (8) that the column vector of Bloch wave satisfies the following eigen-value problem 11 12 The phase factor e iKΛ is the eigen-value of the unit translation matrix (m 11 , m 12 , m 21 and m 22 ) and is given by: So using Equation (11), the dispersion relation for proposed structure can be written as: ements for exponential density profile and linear density are given in the Appendix 1 and Appendix 2, respectively.In case of uniform (constant) density, the unit translation element ily obtained [6].

3.
dth a.Here the thickness of inhomoogeneous dielectric (12) Similarly, we can also find the unit translation elements: m 11 , m 12 , m 21 and m 22 for exponential density profiles.The elements of unit translation el can be eas-

Results and Discussion
The Equation ( 12) is known as dispersion relation which will give all the information regarding PBGs of the considered (linear plasma density and exponential plasma density) profiles.The unit cell of proposed PPCs structure contains inhomogeneous plasma of width b and homogeneous dielectric of wi geneous plasma b is related with hom layer width a by the relation Here d is the ratio of thickness of two layers.The permittivity profile in inhomogeneous plasma layer is and homogeneous dielectric layer is glass with dielectric constant  1 = 2.25.In the case of homogeneous plasma layer, the plasma frequency ω p is constant and its value is  .Also, the width of band gap is large for exponential plasma density profile than the linear plasma density profile.By comparing these results with the case of homogeneous plasma having uniform density profile (dotted line), it is observed th ap of unifo profile and linear density profile is approximately equal dence in at the width of band g rm density for first band gap but differ in higher order band gaps.Figure 2 show that as the angle of inci creases from π 10 , the lower edge frequencies of first and second band gaps of exponential density profile (solid line) remain intact with linear density profile (dash line) and are shifted towards higher frequency.Such types of behaviors are not observed in the case of uniform density profile (dotted line).This is a remarkable result of present investigation.It is also observed that with increase in angle of incidence, the width of band gaps observed that as the thickness of plasma layer changes u ed nd gaps and phase velocities wards lower frequency.At d = 2.0 (shown i are always larger for exponential density profile than corresponding band gaps and phase velocities of the linear the lower edge frequency of first band gap become different for all profiles.From the analysis of Figure 3, it is ered frequency range.The reason is that in the case of exponential density profile, the permittivity contrast changes rapidly than the linear density profile and uniform density, hence produces more Bragg's reflections. Figure 3 are plotted to analyze the effect of thickness of inhomogeneous plasma layer on the dispersion characteristics.Figure 3(a) shows that the lower edge frequency of first and higher order band gaps for exponential density profile (solid line) are different from linear densit e (dash line) at d = 0.5.The interesting feature of this graph is that the lower edge frequency, of first band gap, of linear density profile is now intact with uniform density profile (dotted line).As the thickness of plasma layer increases from d = 0.5 to 1.0, as shown in Figures 3 ency of first band gap for linear density profile and from d = 0.5 to d = 2.0, number of band gaps get increased.Also, it is found that widths of first and higher order band gaps are different for all three density profiles and larger for exponential density profile., ε 1 = 2.25.In this case, the considered volume average permittivity is 0.8394.By comparing Figure 4 from the Figure 2(a) which is plotted at the same parameter with volume average permittivity 0.5.It is observed that as the volume average permittivity increases the widths of band gaps also increase for all profiles and the lower edge frequency of first band gap of all considered profiles remains intact but shifted to wards lower frequency.
From this study, it is found that position and width of The onic crysd because homogeneous plasma having uniform density is rarely realized in the laboratory plasma.A comparison is also made with 1D-PPCs having homogeneous plasma with uniform density profile.It is observed that, for a constant volume average permittivity, the exponential density profile will always give larger band gaps widths than linear density profile and uniform density profile.With increase in angle of incidence, the widths of band gaps increase and at the same time band gaps are shifted towards higher frequency for all considered profiles.The number of band gaps also increases with increase of the thickness of inhomogeneous/homogeneous plasma layer.This analysis shows hat with increase in the volume average permittivity, the w ths dge frequency of first band gap sity profile and linear density profile dispersion characteristic of a plasma phot s tal having inhomogeneous plasma density profile have been studied at the first time in our knowledge.The exponential and linear types of inhomogeneity in the plasma density profile are considere t id of band gaps are increased and the lower edge frequency of first band gap remains intact but shifted towards lower frequency for all profiles.The important finding of this study is that the lower e for exponential den remains intact and shifted towards higher frequency with increase in incident angle.
Finally, it is concluded that the width and position of band gaps strongly depend on the type of density profile of the plasma layer.Therefore, for correct prediction about the width and position of band gaps, the type of density profile should be ensured in PPCs.

Acknowledgements
The authors are grateful to Dr. B. Prasad and Dr. R. D. S. Yadava for their continuous encouragement and supports in many ways.band gaps strongly depend on the density profile and by choosing suitable profile; one can tune the band gaps in desired frequency range.qJ q qt J q q p J q pt J q J t t J t pt J t e t pt p J q qt p J q m pt J q J q qt J q J q qt J q J q qJ q pt J q qt J q pt J q J q qt J q J q qt J q J q t pt p J q qt p J q m e pt J q J q qt J q J q qt J q J q Y t t Y t pt Y t p q J qp t

Appendix: 1 Unit Translation Matrix Elements for 1D-PPCs with Exponential Plasma Density
x p p ik a p p p p p p t J q qt J q e pt J q J q qt J q J q qt J q J q qJ q qt J q q p J q pt J q J t t J t pt J t e t pt p J q qt p J q m pt J q J q qt J q J q qt J q J q Y qJ q pt J q qt J q pt J q J q qt J q J q qt J q J q qt J q qJ q q p J q pt J q J t t J t pt J t t pt p J q qt p J q m e pt J q J q qt J q J q qt J q J q Y t t Y t pt Y t p q J qp t x p p ik a p p p p p p t J q qt J q e pt J q J q qt J q J q qt J q J q    

Figure 1 .
Figure 1.Schematic representation of the unit cell of binary 1D-PPCs having inhomogeneous plasma density layers with homogeneous dielectric layers.
functions of  and .In terms of column vec esentation and usi the peri tor repr incident angles θ, gradation parameter p, a, d and  1 involved in the numerical calculations.

Figure 2 .
Figure 2 are plotted to observe the effect o

Figure 4
is plotted to analyze the effect volume average permittivity on the dispersion characteristics for a = 5000 µm, d = 0.1, π 10 