Next Article in Journal
Crystal Structures of the 43 kDa ATPase Domain of Xanthomonas Oryzae pv. Oryzae Topoisomerase IV ParE Subunit and its Complex with Novobiocin
Next Article in Special Issue
Spatial Beam Filtering with Autocloned Photonic Crystals
Previous Article in Journal
Energy Storage Performance of Sandwich Structured Pb(Zr0.4Ti0.6)O3/BaZr0.2Ti0.8O3/Pb(Zr0.4Ti0.6)O3 Films
Previous Article in Special Issue
Magneto-Optical Isolator Based on Ultra-Wideband Photonic Crystals Waveguide for 5G Communication System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Analysis of Photon Density of States for One-Dimensional Photonic Crystals in a Rectangular Waveguide

1
School of Information Technology, Guangdong Industry Polytechnic, Guangdong 510300, China
2
Multidisciplinary Computational Laboratory, Department of Electrical and Biomedical Engineering, Hanyang University, Seoul 04763, Korea
*
Author to whom correspondence should be addressed.
Crystals 2019, 9(11), 576; https://doi.org/10.3390/cryst9110576
Submission received: 1 August 2019 / Revised: 26 October 2019 / Accepted: 30 October 2019 / Published: 4 November 2019
(This article belongs to the Special Issue Sonic and Photonic Crystals)

Abstract

:
Light propagation in one-dimensional (1D) photonic crystals (PCs) enclosed in a rectangular waveguide is investigated in order to achieve a complete photonic band gap (PBG) while avoiding the difficulty in fabricating 3D PCs. This work complements our two previous articles (Phys. Rev. E) that quantitatively analyzed omnidirectional light propagation in 1D and 2D PCs, respectively, both showing that a complete PBG cannot exist if an evanescent wave propagation is involved. Here, we present a quantitative analysis of the transmission functions, the band structures, and the photon density of states (PDOS) for both the transverse electric (TE) and transverse magnetic (TM) polarization modes of the periodic multilayer heterostructure confined in a rectangular waveguide. The PDOS of the quasi-1D photonic crystal for both the TE and TM modes are obtained, respectively. It is demonstrated that a “complete PBG” can be obtained for some frequency ranges and categorized into three types: (1) below the cutoff frequency of the fundamental TE mode, (2) within the PBG of the fundamental TE mode but below the cutoff frequency of the next higher order mode, and (3) within an overlap of the PBGs of either TE modes, TM modes, or both. These results are of general importance and relevance to the dipole radiation or spontaneous emission by an atom in quasi-1D periodic structures and may have applications in future photonic quantum technologies.

1. Introduction

Photonic crystals (PCs), also known as artificial materials, have attracted much attention in the past three decades due to the tremendous needs of gaining complete control over light propagation and emission [1,2]. PCs, according to the dimension of the periodicity, are divided into three categories: one- (1D), two- (2D), and three-dimensional (3D) crystals. Due to the periodicity, the stop band and the pass band are formed, according to the Bloch theorem [3]. Therefore, periodic dielectric materials are characterized by photonic band gaps (PBGs). As an analogy to the electronic band gaps in solid state materials [4], a PBG in PCs can prohibit the propagation of electromagnetic (EM) waves whose frequencies fall within the band gap region. The continuing success of the semiconductor industry in controlling electric properties of materials from the last century has encouraged us to manipulate the flow of light in PCs and control their optical properties. These PC-based materials are expected to have many applications in optoelectronics, optical communications, and photonic quantum technologies in the next decades [5]. In optical range, PCs have been extensively studied. It was proposed that the emission of EM radiation can be modified by the environment [6,7]. Several environments such as metallic cavities [8], dielectric cavities [9], superlattices [10,11,12,13,14,15], and 2D PCs [16] have been studied. The environmental effects have been described by the photon density of states (PDOS), which is related to the transition rate of the Fermi golden rule. In principle, a PC with a complete PBG can be best realized in a 3D system to prohibit the propagation of electromagnetic waves of any polarization traveling in any direction from any source [1]. However, the difficulty in fabricating such 3D crystals with PBGs in the optical regime prohibits the progress of many applications.
On the other hand, there has been a lot of interest in microwave and millimeter wave applications of PBG, such as the significant progress in the design of filters [17,18], microstrip antennas [19], and slow wave structures [20,21], and so on. However, the design of PBG in this frequency range is still difficult due to complexities of the modeling. There are too many parameters affecting the PBG properties, such as the number of lattice periods [22], lattice shapes [23,24], lattice spacing [25], and relative volume fraction [26,27,28,29,30]. Since the actual fabrication of 3D PCs remains difficult, another simpler choice is periodic dielectric or PC waveguides, which have only a one-dimensional periodic pattern [1]. The rigorous study of the PC waveguide can be traced back to the 1970s, where a more detailed review can be found [31,32]. Recently, it was demonstrated that, by considering a quantum dot spin coupled to a PC waveguide mode, the light–matter interaction can be asymmetric, leading to unidirectional emission and a deterministic entangled photon source, which might have application in future optical quantum devices [5,33]. One interesting feature of electromagnetism in dielectric media is that there is no fundamental length scale, namely the scaling properties of Maxwell’s equations, i.e., the solution of problem at one length scale determines the solutions at all other length scale [1]. In a previous work [34], a multilayer dielectric window in a rectangular waveguide had been studied to achieve a wide-bandwidth transmission of a millimeter wave. A transfer matrix approach was successfully employed to discretize the dielectric function profile and the transmission functions could be calculated efficiently. In principle, the approach can be extended to study a quasi-1D PC, a PC confined in a waveguide. However, the transmission method is limited to study radiation modes in a finite-length system. In order to study the PBG phenomena such as the suppression of spontaneous emission [35] in a quasi-1D PC, the calculation of the dispersion relations or band structures (BS) and the PDOS are needed. Metallic waveguides and cavities are widely used to control microwave propagation. One of the main concerns is visible light energy is quickly dissipated within the metallic components, which makes this method of optical control almost impossible to generalize. Recently, an unconventional superconductivity in magic-angle graphene superlattices had been discovered and studied [36]. The superconductivity might help realize the metallic waveguide confinement of optics in the near future.
In this work, light propagation in 1D PCs enclosed in a rectangular waveguide or quasi-1D PCs is investigated in order to achieve a complete PBG while avoiding the difficulty in fabricating 3D PCs. This work complements two previous articles [15,16] that quantitatively analyzed omnidirectional light propagation in 1D and 2D PCs, respectively, both showing that a complete PBG cannot exist if an evanescent wave propagation is involved. The transfer matrix method is extended to study the transmittance of the quasi-1D PCs for both TE and TM polarization modes [34]. The corresponding BS are obtained by solving the eigenvalue equations with proper periodic boundary conditions following the Bloch theorem [3,4]. The formulas for evaluating the PDOS of the quasi-1D PCs for TE and TM modes are derived, respectively, for determining the PBGs. The contributions of the PDOS from each modes can be distinguished. The model is formulated in Section 2. The calculated results and discussion are presented in Section 3. The conclusions are given in Section 4.

2. Formulations

A transfer matrix approach is employed to discretize the dielectric function profile of the dielectric multilayer heterostructures and the transmission functions are calculated by matching the boundary conditions at each interfaces. In order to solve the PDOS, it is necessary to calculate the dispersion relation, and the corresponding band structures are obtained by solving the eigenvalue equations with proper periodic boundary conditions.

2.1. Transfer Matrix Method

Let us consider a waveguide with its rectangular cross section of sides a and b, and the enclosed multilayer dielectric slab with thickness, ( t 1 , t 2 , t 1 , t 2 , ) and dielectric function, ( ε 1 , ε 2 , ε 1 , ε 2 , ) , as shown in Figure 1. The TE mode (H mode) and TM mode (E mode) in the rectangular waveguide are characterized by the z components of the magnetic field and the electric field, H z and E z , respectively. By definition, these components are never absent in the corresponding modes. The z components of the Helmholtz’s equations for the inhomogeneous media are
ε ( z ) × 1 ε ( z ) × H z + ω 2 ε ( z ) μ ( z ) H z = 0
and
μ ( z ) × 1 μ ( z ) × E z + ω 2 ε ( z ) μ ( z ) E z = 0 .
In these cases, the effect of losses of the medium inside the waveguide is characterized by the complex permittivity ε ( z ) and permeability μ ( z ) . Thus, Equations (1) and (2) can be rearranged as:
2 x 2 + 2 y 2 + ε ( z ) z 1 ε ( z ) z H z + ω 2 ε ( z ) μ ( z ) H z = 0
and
2 x 2 + 2 y 2 + μ ( z ) z 1 μ ( z ) z E z + ω 2 ε ( z ) μ ( z ) E z = 0 .
By symmetry, using separation of variables, Equations (3) and (4) can be split into transverse and longitudinal parts, and the problem can be simplified as solving the one-dimensional Helmholtz’s equation along the z direction, the longitudinal parts,
ε ( z ) z 1 ε ( z ) z ψ ( z ) T E + ω 2 ε ( z ) μ ( z ) k c 2 ψ ( z ) T E = 0
and
μ ( z ) z 1 μ ( z ) z ψ ( z ) T M + ω 2 ε ( z ) μ ( z ) k c 2 ψ ( z ) T M = 0 ,
with eigenvalues k c which are determined by the following eigenvalue equations, the transverse parts,
2 x 2 + 2 y 2 + k c 2 ϕ ( x , y ) T E = 0
and
2 x 2 + 2 y 2 + k c 2 ϕ ( x , y ) T M = 0 .
The corresponding boundary conditions for ϕ ( x , y ) T E and ϕ ( x , y ) T M are
x ϕ ( x , y ) T E | x = 0 , a = 0 a n d y ϕ ( x , y ) T E | y = 0 , b = 0 ,
and
ϕ ( x , y ) T M | x = 0 , a = 0 a n d ϕ ( x , y ) T M | y = 0 , b = 0 .
It then follows that
k c = 2 π λ c = m π a 2 + n π b 2 .
Applying the eigenvalue equations, Equations (7) and (8), and the boundary conditions, Equations (9) and (10), to the general solutions of Equations (1) and (2), the following particular solutions can be found by separation of variables,
H z ( x , y , z ) = H 0 cos m π x a cos n π y b ψ ( z ) T E
and
E z ( x , y , z ) = E 0 sin m π x a sin n π y b ψ ( z ) T M ,
where H 0 and E 0 are determined by the energy of electromagnetic waves propagating inside the waveguide, and m and n are integers. The function ψ ( z ) is chosen for a particular solution since it represents propagating waves in the z-direction. A transfer matrix approach is employed to discretize the dielectric function profile of the heterostructure. For an N-layer dielectric-filled waveguide, ε ( z ) and μ ( z ) can be divided into p = 1 , 2 , , ( N + 2 ) layers with a piecewise constant permittivity ε p and constant permeability μ p , respectively. The discretized one-dimensional Helmholtz’s equation, for the pth region with constant permittivity ε p and constant permeability μ p can be written as
d 2 d z 2 ψ p ( z ) + k p 2 ψ p ( z ) = 0 for z p 1 z z p
with
k p = ω 2 ε p μ p k c 2 = 2 π λ ε r p μ r p λ λ c 2 ,
where ψ p ( z ) represents the wave function in the pth layer, and k p defines the complex wavevector in the same layer along the z-direction, λ represents the wavelength in free space at the operating angular frequency ω , ε r p and μ r p are the relative dielectric constant and permeability of the medium, respectively, and λ c is the cutoff wavelength. The solutions of Equations (12) and (13) can be written as a superposition of the forward and backward traveling wave functions:
ψ p = a p exp ( j k p z ) + b p exp ( j k p z ) f o r z p 1 z z p .
The boundary conditions for ψ ( z ) at the interface between layers p and ( p + 1 ) at position z = z p where p = 1 , 2 , , ( N + 1 ) are (for TE modes)
μ p ψ p ( z p ) = μ p + 1 ψ p + 1 ( z p ) a n d d d z ψ p ( z p ) = d d z ψ p + 1 ( z p ) ,
and (for TM modes)
ε p ψ p ( z p ) = ε p + 1 ψ p + 1 ( z p ) a n d d d z ψ p ( z p ) = d d z ψ p + 1 ( z p ) .
By matching the boundary conditions at each discontinuity, we arrive at
a N + 2 b N + 2 = M N + 1 · · · M p · · · M 1 a 1 b 1 = M 11 M 12 M 21 M 12 a 1 b 1 ,
where for TE modes:
M p = 1 2 exp ( j k p + 1 z p ) 0 0 exp ( j k p + 1 z p ) · μ p μ p + 1 + k p k p + 1 μ p μ p + 1 k p k p + 1 μ p μ p + 1 k p k p + 1 μ p μ p + 1 + k p k p + 1 · exp ( j k p + 1 z p ) 0 0 exp ( j k p + 1 z p ) ,
and for TM modes:
M p = 1 2 exp ( j k p + 1 z p ) 0 0 exp ( j k p + 1 z p ) · ε p ε p + 1 + k p k p + 1 ε p ε p + 1 k p k p + 1 ε p ε p + 1 k p k p + 1 ε p ε p + 1 + k p k p + 1 · exp ( j k p + 1 z p ) 0 0 exp ( j k p + 1 z p ) .
Using Equation (19) with a 1 = 1 , b 1 = r , a N + 2 = t , and b N + 2 = 0 , the reflection and transmission amplitudes, r and t, can be obtained, respectively, by
r = M 21 M 22
and
t = M 11 · M 22 M 12 · M 21 M 22 .
The reflection and transmission coefficients, R and T, can be implicitly represented by S-parameters, S 11 ( d B ) and S 12 ( d B ) , as a function of the operating frequency, for TE modes
S 11 ( d B ) = 10 l o g M 21 M 22 2 a n d S 12 ( d B ) = 10 l o g μ 1 k N + 2 μ N + 2 k 1 M 11 · M 22 M 12 · M 21 M 22 2 ,
and for TM modes
S 11 ( d B ) = 10 l o g M 21 M 22 2 a n d S 12 ( d B ) = 10 l o g ε 1 k N + 2 ε N + 2 k 1 M 11 · M 22 M 12 · M 21 M 22 2 .

2.2. Band Structures

In this model, E z and H z both are periodic with period Λ . According to the Bloch theorem, the electric and magnetic fields in a periodic layered medium are E z ( z ) = E z ( z + Λ ) and H z ( z ) = H z ( z + Λ ) , respectively. The column vector of the Bloch wave satisfies the following eigenvalue equation for consistency
a 3 b 3 = 1 4 A T E / T M B T E / T M C T E / T M D T E / T M · a 1 b 1 = e j k B Λ a 1 b 1 .
For the TE mode:
A T E = e j ( k 1 t 1 k 2 t 2 ) μ 1 μ 2 k 1 k 2 μ 2 μ 1 k 2 k 1 + e j ( k 1 t 1 + k 2 t 2 ) μ 1 μ 2 + k 1 k 2 μ 2 μ 1 + k 2 k 1 ,
B T E = e j [ k 1 t 1 k 2 ( 2 t 1 + t 2 ) ] μ 1 μ 2 + k 1 k 2 μ 2 μ 1 k 2 k 1 + e j [ k 1 t 1 + k 2 ( 2 t 1 + t 2 ) ] μ 1 μ 2 k 1 k 2 μ 2 μ 1 + k 2 k 1 ,
C T E = B T E * ,
and
D T E = A T E * .
For the TM mode:
A T M = e j ( k 1 t 1 k 2 t 2 ) ε 1 ε 2 k 1 k 2 ε 2 ε 1 k 2 k 1 + e j ( k 1 t 1 + k 2 t 2 ) ε 1 ε 2 + k 1 k 2 ε 2 ε 1 + k 2 k 1 ,
B T M = e j [ k 1 t 1 k 2 ( 2 t 1 + t 2 ) ] ε 1 ε 2 + k 1 k 2 ε 2 ε 1 k 2 k 1 + e j [ k 1 t 1 + k 2 ( 2 t 1 + t 2 ) ] ε 1 ε 2 k 1 k 2 ε 2 ε 1 + k 2 k 1 ,
C T M = B T M * ,
and
D T M = A T M * .
The phase factor e j k B Λ is thus the eigenvalue and satisfies the secular equation
A T E / T M e j k B Λ B T E / T M C T E / T M D T E / T M e j k B Λ = 0 .
Finally, the dispersion relation for the Bloch wave function is
k B ( k z , ω ) = 1 t 1 + t 2 cos 1 [ 1 2 ( A T E / T M + D T E / T M ) ] .

2.3. Photon Density of States

The quasi-one-dimensional photonic crystal has been confined in the x y -plane, so the wave vectors k x and k y are determined according to the guiding modes. To perform the PDOS calculation, it is required to use the formal definition which is the number of available photon modes per unit frequency range. Then, we construct two frequencies, namely, ω ( k B ) = ω and ω ( k B ) = ω + Δ ω , where Δ ω is a small increment. We calculate the line therein, and divide it by the line segment occupied by a single mode. The differential line element in K space within the frequency range, is given by Δ L k = Δ k B . Now, Δ k B is defined as
Δ k B = Δ ω | k ω | .
Integrating over the frequency increment, we have that the total phase space line segment contributing to the frequency range ( ω , ω + d ω ) is
ω k d L k = ω k 1 | k ω | d ω ,
where we take the limit of infinitesimal increments. The number of modes within the range ( ω , ω + d ω ) is obtained by dividing the length calculated in Equation (38) by the line segment corresponding to one mode, 2 π / Λ in the phase space. This yields
d N ( ω ) = Λ 2 π ω k 1 | k ω | d ω D ( ω ) d ω .
Because ω is a function of k, we can write
k ω = d ω d k B z ^ .
For the TE mode:
k ω T E = ( t 1 + t 2 ) sin [ 0.5 k B ( t 1 + t 2 ) ] α 1 + α 2 + α 3 + α 4 + α 5 + α 6 z ^
and for the TM mode:
k ω T M = ( t 1 + t 2 ) sin [ 0.5 k B ( t 1 + t 2 ) ] β 1 + β 2 + β 3 + β 4 + β 5 + β 6 z ^ ,
where the functions, α i and β i , with i = 1 , 2 , , 6 , are listed in the Appendix A. Finally, the formula for evaluating the PDOS can be expressed as
D ( ω ) T E / T M = Λ 2 π ω k 1 | k ω T E / T M | d ω .

3. Results and Discussion

All of macroscopic electromagnetism, including the propagation of light in a photonic crystal, is governed by the four macroscopic Maxwell’s equations with no free charges or currents. One interesting feature of electromagnetism in PCs is that there is no fundamental length scale other than the assumption that the system is macroscopic [1]. Therefore, to study physical phenomena in PCs, one may scale a system from the optical frequency range to the microwave one and vice versa if suitable conditions are fulfilled. For the purposes of demonstration and easier verification by experimentalists, a WR28 (7.11 mm × 3.555 mm) rectangular waveguide, usually used for Ka-band millimeter waves, is chosen. The periodic dielectric heterostructure is arranged along the longitudinal ( z ) direction to form the quasi-1D PCs [34]. In the microwave or millimeter wave frequency ranges, the PC experiments are very popular and more affordable compared to those in the optical frequency ones. Nevertheless, in order to extend the results for general use, the data are normalized for solutions at all length scale. Consider the quasi-1D PCs with the arrangement shown in Figure 1, 15 double-layer stacks (n = 30) have been included in the waveguide for investigation. The transmittance is calculated using the transfer matrix approach mentioned above. In order to demonstrate the stop band and pass band, the transmittances, expressed as the S-parameters, S 12 , in dB, of the lowest TE and TM modes for three quasi-1D PCs have been calculated and plotted in Figure 2. For the three cases, the dielectric constants used are the same as ε 1 = 3.8 (quartz) and ε 2 = 1.0 (air), while the thicknesses are varied as ( t 1 , t 2 ) = (1.00, 3.30), (1.00, 3.60), and (1.00, 3.90) mm, corresponding to the filling ratios t 1 / Λ = 23.26 % , 21.74 % , and 20.41 % for the periods of the stacks Λ = 4.3 , 4.6 , and 4.9 mm, respectively. As one can see, the central frequencies of the PBGs for both the T E 10 and T M 11 modes, as shown in Figure 2a,b respectively, are shifted to lower values as the filling ratio decreases. As mentioned above, the frequency axes in the plots are normalized to the cutoff frequency of T E 10 for general use.
In order to look for a complete PBG for all possible modes in the quasi-1D PCs, one may want to enlarge a PBG for a specific mode of choice. For the following four cases, the dielectric constants used are ( ε 1 , ε 2 ) = (2.3, 1.0), (3.8, 1.0), (4.9, 1.0), and (11.4, 1.0), while the thicknesses are varied as ( t 1 , t 2 ) = (2.15, 2.15), (1.00, 3.30), (0.72, 3.58), and (0.27, 4.03) mm, corresponding to the filling ratios t 1 / Λ = 50.00 % , 23.26 % , 16.74 % , and 6.28 % for the fixed period of the stacks Λ = 4.3 mm, respectively. The dielectric constants of 2.3, 3.8, 4.9 and 11.4 used in the chosen microwave frequency range correspond to the dielectric materials, polyethylene, quartz, phenolic resin, and barium sulfate, respectively. As one can see, the calculated transmittances, S 12 (dB), of the four quasi-1D PCs with different 15 double-layer stacks for the T E 10 and T M 11 modes are plotted in Figure 3a,b, respectively. The width of PBGs for both the T E 10 and T M 11 modes in the quasi-1D PCs is widened to larger sizes as decreasing the filling ratio while increasing the dielectric contrast between ε 1 and ε 2 . One may notice that the configurations of these quasi-1D PCs have been specially arranged so that the central frequencies of the PBGs of the T E 10 mode are kept the same while tuning the PBG sizes. However, those of the PBGs of the T M 11 mode are not centralized. This indicates that a tremendous computational effort is still needed to find a complete PBG in spite of the approach with a closed form developed is very efficient for the quasi-1D PCs.
A better way to identify PBGs more accurately and efficiently is to perform the BS calculation and further compute the PDOS [15,16]. Figure 4 shows the band structures of (a) TE10 and (b) TM11 modes for the quasi-1D PC with the dielectric constants ε 1 = 3.8 , ε 2 = 1 and the thicknesses t 1 = 1.00 mm, t 2 = 3.30 mm, corresponding to the case presented in Figure 2 and Figure 3 in a black solid line with a filling ratio of 23.26 % . The PBGs of the TE10 mode show no overlap with those of TM11 mode, all marked in gray stripes in the figure. However, from Figure 4, one may notice that there are no photon states in the frequency ranges below the cutoff frequency of the T E 10 mode nor within the first PBG of T E 10 which is under the cutoff frequency of the T M 11 mode. Note that only the lowest TE amd TM modes are plotted in Figure 4. As there might be other modes involved within the frequency range of the PBG of interest, it is easier to identify a complete PBG from the PDOS plots compared to the BS ones. Figure 5 shows the PDOS of TE10, TE01, TE11, and TM11 modes for the same quasi-1D PC. As one can see, the PDOS of the TE10 and TM11 modes are consistent with the BS calculations shown in Figure 4. The PDOS contributed from the TE01 and TE11 modes tend to fill up the first PBG of the TE10 mode. Other higher order modes are not considered as their cutoff frequencies are too high to contribute any PDOS in the frequency range of the PBG. Finally, the combined PDOS of TE10, TE01, TE11, and TM11 modes shows no photon states in some frequency ranges. The first one is below the cutoff frequency of the T E 10 mode, (0–0.77)fc, the second one is within the first PBG of T E 10 but under the cutoff frequency of the T E 01 mode, (1.26–1.48)fc, and the third one is the overlap of the PBGs of TE10, TE01, and TE11 modes, (1.79–1.87)fc. Therefore, a “complete PBG” can be obtained for some frequency ranges and categorized into three types: (1) below the cutoff frequency of the fundamental TE mode, (2) within the PBG of the fundamental TE mode but below the cutoff frequency of the next higher order mode, and (3) within an overlap of the PBGs of either TE modes, TM modes, or both.

4. Conclusions

In summary, light propagation in quasi-1D PCs has been investigated quantitatively. The transmittances for both the TE and TM modes through a periodic multilayer heterostructure in a rectangular waveguide are calculated using the transfer matrix method. The corresponding band structures are obtained by solving the eigenvalue equations with proper periodic boundary conditions following the Bloch theorem. The formulas for determining the PDOS have been obtained to facilitate identifying the photonic band gaps for all the modes residing in the system. With our approach, the quantitative determination of PDOS can be performed very accurately and efficiently. A complete PBG can exist in these quasi-1D PCs, but the determination must be carefully conducted and verified. It is demonstrated that three types of “complete PBG” can be found and categorized as follows. The first type is the frequency range within which the TE and TM modes are both cutoff, the second type is for which the fundamental TE mode has a PBG while other higher order modes are cutoff, and the last type is an overlap of the PBGs of either TE modes, TM modes, or both. The model might be easier for an experimental validation in a millimeter wave frequency range while the optical counterpart might be possibly pursued as well. We believe these results are of general importance and relevance to the dipole radiation or spontaneous emission by an atom in quasi-1D periodic structures and may have applications in future photonic quantum technologies.

Author Contributions

Conceptualization, M.-C.L.; methodology, R.-F.J. and M.-C.L.; software, R.-F.J. and M.-C.L.; validation, R.-F.J. and M.-C.L.; formal analysis, R.-F.J.; investigation, R.-F.J. and M.-C.L.; resources, R.-F.J. and M.-C.L.; data curation, R.-F.J.; writing—original draft preparation, R.-F.J. and M.-C.L.; writing—review and editing, R.-F.J. and M.-C.L.; visualization, R.-F.J.; supervision, M.-C.L.; project administration, M.-C.L.; funding acquisition, R.-F.J. and M.-C.L.

Funding

This work was partially supported by Guangdong Industry Polytechnic, P. R. China, under Grant Nos. KYRC2018-001 and RC2016-005, the research fund of Hanyang University (HY- 201400000002393), National Research Foundation of South Korea (2015R1D1A1A01061017), and the Alexander von Humboldt Foundation of Germany.

Acknowledgments

The authors would like to thank the late B. Y. Gu at the Institute of Physics, CAS and C. T. Chan at the Department of Physics, HKUST for the helpful comments, discussions, and encouragement.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
1Done-dimensional
2Dtwo-dimensional
3Dthree-dimensional
BSband structure
EMelectromagnetic
PBGphotonic band gap
PCphotonic crystal
PDOSphoton density of states
TEtransverse electric
TMtransverse magnetic

Appendix A. FORMULAS

Here, we give the formulas of the functions employed in Equations (41) and (42).
α 1 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] exp ( j 2 k 1 t 1 ) + exp ( j 2 k 2 t 2 ) ω k 1 ε 1 μ 1 k 2 k 1 2 ε 2 μ 2 k 2 k 1 k 2 + μ 1 μ 2 ,
α 2 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] { 1 + exp [ j 2 ( k 1 t 1 + k 2 t 2 ) ] } ω k 1 ε 1 μ 1 k 2 k 1 2 + ε 2 μ 2 k 2 k 1 k 2 + μ 1 μ 2 ,
α 3 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] exp ( j 2 k 1 t 1 ) + exp ( j 2 k 2 t 2 ) ω k 2 ε 1 μ 1 k 1 + ε 2 μ 2 k 1 k 2 2 k 2 k 1 + μ 2 μ 1 ,
α 4 = j exp [ j ( k 1 t 1 + k 2 t 2 ) ] exp ( j 2 k 1 t 1 ) exp ( j 2 k 2 t 2 ) ω ε 1 μ 1 t 1 k 1 + ε 2 μ 2 t 2 k 2 k 1 k 2 + μ 1 μ 2 k 2 k 1 + μ 2 μ 1 ,
α 5 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] { 1 + exp [ j 2 ( k 1 t 1 + k 2 t 2 ) ] } ω k 2 ε 1 μ 1 k 1 ε 2 μ 2 k 1 k 2 2 k 2 k 1 + μ 2 μ 1 ,
α 6 = j exp [ j ( k 1 t 1 + k 2 t 2 ) ] { 1 + exp [ j 2 ( k 1 t 1 + k 2 t 2 ) ] } ω ε 1 μ 1 t 1 k 1 ε 2 μ 2 t 2 k 2 k 1 k 2 + μ 1 μ 2 k 2 k 1 + μ 2 μ 1 ,
β 1 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] exp ( j 2 k 1 t 1 ) + exp ( j 2 k 2 t 2 ) ω k 1 ε 1 μ 1 k 2 k 1 2 ε 2 μ 2 k 2 k 1 k 2 + ε 1 ε 2 ,
β 2 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] { 1 + exp [ j 2 ( k 1 t 1 + k 2 t 2 ) ] } ω k 1 ε 1 μ 1 k 2 k 1 2 + ε 2 μ 2 k 2 k 1 k 2 + ε 1 ε 2 ,
β 3 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] e x p ( j 2 k 1 t 1 ) + e x p ( j 2 k 2 t 2 ) ω k 2 ε 1 μ 1 k 1 + ε 2 μ 2 k 1 k 2 2 k 2 k 1 + ε 2 ε 1 ,
β 4 = j exp [ j ( k 1 t 1 + k 2 t 2 ) ] exp ( j 2 k 1 t 1 ) exp ( j 2 k 2 t 2 ) ω ε 1 μ 1 t 1 k 1 + ε 2 μ 2 t 2 k 2 k 1 k 2 + ε 1 ε 2 k 2 k 1 + ε 2 ε 1 ,
β 5 = exp [ j ( k 1 t 1 + k 2 t 2 ) ] { 1 + exp [ j 2 ( k 1 t 1 + k 2 t 2 ) ] } ω k 2 ε 1 μ 1 k 1 ε 2 μ 2 k 1 k 2 2 k 2 k 1 + ε 2 ε 1 ,
and
β 6 = j exp [ j ( k 1 t 1 + k 2 t 2 ) ] { 1 + exp [ j 2 ( k 1 t 1 + k 2 t 2 ) ] } ω ε 1 μ 1 t 1 k 1 ε 2 μ 2 t 2 k 2 k 1 k 2 + ε 1 ε 2 k 2 k 1 + ε 2 ε 1 .

References

  1. Joannopoulos, J.D.; Johnson, S.G.; Winn, J.N.; Meade, R.D. Photonic Crystals: Molding the Flow of Light; Princeton University Press: Princeton, NJ, USA, 2008. [Google Scholar]
  2. Joannopoulos, J.D.; Villeneuve, P.R.; Fan, S. Photonic crystals: Putting a new twist on light. Nature 1997, 386, 143–149. [Google Scholar] [CrossRef]
  3. Yariv, A.; Yeh, P. Optical Waves in Crystals; Wiley: New York, NY, USA, 1984. [Google Scholar]
  4. Kittel, C. Introduction to Solid State Physics; Wiley: New York, NY, USA, 1976. [Google Scholar]
  5. O’Brien, J.L.; Furusawa, A.; Vučković, J. Photonic quantum technologies. Nat. Photonics 2009, 3, 687–691. [Google Scholar] [CrossRef]
  6. Purcell, E.M. Spontaneous emission probabilities at radio frequencies. Phys. Rev. 1946, 69, 681. [Google Scholar]
  7. Kleppner, D. Inhibited Spontaneous Emission. Phys. Rev. Lett. 1981, 47, 233–236. [Google Scholar] [CrossRef]
  8. Barut, A.O.; Dowling, J.P. Quantum electrodynamics based on self-energy: Spontaneous emission in cavities. Phys. Rev. A 1987, 36, 649–654. [Google Scholar] [CrossRef] [PubMed]
  9. Rigneault, H.; Monneret, S. Modal analysis of spontaneous emission in a planar microcavity. Phys. Rev. A 1996, 54, 2356–2368. [Google Scholar] [CrossRef]
  10. Dowlin, J.P.; Bowden, C.M. Atomic emission rates in inhomogeneous media with applications to photonic band structures. Phys. Rev. A 1992, 46, 612–622. [Google Scholar] [CrossRef]
  11. Suzuki, T.; Yu, P.K.L. Emission power of an electric dipole in the photonic band structure of the fcc lattice. Opt. Soc. Am. B 1995, 12, 570–582. [Google Scholar] [CrossRef]
  12. Kamli, A.; Babiker, M.; Al-Hajry, A.; Enfati, N. Dipole relaxation in dispersive photonic band-gap structures. Phys. Rev. A 1997, 55, 1454–1461. [Google Scholar] [CrossRef]
  13. Sánchez, A.S.; Halevi, P. Spontaneous emission in one-dimensional photonic crystals. Phys. Rev. E 2005, 72, 056609. [Google Scholar] [CrossRef]
  14. Halevi, P.; Sánchez, A.S. Spontaneous emission in a high-contrast one-dimensional photonic crystal. Opt. Commun. 2005, 251, 109–114. [Google Scholar] [CrossRef]
  15. Lin, M.C.; Jao, R.F. Quantitative analysis of photon density of states for a realistic superlattice with omnidirectional light propagation. Phys. Rev. E 2006, 74, 046613. [Google Scholar] [CrossRef] [PubMed]
  16. Jao, R.F.; Lin, M.C. Efficient and quantitative analysis of photon density of states for two-dimensional photonic crystals with omnidirectional light propagation. Phys. Rev. E 2018, 98, 053306. [Google Scholar] [CrossRef] [Green Version]
  17. Rumsey, I.; Melinda, P.M.; Kelly, P.K. Photonic bandgap structures used as filters in microstrip circuits. IEEE Microw. Guided Wave Lett. 1998, 8, 336–338. [Google Scholar] [CrossRef]
  18. Kim, T.; Seo, C. A novel photonic bandgap structure for low-pass filter of wide stopband. IEEE Microw. Guided Wave Lett. 2000, 10, 13–15. [Google Scholar]
  19. Iluz, Z.; Shavit, R.; Bauer, R. Microstrip antenna phased array with electromagnetic bandgap substrate. IEEE Trans. Antennas Prop. 2004, 52, 1446–1453. [Google Scholar] [CrossRef]
  20. Yang, F.R.; Qian, Y.; Coccioli, R.; Itoh, T. A novel low slow-wave microstrip structure. IEEE Microw. Guided Wave Lett. 1998, 8, 372–374. [Google Scholar] [CrossRef]
  21. Xue, Q.; Shum, K.; Chan, C. Novel 1D microstrip PBG cells. IEEE Microw. Wirel. Compon. Lett. 2000, 10, 403–405. [Google Scholar]
  22. Yang, H.W.; Zhong, W.X.; Sui, Y.K. Analysis of dielectric layer PBG structure using precise integration. In Proceedings of the International Symposium on Computational Mechanics, Beijing, China, 30 July–1 August 2007; pp. 1239–1245. [Google Scholar]
  23. Kshetrimayum, R.S.; Zhu, L. Guided-wave characteristics of waveguide based periodic structures loaded with various FSS strip layers. IEEE Trans. Antennas Prop. 2005, 53, 120–124. [Google Scholar] [CrossRef]
  24. Marini, S.; Covers, Á.; Boria, V.E.; Gimeno, B. Full-wave modal analysis of slow-wave periodic structures loaded with elliptical waveguides. IEEE Trans. Electron Devices 2010, 57, 516–524. [Google Scholar] [CrossRef]
  25. Razavizadeh, S.M.; Sadeghzadeh, R.; Navarri-Cía, M.; Ghattan, Z. Compact THz waveguide filter based on periodic dielectric-gold rings. In Proceedings of the 5th International Conference on Millimeter-Wave and Terahertz Technologies, Tehran, Iran, 18–20 December 2018; pp. 42–44. [Google Scholar]
  26. Pelster, R.; Gasparian, V.; Nimtz, G. Propagation of plane waves and of waveguide modes in quasiperiodic dielectric heterostructures. Phys. Rev. E 1997, 55, 7645–7655. [Google Scholar] [CrossRef] [Green Version]
  27. Amari, S.; Vahldieck, R.; Bornemann, J.; Leuchtmann, P. Propagation in a circular waveguide periodically loaded with dielectric disks. IEEE MTT-S Digest. 1998, 3, 1535–1538. [Google Scholar]
  28. Kesari, V.; Basu, B.N. Analysis of some periodic structures of microwave tubes: Part II: Analysis of disc-loaded fast-wave circular waveguide structures for gyro-travelling-wave tubes. J. Electomagn. Waves Appl. 2018, 32, 1–36. [Google Scholar] [CrossRef]
  29. Christie, L.; Erabati, G.; Jana, M. Analysis of propagation characteristics of circular waveguide loaded with dielectric disks using coupled integral equation technique. In Proceedings of the 5th International Conference on Advances in Computing and Communications, Kochi, India, 2–4 September 2015; pp. 231–234. [Google Scholar]
  30. Miller, R.D.; Jones, T.B. On the effective dielectric constant of columns or layers of dielectric spheres. J. Phys. D Appl. Phys. 1988, 21, 527–532. [Google Scholar] [CrossRef]
  31. Peng, S.T.; Tamir, T.; Bertoni, H.L. Theory of periodic dielectric waveguides. IEEE Trans. Microw. Theory Technol. 1975, MTT-23, 123–133. [Google Scholar] [CrossRef]
  32. Villeneuve, P.R.; Fan, S.; Johnson, S.G.; Joannopoulos, J.D. Three-dimensional photonic confinement in photonic crystals of low-dimensional periodicity. IEE Proc. Optoelectron. 1998, 145, 384–390. [Google Scholar] [CrossRef]
  33. Young, A.B.; Thijssen, A.C.T.; Beggs, D.M.; Androvitsaneas, P.; Kuipers, L.; Rarity, J.G.; Hughes, S.; Oulton, R. Polarization engineering in photonic crystal waveguides for spin-photon entanlers. Phys. Rev. Lett. 2015, 115, 153901. [Google Scholar] [CrossRef] [PubMed]
  34. Lin, M.C. A multilayer waveguide window for wide-bandwidth millimeter wave tubes. Int. J. Infrared Millim. Waves 2007, 28, 355–362. [Google Scholar] [CrossRef]
  35. Yablonovitch, E. Inhibited spontaneous emission in solid-state physics and electronics. Phys. Rev. Lett. 1987, 58, 2059–2062. [Google Scholar] [CrossRef]
  36. Cao, Y.; Fatemi, V.; Fang, S.; Watanabe, K.; Taniguchi, T.; Kaxiras, E.; Jarillo-Herrero, P. Unconventional superconductivity in magic-angle graphene superlattices. Nature 2018, 556, 43–50. [Google Scholar] [CrossRef]
Figure 1. (a) 3D schematic of periodic multilayer heterostructure along the z-direction confined in a rectangular waveguide with a width a and a height b and (b) the corresponding dielectric function profile ε ( z ) in the cross-sectional view.
Figure 1. (a) 3D schematic of periodic multilayer heterostructure along the z-direction confined in a rectangular waveguide with a width a and a height b and (b) the corresponding dielectric function profile ε ( z ) in the cross-sectional view.
Crystals 09 00576 g001
Figure 2. Calculated transmittance, S 12 (dB), of the quasi-1D PCs with 15 double-layer stacks enclosed in the rectangular waveguide for (a) T E 10 and (b) T M 11 modes of light propagation. The dielectric constants used are ε 1 = 3.8 and ε 2 = 1.0 , and the corresponding filling ratios are t 1 / Λ = 23.26 % , 21.74 % , and 20.41 % in the three cases. The frequency axes are normalized to the cutoff frequency of T E 10 .
Figure 2. Calculated transmittance, S 12 (dB), of the quasi-1D PCs with 15 double-layer stacks enclosed in the rectangular waveguide for (a) T E 10 and (b) T M 11 modes of light propagation. The dielectric constants used are ε 1 = 3.8 and ε 2 = 1.0 , and the corresponding filling ratios are t 1 / Λ = 23.26 % , 21.74 % , and 20.41 % in the three cases. The frequency axes are normalized to the cutoff frequency of T E 10 .
Crystals 09 00576 g002
Figure 3. Calculated transmittance, S 12 (dB), of the quasi-1D PCs with 15 double-layer stacks enclosed in the rectangular waveguide for (a) T E 10 and (b) T M 11 modes of light propagation. The dielectric constants used are ε 1 = 2.3, 3.8, 4.9, and 11.4 and ε 2 = 1 , and the corresponding filling ratios are t 1 / Λ = 50.00 % , 23.26 % , 16.74 % , and 6.28 % in the four cases. The frequency axes are normalized to the cutoff frequency of T E 10 .
Figure 3. Calculated transmittance, S 12 (dB), of the quasi-1D PCs with 15 double-layer stacks enclosed in the rectangular waveguide for (a) T E 10 and (b) T M 11 modes of light propagation. The dielectric constants used are ε 1 = 2.3, 3.8, 4.9, and 11.4 and ε 2 = 1 , and the corresponding filling ratios are t 1 / Λ = 50.00 % , 23.26 % , 16.74 % , and 6.28 % in the four cases. The frequency axes are normalized to the cutoff frequency of T E 10 .
Crystals 09 00576 g003
Figure 4. Band structures of (a) TE10 and (b) TM11 modes for the quasi-1D PC with the dielectric constants ε 1 = 3.8 , ε 2 = 1 and the thicknesses t 1 = 1.00 mm, t 2 = 3.30 mm, corresponding to a filling ratio of 23.26 % . The PBGs of TE10 mode show no overlap with those of TM11 mode.
Figure 4. Band structures of (a) TE10 and (b) TM11 modes for the quasi-1D PC with the dielectric constants ε 1 = 3.8 , ε 2 = 1 and the thicknesses t 1 = 1.00 mm, t 2 = 3.30 mm, corresponding to a filling ratio of 23.26 % . The PBGs of TE10 mode show no overlap with those of TM11 mode.
Crystals 09 00576 g004
Figure 5. PDOS of TE10, TE01, TE11, and TM11 modes for the quasi-1D PC with the dielectric constants ε 1 = 3.8 , ε 2 = 1 and the thicknesses t 1 = 1.00 mm, t 2 = 3.30 mm, corresponding to a filling ratio of 23.26 % . The combined PDOS of TE10, TE01, TE11, and TM11 modes shows no photon states in some frequency ranges.
Figure 5. PDOS of TE10, TE01, TE11, and TM11 modes for the quasi-1D PC with the dielectric constants ε 1 = 3.8 , ε 2 = 1 and the thicknesses t 1 = 1.00 mm, t 2 = 3.30 mm, corresponding to a filling ratio of 23.26 % . The combined PDOS of TE10, TE01, TE11, and TM11 modes shows no photon states in some frequency ranges.
Crystals 09 00576 g005

Share and Cite

MDPI and ACS Style

Jao, R.-F.; Lin, M.-C. Quantitative Analysis of Photon Density of States for One-Dimensional Photonic Crystals in a Rectangular Waveguide. Crystals 2019, 9, 576. https://doi.org/10.3390/cryst9110576

AMA Style

Jao R-F, Lin M-C. Quantitative Analysis of Photon Density of States for One-Dimensional Photonic Crystals in a Rectangular Waveguide. Crystals. 2019; 9(11):576. https://doi.org/10.3390/cryst9110576

Chicago/Turabian Style

Jao, Ruei-Fu, and Ming-Chieh Lin. 2019. "Quantitative Analysis of Photon Density of States for One-Dimensional Photonic Crystals in a Rectangular Waveguide" Crystals 9, no. 11: 576. https://doi.org/10.3390/cryst9110576

APA Style

Jao, R. -F., & Lin, M. -C. (2019). Quantitative Analysis of Photon Density of States for One-Dimensional Photonic Crystals in a Rectangular Waveguide. Crystals, 9(11), 576. https://doi.org/10.3390/cryst9110576

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop