In this paper, four apodization functions are proposed for silicononinsulator (SOI) strip waveguides with sidewallcorrugated gratings. The effects of apodization functions on the full width at half maximum (FWHM), the sidelobe level, and the reflectivity of the reflection spectrum are studied using the coupledmode theory (CMT) and the transfermatrix method (TMM). The results show that applying proposed apodization functions creates very good filtering characteristics. Among investigated apodized waveguides, the apodization functions of
Polynomial
and
zpower
have the best performance in reducing sidelobes, where the sidelobe oscillations are entirely removed. Four functions are also used for precise adjustment of the bandwidth. Simulation results show that the minimum and maximum values of the FWHM are 0.74 nm and 8.48 nm respectively. In some investigated functions, changing the apodization parameters decreases the reflectivity which is compensated by increasing the grating length.
I. INTRODUCTION
Bragg gratings are used in many integrated optical components for a large variety of applications including optical signal processing, highspeed optical communications, sensing systems and networking
[1

4]
. Examples of prevalent gratingbased components are add/drop filters for wavelengthdivisionmultiplexing (WDM) communication systems, gratingassisted couplers, dispersion compensators, distributedfeedback and distributedBraggreflector lasers
[5

8]
.
Great compatibility with CMOS structures is made by silicononinsulator (SOI) technology, a promising candidate for largescale integration of optics and electronics on a single silicon platform
[5
,
7
,
9]
. Therefore, in the past decades, SOI waveguides have been attractive choices for integration with Bragg grating structures
[8
,
10
,
11]
. A main group of SOI structures are strip waveguides which are used in many applications, including WDM add/drop filters and grating couplers
[4
,
12
,
13]
. In strip waveguides, or photonic wires, usually submicron cross sections are used
[9
,
10]
. Moreover, the strong confinement of light in the core due to the high index contrast between Si and SiO
_{2}
is a beneficial feature of strip waveguides. So, as a result of small waveguide dimensions and high confinement of light, even small perturbations on the sidewalls can result in high coupling strength
[10]
.
Various methods are suggested for implementing Bragg gratings on the SOI waveguides. Two main categories are: (1) using ion implantation for refractive index modulation
[14]
, and (2): physically corrugating the top surface or the sidewalls of the waveguides
[9
,
11
,
15]
.
light rays reflects back around the Bragg wavelength when it propagates through the Bragg gratings. The Bragg condition is satisfied at the Bragg wavelength,
λ _{D}
, which is defined by
[9]
:
where
n_{eff}
and Λ are the effective refractive index of the fundamental mode of the waveguide and the grating period, respectively. However, due to the presence of some sidelobes at the vicinity of the Bragg wavelength, uniform gratings do not provide the proper performance needed in the integratedoptical components and WDM communication systems, frequently. Grating filters and multiplexers, as fundamental components in the highspeed WDM communication systems, aim to reflect a single channel in the reflection spectrum. As a result, the reduction of the sidelobes is essential for them. On the other hand, the presence of the sidelobes leads to the creation of high cross talk levels which is considered as a drawback for some applications in WDM communication systems
[6]
. Subsequently, the attendance of the sidelobes degrades the functionality of the communication systems. The reduction of the sidelobes is done by using the apodization concept. It is performed by gradually changing the grating strength and therefore the coupling coefficient along the waveguide
[6]
. By using some apodization functions, one can completely remove the sidelobe oscillations, thus resulting in smooth reflection spectra compared with reported results
[9]
for unapodized waveguide. Realization of the apodization can be fulfilled by varying grating duty cycle
[16]
or dimensions of the waveguide
[6]
.
In this paper, four apodization functions are introduced for SOI strip waveguides with sidewall corrugated gratings. Here, the analysis is done using the coupledmode theory (CMT) and the transfermatrix method (TMM). The effects of apodization functions on the sidelobe level, the full width at half maximum (FWHM), and the reflectivity of the grating reflection spectrum are investigated. Improvement of the apodization functions, in order to achieve better sidelobe reduction and filtering characteristics, is accomplished by changing the parameters of the apodization functions.
II. THEORY
Physical perturbation of the waveguide causes the refractive index modulation and brings the coupling between the forward and backward propagation modes. It occurs when the phasematching condition of Eq. (1) is satisfied or simply when
β
≡
π
/Λ .
Consider Ψ
_{1}
(
x,y
) as the transverse mode of the unperturbed waveguide that can propagate in both positive and negative directions. The electric field in the grating section can be expressed by
[19]
:
where
ω
is the angular frequency of the light,
β
is the propagation constant of the mode that is defined by
β
=2
πn_{eff}
/
λ_{0}
where
λ_{0}
is the freespace wavelength, and
F
_{1}
and
G
_{1}
are the field amplitudes of the forward and the backward propagating modes, respectively. The electric field in Eq. (2) satisfies the wave equation
[19]
:
where the relative permittivity,
ε_{r}
, is independent of the grating length while the perturbation in the relative dielectric permittivity, Δ
ε_{r}
, is a periodic function of
z
and can be expanded by Fourier series
[19]
:
Since no gain or loss is assumed,
ε_{r}
and Δ
ε_{r}
are considered to be real and so Δ
ε_{m}
= Δ
ε*_{−m}
, where Δ
ε*_{m}
is the complex conjugate of Δ
ε_{m}
[19]
. Here m =1 is considered to obtain the coupled equations
[20]
. CMT can be applied for orthogonal modes [19]. The orthogonality relation is represented by the following equation
[19]
:
where
μ
_{0}
is the permeability of the vacuum. By using Eq. (2)(5) the two coupledequations can be obtained as follows:
Here Δ
β
_{1}
≡
β
_{1}
− π /Λ and the coupling coefficient is defined as:
where
ε
_{0}
is the permittivity of the vacuum and Δ
ε
_{1}
is the firstorder Fourierexpansion coefficient of the relative dielectric perturbation, Δ
ε_{r}
. Often, grating structures can be analyzed by two methods; the direct numerical integration method and the piecewise uniform approach
[21]
. To solve coupled equations, as a twopoint boundaryvalue problem, the shooting method and TMM can be used
[21]
. The shooting method requires many iterations and therefore is time consuming. However, TMM is an accurate method for analyzing grating structures
[21]
. In this paper, TMM, developed from the piecewise uniform approach, is used to analyze the gratingsassisted waveguides. In this method, the grating with the length of
L
is divided into
N
segments and the grating specifications in each segment are considered to be uniform. For each segment, the CMT is applied and then the transfer matrices of all segments are multiplied by each other.
Equations (6) and (7) can be rewritten in the following form
[19]
:
where
S
is a 2× 2 matrix. The relation between the fields at
z
_{0}
and
z
_{1}
can be presented as
[7]
:
where
C
is the transfer matrix and
z
_{0}
and
z
_{1}
are the initial and ending positions of each segment with uniform gratings, respectively. The analytical solution of Eq. (10) can be expressed by
[7]
:
Matrix exponentials in Eq. (11) can be solved by the Pade approximation
[19]
and the matrices of
S
_{1}
and
S
_{2}
are represented by
[7]
:
Thus, the total transfer matrix can be described as follows:
By using Eq. (14), the reflectivity and the transmittivity of the grating can be represented by:
The relation between the reflectivity
R
, the coupling coefficient
κ
_{11}
, and the grating length
L
, can be expressed as
[6]
:
As it was mentioned, by using the CMT the relation between the forward and the backward propagating modes can be described through a set of equations. On the other hand, by utilizing the TMM, the grating structure is divided into N segments with uniform gratings. For each segment, applying the CMT results in a transfer matrix which can relate the fields at two ends of the segment. Therefore, overall changes in the grating (along the structure as a result of the apodization functions) can be modeled by multiplying the transfer matrices of segments.
III. SIMULATION OF SOI STRIP WAVEGUIDES WITH UNIFORM AND APODIZED GRATINGS
Typical structures of sidewall corrugated SOI strip waveguide (SWCSOISW) with uniform gratings and an exemplary structure with apodized gratings are shown in
Fig. 1
.
(a) Cross section, (b) Top views of the SOI strip waveguide with sidewall corrugated gratings and (c) an exemplary structure for apodized gratings.
The waveguide consists of a silicon layer over the surface of a buried oxide layer with the thicknesses of 220 nm and 2 μm, respectively. For all waveguides reported in this paper, the strip width, W, the grating period, Λ, and the duty cycle are considered to be 500 nm, 310 nm and 50%, respectively. The corrugation width, ΔW, and the grating length L are changed for different waveguides. The phasematching condition, given in Eq. (1), and the effective refractive index of the fundamental mode of the SWCSOISW are depicted in
Fig. 2(a)
. It was assumed that the corrugation width and the grating length are 5 nm and 620 μm, respectively, The crossing point of the two diagrams in
Fig. 2(a)
shows the resonant wavelength of the grating structure. The transmission and the reflection spectra of the waveguide are shown in
Fig. 2(b)
.
The effective index of the fundamental mode and the phasematching condition versus the wavelength and (b) Transmission and reflection spectra of the SOI strip waveguide with 5 nm uniform gratings and the length of 620 μm.
As is clear, the resonant peak of the transmission spectrum is compatible with the Bragg wavelength of Eq. (1). The response exhibits the FWHM of 3.31 nm, the high sidelobe level of 1.8 dB, and the reflectivity of 100%. Suppression of these high sidelobes can be performed by using the apodized gratings. To do this, many types of grating apodization can be utilized. The apodization functions introduced in this paper are:
Exponential 1:
Exponential 2:
Polynomial:
zpower
:
where
a, b, c, e, f, j, k, m, n, p
and
q
are the apodization parameters that may be varied to optimize the filtering characteristics of the structure and
z
′ =
z
/
L
. The apodization functions are depicted in
Fig. 3
.
Apodization functions versus the normalized length (a) Exponential 1, (b) Exponential 2, (c) Polynomial, and (d) zpower.
The schematic of the apodized waveguide with the apodization function of
Exponential 1
, and the simulation results for different apodization parameters with the corrugation width of 5 nm and the grating length of 620 μm, are shown in
Fig. 4
and summarized in
Table 1
.
(a) Schematic of the apodized waveguide with the apodization function of Exponential 1 for a=0.6 and b=c=1. Reflection spectra of the apodized waveguide for (b) b=c= 1, a=0.6 and a=0.8, (c) a=0.7, b=1, c=0.08 and c=1.3, and (d) a=c=1, b=1.2 and b=1.4.
Spectral features of the apodized waveguide with the apodization function ofExponential1
Spectral features of the apodized waveguide with the apodization function of Exponential 1
As it is evident, by varying the apodization parameters, the value of the FWHM, the sidelobe levels and the reflectivity would be changed dramatically. An inspection in the results reveals that the FWHM is lower than the unapodized or uniform waveguide. Therefore, this function can be used for applications with the requirement of the low bandwidth grating structures such as WDM filters
[3
,
10]
. By comparing the results, it is found that, increasing the parameters
a
and
c
leads to a decrease in the FWHM, the sidelobe level and the reflectivity, while increasing the parameter
b
enhances them.
On the other hand, among the reported results in
Table 1
, the lowest values of the FWHM and the sidelobe are 1.09 nm and 15.7 dB, respectively, which can be obtained for
b
=
c
=1 and
a
=0.8. However, the reflectivity is relatively low for this case. According to Eq. (17), it is possible to increase the reflectivity by increasing the grating length. The reflection spectra of the unapodized and the apodized waveguides (with apodization function of
Exponential 1
for
b
=
c
=1,
a
=0.8 and various lengths) are plotted in
Fig. 5
. These results are summarized in
Table 2
.
Comparison between the reflection spectra of the unapodized and the apodized waveguides with the apodization function of Exponential 1 for b=c= 1, a=0.8 and different grating lengths.
Spectral features of the apodized waveguide with the apodization function ofExponential 1forb=c= 1,a=0.8 and different lengths
Spectral features of the apodized waveguide with the apodization function of Exponential 1 for b=c= 1, a=0.8 and different lengths
As shown in
Table 2
, grating length increment leads to an increase in the reflectivity and the sidelobe level while the FWHM decreases and reaches below 1nm. Here, the minimum value of the FWHM is 0.92 nm which is obtained for L=1000 μm.
In
Fig. 6
the schematic and the spectra of the apodized waveguide with the apodization function of
Exponential 2
are shown and the simulation results are presented in
Table 3
in summary.
(a) Schematic of the apodized waveguide with the apodization function of Exponential 2 for e=1, f=1.6, (b) Reflection spectra of the apodized waveguide for: e=1, f=2, f=2.2 and f=2.8, and (c) Comparison between the reflection spectra of the apodized and unapodized waveguides with 2 nm and 5 nm gratings, respectively.
Spectral features of the apodized waveguide with the apodization function ofExponential 2
Spectral features of the apodized waveguide with the apodization function of Exponential 2
The simulations of apodized waveguides are performed using the structures with the grating length of 620 μm and the corrugation width of 2 nm. As it is clear, compared to the unapodized waveguide with the corrugation width of 5 nm, in some cases, the FWHM of the apodized waveguide is increased considerably. From
Table 3
, it can be seen that by enhancing the apodization fparameter, the FWHM decreases dramatically. Moreover, the sidelobe level decreases when this parameter increases from
f
=1.6 to
f
=2.8 and then increases for the higher values. From Table 3, it can be seen that the high and low values of the FWHM with the desired values of the sidelobe level and the reflectivity can be attained by adjusting the apodization parameters. Among these functions, the maximum FWHM of 8.48 nm is obtained by applying the apodization function of
Exponential 2
for
e
=1 and
f
=1.6.
The schematic of the apodized waveguide with the apodization function of
Polynomial
is plotted in
Fig. 7(a)
. The simulation results for the grating length of 620 μm and the corrugation width of 5 nm are shown in
Figs. 7(b), (c)
. These results are summarized in
Table 4
.
(a) Schematic of the apodized waveguide with the apodization function of Polynomial for j=1, k=2. Reflection spectra of the apodized waveguide for: (b) j=1, k=1/4, k=1/2, k=2 and k=4, (c) k=2, j=1 and j=1.3.
Spectral features of the apodized waveguide with the apodization function ofPolynomial
Spectral features of the apodized waveguide with the apodization function of Polynomial
As shown in
Fig. 7
, the apodized waveguide presents extremely good filtering behavior and the sidelobe oscillations of the reflection spectra are entirely removed. The calculated smooth spectra can be utilized in applications such as WDM communication systems where the presence of sidelobes is considered to be a drawback.
Fig. 7(b)
demonstrates the impact of changing the apodization kparameter on the amplitude of the sidelobes, for apodization function of
Polynomial
type. As this parameter increases, the amplitudes of the sidelobes gradually diminish and reach to zero for
k
=2, and the reflectivity is 96%. By further increasing
k
, the spectrum remains smooth, but the reflectivity is reduced. As shown in
Fig. 7(c)
, by choosing
k
=2 and changing the
j
parameter, the FWHM of the spectra can be adjusted. For low values of
j
, the FWHM becomes narrower than that of the unapodized waveguide while the higher values of
j
lead to wider spectra compared to the unapodized one. Based on the results listed in
Table 4
, it can be seen that the minimum FWHM of 1.24 nm is achieved for the parameters
k
=2 and
j
=0.6 with the completely smooth spectrum. However, in this case, the reflectivity is small. It can be increased by increasing the grating length.
The impacts of enhancing the grating length on the reflection spectra of the apodized structure with the apodization function of
Polynomial
are shown in
Fig. 8
and brief results are given in
Table 5
.
Comparison between the reflection spectra of the unapodized and the apodized waveguides with the apodization function of Polynomial for k=2, j=0.6 and different lengths.
Spectral features of the apodized waveguide with the apodization function ofPolynomialforj=0.6,k=2 and different lengths
Spectral features of the apodized waveguide with the apodization function of Polynomial for j=0.6, k=2 and different lengths
Comparing the simulation results of unapodized structure, it is clear that the completelyflat apodized spectra are much narrower than that of the unapodized one. For
L
=1550 μm, the reflectivity reaches above 90% and the FWHM becomes 0.87 nm, that is too small in comparison with the reported results of the unapodized waveguide
[9]
.
The schematic of the apodized waveguide with the last apodization function of
zpower
is shown in
Fig. 9(a)
. Moreover, the simulation results of applying the apodization on the waveguide with the length of 620 μm and the corrugation width of 5 nm are presented in
Fig. 9(b)
and
9(c)
, with a summarization reported in
Table 6
.
(a) Schematic of the apodized waveguide with the apodization function of zpower for m=1, n=1, p=1.5, q=0. Reflection spectra of the apodized waveguide for: (b) m=1, n=1, q=0, p=0.2, p=0.5 and p=1.5 and (c) Comparison between the unapodized and the apodized waveguides for m=1, n=1, p=3 and q=0.
Spectral features of the apodized waveguide with the apodization function ofzpower
Spectral features of the apodized waveguide with the apodization function of zpower
As shown in
Fig. 9(b)
, by increasing the apodization
p
parameter, the sidelobe oscillations gradually diminish until they are totally removed for
p
=1.5. By further increasing of this parameter, the spectrum remains flat, but the FWHM and the reflectivity are reduced.
A comparison between the unapodized and apodized waveguides for
m
=1,
n
=1,
p
=3 and
q
=0 is illustrated in
Fig. 9(c)
. As it is clear, utilizing the apodization has considerable impact on improving the filtering characteristics of the structure and a totallyflat reflection spectrum can be obtained. Absence of the sidelobes enhances the performance of the structure and makes it an attractive choice to be used in applications such as WDM communication systems, where high level sidelobes lead to cross talk between adjacent channels.
On the other hand, by decreasing the apodization
q
 and
n
parameters, the FWHM, the sidelobe level and the reflectivity would be decreased. Based on the results of
Table 6
, the minimum value of the FWHM is 1.05 nm for
m
=1,
n
=−1,
p
=0 and
q
=−1. However, in this case the reflectivity is small.
Figure 10
shows the spectra of the apodized waveguide for
m
=1,
n
=−1,
p
=0,
q
=−1, with different lengths. Also, the results of length increment are listed in
Table 7
.
Comparison between the reflection spectra of the unapodized and the apodized waveguides with the apodization function of zpower for m=1, n=−1, p=0, q=−1 and different lengths.
Spectral features of the apodized waveguide with the apodization function ofzpowerform=1,n=−1,p=0,q=1 and different lengths
Spectral features of the apodized waveguide with the apodization function of zpower for m=1, n=−1, p=0, q=1 and different lengths
As can be deduced from
Fig. 10
and
Table 7
, higher reflectivity and lower FWHM are obtained by increasing the grating length. By comparing the spectra of unapodized and the apodized waveguides, it is found that much narrower spectra with lower sidelobes can be achieved by using the apodization. Moreover, the FWHM of 0.74 nm for the grating length of 1300 μm is the minimum value among the apodized waveguides proposed in this paper.
IV. CONCLUSIONS
In this paper, four apodization functions, applied to SOI strip waveguide with sidewall corrugated gratings are proposed. The effects of apodization functions on the FWHM, the sidelobe level, and the reflectivity are studied. Compared to the reported results for unapodized waveguide by Wang
et al
.
[9]
, proposed apodization functions bring considerable enhancement in filtering behavior of the structure. To reduce the sidelobes, the best performed apodization functions are
Polynomial
and
zpower
, which result in the completely smooth reflection spectra, without any sidelobe oscillations. By applying the apodization of
Exponential 1
and
Exponential 2
on the 620μmlong waveguides, the lowvalue sidelobes of 15.7 dB with too small fluctuations and 20.18 dB are obtained, respectively. In addition to removing the oscillations and reducing the sidelobe level, four functions are also used to adjust the bandwidth. The apodization functions of
Exponential 1
,
Polynomial
and
zpower
, substantially decrease the FWHM while
Exponential 2
can be applied to increase the FWHM. By comparing the results, it is clear that the minimum value of the FWHM is 0.74 nm that is obtained for the 1300μmlong waveguide with the apodization function of
zpower
. On the other hand, using the apodization function of
Exponential 2
, leads to the maximum value of the FWHM of 8.48 nm, for the 620μmlong waveguide. In order to increase the reflectivity of the apodized waveguides with different functions, the simulations are also performed for various grating lengths.
Prabhathan P.
,
Murukeshan V.
,
Jing Z.
,
Ramana P. V.
2009
“Compact SOI nanowire refractive index sensor using phase shifted Bragg grating,”
Opt. Express
17
15330 
15341
Vivien L.
,
Pascal D.
,
Lardenois S.
,
MarrisMorini D.
,
Cassan E.
,
Grillot F.
,
Laval S.
,
Fédéli J.M.
,
El Melhaoui L
2006
“Light injection in SOI microwaveguides using highefficiency grating couplers,”
J. Lightwave Technol
24
3810 
3815
Ozoliņš O.
,
Ivanovs Ģ.
2009
“Realization of optimal FBG bandpass filters for high speed HDWDM,”
Lithuanian J. Electron. Elec. Eng.
4
41 
44
Riziotis C.
,
Zervas M. N.
2003
“Novel fullcyclecouplerbased optical adddrop multiplexer and performance characteristics at 40Gb/s WDM networks,”
J. Lightwave Technol
21
1828 
1837
Shi W.
,
Wang X.
,
Zhang W.
,
Chrostowski L.
,
Jaeger N. A. F.
2011
“Contradirectional couplers in silicononinsulator rib waveguides,”
Opt. Lett.
36
3999 
4001
Hastings J. T.
,
Lim M. H.
,
Goodberlet J. G.
,
Smith H. I.
2002
“Optical waveguides with apodized sidewall gratings via spatialphaselocked electronbeam lithography,”
J. Vac. Sci. Technol. B
20
2753 
2757
Shi W.
,
Wang X.
,
Lin C.
,
Yun H.
,
Liu Y.
,
BaehrJones T.
,
Hochberg M.
,
Jaeger N. A. F.
,
Chrostowski L.
2013
“Silicon photonic gratingassisted, contradirectional couplers,”
Opt. Express
21
3633 
3650
Murphy T. E.
,
Hastings J. T.
,
Smith H. I.
2001
“Fabrication and characterization of narrowband Braggreflection filters in silicononinsulator ridge waveguides,”
J. Lightwave Technol
19
1938 
1942
Wang X.
,
Shi W.
,
Vafaei R.
,
Jaeger N. A. F.
,
Chrostowski L.
2011
“Uniform and sampled Bragg gratings in SOI strip waveguides with sidewall corrugations,”
IEEE Photon. Technol. Lett.
23
290 
292
Wang X.
,
Shi W.
,
Yun H.
,
Grist S.
,
Jaeger N. A. F.
,
Chrostowski L.
2012
“Narrowband waveguide Bragg gratings on SOI wafers with CMOScompatible fabrication process,”
Opt. Express
20
15547 
15558
Wang X.
,
Shi W.
,
Vafaei R.
,
F Jaeger N. A.
,
Chrostowski L.
2010
“Silicononinsulator Bragg gratings fabricated by deep UV lithography,”
in Proc. ACP
501 
502
Tan D. T. H.
,
Ikeda K.
,
Zamek S.
,
Mizrahi A.
,
Nezhad M. P.
,
Krishnamoorthy A. V.
,
Raj K.
,
Cunningham J. E.
,
Zheng X.
,
Shubin I.
,
Luo Y.
,
Fainman Y.
2011
“Wide bandwidth, low loss 1 by 4 wavelength division multiplexer on silicon for optical interconnects,”
Opt. Express
19
2401 
2409
Taillaert D.
,
Bienstman P.
,
Baets R.
2004
“Compact efficient broadband grating coupler for silicononinsulator waveguides,”
Opt. Lett.
29
2749 
2751
Bulk M. P.
,
Knights A. P.
,
Jessop , P. E.
,
Waugh P.
,
Loiacono R.
,
Mashanovich G. Z.
,
Reed G. T.
,
Gwilliam R. M.
2008
“Optical filters utilizing ion implanted Bragg gratings in SOI waveguides,”
Adv. Opt. Technol.
1
1 
6
Honda S.
,
Wu Z.
,
Matsu J.
,
Utaka , K.
,
Edura T.
,
Tokuda M.
,
Tsutsui K.
,
Wada Y.
2007
“Largelytunable wideband Bragg gratings fabricated on SOI rib waveguides employed by deepRIE,”
Electron. Lett.
43
630 
631
Wiesmann D.
,
David C.
,
Germann R.
,
Emi D.
,
Bona G.
2000
“Apodized surfacecorrugated gratings with varying duty cycles,”
IEEE Photon. Technol. Lett.
12
639 
641
He J.
,
Sun X.
,
Zhang M.
2000
“A novel add/drop multiplexer architecture for DWDM network,”
J. Infrared and Millimeter Waves
21
57 
62
De Heyn P.
,
Verstuyft S.
,
Keyvaninia S.
,
Trita A.
,
Van Thourhout D.
2012
“Tunable 4channel ultradense WDM demultiplexer with IIIV photodiodes integrated on silicononinsulator,”
in Proc. ACP
1 
3
Weber J. P.
1993
“Spectral characteristics of coupledwaveguide Braggreflection tunable optical filter,”
Optoelectronics, IEE Proceedings J.
140
275 
284
Riziotis Ch.
,
Zervas M. N.
2001
“Design considerations in optical add/drop multiplexers based on gratingassisted null couplers,”
J. Lightwave Technol.
19
92 
104
Liu G. J.
,
Li Q.
,
Jin G. L.
,
Liang B. M.
2004
“Transfer matrix method analysis of apodized grating couplers,”
Opt. Commun.
235
319 
324