^{1}

^{*}

^{1}

^{2}

^{1}

This paper studied the transmission coefficient and escape frequencies in a system of planar dielectric layers where the refractive index changes from one layer to another through a Gaussian function. The wave equation with normal incidence is analyzed. For the calculations, the transfer matrix formalism is used. In a previous work, the transmission and escape problem for Gaussian electronic superlattices is investigated. Now it studied the electromagnetic modes for a system formed by layers where the refractive index of the structure is modulated by a Gaussian function. The system presents transparency bands of transmission and gaps without transmission. The escape frequencies are situated near these transparency bands but they do not coincide with them. is the frequency (mode) and Γ describes the width of the states. For these systems, the escape states are very wide. A non Gaussian system presents resonance peaks in the transmission and the escape states are narrow. The formation of transparency bands in the transmission for a Gaussian system is attributed to the widening of the escape states.

When we talk about photonic crystals, we refer to periodic structures as systems used to manipulate light. At the 80’s years, Yablonovich and Sajeev proposed to make artificial structures with dielectric material for the study of the electro- magnetic properties of crystal photonics. The first one was interested in inhibiting spontaneous emission, while the second was interested in the study of light introducing a small localized disorder in the periodic system [

The study for energy, electronic, optical and acoustic filters is an interesting and active field. Formerly frequency band pass filters where light can be transmitted in full when the frequency of the incident photons corresponds to the passband and total reflection occurs when the incident energy corresponds to the stopband, were studied. For example, the study of periodic profile systems, i.e. by bilayers formed with refractive indices n 1 and n 2 systems correspondingly [_{1} and n_{2}. In regular systems there are transmissions bands which present resonance peaks inside the bands with high values of transmittance but there is not transparency bands.

On the other hand, we mean by escape frequencies the resonance frequencies for an open system which is separated by the exterior by a partly transparent boundary surface. This structure losses energy to the exterior via radiation and we do not consider radiation coming into the system from the exterior. In this case the frequencies obtained are complex eigenvalues. The formation of transparency bands in a Gaussian system is an outstanding fact. We think that these transparency bands are due to the smoothness of the Gaussian change of refractive index which facilitates the transmission of the waves. On the other hand, this kind of variation of refractive index, with high value of index in the middle of the structure and progressive reduction to the extremes of the system, would make the waves to escape more easily from the structure. This means that there would be a reduction of the states lifetime and a widening in frequency. Then, the transparency bands would be the envelope of the wide escape states. Thus, our purpose in this work is to compare the transparency bands with the escape states.

For electromagnetic waves, in [

The paper is organized as follows. In Section 2, we present the theoretical model. We use the transfer matrix formalism. In Section 3, we describe the structures we are interested in. In Section 4, we present our results. Finally, in Section 5 we formulate our conclusions.

We can follow the same idea as in [

In the case of a polarization wave s where the electric field vector E is transverse to the plane of incidence can be determined by applying the conditions of continuity in the interface. All electric field are perpendicular to the incident plane, with outward direction of the sheet, while two magnetic field vectors are taken so that the flow energy is positive in the direction of propagation. We describe the propagation light with variable refractive index by the equation

d d x [ μ 0 μ ( x ) d E z ( x ) d x ] + [ ω 2 μ 0 ε ( x ) − μ 0 μ ( x ) κ 2 ] E z ( x ) = 0 , (1)

where E z ( x ) and μ 0 μ ( x ) d E z ( x ) d x are continuous functions. The equation of

motion is derived from the Maxwell equations in the substance,

∇ ⋅ D ( r , t ) = ρ ( r , t ) , (2)

∇ × E ( r , t ) = − ∂ B ( r , t ) ∂ t , (3)

∇ ⋅ B ( r , t ) = 0 , (4)

∇ × H ( r , t ) = j ( r , t ) + ∂ D ( r , t ) ∂ t . (5)

If there are no sources, the fields are of the form

E t ( ρ , x , t ) = E 0 t ( x ) e i ( κ ⋅ ρ − ω t ) , (6)

H t ( ρ , x , t ) = H 0 t ( x ) e i ( κ ⋅ ρ − ω t ) , (7)

E z ( ρ , x , t ) = E 0 z ( x ) e i ( κ ⋅ ρ − ω t ) , (8)

H z ( ρ , x , t ) = H 0 z ( x ) e i ( κ ⋅ ρ − ω t ) , (9)

with

ρ = y e y + z e z , (10)

κ = κ y e y + κ z e z , (11)

E 0 t ( x ) = E 0 y ( x ) e y + E 0 z ( z ) e z , (12)

H 0 t ( x ) = H 0 y ( x ) e y + H 0 z ( x ) e z . (13)

With a little algebraic manipulation we can obtain (1) considering that D x ≡ 0 . We study one-dimensional systems where the profile for the refractive index n s y s ( x ) is described as Section 3.

We concentrate in a structure of N plane dielectric layers perpedicular to the x axis. The width of each j layer are denoted by d j and each refractive index by n j . We also consider the study in medium homogeneous, isotropic and lossless. In other words, we consider a system where at its ends, the refractive indices are constant ( n L and n R ), while in the intermediate zone the refractive index n ( z ) is dependent on the position (the growth direction of the heterostructure). We use the formalism of transfer matrix conformed by the dynamical matrix D j and the propagation matrix P j [

We know that the energy flux is given by the Poynting vector S . Imposing the condition of continuity in the interface ( x = 0 ) for E y and H z we have

E 1 + E ′ 1 = E 2 + E ′ 2 , (14)

ϵ 1 μ 1 ( E 1 − E ′ 1 ) cos θ 1 = ϵ 2 μ 2 ( E 2 − E ′ 2 ) cos θ 2 . (15)

θ 1 and θ 2 are the angles of the wave vectors k 1 and k 2 , respectively, with respect to the normal to the interface. These two equations can be rewritten in matrix form as

D 1 ( E 1 E ′ 1 ) = D 2 ( E 2 E ′ 2 ) . (16)

Then we can rewrite the continuity condition in N layers for every interface through a matrix

D j = ( 1 1 ϵ j μ j cos θ j − ϵ j μ j cos θ j ) . (17)

D j is so-called dynamic matrix s for j = 1 , 2 , ⋯ , N .

Analogously we can the same methodology followed for polarization p [

[ A 0 B 0 ] = [ M 11 M 12 M 21 M 22 ] [ A s B s ] , (18)

A 0 , B 0 , A s , B s corresponding to coefficients incoming and outgoing coefficients media respectively. We define the matrix P j which is so-called propagation matrix, it considers the phase change due to the thickness of each layer

P j = [ e − i ϕ j 0 0 e i ϕ j ] . (19)

Φ j is given by ϕ j = k j d j = n j ω d j c cos θ j where k j = n j ω c cos θ j . It is

possible to show that a layered system such as that shown in

D 0 − 1 [ ∏ i = 1 N D j P j D j − 1 ] D s = [ M 11 M 12 M 21 M 22 ] . (20)

The system is composed of N layers, where D 0 and D s refers to matrices of the corresponding incoming and outgoing media. Details can be seen in [

We are interested in studying the scattering problem of (18). In this problem we consider the incoming, reflected and transmitted waves, then the coefficient B s is null. The transmittance T is given by the ratio of the Poynting power flow of the transmitted wave to that of the incidente wave, then

T = n s n 0 | 1 M 11 | . (21)

The reflection coefficient R is obtained in a similar way, or taking into account that R + T = 1 .

The problem of escape is related to having photons confined in the crystal and they leave the system. We consider only outgoing waves, and the coefficients A 0 = 0 and B s = 0 out of the interval ( x 0 , x s ) . Then (18) gives the transcendental equation,

M 11 = 0. (22)

The solutions to (22) have the form ω = ω r − i Γ , where the real part ω r represents the frequency while the imaginary part Γ describes the fact that the modes have finite lifetime and they decay. For simplicity we make a change of variable, in the Section 3 we will explain in detail.

We are interested in a photonic crystal where the refractive index of the structure is modulated by the Gaussian function (see

proposal in [

where n max is the refractive index of the value of the central layer, n min the minimum value of the initial and final layer and σ is the standard deviation. To guarantee that the initial and final layer take the value 1.2 we take

n min = 1.2 − n max exp ( − x / σ ) 2 1 − exp ( − ( x / σ ) 2 ) (23)

We do this in order that the traveling wave traverses the structure in a smooth way as possible. In all cases studied we take n max = 2 . The even layers are constant and for our case we consider them with a value equal to the exterior of the structure. We show a result for system where the odd layers are modulated by a Gaussian function so that the maximum refractive index is at the ends and the minimum at the center of the structure to compare them (Section 4).

The calculations are performed at normal incidence with non-magnetic materials. The ratio of the widths of the layers d e v e n / d o d d = 8 . It was normalized in the frequency range and instead of plotting in function of ω the results are presented in function of f Δ / c . In this case we consider the total width of the structure as Δ = 1 and f is the frequency^{1}.

We refer to regular system to structures formed by refractive indices n 1 and n 2 . The refractive indices are interleaved with odd number of layers N c = 3 , 5 , 7 , ⋯ . For example, for N c = 7 we have refractive indices n 1 , n 2 , n 1 , n 2 , n 1 , n 2 , n 1 and with its corresponding widths d 1 and d 2 .

^{1 ω = 2 π f }.

In this section we present the transmission coefficients and escape frequencies for systems in which the profiles are formed for layers modulated by a Gaussian function. We show it for 17, 25, 37, 45 layers with Gaussian modulation (see

system the odd layers are modulated by a Gaussian function so that the maximum refractive index is at the ends and the minimum at the center of the structure. The regular system case consist in that the value of the refractive index in the odd layers is constant and equal.

We use allowable values for porous silicon. We can observe that the greater the number of layers the passband is better defined.

We made the calculations using Mathematica 10.4. It is important to stress that it was difficult to solve the escape problem in (10) for more than 11 layers. This was overcome separating the element matrix M 11 equals to zero in its real and imaginary part and using an interpolation with known numerical methods.

From Figures 2-4 we see that the Gaussian systems present transparency bands with gaps of no transmission. For these systems the escape states does not necessarily coincide with the resonaces of the passbands and the escape states

are very wide. By contrast, the reversed Gaussian and the regular structures present passbands without transparency, only peaks of resonances of transmission for the eigenstates of the system. At the same time, for these structures, the coincidence of escape states and resonances of transmission is much better and the escape states are narrow. This means that the escape states lifetime for Gaussian systems is shorter than for a reversed Gaussian and regular systems, which is what we expected because the photons are more difficult to escape from these last systems. By contrast, the progressive reduction of refractive index from the middle of the structure to the extremes, facilitates the escape of photons in a Gaussian system . Simultaneously, from the transmission bands and the position of the wide escape states for Gaussian structures, we deduce that the transparency bands are formed by the envelope of the wide escape states, which was the expectation mentioned in the Introduction. Mathematically speaking we are facing two different boundary problems on the same equation. The solutions need not be the same. In fact, the transmission problem has a continuum spectrum while the escape problem has a discrete spectrum with complex values. However, in certain situations, as the case of the regular structure, the two problems in some extent appear to be closely related. In the case of the Gaussian structure the relationship between the two problems is less clear.

Using the Maxwell’s equation and the transfer matrix formalism, we obtained transmission coefficients and the solution of the problem of escape in systems where the refractive index of the layers is modulated by a Gaussian function, and for a reversed Gaussian and a regular systems. The variation of parameters such as the refractive index and the number of layers in Gaussian structures, allow passbands with very good transmission, i.e., very good transparency. This type of systems has potential applications as frequency filters for photons.

The structure has broad intervals of frequency or passbands where there is almost total transmission, separated by stopbands where there is no propagation of photons. The escape frequencies have very wide linewidths and are situated inside or near the passbands, but they do not necessarily coincide with the passbands. For Gaussian structures, the escape frequencies have very wide linewidths and are situated inside or near the passbands, but they do not necessarily coincide with the passbands. For a regular and and inverted Gaussian structures the resonance and the escape frequencies are very close. We associate the large width of the escape frequencies to the formation of transparency bands in the Gaussian structures.

MLSV and DACS thank the support of PRODEP SEP-SES. CIC wishes to thank the support of COZCYT and also of CONACYT (Grant Number 337137). RPA acknowledges hospitality at Autonomous University of Zacatecas.

de la Luz Silba- Vélez, M., Contreras-Solorio, D.-A., Pérez- Álvarez, R. and Cabrera, C.I. (2017) Finite One-Dimensional Photonic Crystal with Gaussian Modulation: Transmission and Escape. Optics and Photonics Journal, 7, 170-180. https://doi.org/10.4236/opj.2017.710017