Resonant transmission of self-collimated beams through coupled zigzag-box resonators : slow self-collimated beams in a photonic crystal

The resonant transmission of self-collimated beams through zigzag-box resonators is demonstrated experimentally and numerically. Numerical simulations show that the flat-wavefront and the width of the beam are well maintained after passing through zigzag-box resonators because the up and the down zigzag-sides prevent the beam from spreading out and the wavefront is perfectly reconstructed by the output zigzag-side of the resonator. Measured split resonant frequencies of twoand three-coupled zigzag-box resonators are well agreed with those predicted by a tight binding model to consider optical coupling between the nearest resonators. Slowing down the speed of self-collimated beams is also demonstrated by using a twelve-coupled zigzag-box resonator in simulations. Our work could be useful in implementing devices to manipulate self-collimated beams in time domain. © 2012 Optical Society of America OCIS codes: (230.4555) Coupled resonators; (260.2030) Dispersion; (260.5950) Selffocusing; (230.5298) Photonic crystals. References and links 1. H. Kosaka, T. Kawashima, A. Tomita, M. Notomi, T. Tamamura, T. Sato, and S. Kawakami, “Self-collimating phenomena in photonic crystals,” Appl. Phys. Lett. 74, 1212–1214 (1999). 2. P. T. Rakich, M. S. Dahlem, S. Tandon, M. Ibanescu, M. Solja čiv́, G. S. Petrich, J. D. Joannopoulos, L. A. Kolodziejski, and E. P. Ippen, “Achieving centimetre-scale supercollimation in a large-area two-dimensional photonic crystal,” Nat. Mater. 5, 93–96 (2006). 3. Z. Lu, S. Shi, J. A. Murakowski, G. J. Schneider, C. A. Schuetz, and D. W. Prather, “Experimental Demonstration of Self-Collimation inside a Three-Dimensional Photonic Crystal,” Phys. Rev. Lett. 96, 173902 (2006). 4. S.-H. Kim, T.-T. Kim, S. S. Oh, J.-E. Kim, H. Y. Park, and C.-S Kee, “Experimental demonstration of selfcollimation of spoof surface plasmons,” Phys. Rev. B 83, 165109 (2011). 5. D. Chigrin, S. Enoch, C. Sotomayor Torres, and G. Tayeb, “Self-guiding in two-dimensional photonic crystals,” Opt. Express11, 1203–1211 (2003). 6. X. Yu and S. Fan, “Bends and splitters for self-collimated beams in photonic crystals,” Appl. Phys. Lett. 83, 3251–3253 (2003). 7. D. W. Prather, S. Shi, D. M. Pustai, C. Chen, S. Venkataraman, A. Sharkawy, G. J. Schneider, and J. Murakowski, “Dispersion-based optical routing in photonic crystals,” Opt. Lett. 29, 50–52 (2004). #161840 $15.00 USD Received 20 Jan 2012; revised 14 Mar 2012; accepted 15 Mar 2012; published 26 Mar 2012 (C) 2012 OSA 9 April 2012 / Vol. 20, No. 8 / OPTICS EXPRESS 8309 8. S.-G. Lee, S. S. Oh, J.-E. Kim, H. Y. Park, and C.-S. Kee, “Line-defect-induced bending and splitting of selfcollimatedbeams in two-dimensional photonic crystals,” Appl. Phys. Lett. 87, 181106 (2005). 9. T.-T. Kim, S.-G. Lee, S.-H. Kim, J.-E Kim, H. Y. Park, and C.-S. Kee, “Asymmetric Mach-Zehnder filter based on self-collimation phenomenon in two-dimensional photonic crystals,” Opt. Express 18, 5384–5389 (2010). 10. T.-T. Kim, S.-G. Lee, S.-H. Kim, J.-E Kim, H. Y. Park, and C.-S. Kee, “Ring-type Fabry-Perot filter based on the self-collimation effect in a 2D photonic crystal,” Opt. Express 18, 17106-17113 (2010). 11. Z. Li, H. Chen, Z. Song, F. Yang, and S. Feng, “Finite-width waveguide and waveguide intersections for selfcollimated beams in photonic crystals,” Appl. Phys. Lett. 85, 4834–4386 (2004). 12. A. Taflove,Computational Electrodynamics : The Finite-Difference Time-Domain Method (Artech House, Boston, 1995). 13. http://ab-initio.mit.edu/wiki/index.php/Meep. 14. S.-G. Lee, J.-S. Choi, J.-E. Kim, H. Y. Park, and C.-S. Kee, “Reflection minimization at two-dimensional photonic crystal interfaces,” Opt. Express 16, 4270–4277 (2008). 15. T.-T. Kim, S.-G. Lee, M.-W. Kim, H. Y. Park, and J.-E. Kim, “Experimental demonstration of reflection minimization at two-dimensional photonic crystal interfaces via antireflection structures,” Appl. Phys. Lett. 95, 011119 (2009). 16. C.-S. Kee and H. Lim, “Coupling characteristics of localized photons in two-dimensional photonic crystals,” Phys. Rev. B67, 073103 (2003). 17. A. Yariv, Y Xu, R. K. Lee, and A. Scherer, “Coupledresonator optical waveguide: a proposal and analysis,” Opt. Lett. 24, 711–713 (1999). 18. M. Bayindir, B. Temelkuran, and E. Ozbay, “Tight-binding description of the coupled defect modes in threedimensional photonic crystals,” Phys. Rev. Lett. 84, 2140–2143 (2000). 19. T. F. krauss, “Why do we need slow light?” Nat. Photon. 2, 448–450 (2008). 20. T. Baba, “Slow light in photonic crystals,” Nat. Photon. 2, 465–473 (2008). 21. J. B. Khurgin, “Slow light in various media: a tutorial,” Adv. Opt. Photon. 2, 287–318 (2010). 22. E. Ozbay, A. Abeyta, G. Tuttle, M. Tringides, R. Biswas, C. T. Chan, C. M. Soukoulis, and K. M. Ho, “Measurement of a three-dimensional photonic band gap in a crystal structure made of dielectric rods,” Phys. Rev. B 50, 1945–1948 (1994). 23. M. Notomi, E. Kuramochi, and T. Tanabe, “Large-scale arrays of ultrahigh-Q coupled nanocavities,” Nat. Photonics2, 741–747 (2008).


Introduction
A self-collimation phenomenon, diffractionless propagation of a light beam in photonic crystals (PCs), has been of great interest in recent years because it could provide a new way to manipulate light propagation in PCs [1][2][3][4][5].It has been experimentally demonstrated that selfcollimated beams can be well guided in PCs without the use of any physical boundary and effectively routed by employing the bends and splitters [6][7][8][9][10].Moreover, self-collimated beams can be crossed without cross talk [11], hence optical devices based on the self-collimated beams have intrinsic potential for high density photonic integrated circuits.
A resonator has been of fundamental and practical interest in optical devices such as cavities, waveguides, filters, couplers, and so on.A mirror is necessary to make a resonator.A 45-degree mirror to give rise to total internal reflection of a self-collimated beam has been proposed [6,8].However, it is difficult to make a resonator composed of 45-degree mirrors because the beams reflected from 45-degree mirrors can hardly make destructive interference.Thus, to realize a resonator of self-collimated beams, it is necessary to make a high quality mirror to reflect strongly self-collimated beams into a backward direction.
Designing a mirror of self-collimated beams can start from reviewing fundamental properties of a self-collimated beam.One of the properties is that a self-collimated beam has a flat-wavefront perpendicular to a propagation direction.One can reasonably conjecture that an appropriate geometry to break the flat-wavefront may cause strong back-reflection of selfcollimated beams.A zigzag-shape line-defect can be one of structures to break the flatwavefront of a self-collimated beam.We will show that a zigzag-shape line-defect can reflect strongly incident self-collimated beams into a backward direction and act as a practical mirror of self-collimated beams.
A self-collimated beam has an uncertainty of a wavevector along a direction perpendicular to the beam propagation direction due to the finite beam size.The uncertain wavevector components will make the beam spread out along the direction perpendicular to the propagation direction in the resonator.However, the resonant phenomenon requires that the beam size of an incident self-collimated beam should be preserved after passing through a resonator.In addition, for no spreading of the beam, an incident self-collimated beam broken by one side of a resonator should be reconstructed by another one.Thus, a resonator composed of zigzag-shape line-defects should be designed carefully and have a two-dimensional shape, like a zigzag-box.
In this paper, we experimentally and numerically demonstrated the resonant transmission of self-collimated beams through a designed zigzag-box resonator.Numerical simulations show that the flat-wavefront and the width of the beam are well maintained after passing through a zigzag-box resonator because the up and the down zigzag-sides prevent the beam from spreading out and the wavefront is perfectly reconstructed by the output zigzag-side of the resonator.We have also investigated resonant transmission characteristics of two-and threecoupled zigzag-box resonators.Split resonant frequencies of the coupled resonators are well agreed with those predicted by a tight-binding (TB) model to consider optical coupling between the nearest resonators.Slowing down the speed of self-collimated beams was also demonstrated by using a twelve-coupled zigzag-box resonator in simulations.

Results and discussion
It has been demonstrated that microwave experiments are very useful in testing properties of newly designed PC devices before applications in infra-red or visible wavelength ranges.In this study, we employ a square lattice PC composed of cylindrical alumina rods with dielectric constant of 9.7 in air.The lattice constant and radius of the rods are a = 5 mm and r = 0.4 a = 2 mm, respectively.Two parallel aluminum plates with periodically drilled holes are used to hold the alumina rods vertically.E-polarized microwaves (the electric-field parallel to the rod axes) of frequencies around 12.5 GHz can propagate with almost no diffraction along the ΓM direction inside the PC [9,10].An HP 8720C network analyzer and two horn antennas are uesed in experiments.Numerical simulations are performed by using a finite-difference timedomain (FDTD) method with a perfectly matched layer absorbing boundary condition [12].A freely available FDTD software package MEEP [13] was employed.The spatial resolutions are △x = △y = a/32.The discrete time step is set to △t = S△x/c, where the Courant factor S is chosen to be 0.5 for stable simulations.To obtain transmission spectra, a Gaussian pulse with a waist of w = 4 a is launched into the PC and the transmitted power P t (ω) is computed at the end of the PC.P t (ω) is normalized to an incident power P i (ω) calculated at near a source plane without the PC structure.Antireflection structures are employed to eliminate unwanted reflection at the PC-air interfaces [14,15].We first investigated transmission properties of self-collimated beams through a zigzag-shape mirror (ZSM) created by missing rods as shown in the inset in Fig. 1.One can clearly see that the flat-wavefront of a self-collimated beam cannot be maintained at the zigzag-shape PC-air interface and the beam is strongly reflected into a backward direction.The flat-wavefront of the reflected beam gets distorted as the height of the zigzag step increases and the reflectance decreases.When the length of the mirror increases, the flat-wavefront of the reflected beam is well maintained and the reflectance increases slightly.Figure 1 represents measured and simulated transmission spectra of self-collimated beams through a ZSM in a frequency range from 12.1 to 12.9 GHz.Measured (simulated) transmittance of self-collimated beams through the ZSM is less than 9 % (3 %), and thus the proposed ZSM can act clearly as a back-reflection mirror for self-collimated beams.
We next designed a two-dimensional zigzag-box resonator as shown in Fig. 2(a).Measured (simulated) transmission spectra plotted in Fig. 2(b) shows the resonant transmission of a selfcollimated beam of the frequency f 0 = 12.575 GHz (12.559GHz).The measured frequency is slightly lower than the simulated one.The discrepancy between them may come from the small uncertainty of dielectric constant of the alumina rod since the simulated frequency was matched to the measured one when the dielectric constant of the alumina rod is 9.68, which is slightly smaller than 9.7.The Q-factor estimated from the relation of f 0 /∆ f in the measurement (the FDTD simulation), where ∆ f is a full width half maximum, is about 607 (718).The simulated electric-field distribution of the beam at f 0 = 12.559 GHz represented in Fig. 2(c) shows that the flat-wavefront and the width of the resonant self-collimated beam are well maintained after passing through the zigzag-box resonator because the up and the down zigzag-sides prevent the beam from spreading out and the wavefront is perfectly reconstructed by the output zigzag-side of the resonator.
We also investigated coupled zigzag-box resonators to exhibit multi-resonant transmission peaks due to the evanescent coupling between individual resonator modes.Figure 3(a) shows the measured and the simulated transmission spectra of self-collimated beams through a twocoupled zigzag-box resonator with the inter-distance a = 15 √ 2 mm.The measured (simulated) resonant frequencies are Ω 1 = 12.518 GHz (12.502GHz) and Ω 2 = 12.630 GHz (12.615GHz).The simulated electric-field distributions of the resonant modes with resonant frequencies of Ω 1 and Ω 2 are represented in Fig. 3 (b) and 3(c), respectively.One can see that the resonant mode with Ω 1 (Ω 2 ) mimics an odd (even) symmetry mode, even though it has been known that a coupling between two identical resonant modes in coupled high dielectric cavities makes their frequency split into a lower frequency of even mode and a higher frequency of odd mode.Kee and Lim have demonstrated that the parity of split resonant modes in two-coupled resonators in a 2D PC could be switched due to the correlation of the inter-distance between resonators and the period of oscillatory decaying evanescent fields of resonant modes.[16].
It is well-known that a tight binding model is very useful in predicting resonant transmission characteristics of coupled resonant systems.According to the TB model under the consideration of the interaction between the nearest resonators only, two resonant frequencies of a two-coupled resonator are given by Ω 2  1,2 ≃ f 2 0 (1 ± β )/(1 ± α), where α and β are TB parameters and f 0 is a resonant frequency of single resonator.The physical meanings of the TB parameters are described in detail in Ref. [17].When the number of resonators increases, a transmission band is formed due to the evanescent coupling of individual resonant modes.Physical quantities such as the bandwidth ∆ f , the dispersion relation f (k) and the group velocity v g (k) of propagation modes are determined by coupling coefficient κ between the nearest resonators, ∆ f ≃ 2 f 0 |κ|, f (k) ≃ f 0 [1 + κ cos(ka)] and v g (k) ≃ −2π f 0 aκ sin(ka), where the coupling coefficient is defined as κ = β − α [17,18].
From the measured (simulated) resonant frequencies f 0 , Ω 1 and Ω 2 , we obtained the TB parameters and coupling coefficient α = −0.0085(0.0022), β = −0.0174(−0.0067), and κ = −0.0089,respectively.From the TB parameters, one can expect that the group velocities of selfcollimated beams should be less than 2π f 0 a|κ| ≃ 0.05 c in the whole transmission band (∆ f ≃ 0.22 GHz from 12.45 to 12.67 GHz) and approach to zero at band edges, provided that the TB model is valid.To check the validity of the TB model of the coupled zigzag-box resonator system for self-collimated beams, we compared theoretical resonant frequencies of a threecoupled zigzag-box resonator predicted by the TB model, ), with the measured (simulated) three resonant frequencies obtained from transmission spectra.Table 1 shows that the measured (simulated) three resonant frequencies well coincide with the frequencies predicted by the TB model.Hence, slowing down self-collimated beams could be possible by using coupled systems composed of numbers of zigzag-box resonators.In recent years, slowing the speed of light down to a remarkably low velocity has been of great interest due to its potential applications such as strong light-matter interactions, optical delay lines, optical buffers, and optical storage [19][20][21].
Table 1.Three resonant frequencies of a three-coupled resonator obtained from measured (simulated) transmission spectra.TB mea.(TB sim. ) frequencies were calculated by the measured (simulated) TB parameters obtained from the measured (simulated) two resonant frequencies of a two-coupled resonator.The unit of frequency is GHz.Measured  To verify the slow propagation of self-collimated beams, we investigated the transmission properties of self-collimated beams through a couple system composed of twelve zigzag-box resonators.The consecutive resonators were made in a PC of length L pc = 270 √ 2 mm and the total length of the resonator region is L res = 195 √ 2 mm.As shown in Fig. 4(a), the simulated transmission band extending from 12.45 GHz to 12.67 GHz agrees well with the band predicted by the TB model.To obtain the group velocity of light in a media, it is essential to find the dispersion relation because the group velocity of light is given by v where k is a wave vector in the media.
In this study, the dispersion relation of the slow light modes has been obtained from the fact that the total phase difference ∆φ between a light propagating through a PC with thickness L and air is given by (k pc −k air )L, where k pc and k air are wave vectors in a PC and air, respectively [22].After phase determinations for three different cases, (1) φ air ( f ) for air, (2) φ pc ( f ) for a PC without a resonator and (3) φ res ( f ) for a PC with the twelve-coupled zigzag-box resonator, Figure 4(b) shows the dispersion relations obtained from the phase calculations (black-thick line) and the TB model (red-thin line) in the transmission band.Overall tendency of the measured dispersion relation is in good agreement with the calculated one by the TB model.The step-like behavior of the measured dispersion relation comes from the finite number of cavities.In the simulations, we first recorded the time varying electric-fields at the end of the PC and then performed Fast Fourier Transforms to obtain the phase values.
The group velocities of self-collimated beams were obtained by taking the derivative of the calculated dispersion relation, v g ( f ) = 2π/[dk( f )/d f ] and the results were plotted in Fig. 4(c) with the theoretical curve from the TB model.The calculated group velocities (black-solid square) oscillate around the theoretical curve (red-thin line) and has local minimum values at resonant frequencies where the transmission exhibits peak values.The oscillatory behavior of the group velocity will disappear and the group velocity curve will be close to the theoretical curve, if the number of cavities becomes infinite [23].The group velocity noticeably decreases to values less than c/100 as a frequency approaches to edges of the transmission band.
It would be challengeable to make zigzag-box resonators in 2D PC slabs operating at optical frequencies, even though we demonstrated their performances in a microwave range.The coupled zigzag-box resonators fabricated in optical wavelength scales could be useful in designing devices to manipulate self-collimated beams in time domain such as delay lines and optical storages in optical communications.

Conclusion
In conclusion, the resonant transmission of self-collimated beams through the zigzag-box resonator was demonstrated experimentally and numerically.Using FDTD simulations, we showed that the flat-wavefront and the size of the resonant self-collimated beam are preserved after passing through the box resonator.We have investigated multi-resonant transmission characteristics of two-and three-coupled zigzag-box resonators and analyzed the split frequencies by using the TB model.Slowing down the speed of self-collimated beams demonstrated by using a twelve coupled zigzag-box resonator could be useful in implementing devices to manipulate self-collimated beams in time domain.

Fig. 1 .Fig. 2 .
Fig.1.Experimental and FDTD simulated transmittances of self-collimated beams through the zigzag-shape mirror in a range of frequency from 12.1 to 12.9 GHz.An inset depicts a zigzag-shape mirror to reflect strongly an incident self-collimated beam with a frequency of 12.5 GHz.

#Fig. 3 .
Fig. 3. Transmission spectra of self-collimated beams through a two-coupled zigzag-box resonator with two resonant frequencies of Ω 1 and Ω 2 (a).Black-thick and red-thin lines indicate experimental and FDTD simulated results, respectively.The simulated electricfield distributions of the resonant modes with resonant frequencies of Ω 1 (b) and Ω 2 (c).

Fig. 4 .
Fig. 4. (a) Transmission spectrum of self-collimated beams through a twelve-coupled zigzag-box resonator.(b) Dispersion relations of the transmission band obtained from the phase calculations (black-thick line) and the TB model (red-thin line).(c) Group velocities obtained from the FDTD simulations (black-solid rectangle) and the TB model (red-thin line).Dashed vertical lines represent the resonant frequencies of a twelve-coupled zigzagbox resonator.