TELKOM NIKA , Vol.13, No .3, Septembe r 2015, pp. 8 20~827   ISSN: 1693-6 930,  accredited  A  by DIKTI, De cree No: 58/DIK T I/Kep/2013   DOI :  10.12928/TELKOMNIKA.v13i3.1733    820      Re cei v ed Ma rch 1 1 , 2015;  Re vised J une  3, 2015; Accepted June 1 5 , 2015   Simple, Easy-use and Low-cost Software for Design of  Single and Cascaded Microring Resonators Using  Semi-numerical Method       Budi Muly an ti 1* , Lilik Has a nah 2 , Tommi Hariy a di 1 , Arjuni B Pantja w a ti 1 ,   Heru Yu w o n o 3 , P.  Susthita Menon 4 , Sahbudin Sha a ri 4    1 Departme n t of Electrical En gi neer ing E ducat ion   2 Departme n t of Ph y s ics E duc ation, Un iv ersit a s Pend idik an  Indon esi a  (UPI),  Jala n Dr. Setia bud hi 2 07, Ban dun g 40 154, In don esia, T e lp/fax: 0 22 2 013 1 63/02 2 20 11 57 3 Arsari Group, Jala n Jen d . Sudirman, Ka v 1 0 - 11, Jakarta Se latan, Ind ones i a   4 Institute of Microen gin eeri ng  and  N a n oel ectronics (IMEN), Univers i ti Keb a ngsa an Mal a ys ia (UKM),   436 00 UKM Ba ngi, Sel a n gor, Mala ysi a   *Corres p o ndi n g  author, e-ma i l : b_mul ya nti@ ya ho o.com       A b st r a ct   Devel o p m ent o f  a simpl e , eas y and  low - cost softw are for desig nin g  of w a vegu ide-c o u p le d sin g l e   and casc ad ed  micror ing r e s onator (MR Rs) using se mi -n umerica l  calc u l atio n bas ed o n  transfer matrix   m e th od (TMM), is presented i n  this pap er. The softw are uses a devic e mode l w h ich is e m b e d ded  on t h e   hig h  i n d e x c o ntrast (HIC) st ructure  of sil i c on-o n -ins ulat or  (SOI) w i th mono moda l c a vi ty for T E -mod e   pol ari z a t i ons, o perati ng  arou n d  15 50  n m  o p tical w a ve le ngt h. T he  ma in  ai of the softw are is to  esti mate   the microri ng r e son a tor p e rfo r ma nce  para m eters, such as  free spectra l  ra nge (F SR) a n d  qua lity factor ( Q - factor). T he so ftw are is very  simple  a nd  ea sy to use.  W i t h  a sta n d a rd  l aptop  co mp ute r , it only  takes  few   secon d s to o b tain tra n s m is sion r e spo n se , F S R and Q - factor of sin g l e MRR  for v a rie d  w a veg u i des   separ ation d i stance a nd rin g  radi us. T he results w e re  then verifie d  usin g simulati on  meth od bas ed on fi ni t e   integr ation tec hni que  usin g 3 D  electro m ag n e tic si mu l a tor,  w h ich ne ed a  hig h  me mory  and  process o r  of   computer  and  take days to  e x ecute the s i mulati on.  W e  fo und  only s m all  discrep ancy,  w h ich in  avera ges   are ab out 4.25 % and 1 0 .80 %  for F S R and Q-factor, respec tively. In gene ral, the results obtai ne d from thi s   softw are are closer to 3D el ec troma g n e tic si mu lati on resu lts.     Ke y w ords : mi croring, se mi-n umerica l  meth od, matrix  tran sfer, free spectral ran ge, qu ali t y factor         Copy right  ©  2015 Un ive r sita s Ah mad  Dah l an . All rig h t s r ese rved .       1. Introduc tion  Microrin g resonators  (M RRs) i s  a ve ry  important  device th at i s  bei ng d e velope nowaday s be cau s of thei r ability to suppo rt the  cr eation of hi g h  den sity integrate d   circu i ts.  Variou s ap pli c ation s  of M R Rs have  be en develo ped , such  as  wa velength filter [1], multiplexing  [2], sen s ors [ 3 ], modulatio n wave, [4]  a nd bio m edi ca l [5] appli c ati ons F o r se n s ing appli c ati on,  there  are ma ny types  of sensor  th at ha ve bee n d e ve loped  such a s  fo Wirel e ss Sensor Network  (WS N) a ppli c ation [6, 7].  MRRs h a se veral advanta ges that of in terest to ma n y  rese arche r s in   this field, that are (1) its  smaller a nd co mpact  size, (2) can be int egrate d  with  other  comp o nents  (light so urce,  detecto r, cou p ler, and m a n y  others), and  (3) very ea sy  to mass p r od uce d  [8].    For de signi n g  single a n d  casca ded  MRRs a s  well as othe r MRRs config uration s several num e r ical m e thod s have been  develop ed, such a s  tra n sf er matrix met hod (T MM) [ 9 ],  finite-differenc e time-domain (FDT D) [1 0], confo r mal  tran sform a tion metho d  [1 1] or mo delin g in  terms of se mi-an a lytic couple d -m ode  theor y (CMT) [12-15], with each  method ha s its  advantag es  a nd di sadva n tage s. The  CMT in  spa c and time  dom ains i s  im ple m ented  usual ly in  two  spa c di mensi o n s  (2 D)  and  ra rel y  in 3D [16]. The  FDT D   method i s  th e mo st po pu lar  simulatio n  to ol for that p u rpo s e, b u t due to  its i n here n t we aknesse s (stai r  ca sing  error and   nume r ical di spersion ), oth e solutio n are l o o k ed  fo r. The  3 D  F D TD an alysi s  of an  add -d rop  filter co nfigu r ation, de scrib ed in  the  stu d y of  Fujii  et al. [17], illu strates high  requireme nts  in  terms of com puter resources.   On  the   othe r hand, TMM h a s sh own  results  t hat are quite  reliabl e for  mo delin g optical   micro-re son a tors an d ju st take s mi nut es in  t he  ca lculatio n [18] . TMM is a  semi -num eri c al  Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM NIKA   ISSN:  1693-6 930       Sim p le, Easy-use and L o w-co st Softwa r e for De sign  of Single and  Ca scade d… (Budi Mulyanti )   821 method that  has  bee n wid e ly kno w n fo r solving  pro b l e ms of M a th ematics an Physics, nam ely  dividing the total system i n to  N su b- sy st em s,  in wh ich a s ub- sy st em only  int e ra ct s wit h  it neigh bors. In the early 1970 s, analytical me tho d  wa s develop ed by Yariv  kno w n a s  CMT   (co uple d  mo de theo ry) u s e s  a mat r ix analytic e n e r gy-cou pling  [19]. Based  on this  wo rk,  it is  possibl e now to formulate nume r ical sol u tions in  ord e r to simulat e  MRRs prop erties  with much   less comp uting co st  co mpared  to FDT D  si mul a tions de scri bed a bove  with  comp arabl e   ac cur a cy .   In this study, we propo se d a simple, ea sy- use a nd lo w co st softwa r e develo ped  usin g a  semi -num eri c al method, n a mely the tra n sfer m a trix method (TM M ) to estimat e  the values of  MRRs’ prope rties, namely FSR  a nd  th e Q-fa ctor.  T h e  re sults a r e th en verifie d  by  a F D T D ’s m o st   popul ar co m m ercial software   CST   Microwave  Studio ,  usin g 3 D   el ectro m ag neti c   simulato b a se d   on finite integration te chn i que, whi c need s a  high  memory an d pro c e s sor  of compute r  and  take s days to  execute the  simulatio n .       2. Rese arch  Metho d   2.1.   Dev i ce Model  MRRs d e vice  model u s e d  in this  study, wa s ba sed  o n  SOI stru ctu r e with  high  refractive   index contra st and op erati ng at  15 50 n m  of wavele ngth a s  sch e m atically sho w n in Fi gu re  1  whi c h co mpri se s of a ring waveg u ide cl osely co uple d  to double straight bus  wa veguide s (Fig ure   1(a )) an d th e cross secti on of the waveguid e   structure is  sh own in Fig u re 1(b ) . The  bus  waveg u ide s   serve a s   evan escent li ght i nput a nd  out put coupl ers,  while  the  ring  wave guide  a c ts  as th e wavel ength  sele cti v e eleme n t.  d epict s the  ring  ra diu s gap  i s  the  se paratio n di sta n ce  betwe en  strai ght and  rin g   waveg u ide,  i s  the  wav eguid e  wi dth  and  i s  the  total waveg u i de  height. Th e f u lly etch ed  waveguid e   structure i s   con s ide r ed  throu ghout thi s   st udy. The  cro s s- se ction of the  MRRs d e vice is 0.30  μ m x  0.55  μ m, on  top of 1 μ m-t h ick bu ried  o x ide (BOX) l a yer .   The  cro s s-se ction of  su ch  waveg u ide s  wa cho s en  to ensure si ngle-mod e  propag ation n e a 1550 nm telecommuni catio n s wavelen g ths [20].          Figure 1.   (a)  Lay out and (b) cro ss  se ction of singl  MRR  with do uble straight  waveg u ide s       2.2.   Calcu t ion of  MRR c h ara c teristics usin g CMT and T MM  Effec t ive refrac trive index,    , can be d e ri ved usin g [19 ]                ( 1 )     Whe r β  is p r opa gation  co nstant. By kn owin g   β  and  the operatin g wavele ngth ,   λ , the value of    can be o b tai ned. The  cou p ling coeffien t can be calculated u s ing  Figure 2.   Evaluation Warning : The document was created with Spire.PDF for Python.
                          ISSN: 16 93-6 930   TELKOM NIKA   Vol. 13, No. 3, September 20 15 :  820 – 827   822   Figure 2. Cou p ling effect b e twee n the st raight waveg u ide an d the ring with ra diu s  R      In Figure 2, it can be  see n   that the cou p ling len g th,  , betwee n  the  straight a nd  ring   waveg u ide s , i n  which the  in put si gnal  of  straig ht   wavel ength i s  co up led  with the  ri ng  wavele ngt h.  The amo unt  of coupli ng is  determi ned b y   the magnitu de of the co u p ling coefficie n t,  k . Therefo r e,  the cou p ling l ength,  , can b e  cal c ulate d  usin g [21]:                                                                                         (2)    Whe r λ  i s  operatin g wav e length,    and    are effective refra c tive ind e x for even and odd   wavele ngth, resp ectively. And by using:                                             (3)    The value of  cou p ling  coef ficient,  , can be cal c ul ated     Figure 3. MRRs  config urati on and re p r e s ente d  co upli ng matrix P and Q       MRRs  config uration  ca n b e  forme d  by  one o r  two  st raight  waveg u ide s  and  on e or m o re   than one ri ng  waveguid e s.  The input si gnal ente r s o ne end of a straig ht wave guide an d go es  forth at the  ot her  end.  Du ri ng its travel i n  a  st rai ght waveguid e mo st  si gnal   will be cou p led a nd  towards the  ring  wavegui d e . If there i s  more than  o ne ri ng, the   sign al will  be  co upled  ag a i n   towards the n e xt ring, or it coul d be cou p led to ward anothe r strai ght waveg u id e. In Figure 3 ,   is a s sumed t o  be a  cou p l i ng matrix be tween th e straight and  rin g  waveg u ide s  or  a coupli n g   matrix bet we en the  two  ri ng  waveg u id es,  while  Q i s  a    cou p ling  matrix of  sig nal p r op agati o n   along  th e wa veguide ring, then  the  tran sfer equ at ion  for MRR  ca be exp r e s sed  in the follo wi ng  equatio n [19]:          , | | | | 1       ( 4 )     Whe r κ  is a  norm a lized couplin g co efficient an d  t is si gnal tra n smitted in the ring  waveg u i de.     Furthe rmo r e,  the cou p ling  matrix can b e  obtained u s i ng [19]:        ,  1 1         ( 5 )     Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM NIKA   ISSN:  1693-6 930       Sim p le, Easy-use and L o w-co st Softwa r e for De sign  of Single and  Ca scade d… (Budi Mulyanti )   823 And pro pag ation matrix  Q is [19,20]:      , 0      0       ( 6 )     Whe r e R i s  ri ng radi us,  an β  is pro pag ation co nsta n t  mentioned i n  Equation (3 ).   Since FSR i s   defined a s  a  distan ce of two  adja c ent pe ak inten s ity [19], therefore:      2               ( 7 )     Whe r e L i s  a  circumfe ren c e of micro r in g circle a nd  is spe ed of li ght in a medi um that can  be  expre s sed a s                ( 8 )     Whe r e    is sp eed of light in  the vacuum.   Based o n  Fig u re 4 an d by usin g bou nda ry conditio n  o f  [19]:    ∝            ( 9 )     Whe r α  i s  a b so rption  co efficient pe unit length,  a c cordi ngly th e tran smitan ce in  the  straight  waveg u ide is  [19]:       | |  | |    | |  | |          ( 1 0 )         Figure 4. Ligh t wave tran smitted in the straig ht wave guide       And the total power in the  microrin g is [19]:    | |  α  | |  | |  | |           ( 1 1 )     If    , then fro m  equatio n (10), full widt h at half ma ximum  (FWHM)  can  be   obtaine d from     | |   ; | | 1       ( 1 2 )      And the Q-fa ctor  can be e x presse d as:          | |   | |         ( 1 3 )       3. Results a nd Discu ssi on  The software develop ed in this study ca n be  used to study device perfo rman ce  su ch a s   freque ncy  re spo n se, FSR, the Q-fa cto r , finesse, an d FWHM fo r singl e a s  well as  ca sca ded  Evaluation Warning : The document was created with Spire.PDF for Python.
                          ISSN: 16 93-6 930   TELKOM NIKA   Vol. 13, No. 3, September 20 15 :  820 – 827   824 MRR  devices, was  re sulte d  in this  stud y. The GUI  p age of the  sof t ware  ca n be  see n  in Fig u re 5.  All of the de sign pa ram e te rs,  su ch  as  ri ng radiu s , ga p, width  of waveguid e , ref r active i ndex es  can  be e a sil y  incorpo r ate d  into the G U I pre c e ded  by cho o si ng  the MRR co nfiguratio n.  The   results  of the  simul a tion, n a mely fre que ncy re spon se , FSR, the Q - facto r , fine sse, and  FWHM   will be soon  shown in the GUI pag e.           Figure 5. The  GUI page of  softwa r e resu lted in this stu d y.      3.1. Single MRR  The sim u latio n s were d one  by varying ring  radi us a n d  wavegui de s eparation di stance of   4-12 µm an d 100-150 nm,  respe c tively at incident  wav e length of 1.5 5  µm. The other pa ram e te rs  are  ke pt co n s tant, incl udi ng SiO 2  an  Si  refra c tive  index of 1.5 277 a nd 3.4 7 77, re sp ectiv e ly,    SiO 2  and Si  height of 1.0 and 0.55 µm,  resp ective ly and the wave guide wi dth o f  0.30 µm. T he  results were  then  verifie d  by  re sults  of sim u lation  usi n g  sim u l a tion m e thod  ba sed  o n  fi nite  integratio n te chni que  usi n g 3D ele c tro m agneti c   sim u lator,  whi c h  need a hig h  memo ry a nd  pro c e s sor  of comp uter  an d take s d a ys to execute  the si mulation .   The fre que ncy pe rforma nce   for  singl e micro r in g for dif f erent value s  of ring radi u s  and 1 00 n m  of  waveg u ide sepa rati on   distan ce, is  shown in Figu re 6.          Figure 6. Single MRR tran smissio n  re sp one of the drop port for dif f erent value s  of ring ra diu s , R  = 6 and 1 0  µm at wavegui de se paration  distan ce, wa veguide  width ,  and waveg u ide heig h t are   100nm, 0.30  µm and 0.55  µm, resp ectiv e ly      Figure 6  sho w s pe aks  of tran smi ssi on i ndicat ed po wer  in   resona n c e co ndition. It  can be   see n   clea rly the p o si sitio n  of two a d j a ce nt pe aks that can  be  tra c ed  at di fferent valu e  of  wavele ngth  d epen ds on  ri ng  radiu s . Su ppo se fo 6 µ m  of  ring  ra di us  and  10 n m  of  sep a ration  Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM NIKA   ISSN:  1693-6 9 30       Sim p le, Easy-use and L o w-co st Softwa r e for De si gn  of Single and  Ca scade d … (Budi Mulyanti )   825 distan ce, resonant pe aks  occur  at  1.55 20 and  1.566 3 of wavele n g th, therefo r e  the FSR valu e of  those p a ra me ters i s  14.29n m.   To comp are   our  TMM  si mulation  re sults, we  cal c ulate the  fre quen cy resp onse u s in g   Finite Integration Te chni que (FIT) u s ing 3D  ele c troma gneti c  simulator  with  the sam e  ring   para m eters and config urations. We use   co mme rc ial softwa r e   from Com p uter  Simul a tion   Tech nolo g y AG named  CST Studio Suite to perform   3D FIT sim u lation. From  simulatio n  using   CST Studio  Suite, we o b t ained FS value at 1 5 .4 0nm for rin g  radi us  6 µm  and  100  nm  of  sep a ratio n  di stan ce. We can se e that the sm aller  th e radi us of th e ring, the g r eater the di st ance  betwe en two  adja c ent re so nant pea ks.   From  our T MM sim u lati on an d CST  simul a tion,  the high est  FSR obtai n e d  at the  smalle st rin g   radiu s , therefore the valu e  of FSR  is in versely p r op o r tional to the  effective leng th   microrin g, L  (= 2 π R). On  the ave r ag e, the differen c e  between  se mi-num eri c al   simulatio n s with  and  CST  si mulation s i s   4.25%, a s   shown in  Fi g u re   7. We  may  con c lud e   that only small   discre pan cy betwe en both  simulation s.           Figure 7. FRS values of Single MRR for different valus of ring radiu s       3.2. Double  and Triple Serially  Casca ded MRR    The analy s is  wa s also don e with dou ble  and triple se rially ca scad ed microri ng  for the  same  value s   of sep a ration  distan ce  an d rin g  r adiu s ,  i.e:  0.1 µm  and 6  µm, re spe c tively. The   device  perfo rmance of the  doubl e an d triple  ca scade d microri n g s   are  simul a ted  and  comp ared   to a si ngle  m i cro r in g an d t he results  are interprete in Figu re  8.  Two  re son a n t  pea ks  occu r at  1.5520  and  1 . 5663 givin g   FSR eq ual to  14.5 nm  reg a rd le ss the n u mbe r  of mi crori n g s  orde r. To  verify the results ,  we c onfirmed the sim u lations results with CST si mulation.       Figure 8. The  FSR values f o r differe nt for the si ngl e, double a nd trip le ca scade d microrin gs at  ring radiu s  an d sep a ratio n  distan ce of 6  µm and 100 n m  respe c tively.       Furthe rmo r e,  to study the perfo rman ce  of mi cro r in g device, the Q - facto r  wa s e s timate d   usin semi -n umeri c al  sim u lation. Fig u re 9  sh ows th e plot s of  Q-f a ctor to  se pa ration  dista n ce i n   the rang e of 0.1  μ m and 0.15  μ m at con s tant ring  radi us of 6  μ m.   Evaluation Warning : The document was created with Spire.PDF for Python.
                          ISSN: 16 93-6 9 30   TELKOM NIKA   Vol. 13, No. 3, September 20 15 :  820 – 827   826   Figure 9. The  plots of Q-fa ctor to s epa ra tion distan ce i n  the rang e o f  0.1  μ m and 0.15  μ m at  con s tant rin g  radiu s  of 6  μ m of semi-nu m eri c al an d CST sim u lations.         It is shown that greate r  Q - facto r  occu rs at  la r g er  sep a r a tion  d i s t an c e . T h e  ob se r v a b le  Q-fa ctor diffe ren c e  in th e  se pa ration   distan ce  of  0.1  μ m and   0.15   μ m i s  375  a n 6 10,  respe c tively. Clea rly, the gap si ze h a s a  signi fi cant effect o n  the Q-fa ct or, yet the FSR   unaffecte d . Despite th e lo w lo ss  pro d u c ed  at the na rro w g ap, it is difficult to f abri c ate a n fine  etchin g i s  re quire d. Fu rth e rmo r e,  we   found  a mi smatch  in th e Q - facto r  v a lue s  for bo th   simulatio n s.  The la rge s t di fference was  found for th e gap  size  of 1 2 0 nm  with 2 0 .5% differen c e.  This ag ain  du e to  different   in effective  re fracti ve i ndex  value s   wh ere the r e i s  no  straig htforward  way of  cal c ulating the   effective ref r active in d e x, espe cially i n  be nd  wav eguid e. The r eby,  approximatio ns a r con s i dere d  in the  MRR  model i n g. Ho weve r, both re sult s experie nced  the  same  tren d,  whi c h the  Q-f a ctor rai s e d   as the   sepa ration  ga ps wi den s.  The ob serva b le Q-fa ctor  differen c e a c ross the sepa ration gap of  0.12  μ m is 1 0 , which cont ribute to 2.6 %  in percent age,  while the ave r age of deviat i on is 10.8 0 %, which  can b e  con s id ere d  relatively low.       4. Conclusio n   The sim p le,  easy-use an d low-cost so ftware for  de sign of  single  and serial  cascad ed  microrin g re sonators  (M RRs) usin th semi-num e r ical  meth od  ha s b een  d e velope d in  this  study. Th so ftware  can  be  used to  e s ti mate th e  M R Rs pe rforman c pa ramete rs, such  as F S R   and Q - fa ctor in only fe w se con d s u s ing a  stan da rd la ptop  co mputer. T h e  re sults  of the  simulatio n   were th en ve ri fied u s ing  si mulation  met hod  ba sed  o n  finite inte g r ation te ch ni que   usin g 3D el e c trom agn etic simulato r. In gene ral,  we found only  small di scre pan cy, which  in   averag es a r about 4.25%  and 10.8 0 % for  FSR a nd  Q-fa ctor, re sp ectively.       Referen ces   [1]  Kazmiercz a k A, Boaerts W ,   T hourho ut DV, Dr ouard E,  Rojo-R ome o  P, G i annon D, G a ffiot F.  Analys is of Si licon  on Ins u l a tor (SO I ) O p tical Micror ing  Add-Dro p  F i l t er Based  on  W a vegu id e   Intersection.  Pr ocee din g s of SPIE. Strasbourg. 2008; 6 996.    [2]  Palocz i G T Hua ng Y, Yar i v A. F r ee-sta ndi ng  all P o l y mer Microri n g  Reso nator O p tical F i lter.   Electron ics Let ter.  2003; 39( 2 3 ): 1650- 16 65.   [3]  Xi ao S, K h a n   SM, Shen  H,  Q i  M. Multipl e -chan nel  sil i co n micro-r e son a tor bas ed  filters for W D M   app licati ons.  Optics Express.  200 7; 15(1 2 ): 7489- 749 8.   [4]  Cha o  CY, G uo LJ. Desig n  and O p timiz a tion of  Micro r ing Res o n a to rs in Bioch e mical Se nsin g   Appl icatio ns.  Journ a l of Li ght w a ve T e chnol o g y.  2006; 2 4 (3) :  1395-1 4 0 2 [5]  Cha o  CY, F u n g  W ,  G uo LJ. Pol y mer Micr orin g Res onat ors for Bioch e m ical Se nsi n g  Applic atio ns.   IEEE Journal  o f  Selected Top i cs in Quantu m   Electron ics.  20 06; 12(1): 1 34- 142   [6]  Hari ya w a n  M Y , G una w a n   A, Putra EH.  W i reless  Se nsor  Net w ork  for F o rest F i re D e tectio n.   T E LKO M NIKA Indon esi an Jou r nal of Electric al Eng i ne eri n g .  2013; 1 1 (3): 5 63-5 74.   [7]  Hase eb K, Ab u Bakar K, Ab dull ah A H , Ah med A. G r id B a sed  Cluster  H ead S e lecti on  Mecha n ism for   W i reless s ens or net w o rk.  T E LKO M NIKA T e leco mmunic a tion  Co mp uti ng El ectron ics and  Co ntrol 201 5; 13(1): 26 9-27 6.  [8]  F i dan bo yl u K,  Efendi ogl u HS.   F i ber O p tic S ensors  and T h eir Ap plic ation s .  5 th  Internatio nal A d vanc ed   T e chnolog ies Simposi u m.  20 10.   Evaluation Warning : The document was created with Spire.PDF for Python.
TELKOM NIKA   ISSN:  1693-6 930       Sim p le, Easy-use and L o w-co st Softwa r e for De sign  of Single and  Ca scade d… (Budi Mulyanti )   827 [9]  Po w e ll O. Sin g le-mo de  Co n d itio n for Si lic on R i b Wav e g u id es J.  Li ght w a ve T e chn o l .  200 2; 20( 10):  185 1-18 55.   [10]  Little BE, Ch u  ST , Haus HA . F o resi, Lai ne  J- P. Microring  Reso nator C h ann el Dr op pin g  F ilter.  J.of   Lig h tw ave T e chno l . 199 7; 15( 6): 998-1 0 0 5 .   [11]  Belar ouci A, H i ll KB, Li u Y,  Xi ong Y, C h a ng  T ,  Craig AE. D e sig n   an d Mo d e lin g of W a ve g u id e-Co upl ed   Microrin g Res o nator.  J. Lumin .  2001; 94( 95): 35-3 8 [12]  Hee bner J, Gr over R, Ibra hi m T .  Optical Micr oreso nator s - T heor y ,  A pplic atio ns a n d  F abric atio n.   Lon do n: Spring er-Verla g Lo nd on Lim i ted, 20 08.   [13]  Boga erts W ,  De He yn  P, Va n  Vaere n b e rgh   T ,   DeVos K, S e lvar aja S K , Cl aes T ,  Dumon  T ,  Bienstman   P, Van T hourhout D, Baets R.  Silico n  microri ng reso nators.  Laser Ph oton ic s Rev . 2012; 6:  47-73.   [14]  Benso n  T ,  Boriskina  S, Se w e ll  P, V u kovic  A, Greedy  S,  Nosic h  A.  Mi cro-Optical Re sonators  f o r   Microlas ers  an d Intergr a ted   Optoelectro n ic s-Rece nt Adva nces  an d F u tur e  C hal len ges . I n  F r onti e rs o f   Plan ar Li ght w a ve Circu it T e chno log y : D e si g n , Simu l a tion   and F a bricati o n. Sprin ger, L ond on. 2 005:  39-7 0 [15]  Hiremath  K,  Hammer M.  Circul a r Inte gr ated Optic a Microres o nator s-Analytic al M e thods  a n d   Co mp utation a l  Aspects . In Photon ic Micror eson ator Res e arch  an d App l i c ations. Spri ng er Series i n   Optical Scie nc es. Spring er Verla g , Lon don.  201 0; 156: 29- 59.   [16]  Stoffer R, Hire m ath KR, Ham m er  M, Prkna L, Ct y r ok J. C y l i ndric al  i n tegr ated o p tical micr oreso nators  -   mode lin g b y   3- D vectoria l cou p le d mode th e o r y Opt. Comm un . 200 5; 25 6: 46-67.   [17]  Fujii M, Fre u de W, Russ er  P.  Efficient  High-S pati a l-O r der F D T D  A nalysis  of  3D  Optical  Ri ng   Reso nator F ilt ers . In Proc. 19th Ann ual R e v. Progress in  Appl. Com put. Electromag ne tics. Montre y ,   CA. 2003: 7 39- 744.   [18]  Masi M. Mod e lin g of S equ ences  of Silic on Micr o-Res onators for O n -Chi p Optica l  Routi ng  a n d   S w itc h i ng.  PhD  T hesis. T r ent o: Univ of  T r ento. 2011.   [19]  Yariv A.  Cou p l ed-Mo de T heo r y  for  Gui ded- W a ve Optics.  I EEE Jour nal  of  Quantu m  El ec tron.  19 73;  9:   919- 933.   [20]  Mul y a n ti B, Hari ya di T ,  Hasana h L, Hazu r a  H, Sahb udi n  S, Menon PS , Abang AE, Mahmu d in D ,   W i ranto G. De sign  and Optim i zatio n  of Co up led Micr orin g R e son a tors (MR Rs) in Si lico n -o n-Insul a tor.   Journ a l Sai n s Malaysi an a.  20 14; 43(2): 2 47- 252   [21]  Palocz i GT , H uan g Y, Yariv A, Mookherj e a  S.  Pol y mer i c Mach-Z eh nderI n terferometer  Using S e ria l l y   Cou p le d Micror ing R e son a tors Optics Express . 2003; 11( 21) : 2666-2 6 7 1 [22]  Mookh e rje a  S. Optical D e vic e s for Com p u t ing:  Mod e li ng  Photon ic Circ u its. Lecture.  San D i eg o :   Univers i t y   of Califor nia. 20 12.                Evaluation Warning : The document was created with Spire.PDF for Python.