Coupling of Photonic Crystal Surface Modes

Yücel, Melike Behiye

doi:https://doi.org/10.1155/2022/8947410

Advances in Condensed Matter Physics

On this page

Abstract Introduction Conclusion Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 8947410 | https://doi.org/10.1155/2022/8947410

Coupling of Photonic Crystal Surface Modes

Melike Behiye Yücel¹

Academic Editor: Ashwini Kumar

Received07 Mar 2022

Accepted06 Jun 2022

Published10 Aug 2022

Abstract

Guiding and evanescent coupling properties of surface modes bound to the interfaces of two-dimensional photonic crystals in close proximity are numerically demonstrated. Interacting photonic crystals are composed of silicon pillars in air, where their outermost layers facing each other are annular. Surface modes are identified through supercell band structure computations, while their excitation by the electromagnetic waves through a perpendicular insertion waveguide is demonstrated using finite-difference time-domain simulations. Lifting the degeneracy between the surface modes as a consequence of bringing two identical photonic crystal surfaces to a sufficient distance results in evanescent coupling in a beating manner whose beat length linearly varies between 10 and 20 periods up to a frequency at which both surface modes travel with the same group velocity. The surface mode coupling phenomenon could be employed either to enhance sensitivity or to reduce device size in bio/chemical sensor applications since the effective travelling length of surface waves increases by about 3.5 times due to evanescent coupling.

1. Introduction

Photonic crystals (PCs) are structures whose refractive indices are periodically modulated. They possess numerous interesting properties, which have been the subject of scientific and technological attention ever since they were first introduced [1]. PCs display photonic band gap (PBG) similar to their electronic counterparts in which no electromagnetic (EM) waves could be transmitted. EM wave transmission for frequencies lying in the PBG of a PC can be achieved by introducing, for instance, linear defects to form waveguides (WGs) in which EM waves are confined in and decay abruptly in the transverse direction [2–5]. Linear defect WGs have found widespread utility in applications such as wavelength division de/multiplexing [6, 7] and unidirectional light transmission [8, 9]. EM waves in two parallel WGs were shown to couple and decouple and thus enable EM power transfer in between if the dielectric region between the WGs is properly designed [10, 11]. Coupling efficiencies of such PC directional couplers have been improved [12, 13] and successfully adopted in applications such as all-optical switching [14], wavelength selective optical filters [12, 15], power splitters [16, 17], and polarizing beam splitters [18].

Yet, another means for EM wave confinement and guiding for frequencies lying in the PBG are realized at the interface between a PC and a homogenous medium in which the PC is embedded [19–26]. EM waves propagating along the PC surface they bound to, namely surface waves (SWs) [19–25], have been extensively studied. SW bands within the PBG should lie below the light line of the homogenous matrix so that SWs do not penetrate either region. Surface modes excited along a corrugated PC surface at the output of a perpendicular insertion WG (IWG) were shown to focus the outgoing EM waves by suppressing diffraction [26, 27] as the modes created at the input surface enhance transmission [28]. The crucial goal of efficiently coupling PC surface and line defect WG modes for all-PC optical circuits is accomplished by modifying the WGs and interfaces [20, 29]. SWs were also employed in achieving directional emission out of a WG in a PC [20]. A 1-to-N beam splitter, which is an important device in all-PC optical circuits, was realized by adjusting the positions of point defects on a modulated surface [30].

Microcavity configurations on PC surfaces are usually employed to utilize SWs in lasing and sensing applications [21–23, 31–35] since the surface modes within the semi-infinite crystal surface in which the translational symmetry along the PC boundary is broken turn into resonant states. Refractive index sensitivities ranging from to , where stands for refractive index unit, are reported for the systems employing microcavity configurations [33–35]. Moreover, dispersion properties of the SWs are highly sensitive to refractive index variations, which could take place as a result of variations in chemical composition and/or concentration, in the homogenous liquid or gas matrix in which the PC is immersed [20–23, 36]. Such properties of SWs have been in use for sensor applications in which sensitivity values of and are achieved on the surfaces of square and triangular PC lattices, respectively [37].

In addition to the studies mentioned above, there are also studies in which different designs are used in sensor applications [38–40]. Wellenzohn and co-workers proposed and designed a two-dimensional photonic crystal defect waveguide biosensor based on complementary metal oxide semiconductor (CMOS)-compatible silicon-on-insulator technology operating in aqueous solutions at 1.34 μm [38]. The authors compared this operation wavelength with the 1.55 μm and reported a significantly smaller propagation loss at the wavelength they studied. On the other hand, Liu et al. proposed a configuration of D-shaped five-hole photonic crystal fiber-based surface plasmon resonance refractive index sensor with a maximum sensitivity of and a maximum detection range of 1.30 to 1.50 [39]. In a study modeling a highly sensitive sensor based on a slot photonic crystal waveguide facilitating the enhancement of the interaction between the analyte and light for highly sensitive evanescent field absorption sensing in fluids, the symmetry of the photonic crystal waveguide modes with and without a slot has been investigated [40]. As a result, it has been shown that there is an enhancement factor of 7.6 in the evanescent field ratio compared with a slab waveguide with a similar thickness by tuning the slot width.

In this article, guiding and evanescent coupling properties of surface modes bound to the interfaces of two semi-infinite PCs with homogenous host when they are brought to close proximity are numerically demonstrated and possible uses of this phenomenon in sensor applications are discussed. Surface modes are identified through supercell band structure (BS) computations, while their excitation by the waves travelling along a perpendicular IWG is demonstrated using two-dimensional (2D) finite-difference time-domain (FDTD) method [41].

The study is organized as follows: Band structures of the PCs, the IWG, a single PC surface, and two interacting PC surfaces are presented in Section 2. Optimization of surface modification to achieve broadband coupling, the resultant modal properties of coupling surface modes with definite parities, and spectral variation of coupling length and group velocities of coupling modes are also discussed in this section. Demonstration of the coupling between surface modes through FDTD simulations is presented and discussed in Section 3. Brief conclusions are provided inSection 4.

2. Design and Dispersion Properties of Photonic Crystal Surfaces

The underlying 2D PC is composed of a square array of silicon pillars in air, whose refractive indices are and , respectively. The PC is laid along the [10] direction (Figure 1). To facilitate surface mode coupling at wavelengths around , lattice constant and radii of the cylinders are set to and , respectively. Such high-aspect ratio pillar-type devices can easily be fabricated through techniques, such as deep reactive ion etching [42–44].

BS of the square PC obtained through the plane-wave expansion (PWE) method [45, 46] with 2¹⁶ () plane waves implemented in RSoft Design Groups BandSOLVE software reveals a complete PBG existing along the direction for a range of wave vector k increased by 0.001 for transverse magnetic (TM) polarized waves, i.e., the electric field being normal to the PC (xy) plane between normalized angular frequencies and . Here, is the angular frequency and c is the speed of light in free space. Corresponding gap width and gap-over-midgap ratio values are and , respectively. Since it is a common practice to couple EM waves into the surface modes through perpendicular IWGs [47], a vertical WG on the left of Figure 1 denoted as IWG is introduced. As it will be discussed in detail below, its width, d_IWG, is chosen such that a sufficiently broad air-like band, which covers the surface bands of interest, can be obtained. TM polarized EM wave is launched from the lower end, while the upper end is left open to prevent backscattering. Calculations of the IWG bands, as well as single/coupled surface bands, are carried out within a supercell approach [48] through the PWE method.

Localized Bloch modes on PC surfaces are obtained by replacing the scatterers at the outmost layer by annuli [37] with inner and outer radii of and , respectively (Figure 1). The center-to-center distance between the annuli on the upper and lower surfaces, , is optimized such that it is not (i) too small to accommodate a single WG band rather than surface bands, or (ii) too large to prevent mode overlap, and thus evanescent coupling. The first 9 elements of the upper PC to the immediate right of the IWG, denoted as the transition region in Figure 1, are not modified such that no surface modes are supported on the upper surface within the transition region to ensure that initially surface modes on only the lower PC surface are excited by the IWG. Thus, only the TM polarized waves from the lower surface can couple to the modes of the upper surface after the transition region.

Projected BSs of the IWG for two different d_IWG values, i.e., and , accompanied by those of a single PC surface interfacing with air host where the surface termination is annular, as well as BSs for two identical such PC surfaces in close proximity, are presented in Figure 2. These BSs are obtained by varying the respective parallel wave vector modulus, , to surface. Figure 2 also presents the projection of bulk bands over the ΓX direction and the light line of air host, i.e., the line. In case of the IWG, 9 rows of scatterers on each side, accompanied by an air gap, are employed. For and , this corresponds to a and supercell, respectively. A single PC supercell with annular termination involves 11 rows of PC elements, 1 row of annular termination, and a thick air layer, resulting in a supercell. Supercell for coupled surfaces involves 6 rows of PC elements terminated by 1 row of annulus on each side and a wide gap, giving rise to a supercell. Geometries of all supercells are shown in Figure 1, although not drawn to scale for compactness. Choice of such supercells ensures sufficient decay of the WG or surface modes into the PC bulk, sufficient decay of the surface modes into air in the single surface case, as well convergence of the obtained bands. This will be discussed in detail below in terms of mode profiles.

Figure 2

(Color online) projected BSs of the IWG for (light blue-colored dashed line with hollow circles) and (purple-colored dashed lines with filled circles), as well as a single PC surface terminated with annular scatterers (magenta-colored dash-dotted lines) and two identical surfaces in proximity (solid lines with triangles where the red-colored filled and the pink-colored empty triangles represent the surface and WG bands, respectively). Inside the PBG, WG bands lie above the light line, while surface bands lie below the light line. The solid shaded areas are the bulk bands, whereas the graded rectangular area represents the full TM PBG. The horizontal thick dash-dotted line (black-colored) indicates the wavelength of interest, i.e., .

Figure 2 shows that all IWG modes lie above the light line for either narrow () or wide () WG width. Only one band for , which covers a small frequency range in the upper portion of the band gap much above the wavelength of interest, i.e., , is observed. In contrast, there are three bands for , one of which is close to the light line and varies in close resemblance with it. In fact, such a broad air-like band spanning the whole band gap can be obtained for covering , as depicted by the horizontal dash-dotted line in Figure 2. Thus, modes of this band can be utilized in exciting surface modes of the adjacent perpendicular surface in Figure 1.

To optimize the single PC surface to cover a large bandwidth around and to lie below the light line, the inner and outer radii of the annular elements on the surface, and , respectively, are treated as optimization parameters. For ease of calculations, is fixed to , while is varied in the range where the case corresponds to a hollow WG. A sufficiently broad surface band (SB), which covers and lies below the light line, can be obtained for . This band can be utilized to obtain surface modes in a frequency range between the lower band gap edges, i.e., and , with a bandwidth of around the central frequency. Air-like band of the IWG covers the whole frequency range of the SB, thus rendering the excitation of each mode within the SB.

Surface modes for adjacent PCs with annular surface elements are investigated using the supercell denoted as the coupling region in Figure 1. Here, the only relevant control parameter is restricted to and it determines the degree of mode overlap between the neighboring PC surfaces. A too large value gives rise to negligible mode overlap between neighboring PC surfaces, whereas a small value results in a significant overlap and thus the WG modes instead of coupling surface modes as discussed above. Figure 2 shows that the SB of the single PC surface lying between and (dash-dotted curve) splits into two SBs (solid lines accompanied by filled triangles) for a moderate of . One of the SBs is closer to the dielectric bands of the PC, whereas the other is to the light line. Although another band lying between and (solid line accompanied by empty triangles) in Figure 2 exists, it is out of scope, since it is a WG band whose modes are air-guided rather than being surface modes.

The discrepancy between the two-split SBs of the neighboring surfaces is primarily determined by . When is much larger than , SBs of the neighboring surfaces collapse into a single degenerate SB of the isolated PC surface. When , however, is sufficiently small, e.g., , degeneracy is lifted and the two-split SBs cover a common frequency range between the top of the bulk dielectric bands, i.e., , and the top of the lower SB at for reasons to be discussed below. The overlap region covers a frequency range of , or equivalently a bandwidth of , which significantly diminishes to as increases to . Moreover, the corresponding WG band depicted by the solid line with hollow triangle in Figure 2 starts overlapping the SBs for . Thus, is the optimal value for SB coupling.

Modal properties of the single PC surface band where and , along with those of the SBs of the adjacent PC surfaces for at , are presented in Figure 3. It is clearly seen that EM energy is concentrated mainly inside the annular region of the isolated PC surface (Figure 3(a)). The field strength rapidly drops toward both the air and PC bulk. Since the perpendicular electric field component, , is negligible inside the PC bulk, as demonstrated by the plot of the cutline along the vertical bisector of the supercell in Figure 3(a), the adopted supercell size provides a well-converged SB in Figure 2.

(a)

(b)

The mode profiles in Figure 3(b) demonstrate that the surface mode splits into an antisymmetric mode and a symmetric mode represented by filled down and up triangles, respectively, in Figure 2. Antisymmetric mode is such that there is a node midway between the neighboring surfaces, whereas there is significant overlap for the symmetric mode (Figure 3(b). This is remarkably similar to the case of coupling two parallel PC WGs [49]. This resemblance can be utilized in obtaining unidirectional light transmission through coupling surface modes. Inspection of the surface modes at a given value reveals that antisymmetric modes are always at higher angular frequencies.

An incident TM polarized wave through the IWG in Figure 1 first couples to the surface mode in the transition region. Thereafter, it couples back and forth between the antisymmetric and symmetric modes until the exit on the right side of the system in Figure 1. A quantitative measure of evanescent coupling between the surface modes is the normalized spatial coupling period, denoted as , where is the coupling length determined by the difference between the values of the modes, , as [10, 49] follows:

Another useful quantity is the beat length, which denotes the full spatial cycle for the energy to couple back and forth once, defined as .

Figure 4(a) demonstrates that varies roughly between 5 and 10 for , and . It first increases monotonically for and then decreases rapidly. at is remarkably short and can lead to design of very compact devices out of the coupling mechanism. varies linearly for in which a linear fit, , reveals that the rate of change with respect to is , or equivalently per THz of frequency shift, or per percent variation around . This is a measurable shift, and tracking the variation of can be used for sensing purposes.

(a)

(b)

A better approach for evaluating the possibility of using surface mode coupling for sensing is considering the first derivative of with respect to , i.e., , the result of which is also presented in Figure 4(a). Although it takes small negative values for frequencies below (equivalently ), it sharply increases for higher frequencies close to the symmetric SB top, or . As also drops sharply around the band edge, the surface mode coupling mechanism can be used in sensing devices.

The coupling behavior of symmetric and antisymmetric surface modes is closely related to group velocities of the modes, which can be determined through the dispersion curves in Figure 2. As the slopes of the two-split surface bands in Figure 2 are generally different, modes with different parity propagate along PC surfaces with different speeds. This, in turn, results in a phase difference between coupling symmetric and antisymmetric modes. Figure 4(b) reveals that (normalized) group velocity of the symmetric mode is higher than that of the antisymmetric mode up to a frequency at which they become equal. For frequencies higher than , the case is reversed (Figure 4(b)). No coupling is expected for where only the antisymmetric mode is excited up to the band top. Thus, three frequencies, i.e., , , and , are chosen for FDTD simulations, which will be discussed in the next section to demonstrate out-of-phase, in-phase coupling, and single-mode excitation, respectively.

3. Demonstration of Surface Mode Coupling

Dynamic visualization of surface mode coupling is carried out through FDTD simulations [41] by RSoft Design Group’s FullWAVE software. Grid size in FDTD computations is chosen as , while the time step is such that to satisfy the Courant–Friedrichs–Lewy stability criterion [50]. The computational domain is surrounded by perfectly matched layer-absorbing boundaries [51]. Surface mode coupling is demonstrated through probing the electric field distribution over the linear patches denoted as and in Figure 1.

FDTD simulation of surface mode coupling at is presented in Figure 5. The TM polarized wave is launched from the bottom of the IWG couples evanescently into the lower surface in the transition region. When the surface wave reaches the coupling region, it couples back and forth between the upper and lower surfaces. The coupling region is such that 2 complete beat cycles are covered (Figure 5).

Figure 5

(Color online) FDTD simulation of surface mode coupling at (top panel). The middle panel is a close-up view of the transition and coupling regions, whereas the bottom panel presents the variation of Poynting vector magnitude along the horizontal dash-dotted cutlines. The plot to right on the middle panel depicts the variation of on the vertical dash-dotted cutline at the output. The thick hollow arrow in the top panel indicates the initially launched wave direction, whereas the shaded area in the bottom panel denotes the coupling region.

The evanescent nature of SWs excited on the lower PC surface in the transition region in Figure 5 reveals that decay ratios are different on the two sides of the lower PC surface, due to significantly different media in either directions, e.g., PC bulk and air. The decay length of the SW in air is considerably larger than that into the PC bulk, as evidenced by the mode profile in Figure 3(a). This, in turn, suggests that SWs on the lower PC could well excite corresponding modes on the upper PC while they travel along the x-axis.

For visual purposes, the area within the dashed rectangle in the top panel of Figure 5, which covers the transition and coupling regions, is enlarged and presented in a different color scale covering a narrower range since the wave components propagating in the IWG have significantly larger magnitude. The middle panel in Figure 5 clearly demonstrates the coupling of surface modes between adjacent surfaces. The bottom panel in Figure 5 is a depiction of the variation of the magnitude of the Poynting vector () along the cutlines SU and SL in Figure 1, denoted by the horizontal arrows to the left of the middle panel of Figure 5. Only the surface modes of the lower PC are excited throughout the transition region, as is almost zero on the upper surface, as depicted in the plot of along the surfaces. In contrast, evanescent coupling rapidly takes place when the coupling region is encountered (Figure 5). The beat length denoted by the arrow between the two vertical dashed lines in Figure 5 is calculated to be which is in very good agreement with the value () predicted from Figure 4(a).

The middle and bottom panels of Figure 5 demonstrate that the coupling symmetric and antisymmetric modes are out of phase, as discussed in the previous section. It is more clearly visible on the plot of that the peaks for the upper and lower PC surfaces never coincide along the surfaces, bottom panel of Figure 5. The plot of at the output of the coupling PC surfaces, denoted by the vertical dash-dotted line in the middle panel of Figure 5, demonstrates that wave output is primarily from the lower surface at , as the coupling region length is such that as two complete beats take place over the two PC surfaces.

In general, the ratio of collected power between the surfaces reveals frequency dependence as depicted in Figure 6, which presents the variation of with vertical position at the output side and frequency. The maximum frequencies of the output power tightly depend on the WG width between the upper and lower PC surfaces with the lattice constant . Maximum power peaks are obtained at y/a = 0.75 for the coupling surface modes. These modes are travelling electromagnetic waves, which are bound to the interface between PCs and uniform host (air). Figure 6 shows that most of the EM energy is output from the upper surface at , while the situation converses at . There is roughly an oscillatory behavior with frequency away from the band edges. For sensor applications, the device can be operated at frequencies for which wave output is almost totally from one of the surfaces such that slight variations in the refractive index of the host result in disruption of this large contrast.

Figure 7 demonstrates SW propagation at and . Both SWs on the two surfaces in Figures 7(a) and 7(b) travel in phase with and without beating (evanescent coupling), respectively. This is clearly visible in the plots of along the PC surfaces calculated along the horizontal dash-dotted lines in Figure 7(a) where peaks of for the upper and lower surfaces occur at the same positions at all values, bottom panels in Figure 7. Such in-phase operation could be employed for sensing either to reduce the sensor size or to enhance the sensitivity since the effective interaction length between the SW and the gas phase is increased by about times.

(a)

(b)

It should be noted from Figures 5 and 7(a) that SWs exit the PC surfaces unevenly for either or . This can easily be deduced from the variation of over the output side (vertical dash-dotted line in Figures 5 and 7(a)). However, SWs at leave in an antisymmetric manner with a zero field strength in the midway (Figure 7(b)) as was pointed out in the discussion of Figure 4(b). There is no coupling at and above this frequency, where the propagating mode is solely the antisymmetric mode in Figure 3(b).

It should be noted that peak positions of the curves for neither upper nor the lower surface can be fitted by a smoothly varying envelope for (Figure 7(a)). In contrast, the corresponding curves in Figure 5 are more likely to be fitted by Gaussian-like curves for . Figure 7(a) suggests that there is interplay of two beating phenomena occurring at different spatial periods. One is obviously due to evanescent coupling. The other mechanism can ascribe to beating due to the finite length of the PC surfaces where two equivalent counter-propagating modes with complementary wave numbers with respect to interfere [49]. Due to this secondary beating, the beat length in Figure 7(a) is calculated as , a value significantly smaller than the one expected from Figure 4(a), i.e.,. The secondary beating length increases as approaches , where of the counter-propagating modes decreases.

4. Conclusion

In conclusion, evanescent coupling between surface modes of parallel two-dimensional photonic crystals is numerically demonstrated. While there exists a single degenerate surface band for two identical surfaces infinitely apart, the degeneracy is lifted as the surfaces are brought together at a length scale compared with lattice periodicity. This, in turn, gives rise to the occurrence of two localized surface bands within the band gap with opposite definite parities. The discrepancy of these modes in the reciprocal space suggests that effective coupling length increases linearly with frequency up to a point where the group velocities of the two interacting modes become equal, after which a rapid fall of the coupling length is observed. Finite-difference time-domain simulations demonstrate surface mode coupling at length scales in agreement with the values predicted from the band structure calculations. Coupling modes propagate along the two surfaces in phase if the group velocities are equal. This can be utilized to achieve highly sensitive compact refractive index sensors, as the effective interaction length increases significantly due to staggering by evanescent coupling.

Data Availability

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The author declares that there are no known conflicts of financial interest or personal relationships that could have appeared to influence the work reported in this study.

References

E. Yablonovitch, “Inhibited spontaneous emission in solid-state physics and electronics,” Physical Review Letters, vol. 58, no. 20, pp. 2059–2062, 1987.
View at: Google Scholar
R. D. Meade, A. Devenyi, J. D. Joannopoulos, O. L. Alerhand, D. A. Smith, and K. Kash, “Novel applications of photonic band gap materials: low‐loss bends and high Q cavities,” Journal of Applied Physics, vol. 75, no. 9, pp. 4753–4755, 1994.
View at: Google Scholar
J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the Flow of Light, Princeton University Press, Princeton, NJ, USA, 2011.
T. Baba, N. Fukaya, and J. Yonekura, “Observation of light propagation in photonic crystal optical waveguides with bends,” in Proceedings of the Integrated Photonics Research, Optical Society of America, 1999.
View at: Google Scholar
S. F. Mingaleev and Y. S. Kivshar, “Nonlinear transmission and light localization in photonic-crystal waveguides,” Journal of the Optical Society of America B, vol. 19, no. 9, pp. 2241–2249, 2002.
View at: Google Scholar
L. Jiu-Sheng, L. Han, and Z. Le, “Compact four-channel terahertz demultiplexer based on directional coupling photonic crystal,” Optics Communications, vol. 350, pp. 248–251, 2015.
View at: Google Scholar
M. Koshiba, “Wavelength division multiplexing and demultiplexing with photonic crystal waveguide couplers,” Journal of Lightwave Technology, vol. 19, no. 12, p. 1970, 2001.
View at: Publisher Site | Google Scholar
M. B. Yucel, A. Cicek, and B. Ulug, “Polarization-independent unidirectional light transmission by an annular photonic crystal prism,” Photonics and Nanostructures-Fundamentals and Applications, vol. 11, no. 3, pp. 270–278, 2013.
View at: Google Scholar
A. O. Cakmak, E. Colak, A. E. Serebryannikov, and E. Ozbay, “Unidirectional transmission in photonic-crystal gratings at beam-type illumination,” Optics Express, vol. 18, no. 21, pp. 22283–22298, 2010.
View at: Google Scholar
S. Boscolo, M. Midrio, and C. G. Someda, “Coupling and decoupling of electromagnetic waves in parallel 2D photonic crystal waveguides,” IEEE Journal of Quantum Electronics, vol. 38, no. 1, pp. 47–53, 2002.
View at: Google Scholar
X. M. Bendana and F. J. García de Abajo, “Power transfer between neighboring planar waveguides,” Optics Express, vol. 20, no. 3, pp. 3152–3157, 2012.
View at: Google Scholar
A. Martinez, F. Cuesta, and J. Marti, “Ultrashort 2-D photonic crystal directional couplers,” IEEE Photonics Technology Letters, vol. 15, no. 5, pp. 694–696, 2003.
View at: Google Scholar
M. Thorhauge, L. H. Frandsen, and P. I. Borel, “Efficient photonic crystal directional couplers,” Optics Letters, vol. 28, no. 17, pp. 1525–1527, 2003.
View at: Google Scholar
F. Cuesta-Soto, A. Martínez, J. Garcia et al., “All-optical switching structure based on a photonic crystal directional coupler,” Optics Express, vol. 12, no. 1, pp. 161–167, 2004.
View at: Google Scholar
J. Zimmermann, M. Kamp, A. Forchel, and R. März, “Photonic crystal waveguide directional couplers as wavelength selective optical filters,” Optics Communications, vol. 230, no. 4–6, pp. 387–392, 2004.
View at: Google Scholar
I. Park, H.-S. Lee, H.-J. Kim et al., “Photonic crystal power-splitter based on directional coupling,” Optics Express, vol. 12, no. 15, pp. 3599–3604, 2004.
View at: Google Scholar
T. B. Yu, M. H. Wang, X. Q. Jiang, Q. H. Liao, and J. Y. Yang, “Ultracompact and wideband power splitter based on triple photonic crystal waveguides directional coupler,” Journal of Optics A: Pure and Applied Optics, vol. 9, no. 1, pp. 37–42, 2006.
View at: Google Scholar
T. Liu, A. R. Zakharian, M. Fallahi, J. V. Moloney, and M. Mansuripur, “Design of a compact photonic-crystal-based polarizing beam splitter,” IEEE Photonics Technology Letters, vol. 17, no. 7, pp. 1435–1437, 2005.
View at: Google Scholar
R. D. Meade, K. D. Brommer, A. M. Rappe, and J. D. Joannopoulos, “Electromagnetic bloch waves at the surface of a photonic crystal,” Physical Review B, vol. 44, no. 19, pp. 10961–10964, 1991.
View at: Google Scholar
A. I. Rahachou and I. V. Zozoulenko, “Waveguiding properties of surface states in photonic crystals,” Journal of the Optical Society of America B, vol. 23, no. 8, pp. 1679–1683, 2006.
View at: Google Scholar
T.-W. Lu, S.-P. Lu, L.-H. Chiu, and P.-T. Lee, “Square lattice photonic crystal surface mode lasers,” Optics Express, vol. 18, no. 25, pp. 26461–26468, 2010.
View at: Google Scholar
S. Xiao and M. Qiu, “Optical microcavities based on surface modes in two-dimensional photonic crystals and silicon-on-insulator photonic crystals,” Journal of the Optical Society of America B, vol. 24, no. 5, pp. 1225–1229, 2007.
View at: Google Scholar
T.-W. Lu, Y. H. Hsiao, W. D. Ho, and P. T. Lee, “High-index sensitivity of surface mode in photonic crystal hetero-slab-edge microcavity,” Optics Letters, vol. 35, no. 9, pp. 1452–1454, 2010.
View at: Google Scholar
R. Moussa, T. Koschny, and C. M. Soukoulis, “Excitation of surface waves in a photonic crystal with negative refraction: the role of surface termination,” Physical Review B, vol. 74, no. 11, Article ID 115111, 2006.
View at: Google Scholar
B. Wang, W. Dai, A. Fang et al., “Surface waves in photonic crystal slabs,” Physical Review B, vol. 74, no. 19, Article ID 195104, 2006.
View at: Google Scholar
E. Moreno, F. J. Garcia-Vidal, and L. Martín-Moreno, “Enhanced transmission and beaming of light via photonic crystal surface modes,” Physical Review B, vol. 69, no. 12, Article ID 121402, 2004.
View at: Google Scholar
P. Kramper, M. Agio, C. M. Soukoulis et al., “Highly directional emission from photonic crystal waveguides of subwavelength width,” Physical Review Letters, vol. 92, no. 11, Article ID 113903, 2004.
View at: Google Scholar
I. Bulu, H. Caglayan, and E. Ozbay, “Beaming of light and enhanced transmission via surface modes of photonic crystals,” Optics Letters, vol. 30, no. 22, pp. 3078–3080, 2005.
View at: Google Scholar
T. Yu, S. Li, N. Liu, T. Wang, Q. Liao, and X. Xu, “Highly efficient coupling between inner and surface fields in photonic crystal waveguides,” IEEE Photonics Technology Letters, vol. 25, no. 15, pp. 1496–1499, 2013.
View at: Google Scholar
Q. Wang, L. Zhang, and Q. Li, “Beam splitting at the output of photonic crystal waveguides with discrete surface point defects,” Optics Express, vol. 18, no. 23, pp. 24245–24257, 2010.
View at: Google Scholar
J.-K. Yang, S.-H. Kim, G.-H. Kim, H.-G. Park, Y.-H. Lee, and S. B. Kim, “Slab-edge modes in two-dimensional photonic crystals,” Applied Physics Letters, vol. 84, no. 16, pp. 3016–3018, 2004.
View at: Google Scholar
A. I. Rahachou and I. V. Zozoulenko, “Light propagation in finite and infinite photonic crystals: the recursive green’s function technique,” Physical Review B, vol. 72, no. 15, Article ID 155117, 2005.
View at: Google Scholar
S. Kita, S. Hachuda, K. Nozaki, and T. Baba, “Nanoslot laser,” Applied Physics Letters, vol. 97, no. 16, Article ID 161108, 2010.
View at: Publisher Site | Google Scholar
E. Chow, A. Grot, L. W. Mirkarimi, M. Sigalas, and G. Girolami, “Ultracompact biochemical sensor built with two-dimensional photonic crystal microcavity,” Optics Letters, vol. 29, no. 10, pp. 1093–1095, 2004.
View at: Google Scholar
A. Di Falco, L. O’faolain, and T. F. Krauss, “Chemical sensing in slotted photonic crystal heterostructure cavities,” Applied Physics Letters, vol. 94, no. 6, Article ID 063503, 2009.
View at: Google Scholar
S. G. Johnson and J. D. Joannopoulos, “Block-iterative frequency-domain methods for Maxwell’s equations in a planewave basis,” Optics Express, vol. 8, no. 3, pp. 173–190, 2001.
View at: Google Scholar
H. Kurt, M. N. Erim, and N. Erim, “Various photonic crystal bio-sensor configurations based on optical surface modes,” Sensors and Actuators B: Chemical, vol. 165, no. 1, pp. 68–75, 2012.
View at: Google Scholar
M. Wellenzohn, E. Melnik, P. Muellner, L. O’Faolain, and R. Hainberger, “Design of a photonic crystal defect waveguide biosensor operating in aqueous solutions at 1.34 μm,” Proceedings, vol. 2, Article ID 1026, 2018.
View at: Publisher Site | Google Scholar
Y. Liu, S. Li, H. Chen, J. Li, and W. Zhang, “Surface plasmon resonance-induced high sensitivity refractive index sensor with adjustable measurement range based on an evanescent field-enhanced D-shaped five-hole photonic crystal fiber,” Journal of Physics D: Applied Physics, vol. 53, no. 11, Article ID 115107, 2020.
View at: Google Scholar
R. Jannesari, G. Pühringer, T. Grille, and B. Jakoby, “Design and analysis of a slot photonic crystal waveguide for highly sensitive evanescent field absorption sensing in fluids,” Micromachines, vol. 11, no. 8, p. 781, 2020.
View at: Google Scholar
K. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Transactions on Antennas and Propagation, vol. 14, no. 3, pp. 302–307, 1966.
View at: Google Scholar
S. Aachboun and P. Ranson, “Deep anisotropic etching of silicon,” Journal of Vacuum Science and Technology A: Vacuum, Surfaces, and Films, vol. 17, no. 4, pp. 2270–2273, 1999.
View at: Google Scholar
N. Jukam and M. S. Sherwin, “Two-dimensional terahertz photonic crystals fabricated by deep reactive ion etching in Si,” Applied Physics Letters, vol. 83, no. 1, pp. 21–23, 2003.
View at: Google Scholar
F. Marty, L. Rousseau, B. Saadany et al., “Advanced etching of silicon based on deep reactive ion etching for silicon high aspect ratio microstructures and three-dimensional micro-and nanostructures,” Microelectronics Journal, vol. 36, no. 7, pp. 673–677, 2005.
View at: Google Scholar
K. M. Ho, C. T. Chan, and C. M. Soukoulis, “Existence of a photonic gap in periodic dielectric structures,” Physical Review Letters, vol. 65, no. 25, pp. 3152–3155, 1990.
View at: Google Scholar
A. David, H. Benisty, and C. Weisbuch, “Fast factorization rule and plane-wave expansion method for two-dimensional photonic crystals with arbitrary hole-shape,” Physical Review B, vol. 73, no. 7, Article ID 075107, 2006.
View at: Google Scholar
E. Moreno, F. J. Garcia-Vidal, and L. Martín-Moreno, “Fast factorization rule and plane-wave expansion method for two-dimensional photonic crystals with arbitrary hole-shape,” Physical Review B, vol. 73, no. 7, Article ID 075107, 2006.
View at: Google Scholar
A. Mekis, S. Fan, and J. D. Joannopoulos, “Bound states in photonic crystal waveguides and waveguide bends,” Physical Review B, vol. 58, no. 8, pp. 4809–4817, 1998.
View at: Google Scholar
D. Zhao, Z. Liu, C. Qiu, Z. He, F. Cai, and M. Ke, “Surface acoustic waves in two-dimensional phononic crystals: dispersion relation and the eigenfield distribution of surface modes,” Physical Review B, vol. 76, no. 14, Article ID 144301, 2007.
View at: Google Scholar
R. Courant, K. Friedrichs, and H. Lewy, “On the partial difference equations of mathematical physics,” IBM Journal of Research and Development, vol. 11, no. 2, pp. 215–234, 1967.
View at: Google Scholar
J. P. Berenger, “A perfectly matched layer for the absorption of electromagnetic waves,” Journal of Computational Physics, vol. 114, no. 2, pp. 185–200, 1994.
View at: Google Scholar

Copyright

Copyright © 2022 Melike Behiye Yücel. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

325

Downloads

536

Citations