Abstract

We investigate coupled-mode theory in designing high index contrast silicon-on-insulator waveguide couplers and arrayed waveguides. We develop and demonstrate a method of solution to the inverse problem of reconstructing the coupling matrix from the modal profiles obtained, in this case, from finite-difference frequency-domain field calculations. We show that whereas supermode theory provides a good approximation of the mode profiles, next-to-nearest-neighbor coupling becomes significant at small separation distances between arrayed waveguides. These distances are quantified for three different silicon-on-insulator material platforms. We also point out the phenomenon of field skewing and deformation at small separations.

© 2009 Optical Society of America

1. Introduction

Silicon-on-insulator (SOI) waveguides and optical on-chip circuitry rely on the high refractive index contrast between core (silicon, n = 3.5) and cladding (silicon dioxide, n = 1.45) materials to guide light in very compact structures and with small bending radii [1, 2]. SOI photonics is one of the most active areas of ongoing research and large-scale integrated circuits are being designed and fabricated. In many of these proposed circuits, one of the most critical waveguide components is the directional coupler between two parallel waveguides, which is used in microring-based filters [3], Mach-Zehnder interferometers and modulators [4, 5, 6], arrayed waveguide structures [7, 8, 9, 10] etc.

Here, we investigate and quantify the limitations of coupled-mode theory (CMT) in designing high index contrast SOI couplers and arrayed waveguides. CMT is a simple and often reliable approach to the design of such structures, in which the coupling coefficient can be written down in terms of overlap integrals of the individual waveguide modes and the refractive index distribution n(x,y) in the cross-sectional plane [11, 12, 13]. More accurate corrections to CMT for slab waveguides have been investigated by Chiang [14], and Payne [15], among others (see references therein).

In this paper, we compare “exact” results of the coupling coefficients of directional couplers calculated using a fully-vectorial finite difference frequency-domain (FDFD) mode-solving computer program with predictions of CMT. We demonstrate how to solve the inverse problem of reconstructing the coupling matrix from the solutions of the FDFD program. The significance of these results lies, firstly, in the resulting simplification of the design process of high-index contrast SOI waveguide directional couplers and coupled-waveguide structures, providing graphs and rule-of-thumb estimates for the validity and invalidity of simple CMT instead of having to rely on time-consuming ab initio numerical simulations in every case. Secondly, our method of solution of the inverse problem can be easily applied to design couplers from a mode solver rather than a time-consuming propagation simulation.

As a test structure which highlights both the applicability and shortcomings of CMT, we will consider the multi-waveguide coupled-array structure [16] which consists of a number of directional couplers parallel to one another, as is often used in arrayed waveguide gratings and multi-element lasers and amplifiers, and which require an accurate estimate of the coupling coefficients to prevent imaging and phase errors [17]. This structure, also called the “multi-slot waveguide,” has been recently fabricated and demonstrated by us in SOI [18]. We show that, because of the narrow waveguide widths allowed by the high index contrast, multi-waveguide structures can reveal significant next-to-nearest-neighbor coupling and other deviations from the conventional picture of modal coupling, and we find the “critical” waveguide-to-waveguide separation distance at which such terms become significant.

2. Coupled-mode theory (CMT) of the modes of multi-slot waveguides

If many slot waveguides are arranged in a parallel array, as would be encountered in the cross-section of an arrayed waveguide grating or coupler, or coupled-waveguide laser, then their modes can often be adequately described by supermode theory [16], which is one of the fundamental predictions of coupled-mode theory (CMT), and hence can be a test of the applicability of CMT to SOI photonics. The first part of this section will briefly describe an ab initio numerical algorithm we have encoded in MATLAB to calculate the modes accurately. The next part of this section presents the analysis of the modes based on supermode theory. In the following sections, we will compare the predictions of CMT with the FDFD calculations.

2.1. Finite-difference frequency-domain (FDFD) algorithm

In the finite-difference frequency-domain (FDFD) algorithm, the dielectric profile of the waveguide’s cross section is discretized on a rectangular grid. As developed by C.L. Xu et al. [20], the vectorial wave equation,

2E+(n2k2β2)E=[·E1n2·(n2E)]

where E (x,y) is the transverse electric field vector, is written in matrix form as

(PxxPxyPyxPyy)(ExEy)=β2(ExEy)

and,

PxxEx=x[1n2(n2Ex)x]+2Exy2+n2k2Ex,PxyEy=x[1n2(n2Ey)y]+2Eyxy,
PyxEx=y[1n2(n2Ex)x]+2Exyx,PyyEy=y[1n2(n2Ex)y]+2Eyx2+n2k2Ey.

In this formulation the displacement vectors n 2 Ex and n 2 Ey are both continuous across any dielectric discontinuity, and with the graded-index approximation, the central difference equations can be applied directly without any special treatment at the boundaries. Fully discretized versions of these operators can be found in Ref. [21]. We have written a set of routines in MAT-LAB to calculate the modes and their effective indices for arbitrarily-shaped index profiles.

 

Fig. 1. (a) Refractive index profile in the transverse plane n(x,y) for an N arrayed-waveguide structure. (b) The same refractive index profile may be decomposed, mathematically, into the sum of parts, Δn2 i = n2 i -n2 s , each of which appears in integrals equation for the coupling coefficients. For the structures considered in this paper, h = 500 nm, w = 200 nm, and s varies over the range 50 nm to 1μm. For these waveguide widths and heights (similar to those in Ref. [18]), the polarization direction of the principal transverse component of the electric field is indicated for the (quasi) TE and TM modes.

Download Full Size | PPT Slide | PDF

2.2. Coupled mode theory and its predictions

To describe the modes of multislot waveguides, we begin with the wave equation [19],

2E+ω2c2n2(x,y)E=0,

and consider each polarization in turn. (TE and TM polarization are defined in terms of as the major component of the electric field, which for the structure in question, are polarized vertically and horizontally as shown in Fig. 1.)

2.2.1. TE Polarization

Consider an array of N single mode waveguides, whose refractive index profile is shown schematically in Fig. 1a. As long as the waveguides are not too close to each other (i.e., greater than a separation distance which we will investigate and quantify in a subsequent section), the transverse mode profile of the multislot waveguide structure can be approximated by an expansion of the individual high index waveguide modes

Ê=ŷE(x,y)eiβz=ŷ[Σl=1NAl𝓔1(x,y)]eiβz.

As shown by Fig. 1, the relative dielectric coefficient distribution of the entire N waveguide structure n 2 (x,y) can be written as a sum of individual waveguide contributions, so that

n2(x,y)=ns2(x,y)+Σl=1NΔnl2(x,y)

where n 2 s (x,y) corresponds to the cladding. Thus, n 2 s (x,y) + Δn 2 l (x,y) would yield the dielectric coefficient profile of the l th waveguide in the absence of the others. Substituting the above two equations into the wave equation, we have,

(2+ω2c2[ns2(x,y)+Σl=1NΔnl2(x,y)]β2)[Σl=1NAl𝓔1(x,y)]=0.

The modes of the individual waveguides satisfy their respective eigenvalue equations,

(2+ω2c2[ns2(x,y)+Δnl2(x,y)]βl2)𝓔1(x,y)=0

and therefore, using Eq. (7), Eq. (6) can be written as,

Σl=1NA1(Δl+ω2c2Σm=1mlNΔnm2(x,y))𝓔1(x,y)=0

where,

Δ1βl2β2.

N equations are formed by multiplying Eq. (8) by 𝓔j * (j = 1,2,…,N), and integrating each of these equations over x and y,

Σl=1NA1(Δl∫∫𝓔j*𝓔1dxdy+ω2c2Σm=1mlN∫∫𝓔j*Δnm2(x,y)𝓔1dxdy)=0,j=1,2,,N.

We define the modal overlap integrals as follows:

Ijl=∫∫𝓔j*𝓔1dxdy,
κjl=ω2c2Σm=1mlN∫∫𝓔j*Δnm2(x,y)𝓔1dxdy,

with the normalization

∫∫𝓔j*𝓔1dxdy=1.

Ijl is the overlap integral of the modes of two waveguides which are not orthogonal to each other (particularly in the case of small waveguide separation), and κ are the self-coupling and cross-coupling (exchange coupling) coefficients familiar from coupled-mode theory [19, p. 362].

Equation (9) can then be written in matrix form as an eigenvalue problem,

(β12+κ11Δ2I12+κ12ΔN1I1,N1+κ1,N1ΔNI1N+κ1NΔ1I21+κ21β22+κ22ΔN1I2,N1+κ2,N1ΔNI2N+κ2NΔ1I31+κ31Δ2I32+κ32ΔN1I3,N1+κ3,N1ΔNI3N+κ3NΔ2IN1+κN1Δ2IN2+κN2ΔN1I1,N1+κN,N1βN2+κNN)(A1A2A3AN)
=(β20000β20000β20000β2)(A1A2A3AN).

(The matrix on the left-hand side of the above equation will be referred to as M.)

 

Fig. 2. TE Polarization Ey: The modes of an N = 5 coupled waveguide array for λ = 1550 nm, calculated using coupled-mode theory (blue solid lines), and a finite-difference frequency-domain algorithm (black crosses). The coupled-mode theory calculations were done by using the effective index method, calculating the overlap integrals, solving Eq. (11), and reassembling the field. Waveguide height = 500 nm, width = 200 nm, separation = 200 nm, ncore = 3.47, and nclad =1.46. Under nearest neighbor coupling, the scaling relationship predicted by Eq. (13) adequately predicts the field amplitudes within each waveguide.

Download Full Size | PPT Slide | PDF

If we assume that only nearest neighbor coupling is significant, then the integrals in Eq. (10) are nonzero only when l = j -1, j, j +1. M takes the tridiagonal form,

M=(β12+κ11Δ2I12+κ1200Δ2I21+κ21β22+κ22000Δ2I32+κ32000ΔN1IN,N1+κN,N1βN2+κNN).

If the waveguides are identical and equally spaced, M can be further simplified by setting β 2 1 = β 2 2 … = β N 1β 2 0 and also, I l,l+1 = I l-1,lI,κ l,l+1 = κ l-1,lκ. However, even if the waveguides are identical and equally spaced, κ 11 and κNN are not equal to κ 22, κ 33,…,κ N-1N-1. In fact, Eq. (10b) shows that for those waveguides at the edges (l = 1 and I = N) there are approximately only half as many contributing terms as the other waveguides: there are no waveguides to the left of the l = 1 waveguide, and there are no waveguides to the right of the l = N waveguide, whereas all the other waveguides have contribution terms from both the left and right halves of their modal profiles.

We define κ selfκ 22, κ 33,…, κ N-1N-1, κ self,edgeκ 11 and κNN, and δκ self = κ self - κ self,edge.

 

Fig. 3. TM Polarization Ex: The modes of an N = 5 coupled waveguide array, calculated using coupled-mode theory (blue solid lines), and a finite-difference frequency-domain algorithm (black crosses). The coupled-mode theory calculations were done by using the effective index method, calculating the overlap integrals, solving Eq. (11), and reassembling the field. Waveguide height = 500 nm, width = 200 nm, separation = 1μm, ncore = 3.47, and nclad =1.46. Under nearest neighbor coupling, the scaling relationship predicted by Eq. (13) adequately predicts the field amplitudes within each waveguide.

Download Full Size | PPT Slide | PDF

To first order in the perturbation δκ self, the eigenvectors are

Al(m)=(2N+1)1/2sinlmπN+1δκself(2N+1)3/2
×Σm=1nmN1β2(m)β2(m)[sinN+1sinN+1+sinmNπN+1sinnNπN+1]sinlnπN+1

where m is the modal number and l indicates which high-index rib waveguide (or low-index slot) is being described. [The expression for the eigenvalues is written later, Eq. (22).]

For large N, the second term in the above expression, is smaller than the first by (N + 1)-1 and can be ignored, yielding a simpler expression. The progression of peak-amplitude values (in the high index regions) {A (m) l},l = 1,…,N matches with the numerical calculations shown in Fig. 2. However, we shall see that the agreement is good only at large separation distances between the individual waveguides.

2.2.2. TM Polarization

For the TM polarization (in which the electric field is normal to the waveguide/slot boundary), we again start with the wave equation only now defined in terms of the magnetic field, which we expand in terms of the individual waveguide modes,

H=ŷH(x,y)eiβz=ŷ[Σl=1NAlH1(x,y)]eiβz.

Under nearest neighbor coupling, the magnetic field will also obey the scaling relationship of Eq. (13). If it can be assumed that ∣∂Hz/y∣ ≪ ∣∂Hy/z∣, then Hy and Ex are related by

Ex=βε(x,y)ωHy,

so that, within the high index ribs, we expect that Eq. (13) describes the scaling relationship of the peak electric field amplitudes. Although the peak electric field amplitudes of the entire modal profile are found not in the high-index regions, but in the low-index regions (just inside the core-cladding boundary), they obviously satisfy the same scaling law, as can be seen in Fig. 3.

3. Numerically-assisted CMT: The “Inverse Problem”

CMT offers valuable physical insight into how waveguides couple—in particular, the structure of matrix M in Eq. (11) is revealing—but the quantitative predictions of CMT are in error in high-index-contrast SOI structures at short separation distances. To obtain a numerically-accurate picture of modal coupling, we propose a new extension of CMT, which we call “numerically assisted” CMT, to use the simulation results of the FDFD algorithm to back-calculate the elements of the coupling matrix M. We can thereby check if the assumption of nearest-neighbor coupling is valid at short separation distances, and identify various other interesting coupling phenomena (e.g., non-Hermiticity of M) which have not been pointed out earlier.

To develop NA-CMT, we use the following mathematical procedure, based on a matrix theorem previously developed for coupled-resonator structures [23].

  1. First, we solve for the supermodes using FDFD, which does not contain any of the limitations of nearest-neighbor CMT under investigation. The propagation constants of the supermodes are also obtained by this algorithm.
  2. Having obtained both the eigenvectors (peak amplitudes) and eigenvalues (propagation constants), we construct the (non-singular) matrix of eigenvectors Amn (whose columns are the linearly-independent supermodes), and the diagonal matrix of eigenvalues, Λ = diag{β 2 m}.
  3. Next, we reassemble M [see Eq. (11)] by using the matrix theorem cited in Ref. [23, Eq. (7), Lem. 1–2]: if the eigenvalues are distinct (which they are in this case), M can be reconstructed as follows: M = AΛA -1. The matrix is unique to within a similarity transformation, which does not affect the following step. An example is shown in Table 1. (Notice that κ 11 and κ 55 are approximately one-half of κ 22, κ 33, or κ 44, as discussed earlier.) The values of the reconstructed M matrix may be useful to design couplers in the strongly-coupled regime from the output of the FDFD mode-solver algorithm itself, without having to carry out time-consuming beam-propagation simulations.

4. Asymptotic accuracy of numerically-assisted CMT

In this paper, the eigenvalues (λk) and eigenvectors (A (k)) are obtained from a computer simulation. However, they may be obtained from measurements on fabricated structures, in order to test whether the intended coupling matrix was successfully obtained in practice. The experimental procedure to measure eigenvalues and eigenvectors could be similar to that used to image the modes of laser resonators.

We assume that the measurements result in some small, uncorrelated errors in the eigenvalues (Δλk) and eigenvectors (Δuk). The inversion algorithm presented in the previous section can also be used with measured data. In this section we study the accuracy of the nearest-neighbor coupling and next-to-nearest-neighbor coupling coefficients in terms of Δλk and Δuk.

First, we will carry out a simple theoretical estimation. We will assume that the coupling matrix is Hermitian. After some algebra, the error in any element of M can be written to first order as

ΔMij=Σk=1N[(uik+ujk)Δukλk+uikujkΔλk].

For simplicity, in this paper, we assume that Δuk is zero, i.e., the errors are only in the measured eigenvalues, since for identical arrayed waveguide structures, successive eigenvectors look quite different from each other and are easily distinguished [23]. We will also use the simpler form of the eigenvectors, retaining only the first term of Eq. (13), so that the errors in the reconstructed nearest-neighbor coupling and next-to-nearest-neighbor coupling coefficients are

ΔMjj+1=Σk=1N2N+1sinjkπN+1sin(j+1)N+1Δλk,
ΔMj1j+1=Σk=1N2N+1sin(j1)N+1sin(j+1)N+1Δλk.

Assuming that Δλk are uncorrelated identically-distributed random variables with mean Eλ] and variance Var[Δλ], the mean and variance of the nearest-neighbor coupling and next-to-nearest-neighbor coupling coefficients can be calculated. Both ΔM jj+1 and ΔM j-1j+1 are zero-mean, since, for example,

E[ΔλMjj+1]=E[Δλ]Σk=1N2N+1sinjkπN+1sin(j+1)N+1=0

because the summation vanishes as a consequence of the orthogonality of the eigenvectors (the sum is equal to u (j).u (j+1) = 0).

To calculate the variance, we see that

Var[ΔλMjj+1]=Var[Δλ](1N+1)2Σk=1N(cos2N+1cos2jkπN+1)2
=Var[Δλ](1N+1),

and the same result is obtained for Var[ΔM j-1j+1].

Summarizing the results,

Nearestneighbor:E[ΔλMjj+1]=0,Var[ΔλMjj+1]=Var[Δλ]/(N+1),
Nexttonearest:E[ΔλMj1j+1]=0,Var[ΔλMj1j+1]=Var[Δλ]/(N+1).

Numerical calculations, shown in Fig. 4 confirm Eq. (;20). (Numerical calculations show that the same relations are seen to hold in the case of non-identical waveguides, in which case the off-diagonal terms of the coupling matrix are not identical along the sub-diagonals, and also for slightly asymmetric matrices.)

 

Fig. 4. Error versus N: Exact eigenvalues of a tridiagonal symmetric matrix of size N were perturbed by values chosen from a uniform random distribution with variance chosen to be ten percent of the first eigenvalue. The variance and mean of the reconstructed nearest-neighbor coupling and next-to-nearest-neighbor coupling coefficients are plotted, calculated from a distribution of coupling matrices generated by 105 iterations, showing that Eq. (20) is a good predictor of the reconstruction accuracy.

Download Full Size | PPT Slide | PDF

These results show that the error in reconstructing the coupling coefficients decreases, rather than increases, as the number of inaccurately-measured eigenvalues increases. This results from (spectral) averaging: each reconstructed coupling coefficient averages over the entire spectrum of eigenvalues, and therefore, benefits from the law of averages. In contrast, directly measuring a coupling coefficient e.g., by a local near-field probe of the field in the coupling region, does not benefit from any ensemble averaging.

5. Discussion

5.1. Next-to-nearest-neighbor coupling

As Table 1 shows (calculated at one specific value of the waveguide separation distance), M contains useful information about non nearest-neighbor coupling. We can read off whichever coupling coefficients are needed: in particular, we calculate the ratio ∣κ 13/κ 12∣, i.e., the ratio of next-to-nearest-neighbor coupling coefficient to the nearest-neighbor coupling coefficient.

Tables Icon

Table 1. An example of a reconstructed coupling (M) matrix from FDFD calculations of eigenmodes and eigenvalues. Si/SiO2, TE polarization, separation s = 350 nm, β 0 = 2.26128 (2π/λ). Although the nearest-neighbor coupling coefficients dominate, the self-coupling and off-tridiagonal coupling terms are non-zero.

First we estimate the expected dependency of this ratio of coupling coefficients to the edge-to-edge separation, s. Using Kuznetsov’s solution for the coupling coefficients of two slab waveguides [22], we observe that κ in both the TE and TM cases varies with s as κ ~ e -ps where p is the field decay length in the cladding. Therefore, the ratio κ 13/κ 12 for both polarizations has the following expression (to leading order),

Δ3I13+κ13Δ2I12+κ13κ13κ12=ep(2s+w)eps=ep(s+w)

i.e., the ratio of next-to-nearest-neighbor coupling coefficient to the nearest-neighbor coupling coefficient should fall off exponentially with increasing separation.

Figure 5 shows the calculations of this ratio using the above algorithm. The exponential fit describes the TE polarization much better than it does the TM polarization, indicating that some of the central assumptions of CMT are starting to fail for the TM polarization at short distances. The next section will describe another symptom of the failure of CMT, obtained by looking at the eigenvalues, i.e., the propagation constants, of the supermodes.

5.2. Eigenvalue fanout: effective index of the supermodes versus separation distance

Another way to evaluate the predictions of CMT is the theory behind the eigenvalues of Eq. (11), which predicts that the effective index of the m-th supermode is given by the equation

β2(m)=β02+κself+2(κ+Δ0I)cosN+12δκselfN+1(sin2N+1+sin2NmπN+1),

where β 0 is the propagation constant of a single waveguide in isolation. Note that for N = 5, the m = 3 supermode has the special property that the right-hand-side of the above expression , i.e., the index of that supermode does not change with the coupling coefficient κ. Hence, n (3) eff is only weakly dependent on the separation distance (through the self-coupling coefficients, κ 11,κ 22,…,κ 55).

To verify this prediction, Fig. 6 shows the effective index calculated by FDFD for each of the five supermodes at various separation distances in three different silicon-based material systems. These values of the effective index take into account coupling-induced frequency shifts (CIFS, [25]) because M itself results from a numerical calculation of the supermodes (and their eigenfrequencies), rather than individual waveguide modes and the propagation constants of isolated waveguides.

In the limit of large separation, the effective indices of all the supermodes tends to that of the single waveguide. As the separation distance is decreased, the coupling coefficients increase, and the effective indexes of the different modes separate [24]. The first three modes remain guided even as s shrinks to zero, since their effective indices are higher than that of the single waveguide. From Fig. 6, one can read off the waveguide separation distance at which conventional CMT is expected to fail, and more accurate design tools, such as FDFD calculations, should be used to accurately predict the coupling coefficients.

 

Fig. 5. Ratio of coupling coefficients for different separation distances extracted from Eq. (11), which was reconstructed using an algorithm described in the text. (a) TE Polarization An exponential fit expected from a simple nearest-neighbor-coupling theory holds throughout this regime. (b) TM Polarization At a separations less than 450 nm, the ratio deviates significantly from the predicted behavior. (c) TE Polarization The ratio of cross coupling coefficients show that the reconstructed coupling matrix M becomes asymmetric as the waveguide separation is reduced. (d) TM Polarization The asymmetry of the coupling matrix begins at a larger separation.

Download Full Size | PPT Slide | PDF

An interesting observation obtains from the m = 3 supermode: at a certain (small) waveguide separation, n (3) eff is no longer independent of s and begins to deviate substantially from a straight line, contrary to the prediction of Eq. (22). This deviation is much more pronounced in the case of the TM polarization.

 

Fig. 6. Left column: TE Polarization, and Right column: TM polarization. Effective index of the five supermodes for different separation distances with ncore = 3.47, and (a,b) nclad =1.46, (c,d) nclad = 1 (e,f) nclad = 2.05. For each case as the separation between the waveguides increases, the effective indexes of the modes converge to that of the single waveguide. These values are (a) neff = 2.36 and (b) neff = 1.66 for oxide cladding, (c) neff = 2.24 and (d) neff = 1.07 for air cladding, (e) neff = 2.56 and (f) neff = 2.26 for nitride cladding. The shaded regions indicate > 5% deviation of neff for the m=3 supermode from its theoretical value, which as discussed in the text, is predicted by CMT to be independent of the separation distance.

Download Full Size | PPT Slide | PDF

5.3. Field skewing and reshaping

At short separation distances, the reconstructed matrix M can become non-symmetric (non-Hermitian), although the eigenvalues remain stricly real as long as the mode is above cut-off. This can be seen in fact in the matrix written in Table 1 and Fig. 5(c,d): κ 12κ 21 and κ 13κ 31, etc.

 

Fig. 7. TM Polarization Ex: The field profile of the fifth eigenmode in the first waveguide. When the separation is decreased below 450 nm, the peak of the field in the high-index rib indicated by the dotted red line in (a) is no longer centered, and the mode shape is considerably altered, thereby changing both κ and n eff. Consequently, CMT can no longer accurately predict the mode coupling.

Download Full Size | PPT Slide | PDF

The reason for this asymmetry is that the fields within the individual waveguides are no longer centered between the dielectric boundaries. As shown in Fig. 7. the modal profile starts to deviate in the location of its maxima and minima. For example, the peaks of the field in the outermost ribs are skewed and no longer centered in the middle of the dielectric boundaries, and can even reach the boundaries of the high-index and low-index regions. It is no longer accurate to read off the peak amplitudes of the supermode in order to write the eigenvectors A (m) in Eq. (13)—doing so would result in asymmetric M matrices.

At a short separation distance of 80 nm, Fig. 8 shows the coupled-mode theory used to reconstruct the m = 1 and m = 5 supermode (plotted with continuous lines) and the supermode calculation of FDFD (with crosses). Note that in both cases the field is asymmetrically centered within the dielectric boundaries of the outer waveguides. Recall that CMT is based on writing the field as a summation of the scaled individual waveguide modes, Fig. 8(a,d), each of which is centered within its own core-cladding boundaries. At short separation distances, when, for example, there is a significant contribution of the (asymmetric) tail from the field in the second waveguide to the (symmetric) mode of the first waveguide, CMT itself predicts a lateral shift of the peak (of the sum) away from the exact center of the waveguide. The scaling relationships from Eq. (13) will enhance this effect for an multi-waveguide arrayed structure compared to a (twin-waveguide) directional coupler.

 

Fig. 8. TE Polarization Ey: Using the exact solution from a FDFD simulation of a single waveguide, the horizontal cross section is extracted and five copies are shifted from one another so that their separation corresponds to a waveguide separation of 80 nm. (a) These individual waveguide modes are scaled in accordance with Eq. (13) for the fundamental mode (m=1). (b) The sumation of the individual waveguide modes; superimposed is the FDFD solution of the entire five waveguide structure. (c) Zoomed in to just the first waveguide. CMT and FDFD show a shift of the mode towards the center of the waveguide structure. (d-e) The fifth mode, both CMT and FDFD show a shift towards the edge of the waveguide structure however FDFD shows a shift of greater magnitude.

Download Full Size | PPT Slide | PDF

For the fundamental mode, Fig. 8(a-c), the summation of the fields associated with the first (blue) and second (green) waveguides results in the peak shifting towards the center of the five waveguide structure, which qualitatively agrees with the FDFD simulation. But the FDFD result for the fifth mode shows a shift of greater magnitude, now towards the outer edge of the waveguide structure, indicating that CMT no longer accurately predicts the modal profile of the supermode. TM polarized modes start to shift at a larger separation, due to the field discontinuities at the waveguide boundaries and electric field enhancement in the cladding regions.

Note also, as shown in Fig. 8(e), that FDFD predicts a different exponential decay constant of the field wings, compared to CMT. This is a fundamental failure of CMT in the sense that the eigenmode of the composite structure can no longer be written as the sum of modes of individual waveguides (in isolation from each other). Until a satisfactory theory of mode evolution in this strong-coupling regime can be obtained, we recommend that designers rely on direct numerical simulations to obtain field profiles.

5.4. Polarization hybridization

Alongside this phenomenon, we also observe that the polarization becomes strongly hybridized as the separation distance is reduced. Figure 9 shows the major and minor E field components at large and small separation distances, along with the cross section of the Poynting vector, which indicates power flow. Notice that at small separation distances, the polarization component that was previously negligible has become, in fact, the dominant one. Furthermore, the power is actually carried above and below the waveguide structure at the outer edges, rather than within the inner slots, contrary to the original intention of slot waveguides.

These calculations suggest that at small separation distances, CMT of high-index contrast waveguides should consider both polarization components of the electric field. Mathematically, the dimensionality of the basis set (of eigenvectors) depends on the separation distance. Practically, one should pay careful attention to where the power flow is concentrated in waveguide couplers at short separation distances, since both polarization components of the electric field (with different coupling lengths due to different modal effective indices) play a role in the exchange of power between the constituent waveguides.

 

Fig. 9. TM Polarization: finite element simulations (COMSOL) of the transverse field amplitudes, Ex and Ey, and Power flow, Pz for a separation of 1 μm (first Row) and 100 nm (second row) for the fifth mode. For small separations the polarization becomes strongly hybridized and the power is carried at the outer edges.

Download Full Size | PPT Slide | PDF

6. Conclusion

In this paper, we have studied the validity of coupled-mode theory CMT for high-index contrast (e.g., silicon-based) optical waveguiding structures, in particular directional couplers and multi-slot waveguides. We have used a fully vectorial finite-difference frequency-domain (FDFD) algorithm to obtain modal profiles and effective indices of the supermodes in a non-perturbative way. When the modal profiles can be “discretized” to read off peak amplitudes within each of the waveguide cores, a theorem from matrix algebra can be employed to solve the “inverse problem”—to reconstruct accurately the matrix of coupling coefficients M (a procedure we have called numerically-assisted coupled-mode theory, NA-CMT). The NA-CMT framework can be used to find out when the nearest-neighbor-coupling approximation breaks down.

Aside from the inversion procedure, the results of FDFD calculations also directly address this problem. We have presented the eigenvalue fanout graphs for three different silicon-based material platforms and indicated the separation distances at which the CMT predictions break down. We have also shown that there are significant deviations in the modal shape as the waveguide separation distance decreases and the state of polarization of the field may change between quasi-TE or quasi-TM to a strongly-hybridized polarization.

We suggest that the terminology “strong coupling” in the context of waveguide directional couplers be used to describe this coupling regime, for which a simple theory is not yet available, rather than simply large values of the length-integrated coupling coefficient, which can approach unity even with small coupling coefficients per unit length for long couplers, as this latter regime has already been well studied.

Acknowledgment

The authors are grateful to the National Science Foundation for support (ECCS-0642603 and ECCS-0723055) and thank Jung S. Park and Mark A. Schneider for useful discussions.

References and links

1. F. Xia, L. Sekaric, and Y. Vlasov, “Ultracompact optical buffers on a silicon chip,” Nat. Photonics 1, 65–71 (2007). [CrossRef]  

2. Q. Xu, D. Fattal, and R. G. Beausoleil, “Silicon microring resonators with 1.5-μm radius,” Opt. Express 16, 4309–4315 (2008). [CrossRef]   [PubMed]  

3. F. Xia, M. Rooks, L. Sekaric, and Y. Vlasov, “Ultra-compact high order ring resonator filters using submicron silicon photonic wires for on-chip optical interconnects,” Opt. Express 15, 11934–11941 (2007). [CrossRef]   [PubMed]  

4. A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004). [CrossRef]   [PubMed]  

5. Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005). [CrossRef]   [PubMed]  

6. W. M. J. Green, M. J. Rooks, L. Sekaric, and Y. A. Vlasov, “Optical modulation using anti-crossing between paired amplitude and phase resonators,” Opt. Express 15, 17264–17272 (2007). [CrossRef]   [PubMed]  

7. X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

8. P. Cheben, J. H. Schmid, A. Delage, A. Densmore, S. Jannz, B. Lamontagne, J. Lapointe, E. Post, P. Waldron, and D. C. Xu, “A high-resolution silicon-on-insulator arrayed waveguide grating microspectrometer with sub-micrometer aperture waveguides,” Opt. Express 15, 2299–2306 (2007). [CrossRef]   [PubMed]  

9. K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).

10. P. Dumon, W. Bogaerts, D. V. Thourhout, D. Taillaert, and R. Baets, “Compact wavelength router based on a Silicon-on-insulator arrayed waveguide grating pigtailed to a fiber array,” Opt. Express 14, 664–669 (2006). [CrossRef]   [PubMed]  

11. H. Kogelnik and C. V. Shank, “Coupled-mode theory of distributed feedback lasers,” Appl. Phys. 43, 2327–2335 (1972). [CrossRef]  

12. A. Hardy and W. Streifer, “Coupled-mode theory of parallel waveguides,” J. Lightwave Technol. LT-3, 1135–1146 (1985). [CrossRef]  

13. W. P. Huang, “Coupled-mode theory for optical waveguides: An overview,” J. Opt. Soc. Am. A 11, 963–983 (1994). [CrossRef]  

14. K. S. Chiang, “Coupled-zigzag-wave theory for guided waves in slab waveguide arrays,” J. Lightwave Technol. 10, 1380–1387 (1992). [CrossRef]  

15. F. P. Payne, “An analytical model for the coupling between the array waveguides in AWGs and star couplers,” Opt. Quantum Electron. 38, 237–248 (2006). [CrossRef]  

16. E. Kapon, J. Katz, and A. Yariv, “Supermode analysis of phase-locked arrays of semiconductor lasers,” Opt. Lett. 10, 125–127 (1984). [CrossRef]  

17. A. Klekamp and R. Munzner, “Calculation of imaging errors of AWG,” J. Lightwave Technol. 21, 1978–1986 (2003). [CrossRef]  

18. S. H. Yang, M. L. Cooper, P. R. Bandaru, and S. Mookherjea, “Giant birefringence in multi-slotted silicon nanophotonic waveguides,” Opt. Express 16, 8306–8316 (2008). [CrossRef]   [PubMed]  

19. P. Yeh, Optical Waves in Layered Media (John Wiley & Sons, New York, 2005).

20. C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994). [CrossRef]  

21. W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993). [CrossRef]  

22. M. Kuznetsov, “Expressions for the coupling coefficient of a rectangular waveguide directional coupler,” Opt. Lett. 8, 499–501 (1983). [CrossRef]   [PubMed]  

23. S. Mookherjea, “Spectral characteristics of coupled resonators,” J. Opt. Soc. Am. B 23, 1137–1145 (2006). [CrossRef]  

24. G. Lenz and J. Salzman, “Eigenmodes of multiwaveguide structures,” J. Lightwave Technol. 8, 1803–1809 (1990). [CrossRef]  

25. M. Popovic, C. Manolatou, and M. Watts, “Coupling-induced resonance frequency shifts in coupled dielectric multi-cavity filters,” Opt. Express 14, 1208–1222 (2006). [CrossRef]   [PubMed]  

26. E. Marcatili, “Improved coupled-mode equations for dielectric guides,” IEEE J. Quantum Electron. QE-22, 988–993 (1986). [CrossRef]  

27. H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987). [CrossRef]  

References

  • View by:
  • |
  • |
  • |

  1. F. Xia, L. Sekaric, and Y. Vlasov, “Ultracompact optical buffers on a silicon chip,” Nat. Photonics 1, 65–71 (2007).
    [Crossref]
  2. Q. Xu, D. Fattal, and R. G. Beausoleil, “Silicon microring resonators with 1.5-μm radius,” Opt. Express 16, 4309–4315 (2008).
    [Crossref] [PubMed]
  3. F. Xia, M. Rooks, L. Sekaric, and Y. Vlasov, “Ultra-compact high order ring resonator filters using submicron silicon photonic wires for on-chip optical interconnects,” Opt. Express 15, 11934–11941 (2007).
    [Crossref] [PubMed]
  4. A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
    [Crossref] [PubMed]
  5. Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
    [Crossref] [PubMed]
  6. W. M. J. Green, M. J. Rooks, L. Sekaric, and Y. A. Vlasov, “Optical modulation using anti-crossing between paired amplitude and phase resonators,” Opt. Express 15, 17264–17272 (2007).
    [Crossref] [PubMed]
  7. X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.
  8. P. Cheben, J. H. Schmid, A. Delage, A. Densmore, S. Jannz, B. Lamontagne, J. Lapointe, E. Post, P. Waldron, and D. C. Xu, “A high-resolution silicon-on-insulator arrayed waveguide grating microspectrometer with sub-micrometer aperture waveguides,” Opt. Express 15, 2299–2306 (2007).
    [Crossref] [PubMed]
  9. K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).
  10. P. Dumon, W. Bogaerts, D. V. Thourhout, D. Taillaert, and R. Baets, “Compact wavelength router based on a Silicon-on-insulator arrayed waveguide grating pigtailed to a fiber array,” Opt. Express 14, 664–669 (2006).
    [Crossref] [PubMed]
  11. H. Kogelnik and C. V. Shank, “Coupled-mode theory of distributed feedback lasers,” Appl. Phys. 43, 2327–2335 (1972).
    [Crossref]
  12. A. Hardy and W. Streifer, “Coupled-mode theory of parallel waveguides,” J. Lightwave Technol. LT-3, 1135–1146 (1985).
    [Crossref]
  13. W. P. Huang, “Coupled-mode theory for optical waveguides: An overview,” J. Opt. Soc. Am. A 11, 963–983 (1994).
    [Crossref]
  14. K. S. Chiang, “Coupled-zigzag-wave theory for guided waves in slab waveguide arrays,” J. Lightwave Technol. 10, 1380–1387 (1992).
    [Crossref]
  15. F. P. Payne, “An analytical model for the coupling between the array waveguides in AWGs and star couplers,” Opt. Quantum Electron. 38, 237–248 (2006).
    [Crossref]
  16. E. Kapon, J. Katz, and A. Yariv, “Supermode analysis of phase-locked arrays of semiconductor lasers,” Opt. Lett. 10, 125–127 (1984).
    [Crossref]
  17. A. Klekamp and R. Munzner, “Calculation of imaging errors of AWG,” J. Lightwave Technol. 21, 1978–1986 (2003).
    [Crossref]
  18. S. H. Yang, M. L. Cooper, P. R. Bandaru, and S. Mookherjea, “Giant birefringence in multi-slotted silicon nanophotonic waveguides,” Opt. Express 16, 8306–8316 (2008).
    [Crossref] [PubMed]
  19. P. Yeh, Optical Waves in Layered Media (John Wiley & Sons, New York, 2005).
  20. C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
    [Crossref]
  21. W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993).
    [Crossref]
  22. M. Kuznetsov, “Expressions for the coupling coefficient of a rectangular waveguide directional coupler,” Opt. Lett. 8, 499–501 (1983).
    [Crossref] [PubMed]
  23. S. Mookherjea, “Spectral characteristics of coupled resonators,” J. Opt. Soc. Am. B 23, 1137–1145 (2006).
    [Crossref]
  24. G. Lenz and J. Salzman, “Eigenmodes of multiwaveguide structures,” J. Lightwave Technol. 8, 1803–1809 (1990).
    [Crossref]
  25. M. Popovic, C. Manolatou, and M. Watts, “Coupling-induced resonance frequency shifts in coupled dielectric multi-cavity filters,” Opt. Express 14, 1208–1222 (2006).
    [Crossref] [PubMed]
  26. E. Marcatili, “Improved coupled-mode equations for dielectric guides,” IEEE J. Quantum Electron. QE-22, 988–993 (1986).
    [Crossref]
  27. H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
    [Crossref]

2008 (2)

2007 (5)

2006 (4)

2005 (1)

Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
[Crossref] [PubMed]

2004 (1)

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

2003 (1)

1994 (2)

W. P. Huang, “Coupled-mode theory for optical waveguides: An overview,” J. Opt. Soc. Am. A 11, 963–983 (1994).
[Crossref]

C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
[Crossref]

1993 (1)

W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993).
[Crossref]

1992 (1)

K. S. Chiang, “Coupled-zigzag-wave theory for guided waves in slab waveguide arrays,” J. Lightwave Technol. 10, 1380–1387 (1992).
[Crossref]

1990 (1)

G. Lenz and J. Salzman, “Eigenmodes of multiwaveguide structures,” J. Lightwave Technol. 8, 1803–1809 (1990).
[Crossref]

1987 (1)

H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
[Crossref]

1986 (1)

E. Marcatili, “Improved coupled-mode equations for dielectric guides,” IEEE J. Quantum Electron. QE-22, 988–993 (1986).
[Crossref]

1985 (1)

A. Hardy and W. Streifer, “Coupled-mode theory of parallel waveguides,” J. Lightwave Technol. LT-3, 1135–1146 (1985).
[Crossref]

1984 (1)

1983 (1)

1972 (1)

H. Kogelnik and C. V. Shank, “Coupled-mode theory of distributed feedback lasers,” Appl. Phys. 43, 2327–2335 (1972).
[Crossref]

Baba, T.

K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).

Baets, R.

Bandaru, P. R.

Beausoleil, R. G.

Bogaerts, W.

Chaudhuri, S. K.

C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
[Crossref]

Cheben, P.

Chen, X.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Chiang, K. S.

K. S. Chiang, “Coupled-zigzag-wave theory for guided waves in slab waveguide arrays,” J. Lightwave Technol. 10, 1380–1387 (1992).
[Crossref]

Cohen, O.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Cooper, M. L.

Dadap, J. I.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Delage, A.

Densmore, A.

Dumon, P.

Fattal, D.

Green, W. M.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Green, W. M. J.

Hardy, A.

A. Hardy and W. Streifer, “Coupled-mode theory of parallel waveguides,” J. Lightwave Technol. LT-3, 1135–1146 (1985).
[Crossref]

Haus, H. A.

H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
[Crossref]

Hsieh, I.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Huang, W. P.

W. P. Huang, “Coupled-mode theory for optical waveguides: An overview,” J. Opt. Soc. Am. A 11, 963–983 (1994).
[Crossref]

C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
[Crossref]

W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993).
[Crossref]

H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
[Crossref]

Jannz, S.

Jones, R.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Kapon, E.

Katz, J.

Kawakami, S.

H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
[Crossref]

Klekamp, A.

Kogelnik, H.

H. Kogelnik and C. V. Shank, “Coupled-mode theory of distributed feedback lasers,” Appl. Phys. 43, 2327–2335 (1972).
[Crossref]

Kuznetsov, M.

Lamontagne, B.

Lapointe, J.

Lenz, G.

G. Lenz and J. Salzman, “Eigenmodes of multiwaveguide structures,” J. Lightwave Technol. 8, 1803–1809 (1990).
[Crossref]

Liao, L.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Lipson, M.

Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
[Crossref] [PubMed]

Liu, A. S.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Liu, X.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Manolatou, C.

Marcatili, E.

E. Marcatili, “Improved coupled-mode equations for dielectric guides,” IEEE J. Quantum Electron. QE-22, 988–993 (1986).
[Crossref]

Mookherjea, S.

Motegi, A.

K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).

Munzner, R.

Nicolaescu, R.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Ohno, F.

K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).

Osgood, R. M.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Paniccia, M.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Panoiu, N. C.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Payne, F. P.

F. P. Payne, “An analytical model for the coupling between the array waveguides in AWGs and star couplers,” Opt. Quantum Electron. 38, 237–248 (2006).
[Crossref]

Popovic, M.

Post, E.

Pradhan, S.

Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
[Crossref] [PubMed]

Rooks, M.

Rooks, M. J.

Rubin, D.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Salzman, J.

G. Lenz and J. Salzman, “Eigenmodes of multiwaveguide structures,” J. Lightwave Technol. 8, 1803–1809 (1990).
[Crossref]

Samara-Rubio, D.

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Sasaki, K.

K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).

Schmid, J. H.

Schmidt, B.

Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
[Crossref] [PubMed]

Sekaric, L.

Shank, C. V.

H. Kogelnik and C. V. Shank, “Coupled-mode theory of distributed feedback lasers,” Appl. Phys. 43, 2327–2335 (1972).
[Crossref]

Stern, M.S.

C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
[Crossref]

Streifer, W.

A. Hardy and W. Streifer, “Coupled-mode theory of parallel waveguides,” J. Lightwave Technol. LT-3, 1135–1146 (1985).
[Crossref]

Taillaert, D.

Takekoshi, M.

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Thourhout, D. V.

Vlasov, Y.

Vlasov, Y. A.

W. M. J. Green, M. J. Rooks, L. Sekaric, and Y. A. Vlasov, “Optical modulation using anti-crossing between paired amplitude and phase resonators,” Opt. Express 15, 17264–17272 (2007).
[Crossref] [PubMed]

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Waldron, P.

Watts, M.

Whitaker, N. A.

H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
[Crossref]

Xia, F.

F. Xia, L. Sekaric, and Y. Vlasov, “Ultracompact optical buffers on a silicon chip,” Nat. Photonics 1, 65–71 (2007).
[Crossref]

F. Xia, M. Rooks, L. Sekaric, and Y. Vlasov, “Ultra-compact high order ring resonator filters using submicron silicon photonic wires for on-chip optical interconnects,” Opt. Express 15, 11934–11941 (2007).
[Crossref] [PubMed]

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Xu, C. L.

C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
[Crossref]

W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993).
[Crossref]

Xu, D. C.

Xu, Q.

Q. Xu, D. Fattal, and R. G. Beausoleil, “Silicon microring resonators with 1.5-μm radius,” Opt. Express 16, 4309–4315 (2008).
[Crossref] [PubMed]

Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
[Crossref] [PubMed]

Yang, S. H.

Yariv, A.

Yeh, P.

P. Yeh, Optical Waves in Layered Media (John Wiley & Sons, New York, 2005).

Appl. Phys. (1)

H. Kogelnik and C. V. Shank, “Coupled-mode theory of distributed feedback lasers,” Appl. Phys. 43, 2327–2335 (1972).
[Crossref]

Elec. Lett. (1)

K. Sasaki, F. Ohno, A. Motegi, and T. Baba,“Arrayed waveguide grating of 70×60 μm2 size based on Si photonic wire waveguides,” Elec. Lett. 41, (2007).

IEE Proc.-Optoelectron. (1)

C. L. Xu, W. P. Huang, M.S. Stern, and S. K. Chaudhuri, “Full-vectorial mode calculations by finite difference method,” IEE Proc.-Optoelectron. 141, 281–286 (1994).
[Crossref]

IEEE J. Quantum Electron. (2)

W. P. Huang and C. L. Xu, “Simulation of three-dimensional optical waveguides by a full-vector beam propagation method,” IEEE J. Quantum Electron. 29, 2639–2649 (1993).
[Crossref]

E. Marcatili, “Improved coupled-mode equations for dielectric guides,” IEEE J. Quantum Electron. QE-22, 988–993 (1986).
[Crossref]

J. Lightwave Technol. (5)

H. A. Haus, W. P. Huang, S. Kawakami, and N. A. Whitaker, “Coupled-mode theory of optical waveguides,” J. Lightwave Technol. LT-5 , 16–23 (1987).
[Crossref]

A. Hardy and W. Streifer, “Coupled-mode theory of parallel waveguides,” J. Lightwave Technol. LT-3, 1135–1146 (1985).
[Crossref]

G. Lenz and J. Salzman, “Eigenmodes of multiwaveguide structures,” J. Lightwave Technol. 8, 1803–1809 (1990).
[Crossref]

K. S. Chiang, “Coupled-zigzag-wave theory for guided waves in slab waveguide arrays,” J. Lightwave Technol. 10, 1380–1387 (1992).
[Crossref]

A. Klekamp and R. Munzner, “Calculation of imaging errors of AWG,” J. Lightwave Technol. 21, 1978–1986 (2003).
[Crossref]

J. Opt. Soc. Am. A (1)

J. Opt. Soc. Am. B (1)

Nat. Photonics (1)

F. Xia, L. Sekaric, and Y. Vlasov, “Ultracompact optical buffers on a silicon chip,” Nat. Photonics 1, 65–71 (2007).
[Crossref]

Nature (2)

A. S. Liu, R. Jones, L. Liao, D. Samara-Rubio, D. Rubin, O. Cohen, R. Nicolaescu, and M. Paniccia, “A highspeed silicon optical modulator based on a metal oxide-semiconductor capacitor,” Nature 427, 615–618 (2004).
[Crossref] [PubMed]

Q. Xu, B. Schmidt, S. Pradhan, and M. Lipson, “Micrometer-scale silicon electro-optic modulator,” Nature 435, 325–327 (2005).
[Crossref] [PubMed]

Opt. Express (7)

P. Cheben, J. H. Schmid, A. Delage, A. Densmore, S. Jannz, B. Lamontagne, J. Lapointe, E. Post, P. Waldron, and D. C. Xu, “A high-resolution silicon-on-insulator arrayed waveguide grating microspectrometer with sub-micrometer aperture waveguides,” Opt. Express 15, 2299–2306 (2007).
[Crossref] [PubMed]

F. Xia, M. Rooks, L. Sekaric, and Y. Vlasov, “Ultra-compact high order ring resonator filters using submicron silicon photonic wires for on-chip optical interconnects,” Opt. Express 15, 11934–11941 (2007).
[Crossref] [PubMed]

W. M. J. Green, M. J. Rooks, L. Sekaric, and Y. A. Vlasov, “Optical modulation using anti-crossing between paired amplitude and phase resonators,” Opt. Express 15, 17264–17272 (2007).
[Crossref] [PubMed]

Q. Xu, D. Fattal, and R. G. Beausoleil, “Silicon microring resonators with 1.5-μm radius,” Opt. Express 16, 4309–4315 (2008).
[Crossref] [PubMed]

S. H. Yang, M. L. Cooper, P. R. Bandaru, and S. Mookherjea, “Giant birefringence in multi-slotted silicon nanophotonic waveguides,” Opt. Express 16, 8306–8316 (2008).
[Crossref] [PubMed]

P. Dumon, W. Bogaerts, D. V. Thourhout, D. Taillaert, and R. Baets, “Compact wavelength router based on a Silicon-on-insulator arrayed waveguide grating pigtailed to a fiber array,” Opt. Express 14, 664–669 (2006).
[Crossref] [PubMed]

M. Popovic, C. Manolatou, and M. Watts, “Coupling-induced resonance frequency shifts in coupled dielectric multi-cavity filters,” Opt. Express 14, 1208–1222 (2006).
[Crossref] [PubMed]

Opt. Lett. (2)

Opt. Quantum Electron. (1)

F. P. Payne, “An analytical model for the coupling between the array waveguides in AWGs and star couplers,” Opt. Quantum Electron. 38, 237–248 (2006).
[Crossref]

Other (2)

P. Yeh, Optical Waves in Layered Media (John Wiley & Sons, New York, 2005).

X. Liu, I. Hsieh, X. Chen, M. Takekoshi, J. I. Dadap, N. C. Panoiu, R. M. Osgood, W. M. Green, F. Xia, and Y. A. Vlasov, “Design and fabrication of an ultra-compact silicon on insulator demultiplexer based on arrayed waveguide gratings,” in Proceedings of the Conference on Lasers and Electro-Optics (CLEO, 2008), paper CTuNN1.

Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (9)


                  Fig. 1.
Fig. 1.

(a) Refractive index profile in the transverse plane n(x,y) for an N arrayed-waveguide structure. (b) The same refractive index profile may be decomposed, mathematically, into the sum of parts, Δn2 i = n2 i -n2 s , each of which appears in integrals equation for the coupling coefficients. For the structures considered in this paper, h = 500 nm, w = 200 nm, and s varies over the range 50 nm to 1μm. For these waveguide widths and heights (similar to those in Ref. [18]), the polarization direction of the principal transverse component of the electric field is indicated for the (quasi) TE and TM modes.

Fig. 2.
Fig. 2.

TE Polarization Ey : The modes of an N = 5 coupled waveguide array for λ = 1550 nm, calculated using coupled-mode theory (blue solid lines), and a finite-difference frequency-domain algorithm (black crosses). The coupled-mode theory calculations were done by using the effective index method, calculating the overlap integrals, solving Eq. (11), and reassembling the field. Waveguide height = 500 nm, width = 200 nm, separation = 200 nm, ncore = 3.47, and nclad =1.46. Under nearest neighbor coupling, the scaling relationship predicted by Eq. (13) adequately predicts the field amplitudes within each waveguide.

Fig. 3.
Fig. 3.

TM Polarization Ex : The modes of an N = 5 coupled waveguide array, calculated using coupled-mode theory (blue solid lines), and a finite-difference frequency-domain algorithm (black crosses). The coupled-mode theory calculations were done by using the effective index method, calculating the overlap integrals, solving Eq. (11), and reassembling the field. Waveguide height = 500 nm, width = 200 nm, separation = 1μm, ncore = 3.47, and nclad =1.46. Under nearest neighbor coupling, the scaling relationship predicted by Eq. (13) adequately predicts the field amplitudes within each waveguide.

Fig. 4.
Fig. 4.

Error versus N: Exact eigenvalues of a tridiagonal symmetric matrix of size N were perturbed by values chosen from a uniform random distribution with variance chosen to be ten percent of the first eigenvalue. The variance and mean of the reconstructed nearest-neighbor coupling and next-to-nearest-neighbor coupling coefficients are plotted, calculated from a distribution of coupling matrices generated by 105 iterations, showing that Eq. (20) is a good predictor of the reconstruction accuracy.

Fig. 5.
Fig. 5.

Ratio of coupling coefficients for different separation distances extracted from Eq. (11), which was reconstructed using an algorithm described in the text. (a) TE Polarization An exponential fit expected from a simple nearest-neighbor-coupling theory holds throughout this regime. (b) TM Polarization At a separations less than 450 nm, the ratio deviates significantly from the predicted behavior. (c) TE Polarization The ratio of cross coupling coefficients show that the reconstructed coupling matrix M becomes asymmetric as the waveguide separation is reduced. (d) TM Polarization The asymmetry of the coupling matrix begins at a larger separation.

Fig. 6.
Fig. 6.

Left column: TE Polarization, and Right column: TM polarization. Effective index of the five supermodes for different separation distances with ncore = 3.47, and (a,b) nclad =1.46, (c,d) nclad = 1 (e,f) nclad = 2.05. For each case as the separation between the waveguides increases, the effective indexes of the modes converge to that of the single waveguide. These values are (a) neff = 2.36 and (b) neff = 1.66 for oxide cladding, (c) neff = 2.24 and (d) neff = 1.07 for air cladding, (e) neff = 2.56 and (f) neff = 2.26 for nitride cladding. The shaded regions indicate > 5% deviation of neff for the m=3 supermode from its theoretical value, which as discussed in the text, is predicted by CMT to be independent of the separation distance.

Fig. 7.
Fig. 7.

TM Polarization Ex : The field profile of the fifth eigenmode in the first waveguide. When the separation is decreased below 450 nm, the peak of the field in the high-index rib indicated by the dotted red line in (a) is no longer centered, and the mode shape is considerably altered, thereby changing both κ and n eff. Consequently, CMT can no longer accurately predict the mode coupling.

Fig. 8.
Fig. 8.

TE Polarization Ey : Using the exact solution from a FDFD simulation of a single waveguide, the horizontal cross section is extracted and five copies are shifted from one another so that their separation corresponds to a waveguide separation of 80 nm. (a) These individual waveguide modes are scaled in accordance with Eq. (13) for the fundamental mode (m=1). (b) The sumation of the individual waveguide modes; superimposed is the FDFD solution of the entire five waveguide structure. (c) Zoomed in to just the first waveguide. CMT and FDFD show a shift of the mode towards the center of the waveguide structure. (d-e) The fifth mode, both CMT and FDFD show a shift towards the edge of the waveguide structure however FDFD shows a shift of greater magnitude.

Fig. 9.
Fig. 9.

TM Polarization: finite element simulations (COMSOL) of the transverse field amplitudes, Ex and Ey , and Power flow, Pz for a separation of 1 μm (first Row) and 100 nm (second row) for the fifth mode. For small separations the polarization becomes strongly hybridized and the power is carried at the outer edges.

Tables (1)

Tables Icon

Table 1. An example of a reconstructed coupling (M) matrix from FDFD calculations of eigenmodes and eigenvalues. Si/SiO2, TE polarization, separation s = 350 nm, β 0 = 2.26128 (2π/λ). Although the nearest-neighbor coupling coefficients dominate, the self-coupling and off-tridiagonal coupling terms are non-zero.

Equations (32)

Equations on this page are rendered with MathJax. Learn more.

2 E + ( n 2 k 2 β 2 ) E = [ · E 1 n 2 · ( n 2 E ) ]
( P xx P xy P yx P yy ) ( E x E y ) = β 2 ( E x E y )
P xx E x = x [ 1 n 2 ( n 2 E x ) x ] + 2 E x y 2 + n 2 k 2 E x , P xy E y = x [ 1 n 2 ( n 2 E y ) y ] + 2 E y x y ,
P yx E x = y [ 1 n 2 ( n 2 E x ) x ] + 2 E x y x , P yy E y = y [ 1 n 2 ( n 2 E x ) y ] + 2 E y x 2 + n 2 k 2 E y .
2 E + ω 2 c 2 n 2 ( x , y ) E = 0 ,
E ̂ = y ̂ E ( x , y ) e iβz = y ̂ [ Σ l = 1 N A l 𝓔 1 ( x , y ) ] e iβz .
n 2 ( x , y ) = n s 2 ( x , y ) + Σ l = 1 N Δ n l 2 ( x , y )
( 2 + ω 2 c 2 [ n s 2 ( x , y ) + Σ l = 1 N Δ n l 2 ( x , y ) ] β 2 ) [ Σ l = 1 N A l 𝓔 1 ( x , y ) ] = 0 .
( 2 + ω 2 c 2 [ n s 2 ( x , y ) + Δ n l 2 ( x , y ) ] β l 2 ) 𝓔 1 ( x , y ) = 0
Σ l = 1 N A 1 ( Δ l + ω 2 c 2 Σ m = 1 m l N Δ n m 2 ( x , y ) ) 𝓔 1 ( x , y ) = 0
Δ 1 β l 2 β 2 .
Σ l = 1 N A 1 ( Δ l ∫∫ 𝓔 j * 𝓔 1 dxdy + ω 2 c 2 Σ m = 1 m l N ∫∫ 𝓔 j * Δ n m 2 ( x , y ) 𝓔 1 dxdy ) = 0 , j = 1,2 , , N .
I jl = ∫∫ 𝓔 j * 𝓔 1 dxdy ,
κ jl = ω 2 c 2 Σ m = 1 m l N ∫∫ 𝓔 j * Δ n m 2 ( x , y ) 𝓔 1 dxdy ,
∫∫ 𝓔 j * 𝓔 1 dxdy = 1 .
( β 1 2 + κ 11 Δ 2 I 12 + κ 12 Δ N 1 I 1 , N 1 + κ 1 , N 1 Δ N I 1 N + κ 1 N Δ 1 I 21 + κ 21 β 2 2 + κ 22 Δ N 1 I 2 , N 1 + κ 2 , N 1 Δ N I 2 N + κ 2 N Δ 1 I 31 + κ 31 Δ 2 I 32 + κ 32 Δ N 1 I 3 , N 1 + κ 3 , N 1 Δ N I 3 N + κ 3 N Δ 2 I N 1 + κ N 1 Δ 2 I N 2 + κ N 2 Δ N 1 I 1 , N 1 + κ N , N 1 β N 2 + κ NN ) ( A 1 A 2 A 3 A N )
= ( β 2 0 0 0 0 β 2 0 0 0 0 β 2 0 0 0 0 β 2 ) ( A 1 A 2 A 3 A N ) .
M = ( β 1 2 + κ 11 Δ 2 I 12 + κ 12 0 0 Δ 2 I 21 + κ 21 β 2 2 + κ 22 0 0 0 Δ 2 I 32 + κ 32 0 0 0 Δ N 1 I N , N 1 + κ N , N 1 β N 2 + κ NN ) .
A l ( m ) = ( 2 N + 1 ) 1 / 2 sin lmπ N + 1 δ κ self ( 2 N + 1 ) 3 / 2
× Σ m = 1 n m N 1 β 2 ( m ) β 2 ( m ) [ sin N + 1 sin N + 1 + sin mNπ N + 1 sin nNπ N + 1 ] sin lnπ N + 1
H = y ̂ H ( x , y ) e iβz = y ̂ [ Σ l = 1 N A l H 1 ( x , y ) ] e iβz .
E x = β ε ( x , y ) ω H y ,
Δ M ij = Σ k = 1 N [ ( u ik + u jk ) Δ u k λ k + u ik u jk Δ λ k ] .
Δ M jj + 1 = Σ k = 1 N 2 N + 1 sin jkπ N + 1 sin ( j + 1 ) N + 1 Δ λ k ,
Δ M j 1 j + 1 = Σ k = 1 N 2 N + 1 sin ( j 1 ) N + 1 sin ( j + 1 ) N + 1 Δ λ k .
E [ Δλ M jj + 1 ] = E [ Δλ ] Σ k = 1 N 2 N + 1 sin jkπ N + 1 sin ( j + 1 ) N + 1 = 0
Var [ Δλ M jj + 1 ] = Var [ Δλ ] ( 1 N + 1 ) 2 Σ k = 1 N ( cos 2 N + 1 cos 2 jkπ N + 1 ) 2
= Var [ Δλ ] ( 1 N + 1 ) ,
Nearest neighbor : E [ Δλ M jj + 1 ] = 0 , Var [ Δλ M jj + 1 ] = Var [ Δλ ] / ( N + 1 ) ,
Next to nearest : E [ Δλ M j 1 j + 1 ] = 0 , Var [ Δλ M j 1 j + 1 ] = Var [ Δλ ] / ( N + 1 ) .
Δ 3 I 13 + κ 13 Δ 2 I 12 + κ 13 κ 13 κ 12 = e p ( 2 s + w ) e ps = e p ( s + w )
β 2 ( m ) = β 0 2 + κ self + 2 ( κ + Δ 0 I ) cos N + 1 2 δ κ self N + 1 ( sin 2 N + 1 + sin 2 Nmπ N + 1 ) ,

Metrics