Bandwidth Improvement for a Photonic Crystal Optical Y-splitter
Bandwidth Improvement for a Photonic Crystal Optical Y-splitter
Journal of the Optical Society of Korea. 2011. Sep, 15(3): 283-288
Copyright ©2011, Optical Society of Korea
  • Received : January 26, 2011
  • Accepted : August 10, 2011
  • Published : September 25, 2011
Export by style
Cited by
About the Authors
Mohammad, Danaie
Hassan, Kaatuzian

In this study, a wide-band photonic crystal Y-splitter for TE modes is proposed. A triangular lattice of air holes etched in a GaAs slab is used as the platform. In order to numerically analyze the structures,plane wave expansion (PWE) and finite difference time domain (FDTD) methods are used. In comparison with the structures reported in the literature, the proposed topology has a less complexity while it provides more than 100nm bandwidth. The simplicity of the design, its high transmission ratio and its wide bandwidth makes it a suitable choice for the implementation of photonic crystal integrated circuits.
Photonic crystals (PCs) are periodic dielectric structures which under certain conditions prohibit the propagation of electromagnetic waves within certain frequency bands. This phenomenon provides the ability of guiding and manipulating the flow of light [1 , 2] . PCs can be periodic in one, two or three dimensions. Slab PCs are periodic in two dimensions and have a finite thickness in the other dimension.
PCs can be used to design passive, active and nonlinear optical devices. Waveguides, filters, resonators, lasers, logical gates, etc. [3 - 8] can be implemented using PCs. Waveguides are realized by introducing a line defect in the PC structure [9] . The most popular method for creating a waveguide for TE modes is to remove a row of holes in a triangular lattice of air holes (which are etched in a dielectric slab). The resultant waveguide is often referred to as a “W1” waveguide.It is worth mentioning that in the same manner, sharp bends and branches can be created in PCs. The quality of optical bends and junctions are very important for many integrated optical circuits. The unwanted reflection from bends and branches can be a nuisance in many instances. Low reflection (high transmittance) devices are required for large scale integrated optical circuits. Various methods have been proposed in the literature to improve the quality of the bends [10 - 14] , Y-junctions and splitters [15 - 20] . Among the methods proposed to achieve higher transmission ratios, the most effective method is to modify the structure of the waveguide in the vicinity of the discontinuity. It can be done by introducing new holes [13 - 14] or by changing the radii of the holes in that region. Due to the complexity of the equations governing PC behavior, analytic optimization of PC structures is very difficult and in some instances almost impossible. Numerical optimization methods are usually used for this purpose. Using optimization techniques a wideband photonic crystal Y-splitter is proposed in this paper. The rest of the paper is organized as follows: In Section 2 we describe the numerical aspects of the work. Section 3 is devoted to the optimization of a Y-junction and consequently a Y-splitter for increasing the bandwidth. Section 4 presents the simulation results regarding the proposed structures,and the final section is devoted to the conclusions.
Different slab materials can be used for creating PCs, among them GaAs and silicon are the most frequently used materials.The refractive index of GaAs, which is used in our simulations, is about 3.4 for wavelengths around 1550 nm.In our designs a photonic crystal slab with a triangular pattern of circular holes is assumed. It can provide a considerable bandgap and can easily be fabricated using integrated circuit manufacturing technology. Assuming that “ a ” is the lattice constant of the structure, we choose the radius of the holes to be 0.3 a . Choosing larger holes will reduce the guiding strength in the vertical direction [14] . The slab thickness is considered to be 0.6 a and the PC is assumed to have an air cladding (which can be easily implemented using an air bridge.).
It is shown that in order to reduce the computational process, the slab structure can be modeled by a 2D structure of infinite height with an effective index [10 - 11] , [14] , [20 - 25] . In [21] a detailed study regarding the replacement of 3D analysis with 2D analysis has been made. It is shown that the 3D FDTD simulation results produce almost the same results as its 2D effective index approximation. The simplest method for calculation of the effective index is to compare the band diagrams of the 3D PC slab and the equivalent infinite 2D PC. The effective index of the 2D PC is chosen such that its band diagram approximately matches the band diagram of the 3D slab below the light cone. In our simulations we chose an effective index equal to 2.76. The same value has been used in [10 - 11] , [20] .Since the Y-splitter is not a straight waveguide, there always exists a vertical loss at the region where the propagation vector is changed (and the leaky modes can be excited).The mentioned loss cannot be modeled using the 2D effective index method. Although the vertical loss usually [10 - 11] does not create a considerable difference in the results but to be precise the transmittance notion in this paper incorporates the mentioned leakage.
The FDTD area is surrounded by a perfectly matched layer (PML) of thickness “ a ” and the mesh size used for analysis is chosen to be a /20. Also, the simulation time is chosen large enough so that increasing it further does not have any noticeable impact on the numerical results.
Using the plane wave expansion (PWE) method, it can be easily seen that the PC slab used in our simulations has a bandgap from 0.260 ( a/λ ) to 0.324 ( a/λ ). When a row of holes is removed from the structure of the PC, two waveguide modes are created in the bandgap region. Using the PWE method we find the band diagram to be as depicted in Fig. 1 . Index guiding and gap guiding are the different mechanisms for wave propagation in a photonic crystal slab.Since the photonic crystal slab is only periodic in two dimensions and has a finite size in the vertical direction, it does not have a complete bandgap. Hence some portion of the waves are radiated out of the slab and some modes use total reflection phenomenon to propagate inside the waveguide.The shaded side of Fig. 1 , which is labeled as the light cone, specifies the wave numbers for the modes that can be coupled to the air (the leaky modes). Therefore the wave-numbers that reside in the mentioned section cannot be used.
As it can be seen, the waveguide associated with the diagram shown in Fig. 1 has a single mode bandwidth between 0.265 ( a/λ ) to 0.287 ( a/λ ). A good PC branch should be able to use all that available W1 waveguide bandwidth. In the case where the central wavelength is chosen at 1550 nm,the lattice constant will be approximately equal to 430 nm.
The simplest geometry which can be used as a branch is shown in Fig. 2 a. The corresponding PC is the commonly used 2D triangular lattice of circular holes embedded in a 0.6a thick GaAs slab. By removing two columns of air holes adjacent to a W1 waveguide, a 120-degree Y-branch can be created. An ideal branch should divide the incident light beam to two distinct sections; with 50% of the incident power
Lager Image
The band diagram of the W1 waveguide used in thispaper.
Lager Image
(a) A simple PC Y-branch structure (b) Our proposedstructure for a Y-junction (c) The transmittance curve for thejunctions.
traveling in each direction (100% transmission ratio). Zero reflection is expected for the entire frequency range of the corresponding W1 waveguide; therefore the bandwidth used in this paper is normalized to the single-mode region of the W1 waveguide. For the branch shown in Fig. 2 a, FDTD simulation results show that the maximum transmission coefficient does not exceed 48% (24% for each arm), which is equal to a 52% reflection of the incident power. The high reflection ratio makes the mentioned structure unsuitable for many integrated optic applications. By introducing modifications to the structure of the branch the transmission spectrum can be improved [16 - 17] , [19] . There have been some detailed topologies reported in the literature which enhance the bandwidth [26 , 27] . The structure implemented in [26] provides a bandwidth of more than 100 nm. However; the structures presented in [26 , 27] are hard to define, manufacture and reproduce. Some holes are not circular; some are merged together; while for our design, the each hole is circular and has a defined radius.
Our proposed Y-junction structure is shown in Fig. 2 b.In comparison to the previous branch, three additional holes are introduced to the structure and the radius of one hole is also increased. An optimization algorithm (using FDTD method for evaluating the structures) is used to calculate the optimum values for the holes. The method used for the optimization of bandwidth is taken from [3] . It uses a neural network optimization engine to find the optimum values for the variables. An evaluation section calculates the transmittance curve for each set of input variables. The next block calculates the percentage of the transmittance curve which is more than 0.9. The obtained number is considered as the bandwidth. In other words, instead of using the conventional 3 dB criteria for calculating the bandwidth,we use the 90% limit which is approximately equal to 0.45 dB. The algorithm mentioned finds the input variable set that maximizes the bandwidth. In order to optimize the junction, the central junction’s hole’s radius is chosen as Var2, The right side hole’s radius is considered to be Var1 and the remaining two left-side holes are assigned an identical radius denoted by Var3. After the optimization phase,the optimum variables are obtained as: Var1=0.4a, Var2=0.31a and Var3=0.12a. Fig. 2 c shows the transmittance curve(summation of the two output arms) for the proposed junction in comparison to a simple junction. It is seen that a bandwidth equal to 94% is obtained.
Yang et. al. used a triangular lattice of dielectric rods to design a Y-Branch [17] . They showed that the Y-branch can be treated as a cavity that couples with the input and output waveguides. In order to enhance the transmittance,additional rods are added to the junction area and the corner rods are displaced. The movement of the corner rods can increase the volume of the cavity and make the cavity mode resonant with the waveguide modes. It is shown in [17] that using the coupled mode theory the reflection coefficient can be expressed as (1).
Lager Image
Where, ω 0 is the resonance frequency, τ i is the time constant regarding the amplitude decay of the resonance into the ith port. It is seen that if the (2) is satisfied, the reflection will be zero for ω = ω 0 .
Lager Image
We already know that due to symmetry τ 2 3 ; therefore the condition mentioned can be simplified as: τ 2 =2τ 1 . It means that the decay rate of resonance to the output ports should be half the decay rate of resonance to the input port. In a simple Y-junction all of the three time constants are equal. The tapered structure of the holes inserted in our proposed design, improves the coupling strength between the resonator and the input waveguide. It results in the reduction of τ 1 . If the tapering is designed so that τ 1 is reduced to half τ 2 , then the zero reflection condition can be satisfied. Numerical optimization is used to design such a tapering. In [17] two rods are placed at the output waveguides to reduce the coupling strength of the output waveguides (to increase τ 2 ), while we decrease τ 1 to obtain similar results.
In order to create a PC Y-splitter the aforementioned optimized Y-junction and two additional optimized PC 60-degree bends are combined. Figure 3 depicts the Y-splitter structure proposed in this paper and the recent bi-periodic Y-splitter proposed in [20] . Five additional holes have been added to the structure of a conventional PC bend to maximize the transmission bandwidth [14] . The radii of four of these holes are 0.11a and the other’s radius is 0.19 a .
For the bend embedded in Fig. 3 b, ΓJ denotes the 0-degree and equivalently 60-degree directions. When light passes through the bend, it can be assumed that, first it transmits along the 0-degree direction then turns in to the transient modes of 30-degree and finally turns back to the 60-degree direction. Assuming that ΓX represents the direction at the bend, namely 30-degrees, it can be easily shown that the mentioned structure has one guiding mode in the ΓJ direction and two guiding modes in the ΓX direction [14] . Therefore the coupling of energy between the modes is not perfect.It results in a higher reflection or equally a lower transmission from the bend. Modifying the topology of the waveguide in the vicinity of the bend can change the profile of the guiding modes in ΓX direction. In [14] , some additional
Lager Image
(a) The bi-periodic structure proposed in [20] (b) Thestructure proposed in this paper (c) The transmittance for thestructures in Fig. 3a and Fig. 3b.
holes with identical sizes are introduced in the structure to move one of the ΓX modes out of the band and enhance the coupling. It is a well-known principle that gradual change or tapering in the vicinity of discontinuity can lead to the reduction of reflection. It is expected that if the radii of the holes change in a gradual manner a lower reflection coefficient can be obtained. The optimization process of the previous section confirms the same concept.
The transmission spectrum of the structures shown in Fig. 3 a and Fig. 3 b are depicted in Fig. 3 c. As it can be seen, regarding the structure in Fig. 3 a, for only a small portion of the available W1 waveguide bandwidth, the transmission ratio exceeds 0.8. The transmittance of the splitter proposed in this paper is also shown in Fig. 3 c, whose normalized bandwidth is 79%, which is a considerable improvement in comparison with the structure reported in [20] . Furthermore;the proposed structure has the benefit of simplicity in regard to the bi-periodic structure in Fig. 3 a. The small holes in the mentioned figure are very close to each other. A slight manufacturing error might cause them to join together;which results in the total failure of the structure.
To have a better visual perception of the proposed Y-Splitter,the time domain simulation results regarding the amplitude of the vertical component of the magnetic field for an input wavelength equal to 1550 nm has been presented in Fig. 4 . As can be seen in this figure, the incident power is neatly divided in to two distinct sections.
We have determined the normalized holes’ radii as:0.4 a , 0.31 a and 0.12 a . However due to the implementation errors the radii of the mentioned holes usually have variations
Lager Image
Contour map of the magnitude of vertical magneticfield component (Hy) for the input wavelength of 1550 nm.
from the optimum values. Some holes might become bigger while the others become smaller than the desired value.These variations normally have a Gaussian distribution. Here each hole’s radius is independently assumed to have a variance of 4 nm. A set of 150 statistically different samples are produced, in which each variable has a Gaussian variation from the desired values. The histogram of the mentioned set is shown in Fig. 5 a.
As it is seen, the first variable varies independently in the 0.38 a -0.42 a range, the second in the 0.29 a -0.33 a range and finally the last variable in the 0.1 a -0.14 a range. In reality the variation is much less than what we have considered here. The normalized bandwidth is calculated for each of the 150 cases and the histogram of the obtained bandwidth is shown in Fig. 5 b. The maximum bandwidth obtained is 0.79 according to our simulations; however if we consider the samples with a bandwidth more than 0.70 as acceptable, then 68% of all our samples have a bandwidth higher than 0.7. The mentioned histogram shows the robustness of our design to harsh implementation variations, where each hole can have a 4 nm variation in its radius.
What we considered before was the worst case scenario;however in reality the hole sizes usually vary in the same manner. They either become a bit wider or narrower. In order to simulate the mentioned effect, we first assume that all holes become 0.01a larger (0.41 a , 0.32 a and 0.13 a and then 0.01 a smaller (0.39 a , 0.30 a and 0.11 a ), and we present the results in Fig. 6 to analyze the effect. As depicted in Fig. 6 , when the radii of the holes are decreased the right side of the transmittance curve moves up and the left side moves down. The opposite happens when the hole sizes are increased. Based on that observation, in order to have the maximum bandwidth, the values regarding the normal state in Fig. 6 should be used.
In order to investigate how the three added holes improve
Lager Image
(a) Histogram for the variables applied in optimizationof Y-junction used for Monte-Carlo simulation (b)Histogram of the Y-splitter bandwidth obtained afterMonte-Carlo simulation.
Lager Image
Transmittance for ±0.01a variation in extra holes’radii.
the transmission spectrum, the transmittance curves for different configurations are presented in Fig. 7 . The notation r 1 -r 4
Lager Image
Transmittance for different junction configurations:(I) r1= 0.4a r2-r4 =0 (II) r1 = 0.4a r2= 0.31a r3-r4=0 (III) r1 =0.4a r2= 0.31a r3 = 0.12a r4 = 0 (IV) r1 = 0.4a r2= 0.31a r3= 0.12a r4 = 0.12a (V) r1= 0.3a r2-r4 =0.
Lager Image
Transmittance versus wavelength for the proposedY-splitter.
specifies the radii of the holes. Removing a hole can be done by assigning a zero radius to it. For case (V) in Fig. 7 the radius of the junction hole (r 1 ) is assumed 0.3a,as the rest of lattice and the rest of additional holes do not exist. A simple junction is obtained in this case. As can be seen in Fig. 7 the worst transmission happens for the simple junction. When the radius of junction hole is increased from 0.3 a to 0.4 a the transmission improves (Case I).When an extra hole is added beside the junction hole (r 2 =0.31 a ), the transmittance curve improves further (Case II).Adding the second and third holes in the tapered structure also improves the transmittance further (Cases III and IV).
Finally in Fig. 8 the transmittance curve is depicted in dB scale versus wavelength to provide a more physical insight on the performance of the Y-splitter.
A PC Y-splitter for TE modes was proposed in this paper.FDTD method was used to analyze the proposed structure.Simulation results confirm that in comparison to the most recent structures reported in literature, the proposed Y-splitter provides an improved transmittance and bandwidth, while it has a much simpler design. For a bandwidth more than 120 nm the transmittance curve has less than 0.5 dB fluctuations.
Kaatuzian H 2007 Photonics (II) Tehran Polytechnic Press Persian edition Tehran Iran
Noda S , Baba T Roadmap on Photonic Crystal Kluwer Academic Publisher Norwell USA
Danaie M , Attari A. R , Mirsalehi M. M , Naseh S “Neuro-Fuzzy optimization of photonic crystal structures” in Proc. IEEE International Conference on Computer as aTool (IEEE EUROCON'07) (Warsaw Poland Sept. 2007)
Danaie M , Attari A. R , Mirsalehi M. M , Naseh S 2008 Optimization of two-dimensional photonic crystal waveguides for TE and TM polarizations Opt. Appl. 38 643 - 655
Song G. H , Kim S , Hwang K.-H 2002 FDTD simulation of photonic-crystal lasers and their relaxation oscillation J. Opt. Soc. Korea 6 87 - 95    DOI : 10.3807/JOSK.2002.6.3.087
Ryu H.-Y , Kim S.-H , Kwon S.-H , Park H.-G , Lee Y.-H 2002 Low-threshold photonic crystal lasers from InGaAsP free-standing slab structures J. Opt. Soc. Korea 6 59 - 71    DOI : 10.3807/JOSK.2002.6.3.059
Shinya A , Takara H , Kawanishi S 2008 All-optical flip-flop circuit composed of coupled two-port resonant tunneling filter in two-dimensional photonic crystal slab Opt. Express 14 1230 - 1235
Park Q.-H 2002 Nonlinear localized modes in photonic crystals J. Opt. Soc. Korea 6 96 - 99    DOI : 10.3807/JOSK.2002.6.3.096
Johnson S. G , Villeneuve P. R , Fan S , Joannopoulos J. D 2000 Linear waveguides in photonic-crystal slabs Phys. Rev.B 62 8212 - 8222
Chutinan A , Noda S 2000 Waveguide and waveguide bend in two dimensional photonic crystal slabs Phys. Rev. B 63 4488 - 4492
Chutinan A , Okano M , Noda S 2002 Wider bandwidth with high transmission through waveguide bends in two-dimensional photonic crystal slabs Phys. Lett. 80 1698 - 1700
Strasser P , Stark G , Robin F , Erni D , Rauscher K , Wüest R , Jackel H 2008 Optimization of a 60° waveguide bend in InP-based 2D planar photonic crystals J. Opt. Soc.Am. B 25 67 - 73    DOI : 10.1364/JOSAB.25.000067
Lavrinenko A. V , Tetu A , Frandsen L. H , Fage-pedersen J , Borel P. I 2007 Optimization of photonic crystal 60 waveguide bends for broadband and slow-light transmission Appl.Phys. B 87 53 - 56    DOI : 10.1007/s00340-006-2507-4
Danaie M , Attari A. R , Mirsalehi M. M , Naseh S 2008 Design of a high efficiency wide-band 60º bend for TE polarization Photon. Nanostruct: Fundam. 6 188 - 193    DOI : 10.1016/j.photonics.2008.08.003
Wilson R , Karle T. J , Moerman I , Krauss T. F 2003 Efficient photonic crystal Y-junctions J. Opt. A: Pure Appl. 5 S76 - S80    DOI : 10.1088/1464-4258/5/4/358
Frandsen L. H , Borel P. I , Zhuang Y. X , Harpøth A , Thorhauge M , Kristensen M 2004 Ultralow-loss 3-dB photonic crystal waveguide splitter Opt. Lett. 29 1623 - 1625    DOI : 10.1364/OL.29.001623
Yang W , Chen X , Shi X , Lu W 2010 Design of high transmission Y-junction in photonic crystal waveguides Physica B 405 1832 - 1835    DOI : 10.1016/j.physb.2010.01.056
Ayre M , Karle T. J , Wu L , Davies T , Krauss T. F 2005 Experimental verification of numerically optimized photonic crystal injector Y-Splitter and bend IEEE J. On Select.Areas in Commun. 23 1390 - 1395    DOI : 10.1109/JSAC.2005.851169
Kaatuzian H , Danaie M , Foghani S “Design of a highefficiency wide-band 60 degree Y-branch for TE polarization” in Proc. Opto-electronics and Communications Conference(OECC) (Hong Kong Jul. 2009) 1 - 2
Ghaffari A , Djavid M , Abrishamian M. S 2009 Bi-periodic photonic crystal Y-splitter Physica E 41 1495 - 1499    DOI : 10.1016/j.physe.2009.04.025
Povinelli M. L , Johnson S. G , Fan S , Joannopoulos J. D 2001 Emulation of two-dimensional photonic crystal defect modes in a photonic crystal with a three-dimensional photonic band gap Phys. Rev. B 64
Liu Q , Ouyang Z , Wu C. J , Liu C. P , Wang J. C 2008 All-optical half adder based on cross structures in two-dimensional photonic crystals Opt. Express 16 18992 - 19000    DOI : 10.1364/OE.16.018992
Zhang Y , Li B 2007 Ultracompact waveguide bends with simple topology in two-dimensional photonic crystal slabs for optical communication wavelengths Opt. Lett. 32 787 - 789    DOI : 10.1364/OL.32.000787
Danaie M , Kaatuzian H 2011 Design of a photonic crystal differential phase comparator for a Mach?Zehnder switch J. Opt. 13    DOI : 10.1088/2040-8978/13/1/015504
Danaie M , Kaatuzian H 2011 Improvement of power coupling in a nonlinear photonic crystal directional coupler switch Photon. Nanostruct: Fundam. Appl. 9 70 - 81    DOI : 10.1016/j.photonics.2010.10.002
Borel P. I , Frandsen L. H , Harpøth A , Kristensen M , Jensen J.S , Sigmund O 2005 Topology optimised broadband photonic crystal Y-splitter Electron. Lett. 41 69 - 71    DOI : 10.1049/el:20057717
Tetu A , Kristensen M , Frandsen L , Harpøth A , Borel P , Jensen J , Sigmund O 2005 Broadband topology-optimized photonic crystal components for both TE and TM polarizations Opt. Express 13 8606 - 8611    DOI : 10.1364/OPEX.13.008606