Nanophotonic wire waveguides play an important role for the realization of highly dense integrated photonic circuits. The miniaturization of optoelectronic devices and realization of ultra-small integrated circuits strongly demand compact waveguide branches. T-shaped versions of nanophotonic wires are the first stage of both power splitting and optical-interconnection systems based on guided-wave optics; however, the acute transitions at the waveguide junctions typically induce huge bending losses in terms of radiated modes. Both 2D and 3D finite-difference time-domain methods are employed to monitor the efficient light propagation. By introducing appropriate combinations of dielectric posts around the dielectric-waveguide junctions within the region, we are able to reduce the bending losses dramatically and increase the transmission efficiency from low values of 18% in the absence of the dielectric posts to approximately 49% and 43% in 2D and 3D cases, respectively. These findings may lead to the implementation of such T-junctions in near-future high-density integrated photonics to deliver optical-clock signals via H-tree network.
© 2011 OSA
Nanophotonic wire waveguide structures play a crucial role for compact photonic integrated circuits (PICs) and are deployed for various applications in the design of optoelectronic devices such as power splitters, ring resonators, optical logic gates, wavelength filters, optical switches, and optical-interconnects in H-tree network [1–8]. High-density integrated circuits require the minimization of the dimensions of the photonic devices. In addition, sharp 90° bends are desirable to reduce the overall system size. This strongly militates for the use of dielectric waveguides with high refractive-index contrast. Although standard dielectric waveguides with low-index contrast provide low loss with respect to the transmission efficiency, their deployment requires relatively small bend angles with consequent large bend radii, which is not desirable for high-density PICs [9–11].
Increasing the index contrast may help us improve the optical confinement and increase the density of circuits. For single-mode operation, the dimensions of the waveguide should be reduced. Although Refs. 12 and 13 have reported the propagation losses for straight photonic wires to be ~0.24 dB/mm, losses become excessive as the bends become sharper, especially for 90° bends [5, 14]. At this stage, it may be helpful to emphasize that the bending losses will be of no consequence when bend radius is adjusted to be high enough. On the other hand, while the bend radius is dramatically decreased to have compact devices, several approaches have been described in the literature to reduce severe bending losses. Corner mirrors, waveguide resonators, air trenches, and topology optimization have been proposed [15–19]. In general, the complexity of the proposed solutions and their limited operating bandwidth are issues of concern. Here, we propose a new solution to enhance transmission efficiency of narrow-waveguide T-junctions that can be easily incorporated in H-tree networks. Most of the time, achieving high transmission only within a limited band is not enough and broadband solutions are preferable. We aimed therefore to design bends and branches which provide small bending losses over high bandwidth and occupy small area.
T-junctions are a basic part of a more complicated network configuration such as H-tree circuits that can be used for power distribution and optical-signal clock delivery. An effective design that promotes high power transmission through sharp corners can be systematically duplicated to build larger structures, such as optical interconnects that may contain many branches . Optical clock signals should be distributed to each receiver (photo-detector) with an equal time delay. To achieve that, high-index contrast waveguides with a special layout that resemble an H-tree shape can be used. The potential advantages of optical interconnects over electrical ones include high speed, high bandwidth, and low power consumption. High-speed processors handling large data loads demand the delivery of optical rather than electrical clock signals. Potential drawbacks of electrical interconnects include possible Joule heating that may degrade the attractive features of them such as faster speed and higher density. This argues for high-transmission junctions to minimize requisite optical power levels as well as stray light that may itself contribute to heating. High-efficiency methods that utilize resonant structures, however, may not be robust with respect to wavelength shifts due to any changes occurring in the environment. Such resonant wavelength shifts may occur with the changes in the chip temperature due to the temperature dependence of refractive index of the optical cavities. Hence, our approach is broadband from the start. The waveguide material should be transparent to over the entire range of potential operating wavelengths; consequently, we select Si as a waveguide material.
In this paper, a hybrid structure with a zero bending radius is discussed. The studies reported in Refs. 5, 9, 10, 11 and 14 necessitate small bend angles with consequence of large bend radii. The combination of photonic crystals (PC) and photonic wires is proposed in order to cope with the problem of high bend losses. A photonic dielectric wire is surrounded by air and the PC is composed of a square lattice of dielectric rods that are only present in the vicinity of the junctions where excess bend losses would otherwise appear. By means of the photonic band gap (PBG) mechanism [21–25], the PC behaves as a photonic insulator preventing the light from penetrating into the surrounding region resulting in higher light confinement. After we demonstrate numerically effective signal delivery of T-junction photonic wires both in 2D and 3D, we finally implemented 1 to 4 fanouts in an H-tree configuration in which 1 to 2n fanouts are possible, where n is an integer. We should emphasize that the proposed solution in Ref. 21 is only for a PC-waveguide structure. However, we here offer a hybrid structure that is composed of a PC and high-index contrast nanowire waveguide. Besides, the PC portion is retained only proximate to the junction itself. The remainder of the structure is Si-waveguide on top of an insulator layer.
The rest of the paper is organized as follows: To address the bend losses in narrow waveguides, we propose optimized waveguide corners designed for H-tree interconnects in Sec. I. Section II describes light transport performance for three cases: standard T-junctions, the same junction surrounded by PC, and the optimized composite version of the structure. Section III outlines 3D investigations of the optimized T-junction. Light transmission in an H-tree network is presented in Sec. IV. We give concluding remarks in Sec. V.
2. Two-dimensional T-junction design and enhanced power transmission
The original T-junction structure is made up of Si () surrounded by air () as shown in Fig. 1 . The width of the waveguide, is represented in terms of a constant, a (). The parameter a is defined shortly below. For numerical studies in 2D configurations, a transverse magnetic (TM) guided mode is assumed, in which case the corresponding nonzero electric- and magnetic-field components are Ez, Hx, and Hy respectively. In order to improve the poor power transmission capability of a standard T-junction, some modifications will be made in the area proximate to the junction.
For our 2D analysis of T-junction structures, we first carry out simulations for the propagation of electromagnetic waves by using the 2D finite-difference time-domain (FDTD) method . The spectral content of the transmitted power is calculated by taking the ratio of the output power to input power. For the correct normalization of the transmission efficiency, we need to detect both incident and output pulses. The input power is calculated in the absence of sharp corners. Output detector planes are placed outside the engineered area. To eliminate back reflections originating from the ends of the finite computational window, we used perfectly matched layer absorbing boundary conditions . In Fig. 2 the normalized transmission coefficient for the T-shaped structure is presented with a solid line. As can be seen from the figure, this structure has poor transmission efficiency of around 18%. In the lossless case, that value should correspond to 50% efficiency.
In Fig. 3 , a snap shot of the pulse propagation for the T-junction structure is shown to illustrate the severe power leakage in the junction area. As can be clearly observed, radiated field leaves the waveguide and consequently the resulting bend loss at the junction of the structure is very high. The 18% transmission efficiency in each branch gives rise to 64% radiation loss in total due to the bend loss.
A hybrid solution composed of a high-index-contrast waveguide and a PC is considered next in order to enhance the transmission efficiency. For that purpose, dielectric rods made of Si with radii r = 0.20a are placed periodically around the junction of the wire as shown in Fig. 4 . The periodicity of the dielectric rods’ locations is indicated by lattice constant a and it is taken to be nm in order to target telecom wavelengths. Such a straightforward placement of the periodic rods suppresses the light scattered at the bend, but only slightly improves the transmission efficiency as compared to standard T-junction. The result is shown with a dashed-green line in Fig. 2.
The small increment in the transmission efficiency can be improved upon as we show in the next step. We analyzed the structure carefully and decided to rearrange the locations and number of the PC rods to obtain better results. For this reason, two additional rods marked as “A” in Fig. 5 are located on the diagonal between the second-nearest neighbors as shown and the positions of the corner rods are displaced by (δx, δy). In a recent study , we performed both 2D and 3D studies on sharp 90° bends. For the optimum solution, we applied the same design procedure, mentioned in Ref. 28, for the rods of the sharp-bend structure to obtain efficient T-junction hybrid structures. We optimized the parameters of the crucial unit cells within the junction area. The radii of the designated Si rods (,) are chosen as (0.18a, 0.26a). The radii of the rest of the rods are kept the same, namely as 0.20a. The power transmission efficiency for the designed structure can be seen in Fig. 2 as the dotted-red line. With the help of that modified hybrid structure, the power transmission efficiency dramatically increases and reaches up to 48.6%. This value gets near to 97%, which implies the system to become highly efficient. We should emphasize that this value is valid under the assumption of absent insertion and scattering losses. As can be seen in Fig. 2, the normalized transmission stays well above 40% (ideal case is 50%) for a bandwidth of nearly nm. In Fig. 6 , one can observe a snapshot of pulse propagation in the structure when the normalized frequency of the continuous pulse is chosen to be. It should be noted that the located source size is kept sufficiently small that it can be considered as a point source. In that case, it isotropically emits and some part of the radiated power hits the PC front surfaces. That is why some undesired electric field appears around the source location. The negligible back reflection can be seen both qualitatively (visually) and quantitatively (transmission coefficient). The qualitative explanation is based on the inspection of the time domain snap shot of the electric field presented in Fig. 6. If there were a strong back reflection that may arise due to the presence of an acute corner, the forward and backward propagating waves would interfere. As a result, the spatial periodicity of the field may be disturbed. However, as we can see from Fig. 6, there is no strong back reflection having an effect on the periodicity of the wave oscillation.
The wavelength of the operating signal is taken to beµm. In that case, the optimal lattice constant is nm as mentioned previously. The efficient division of the optical pulse into the two branches can be observed in Fig. 6.
Even though surrounding the photonic wire waveguide directly by PCs provides slightly better results, the carefully adjusted displacement and incorporation of additional dielectric rods around the junction area dramatically improves the power transmission capability. The efficiency of 40% is the total transmission of the two branches. That means 20% transmission per branch. This can also be inferred from the transmission efficiency vs. frequency graph, shown in Fig. 2 (dashed-green line). However, in the optimized hybrid structure, the total transmission efficiency reaches near to 98% because the numerically reported efficiency at each branch is approximately 49%. The origin of the optimization stage was partially borrowed from the recent work . The starting point of the design is to select the parameters of the periodic dielectric medium that surround the corner region of a standard nano-wire waveguide of width w. The lattice constant, a of the periodic rods that should surround the corner region can be determined according to desired operating wavelength. The radii of dielectric rods, r can be selected in such a way that the band gap of the periodic structure should cover the waveguide mode. After these selections, the critical rods closest to the corner region should be manipulated in terms of their locations and radii to enhance the amount of light transportation. Besides, additional rods can be incorporated in the diagonal directions, viz. the cylinder rods represented by “A” in Fig. 5. The confinement mechanism is replaced by band gap effect instead of total internal reflection due to periodic arrangement of rods around the corner area. The adjustment of locations and radii of rods provide an effective 45° turn while implementing 90° bend. In the same study, we have proved that moving from 2D to 3D FDTD analysis, small deviations occur in the overall performance. The reason for that is the high-index contrast waveguide surrounded by lower refractive indices both above and below the core region. In the next section, 3D analysis of T-junction will be presented.
To check the robustness of the proposed design against some imperfections that may occur during the fabrication procedure we have made some modifications on the upper part of the designed structure as shown in Fig. 5. Firstly, the rod marked “B” as shown in Fig. 5 is shifted back and forth by 0.20a steps. An approximate 1% variation in the power transmission efficiency of the structure was observed. Secondly, the radius of rod B was changed within the interval 0.16a and 0.26a. Again, we obtained only 0.5% variations in the transmission efficiency. Lastly, we rearranged the radii of the rods that are placed next to rod “B”. In this case, the transmission efficiency is only affected by ~0.3%. These observations indicate that the top part of the T-junction is robust against changes in the parameters of the rods and inadvertent displacements from their ideal locations.
3. Three-dimensional analysis of T-junction
While the 2D analysis given above is indicative of performance, a 3D analysis is desirable to obtain more accurate results  and to understand the role of out-of-plane and substrate losses occur. A similar 3D data-extraction procedure as in 2D case is conducted in this section. The optimization steps carried out in 2D are repeated in 3D to finely tune and enhance the transmission efficiency of T-junctions. The schematic drawing of the 3D structure is shown in Fig. 7 . The T-junction photonic wire waveguide is surrounded by dielectric Si rods, with and laid upon a SiO2 substrate layer. Above the core region, air is used as cladding. The core region is sandwiched (along the z-direction) by SiO2 from the bottom and air from the top. In addition to that, the corner regions are accompanied by periodic dielectric rods around corner region. The radii of the silicon rods are selected as (117.6 nm). This change was made to obtain optimal results in 3D concerning the transmission efficiency. The height of the waveguide is taken to be h = 3a, which is nearly equal to 1.535 µm. The substrate beneath the Si rods and the wire is composed of silica with refractive index. The width of the waveguide is 0.20a. The parameters of the structure (waveguide width, radii of rods and slab height) are selected to have large bandwidth for the fundamental guided TM mode. It is known that the slab thickness is usually higher for TM mode than TE mode if one wants to obtain wide bandwidth . Hence, we selected the height of the slab for TM mode according to 3D FDTD analysis that ensures both high transmission for single 90° waveguide bend and single mode operation spanning whole band gap region [25, 28, 30]. The band-gap guided mode appears in the corner region due to surrounding periodic dielectric rods. The other regions are composed of uniform dielectric wire waveguide. As a result, the confinement mechanism is only due to the total-internal reflection. It can be feasible to fabricate such a waveguide with a high aspect ratio, which corresponds to a value of in our study, thanks to the recent advancements in the semiconductor fabrication technology .
The transmission efficiency in 3D configuration was calculated and the result is presented in the solid-blue line in Fig. 8 . As compared to 2D analysis, the transmission efficiency drops 4-5% near the central region of the spectrum, which is actually due to the out-of-plane loss. The finding dictates that the leading loss mechanism in 3D is, again, due to bend loss. For high-index contrast dielectric waveguides, the out-of-plane and substrate to leakage losses are not a major concern as we see from the 3D FDTD results. This is mainly because of the fact that the index contrast between core and cladding/substrate is high enough for the total internal reflection to give strong vertical field confinement. 4-5% efficiency drop corresponds to a power level change that occurs when 3D analysis is performed instead of 2D one. The out-of-plane and leakage to substrate losses are only considered in 3D. As a result, the comparison between 4 and 5% additional loss due to out-of-plane and leakage to substrate losses versus 64% bending loss shows the fact that the main loss contribution in these structures arises because of the sharp corners. It is again important to mention that we excluded insertionand propagation losses in the current work. On the other hand, we should emphasize the importance of the scattering losses that may occur due to sidewall roughness. It will be very crucial to determine the overall performance of fabricated systems; however, we do not address such loss being outside the scope of the present work. The other critical aspect of deployed systems is the input and output coupling. The size mismatch between various waveguides necessitates the implementation of spot-size converters. Novel solutions achieving efficient couplers were proposed in the literature [32, 33].
The efficient optical pulse propagation in T-junction both in 2D and 3D cases motivates us move to the next stage in which we designed a 1x4 H-tree network to deliver optical signal.
4. 1x4 optical pulse distribution in H-tree network
This section outlines a prototype for optical signal distribution in an H-tree optical waveguide network. It’s well-known that optical clock signal distribution for routing timing signals to all parts of a chip with the same propagation delays plays a key role to reduce skew of clocking in VLSI design [20, 34]. Therefore, the use of H-tree interconnections is inevitable for high-density photonic integrated circuits. Below, we demonstrate numerically transmission of light to 4 fanouts. By duplicating the optimized T-junction in a systematic manner, we obtain an H-tree configuration that is represented in Fig. 9 . This topology ensures the delivery of the input signal to output ports at equal propagation delays. The parameters of the structure are kept the same as before. The source is placed in the input waveguide. The output signals are detected at the end of each branch.
The calculated transmission efficiency of the 1x4 branches is around 23%, as can be seen in Fig. 10 . We observed narrowing in the bandwidth despite the high transmission. Such an observation can be explained as follows: Even though we do not directly incorporate a cavity configuration at the junction, surrounding the waveguide by the dielectric rods gives rise to a very low-Q cavity (Q<10) type of configuration. Each additional T-junction included in the structure contributes to the enhancement of Q associated with the junctions. As a result, the slightly larger Q diminishes the overall bandwidth of the nanowire waveguide.
In order to demonstrate the high transmission capability of the designed branches in the H-tree network configuration, we prepared a snap shot of the pulse propagation down the waveguide and also provided a movie (Media 1). The effective power distribution is apparent in Fig. 11 . Each time the input pulse is divided equally, the color of the electric field in the respective waveguide branch fades. Even though we have achieved 4 fanouts in the H-tree in Fig. 11, the idea is easily generalize to fanouts, with n an integer.
We summarize the best values of the transmission efficiencies achieved in the present work. It provides us a concise way to judge the performance of the proposed idea. The difference between the values of transmission efficiencies in 2D and 3D configurations corresponds to out-of-plane and leakage to substrate losses. (43x2)% efficiency represents 3D result. The 14% loss includes bending, out-of-plane and leakage to substrate losses for one-junction. On the other hand, (23x4)% efficiency represents the findings of 2D analysis. The 8% loss is due to the 1x4 H-tree network configuration. As a result, loss per junction is 4%. The difference between the 2D and 3D results is 10% that encloses out-of-plane and leakage to substrate losses for one-junction, under the assumption that the bending losses in 2D and 3D structures are equal. When there is no manipulation made to the corner region, the loss per junction is 64%. That means we report a sixteen times improvement in terms of bending loss reduction (64% loss reduces down to 4%).
The flow of the input power at each junction can be monitored as shown in Fig. 11. If we define the power split at first junction as T1 then at the second junction, the wave splits again. This efficiency is defined as T2. There exists a relationship between T1 and T2 as. When the value of T1 is taken as, as reported in Sec. 2 then T2 turns out to be. This number is very close to 23% efficiency value as reported in Sec. 4. We have 0.12 dB bending loss for 1x2 T-junction and 0.377 dB bending loss for 1x4 H-tree network case according to 2D study. On the other hand, 3D analysis gives us a loss value of 0.49 dB that includes other loss mechanisms in addition to bending loss. Ref. 8 experimentally reported 0.15 dB and 2.41 dB excess losses for 1x2 and 1x16 splitter cases, respectively that excluded insertion and propagation losses.
The low transmission through T-junctions made of high-index contrast narrow waveguides associated with bend losses is addressed and a novel solution is proposed. Transmission efficiencies extending up to 43% (ideal case corresponds to 50%) are predicted based on a 3D FDTD analysis. Sharp corners appearing along the waveguide branches were optimized for low-loss and uniform light propagation with a balanced propagation delay. Taking advantage of the optimized T-junction, we have shown efficient distribution of an input pulse in a 1x4 H-tree network configuration. The proposed solution in the present study may contribute to the widespread usage of photonic integrated circuits based on high-index Si nanophotonic wires.
I.H.G and H.K. gratefully acknowledge the financial support of the Scientific and Technological Research Council of Turkey (TUBITAK), Project numbers 108T717 and 110T306. H.K. also acknowledges support from the Turkish Academy of Sciences Distinguished Young Scientist Award (TUBA-GEBIP). D.S.C. acknowledges the support of the National Science Foundation under grant ECCS 0925713.
References and links
1. J. Gamet and G. Pandraud, “Ultralow-loss 1 x 8 splitter based on field matching Y junction,” IEEE Photon. Technol. Lett. 16, 2060–2062 (2004). [CrossRef]
2. B. E. Little, J. S. Foresi, G. Steinmeyer, E. R. Thoen, S. T. Chu, H. A. Haus, E. P. Ippen, L. C. Kimerling, and W. Greene, “Ultra-compact Si-SiO2 microring resonator optical channel dropping filters,” IEEE Photon. Technol. Lett. 10(4), 549–551 (1998). [CrossRef]
3. D. R. Lim, B. E. Little, K. K. Lee, M. Morse, H. H. Fujimoto, H. A. Haus, and L. C. Kimerling, “Micro-sized channel dropping filters using silicon waveguide devices,” Proc. SPIE 3847, 65–71 (1999). [CrossRef]
4. T. Fujisawa and M. Koshiba, “All-optical logic gates based on nonlinear slot-waveguide couplers,” J. Opt. Soc. Am. B 23(4), 684–691 (2006). [CrossRef]
5. W. Bogaerts, P. Dumon, D. V. Thourhout, D. Taillaert, P. Jaenen, J. Wouters, S. Beckx, V. Wiaux, and R. G. Baets, “Compact Wavelength-Selective Functions in Silicon-on-Insulator Photonic Wires,” IEEE J. Sel. Top. Quantum Electron. 12(6), 1394–1401 (2006). [CrossRef]
6. C. A. Barrios, “High-performance all-optical silicon microswitch,” Electron. Lett. 40(14), 862–863 (2004). [CrossRef]
7. G. T. Reed and A. P. Knights, Silicon Photonics: An Introduction (Wiley, 2004).
8. S. H. Tao, Q. Fang, J. F. Song, M. B. Yu, G. Q. Lo, and D. L. Kwong, “Cascade wide-angle Y-junction 1 x 16 optical power splitter based on silicon wire waveguides on silicon-on-insulator,” Opt. Express 16(26), 21456–21461 (2008). [CrossRef]
11. B. Chen, T. Tang, Z. Wang, H. Chen, and Z. Liu, “Flexible optical waveguides based on the omnidirectional reflection of one-dimensional photonic crystals,” Appl. Phys. Lett. 93(18), 181107 (2008). [CrossRef]
12. W. Bogaerts, R. Baets, P. Dumon, V. Wiaux, S. Beckx, D. Taillaert, B. Luyssaert, J. Van Campenhout, P. Bienstman, and D. Van Thourhout, “Nanophotonic waveguides in silicon-on-insulator fabricated with CMOS technology,” J. Lightwave Technol. 23(1), 401–412 (2005). [CrossRef]
13. P. Dumon, W. Bogaerts, V. Wiaux, J. Wouters, S. Beckx, J. Van Campenhout, D. Taillaert, B. Luyssaert, P. Bienstman, D. V. Thourhout, and R. Baets, “Low-loss SOI photonic wires and ring resonators fabricated with deep UV lithography,” IEEE Photon. Technol. Lett. 16(5), 1328–1330 (2004). [CrossRef]
15. M. Popovic, K. Wada, S. Akiyama, H. A. Haus, and J. Michel, “Air trenches for sharp silica waveguide bends,” J. Lightwave Technol. 20(9), 1762–1772 (2002). [CrossRef]
16. C. Manolatou, S. G. Johnson, S. Fan, P. R. Villeneuve, H. A. Haus, and J. D. Joannopoulos, “High-density integrated optics,” J. Lightwave Technol. 17(9), 1682–1692 (1999). [CrossRef]
17. B. Jalali and S. Fathpour, “Silicon Photonics,” J. Lightwave Technol. 24(12), 4600–4615 (2006). [CrossRef]
19. J. S. Jensen and O. Sigmund, “Topology optimization of photonic crystal structures: a high-bandwidth low-loss T-junction waveguide,” J. Opt. Soc. Am. B 22(6), 1191–1198 (2005). [CrossRef]
20. L. Pavesi and G. Guillot, Optical Interconnects: The Silicon Approach (Springer, 2006).
21. A. Mekis, J. C. Chen, I. Kurland, S. Fan, P. R. Villeneuve, and J. D. Joannopoulos, “High Transmission through Sharp Bends in Photonic Crystal Waveguides,” Phys. Rev. Lett. 77(18), 3787–3790 (1996). [CrossRef]
22. M. Notomi, “Manipulating light with strongly modulated photonic crystals,” Rep. Prog. Phys. 73(9), 096501–096557 (2010). [CrossRef]
23. H. Kurt and D. S. Citrin, “Photonic-crystal heterostructure waveguides,” IEEE J. Quantum Electron. 43(1), 78–84 (2007). [CrossRef]
26. A. Taflove and S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method (Artech House, 2005).
27. J. P. Berenger, “A perfectly matched layer for the absorption of electromagnetic waves,” J. Comput. Phys. 114(2), 185–200 (1994). [CrossRef]
28. H. Kurt, I. H. Giden, and K. Ustun, “Highly efficient and broadband light transmission in 90° nano-photonic wire waveguide bends,” J. Opt. Soc. Am. B 28(3), 495–501 (2011). [CrossRef]
29. Z. Sheng, D. Dai, and S. He, “Comparative Study of Losses in Ultrasharp Silicon-on-Insulator Nanowire Bends,” IEEE J. Sel. Top. Quantum Electron. 15(5), 1406–1412 (2009). [CrossRef]
30. S. G. Johnson, S. Fan, P. R. Villeneuve, J. D. Joannopoulos, and L. Kolodziejski, “Guided modes in photonic crystal slabs,” Phys. Rev. B 60(8), 5751–5758 (1999). [CrossRef]
31. G. Si, A. J. Danner, S. Lang Teo, E. J. Teo, J. Teng, and A. A. Bettiol, “Photonic crystal structures with ultrahigh aspect ratio in lithium niobate fabricated by focused ion beam milling,” J. Vac. Sci. Technol. B 29(2), 021205 (2011). [CrossRef]
32. H. Kurt and D. S. Citrin, “A novel optical coupler design with graded-index photonic crystals,” IEEE Photon. Technol. Lett. 19(19), 1532–1534 (2007). [CrossRef]
33. O. Cakmak, E. Colak, H. Caglayan, H. Kurt, and E. Ozbay, “High efficiency of graded index photonic crystal as an input coupler,” J. Appl. Phys. 105(10), 103708 (2009). [CrossRef]
34. T. Fukazawa, A. Sakai, and T. Baba, “H-tree-type optical clock signal distribution circuit using a Si photonic wire waveguide,” Jpn. J. Appl. Phys. 41(Part 2, No. 12B), L1461–L1463 (2002). [CrossRef]