Open access peer-reviewed chapter

Propagation of Electromagnetic Waves in Anisotropic Photonic Structures

By V.I. Fesenko, I.A. Sukhoivanov, S.N. Shul’ga and J.A. Andrade Lucio

Submitted: May 31st 2012Reviewed: November 5th 2012Published: February 13th 2013

DOI: 10.5772/54847

Downloaded: 1906

1. Introduction

In this chapter we discuss plane-wave propagation in a layered arbitrarily anisotropic media.One-dimensional (1D) arbitrary layered structure is usually formed by stacking together layers of several different materials with some specific thickness djand refractive index njas depicted on Fig.1.

Figure 1.

Arbitrary layered structure

Nowadays layered photonic structures (LPSs) are key of optoelectronic and microwave devices such as, phased-array antennas, microcavities and mirrors [1,2], filters of xWDM systems [3], waveguide structures, photodetectors, sensors and others. In case of active devices the layered structures are usually used in form of superlattices [4-6], multiple quantum wells [7,8] and asymmetric multiple quantum wells [9].

Different kinds of materials are used today for design LPSs, such as linear and nonlinear [3,10] dielectric materials, anisotropic or bi-anisotropic materials [11], chiral media [12], metamaterials [13], etc. If one use these materials one can effectively control emission, propagation and detection of the electromagnetic waves, and develop new designs of the photonic devices and those one for other parts of the electromagnetic spectrum. For instance, very active research worldwide is concentrated currently on the THz range, attempting to overcome the so called problem of the terahertz gap [14].

The analysis of propagation of the electromagnetic waves in periodic [3,11,13,15], quasiperiodic [12,16,17] and random [18] layered media is a problem which extends over the all fields in the modern physics. Optics is the area where it is crucial to calculate spectral characteristics, absorbance coefficients, polarization properties and other features of the multilayer structures in a wide spectral range and at various thicknesses or material properties of constituents. Even in the fiber optics, where usually propagation characteristics of the optical pulses [19, 20] are considered, the spectral characteristics become of the principal interest when optical channel incorporates such inhomogeneities as fiber Bragg gratings and fiber knots, and microresonators based on them.

In addition the use of the optical control techniques for phased-array antennas [21-22] promises to alleviate many of the problems associated with traditional electronic steering systems. The unique properties of layered anisotropic photonic structures (for example see [11]) are suitable for this application.

Here we discuss the optical properties of one-dimensional arbitrarily anisotropic photonic multilayers. The main objective of the chapter is the obtaining of a solution to the numerical problem of the electromagnetic plane wave interaction with arbitrarily anisotropic and arbitrarily inhomogeneous one-dimensional photonic structures. It is well known that many of novel technological designs have resulted from analysis of the properties of materials and creation of new structural configurations for them. In order to develop a new structural configuration with unique properties, one needs to thoroughly understand the characteristics of the structure. This can be accomplished by applying an advanced computational engine.

Today several numerical techniques are commonly used to compute the spectral characteristics of the layered photonic structures and electromagnetic field distribution in the interiors [3,12-15,23-31]. The most known and, probably, most usable are the finite element method (FEM) [3,24], the transfer matrix method (TMM) [25,26], the finite difference time-domain method (FDTD) [15,27] and the beam propagation method (BPM) [28]. Unfortunately, some methods are entirely disregarding the anisotropy and the inhomogeneity of the constituent materials. The most of methods which do account for material anisotropy require that the permittivity tensor be diagonal. Others allow for nondiagonal tensors, but require that the off-diagonal elements be small in comparison to the diagonal terms [28]. Although these techniques are adequate for many layered structures, they cannot be easily applied to multilayers in which the anisotropy is arbitrary oriented along an oblique axis on the random layers. In contrast to these numerical methods the method discussed here makes it possible the analyzing of inhomogeneousone-dimensional anisotropic multilayers with an arbitrary permittivity tensor and the optical axis arbitrarily oriented on any layer of LPS.

A general theory of electromagnetic propagation in periodic anisotropic layered media has been treated by a number of authors [26, 29-31]. The present chapter describes the efficient physico-mathematical model pertinent to one-dimensional optical-range microstructures based upon anisotropic materials. The electromagnetic field scalarization procedure [29-30] is used after the initial vector electromagnetic diffraction problem is reduced to the boundary problem for two scalar potentials. As a result, a set of linear algebraic equations are obtained.By solving them we find the unknown transmission and reflection spectra for the structure under study. The major advantage of the proposed method is that homogeneous, piecewise homogeneous and continuously inhomogeneous flat-layered anisotropic media can be analyzed on the same footing.

The chapter is organized as follows: in Section 2 we present the method of calculation employed in the chapter, which is based on electromagnetic field scalarization procedure in conjunction with the finite-difference method; Section 3 is devoted to the presentation of the numerical results, together with the discussion of their main features; we summarize our study and conclude the chapter in Section 4.

2. Theoretical model

In this section we present the mathematical background for calculation spectral characteristics of anisotropic layered media. The presented theory is applicable to any anisotropic layers with arbitrary orientation of the optical axes on each layer and for arbitrary angle of incidence.

The structure under consideration is schematically depicted in Fig. 2. Let’s introduce the Cartesian coordinate system x, y, z such that the z axis is directed vertically upward. In this coordinate system, an inhomogeneous anisotropic layered medium is represented by the single layer that occupies the domain b<z<0, <x,y<+. The upper free half-space z>0and substrate z<bare homogeneous and isotropic and have permittivities ε0,μ0and εc,μc, respectively. In general case, all layers in this geometry are lossy media.

Figure 2.

The benchmark photonic structure

The layered character of the media consists in next conditions or combination of them:

  • permeability μ^and permittivity ε^are continuous functions of the variable z. These complex-valued tensors η^=ε^,μ^can be expressed in Cartesian coordinates as:

η^=[ηxxηxyηxyηyxηyyηyzηzxηzyηzz];E1
  • the media are piecewise continuous, i.e. there are boundary surfaces z=zj=const,(const>0), where properties of anisotropic media are varying stepwise;

  • in the points z=0and z=b, (b>0) the medium is bounded by homogeneous conducting planes or the planes permeable for the electromagnetic field. In our case the layered medium is confined by the impedance planes. These planes are characterized by impedance dyads L^(a)and L^(u):

L^(a,u)=[L^xx(a,u)L^xy(a,u)L^yx(a,u)L^yy(a,u)].E2

Here indices a (above) and u (under) correspond to upper (z=0) and lower (z=b) boundaries, respectively.

We will start with Maxwell’s equations, for complex field vectors E(R)and H(R), in the form:

×E(R)ik0μ^(z)H(R)=(4π/c)M(R),×H(R)+ik0ε^(z)E(R)=(4π/c)J(R),E3

where k0 is the wave number in free space; c is the velocity of light; J(R)and M(R)are electric and magnetic volume current densities, respectively; R=(x,y,z)is the radius vector. Here we assume the harmonic time dependence exp(iωt)of the fields. Equations (3) are satisfied everywhere within the medium except the interfaces.

The electric and magnetic fields must satisfy the suitable boundary condition at the interfaces of the layered media:

  • at the first, the tangential components of the electromagnetic field must be continuous at the all boundaries of the layered media:

{z0×E}=0,{z0×H}=0,(z=zj,j=1,2,...,N).E4

where z0is unit vector along zaxis; EE(R),HH(R)are electric and magnetic field components that are orthogonal to z0. Throughout the chapter we use braces {}for next operator designation {f(z)}f(z+0)f(z0);

  • at the second, we introduce impedance boundary conditions that are desired on the above and the bottom boundaries of the inhomogeneous anisotropic structure:

E+L^(a)z0×H=0,(z=0);EL^(u)z0×H=0,(z=b).E5

The next one, in what follows we assumed that external sources and electromagnetic field components are represented by spatial harmonics with wave vector κ=(κx,κy,0)

J(R)J(κ,z)exp(iκr),M(R)M(κ,z)exp(iκr);E6
E(R)E(κ,z)exp(iκr),H(R)H(κ,z)exp(iκr).E7

In expressions (6) – (7) κx,yare the arbitrary complex constants, J(κ,z),M(κ,z)and E(κ,z),H(κ,z)are the vector amplitudes of the sources and the fields, respectively.

Now let’s try to obtain general expressions for the transmittance and reflectance of a layered medium.

2.1. Reduction of the electromagnetic field diffraction problem to a boundary value problem for scalar potentials

Now let’s consider in details the solving of the initial electromagnetic field diffraction problem. At first we should to introduce right-hand basis of vectors az,al,at:

az=z0,al=n,at=z0×n.E8

In writing eq (8) we have assumed that: n=κ/κis the unit vector, that is situated in the plane z=0; κ=κx2+κy2is the branch of the square root, which is chosen such that condition 0arg()πshall be satisfied. These unit vectors obey the orthogonality relations:

aσaτ=0,(στ).E9

Then after scalar multiplication of the Maxwell’s equations (3) by unit vectors az,aland using simple conversion of vector algebra we obtain Ez,l=az,lE(k,z),Hz,l=az,lH(k,z)in terms of scalar potentials:

e(κ,z)=atE(κ,z),h(κ,z)=atH(κ,z).E10

In the basis of orthogonal vectors (8), vector amplitudes E(κ,z),H(κ,z)can be written down as:

E(κ,z)=Vε(n,z)e(κ,z)Wε(κ)h(κ,z)+(4πi/k0c)α^ε(n,z)J(κ,z),H(κ,z)=Vμ(n,z)h(κ,z)Wμ(κ)e(κ,z)+(4πi/k0c)α^μ(n,z)M(κ,z).E11

The following symbols are used in (11): Vη(n,z)and α^η(n,z), (η=ε,μ)are the vector functions and the dyad functions, respectively; Wη(κ)– vector differential operator:

Vη(n,z)=z0×n+1/aη(n,z)[bη(n,z)n+cη(n,z)z0],E12
aε(n,z)a^ε(n,z)at×ε(z)×at=(εlznεllz0)z0+(εzlz0εzzn)n,E13
Wε(κ)=(1/aη(n,z)k0)×[(εzznεzlz0)iz+(εzlz0εzzn)k0],E14
aεaμ,a^εa^μ,W^εW^μ,(εμ).E15

In expressions (12) – (15), scalars εστεστ(n,z),μστμστ(n,z),(σ,τ=z,l,t)are components of the relevant dyads ε,μ(1):

εστ(n,z)=aσε(z)aτ,μστ(n,z)=aσμ(z)aτ.E16

The scalar functions aη(n,z),bη(n,z),cη(n,z),(η=ε,μ)that are used in (12) – (14) have the form:

aε=εzzεllεzlεlz,bε=εztεlzεzzεlt,cε=εltεzlεllεzt,E17

aεaμ,bεbμ,cεcμ,(εμ).

On the next step after scalar multiplication of the Maxwell’s equations (3) by the unit vector at, and using the expressions (11) we obtain:

Dss(κ)h(κ,z)+Dsp(κ)e(κ,z)=(4π/c)qs(κ,z),Dpp(κ)e(κ,z)+Dps(κ)h(κ,z)=(4π/c)qp(κ,z).E18

Expressions (18) represent the system of coupled ordinary differential equations for two scalar potentials e(κ,z)and h(κ,z)within the interval -b<z<0; the external sources are entered into the quantities qv(κ,z); Dvξ(κ)are the scalar operators that depend on κ(v,ξ=s,p). These operators in explicit form are written as follows:

qs(κ,z)=ik0[z0×n+1aμ(n,z)(eμz0+dμn)]M(κ,z)++[(zεlz+iκεll)z0(zεzz+iκεzl)n]J(κ,z)aε(n,z),qsqp(εμ,JM,MJ);E19
Dss(κ)=zεzzaε(n,z)z+iκ[zεlzaε(n,z)+εzlaε(n,z)z]+k02|μ^(z)|aμ(n,z)κ2εllaε(n,z);E20
ik01Dsp(κ)=iκ[cε(n,z)aε(n,z)+eμ(n,z)aμ(n,z)]zbε(n,z)aε(n,z)zeμ(n,z)aμ(n,z);E21
DssDpp,DspDps,(εμ)E22

where:

dε=εzzεtlεzlεtz,eε=εllεtzεlzεtl,dεdμ,eεeμ,(εμ)E23

After the substituting expressions (11) into condition (4) we obtain boundary conditions for scalar potentials e(κ,z), h(κ,z)that satisfied on all boundaries z=zj:

{e(κ,z)}=0,{h(κ,z)}=0,{1/aμ(n,z)×[(μzzziκμlz)e(κ,z)ik0bμ(n,z)h(κ,z)]}=0,{1/aε(n,z)×[(εzzziκεlz)h(κ,z)ik0bε(n,z)e(κ,z)]}=0.E24

Substituting expressions (11), for vector amplitudes E(κ,z),H(κ,z), into impedance boundary conditions (5) we obtain two pair of equations for scalar potentials and their derivates. The first pair of equations for lower boundary z=b+0has the form:

[ik0app(u)(κ)+bpp(u)(κ)z]e(κ,z)+ik0aps(u)(κ)h(κ,z)=0,[ik0asp(u)(κ)+bsp(u)(κ)z]e(κ,z)+[ik0ass(u)(κ)+bss(u)(κ)z]h(κ,z)=0,E25

where:

app(u)(κ)=1+κμlz(κ,z)k0aμ(n,z)Ltt(u)(κ),bpp(u)(κ)=μzz(z)aμ(n,z)Ltt(u)(κ),aps(u)(κ)=Ltl(u)(κ)+bμ(n,z)aμ(n,z)Ltt(u)(κ),bss(u)(κ)=εzz(z)aε(n,z),bsp(u)(κ)=μzz(z)aμ(n,z)Llt(u)(κ),asp(u)(κ)=bε(n,z)aε(n,z)+κμlz(κ,z)k0aμ(n,z)Llt(u)(κ),ass(u)(κ)=κεlz(κ,z)k0aε(n,z)bμ(n,z)aμ(n,z)Llt(u)(κ)+Lll(u)(κ).E26

In expressions (26) the equality is used Lστ(u)(κ)=aσL^(u)aτ, (σ,τ=z,l,t).

The boundary conditions for top plane (z=0)are similar to (25). They are resulting from the following replacements in (25) – (26): bpp(u)(κ)bpp(a)(κ),aps(u)(κ)aps(a)(κ),bsp(u)(κ)bsp(a)(κ),bss(u)(κ)bss(a)(κ);in formulas for app(u)(κ),asp(u)(κ),ass(u)(κ)upper index u → a and κκ.

2.2. Numerical solutions by finite-difference method

Now let’s build a numerical solution of the problem of monochromatic plane wave diffraction on the inhomogeneous anisotropic layered structure. We will assume that the structure is piecewise homogeneous along the axis z and within each homogeneous layer the anisotropic material is gyrotropic one, or, particularly, an uniaxial material with arbitrary orientation of the optical axes.

The electromagnetic properties of the benchmark structure in a fixed point of the space are defined by permeability and permittivity dyads:

ε^(z)=ε(z)I^+(ε||(z)ε(z))aaif(z)a×I^,μ^(z)=μ(z)I^+(μ||(z)μ(z))bbig(z)b×I^.E27

Here, ε(z),ε||(z),f(z)and μ(z),μ||(z),g(z)are twice differentiable functions of the variable z; I^is the identity dyad; aand bare the unit vectors in the direction of the optical axes which have the following components in the Cartesian coordinate system:

aa(z)=(cosθasinφa,cosθacosφa,sinθa),bb(z)=(cosθbsinφb,cosθbcosφb,sinθb).E28

For the sake of clarity, the sloping angles θa,θband the azimuthal angles φa,φb, which determine the optical axes direction, are shown in Fig. 3 and therewith we have:

π/2θa,θbπ/2,0φa,φb2π.E29

Figure 3.

The optical axis orientation in the arbitrary layer

At the same time we consider that an incident (s- or p-polarized) plane wave arrives from the free half-space (z>0)in the direction of the unit vector linwhich is determined by the sloping angle θand the azimuthal angle φ, as depicted in Fig.2. Its components in the Cartesian coordinate system are as follows:

lin=(cosθcosφ,cosθsinφ,sinθ),E30
π/2θπ/2,0φ2π.E31

Let us assume the inhomogeneous anisotropic structures under consideration with thickness b are placed on isotropic (or anisotropic) homogeneous substrate with permittivity and permeability εc,μc. In general case, εcand μcare complex values. We assume also that anisotropic layer is bounded above (z>0) by free half-spice with ε0and μ0. In the case presented here, the anisotropic layer is inhomogeneous in the thickness; it means, that structure’s parameters are depended on the coordinate z, but they are invariable along the axes x and y. Generally these parameters are piecewise continuous functions of z.

Then scalar potentials e(κ,z), h(κ,z)beyond the anisotropic layer will be presented by the following expressions:

  • in the free half-spice (0<z<+):

e(κ,z)=(eiγ0z+Rppeiγ0z)Ap+RpsAseiγ0z,h(κ,z)=(eiγ0z+Rsseiγ0z)As+RspApeiγ0z;E32
  • in the substrate (<z<b):

e(κ,z)=[TpsAs+TppAp]eik0nc(z+b),h(κ,z)=[TssAs+TspAp]eik0nc(z+b).E33

Here γ0=k0sinθ;nc=εcμccosθ,(0argncπ); complex values Apand Ascharacterize electromagnetic (s or p –polarized) wave’s components. The s- (p-) polarized light corresponds to an electric (magnetic) field be parallel to the layers.

In the expressions (32) – (33) we have introduced the complex reflection Rvξ(v,ξ=p,s)and transmission Tvξcoefficients, which depend on: wave number in the free space k0; angles θ, φand other geometrical and electrodynamical parameters of the problem. The coefficients with similar lower indices (v=ξ)describe conversion of the incident wave into the wave with the same polarization. By analogy, reflection and transmission coefficients with dissimilar lower indexes (vξ)describe conversion of the incident wave into the wave with the orthogonal polarization. In that notation, the left lower index νcorresponds to the polarization of the reflected/transmitted wave; the right lower index ξ corresponds to the polarization of the incident wave.

Notice, that the presence of the “crossed” reflection (Rsp,Rps)and transmission coefficients (Tsp,Tps), which are responsible for incident plane wave depolarization, is the specific properties of anisotropic media (see, for example [31]).

As it follows from the expressions (32) – (33):

  • for the case of s– polarized incident plane wave (As=1,Ap=0):

Rss=h(κ,0)1,Rps=e(κ,0),Tss=h(κ,b),Tps=e(κ,b);E34
  • for the case of p– polarized incident plane wave (As=0,Ap=1):

Rpp=e(κ,0)1,Rsp=h(κ,0),Tpp=e(κ,b),Tsp=h(κ,b).E35

The impedance boundary conditions (25) for the scalar potentials can be rewritten in the following form:

h/z=ik0(λsp(a)e+λss(a)h)+fs,(z=0);e/z=ik0(λpp(a)e+λps(a)h)+fp,E36
h/z=ik0(λsp(u)e+λss(u)h),(z=b+0).e/z=ik0(λpp(u)e+λps(u)h),E37

The values λvξ(a,u),fvare depending on the angles θand φand have the form:

λpp(a)=1μzz(0)[aμ(0)sinθμlz(0)cosθ],λps(a)=bμ(0)μzz(0),λsp(a)=bε(0)εzz(0),λss(a)=1εzz(0)[aε(0)sinθεlz(0)cosθ];E38
λpp(u)=1μzz(b)[aμ(b)ncμc+μlz(b)cosθ],λps(u)=bμ(b)μzz(b),λsp(u)=bε(b)εzz(b),λss(u)=1εzz(b)[aε(b)ncμc+εlz(b)cosθ];E39
fs=2ik0aε(0)sinθAs/εzz(0),fp=2ik0aμ(0)sinθAp/μzz(0).E40

Let us build the finite-difference procedure that approximately describes the system of coupled differential equations (18) and the boundary conditions (36) – (37). Taking into account the fact, that the external sources are absent inside the anisotropic layer, the substitutions qp=qs0are required in formulas (18).

At first, we divide the segment b<z<0on the N equal parts. After that the grid step is immediately obtained as Δb=b/N, and then the grid’s knot set z0,z1,...,zNis defined by formula zj=jΔbb,(j=0,1,...,N). Notice, that this knot set includes boundary points z0=band zN=0as well. The partial derivatives in the differential equations (18) are approximated by the central difference; in the boundary conditions (36) – (37), they are approximated by the left-hand difference and the right-hand difference, correspondingly. Then, we obtain the system of linear algebraic equations with dimension 2N+2for the unknown complex variables xkh(κ,zk),yke(κ,zk),(k=0,1,...,N). For the case of s-polarized wave this system has the form:

{xN(1ik0bλss(a))xN1ik0bλsp(a)yN=fss,ik0bλps(a)xN+yN(1ik0bλpp(a))yN1=0;E41
{Ajxj+1+Bjxj+Cjxj1+Djyj+1+Fjyj+Gjyj1=0,Pjxj+1+Qjxj+Rjxj1Kjyj+1LjyjMjyj1=0,(j=1,2,...,N1);E42
{x1x0(1+ik0bλss(u))ik0bλsp(u)y0=fss,y1y0(1+ik0bλpp(u))ik0bλps(u)y0=0.E43

Here, the first equation set (41) is corresponding to the finite-difference approximation of the boundary conditions (36) imposed at z=0, As=1, Ap=0. The next equation set (42) is the finite-difference approximation of the differential equations (18) and the next one equation set (43) is the finite-difference approximation of the boundary conditions (37) imposed z=b;fss=Δbfs.

The coefficients Aj,...,Mjthat enter the equations (42) are given by the expressions:

Aj=εzz(zj)aε(zj)+14[εzz(zj+1)aε(zj+1)εzz(zj1)aε(zj1)]+ik0bcosθ2[εzl(zj)aε(zj)+εlz(zj+1)aε(zj+1)],Bj=2εzz(zj)aε(zj)+(k0b)2[|μ^(zj)|aμ(zj)cos2θεll(zj)aμ(zj)],Dj=ik0b2[dμ(zj)aμ(zj)+bε(zj+1)aε(zj+1)],Cj=εzz(zj)aε(zj)14[εzz(zj+1)aε(zj+1)εzz(zj1)aε(zj1)]ik0bcosθ2[εzl(zj)aε(zj)εlz(zj1)aε(zj1)],Fj=(k0b)2cosθ[cε(zj)aμ(zj)eμ(zj)aμ(zj)],Gj=ik0b2[dε(zj)aμ(zj)bε(zj1)aε(zj1)],AjKj,BjLj,CjMj,DjPj,FjQj,GjRj,(εμ).E44

The system of equations for the p-polarized incident plane wave, can be received from (41) – (43) by the substitutions fss0and 0fpp=Δbfpin the systems (41), (43). For the case, when anisotropic layer is uniaxial media: f(z)=g(z)0.

As a result, a set of linear algebraic equations (41) – (43) with dimension 2N+2is derived. Obtained linear system of equations can be solved by standard techniques such as Gauss method. By solving it we find the unknown transmission and reflection factors for the structure under study.

3. Results and discussion

In this section we present the results of numerical simulations that illustrate the influence of anisotropy and inhomogeneities of the materials composed of layers in the structure being studied upon the mechanisms of diffraction of the incident plane electromagnetic wave.

All structures presented here are based on the porous silicon (PSi). Today, PSi plays an important role in a number of applications. These include microcavities [11,32], photonic crystals [33], waveguide structures [34], photodetectors [35], sensors [36], etc. Besides, PSi has the potential to be an optically active material in the case when an acceptable electro- or thermo-optic media is infiltrated into the pores [32]. Therefore, porous silicon is an excellent candidate for tunable optical interconnects and optical switches. For all these applications a strict control over the reflectance and transmission properties of PSi layers is required.

Today porous silicon attracts a great deal of attention because it’s a material with great technological promise. The main advantages of PSi may be summarized as follows:

  • PSi is a simple and low cost dielectric material that can be easy prepared;

  • PSi is a promising material for photonic applications due to its excellent thermal and mechanical properties, obvious compatibility with standard Si-based technologies;

  • PSi is a suitable material for the formation of arbitrary multilayers. PSi multilayers are almost arbitrary combination of layers with different thickness and porosity (refractive index), because these two parameters can be relatively easily controlled during the formation process of porous silicon.

It is well known that the PSi films can be produced by anode electrochemical etching of the monocrystalline silicon plates [37]. The nanometer-size pores tend to grow in the direction of electrochemical etching and, accordingly, nanocrystal formation sets in. The porosity P and the effective refractive index neff=εeffof PSi are controlled by the current density under electrochemical etching, because the effective refractive index of PSi is determined by the porosity and refractive index of the medium inside the pores. Thus, by periodically varying the magnitude of current density we are able to obtain the structure with alternating layers of different porosity and, consequently, with different refractive indices.

The bulk silicon crystal is not birefringent due to its cubic crystal symmetry. However, porous modification of silicon can exhibit strong in-plane anisotropy of the refractive index [37,38]. The observed birefringence depends on the porosity, the size of Si nanocrystals, the spacing between them and the dielectric properties of surrounding medium. As was demonstrated in [39], PSi layers, with dimension of the pore about 10-30 nm, have properties of the negative uniaxial crystal with diagonal permittivity tensor

ε=[n02000n02000ne2],E45

whose birefringence magnitude Δn=noneis up to 0.24. In expression (45) n0is the index of refraction for the waves polarized perpendicularly to the optical axis, which are called “ordinary” or “o – waves”; neis the index of refraction for the waves polarized parallel to the optical axis, which are called “extraordinary” or “e – waves”. It is important to note that if ne>n0the crystal is said to be positively uniaxial, in opposite case if ne<n0the crystal is said to be negatively uniaxial. When a linearly polarized wave of arbitrary polarization direction enters an anisotropic medium, it will be split into two components polarized along the two allowed polarization directions which are determined by the direction of the wave vector relative to the axes of the indicatrix. As a result, each s- or p- polarized plane wave incident on such an anisotropic photonic structure will generate two reflected and two transmitted plane waves containing both s- and p- polarized planewaves. For the special cases when the principal axes of the layers are parallel or perpendicular to the fixed axes, the s-and p-polarized waves remain uncoupled [26].

Two types of PBG structures are investigated, namely, a distributed Bragg reflector (DBR) in other words photonic crystal (PhC) and a microcavity.

Figure 4.

Schematic presentation of PSi-based PhC: the dark layers have high porosity (low refractive index) and the bright layers are of low porosity (high refractive index)

Figure 5.

Reflectance spectra Rss of the PhC: solid line – numerical simulation; scatter – experimental results [39]. The calculations involved the value of sloping and azimutal angles: θ=700, θa=450 and φ=φa=00

Figure 6.

Reflectance spectra Rss as function of incident sloping angle θ. In this case θa=450, λc=800 nm

The simplest multilayered PBG structure is one-dimensional photonic crystal as depicted on Fig. 4. It’s well known that PhCs are class of optical media represented by the natural or artificial structures with periodic modulation of the refractive index. Such optical media have some peculiar properties which gives an oportunity for a number of applications to be implemented on their basis. The most important property which determines practical significance of the PhC is the presence of the omnidirectional photonic band gap. The PBG refers to the energy or frequency range where the light propagation is prohibited inside the PhC. As an example of such a PhC one can give a Bragg grating which is widely used as a distributed reflector in vertical cavity surface emitting lasers. Besides, such structures are widely used as antireflecting coatings which allow dramatically decrease the reflectance from the surface and are used to improve the quality of lenses, prisms and other optical components.

Figure 7.

Reflectance spectra Rss as function of the incident azimuthal angle φ

Figure 8.

Reflectance spectra Rps as function of the incident azimuthal angle φ

First, we study light propagation in a one-dimensional photonic crystal, which was identical to PhC that was experimental investigated by Aktsipetrov et al. [39]. It consists of 25 pairs (building blocks) of lossy anisotropic layers with refractive indices: no1=1.39+0.008i,ne1=1.32+0.008i,no2=1.58+0.008i,ne2=1.5+0.008i.The materials are assumed to be nonmagnetic so that μ=1throughout the whole layered medium. The physical thickness of the layers was chosen such that the optical thickness of layers was equal to λc/4, where λc=800nmis the Bragg wavelength corresponding to the photonic band gap (PBG) centre at the normal incidence (θ=90). The numerical modeling were made for the case where the optical axes of all structure layers were oriented in one and the same direction.

The comparison of the results of the proposed numerical scheme with results of the cited experimental work [39] is presented on Fig. 5. The reflectance spectra demonstrate existence of a PBG with a reflectance of about 0.9 in the wavelength region of 770–860 nm. As we see theoretical result is in very good agreement with the experimental one. So proposed method can be successfully used for computation of the spectral characteristics of 1D anisotropic layered structures.

For all the structures discussed here, the reflection/transmission spectrum is strongly dependent on the incident sloping θand azimuthal φangles. The dependences of the reflection spectra of PSi-based distributed Bragg reflectors from variation of the sloping angle θand the azimuthal angle φof the incident s-polarized plane wave are shown in Fig. 6. – Fig. 8. In all cases we take λc=800nm.

It may be seen from the Fig. 6 that in agreement with the theory the width of the high reflectance region (width of the photonic band gap) is decreased with decreasing the sloping angle θ, and the central wavelength of PBG is shifted to the short wavelength region. In this case we assume that φ=φa=00in result the conversion incident wave to wave with orthogonal polarization is absent.

The influence of the azimuthal angle φon the conversion of the incident plane wave is show in Fig.7-8. The curves presented are calculated at θ=800, θa=450and φa=00. Note that in a uniaxial crystal the maximum angular separation of the “o” and “e” waves, in other words the maximum conversion of the incident linearly s- or p- polarized plane wave into the cross-polarized wave, occurs when the wave vector has the angle φmax450with the optic axis [31]. Also we should note that the value of φmaxis proportional to |none|. Analyzing the figures 7-8, we can see that when the incident azimuthal angle φvaries:

  • the reflectance spectra demonstrate conversion of the incident s-polarized plane wave into the wave with same (Fig. 7) and orthogonal (Fig. 8) polarization. It occurs within the wide range of the incident azimuthal angles;

  • the spectra also show a shift in the PBG spectral position when the azimuthal angle is changed. The largest shift, about 20 nm, is observed under variation of the incident azimuthal angle φfrom 00 to 900. This results are in agreement with [39];

  • the maximum conversion of the incident s-polarized wave into the wave with orthogonal polarization occurs at the incident azimuthal angle φ450(see Fig. 8(b)), that is in conformity with [31]. In this case two reflection peaks with magnitude about 0.3 are clearly observed in the reflectance spectra Rps.The first reflection peak with the central wavelength 780 nm is corresponding to the short-wavelength PBG edge for case φ=900. The second one with the central wavelength 835 nm – to the long-wavelength PBG edge for case φ=00

The influence of the sloping and the azimuthal angles θa,φa(which determine direction of the optical axes) on the conversion incident s-polarized plane wave into wave with orthogonal polarization at the Bragg wavelength (λc=800nm)is shown in Fig. 9(a,b). Reflection coefficients Rpsis calculated at incident angles θ=800and φ=00. It is clear from this graph that:

  • the maximum value of Rpsis around 0.4 (see Fig. 9(b)) and corresponds to the sloping angle θa=00and the azimuthal angles φa=mπ/4, (m=1,3,5,7)that is in accordance with [31];

  • in the case φa=mπ/2, (m=0,1,2,3,4)the conversion of the incident wave into the wave with orthogonal polarizations is virtually absent Rps0.

Finally, we can conclude that reflection coefficients of the investigated structures are very sensitive to the azimuthal angle φaand the incidence angles θand φ. When the angles of incidence θ=900φ=00and the azimuthal angle φapossesses values 0° or 90°, the s- and p-modes are almost uncoupled, and the value of the Rpsremains very small. The situation is similar to that one in the case of isotropic materials. With the exception of these particular values of φa, there is a great mixing between s-and p-modes that gives rise to the emergence of a large reflection coefficient Rps.

Figure 9.

Reflection coefficient Rps as function of the sloping and the azimuthal angles

The PSi-based mirrors studied above are periodic structures, but one important property of PhC is the presence of narrow resonance (localized defect modes) in the PBG region when a disorder is introduced in their periodic structure. Usually “defect” is a layer with half-wavelength optical thickness that inserted in the middle of the dielectric stack. Fig. 10 shows a scheme of a typical PSi-based microcavity structure consisting in an active layer sandwiched between two distributed Bragg reflectors (PSi mirrors). At that, the reflectance spectrum of the whole structure is changed. As a rule, when the parameters of defect layer are properly selected, sharp transmittance peak appears within the main reflectance band. The wavelength of the transmittance peak corresponds to the resonant wavelength of the defect. For instance, if the thickness of the defect layer is twice larger than it was in defectless structure, the transmittance peak appears at the Bragg wavelength of the corresponding defectless structure. Thus, if the radiation with wavelength equal to Bragg wavelength of defectless reflector falls at this structure it will pass the structure almost without the reflectance. If the defect thickness will be slightly different, the wavelength of the transmittance peak will be different as well.

Figure 10.

Schematic presentation of PSi-based microcavity

The advantage of using PSi microcavities is that the position of the transmittance peak is completely tunable by changing the properties of the central layer (i.e. porosity and thickness) during the electrochemical etching condition and by the infiltration of organic molecules (for example liquid crystals [32]).

Now, let’s discuss the influence of anisotropy of layers on the optical properties of the microcavitiy. In cited case the PSi-based microcavity consist of two mirrors (each of them consists of 12 building blocks) separated by active layer with optical thickness λc/2and low porosity.

The Fig. 11 shows the reflectance spectra of the microcavity for both s- (solid curve) and p- (dashed curve) polarized incident plane waves. It should be noted that spectra reveal the presence of photonic band gap with a reflectance of about 0.85-0.9 and microcavity mode within the PBG for both polarizations. Analyzing both spectra together, we can observe a spectral shift of the microcavity mode when the incident polarization is changed. This particular effect is also the manifests an attribute of birefringence and can be used in dividing the incident radiation polarization. In our case, the central wavelengths of resonance peaks are λsc=796nmand λpc=817nmfor s- and p- polarization, respectively. Taking into account this result, we can consider that the variation in the incident radiation polarization brings about the shift Δλ=|λpcλsc|of the microcavity mode, and this shift may be as high as 21 nm for presented case.

Figure 11.

s - and p- polarized reflectance spectra of microcavity. In this case λc=800 nm, θ=800, φ=φa=00, θa=450

Figure 12 plots the spectral shift of the microcavity mode for p-and s-polarized light, as a function of the of azimuthal angle φa. From this figure, some particular effects can be highlighted:

  • the orientation of the optical axes of the layers with respect to incidence plane of s- or p-polarized waves influences on the quantity of resonant peaks within PBG area and its location on the wavelength scale. Particularly, with the increasing of azimuthal angle, the “additional” well recognized resonant peak is appeared in the stopbands for both polarizations of the incident plane waves. It correlates with results of the paper [40];

  • as a general trend, the value of the “main” resonant peak is decreasing with increasing the azimuthal angle φaand it completely disappears when φa=900. There against, the magnitude of the “additional” resonance peak is increasing with increasing φaand it has maximum value in the mentioned above case;

  • both transmission peaks have identical magnitudes when azimuthal angle φa480for both s- and p- polarization;

  • the extreme cases, namely φa=00and φa=900, correspond to presence of single-defect mode within PBG region, while in the other cases this mode splits into two defect modes. The comparison of the results presented in Fig. 12(a) and 12(b) shows that λs0=λp90=796nmand λs90=λp0=817nm. Here we used next definition for central wavelengths – λxy. In this notation, the left subscript index x=s,pis corresponding to polarization of incident wave; the right subscript index y=0,90is corresponding to quantity of azimuthal angle φa;

  • also, in the extreme cases, we observed the shift of the resonance peaks Δλ=|λp0λp90|=|λs0λs90|about 21 nm to the short-wavelength and long-wavelength regions for p- and s-polarization, respectively. This result is in good agreement with the experimental results [39].

Figure 12.

Microcavity resonance shift based on the azimuthal angle φa change. Reflectance spectra for p- and s- polarized incident plane wave are depicted on Fig. 12(a) and Fig. 12(b), respectively.

The Fig. 13 shows characteristics similar to those one presented on Fig.12, but for cross-polarized components Rspand Rps. In these figures we can see that within the PBGs areas the reflection coefficients Rspand Rpshave local minimum that occur around the wavelength λ(λp0+λp90)/2(λs0+λs90)/2.

Figure 13.

Reflectance spectra Rsp (a) and Rps (b) as function of the azimuthal angle φa. In this case λc=800 nm, θ=800, φ=00, θa=450.

4. Conclusions

In conclusion, the present chapter describes the mathematical background for calculation of spectral characteristics of the birefringent layered media. On the basis of the presented theoretical description and its numerical approximation the influence of material anisotropy of PSi-based layered photonic structures on their optical properties has been theoretically investigated. All these multilayer PBG structures have been designed for 0.8 μm applications.

The reflectance spectra of some photonic structures for both polarization of the incident plane wave are calculated. The agreement between the numerical calculations and the experiments [39] are obtained. It was shown numerically that anisotropy of layered media reduces to polarization transformation of the incident plane wave. Notably, the maximum conversion of incident plane wave into the wave with orthogonal polarization occurs when the optical axis of the structure have angles φa=mπ/4, (m=1,3,5,7)with respect to the incidence plane of s- or p-polarization incident plane. In contrast to this, the conversion of the incident wave into the wave with orthogonal polarizations is absent when φa=mπ/2, (m=0,1,2,3,4).

Also from our numerical calculations, we conclude that the location of the PBG edges and location of microcavity modes within PBG region are different for s- and p- polarized waves and their spectral positions change under rotation of the optical axis of the structure with respect to the wave vector of the incident electromagnetic wave. The variation of polarization of the incident radiation brings about the shift of the microcavity mode, and this shift may be as high as 21 nm for discussed case. This particular effect can be used in the devices that divide the incident radiation according to polarization.

In addition the obtained results can be used in the designing of PSi-based photonic devices:

  • spectral- and polarization-selective elements (filters);

  • electrically tunable filters and optical switches. In this instance it is required that the pore be filled with liquid crystals;

  • polarization converters;

  • time-delay elements for optically controlled phased-array antenna system.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

V.I. Fesenko, I.A. Sukhoivanov, S.N. Shul’ga and J.A. Andrade Lucio (February 13th 2013). Propagation of Electromagnetic Waves in Anisotropic Photonic Structures, Advances in Photonic Crystals, Vittorio M.N. Passaro, IntechOpen, DOI: 10.5772/54847. Available from:

chapter statistics

1906total chapter downloads

1Crossref citations

More statistics for editors and authors

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

Access personal reporting

Related Content

This Book

Next chapter

Threshold Mode Structure of Square and Triangular Lattice Gain and Index Coupled Photonic Crystal Lasers

By Marcin Koba

Related Book

First chapter

Gaseous Scanning Electron Microscope (GSEM): Applications and Improvement

By Lahcen Khouchaf

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

More About Us